点在多边形内算法——判断一个点是否在一个复杂多边形的内部

新页面(new page)介绍了将样条曲线添加到此技术的内容。也可以访问多边形内最短路径页(shortest-path-through-polygonpage)!


图 1

图1显示了一个具有14条边的凹多边形。我们要判断红色点是否在多边形内。

解决方案是将测试点的Y坐标与多边形的每一个点进行比较,我们会得到一个测试点所在的行与多边形边的交点的列表。在这个例子中有8条边与测试点所在的行相交,而有6条边没有相交。如果测试点的两边点的个数都是奇数个则该测试点在多边形内,否则在多边形外。在这个例子中测试点的左边有5个交点,右边有三个交点,它们都是奇数,所以点在多边形内。

(注意:这个算法适用于顺时针和逆时针绘制的多边形。)


图 2

图2显示了多边形自交的情况。在这个例子中多边形的10条边有些互相交叉。这种情况很像汇编语言中的“异或”(XOR)。多边形中重叠的部分剔除。因此测试点在多边形的外面,我们从它的左右两边各有两个交点也可以看出来。


图 3

图3中多边形没有重叠,但是有两条边相交。这种情况下算法也没有问题,任然可以正常工作。


图 4

图4显示了当我们要测试的点所在的扫描行正好穿过多边形的一个顶点,因此扫描行与边a有一个交点,与边b也有一个交点,一共有两个角点,测试点的右边也是同样的情况,那按照我们的算法得到:测试点在多边形外的结论,但这显然是错误的!

要解决这种情况遇到的问题非常简单。边上的点是否与扫描行相交,我们要看边的两个端点是否是在扫描行的两侧,在扫描行上或上方的端点我们把它认为是同一种情况,那图4中边a的一个端点在扫描线的下方,另一个点在扫描线上或上方,所以边a与扫描线相交,而边b的两个端点都在扫描线上或上方,所以边b与扫描线不相交。


图 5

图5显示的多边形上一条边完全与扫描行重合的情况。根据图4中具体描述的边c的一个端点在扫描线的下方另一个端点在扫描线上或上方,所以边c与扫描线有一个交点,而边d的两个端点都在扫描线上或以上,所以无交点,边e也是两个端点都在扫描线上或以上,所以也没交点。


图 6

图6说明了一种特殊的情况(由加州州立理工大学的John David Munch指出)。多边形的一个角刚好落在扫描线上。其实这也没问题,上面的图中只有红色的边与扫描线相交产生交点,所以第一张图有1个交点第二张图有3个交点,交点的个数任然是奇数个,所以测试点在多边形内部。

多边形的边

如果测试点刚好在多边形的边上,则这种算法得到的结果是不确定的;例如结果可能是“在里面”或“不在里面”,这取决于很多不定的因素例如多边形在坐标系统中的方向。(不过这也没啥影响,因为多边形的边是无限小的,and points that fall right on the edge can go either way withouthurting the look of the polygon)

 

C代码例子

// Globals which should be set before calling this function:
//
// int    polySides  =  how many cornersthe polygon has
// float  polyX[]    =  horizontalcoordinates of corners
// float  polyY[]    =  verticalcoordinates of corners
// float  x,y       =  point to be tested
//
// (Globals are used in this example for purposes of speed.  Change as
// desired.)
//
//  Thefunction will return YES if the point x,y is inside the polygon, or
//  NOif it is not.  If the point is exactly on the edge of the polygon,
// then the function may return YES or NO.
//
// Note that division by zero is avoided because the division is protected
//  bythe "if" clause which surrounds it.

bool pointInPolygon() {

  int   i,j=polySides-1 ;
  bool  oddNodes=NO     ;

  for (i=0;i
    if(polyY[i]=y
    ||  polyY[j]=y) {
      if(polyX[i]+(y-polyY[i])/(polyY[j]-polyY[i])*(polyX[j]-polyX[i])
        oddNodes=!oddNodes;}}
    j=i;}

  returnoddNodes; }

 

