metpy - 将图像保存为 geotiff（或 png）

Question

我正在使用metpy 和siphon 来下拉和显示雷达数据。理想情况下，我只想保存雷达层以在 LeafletJS 中显示。

GDAL 似乎提供了一个使用 gdal_translate 的解决方案：

gdal_translate -of PNG -ot Byte -b 1 -scale 10 100 1 255 mygribfile.grib2 out.png

不过，如果我已经通过metpy和cartopy库创建了数据层，我假设有一种方法可以只获取雷达层，包括相关的颜色表，并以地理参考的方式将其保存为图像这样我就可以将它作为一个图层包含在其他地方。（我认为 matplotlib.pyplot.savefig() 不会保留地理信息。）

是否有人可以向我指出一些可以实现这些方法的示例代码？我的谷歌搜索（到目前为止）是空的。

谢谢！

Answer 1

matplotlib 或metpy 中没有任何东西可以本地处理这个问题。您最好的选择可能是尝试使用rasterio创建 geotiff，这将允许您指定地理位置信息。您应该能够将 rasterio 传递给与渲染图形相对应的 numpy 数组，如本示例中所做的那样。