Python 使用光栅和Geopandas显示感兴趣区域的光谱

我相信这是一个非常简单的问题，但出于某种原因，我真的很难找到这样做的方法

我绘制了一幅卫星图像，有13个波段和位于其顶部的形状文件，使用光栅和Geopandas，如下所示：

shapefile = gpd.read_file('ROI.shp')


# plots both elements
fig, ax = plt.subplots(figsize=(20,10))

ax = rasterio.plot.show(src.read([4,3,2]),transform=src.transform,title='Image- bands 4,3,2 ',vmin=0,vmax=30000,with_bounds=True, ax=ax)
shapefile.plot(ax=ax)

这是shapefile的表格：

id  ROI geometry
0   4   Urban   POLYGON ((34.45501 31.53714, 34.48213 31.52649...
1   2   Sea POLYGON ((34.41933 31.54509, 34.40629 31.50923...
2   3   Agriculture MULTIPOLYGON (((34.49357 31.36081, 34.49463 31

我的问题是，我有没有办法描绘出我所创建的感兴趣区域的全景？例如海洋、农业和城市区域？ 我想画一张有光谱的端线图。 我能想到的唯一方法是光栅化形状文件，然后是熊猫并绘制它，我相信这是更好的方法来获得与光栅相同的结果

---

当X轴是图像的波段（或波长），Y轴是该波段的值时，我希望的输出是折线图。

您好！我不太确定你的产量。我理解如下：对于组合的ROI（Urban+…），您希望有一个13条（每个波段一条）的直方图，每个条显示该波段ROI覆盖的所有像素的总和（？）？对吗？嗨！我不太确定你的产量。我理解如下：对于组合的ROI（Urban+…），您希望有一个13条（每个波段一条）的直方图，每个条显示该波段ROI覆盖的所有像素的总和（？）？这是正确的吗？