File: poly2ccw.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 (68 lines) | stat: -rw-r--r-- 2,221 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
## 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{xccw},@var{yccw}] = poly2ccw (@var{x},@var{y})
## Convert Polygons to counterclockwise 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{poly2cw,ispolycw}
## @end deftypefn

## Created: 2017-07-18

function [xccw, yccw]=poly2ccw(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);

  [xccw, yccw] = orientPolygon(x,y,"ccw");

  if(xccw(1) == xccw(2) && yccw(1) == yccw(2))
    xccw = circshift(xccw,-1);
    yccw = circshift(yccw,-1);
  endif

endfunction

%!test
%! x = [0 0 1 1 0]; y = [0 1 1 0 0];
%! [xccw,yccw] = poly2ccw(x,y);
%! xexp = [0; 1; 1; 0; 0]; yexp = [0; 0; 1; 1; 0];
%! assert (xccw,xexp);
%! assert (yccw,yexp);

%!test
%! x=[0 0 2 2 NaN 0 2 0]; y=[0 2 2 0 NaN 0 0 3];
%! [xccw,yccw]=poly2ccw(x,y);
%! xexp=[0; 2; 2; 0; NaN; 0; 2; 0];
%! yexp=[0; 0; 2; 2; NaN; 0; 0; 3];
%! assert (xccw,xexp);
%! assert (yccw,yexp);