爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8848|回复: 17

[程序设计] 利用matlab做双线性插值遇到的问题

[复制链接]

新浪微博达人勋

发表于 2015-12-7 15:16:02 | 显示全部楼层 |阅读模式

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

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

x
   使用matlab进行双线性插值,原数据分辨率为1.4°*1.4°,为模式输出结果,想插值成0.5*0.5的数据,代码如下

       fil1 = ('G:\IAP_hindcast_output\198003\1980030100\FD14IAP5_lh23_5_Hindcast.cam2.h0.1980-06.nc');
       prep=ncread(fil1,'PRECT');
       a =prep*1000*60*60*24
       lon=ncread(fil1,'lon');
       lat=ncread(fil1,'lat');
        [xx,yy]=meshgrid(lon,lat);
        [xi,yi]=meshgrid(0.25:0.5:359.75,0.25:0.5:179.75);
        c=interp2(xx,yy,a,xi,yi);



却显示如下错误
       Error using griddedInterpolant
         The size of the grid point arrays should match the size of the value arrays.
       Error in interp2/makegriddedinterp (line 220)
         F = griddedInterpolant(varargin{:});
          Error in interp2 (line 133)
       F = makegriddedinterp(X, Y, V, method);
     Error in chazhi1 (line 17)
           c=interp2(xx,yy,a,xi,yi);


请大神帮忙,实在不知道哪里出问题,多谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-7 17:20:25 | 显示全部楼层
显示的错误是你的xx yy与a的大小不符,是不是你a的时次不止一个?
建议检查下xx yy a的大小
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-7 17:22:10 | 显示全部楼层
泰达巴神 发表于 2015-12-7 17:20
显示的错误是你的xx yy与a的大小不符,是不是你a的时次不止一个?
建议检查下xx yy a的大小

看到二楼的头像,好熟悉,话说高中买的利物浦球衣还在
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-7 18:07:57 | 显示全部楼层
泰达巴神 发表于 2015-12-7 17:20
显示的错误是你的xx yy与a的大小不符,是不是你a的时次不止一个?
建议检查下xx yy a的大小

找到问题了,是xx,yy对应的经纬反了,多谢多谢!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-7 18:09:03 | 显示全部楼层
rs06 发表于 2015-12-7 17:22
看到二楼的头像,好熟悉,话说高中买的利物浦球衣还在

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

新浪微博达人勋

发表于 2015-12-7 20:27:40 | 显示全部楼层
橙子鲜 发表于 2015-12-7 18:07
找到问题了,是xx,yy对应的经纬反了,多谢多谢!

你们师门都做模式相关的东西啊?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-12-7 20:28:41 | 显示全部楼层
rs06 发表于 2015-12-7 17:22
看到二楼的头像,好熟悉,话说高中买的利物浦球衣还在

哈哈,我名字加头像是我喜欢的两个球队加一个球员
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-12-8 09:43:21 | 显示全部楼层
泰达巴神 发表于 2015-12-7 20:27
你们师门都做模式相关的东西啊?

有些人是做模式的,不过我是处理他们转出来的结果
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-7 03:43:37 | 显示全部楼层
橙子鲜 发表于 2015-12-7 18:07
找到问题了,是xx,yy对应的经纬反了,多谢多谢!

您好。最近也在处理模式数据,参考了您的方法,得到结果了。但是有个疑问,现在是二维数据,如果是三维数据,那类似的,meshgrid就生成三维网格,用interp3函数,但是原始的nc文件读出的时间是字符型,请问这个您是怎么处理的?希望得到您的帮助。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-5-19 10:04:38 | 显示全部楼层
那个a是什么意思?用到了吗?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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