引用Basemap实例的坐标系

Question

我使用matplotlib底图工具包（mpl_basemap.Basemap）使用pcolor()或pcolormesh()在地图背景上绘制GeoTIFF图像。为此，我需要将栅格数据的x和y坐标的网格网格转换为绘图的坐标系，例如极坐标立体图。目前，我就是这样做的：

• 创建一个osr.SpatialReference()对象，并使用它通过ImportFromWkt()方法获取GeoTiff图像的坐标，并以GeoTIFF的投影作为参数进行调用。
• 通过.ExportToProj4()将此空间参考导出到proj4，并使用它创建pyproj.Proj()对象。
• 在GeoTIFF图像的东向和北向值的网格网格上使用inverse=True调用此类实例以创建纬度/长度网格颜色
• 在此lat / long meshgrid上调用Basemap类实例，以在Basemap对象的坐标系中创建meshgrid

是否有更简单的方法不通过纬度/经度表示？如果我可以获得对Basemap实例的坐标系的引用，我可以直接在GeoTIFF数据集和地图之间设置坐标转换。