判断点在多边形内

1. 角度和判别法（适用于任意多边形）

double angle = 0;

realPointList::iterator iter1 = points.begin();

for (realPointList::iterator iter2 = (iter1 + 1); iter2 < points.end();++iter1, ++iter2)

{

double x1 = (*iter1).x - p.x;

double y1 = (*iter1).y - p.y;

double x2 = (*iter2).x - p.x;

double y2 = (*iter2).y - p.y;

angle += angle2D(x1, y1, x2, y2);

}

if (fabs(angle - span::PI2) < 0.01) return true;

else return false;

另外，可以使用bounding box来加速。

if (p.x < (*iter)->boundingBox.left ||

p.x > (*iter)->boundingBox.right ||

p.y < (*iter)->boundingBox.bottom ||

p.y > (*iter)->boundingBox.top) 。。。。。。

对于多边形来说，计算bounding box非常的简单。只需要把水平和垂直方向上的最大最小值找出来就可以了。

1. 角度和判别法（适用于任意多边形）

double angle = 0;
realPointList::iterator iter1 = points.begin();
for (realPointList::iterator iter2 = (iter1 + 1); iter2 < points.end();++iter1, ++iter2)
{
double x1 = (*iter1).x - p.x;
double y1 = (*iter1).y - p.y;
double x2 = (*iter2).x - p.x;
double y2 = (*iter2).y - p.y;
angle += angle2D(x1, y1, x2, y2);
}

if (fabs(angle - span::PI2) < 0.01) return true;
else return false;

另外，可以使用bounding box来加速。
if (p.x < (*iter)->boundingBox.left ||
p.x > (*iter)->boundingBox.right ||
p.y < (*iter)->boundingBox.bottom ||
p.y > (*iter)->boundingBox.top) 。。。。。。

对于多边形来说，计算bounding box非常的简单。只需要把水平和垂直方向上的最大最小值找出来就可以了。

对于三角形：第四点分别与三角形的两个点的交线组成的角度分别设为j1,j2,j3，只要j1+j2+j3>360就不在三角形范围中。

2. 水平/垂直交叉点数判别法（适用于任意多边形）

注意到如果从P作水平向左的射线的话，如果P在多边形内部，那么这条射线与多边形的交点必为奇数，如果P在多边形外部，则交点个数必为偶数（0也在内）。所以，我们可以顺序考虑多边形的每条边，求出交点的总个数。还有一些特殊情况要考虑。假如考虑边(P1,P2)，
1)如果射线正好穿过P1或者P2,那么这个交点会被算作2次，处理办法是如果P的从坐标与P1,P2中较小的纵坐标相同，则直接忽略这种情况
2)如果射线水平，则射线要么与其无交点，要么有无数个，这种情况也直接忽略。
3)如果射线竖直，而P0的横坐标小于P1,P2的横坐标，则必然相交。
4)再判断相交之前，先判断P是否在边(P1,P2)的上面，如果在，则直接得出结论：P再多边形内部。

对于三角形：第四点分别与三角形的两个点的交线组成的角度分别设为j1,j2,j3，只要j1+j2+j3>360就不在三角形范围中。

2. 水平/垂直交叉点数判别法（适用于任意多边形）

注意到如果从P作水平向左的射线的话，如果P在多边形内部，那么这条射线与多边形的交点必为奇数，如果P在多边形外部，则交点个数必为偶数（0也在内）。所以，我们可以顺序考虑多边形的每条边，求出交点的总个数。还有一些特殊情况要考虑。假如考虑边(P1,P2)，

1)如果射线正好穿过P1或者P2,那么这个交点会被算作2次，处理办法是如果P的从坐标与P1,P2中较小的纵坐标相同，则直接忽略这种情况

2)如果射线水平，则射线要么与其无交点，要么有无数个，这种情况也直接忽略。

3)如果射线竖直，而P0的横坐标小于P1,P2的横坐标，则必然相交。

4)再判断相交之前，先判断P是否在边(P1,P2)的上面，如果在，则直接得出结论：P再多边形内部。

==================================================================================

int CUtility::IsPtOnGon(const AcGePoint2dArray& ptArr,AcGePoint2d& pt)
{
if(ptArr.length()<=2) return 0;
//判是否在线上===================================================
AcGeCurve2d* pCurve = (AcGeCurve2d*) new AcGePolyline2d(ptArr);
if(pCurve->isOn(pt))
{
delete pCurve;
pCurve= NULL;
return 1;
}
delete pCurve;
pCurve = NULL;
//三角形===========================================================
if(ptArr.length() == 3)
{
AcGeVector2d a,b,c;
a = ptArr[0] - pt;
b = ptArr[1] - pt;
c = ptArr[2] - pt;
double angle = a.angleTo(b);
angle += b.angleTo(c);
angle += c.angleTo(a);

if(fabs(angle - PI*2.0) < 0.001) return 1;
return 0;
}

//多边形===============================================================
//判是否在Extend外
AcGePoint2d ptMin,ptMax;
FindExtInPtArr(ptArr,ptMin,ptMax);
if(pt.x - ptMin.x < 0.0001 || pt.x - ptMax.x > 0.0001) return 0;
if(pt.y - ptMin.y < 0.0001 || pt.y - ptMax.y > 0.0001) return 0;

//点到各顶点的夹角和=2PI
AcGeVector3dArray vArr;
AcGeVector3d pRef;
for(int i=0;i<ptArr.length();i++)
{
AcGeVector3d pVec = AcGePoint3d(ptArr[i].x,ptArr[i].y,0) - AcGePoint3d(pt.x,pt.y,0);
vArr.append(pVec);
}
double angle = 0.0;
int j=0;
for(int i=0;i<vArr.length()-1;i++)
{
pRef = vArr[i].crossProduct(vArr[i+1]);
pRef.normalize();
angle += -pRef.z * vArr[i].angleTo(vArr[i+1]);
j++;
}
pRef = vArr[j].crossProduct(vArr[0]);
pRef.normalize();

angle += -pRef.z * vArr[j].angleTo(vArr[0],pRef);
double sumAngle = PI*2.0;
if(fabs(fabs(angle) - sumAngle) < 0.001) return 1;
else return 0;
}