那些个多变的缺测值

nnlovess

对于我这个菜鸟新手来说，真是每天都有新问题，还好咱们的论坛很强大，总能找到解决办法或者解决问题的思路！下面是一个找到缺测值的方法，发来希望能给同样菜鸟的新手一点帮助吧！今天在用grads提取nc文件的时候，出图变成了以下的状态：

  

查了一打的帖子发现应该是missing data的问题，也有很多的方法去判断缺测值应该是多少，可是我查了了以后依然不对（有-9.99e33，-32767等），最后是通过MeteoInfo解决了问题。
步骤也很简单，首先将转换好的grd文件的ctl写好，主要要保证别的地方没有问题，例如：
dset d:\20002009sst\mn20001.grd
undef -9.99e33
xdef 241 linear  -59.875  0.25
ydef 121 linear  10.125  0.25
zdef  1 linear 0 1
tdef  1 linear 00z01mar2000 1mo
vars 1
sst=>sst  0 t,y,x Daily Sea Surface Temperature
endvars

这个的ctl文件缺测值就是错误的，我们不去管它，直接用MeteoInfo打开，用view data可以看到其实是有数据的，而这里面就可以看到我们需要找的missing data，原来是-9.99e8，而不是之前查到的那些值。

将这个值改正过后我们就能正确的显示结果了。

river

这个确实比较麻烦的，我原来是用forrtran读出来看看缺测值

mofangbao

如果你使用GrADS的fwrite输出的话，十有八九会变成-9.99e8的，下次直接用这个值试就行了。
不过楼主的方法在不知道缺测值的时候测试一下也是挺有用的。
类似的方法除了楼上说的fortran输出，还可以直接观察你d一个变量的时候GrADS给出的提示信息来猜一下，或者用兰溪的grads2ascii脚本输出为txt看看。
总之，谢谢楼主分享的经验

凌寒

如果没有缺测值那undef应该怎么设置呢？

mofangbao

凌寒 发表于 2013-7-18 14:20

                               
登录/注册后可看大图

如果没有缺测值那undef应该怎么设置呢？

那就随便写个你数据中不可能出现的值 比如很大或者很小的数

nnlovess

mofangbao 发表于 2013-7-18 08:20

                               
登录/注册后可看大图

如果你使用GrADS的fwrite输出的话，十有八九会变成-9.99e8的，下次直接用这个值试就行了。
不过楼主的方法 ...

嘿嘿，第一次发，谢谢肯定啊！以后会继续努力的！

麦田_smile

没看到哪里能view data

奋豆儿

mofangbao 发表于 2013-7-18 08:20

                               
登录/注册后可看大图

如果你使用GrADS的fwrite输出的话，十有八九会变成-9.99e8的，下次直接用这个值试就行了。
不过楼主的方法 ...

谢谢清风，今天用fwrite提取变量生成.grd文件后，怎么画图显示的值都不对，不知道该怎么解决，在论坛浏览相关帖子时，看到你的回复，一言惊醒梦中人啊。按照你的提示修改缺省值为-9.99e8马上解决问题。但是不知道清风能否解释一下，使用GrADS的fwrite输出的话，缺省值会变成-9.99e8的原因是什么呢？

奋豆儿

论坛确实很强大，高了半天弄不出来的问题，在这里找到解决方法！

mofangbao

奋豆儿 发表于 2013-12-19 15:29

                               
登录/注册后可看大图

谢谢清风，今天用fwrite提取变量生成.grd文件后，怎么画图显示的值都不对，不知道该怎么解决，在论坛浏览 ...

没有为什么。。。GrADS内部缺测值就是这个