


ISINPOLYLARGE - indicate if points are inside a polygon
Usage:
isin = isinpolylarge(x,y,xp,yp)
See also ISINPLOY

0001 function isin = isinpolylarge(x,y,xp,yp) 0002 %ISINPOLYLARGE - indicate if points are inside a polygon 0003 % 0004 % Usage: 0005 % isin = isinpolylarge(x,y,xp,yp) 0006 % 0007 % See also ISINPLOY 0008 0009 size=10000; 0010 num=floor(length(x)/size); 0011 isin=zeros(length(x),1); 0012 0013 for i=1:num, 0014 disp(i); 0015 isin((i-1)*size+1:i*size)=isinpoly(x((i-1)*size+1:i*size),y((i-1)*size+1:i*size),xp,yp); 0016 end 0017 isin(num*size+1:length(x))=isinpoly(x(num*size+1:length(x)),x(num*size+1:length(x)),xp,yp); 0018