如何用matlab 插值NC数据调整分辨率

veronicasyd

前两天刚发过一个帖子询问关于将CMIP5中的模式，插值成统一分辨率的方法。经过Aires的指导，初步确定用matlab先读取NC数据后，在进行插值。本来打算进行双线性插值，但是读取数据后，发现那个降水率的数据PrData是三维的[lat,lon,time]，数据如下：

所以后来打算用三维差值下看看，程序如下：
PrData  = ncread(filename,'pr');
TimeData  = ncread(filename,'time');
LonData  = ncread(filename,'lon');
LatData  = ncread(filename,'lat');
[xx,yy,zz]=meshgrid(LatData,LonData,TimeData);
[xi,yi,zi]=meshgrid(128:1:135,43:1:48,TimeData);%我只想要得到128~135E，43~48N这个范围的数据~
vi=interp3(xx,yy,zz,PrData,xi,yi,zi);

但是出来的结果是vi都是NaN, 不知道到底是哪里出了问题。

PS：如果想用双线性插值的话，这样的数据要怎么处理。
（因为初步接触这些程序不久，所以好多都不懂，百度和论坛上的资料看了好多，还是很迷茫，希望能有大神解惑~）

沙颖凯

veronicasyd 发表于 2014-5-4 18:54
好像不是的感觉，改了之后提示：
??? Error using ==> interp3 at 128
Matrices X,Y and Z must be the ...

PrData=permute(PrData, [2 1 3]);
[xx,yy,zz]=meshgrid(LonData,LatData,TimeData);

沙颖凯

改成
[xx,yy,zz]=meshgrid(LonData,LatData,TimeData);

veronicasyd

沙颖凯 发表于 2014-5-4 18:47
改成
[xx,yy,zz]=meshgrid(LonData,LatData,TimeData);

好像不是的感觉，改了之后提示：
??? Error using ==> interp3 at 128
Matrices X,Y and Z must be the same size as V.

veronicasyd

沙颖凯 发表于 2014-5-4 18:59
PrData=permute(PrData, [2 1 3]);
[xx,yy,zz]=meshgrid(LonData,LatData,TimeData);

哇！！成功了。点无数个赞~不过再问个小问题~~如何导出三维数据呢~~导出成Excel是不是有利于之后的分析。。

沙颖凯

veronicasyd 发表于 2014-5-4 19:37
哇！！成功了。点无数个赞~不过再问个小问题~~如何导出三维数据呢~~导出成Excel是不是有利于之后的分析。 ...

xlswrite(FILE, ARRAY, SHEET, RANGE)

veronicasyd

沙颖凯 发表于 2014-5-4 20:04
xlswrite(FILE, ARRAY, SHEET, RANGE)

没办法用xls导出，好像是因为是三维的，错误提示：
??? Error using ==> xlswrite at 154
Dimension of input array cannot be higher than two.

沙颖凯

veronicasyd 发表于 2014-5-4 20:10
没办法用xls导出，好像是因为是三维的，错误提示：
??? Error using ==> xlswrite at 154
Dimension of ...

想想也是，不过128*64*1872对EXCEL来讲太大了，感觉要用EXCEL就没有什么好办法了

veronicasyd

沙颖凯 发表于 2014-5-4 20:14
想想也是，不过128*64*1872对EXCEL来讲太大了，感觉要用EXCEL就没有什么好办法了

数据在插值的时候只算了，128~135E，43~48N这个范围的数据~所以数据量也不是很大，是6*8*1872~但是就是不知道怎么样导出这样的三维数据~