nc文件转为2进制文件出错

yx88229 · 发表于 2012-12-20 11:26:55

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

nc文件的ctl为

登录/注册后可看大图

nc文件画出的psl全球图为
QQ截图20121219230647.png

nc文件转为2进制的gs为：
'reinit'
'sdfopen f:/psl_mh.nc'
'set gxout fwrite'
'set fwrite f:/psl_mh.dat'
'set lon 0 358.125'
'set lat -88.5722 88.5722'
t0=1
while(t0<=1200)
'set t 't0
'd psl'
t0=t0+1
endwhile
'disable fwrite'

而用2进制文件画出来psl的图却是

等值线交叉了，请问是怎么回事？
编写的ctl为：
dset f:/psl_mh.bin
title MPI-ESM-P model output prepared for CMIP5 mid-Holocene
undef 9.96921e+36
xdef 192 linear 0 1.875
ydef 96 levels -88.5722 -86.7225 -84.862 -82.9989 -81.135 -79.2706 -77.4059 -75.5411 -73.6761 -71.8111 -69.9461 -68.081 -66.2159 -64.3507 -62.4856 -60.6204 -58.7552 -56.89 -55.0248 -53.1596 -51.2944 -49.4292 -47.5639 -45.6987 -43.8335 -41.9682 -40.103 -38.2377 -36.3725 -34.5072 -32.642 -30.7767 -28.9115 -27.0462 -25.181 -23.3157 -21.4505 -19.5852 -17.72 -15.8547 -13.9894 -12.1242 -10.2589 -8.39367 -6.52841 -4.66315 -2.79789 -0.93263 0.93263 2.79789 4.66315 6.52841 8.39367 10.2589 12.1242 13.9894 15.8547 17.72 19.5852 21.4505 23.3157 25.181 27.0462 28.9115 30.7767 32.642 34.5072 36.3725 38.2377 40.103 41.9682 43.8335 45.6987 47.5639 49.4292 51.2944 53.1596 55.0248 56.89 58.7552 60.6204 62.4856 64.3507 66.2159 68.081 69.9461 71.8111 73.6761 75.5411 77.4059 79.2706 81.135 82.9989 84.862 86.7225 88.5722
zdef 1 linear 0 1
tdef 1200 linear 12Z16JAN1850 1yr
vars 1
psl 0 -999
endvars

急求啊！！！

yx88229 · 发表于 2012-12-20 11:34:34

这是转为2进制文件画出来psl的图

mofangbao · 发表于 2012-12-20 14:02:46

'set lon 0 358.125'
'set lat -88.5722 88.5722'
你的y方向好像是列出来的，并不是线性的，你把这两句话去掉看看

yx88229 · 发表于 2012-12-20 14:41:32

mofangbao 发表于 2012-12-20 14:02

登录/注册后可看大图

'set lon 0 358.125'
'set lat -88.5722 88.5722'
你的y方向好像是列出来的，并不是线性的，你把这两句话 ...

把那两句去掉后也是错的
图出来是 QQ截图20121220143943.png

river · 发表于 2012-12-20 18:54:33

'set lon 0 358.125'
'set lat -88.5722 88.5722'
还是改成'set x 1 192'
'set y 1 96' 吧，这样比较保险

yx88229 · 发表于 2012-12-20 21:37:22

river 发表于 2012-12-20 18:54

登录/注册后可看大图

'set lon 0 358.125'
'set lat -88.5722 88.5722'
还是改成'set x 1 192'

改了之后还是不行，画出来跟上面一样

river · 发表于 2012-12-20 22:13:21

yx88229 发表于 2012-12-20 21:37

登录/注册后可看大图

改了之后还是不行，画出来跟上面一样

那就是ctl的问题了吧，你看看是不是ctl里的时间那一项有问题啊，你设置的是1yr 原始文件应该是1mo啊

yx88229 · 发表于 2012-12-21 15:04:28

river 发表于 2012-12-20 22:13

登录/注册后可看大图

那就是ctl的问题了吧，你看看是不是ctl里的时间那一项有问题啊，你设置的是1yr 原始文件应该是1mo啊

哦，这个是笔误，后来改了之后还是一样的不行。。。

tulipwaugh · 发表于 2012-12-21 16:07:37

试试用！ncdump -c 文件路径查看下原nc的缺省值，有可能是缺省值有问题...

小草 · 发表于 2012-12-21 16:10:10

有可能是x方向上抱点的问题，即0和360为同一个点，fwrite时应该写'set lon 0 360'，你可以试一下。

		自动登录	找回密码
密码			立即注册

[分享资料] nc文件转为2进制文件出错

登录后查看更多精彩内容~