python - 如何有效地仅绘制那些落在 map 边界框中的点？

Question

我有以下代码，其中 mp 是 Fiona 框的 Shapely Multipolygon:

ax = fig.add_subplot(111)
minx, miny, maxx, maxy = mp.bounds
w, h = maxx - minx, maxy - miny
ax.set_xlim(minx, maxx)
ax.set_ylim(miny, maxy)
ax.set_aspect(1)
patches = []
for idx, p in enumerate(mp):
    patches.append(PolygonPatch(p, fc='#ffffff', ec='#000000', alpha=1., zorder=1))
ax.add_collection(PatchCollection(patches, match_original=True))
plt.show()

我有一个 pandas 数据框，其中包含许多带标签的纬度/经度对。

> 0 -34.026810  -71.101670  A
> 1 -32.986970  -71.501580  G
> 2 -40.931400  -73.027190  C

但是，并非 df 中的所有点都“落入” map 边界内。问题是:如何有效地仅映射 map 边界框中的那些点？

Answer 1

我想在 pandas 级别上过滤点是最有效的。为此，使用其 bounds 属性获取 Multipolygon 的边界，然后使用

过滤数据帧

df[(df.lat > latmin) & (df.lat < latmax) 
   & (df.long > longmin) & (df.long < longmax)]

我假设您的纬度列名为 lat，经度列名为 long。 latmin、latmax 等是从多重多边形边界中获得的值。