postgresql - Postgis – 与 ST_Within 相反

Question

我正在尝试获取不在几个多边形内的所有点:

SELECT pt.geom 
FROM points pt, polygons po
WHERE ST_Within(pt.geom, po.geom) = false;

但没有用——我得到了所有的分数。

Answer 1

您看到“所有点”的原因是您有不止一个点和不止一个多边形。如果您只有其中之一，它就可以工作。

其实有几种方法可以解决这个问题，但这可能是最简单、最快的。

如果您有主键或唯一键(例如 gid):

SELECT pt.*
FROM points pt
WHERE pt.gid NOT IN (
    SELECT pt.gid FROM points pt, polygons po
    WHERE ST_Within(pt.geom, po.geom)
);

请注意，它还会返回任何具有 NULL 几何图形的点。

另一种方法是使用EXCEPT:

SELECT pt.*
FROM points pt
EXCEPT SELECT pt.*
FROM points  pt, polygons po
WHERE ST_Within(pt.geom, po.geom);

还有一种方法是使用 ST_Union 将所有多边形融合成一个大几何体，但这种方法会慢得多。