octave-cvsupdate Mailing List for Octave Forge (Page 27)

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 422-6466

Revision: 10448
          http://octave.svn.sourceforge.net/octave/?rev=10448&view=rev
Author:   jpicarbajal
Date:     2012-05-16 08:00:43 +0000 (Wed, 16 May 2012)
Log Message:
-----------
geometry: Adding file from Simeneon

Modified Paths:
--------------
    trunk/octave-forge/main/geometry/DESCRIPTION
    trunk/octave-forge/main/geometry/INDEX
    trunk/octave-forge/main/geometry/NEWS

Added Paths:
-----------
    trunk/octave-forge/main/geometry/inst/geom2d/closed_path.m

Modified: trunk/octave-forge/main/geometry/DESCRIPTION
===================================================================

--- trunk/octave-forge/main/geometry/DESCRIPTION	2012-05-15 20:40:30 UTC (rev 10447)
+++ trunk/octave-forge/main/geometry/DESCRIPTION	2012-05-16 08:00:43 UTC (rev 10448)
@@ -1,7 +1,7 @@
 Name: Geometry
 Version: 1.5.0
 Date: 2012-04-16
-Author: David Legland <dav...@gr...>, José Luis García Pallero <jgp...@gm...>, Juan Pablo Carbajal <car...@if...>
+Author: David Legland <dav...@gr...>, José Luis García Pallero <jgp...@gm...>, Juan Pablo Carbajal <car...@if...>, Simeon Simeonov <ss...@ca...>
 Maintainer: Juan Pablo Carbajal <car...@if...>
 Title: Computational Geometry
 Description: Library for geometric computing extending MatGeom functions. Useful to create, transform, manipulate and display geometric primitives.

Modified: trunk/octave-forge/main/geometry/INDEX
===================================================================
--- trunk/octave-forge/main/geometry/INDEX	2012-05-15 20:40:30 UTC (rev 10447)
+++ trunk/octave-forge/main/geometry/INDEX	2012-05-16 08:00:43 UTC (rev 10448)
@@ -14,6 +14,7 @@
  polygons2d
 2D Points
  centroid
+ closed_path
  distancePoints
  drawPoint
  isCounterClockwise

Modified: trunk/octave-forge/main/geometry/NEWS
===================================================================
--- trunk/octave-forge/main/geometry/NEWS	2012-05-15 20:40:30 UTC (rev 10447)
+++ trunk/octave-forge/main/geometry/NEWS	2012-05-16 08:00:43 UTC (rev 10448)
@@ -30,6 +30,9 @@
    expandPolygon.m, medialAxisConvex.m, polygonLoops.m, polygonSelfIntersections.m
    polylineSelfIntersections.m, splitPolygons.m
 
+ - close_path.m : given a set of points in the plane calculate a piecewise linear
+                  simple path that passes through all points.
+
 * Changed functions:
  - distancePointEdge.m: Now the function computes the distance between all points
                         and all edges. A third optional argument provides

Added: trunk/octave-forge/main/geometry/inst/geom2d/closed_path.m
===================================================================
--- trunk/octave-forge/main/geometry/inst/geom2d/closed_path.m	                        (rev 0)
+++ trunk/octave-forge/main/geometry/inst/geom2d/closed_path.m	2012-05-16 08:00:43 UTC (rev 10448)
@@ -0,0 +1,105 @@
+%% Copyright (C) 2012 Simeon Simeonov <sim...@gm...>
+%%
+%% This program is free software; you can redistribute it and/or modify
+%% it under the terms of the GNU General Public License as published by
+%% the Free Software Foundation; either version 3 of the License, or
+%% (at your option) any later version.
+%%
+%% This program is distributed in the hope that it will be useful,
+%% but WITHOUT ANY WARRANTY; without even the implied warranty of
+%% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+%% GNU General Public License for more details.
+%%
+%% You should have received a copy of the GNU General Public License
+%% along with this program; if not, see <http://www.gnu.org/licenses/>.
+%%
+%% -*- texinfo -*-
+%% @deftypefn {Function File} {@var{y} =} polygon (@var{x})
+%% Returns a simple closed path that passes through all the points in x.
+%% x is a vector containing 2D coordinates of the points.
+%%
+%% @end deftypefn
+
+%% Author: Simeon Simeonov <sim...@gm...>
+
+function y = closed_path(x)
+
+  if(size(x,1) > 1 && size(x,1) < size(x,2))
+      x = x';
+  end
+
+  N = size(x,1); % Number of points
+  idx = zeros(N, 1); % ind contains the indices of the sorted coordinates
+
+  a = find(x(:,2)==min(x(:,2)));
+
+  if(size(a,1) > 1)
+      [~, i] = sort(x(a,1));
+      a = a(i(1));
+  end
+
+  x_1 = x(x(:,2)==x(a,2),:); % find all x with the same y coordinate
+
+  if(x(a,1) == min(x(:,1)))
+      x_2 = x(x(:,1)==x(a,1),:); % find all x with the same x coordinate
+  else
+     x_2 = x(a,:);
+  end
+
+  if(size(x_1,1) > 1 || size(x_2,1) > 1)
+      if(size(x_1,1) > 1)
+          x_1 = sort(x_1); % Sort by x coordinate
+          y(1,:) = x(a,:); % original starting point
+      end
+
+      if (size(x_2,1) > 1)
+          x_2 = sort(x_2, 'descend');
+      end
+
+      x_not = [x_1; x_2];
+      i = ismember(x,x_not,'rows');
+      x(i, :) = [];
+      x = [x_1(size(x_1,1),:); x];
+      x_1(size(x_1, 1),:) = [];
+      N = size(x,1);
+      a = 1;
+  else
+      x_1 = [];
+      x_2 = x(a,:);
+  end
+  d = x - repmat(x(a,:), N, 1);
+  th = d(:,2)./(d(:,1) + d(:,2));
+
+  [~, idx0] = ismember(sort(th(th==0)), th);
+  [~, idx1] = ismember(sort(th(th>0)), th);
+  [~, idx2] = ismember(sort(th(th<0)), th);
+
+  idx = [a; idx0; idx1; idx2];
+  % I contains the indices of idx in a sorted order. [v i] = sort(idx) then
+  % i==I.
+  [~,I,J]= unique(idx);
+  if(size(I,1) ~= size(J,1))
+      R = histc(J, 1:size(I,1)); % R(1) will always be 1?
+      idx_sorted = idx(I);
+      r = find(R>1);
+      for ri = r'
+          idx_repeated = idx_sorted(ri);
+          idx(idx==idx_repeated) = find(th==th(idx_sorted(ri)));
+      end
+  end
+
+  y = [x_1; x(idx,:); x_2;];
+
+endfunction
+
+%!demo
+%! maxInt = 100;
+%! N = 25;
+%!
+%! for i = 1:5
+%!   x = randi(maxInt, N, 2);
+%!   y = closed_path(x);
+%!   plot(y(:,1), y(:,2), '*');
+%!   hold on;
+%!   plot(y(:,1), y(:,2));
+%! end

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.





