MATLAB关于四维EC资料插值问题求教！初学者

quentinwan

我下载的EC的资料'Temperature是四维资料（lon,lat,lev,time），经纬度是从33N-39N,115E-125E，0.25*0.25，0.0KM--10KM30层，时间维数为180天，请问我要是想插值到36.7N,117.3E,垂直层为60层，我该怎么写程序，我是matlab的初学者，请大家帮帮我！万分感谢！

quentinwan

盼望着，盼望着，盼望着您能指点我这段代码

我是柿子

为楼主工作量堪忧，目测俺木有简单快捷的方法，坐等大神的出现

二爷名声在外

差值的算法就有很多种，主要是找算法，剩下的编程其实本身不难……

quentinwan

我今早刚睡醒，想了想是不是这样处理，请各位高手给看看：
%i=经度维数计数器；j=维度计数器；k=垂直层数计数器；t=时间维数
station=double(ones(length(lev),length(time));
for t=1:length(time)
   for k=1:length(lev)
       station(k,t)=interp2(lon,lat,temperature(:,:,k,t),117.3,36.7,'spline')；
   end
end


不知道这样写对不

quentinwan

再附上我的变量属性，大家帮我解答下啊


Variables:
    longitude
           Size:       41x1
           Dimensions: longitude
           Datatype:   single
           Attributes:
                       units     = 'degrees_east'
                       long_name = 'longitude'
    latitude
           Size:       25x1
           Dimensions: latitude
           Datatype:   single
           Attributes:
                       units     = 'degrees_north'
                       long_name = 'latitude'
    level   
           Size:       31x1
           Dimensions: level
           Datatype:   int32
           Attributes:
                       long_name = 'model_level_number'
    time     
           Size:       184x1
           Dimensions: time
           Datatype:   int32
           Attributes:
                       units     = 'hours since 1900-01-01 00:00:0.0'
                       long_name = 'time'
    t        
           Size:       41x25x31x184
           Dimensions: longitude,latitude,level,time
           Datatype:   int16
           Attributes:
                       scale_factor  = 0.00153
                       add_offset    = 257
                       _FillValue    = -3.28e+04
                       missing_value = -3.28e+04
                       units         = 'K'
                       long_name     = 'Temperature'
                       standard_name = 'air_temperature'

我是柿子

你不是还要把垂直方向差值成60层吗

quentinwan

我想先水平插值，再垂直插值！有好心人请指导下我，我QQ：123744035

quentinwan

是不是得用griddata函数啊

quentinwan

for i=1:t
dyt=griddata(double(lon'),double(lat'),h_ec,double(temp(:,:,:,i)),118.67,37.43,h_rm,'linear');
end
h_ec=[11.29 10.42 9.59 8.81 8.08 7.38 6.71 6.09 5.50 4.94 4.42 3.94 3.48 3.06 2.67 2.31 1.98 1.68 1.41 1.17 0.95 0.76 0.60 ...
     0.46 0.34 0.24 0.16 0.10 0.06 0.03 0.01].*1000;

h_rm=[0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.60 0.70 0.80 0.90 1.00 1.10 1.20 1.30 1.40 1.50 1.60 1.70 1.80 1.90 2.00 ...
    2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 4.50 4.75 5.00 5.25 5.50 5.75 6.00 6.25 6.50 6.75 7.00 7.25 7.50 7.75 8.00 8.25 8.50 ...
    8.75 9.00 9.25 9.50 9.75 10.00].*1000;