5.数据裁剪

为了减少计算量，先判定数据是否在该区域的最小外接矩形内，再执行点是否在区域内的算法（在本次毕业设计中，效果并不大，可以不用进行）

方法一：射线法

使用射线法
https://www.cnblogs.com/guogangj/p/5127527.html
https://blog.csdn.net/weixin_42943114/article/details/124645273 (更多的方法，更好的讲解)
一个点引出一根“射线”，与多边形的任意若干条边相交，累计相交的边的数目，如果是奇数，那么点就在多边形内，否则点就在多边形外。

方法二：调用matlab函数inpolygon

虽然matlab里的示例给出的是凸多边形，但 https://blog.csdn.net/weixin_42943114/article/details/124645273 试验了可以凹多边形也可以使用

clear
clc
close all
%多边形定义（连线必须按照首尾相接的顺序）
BD=[0,0;0,-1;1,0;2,-1;2,0;2,1;2,2;1,2;1,1;0,1];%凹多边形
%BD=[0,-1;0,-2;1,-1;2,-2;2,-1;2,0;2,3;-1,3;-1,1;...
%    3,1;3,1.5;0,1.5;0,2.5;1,2.5;1,0;0,0];%凹多边形，而且自相交
%BD=[0,1;-1,0;-1,-1;0,-1;1,-1;1,0];%凸多边形
%要判断的点
[X,Y]=meshgrid(-5:0.1:5,-5:0.1:5);
xy2=[X(:),Y(:)];

%方法6 matlab自带方法
[IsInPoly,IsOnBD]=inpolygon(xy2(:,1),xy2(:,2),BD(:,1),BD(:,2));

figure()
hold on
plot(BD(:,1),BD(:,2))
scatter(xy2(:,1),xy2(:,2),24,IsInPoly+2*IsOnBD,'Marker','.')

tic
% 第二种方法，用内置函数 inpolygon
flags = inpolygon(P(:,1),P(:,2),Poly(:,1),Poly(:,2));
In_P = P(flags== 1,:)            % 根据判定，找出在判定区域内的点坐标
toc                              %

补充

试验了matlab的inploygon函数，在边界点数少和完全不沾边的文件时，速度很快，其余时候都很慢，大概30~70s，简直无法忍受，现在我打算先用最小外接矩形把原始数据剪切，然后在交叉点时再计算是否在准确边界内。