2001	Jan	Feb	Mar	Apr	May	Jun	Jul (10)	Aug (5)	Sep (3)	Oct (41)	Nov (41)	Dec (33)
2002	Jan (75)	Feb (10)	Mar (170)	Apr (174)	May (66)	Jun (11)	Jul (10)	Aug (44)	Sep (73)	Oct (28)	Nov (139)	Dec (52)
2003	Jan (35)	Feb (93)	Mar (62)	Apr (10)	May (55)	Jun (70)	Jul (37)	Aug (16)	Sep (56)	Oct (31)	Nov (57)	Dec (83)
2004	Jan (85)	Feb (67)	Mar (27)	Apr (37)	May (75)	Jun (85)	Jul (160)	Aug (68)	Sep (104)	Oct (25)	Nov (39)	Dec (23)
2005	Jan (10)	Feb (45)	Mar (43)	Apr (19)	May (108)	Jun (31)	Jul (41)	Aug (23)	Sep (65)	Oct (58)	Nov (44)	Dec (54)
2006	Jan (96)	Feb (27)	Mar (69)	Apr (59)	May (67)	Jun (35)	Jul (13)	Aug (461)	Sep (160)	Oct (399)	Nov (32)	Dec (72)
2007	Jan (316)	Feb (305)	Mar (318)	Apr (54)	May (194)	Jun (173)	Jul (282)	Aug (91)	Sep (227)	Oct (365)	Nov (168)	Dec (18)
2008	Jan (71)	Feb (111)	Mar (155)	Apr (173)	May (70)	Jun (67)	Jul (55)	Aug (83)	Sep (32)	Oct (68)	Nov (80)	Dec (29)
2009	Jan (46)	Feb (18)	Mar (95)	Apr (76)	May (140)	Jun (98)	Jul (84)	Aug (123)	Sep (94)	Oct (131)	Nov (142)	Dec (125)
2010	Jan (128)	Feb (158)	Mar (172)	Apr (134)	May (94)	Jun (84)	Jul (32)	Aug (127)	Sep (167)	Oct (109)	Nov (69)	Dec (78)
2011	Jan (39)	Feb (58)	Mar (52)	Apr (47)	May (56)	Jun (76)	Jul (55)	Aug (54)	Sep (165)	Oct (255)	Nov (328)	Dec (263)
2012	Jan (82)	Feb (147)	Mar (400)	Apr (216)	May (209)	Jun (160)	Jul (86)	Aug (141)	Sep (156)	Oct (6)	Nov	Dec
2015	Jan	Feb	Mar	Apr	May	Jun	Jul (1)	Aug	Sep (1)	Oct	Nov (1)	Dec (2)
2016	Jan	Feb (2)	Mar (2)	Apr (1)	May (1)	Jun (2)	Jul (1)	Aug (1)	Sep	Oct	Nov (1)	Dec
2019	Jan	Feb	Mar (1)	Apr	May	Jun	Jul	Aug	Sep	Oct	Nov	Dec
2021	Jan	Feb	Mar	Apr (3)	May (4)	Jun (8)	Jul (2)	Aug (5)	Sep (9)	Oct	Nov	Dec

octave-cvsupdate Mailing List for Octave Forge (Page 27)

A collection of packages providing extra functionality for GNU Octave

octave-cvsupdate — CVS update notifications