我有一个像这样的Pandas数据帧：

 id   lat  lon
 jhg  2.7  3.5
 ytr  3.1  3.5
 ...

我还有一个带有一些多边形的Geopandas数据帧多边形.现在,我想只绘制一些多边形内的df点.所以我应该能够做像poly.intersects(p)这样的事情,其中​​p是一个Shapely Point.但我做错了;

from shapely.geometry import Point
for index, row in df.iterrows():
    t = poly.intersects(Point(row.lon, row.lat))

传递具有lat / lon点的数据帧并将其覆盖到poly的最佳方法是什么？请注意,我可以定义最小/最大纬度/经度范围,但也可以在聚合物外部打印点,但在(较大的)边界框内打印.

最佳答案 你的出发点：

import pandas as pd
from shapely.geometry import box
import matplotlib.pyplot as plt

from matplotlib.collections import PatchCollection
from matplotlib.patches import Polygon
from shapely.geometry import Point
import seaborn as sns
import numpy as np

# some pretend data
data = {'lat':[2.7,3.5,1.4,2.3,.9,1.9], 'lon':[1.2,.9,1.9,2.2,3,1.1]}
df = pd.DataFrame(data)

# the 'bounding' polygon
poly = box(1,1,2,2)
patches  = PatchCollection([Polygon(poly.exterior)], facecolor='red', linewidth=.5, alpha=.5)


# plot the bounding box 
fig, ax = sns.plt.subplots(1, figsize=(4,4))
ax.add_collection(patches, autolim=True)

# plot the lat/lon points
df.plot(x='lat',y='lon', kind='scatter',ax=ax)
plt.show()

这些数字看起来像这样：

摆脱不需要的点的一种方法是使用布尔掩码：

#probably more efficient ways to do this, but this works
mask = [poly.intersects(Point(lat,lon)) for lat,lon in zip(df.lat,df.lon)]
df = df[mask]

# make new plot (must make a new 'patch' object)
patches1  = PatchCollection([Polygon(poly.exterior)], facecolor='red', linewidth=.5, alpha=.5)
fig1, ax1 = sns.plt.subplots(1, figsize=(4,4))
ax1.add_collection(patches1, autolim=True)

# make the axis bounds the same
ax1.set_xlim(ax.get_xlim())
ax1.set_ylim(ax.get_ylim())

# plot the lat/lon points
df.plot(x='lat',y='lon', kind='scatter',ax=ax1)
plt.show()

给我这个形象.

请注意,您可以使用其他更快的方式创建布尔掩码,例如lat是否位于多边形中的最高点之上.那些可能不是完美的,但可以减少问题,所以你没有多次调用intersects().

[编辑：如果您的多边形是一个矩形,]另一种方式(正如您在问题中所建议的那样)将只是围绕边界多边形“裁剪”图像.这是一个更快的解决方案,因为您不必一遍又一遍地调用intersects()函数.要根据边界多边形修剪图像,可以在plt.plot()之前插入此权限：

ax.set_xlim((np.min(poly.exterior.xy[0]),np.max(poly.exterior.xy[0])) )
ax.set_ylim((np.min(poly.exterior.xy[1]),np.max(poly.exterior.xy[1])) )

给出以下内容：