Zwei nächste Punkte auf der Grenze der Postgis-Geometrie

Question

Ich habe eine Tabelle geofences, die geometry von Polygon speichert.

Ich habe auch einen Punkt A, die in der Geometrie ist. Was ich zu tun habe, ist die beiden nächsten Punkte von Punkt A zu finden, die auf der Oberfläche der Polygongeometrie liegen.

Funktion in 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; 

Ich möchte nur zwei Punkte auf der Grenze der Polygongeometrie zu finden. Genau wie ich habe eine von der Frage gefunden:

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; 

Andere dann diesem Punkt Ich habe noch eine finden, mit Abstand größer ist dann der Punkt oben finden, aber kleiner als die übrigen Punkte.

Answer 1

Verwenden Sie ST_DumpPoints(), um die Punkte des Polygons auszugeben, und wählen Sie dann aus dieser Reihenfolge mit ST_Distance bis A limit 2.?

So ist es so etwas wie

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

(davon ausgegangen, dass dies eine innere select ist, wo Poly das Polygon ist, A der Punkt, zu vergleichen und geom ist die geom Spalte einer der Punkte in der Poly wird verglichen)

Answer 2

Im Allgemeinen gibt es keinen zweiten Punkt auf dem Grenzpolygon, wenn Sie die Linien einschließen. Genau wie es keine echte Nummer gibt, die am nächsten an Null ist. Entweder möchten Sie nur die Punkte an den Ecken berücksichtigen, wie Markus es vorschlägt. Oder Sie haben nur einen nächsten Punkt.

Answer 3

1) Art eines Links-Feld-Idee, aber die Sekunde -closest Punkt zu Ihrem Ziel zu finden, warum nicht den nächsten Punkt zu dem Punkt finden, die Sie bereits gefunden?

2) Oder mehr Germaine auf Ihre spezielle Frage,

• die Menge von Punkten innerhalb eines vernünftigen Bereichs des Punktes finden,
• den Schnittpunkt dieses Satzes mit dem Satz von Punkten finden, die auf das Polygon Grenze; ​​

3) Noch weiter in die linke Feld, Dump einige Ihrer Daten-Sets in Mongo (die ich eine andere PostGIS Funktion sein bin zu raten, kann nicht PostG in eine Weile so bin ich nicht sicher verwendet) und Verwenden Sie die Funktion $ near ... http://docs.mongodb.org/manual/reference/operator/near/

Answer 4

Ich gehe davon aus Sie die Kante des Polygons finden möchten, dass die am nächsten an der betreffenden Stelle geht

Um den Abstand ‚d‘ von Punkt ‚C‘ aus der Leitung zu erhalten [A , B] zunächst alle Punkte so übersetzen A bei 0,0

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

Dann normalisieren B ist, so dass es der Länge 1 ist.0

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

Suche Länge von A, die erhalten zu C

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

nun der Entfernung im rechten Winkel

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

nicht sicher, wie das in SQL zu tun. Sie müssen das oben für jede Kante auf Ihrem Vieleck tun, und dann finden Sie den kleinsten Wert 'd'

Übrigens wird Ihnen das Zeichen des Kreuzprodukts von B und C jetzt sagen, ob der Punkt eingeschaltet ist das Innere des Polygons oder nicht