高斯格点转换为2*2均匀格点，这样做可以吗？

麦田_smile

下面是原始资料的ctl：
DSET e:\wy\ocean.137.gdat
OPTIONS TEMPLATE
TITLE mom2 output
UNDEF 0.0000000E+00
XDEF 360 LINEAR 0.5 1.
YDEF  132 LEVELS
-73.5000  -70.5092  -67.5459  -64.6277  -61.7709  -58.9906  -56.2998  -53.7093  -51.2275
   -48.8602  -46.6102  -44.4775  -42.4593  -40.5498  -38.7406  -37.0209  -35.3777  -33.7959  -32.2592
   -30.7500  -29.2500  -27.7550  -26.2749  -24.8193  -23.3974  -22.0180  -20.6888  -19.4170  -18.2087
   -17.0688  -16.0012  -15.0084  -14.0917  -13.2512  -12.4855  -11.7920  -11.1670  -10.6055  -10.1013
    -9.6474   -9.2359   -8.8582   -8.5050   -8.1667   -7.8333   -7.5000   -7.1667   -6.8333   -6.5000
    -6.1667   -5.8333   -5.5000   -5.1667   -4.8333   -4.5000   -4.1667   -3.8333   -3.5000   -3.1667
    -2.8333   -2.5000   -2.1667   -1.8333   -1.5000   -1.1667   -0.8333   -0.5000   -0.1667    0.1667
     0.5000    0.8333    1.1667    1.5000    1.8333    2.1667    2.5000    2.8333    3.1667    3.5000
     3.8333    4.1667    4.5000    4.8333    5.1667    5.5000    5.8333    6.1667    6.5000    6.8333
     7.1667    7.5000    7.8333    8.1667    8.5050    8.8582    9.2359    9.6474   10.1013   10.6055
    11.1671   11.7921   12.4855   13.2512   14.0917   15.0084   16.0012   17.0688   18.2087   19.4171
    20.6888   22.0180   23.3974   24.8193   26.2749   27.7550   29.2500   30.7500   32.2644   33.8215
    35.4479   37.1676   39.0009   40.9639   43.0676   45.3176   47.7139   50.2509   52.9176   55.6979
    58.5715   61.5144   64.5000
ZDEF 32 LEVELS
5.00   15.00   25.00   35.00   45.00   55.00   65.00   70.00   85.00   95.00
   105.00  115.4894  127.3993  141.5214  158.4312  178.4312  201.5214  227.3993  255.4895  285.00
   315.00  347.6795  387.6795  437.6795  497.6795  565.00  635.00  789.4041 1174.4041 1874.4041
  2889.4041 4135.00
TDEF 4 LINEAR sep0137 1mo
VARS 7
tt    32 99 temperature
s    32 99 salinity
u    32 99 u current
v    32 99 v current
hflx 1  99 heat flux
tx   1  99 taux
ty   1  99 tauy
ENDVARS

我现在只需要level 5.00的tt资料，我用的方法是grads中的re命令，gs如下：
'reinit'
'open e:\wy\ocean_1.ctl'
'set gxout fwrite'
'set fwrite e:\wy\ocean_cz.grd'
'set z 1'
tt=1
while(tt<=4)
'set t ' tt
'set x 1 360'
'set y 1 132'
'd re(tt,180,linear,0,2,69,linear,-72,2,ba)'
*'d re(v,180,linear,0,2,89,linear,-88,2,ba)'
tt=tt+1
endwhile
'disable fwrite'
'reinit'
;

下面给出插值前后对比，两者相似，就是不懂这种命令能否这么使用，所以上来咨询一下：
下图是插值前：

下图是插值后：

mofangbao

图形相似应该没啥问题

麦田_smile

mofangbao 发表于 2013-9-5 09:23

                               
登录/注册后可看大图

图形相似应该没啥问题

是的呢 就是不懂这样做是否可以 如果是均匀插成均匀肯定没问题

lxlxllx89

高斯插成均匀格点一般也用双线性插值

麦田_smile

lxlxllx89 发表于 2013-9-5 09:52

                               
登录/注册后可看大图

高斯插成均匀格点一般也用双线性插值

那这样做可以吗亲

lxlxllx89

麦田_smile 发表于 2013-9-5 09:54

                               
登录/注册后可看大图

那这样做可以吗亲

我么用过不敢给定论，但是估计是可以的

麦田_smile

lxlxllx89 发表于 2013-9-5 09:55

                               
登录/注册后可看大图

我么用过不敢给定论，但是估计是可以的

谢谢 我也举得可以  因为插值前后还是很相似的

river

基本上是一样的，我觉得可以用

麦田_smile

river 发表于 2013-9-6 20:35

                               
登录/注册后可看大图

基本上是一样的，我觉得可以用

有一个问题 这个放在1.9版本做是没问题的 但是2.0的反而不行 奇怪jng

river

麦田_smile 发表于 2013-9-6 20:36

                               
登录/注册后可看大图

有一个问题 这个放在1.9版本做是没问题的 但是2.0的反而不行 奇怪jng

你可以再看一下官方文档里re()命令在新版中有什么变化