檢查經緯度是否在諾基亞HERE地圖多邊形中

Question

有沒有方法可以確定經緯度是否在nokia.maps.map.Polygon？

我可以使用方法nokiaMap.getObjectsAt來檢查給定地圖像素位置下是否有多邊形，但這對經度和緯度沒有幫助。

有一種方法nokiaMap.geoToPixel我想用得到的像素位置，然後調用nokiaMap.getObjectsAt，但該文件指出「可見地圖區域外結果值可能是非常不可靠的「這是我見過的情況，所以不能使用這種方法。

nokiaMap.geoToPixel文檔 http://developer.here.com/docs/maps_js/topics_api_pub/nokia.maps.map.Display.html#topic-apiref__Display-geoToPixel-method

Answer 1

這需要從計算幾何的point in polygon算法的實現。這方面的例子可以發現如下：

nokia.maps.map.Polygon.prototype.containsCoordinate = function(arg) { 
    var latitude, longitude; 

    if (arg instanceof nokia.maps.geo.Coordinate){ 
     latitude = arg.latitude; 
     longitude = arg.longitude; 
    } else if (arg instanceof Array && arg.length === 2 && !isNaN(arg[0]) && !isNaN(arg[1])){ 
     latitude = arg[0]; 
     longitude = arg[1]; 
    } 

    // Fail fast if not in the bounding box - don't bother with Ray Cast method 
    if (!poly.getBoundingBox().contains(
     nokia.maps.geo.BoundingBox.coverAll(
     [new nokia.maps.geo.Coordinate(latitude, longitude)]))){ 
     return false; 
    } 

    var inPoly = false, 
     path = this.path, 
     numPoints = path.getLength(), 
     j = numPoints-1; 

     for(var i=0; i < numPoints; i++) { 
     var vertex1 = path.getLatLng(i), 
     vertex2 = path.getLatLng(j); 

     if (vertex1[1] < longitude && vertex2[1] >= longitude || 
      vertex2[1] < longitude && vertex1[1] >= longitude) { 
      if (vertex1[0] + (longitude - vertex1[1])/(vertex2[1] - vertex1[1]) 
      * (vertex2[0] - vertex1[0]) < latitude) { 
      inPoly = !inPoly; 
      } 
     } 
     j = i; 
     } 
    return inPoly; 
    }; 

您可以通過調用它無論是Coordinate或一對數字。 e.g：

poly.containsCoordinate(coordinate) 

Working Example

Answer 2

我寫了一個使用光線投射算法來計算，如果該點在多邊形或不是一個PL/SQL程序。對多邊形的頂點數沒有限制。在這裏，頂點的經度和緯度以及要檢查的點以十進制格式提交，並且是分離的。

create or replace function PIP(verLong in varchar2, 
           verLat in varchar2, 
           xLong in number, 
           yLat in number) return boolean 

is 

    i number := 0; 

    j number := 0; 

    c boolean := false; 

    tempLong number; 

    tempLat number; 

    latDegree number; 

    minLength number; 

    pi NUMBER := 3.1415926; 

    x number := xLong; 

    y number := yLat; 

    TYPE coodArray IS TABLE OF VARCHAR2(1000) INDEX BY BINARY_INTEGER; 

    longArray coodArray; 

    latArray coodArray; 

begin 

    --Convert reference point in meters 

    --Round number to 4 places 

    x := round(x, 4); 

    y := round(y, 4); 

    --convert latitude to degrees 

    latDegree := y/60; 

    --convert latitude from degrees to radians 

    latDegree := (latDegree * pi)/180; 

    --compute minute or arc in meters given the latitude value 

    minLength := 1852.3 - (9.4 * Cos(2 * latDegree)); 

    --Multiply the latitude and longitude pair with the length factor to get co-ordinates in meters 

    x := x * minLength; 

    y := y * minLength; 

    loop 

    --Retrieve the longitudes and latitudes from the string of latitudes or longitudes 

    tempLong := to_number(get_token(verLong, i + 1, '|')); 

    tempLat := to_number(get_token(verLat, i + 1, '|')); 

    i := i + 1; 

    if (tempLong is null) then 

     exit; 

    end if; 

    --Round the numbers upto four decimals 

    tempLong := round(tempLong, 4); 

    tempLat := round(tempLat, 4); 

    --convert latitude to degrees 

    latDegree := tempLat/60; 

    --convert latitude from degrees to radians 

    latDegree := (latDegree * pi)/180; 

    --compute the minute of an arc in meters given the latitude value 

    minLength := 1852.3 - (9.4 * Cos(2 * latDegree)); 

    --Multiply the latitude and longitude pair with the length factor to get co-ordinates in meters 

    tempLong := tempLong * minLength; 

    tempLat := tempLat * minLength; 

    longArray(i) := tempLong; 

    latArray(i) := tempLat; 

    end Loop; 

    j := latArray.COUNT; 

    --Check whether the point lies inside the polygon formed by the vertices (Ray casting algorithm) 

    for i in 1 .. j loop 

    if (((latArray(i) > y) != (latArray(j) > y)) and 
     (x < (Longarray(j) - Longarray(i)) * (y - Latarray(i))/
     (latarray(j) - Latarray(i)) + longarray(i))) then 

     c := not c; 

    end if; 

    j := i; 

    end loop; 

    return c; 

end PIP; 

感謝， Noorul