<问题已解决>求助，Grads提取nc形成dat后图形与原图不一致

xiaoting

因要提取nc的部分数据进行处理，先前看到“http://bbs.06climate.com/forum.p ... ge%3D11&page=1grads”和“http://bbs.06climate.com/forum.p ... 670&extra=page%3D21”帖子中提到提取后的数据会有一些错误，所以在进行数据处理前先进行了下检验，果真图形不一样，麻烦高手指点下，看下gs和ctl哪儿出问题了，谢谢！
提取nc的gs:
'reinit'
'set gxout fwrite'
'set fwrite e:/nc/hgt.200801.dat'
'sdfopen e:/nc/hgt.2008.nc'
'set lon 50 160'
'set lat 10 70'
'set lev 700'
t=1
while(t<=31)
'set t 't
'd hgt'
t=t+1
endwhile
'close 1'
'disable fwrite'
'reinit'

dat配套的ctl:
DSET E:/nc/hgt.200801.dat
TITLE Upper Data
undef -9.99E33
xdef 44 linear 50 2.5
ydef 24 linear 10 2.5
zdef 1 levels 700
tdef 31 linear 01jan2008 1dy
vars 1
h 1 99 geo.height.
Endvars

正确的图

转换成.dat后再画出来的图

mofangbao

把你原来的nc文件 q ctlinfo 看看出来的信息，看图感觉是你的经纬度设置有问题，还有转换的时候建议换算成 set x  set y的形式，因为你的经纬度范围不一定刚好落在格点上，所以你的转换后的ctl应该根据实际的格点数来写，用 set x set y的形式就能很方便的看出来

叩容

你的问题和他的不一样。
他的经纬度没错，你的这个看样子是经纬度有问题。。。。

叩容

孩子，你的grads课肯定没好好听。
真是的。自己去检查下ctl！
x的点数应该是(160-50)/2.5+1，你少了一个点。y同理。

修远

整个问题很明显，你用grads提取数据是set lon 50 160 因为是隔2.5取一个，所以取了45个（不信自己数数），set lat 也是一样。因而需要修改ctl，x 45；y25；即可

xiaoting

叩容 发表于 2012-5-3 22:19

                               
登录/注册后可看大图

孩子，你的grads课肯定没好好听。
真是的。自己去检查下ctl！
x的点数应该是(160-50)/2.5+1，你少了一个点 ...

恩，谢谢，问题解决啦，我的grads是在这个论坛学的，呵呵，有很多基础知识还需一点点完善

xiaoting

修远 发表于 2012-5-3 22:27

                               
登录/注册后可看大图

整个问题很明显，你用grads提取数据是set lon 50 160 因为是隔2.5取一个，所以取了45个（不信自己数数），s ...

谢谢，问题已圆满解决

-逆光年年-

想请问下，我想画省级或市级的图，分辨率应该怎样设置呢？是1*1度吗？谢谢