爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4507|回复: 2

[程序设计] Matlab做线性插值时出错

[复制链接]

新浪微博达人勋

发表于 2019-8-23 19:44:56 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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);

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-8-24 14:28:50 | 显示全部楼层
检查一下维度是不是搞错了?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-9-26 14:55:28 | 显示全部楼层
楼主  我也是准备对day数据进行插值  qq:1355119420  交流一下啊
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表