我使用 matplotlib 底图工具包 ( mpl_basemap.Basemap) 在地图背景上使用pcolor()或绘制 GeoTIFF 图像pcolormesh()。为此，我需要将栅格数据的 x 和 y 坐标的网格网格转换为绘图的坐标系，例如极地立体图。目前，这就是我完成此任务的方式：
- 创建一个
osr.SpatialReference()对象并使用它通过其ImportFromWkt()方法获取 GeoTiff 图像的坐标，该方法以 GeoTIFF 的投影作为参数调用。 - 将此空间参考导出到 proj4 via
.ExportToProj4()并使用它来创建pyproj.Proj()对象。 inverse=True在 GeoTIFF 图像的东向和北向值的网格网格上调用此类实例以创建纬度/经度网格网格- 在这个 lat/long meshgrid 上调用 Basemap 类实例以在 Basemap 对象的坐标系中创建一个 meshgrid
有没有更简单的方法不通过纬度/经度表示？如果我可以获得对 Basemap 实例坐标系的引用，我可以直接在 GeoTIFF 数据集和地图之间建立坐标变换。