我正在使用metpy和siphon下拉并显示雷达数据。理想情况下，我想仅保存雷达层以显示在（例如）中。
GDAL似乎使用gdal_translate提供解决方案：
gdal_translate -of PNG -ot Byte -b 1 -scale 10 100 1 255 mygribfile.grib2 out.png
不过，如果我已经通过Metpy和Cartopy库创建了数据层，我假设只有一种仅摄取雷达层的方法，包括关联的colortables，并将其作为图像保存在地理上 - 引用的方式，以便我可以将其包含在其他地方。（我认为matplotlib.pyplot.savefig（）不会维护地理信息。）
是否有人可以将我指向一些可以实现这些手段的示例代码？到目前为止，我的谷歌搜索已经空了。
谢谢！
Matplotlib或Metpy中没有任何东西可以在本地处理此操作。最好的选择可能是尝试使用rasterio创建Geotiff，这将使您可以指定地理位置信息。您应该能够像在此示例中一样，将rasterio射击与渲染图相对应的numpy数组。