关于绘制等经纬度netcdf格式雷达图回波区位置偏移的问题

樱桃小丸子

最近，正在应领导需求，将一个二进制雷达数据产品转换成nc格式，并绘图。在顺利转换格式后，绘图时发现python基于basemap绘图时，雷达数据稍有偏移。如下图1：

图1

，图2

图2

是雷达处理程序自动画的图，应该是grads画的，虽然colorbar有点不一样，但是明显可以看到，红笔圈出的地方的存在一定的位置偏移，用雷达中心经纬度对比了一下，图1的经纬度确实是错的，我开始以为是投影的问题，换了几个不同的投影方式，都不对，后来仔细研究一下，查了很多资料，但是basemap关于makegrid的资料非常少，我尝试解释一下，不过理解的不一定对，雷达数据程序生成的格点数据是等经纬度的数据，但是basemap中makegrid的数据应该是等距的，所以应该先将经纬度坐标转换成投影坐标，之后位置就对上了，如图3

图3

。我之前的代码是在家园下载的其他大神的，如图4

图4

，但是用同样方法画出来的就不对，可能是源代码是读wrfout数据的，本身数据就已经是投影之后的坐标，所以不需要转换，因此把代码改为图5，

图5

。

Masterpiece

basemap.makegrid按照文档说明是
return arrays of shape (ny,nx) containing lon,lat coordinates of an equally spaced native projection grid.
的确和你猜测的一样，返回的是一个等距离网格的经纬度。如果等经纬度网格的数据，用这个等距离网格的坐标去画图，那南北方向就会出现位置的偏移。

一般我们都不用这个makegrid，都是用一维的lon和lat，直接numpy.meshgrid结合成二维，然后再用basemap投影一次绘图坐标

樱桃小丸子

Masterpiece 发表于 2020-4-21 13:37
basemap.makegrid按照文档说明是
return arrays of shape (ny,nx) containing lon,lat coordinates of an  ...

恩，是的，之前也是在群里下载的其他大神画图的代码，在他的基础上只是改了路径和一些变量名而已，没想到出现了这样的问题，看来还是需要自己多钻研，不能投机取巧。感谢解答。

ququ

请问下能共享下代码吗？或者源代码的地址？