postgresql - instalar - postgis install

Dos puntos más cercanos en el límite de la geometría de Postgis

Question

Tengo una tabla de geofences que almacena geometry de polígono.

También tengo un punto A que está dentro de la geometría. Lo que tengo que hacer es encontrar los dos puntos más cercanos desde el punto A que se encuentran en la superficie de la geometría del polígono.

Función en 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;

Solo quiero encontrar dos puntos en el límite de la geometría del polígono. Al igual que he encontrado uno de la consulta:

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;

Aparte de este punto, tengo que encontrar uno más con una distancia mayor que el punto que se encuentra arriba pero más pequeño que el resto de puntos.

Answer 1

1) Tipo de idea de campo izquierdo, pero para encontrar el segundo punto más cercano a su destino, ¿por qué no encontrar el punto más cercano al punto que ya encontró?

2) O, más germaine a su pregunta específica,

• encontrar el conjunto de puntos dentro de un rango razonable del punto,
• encuentre la intersección de ese conjunto con el conjunto de puntos que se encuentran en el borde del polígono (que supongo que puede ser otra función de PostGIS; no he usado postG desde hace tiempo, así que no estoy seguro)

3) Más lejos en el campo izquierdo, descargue parte de su conjunto de datos en Mongo y use la función $ near ... http://docs.mongodb.org/manual/reference/operator/near/

Answer 2

Generalmente, no hay un segundo punto más cercano en el polígono de límites, si incluye las líneas. Al igual que no hay un número real segundo más cercano a cero. O solo desea considerar los puntos en las esquinas, como sugiere Markus. O solo tienes un punto más cercano.

Answer 3

Supongo que desea encontrar el borde del polígono que pasa más cerca del punto en cuestión

Para obtener la distancia ''d'' del punto ''C'' de la línea [A, B] Primero traduzca todos los puntos, de modo que A esté en 0,0

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

Luego normaliza B por lo que es de longitud 1.0

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

Encuentra la longitud de A que está en ángulo recto a C

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

Ahora toma la distancia

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

No estoy seguro de cómo hacer eso en SQL. Deberá hacer lo anterior para cada borde de su polígono, y luego encontrar el valor más pequeño ''d''

Por cierto, el signo del producto cruzado de B y C ahora te dirá si el punto está dentro del polígono o no.

Answer 4

Use ST_DumpPoints () para volcar los puntos del polígono, luego seleccione de esa orden por ST_Distance a un límite 2.?

Así que es algo como

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

(asume que esta es una selección interna donde poli es el polígono, A es el punto para comparar y geom es la columna de geom de uno de los puntos en el que se está comparando poli)