How to return an array without points in the edge?

Dear guys! How to determine points which are not on the edge?

I used this function:

[in,on] = inpolygon (kn, xq ,xv, yv);

But then how to remove from in points contained in on?

[Adam Danz]

The second output "on" it true for points that are on the edges.

https://www.mathworks.com/help/matlab/ref/inpolygon.html#buanwyg-1-on

Answer 1

[Adam Danz]

Use the logical outputs to the inpolygon function.

[in,on] = inpolygon (kn, xq ,xv, yv);

• on indicates all data that are on the edges of the polygon.
• ~on indicates all data that are inside and outside of the polygon but not on the edges.
• in indicates all data that are in the polygon but not on the edges.
• ~in & ~on indicates all data not in the polygon nor on the edges.

Demo:

L = linspace(0,2*pi,6);
xv = cos(L)';
yv = sin(L)';
rng default
xq = [randn(60,1); xv(1:2)]; 
yq = [randn(60,1); yv(1:2)]; 


[in,on] = inpolygon(xq,yq,xv,yv);

clf
plot(xv, yv, 'k-', 'LineWidth', 1, 'DisplayName', 'Polygon')
hold on 
plot(xq(in),yq(in),'c*','DisplayName','in')                 %plot data IN polygon
plot(xq(on),yq(on),'ro','LineWidth',2,'DisplayName','on')   %plot data on border  
plot(xq(~in&~on),yq(~in&~on),'k.','DisplayName','out')      %plot data outside of polygon
axis equal      
legend()

To remove points points on the edge, you can replace thier values with NaN (1st line below) or remove them (2nd line below).

xq(on) = nan;  yq(on) = nan;    % replace with NaN
xq(on) = [];  yq(on) = [];      % -OR- remove