如何解决使用Rasterio和Geopandas显示感兴趣区域的光谱
我相信这是一个非常简单的问题,但是由于某些原因,我确实很难找到解决方法。
我使用Rasterio和Geopandas绘制了一个卫星图像,该图像具有13个波段和位于其顶部的shapefile,如下所示:
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,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
我的问题是,有什么方法可以绘制我创建的感兴趣区域的光谱?例如海洋,农业和城市区域? 我想获得具有频谱的最终折线图。 我唯一能想到的方法是先将shapefile栅格化,然后将熊猫栅格化并绘制出来,我相信Emsut是与Rasterio取得相同结果的更好方法。
当X轴是图像的波段(或波长)而Y轴是该波段上的值时,我想要的输出是折线图。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。