下面是由Nathan Mercer提供的效力比较的的版本,蓝色的代码eliminatescalculations on sides that are entirely to the right of the test point.。虽然对某些多边形来说可能不是最快的方法,但是大多数情况下是最快的。

// Globals which should be set before calling this function:
//
// int    polySides  =  how many cornersthe polygon has
// float  polyX[]    =  horizontalcoordinates of corners
// float  polyY[]    =  verticalcoordinates of corners
// float  x,y       =  point to be tested
//
// (Globals are used in this example for purposes of speed.  Change as
// desired.)
//
//  Thefunction will return YES if the point x,y is inside the polygon, or
//  NOif it is not.  If the point is exactly on the edge of the polygon,
// then the function may return YES or NO.
//
// Note that division by zero is avoided because the division is protected
//  bythe "if" clause which surrounds it.

bool pointInPolygon() {

  int   i,j=polySides-1 ;
  bool  oddNodes=NO     ;

  for (i=0;i
    if((polyY[i]< y && polyY[j]>=y
    ||   polyY[j]=y)
    &&  (polyX[i]<=x || polyX[j]<=x)){
      if(polyX[i]+(y-polyY[i])/(polyY[j]-polyY[i])*(polyX[j]-polyX[i])
        oddNodes=!oddNodes;}}
    j=i;}

  returnoddNodes; }

这里是由提供了另一个高效的版本。内部的if语句消除或替代了异或操作。

// Globals which should be set before calling this function:
//
// int    polySides  =  how many cornersthe polygon has
// float  polyX[]    =  horizontalcoordinates of corners
// float  polyY[]    =  verticalcoordinates of corners
// float  x,y       =  point to be tested
//
// (Globals are used in this example for purposes of speed.  Change as
// desired.)
//
//  Thefunction will return YES if the point x,y is inside the polygon, or
//  NOif it is not.  If the point is exactly on the edge of the polygon,
// then the function may return YES or NO.
//
// Note that division by zero is avoided because the division is protected
//  bythe "if" clause which surrounds it.

