Mercurial > hg > octave-lyh
annotate scripts/geometry/inpolygon.m @ 17157:73a3c1580974
Correctly identify PNG image types (bug #32986, bug #36820 and bug #37689)
* __magick_read_.cc (read_images): check original color type from PNG image
attribute because GraphicsMagick changes it for optimization.
author | Carnë Draug <carandraug@octave.org> |
---|---|
date | Fri, 02 Aug 2013 19:05:42 +0100 |
parents | 5d3a684236b0 |
children | f4c8c66faf34 |
rev | line source |
---|---|
14138
72c96de7a403
maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents:
13057
diff
changeset
|
1 ## Copyright (C) 2006-2012 Frederick (Rick) A Niles |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
2 ## and S�ren Hauberg |
6847 | 3 ## |
4 ## This file is part of Octave. | |
5 ## | |
6 ## Octave is free software; you can redistribute it and/or modify it | |
7 ## under the terms of the GNU General Public License as published by | |
7016 | 8 ## the Free Software Foundation; either version 3 of the License, or (at |
9 ## your option) any later version. | |
6847 | 10 ## |
11 ## Octave is distributed in the hope that it will be useful, but | |
12 ## WITHOUT ANY WARRANTY; without even the implied warranty of | |
13 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
14 ## General Public License for more details. | |
15 ## | |
16 ## You should have received a copy of the GNU General Public License | |
7016 | 17 ## along with Octave; see the file COPYING. If not, see |
18 ## <http://www.gnu.org/licenses/>. | |
6847 | 19 |
20 ## -*- texinfo -*- | |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
21 ## @deftypefn {Function File} {[@var{in}, @var{on}] =} inpolygon (@var{x}, @var{y}, @var{xv}, @var{yv}) |
6847 | 22 ## |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
23 ## For a polygon defined by vertex points @code{(@var{xv}, @var{yv})}, determine |
6847 | 24 ## if the points @code{(@var{x}, @var{y})} are inside or outside the polygon. |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
25 ## The variables @var{x}, @var{y}, must have the same dimension. The optional |
6847 | 26 ## output @var{on} gives the points that are on the polygon. |
27 ## | |
28 ## @end deftypefn | |
29 | |
30 ## Author: Frederick (Rick) A Niles <niles@rickniles.com> | |
31 ## Created: 14 November 2006 | |
32 | |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
33 ## Vectorized by S�ren Hauberg <soren@hauberg.org> |
6847 | 34 |
35 ## The method for determining if a point is in in a polygon is based on | |
36 ## the algorithm shown on | |
37 ## http://local.wasp.uwa.edu.au/~pbourke/geometry/insidepoly/ and is | |
38 ## credited to Randolph Franklin. | |
39 | |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
40 function [in, on] = inpolygon (x, y, xv, yv) |
6847 | 41 |
7125 | 42 if (nargin != 4) |
43 print_usage (); | |
44 endif | |
45 | |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
46 if (! (isreal (x) && isreal (y) && ismatrix (y) && ismatrix (y) |
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
47 && size_equal (x, y))) |
6847 | 48 error ("inpolygon: first two arguments must be real matrices of same size"); |
7125 | 49 elseif (! (isreal (xv) && isreal (yv) && isvector (xv) && isvector (yv) |
10549 | 50 && size_equal (xv, yv))) |
6847 | 51 error ("inpolygon: last two arguments must be real vectors of same size"); |
52 endif | |
53 | |
7125 | 54 npol = length (xv); |
6847 | 55 do_boundary = (nargout >= 2); |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
56 |
14868
5d3a684236b0
maint: Use Octave coding conventions for cuddling parentheses in scripts directory
Rik <octave@nomad.inbox5.com>
parents:
14237
diff
changeset
|
57 in = zeros (size (x), "logical"); |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
58 if (do_boundary) |
14868
5d3a684236b0
maint: Use Octave coding conventions for cuddling parentheses in scripts directory
Rik <octave@nomad.inbox5.com>
parents:
14237
diff
changeset
|
59 on = zeros (size (x), "logical"); |
6847 | 60 endif |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
61 |
6847 | 62 j = npol; |
63 for i = 1 : npol | |
64 delta_xv = xv(j) - xv(i); | |
65 delta_yv = yv(j) - yv(i); | |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
66 ## distance = [distance from (x,y) to edge] * length(edge) |
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
67 distance = delta_xv .* (y - yv(i)) - (x - xv(i)) .* delta_yv; |
6847 | 68 ## |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
69 ## is y between the y-values of edge i,j |
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
70 ## AND (x,y) on the left of the edge ? |
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
71 idx1 = (((yv(i) <= y & y < yv(j)) | (yv(j) <= y & y < yv(i))) |
10549 | 72 & 0 < distance.*delta_yv); |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
73 in (idx1) = !in (idx1); |
6847 | 74 |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
75 ## Check if (x,y) are actually on the boundary of the polygon. |
6847 | 76 if (do_boundary) |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
77 idx2 = (((yv(i) <= y & y <= yv(j)) | (yv(j) <= y & y <= yv(i))) |
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
78 & ((xv(i) <= x & x <= xv(j)) | (xv(j) <= x & x <= xv(i))) |
10549 | 79 & (0 == distance | !delta_xv)); |
11469
c776f063fefe
Overhaul m-script files to use common variable name between code and documentation.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
80 on (idx2) = true; |
6847 | 81 endif |
82 j = i; | |
83 endfor | |
7125 | 84 |
6847 | 85 endfunction |
86 | |
87 | |
88 %!demo | |
14237
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
89 %! xv = [ 0.05840, 0.48375, 0.69356, 1.47478, 1.32158, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
90 %! 1.94545, 2.16477, 1.87639, 1.18218, 0.27615, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
91 %! 0.05840 ]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
92 %! yv = [ 0.60628, 0.04728, 0.50000, 0.50000, 0.02015, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
93 %! 0.18161, 0.78850, 1.13589, 1.33781, 1.04650, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
94 %! 0.60628 ]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
95 %! xa = [0:0.1:2.3]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
96 %! ya = [0:0.1:1.4]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
97 %! [x,y] = meshgrid (xa, ya); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
98 %! [in,on] = inpolygon (x, y, xv, yv); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
99 %! inside = in & !on; |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
100 %! |
14237
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
101 %! clf; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
102 %! plot (xv, yv); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
103 %! hold on; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
104 %! plot (x(inside), y(inside), "@g") |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
105 %! plot (x(!in), y(!in), "@m"); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
106 %! plot (x(on), y(on), "@b"); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
107 %! hold off; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
108 %! disp ("Green points are inside polygon, magenta are outside,"); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
109 %! disp ("and blue are on boundary."); |
6847 | 110 |
14237
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
111 %!demo |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
112 %! xv = [ 0.05840, 0.48375, 0.69356, 1.47478, 1.32158, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
113 %! 1.94545, 2.16477, 1.87639, 1.18218, 0.27615, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
114 %! 0.05840, 0.73295, 1.28913, 1.74221, 1.16023, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
115 %! 0.73295, 0.05840 ]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
116 %! yv = [ 0.60628, 0.04728, 0.50000, 0.50000, 0.02015, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
117 %! 0.18161, 0.78850, 1.13589, 1.33781, 1.04650, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
118 %! 0.60628, 0.82096, 0.67155, 0.96114, 1.14833, ... |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
119 %! 0.82096, 0.60628]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
120 %! xa = [0:0.1:2.3]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
121 %! ya = [0:0.1:1.4]; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
122 %! [x,y] = meshgrid (xa, ya); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
123 %! [in,on] = inpolygon (x, y, xv, yv); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
124 %! inside = in & !on; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
125 %! |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
126 %! clf; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
127 %! plot (xv, yv); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
128 %! hold on; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
129 %! plot (x(inside), y(inside), "@g"); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
130 %! plot (x(!in), y(!in), "@m"); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
131 %! plot (x(on), y(on), "@b"); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
132 %! hold off; |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
133 %! disp ("Green points are inside polygon, magenta are outside,"); |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
134 %! disp ("and blue are on boundary."); |
13057
85c77dd2e85b
codesprint: tests for inpolygon
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
135 |
85c77dd2e85b
codesprint: tests for inpolygon
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
136 %!test |
85c77dd2e85b
codesprint: tests for inpolygon
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
137 %! [in, on] = inpolygon ([1, 0], [1, 0], [-1, -1, 1, 1], [-1, 1, 1, -1]); |
85c77dd2e85b
codesprint: tests for inpolygon
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
138 %! assert (in, [false, true]); |
85c77dd2e85b
codesprint: tests for inpolygon
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
139 %! assert (on, [true, false]); |
14237
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
140 |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
141 %% Test input validation |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
142 %!error inpolygon () |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
143 %!error inpolygon (1, 2) |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
144 %!error inpolygon (1, 2, 3) |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
145 %!error inpolygon (1, [1,2], [3, 4], [5, 6]) |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
146 %!error inpolygon ([1,2], [3, 4], [5, 6], 1) |
11949c9795a0
Revamp %!demos in m-files to use Octave coding conventions on spacing, etc.
Rik <octave@nomad.inbox5.com>
parents:
14222
diff
changeset
|
147 |