我对 python 很熟悉，但对地理空间编码很陌生。
我有一个坐标列表。我也是人口普查区块组的 shapefile。
我想知道如何找到每个坐标 1 英里内的人口普查区块组。即使只有一小部分块组在 1 英里范围内，我仍然希望该组出现在结果中。我还希望附近块组的集合与它们附近的坐标相关联。
如果有任何不清楚的地方，请告诉我，我会尽量解释得更清楚。
谢谢！
PS我认为这个问题正在解决我正在努力解决的一些问题（https://gis.stackexchange.com/questions/11987/polygon-overlay-with-shapely）。这意味着我应该使用 Shapely，为我拥有的每组坐标创建一个缓冲区，然后以某种方式找到哪些人口普查区域被联合起来。