Accelerating the pace of engineering and science

• Trials

# inpolygon

Points inside polygonal region

## Syntax

IN = inpolygon(X,Y,xv,yv)
[IN ON] = inpolygon(X,Y,xv,yv)

## Description

IN = inpolygon(X,Y,xv,yv) returns a matrix IN the same size as X and Y. Each element of IN is assigned the value 1 or 0 depending on whether the point (X(p,q),Y(p,q)) is inside the polygonal region whose vertices are specified by the vectors xv and yv. In particular:

 IN(p,q) = 1 If (X(p,q),Y(p,q)) is inside the polygonal region or on the polygon boundary IN(p,q) = 0 If (X(p,q),Y(p,q)) is outside the polygonal region

[IN ON] = inpolygon(X,Y,xv,yv) returns a second matrix ON the same size as X and Y. Each element of ON is assigned the value 1 or 0 depending on whether the point (X(p,q),Y(p,q)) is on the boundary of the polygonal region whose vertices are specified by the vectors xv and yv. In particular:

 ON(p,q) = 1 If (X(p,q),Y(p,q)) is on the polygon boundary ON(p,q) = 0 If (X(p,q),Y(p,q)) is inside or outside the polygon boundary

## Examples

```L = linspace(0,2.*pi,6); xv = cos(L)';yv = sin(L)';
xv = [xv ; xv(1)]; yv = [yv ; yv(1)];
x = randn(250,1); y = randn(250,1);
in = inpolygon(x,y,xv,yv);
plot(xv,yv,x(in),y(in),'r+',x(~in),y(~in),'bo')```