求助：用fwrite读取出的数据所绘图与用原始数据出图不一致

leeol

原始数据：NCEP再分析资料：uwnd.day.1968-1996.ltm.nc

使用fwrite读取数据（GS如下）：
'set fwrite d:\data\850hPa.grd'
'set gxout fwrite'
'sdfopen d:\data\uwnd.day.1968-1996.ltm.nc'
'set lev 850'
'set t 1'
'd uwnd'
'disable fwrite'
;


自己编写的CTL：
dset  d:\data\850hPa.grd                 
TITLE 850hPa uwnd
UNDEF  -9.9E33                                 
XDEF 144 LINEAR  0 2.5           
YDEF 73 LINEAR  -90 2.5         
ZDEF 1 LEVELS 850
TDEF 1 LINEAR 00Z01JAN0001 1440mn      
vars 1                                                     
uwnd     1   t,z,y,x   U Winds
ENDVARS  


用850hPa.grd文件画图：
'open d:\data\850hPa.ctl'
'set lat 0 90'
'set lon 80 160'
'd uwnd'
'printim d:\nc\1.png white'

；



用原始数据uwnd.day.1968-1996.ltm.nc画图：
'sdfopen d:\data\uwnd.day.1968-1996.ltm.nc'
'set lat 0 90'
'set lon 80 160'
'set lev 850'
'set t 1'
'd uwnd'
'printim d:\nc\2.png white'
;



现在存在疑惑，为什么这两个图不一样，望大家帮忙解惑，。

river

提取资料的gs里面把具体的经纬度什么的都设置好看看，即使全球的提取也最好加上set x 1 144和set y 1 73的设置

风儿飘飘

用fwrite提取的二进制数据，缺测为-9.99E+08，你将缺测设置更改一下试试

传说中的谁

原来的X是从0开始还是-180开始？最好一致。这段时间很多这个问题的帖子，楼主去翻翻看别人怎么解决的吧

leeol

river 发表于 2013-8-8 18:11

                               
登录/注册后可看大图

提取资料的gs里面把具体的经纬度什么的都设置好看看，即使全球的提取也最好加上set x 1 144和set y 1 73的设 ...

谢谢，果然是要在fwrite的时候设置好维度，加上set x 1 144和set y 1 73后一致了

leeol

传说中的谁 发表于 2013-8-8 22:23

                               
登录/注册后可看大图

原来的X是从0开始还是-180开始？最好一致。这段时间很多这个问题的帖子，楼主去翻翻看别人怎么解决的吧

管理员你提到的是这个帖子：http://bbs.06climate.com/forum.p ... &extra=page%3D5

我用前后数据绘制全球图确实发现了如上贴中的数据偏移现象，但是个别经纬度的数值有差异，故不能以修改ctl解决。
再用了1L的方法后解决了问题后，应该是之前使用fwrite的维度设定的不够具体。

leeol

http://bbs.06climate.com/forum.p ... page%3D9&page=1

http://bbs.06climate.com/forum.p ... page%3D9&page=1

这2个帖子也是用grads将nc转成grd出现问题，但与我的问题不一样，标记下以供参考

leeol

风儿飘飘 发表于 2013-8-8 19:32

                               
登录/注册后可看大图

用fwrite提取的二进制数据，缺测为-9.99E+08，你将缺测设置更改一下试试

这个不一定，要看原文件的缺测值

奋豆儿

过来学习，fwrite很有用！