- 积分
- 251
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-6-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
用Matlab2014a对下载的全球模式数据进行双线性插值,遇到了一个问题,Matlab小白,不会解决,求大神帮忙。
代码如下:
clear all;close all;clc
pr=ncread('D:\CMIP5\MIROC5\pr_day_MIROC5_rcp45_r1i1p1_20200101-20291231.nc','pr')-273.15;%读取。nc数据中的pr
lat=ncread('D:\CMIP5\MIROC5\pr_day_MIROC5_rcp45_r1i1p1_20200101-20291231.nc','lat');%读取读取NC数据的纬度数据
lon=ncread('D:\CMIP5\MIROC5\pr_day_MIROC5_rcp45_r1i1p1_20200101-20291231.nc','lon');%读取读取NC数据的经度数据
pr=pr(:,:,613:1872);
lon=[lon(129:256,:)-128;lon(1:128,:)];
for n=1:3650
pr(:,:,n)=[pr(129:256,:,n);pr(1:128,:,n)];
end;
[xx,yy]=meshgrid(lon,lat);%生成格点
[xi,yi]=meshgrid(-179.75:0.5:179.75,-89.75:0.5:89.75);
pr=permute(pr,[2,1,3]);%交换第一维和第二维的位置
for i=1:3650
pr_MIROC5(:,:,i)=interp2(xx,yy,pr(:,:,i),xi,yi,'linear');
end;
错误如下:
出错 interp2>makegriddedinterp (line 226)
F = griddedInterpolant(varargin{:});
出错 interp2 (line 136)
F = makegriddedinterp(X, Y, V, method,extrap);
|
|