求助！关于re函数的使用

寄居蟹

使用re函数出现问题。我想把1°*1°海温资料使用re函数转化成2°*2°的资料海温原始数据资料的ctl描述文件
dset g:\sst_1.dat
undef -999.0
title hsst
xdef 360 linear -180.0  1.0
ydef 180 linear  -90.  1.0
zdef 1 levels 1000
tdef 612 linear 00z01jan1950  1mon
vars 1
sst 0 99  globle ocean sst
endvars
我读取时写的re的gs文件
'reinit'
'open G:\sst1.ctl'
'set gxout fwrite'
'set fwrite G:\sst2.grd'
'set t 1'
'set lev 1000'
'set x 1 360'
'set y 1 180'
*'set lat -90 89'
*'set lon 0 360'
*'define s1=re(sst,2,2,bl)'
'define s1=re(sst,180,linear,0,2,91,linear,-90,2,bl)'
'd s1'
'disable fwrite'
;
数据读出后搭配的ctl描述文件

dset g:\sst2.grd
undef -999.0
title  sst
xdef 180 linear 0  2.0
ydef 91 linear  -90.  2.0
zdef 1 levels 1000
tdef 1 linear 00z01jan1950  1mon
vars 1
sst 0 99  globle ocean sst
endvars

同一时次画出的图

原数据图

读出图

请教，从图中看感觉数据应该是读入进去了，就是经纬度偏移太厉害，海洋上的数据都到陆地上了。不管是re函数的简写还是所有参数都设置的写法都用了，不管是设置经纬度set lat,lon还是设置格点set x,y都试了，结果都是这种。
请问这种re函数应该怎么使用？我的gs和ctl有什么问题？谢谢！！

kongfeng0824

xdef和ydef前后的ctl不一致呢？

传说中的谁

原来的x是-180到180，你的ctl里面是0 到360，是不是这个不一样造成的？

寄居蟹

kongfeng0824 发表于 2013-7-24 10:05

                               
登录/注册后可看大图

xdef和ydef前后的ctl不一致呢？

re函数是改变数据分辨率的，原来的1*1变成2*2，所以Ctl不同

寄居蟹

传说中的谁 发表于 2013-7-24 10:20

                               
登录/注册后可看大图

原来的x是-180到180，你的ctl里面是0 到360，是不是这个不一样造成的？

不是这个原因哪，改成-180起始后图形又变成这个样子

whmooncool

顶起！也有相关问题呐~~

寄居蟹

whmooncool 发表于 2013-7-25 08:45

                               
登录/注册后可看大图

顶起！也有相关问题呐~~

有木有人知道啊~~

风儿飘飘

刚搜了下相关资料，发现定义'define s1=re(sst,180,linear,0,2,91,linear,-90,2,bl)'，最后一个参数应该为''ba'',而非'bl',试一下，不知解决问题了没

风儿飘飘

回信还减少金钱呢

风儿飘飘

我也会遇到相同的问题