当计算多边形的面积时,使用shapely/pyproj的结果与谷歌地球引擎或geojson.io相比有很大的不同。我认为这与投影有关,当我研究它时,我发现了一些示例,例如(How do I get the area of a GeoJSON polygon with Python),它们建议使用阿尔德斯等面积投影,我也想知道谷歌地球引擎和geojson.io是否有一种自动的方式来选择区域适当的区域投影,给定多边形</
我用python和geojson来做这件事,我想指定一个点,这个点将是一个正方形的中心,假设平方是1英里×1英里,我想列出所有在广场上找到的点和圆点,包括大于正方形的圆点。我有多个geojson文件,所以需要检查几次,这很好。我一直在使用下面的代码来检查细胞中心是否靠近正方形的中心,但是会有奇怪形状的多边形的问题。我真的很想知道在广场上找到的所有物品/功能。import json
from shapely.geometry import sh
现在我想有效地确定哪些多边形是孔,并减去孔,即减去完全位于另一个多边形内部的较小的多边形,并返回一个“清理”多边形的列表。一对洞和父多边形应该像这样变换(所以基本上是从父多边形中减去洞):
在Stackoverflow和gis.stackexchange.com上有很多类似的问题,但我还没有找到一个真正解决这个问题的from shapely.geometry import Point
from