关于几何（CAD二次开发实现判断点与多边形相对关系）

 一、本节课程
 

C++ ARX二次开发点和闭合多段线的位置关系

二、本节要讲解的知识点

结合自己的业务场景，想想自己开发中可能碰到的点和闭合多段线相对关系的需求：统计多段线内部的所有图块；还可以拓展判断任意曲线跟多段线的相对关系。

三、具体内容

1. 计算机图形学、计算几何。判断点是否在多边形内部，一般有以下方法：

1. 叉乘判断（适合凸多边形）：如果将多边形的所有顶点按逆时针排序，那么判断点和每一条边的位置关系，如果点在多边形每一条边的左侧，那么点就在我们的多边形内部。
2. 射线法：从给定点出发，沿着X轴正方向或者负方向做一条射线（射线可能跟多边形没有交点），计算射线跟多边形的交点数量，如果是奇数个交点，在内部；偶数个交点在外部。处理下点就在多边形的顶点上的特例。

1. 本节提供的算法的思路是：

1. 计算点到多段线的最近点，如果两点之间的距离少于容差，则认为点就在多段线上。
2. 从给定点（要判断的点）出发，沿最近点到给定点（要判断的点）的方向 做一个射线（之所以这么来做射线，是为了尽可能早可能早判断出点不在多边形内部的情况），计算射线和多段线的交点。如果交点中包含了多段线的顶点，就将射线旋转2度继续进行判断；如果不包含多段线的顶点，则根据交点个数来判断点和多段线的位置关系。

工具函数的添加：

1. MathUtil.h头文件中添加：

static bool IsEqual(double a,double b,double tol=1.0*10E-7); static int GetRand(int minValue,int maxValue); static double Round(double a,int precision); static int Round(double a);

1. MathUtil.cpp文件中添加：

bool CMathUtil::IsEqual(double a,double b,double tol/*=1.0*10E-7*/) { return  (fabs(a-b)<tol); } //对一个数按照指定的小数位数进行四舍五入 double CMathUtil::Round(double a,int precision) { return (floor(a*pow(10.0,precision)+0.5))/pow(10.0,precision); }  int CMathUtil::Round(double a) { return (int)(a+0.5); } int CMathUtil::GetRand(int minValue,int maxValue)  { assert(maxValue-minValue>0); int value=rand(); int rc=minValue+(int)CMathUtil::Round(((double)value)/RAND_MAX*(maxValue-minValue));//将生成的随机数映射到区间上。 return rc; }

 

1. GePointUtil.h头文件中添加：

static bool IsEqual(const AcGePoint3d &firstPoint,const AcGePoint3d &secondPoint,double tol=1.0E-7); static bool IsEqual(const AcGePoint2d &firstPoint,const AcGePoint2d &secondPoint,double tol=1.0E-7); static int FindPoint(const AcGePoint3dArray &points,const AcGePoint3d &point,double tol=1.0E-7); static int FindPoint(const AcGePoint2dArray &points,const AcGePoint2d &point,double tol=1.0E-7); static void FilterEqualPoints(AcGePoint3dArray &points,double tol=1.0E-7); static void FilterEqualPoints(AcGePoint3dArray &points,const AcGePoin