@Daniele,
I recommend using meshgrid() to make the 2D mesh that covers and extends beyond the plygon of interest, the use reshape() to turn the grids into vectors, and inpolygon() to exclude points that fall outside the polygon of interest.
Example: Create a set of points with 0.1 degree spacing inside the hexagon with corners at [-75,41; -74,40; -74,39; -75,38; -76,39; -76,40]
[X,Y]=meshgrid([-76.2:0.1:-73.8], [37.8:0.1:41.2]);
fprintf('Size of X: %d x %d.\n',size(X))
X=reshape(X,[],1); % reshape X to have 1 column
Y=reshape(Y,[],1); % reshape Y to have 1 column
fprintf('Size of X: %d x %d.\n',size(X))
Define the polygon:
xv=[-75,-74,-74,-75,-76,-76]; % x vertices of polygon
yv=[41,40,39,38,39,40]; % y vertices of polygon
Find who is in and who is out:
in=inpolygon(X,Y,xv,yv);
Xin=X(in); Yin=Y(in); % points inside and on the polygon
Xout=X(~in); Yout=Y(~in); % points outside the polygon
Plot the in and out points in green and red
plot(Xin,Yin,'go');
hold on; axis equal
plot(Xout,Yout,'rx');
OK
