关于CTL中undef设置的问题

凌寒

我的理解是，undef应该设置成一个无效值。
在卫星降水估计中，
1）如果将undef设置为0，则画出来的图为：
根据这个图示是不是可以认为网格点中的数据都为0？
2）可是将undef改为-999后，图为：
说明并不是所有的网格点数据都为0

这是为什么呢？

尽头的尽头

undef是缺测值啊，跟你数据对应起来就好了

凌寒

尽头的尽头 发表于 2013-7-18 10:11

                               
登录/注册后可看大图

undef是缺测值啊，跟你数据对应起来就好了

undef应该是grads不用去理会的值，我想问的是，为什么我选用不同的undef值，会得出矛盾的结论呢

mofangbao

传个测试数据上来

凌寒

mofangbao 发表于 2013-7-18 10:32

                               
登录/注册后可看大图

传个测试数据上来

请见附件。

mofangbao

凌寒 发表于 2013-7-18 10:42

                               
登录/注册后可看大图

请见附件。

ctl？

凌寒

mofangbao 发表于 2013-7-18 10:52

                               
登录/注册后可看大图

ctl？

o 不好意思
CTL：

dset FY2E_PRE_003_OTG_20110808_0800.AWX
title FY2E_PRE
options yrev
undef -999
xdef 1001 linear 55 0.1
ydef 501 linear 0 0.1
zdef 1 linear 300 1
tdef 1 linear 03z30jun2013 1hr
vars 1
pre 0 -1,40,1 PRE(MM/HOUR)
endvars


gs：
'reinit'
'open /home/FY2E_PRE/PRE.ctl'

'set gxout shaded'
'set grads off'

'set mpdset cnworld'
'set map 1'
'set clevs 0 10 20 30'
'set ccols 0 2 3 4 9'

'd pre'
'cbarn 1.0 0 5.5 0.2'

'q file'
t=sublin(result,3)
FulFile=subwrd(t,2)
length=strlen(FulFile)
PostFileName=substr(FulFile,35,length)
length2=strlen(PostFileName)
PreFileName=substr(PostFileName,1,length2-5)

'q ctlinfo'
TitleString=sublin(result,2)
Title=subwrd(TitleString,2)

'draw title 'Title''

尽头的尽头

凌寒 发表于 2013-7-18 10:19

                               
登录/注册后可看大图

undef应该是grads不用去理会的值，我想问的是，为什么我选用不同的undef值，会得出矛盾的结论呢

没有矛盾啊，是因为你缺测值设定的和数据不匹配，画的图就是错的，结果是错的，不是矛盾的

river

设置了正确的缺测值，就能画出正确的图了

凌寒

mofangbao 发表于 2013-7-18 10:52

                               
登录/注册后可看大图

ctl？

我发现这类数据（降水估计）是用两个字节来表示一个数据，GrADS如何区分用多个字节表示一个数据和用1个字节来表示一个数据的情况呢？