地图上点与多边形的关系判断

大家都知道，点与多边形的关系无非是点在多边形内、点在多边形外、点在多边形的边界上。那么，如何用代码是实现肉眼就能辨别的地理位置关系呢？~

---

多边形的外包矩形

> 如果一个点不在多边形的外包矩形内，那它肯定不在这个多边形内

射线法判断奇偶数

> 计算射线与多边形各边的交点，如果是偶数，则点在多边形外，否则在多边形内

  //判断点是否在多边形内：
       function isPointInPolygon(point,polygon){

        //如果是百度的点，用百度已有的方法判断该点是否在多边形的外包矩形内
        if((point instanceof BMap.Point)&&(polygon instanceof BMap.Polygon)){

              var polygonBounds = polygon.getBounds();
              var sw = polygonBounds.getSouthWest(); //西南脚点
              var ne = polygonBounds.getNorthEast(); //东北脚点

              var inRect = point.lng >= sw.lng && point.lng <= ne.lng && point.lat >= sw.lat && point.lat <= ne.lat ;

              if(!inRect){

               return false ; //如果不在外包矩形内，更不会在多边形内。直接返回false 

              }else{


                  var pts = polygon.getPath();//获取多边形的点集合。
                  /**
                  基本思想是利用射线法，计算射线与多边形各边的交点，如果是偶数，则点在多边形外，否则在多边形内
                  */
                  var N = pts.length;
                  var boundOrVertex = true; //如果点位于多边形的顶点或边上，也算做点在多边形内，直接返回true
                  var intersectCount = 0;//cross points count of x 
                  var precision = 2e-10; //浮点类型计算时候与0比较时候的容差
                  var p1, p2;//neighbour bound vertices
                  var p = point; //测试点
                  p1 = pts[0];//left vertex 

                  for(var i = 1; i <= N; ++i){//check all rays            

                   if(p.equals(p1)){
                        return boundOrVertex;//p is an vertex
                    }

                   p2 = pts[i % N];//right vertex            

                   if(p.lat < Math.min(p1.lat, p2.lat) || p.lat > Math.max(p1.lat, p2.lat)){//ray is outside of our interests                
                    p1 = p2; 
                    continue;//next ray left point
                   }

                   if(p.lat > Math.min(p1.lat, p2.lat) && p.lat < Math.max(p1.lat, p2.lat)){//ray is crossing over by the algorithm (common part of)
                       if(p.lng <= Math.max(p1.lng, p2.lng)){//x is before of ray                    
                          if(p1.lat == p2.lat && p.lng >= Math.min(p1.lng, p2.lng)){//overlies on a horizontal ray
                             return boundOrVertex;
                    }

                       if(p1.lng == p2.lng){//ray is vertical                        
                         if(p1.lng == p.lng){//overlies on a vertical ray
                            return boundOrVertex;
                          }else{//before ray
                            ++intersectCount;
                          } 
                        }else{//cross point on the left side                        
                       var xinters = (p.lat - p1.lat) * (p2.lng - p1.lng) / (p2.lat - p1.lat) + p1.lng;//cross point of lng                        
                       if(Math.abs(p.lng - xinters) < precision){//overlies on a ray
                            return boundOrVertex;
                        }

                        if(p.lng < xinters){//before ray
                            ++intersectCount;
                        } 
                    }
                 }
                 }else{//special case when ray is crossing through the vertex                

                    if(p.lat == p2.lat && p.lng <= p2.lng){//p crossing over p2                    
                    var p3 = pts[(i+1) % N]; //next vertex                    
                      if(p.lat >= Math.min(p1.lat, p3.lat) && p.lat <= Math.max(p1.lat, p3.lat)){//p.lat lies between p1.lat & p3.lat
                        ++intersectCount;
                      }else{
                        intersectCount += 2;
                     }
                 }
               }            

               p1 = p2;//next ray left point
               }

              if(intersectCount % 2 == 0){//偶数在多边形外
                 return false;
              } else { //奇数在多边形内
                 return true;
               }  

              }
         }
       }

其中3个点的坐标如下：

  var pointA = new BMap.Point(121.333777,31.206659);

  var pointB  = new BMap.Point(121.307313,31.203956);

  var pointC  = new BMap.Point(121.317967,31.224201);

多边形的可参考上一篇计算面积的点集合。

   //调用计算面积的方法
       setTimeout(function(){

        alert("A是否在范围内："+isPointInPolygon(pointA,polygon) );
        alert("B是否在范围内："+isPointInPolygon(pointB,polygon) );
        alert("C是否在范围内："+isPointInPolygon(pointC,polygon) );


       },1000);

因为在代码里把在多边形上的点默认值为boundOrVertex =true 。所以C点也属于范围内的

通过计算点与范围的关系，可以拓展计算范围与范围的交叉重叠计算，在实际业务中很有用喔。