grads处理nc雷达数据，谢谢

river

lxe3399 发表于 2021-4-12 16:48
谢谢尝试！非常感谢
的确很困惑

问题解决了，这个一定要注意数据格式，不同的数据格式在写描述文件时有些细微的差别。

dset G:/test/1/1606290700.nc
title radar
dtype netcdf
undef  9.96921e+36
unpack scale_factor add_offset 0.5 -33
*options yrev
xdef 1019 linear 111.69 0.0099999998
ydef 881 linear 35.348999 0.0099999998
zdef 20 levels 500 1000 1500 2000 2500 3000 3500 4000
4500 5000 5500 6000 7000 8000 9000 10000 12000 14000
17000 19000
tdef 1 linear 00Z01JAN0001 1mn
vars 1
reflectivity=>reflectivity 20 z,y,x Reflectivity   *这句里面“=>”前面的名字必须严格和文件里的标变量名一致，后面就可以自己定义，比如refl
endvars

lxe3399

river 发表于 2021-4-13 01:02
问题解决了，这个一定要注意数据格式，不同的数据格式在写描述文件时有些细微的差别。

dset G:/test/1 ...

依旧谢谢大神的尝试！
这样写ctl 画出的图经纬度是没有偏移，可是dbz值依然不对呢，读出来的dbz值为非编码值，您注意到了吗？谢谢！


所以
还是回到的原点，即：
读对了编码值  经纬度有偏移，
经纬度无偏移，编码值没有读对
是不是很崩溃/哭

river

lxe3399 发表于 2021-4-13 14:13
依旧谢谢大神的尝试！
这样写ctl 画出的图经纬度是没有偏移，可是dbz值依然不对呢，读出来的dbz值为非编 ...

你指的编码值是0-255吗？我觉得是你误会了，资料详情里面提示说有效值为 0S-255S，所以不在这个范围的值就是无效值，你在画图的时候限制一下数值，不画无效值就可以了啊。

dset G:/test/1/1606290700.nc
title radar
dtype netcdf
undef  9.96921e+36
unpack scale_factor add_offset 0.5f -33s
*options yrev
xdef 1019 linear 111.69 0.0099999998
ydef 881 linear 35.348999 0.0099999998
zdef 20 levels 500 1000 1500 2000 2500 3000 3500 4000
4500 5000 5500 6000 7000 8000 9000 10000 12000 14000
17000 19000
tdef 1 linear 00Z01JAN0001 1mn
vars 1
reflectivity 20 z,y,x Reflectivity
endvars


看下我画的，你再研究一下，和实况对比看看对不对。

radar2

画图的脚本
'reinit'
'open g:/test/1/radar.ctl'
'set grads off'
'set mpdset cnhimap'
'set z 8'
'set cmin 0'
'set gxout shade2b'
'g:/test/1/15colors.gs'
'set clevs   0  5 10 15 20 25 30 35 40 45 50 55 60 65 70'
'set ccols  -1 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34'
'd reflectivity'
'cbar'
'gxprint g:/test/1/radar.png white'
;


  

lxe3399

river 发表于 2021-4-14 00:52
你指的编码值是0-255吗？我觉得是你误会了，资料详情里面提示说有效值为 0S-255S，所以不在这个范围的值 ...

再次感谢您的尝试！
是的编码值是0-255，读出来这个编码值之后，再除以2减去33即为正确的雷达回波dbz值。（要读出正确的编码值，我感觉好像ctl里vars那里  一定得写-1，40，1而不是写成z，y，x，（不确定，供参考））
您画的图dbz值与实况相比  强回波没有出来
我发帖的一楼第二个图即为同时次的实况图，这是一次强对流天气过程，此时次北京的西北和西南部均被强回波覆盖。
您可以缩小区域画一下北京周边，就能看到您图中北京西北和西南地区强回波是空白。或者不人为设dbz值，让它自动出图会看的很明显。

lxe3399

river 发表于 2021-4-14 23:44
我修改了你的描述文件，增加了fileheader项，基本能画出你要的效果，和原图在纬度上稍微差一点点。filehe ...

好的  多谢您的尝试！
我再看一看
非常感谢！

lxe3399

请问大神，为何设置为问题已解决，有了最佳答案之后，您这条回复就看不到了呢？谢谢！