Calculando la distancia entre un punto y la línea más cercana, PostGIS

Question

Tengo que agregar una columna a mi tabla de puntos almacenando la distancia (metros) entre cada punto y la línea más cercana. Ambos están en el mismo sistema de coordenadas proyectadas. ¿Cómo puedo hacerlo? He intentado esto pero no funciona:

ALTER TABLE st09_2019_01 ADD coast_dist_M float8;
INSERT INTO st09_2019_01 (coast_dist_M)
VALUES
(SELECT ST_Distance(st09_2019_01.geom, coastline_ln.shape)
FROM st09_2019_01, coastline_ln
ORDER BY ST_Distance(st09_2019_01.geom, coastline_ln.shape)
LIMIT 1)

Answer 1

No puedes INSERTAR en columnas, solo añadir filas a la tabla.

En su lugar, ejecuta un UPDATE:

UPDATE st09_2019_01 AS pt
  SET  coast_dist_m = (
    SELECT ST_Distance(pt.geom, ln.shape)
    FROM   coastline_ln AS ln
    ORDER BY
           pt.geom <-> ln.shape
    LIMIT  1
  )
;

---

Relacionado:

• ¿Cómo hacer un UPDATE con una consulta de vecino más cercano LATERAL?

Información adicional sobre (espacial) (K)NN y el operador <->:

• Forma eficiente de encontrar la característica más cercana entre tablas enormes de Postgres
• Vecino más cercano
• PostGIS: Error de unión ST_ClosestPoint
• Optimizando mi script; tarda mucho tiempo; distancia entre muchos objetos PostGIS

Answer 2

Ambas tablas necesitan un id. El término "&& ST_Expand(p.geom, 10000))" establece el tamaño de la ventana para el cálculo de la distancia (10km). Puedes cambiarlo o quitarlo. Este código calculará la distancia mínima entre cada punto en st09_2019_01 y el punto más cercano en una característica en la capa coastline_ln.

update st09_2019_01 set coast_dist_M = tmp1.distance from (select tmp.pid, min(ST_distance(ST_closestpoint(tmp.r_geom, tmp.p_geom),tmp.p_geom)) as distance from (select p.id pid, r.id rid, p.geom p_geom, r.geom r_geom from st09_2019_01 as p, coastline_ln as r where ST_closestpoint(r.geom,p.geom) && ST_Expand(p.geom, 10000)) as tmp group by pid ) tmp1 where tmp1.pid=st09_2019_01.id;