Postgis 空间查询

Question

我正在尝试创建一个查询，该查询将返回不与线相交的点列表。在大多数情况下，它工作得很好而且非常快。不过我想知道最近的线路有多远。

我尝试使用 st_distance 函数，但它返回 null。有没有人有任何提示或指导来帮助我解决这个问题？

这是我目前的 SQL 语句：

SELECT p.objectid, p.voidentifier, p.featuretypecd, p.heightagl, st_distance(p.geom::geography, l.geom::geography), p.geom 
FROM vo_point p 
  LEFT JOIN vo_line l ON ST_Intersects(p.geom, l.geom) 
WHERE l.featureid IS NULL 
  and p.featuretypecd in ('537','538','540','541','542','543','544');

Answer 1

我相信您正在寻找子查询：

SELECT p.objectid, p.voidentifier, p.featuretypecd, p.heightagl, p.geom,
  (SELECT min(ST_Distance(p.geom::geography, l2.geom::geography)) FROM vo_line l2)
FROM vo_point p 
LEFT JOIN vo_line l ON ST_Intersects(p.geom, l.geom) 
WHERE l.featureid IS NULL 
  AND p.featuretypecd IN ('537','538','540','541','542','543','544');

这个查询应该返回点和线相交的所有记录，它还会返回到最近线的距离。

Answer 2

所以一个多小时后，我决定中止并尝试不同的方法。 相反，我创建了两个视图：

create view pylon_temp as SELECT p.objectid, p.voidentifier, p.featuretypecd, p.heightagl, p.geom
FROM vo_point p LEFT JOIN vo_line l ON ST_Intersects(p.geom, l.geom)
WHERE l.featureid IS NULL and p.featuretypecd in ('537','538','540','541','542','543','544');

create view pylon_errors as
select p.objectid, p.voidentifier, p.featuretypecd, p.heightagl, min(st_distance(p.geom::geography, l.geom::geography))as distance, p.geom
from pylon_temp p, vo_line l
where st_dwithin(p.geom,l.geom,0.001)
group by p.objectid, p.voidentifier, p.featuretypecd, p.heightagl, p.geom
order by distance desc;

这可能不是最优雅的方式，但它非常快，我可以在我的网络应用程序/GIS 应用程序中使用视图。