Postgis幾何邊界上的兩個最近點

Question

我有一個表geofences，它存儲了多邊形的geometry。

我也有一個點A這是幾何內。我需要做的是找到位於多邊形幾何體表面上的點A的兩個最近點。

功能在PostGIS的：

CREATE OR REPLACE FUNCTION accuracyCheck(Polygon geometry 
             ,decimal lat 
             ,decimal lon) 
    RETURNS VARCHAR AS 
$BODY$ 

DECLARE height DECIMAL; 
DECLARE accuracy VARCHAR(250); 

BEGIN 

CREATE TEMPORARY TABLE closePointStorage AS 
SELECT ST_AsText(ST_ClosestPoint(geometry 
           ,ST_GeomFromText('POINT(lat lon)',0) 
           ) 
       ) AS closestPoint 
FROM (
    SELECT ST_GeomFromText(geometry) as geometry 
    FROM gfe_geofences 
    WHERE is_active=true 
    ) As tempName; 

CREATE TEMPORARY TABLE areaStorage ON COMMIT DROP AS 
SELECT ST_Area(ST_GeomFromText('Polygon((23.0808622876029 96.1304006624291 
             ,28.0808622876029 99.1304006624291 
             ,100    200 
             ,23.0808622876029 96.1304006624291 
             ))' 
           ,0) 
      ) AS area; 

CREATE TEMPORARY TABLE distanceStorage ON COMMIT DROP AS 
SELECT ST_Distance(
      ST_GeomFromText('POINT(23.0808622876029 96.1304006624291)',-1) 
     ,ST_GeomFromText('POINT(28.0808622876029 99.1304006624291)',-1) 
     ) AS distance; 

height = (SELECT area FROM areaStorage) 
     /(0.5*(SELECT distance FROM distanceStorage)); 

IF height < (SELECT radius_meters 
      FROM gfe_geofences Where is_active=true) THEN 
    accuracy = "FullConfirm"; 
    RETURN accuracy; 
ELSE 
    accuracy = "PartiallyConfirm"; 
    RETURN accuracy; 
END IF; 

END; 
$BODY$ LANGUAGE plpgsql; 

我只是想找到多邊形形狀的邊界兩分。就像我已經找到一個從查詢：

CREATE TEMPORARY TABLE closePointStorage AS 
SELECT ST_AsText(ST_ClosestPoint(geometry 
           ,ST_GeomFromText('POINT(lat lon)',0) 
           ) 
       ) AS closestPoint 
FROM (
    SELECT ST_GeomFromText(geometry) as geometry 
    FROM gfe_geofences 
    WHERE is_active=true 
    ) 
AS tempName; 

其他然後這一點我必須找到一個更隨距離大於點發現以上，但隨後小點的休息。

Answer 1

使用ST_DumpPoints（）轉儲多邊形的點，然後通過ST_Distance選擇該順序爲限制2。

因此，它是像

SELECT * from ST_DumpPoints(poly) order by ST_Distance(A,geom) asc limit 2; 

（假定這是一個內部的選擇，其中多是多邊形，是比較點和GEOM是在該點之一的GEOM列poly被比較）

Answer 2

如果包含線，那麼通常在邊界多邊形上沒有第二近點。就像沒有第二個接近於零的實數一樣。 要麼你只是想考慮在角落的點，就像馬庫斯建議的那樣。 或者你只有一個最近的點。

Answer 3

1）一種左場的想法，但找到第二個 - 最接近您的目的地，爲什麼不找到最接近你已經找到的點？

2）或者，更傑曼到您的具體問題，

• 找到一點有一些合理的範圍之內的點的集合，
• 發現集合的交集與設定躺在點多邊形邊界（我猜可能是另一個PostGIS的功能;沒有使用postG在一段時間，所以我不知道）

3）到更遠的左場，傾倒你的一些數據集到蒙戈和使用$ near功能... http://docs.mongodb.org/manual/reference/operator/near/

Answer 4

我假設你想找到一個通過了問題

要獲得線點「C」的距離「d」 [A最接近點的多邊形的邊，B] 首先翻譯所有點使A爲0,0

B -= A //vector subtraction 
C -= A 

然後正常化乙所以它是長度爲1。從A 0

len = sqrt(B . B) //dotproduct of two vectors is the length squared 
B /= len //scalar divide by length 

查找長度成直角到C

dotp = B . C //dot product again 
closestPointOnLine = B * dotp //scalar multiply 

現在得到的距離

diff = (C - ClosestPointOnLine) 
d = sqrt(diff . diff) 

不知道如何做，在SQL。你需要爲你的多邊形上的每個邊做上述操作，然後找到最小值'd'。

順便說一句，B和C的叉積的符號現在會告訴你該點是否開啓多邊形的內部或不是