查找与点相交的所有区域，反之亦然 - PostGIS

Question

我正在开发一个使用 PostGIS 来保存空间数据的项目，其中有保存几何点数据的记录和保存几何区域数据的记录。为了解决我的问题，我正在寻找两个可以采用地理形状而不是几何形状作为参数的查询。

对于查询 AI 需要它返回与给定区域相交的所有点。

对于查询 BI 需要它返回与给定点相交的所有区域。

Answer 1

无论您走向哪个方向，查询都几乎相同，它仅取决于您的选择标准，假设您有下表：

---------------------------------------------------
| id | name  | otherdata | geometry(MULTIPOLYGON) |
---------------------------------------------------
| 1  | poly1 | blah      | 1020304050......       |
---------------------------------------------------

使用以下内容：

SELECT * FROM poly_table

或

SELECT id,name,otherdata,AsText(geometry)

将按原样显示该表（后者将提供更英语可读的输出）

现在，如果您也有一个类似的表，但您有点而不是几何图形中的多边形，要对它们进行交叉过滤，您可以使用简单的连接。

EG：

SELECT points.* FROM points_table points
INNER JOIN polygon_table polys ON ST_Within(points.geometry,polys.geometry)
WHERE polys.id = 1 -- Or some other where clause to restrict the polygon area to the one you want.

这将返回点表中落在代表您感兴趣的区域的所选多边形内的所有点。

注意：我在这里使用 ST_Within，但 PostGIS 中可能还有其他几何运算符可以执行许多其他类型的检查这将使您的相交和/或边界测量更加准确，但对于所有这些，原理都是相同的。

对于您的第二个查询，所需的 GIS-SQL 将非常相似，但您不需要返回给定多边形的点，而是需要反转您的条件以返回给定点的多边形：

SELECT polys.* FROM polygon_table polys
INNER JOIN points_table points ON ST_Within(points.geometry,polys.geometry)
WHERE points.id = 1 -- Or some other where clause to restrict the point to the one you want.

注意：我们将参数保留在 ' ST_Within 的顺序完全相同，因为 ST_Within 使用以下规则进行处理：

returns true IF geometry A is within geometry B

由于该点仍然是两个实体中较小的一个，因此我们仍然需要检查相同的空间顺序。