提取NC数据画出来的图有问题

jeffggg

我提取NC降水数据 然后你配ctl 画出来的图 是这样  我噶虐很奇怪啊 难道是缺测值出问题了？我gs ctl 里缺测值都是q ctlinfo NC文件显示的缺测值啊 怎么会有问题呢？

善人/jw

是很奇怪，但LZ你就贴个图出来，别人想帮也帮不了你~~

jeffggg

善人/jw 发表于 2013-4-29 16:07

                               
登录/注册后可看大图

是很奇怪，但LZ你就贴个图出来，别人想帮也帮不了你~~

这就是你个很简单的 提取 把NC资料的数据提取出来生成grd 我自己配CTL 仅此而已

dzxconan

应该gs写出来的格点有问题或者经纬度问提吧？

jeffggg

dzxconan 发表于 2013-4-29 17:02

                               
登录/注册后可看大图

应该gs写出来的格点有问题或者经纬度问提吧？

我很注意这个了！我仔细检查过了 个点数目和NC里一模一样

dzxconan

jeffggg 发表于 2013-4-29 17:06

                               
登录/注册后可看大图

我很注意这个了！我仔细检查过了 个点数目和NC里一模一样

缺测呢 几个维数都set了吗？

kongfeng0824

可以把你的ctl内容贴出来看看。

jeffggg

kongfeng0824 发表于 2013-4-29 21:10

                               
登录/注册后可看大图

可以把你的ctl内容贴出来看看。

'reinit'
'sdfopen I:\CMAP\cmap\3\cmap\precip.pentad.mean.nc'
'set x  1 144'
'set y  1 72'
'set lev 0 '
'set gxout fwrite'
'set fwrite I:\CMAP\cmap\3\cmap\p2338.grd'
p=1
while(p<=2338)
'set t 'p''
'd precip'
p=p+1
endwhile
'disable fwrite'
;

jeffggg

kongfeng0824 发表于 2013-4-29 21:10

                               
登录/注册后可看大图

可以把你的ctl内容贴出来看看。

这个是画图的ctl  是29年每年33-40后的平均
DSET I:\CMAP\cmap\3\cmap\3340.grd
UNDEF -9.96921e+36
TITLE cmap patend
xdef  144 linear  1.25  2.5
ydef  72 linear  -88.75  2.5
ZDEF  1  LEVELS 0
TDEF  29 LINEAR  00z10JUN1979  365dy
vars  1
precip   0 99 precippentad                     
ENDVARS

mofangbao

提取之后的缺测为-9.99e+08