File: poly2cw.m

package info (click to toggle)
octave-geometry 4.1.0-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 720 kB
  • sloc: cpp: 5,358; python: 379; objc: 328; makefile: 25
file content (70 lines) | stat: -rw-r--r-- 2,198 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
## Copyright (C) 2016 - Amr Mohamed
## Copyright (C) 2017 - Piyush Jain
##
## This file is part of Octave.
##
## Octave 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.
##
## Octave 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 Octave; see the file COPYING.  If not,
## see <http://www.gnu.org/licenses/>.

## -*- texinfo -*-
## @deftypefn {Function File} [@var{xcw},@var{ycw}] = poly2cw (@var{x},@var{y})
## Convert Polygons to clockwise contours(polygons).
##
## @var{x}/@var{y} is a cell array or NaN delimited vector of polygons, representing the x/y coordinates of the points. If x1 and y1 can contain multiple contours, represented either as NaN-separated vectors or as cell arrays, then each contour is converted to counter-clockwise ordering.
## @var{xccw}/@var{yccw} has the same format of the input.
##
## @seealso{poly2ccw,ispolycw}
## @end deftypefn

## Created: 2017-07-18

function [xcw, ycw]=poly2cw(x,y);
  if (nargin != 2)
  #case of wrong number of input arguments
   print_usage();
  endif

  if(isempty(x) || isempty(y))
  	error ('Octave:invalid-input-arg', ...
            "poly2ccw: Empty arguments");
  endif

   x = reshape(x, numel(x), 1);
   y = reshape(y, numel(y), 1);

  [xcw, ycw] = orientPolygon(x,y,"cw");

  if(xcw(1) == xcw(2) && ycw(1) == ycw(2))
    xcw = circshift(xcw,-1);
    ycw = circshift(ycw,-1);
  endif

endfunction


%!test
%! x = [0; 1; 1; 0; 0]; y = [0; 0; 1; 1; 0];
%! [xcw,ycw] = poly2cw(x,y);
%! xexp = [0; 0; 1; 1; 0]; yexp = [0; 1; 1; 0; 0];
%! assert (xcw,xexp);
%! assert (ycw,yexp);


%!test
%! x=[0 0 2 2 NaN 0 2 0]; y=[0 2 2 0 NaN 0 0 3];
%! [xcw,ycw]=poly2cw(x,y);
%! xexp=[0; 0; 2; 2; NaN; 0; 0; 2];
%! yexp=[0; 2; 2; 0; NaN; 0; 3; 0];
%! assert (xcw,xexp);
%! assert (ycw,yexp);