bool pointInPolygon() {

  int   i,j=polySides-1 ;
  bool  oddNodes=NO     ;

  for (i=0;i
    if((polyY[i]< y && polyY[j]>=y
    ||   polyY[j]=y)
    && (polyX[i]<=x || polyX[j]<=x)) {
      oddNodes^=(polyX[i]+(y-polyY[i])/(polyY[j]-polyY[i])*(polyX[j]-polyX[i])
    j=i;}

  returnoddNodes; }

 

整型问题

假如你试图画一个像下图中蓝色的多边形,但是出来的却是有横线和竖线组成的边的多边形,如下图中的红色的多边形。出现这种情况可能是将多边形的顶点坐标变量定义成了整型而非浮点型,仔细地检查你的代码确保多边形的顶点是以浮点型定义且传递的。




转载:http://blog.csdn.net/hjh2005/article/details/9246967


1. 叉乘判别法(只适用于凸多边形)

想象一个凸多边形,其每一个边都将整个2D屏幕划分成为左右两边,连接每一边的第一个端点和要测试的点得到一个矢量v,将两个2维矢量扩展成3维的,然后将该边与v叉乘,判断结果3维矢量中Z分量的符号是否发生变化,进而推导出点是否处于凸多边形内外。这里要注意的是,多边形顶点究竟是左手序还是右手序,这对具体判断方式有影响。

2. 面积判别法(只适用于凸多边形)

第四点分别与三角形的两个点组成的面积分别设为S1,S2,S3,只要S1+S2+S3>原来的三角形面积就不在三角形范围中.可以使用海伦公式 。推广一下是否可以得到面向凸多边形的算法?(不确定)

3. 角度和判别法(适用于任意多边形)

[cpp]  view plain copy
  1. double angle = 0;  
  2. realPointList::iterator iter1 = points.begin();  
  3. for (realPointList::iterator iter2 = (iter1 + 1); iter2 < points.end(); ++iter1, ++iter2)  
  4. {  
  5. double x1 = (*iter1).x - p.x;   
  6. double y1 = (*iter1).y - p.y;   
  7. double x2 = (*iter2).x - p.x;  
  8. double y2 = (*iter2).y - p.y;   
  9. angle += angle2D(x1, y1, x2, y2);  
  10. }  
  11. if (fabs(angle - span::PI2) < 0.01) return true;  
  12. else return false;  
  13. 另外,可以使用bounding box来加速。  
  14. if (p.x < (*iter)->boundingBox.left ||  
  15. p.x > (*iter)->boundingBox.right ||  
  16. p.y < (*iter)->boundingBox.bottom ||  
  17. p.y > (*iter)->boundingBox.top) 。。。。。。  


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

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

4. 水平/垂直交叉点数判别法(适用于任意多边形)

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

[cpp]  view plain copy
  1. BOOL PtInPolygon (POINT p, LPPOINT ptPolygon, int nCount)  
  2.   
  3. {  
  4.   
  5.         // 交点个数  
  6.   
  7.        int nCross = 0;  
  8.   
  9.   
  10.        for (int i = 0; i < nCount; i++)   
  11.   
  12.        {  
  13.   
  14.               POINT p1 = ptPolygon[i];  
  15.   
  16.               POINT p2 = ptPolygon[(i + 1) % nCount];// 最后一个点与第一个点连线  
  17.   
  18.   
  19.               if ( p1.y == p2.y )  
  20.   
  21.                      continue;  
  22.   
  23.               if ( p.y < min(p1.y, p2.y) )  
  24.   
  25.                      continue;  
  26.   
  27.               if ( p.y >= max(p1.y, p2.y) )  
  28.   
  29.                      continue;  
  30.   
  31.   
  32.               // 求交点的x坐标  
  33.   
  34.               double x = (double)(p.y - p1.y) * (double)(p2.x - p1.x) / (double)(p2.y - p1.y) + p1.x;  
  35.   
  36.   
  37.               // 只统计p1p2与p向右射线的交点  
  38.   
  39.               if ( x > p.x )  
  40.   
  41.               {  
  42.   
  43.                      nCross++;  
  44.   
  45.               }  
  46.   
  47.        }  
  48.   
  49.   
  50.        // 交点为偶数,点在多边形之外  
  51.   
  52.        return (nCross % 2 == 1);  
  53.   
  54. }  
转载:http://blog.csdn.net/orange_xxx/article/details/7445803

再经典不过的算法了: 

 

// 功能:判断点是否在多边形内 

// 方法:求解通过该点的水平线与多边形各边的交点 

// 结论:单边交点为奇数,成立

 

//参数: 

// POINT p 指定的某个点 

// LPPOINT ptPolygon 多边形的各个顶点坐标(首末点可以不一致) 

// int nCount 多边形定点的个数 

 

 

BOOL PtInPolygon (POINT p, LPPOINT ptPolygon, int nCount) 

  int nCross = 0; 

 

  for (int i = 0; i < nCount; i++) 

  

    POINT p1 = ptPolygon[i]; 

    POINT p2 = ptPolygon[(i + 1) % nCount]; 

 

    // 求解 y=p.y 与 p1p2 的交点 

 

    if ( p1.y == p2.y ) // p1p2 与 y=p0.y平行 

      continue; 

 

    if ( p.y < min(p1.y, p2.y) ) // 交点在p1p2延长线上 

      continue; 

    if ( p.y >= max(p1.y, p2.y) ) // 交点在p1p2延长线上 

      continue; 

 

    // 求交点的 坐标 -------------------------------------------------------------- 

    double x = (double)(p.y - p1.y) * (double)(p2.x - p1.x) / (double)(p2.y - p1.y) + p1.x; 

 

    if ( x > p.x ) 

      nCross++; // 只统计单边交点 

  

 

   // 单边交点为偶数,点在多边形之外 --- 

   return (nCross % 2 == 1); 




你可能感兴趣的:(算法)