python - PostGIS ST_intersects python 等价物

Question

我目前使用 PostGIS 作为我在 python 脚本中执行的许多空间函数的主干。特别是获取几个 shapefile 几何形状并查看它们是否相交，然后将它们分类到单独的目录中。我使用 shp2pgsql 上传 shapefile，然后使用 ST_Intersects 将它们关联起来，然后使用脚本中的 os/shutil 函数对它们进行排序。

我的问题是我们的一个团队只在政府网络上工作，无法让他们的系统管理员批准 postgres/postGIS。是否有一个 python 模块/函数可以在不需要 postgres 的情况下执行与 ST_Intersects 相同的几何相关性？或者，如果我需要自己编写，是否有与几何有关的算法的站点。例如，如果我有一个左上角和一个右下角坐标，我该如何计算其他两点。我不要求任何人为我编写代码，只是一些帮助指向正确的方向。

还有 WGS 1984 中执行的所有数据

Answer 1

有许多读取 Shapefile 的工具，您可以使用它们来获取它们的范围或界限。这些可用于使用 Rtree 构建 R 树索引pacakge，其中有一些good examples in the documentation .使用 R 树索引，您可以使用 intersection 查看边界框相交的位置。这类似于 PostGIS 的 GiST 索引，但根据我的经验，它的构建和使用速度要快得多。如果/当您需要对几何进行详细交集时，您可以使用 Shapely，后者又使用 GEOS，它与 PostGIS 使用的库相同。所以它们都以相似的方式相关。

查看这些相关问题:

• Looking for a fast way to find the polygon a point belongs to using Shapely
• Faster way of polygon intersection with shapely