python - 跨多个地理围栏的地理围栏半径分布/比率

Question

我知道这个问题可能是简单的几何问题，但我正在为一个地理区域使用特定的 shapefile 并寻找特定的 gis 解决方案(希望是代码)。

本质上，我想找到覆盖每个独立地理边界的黑色区域。覆盖红色区域、蓝色区域和绿色区域的百分比是多少？

我有中心的坐标，以及形状文件 (.shp)。半径已知。除了开源(R 和 Python)之外，对库没有偏好。

Answer 1

如果您可以考虑这样的选择，则可以将这些几何图形加载到 PostGIS 数据库中。之后，您可以执行简单的 SQL 查询，使用内置函数 st_buffier() 和 st_intersection() 来获得所需的输出。否则，您可以使用 GDAL 库的 Python 实现。