雷达数据绘制经纬度图

channam

请问  雷达数据画这种经纬度平面图  仰角是怎么处理的呢？ 这两张图片都没有说明是不是画的某一个仰角的图。

第一张图是来自这个帖子 http://bbs.06climate.com/forum.p ... 59768&extra=&page=1
里面有完整的代码

从他的第二段程序来看，x、y坐标是二维坐标，第一维代表了仰角，插值的时候是怎么操作呢？
这一步没太看明白，rf,x,y = griddata((x,y), rf, method='nearest')，等号左边rf，x，y是什么意思？是对这三个要素都插值吗？
从网上查到这个函数griddata(points, values, xi, method='linear', fill_value=numpy.nan, rescale=False)   xi是需要插值的空间，一般用 numpy.mgrid 函数生成后传入；而这里是没有设置 xi 的 。

小其其格

scan = 1
rf = f['Reflectivity'][scan,:,:]

这里scan就是仰角序号，1是第一个仰角

channam

小其其格 发表于 2022-1-12 21:19
scan = 1
rf = f['Reflectivity'][scan,:,:]

我以为scan是从第一行开始读的意思，按我理解是这句是提取仰角值的   ele = f['elevationR'][scan,:]  

下面这几句，j 是仰角的总个数，循环了每个仰角，x和y的第一维代表仰角，第二维代表距离中心点的距离
for j in xrange (len(e)):
    x[j,:] = dis * cos(e[j]) * cos(a[j])
    y[j,:] = dis * cos(e[j]) * sin(a[j])

就是后面插值的时候不知道仰角是怎么处理的

请问您之前处理过雷达数据吗？我是想处理CB型的雷达数据，这个程序是处理SA型数据的，像这些'Reflectivity'、'elevationR'都是SA数据自带的列标签吗？