分离数据帧纬度/纬度对，并根据列值绘制多个数字

Question

我有一个形成多边形的纬度/纬度对的数据框。纬度/纬度列中带有"GAP"和 NaN 的行是分隔符。所以在这种情况下，我有 4 个多边形，具有多个纬度/纬度位置。我的目标是将这些多边形彼此分开，然后使用卡托皮绘制。

0      1           2
0   POINT  87.6298     9.397332
1   POINT  87.8435     9.842206
2   POINT  87.2354     9.472004
4     GAP         NaN       NaN
5   POINT  87.8354     9.397332
6   POINT  87.9544     9.472004
7   POINT  87.9632     9.191509
8   POINT  87.6244     9.221509
9   POINT  87.4554     9.397332
10    GAP         NaN       NaN
11  POINT  87.6249     9.397332
12  POINT  87.7556     9.221509
13  POINT  87.5567     9.086767
14  POINT  87.3222     9.397332
15    GAP         NaN       NaN
16  POINT  87.6554     9.221509
17  POINT  87.9667     9.191509
18  POINT  87.8854     9.056767
19  POINT  87.4452     9.086767

假设只要运行此命令，每个多边形中的多边形数量和纬度/纬度对的数量就会发生变化。

设置的示例代码如下：

df = pd.read_excel(xl, sheet_name=0, header=None)
#change column names
df.rename(columns={1:'lon', 2:'lat'},inplace=True)
#replace GAP with nan so i can separate by group numbers with each line of nans (didnt work with 'GAP')
df.replace('GAP',np.nan, inplace=True)
df['group_no'] = df.isnull().all(axis=1).cumsum()
#define amount of unique numbers and put into list for looping
numbers = df['group_no'].unique()
aa = list(numbers)

这是我迷路的地方(设置后和绘制代码之前的区域(，如下所示。

a_lon, a_lat = 87.8544, 8.721576
b_lon, b_lat = 87.6554, 8.585951

fig,ax = plt.subplots()
plt.figure(figsize=(10.6,6))
proj = ccrs.PlateCarree()
ax = plt.axes(projection=proj)
ax.stock_img()
ax.set_extent([90, 85, 7, 11], crs=ccrs.PlateCarree()) 

proj = ccrs.Geodetic()
plt.plot([a_lon, b_lon], [a_lat, b_lat], linewidth=1, color='blue', transform=proj)
#plt.show()

如您所见，我用 NaN 替换了"GAP"，然后用"group_no"分隔行，这是一个新列。然后把楠的拿走了。生成的数据帧：

0         lon       lat  group_no
0   POINT  87.6298     9.397332         0
1   POINT  87.8435     9.842206         0
2   POINT  87.2354     9.472004         0
4   POINT  87.8354     9.397332         1
5   POINT  87.9544     9.472004         1
6   POINT  87.9632     9.191509         1
7   POINT  87.6244     9.221509         1
8   POINT  87.4554     9.397332         1
10  POINT  87.6249     9.397332         2
11  POINT  87.7556     9.221509         2
12  POINT  87.5567     9.086767         2
13  POINT  87.3222     9.397332         2
15  POINT  87.6554     9.221509         3
16  POINT  87.9667     9.191509         3
17  POINT  87.8854     9.056767         3
18  POINT  87.4452     9.086767         3

我已经尝试了一些事情，但似乎无法完成交易。我尝试使用group_by在字典中将它们分开，键是group_no，值是纬度/纬度对，但我对字典的了解不足以操纵它们并为每个"键"绘制图。

我还尝试将每个数据帧分离到一个新的数据帧中。以为我可以循环并创建一个 df0、df1 等，然后使用 for 循环进行绘图，但也无法弄清楚。

任何帮助将不胜感激，如果需要更多详细信息，请询问。

Answer 1

你快到了，如果你在组号上调用 groupby，你可以拉出每个组并获取纬度/纬度对。 确保您也设置了适当的投影。

from shapely.geometry import Polygon
for group_no,group_data in df.groupby('group_no'):
poly_coords = group_data[['lon','lat']].values
# Whatever function you are using to create shape with the 'poly_coords e.g.'
polygon = Polygon(poly_coords)
#add to map ...