求助：利用matlab，站点数据插值到格点数据遇到的问题

kc121 · 发表于 2012-6-11 21:26:27

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

目前手里有320*160和不等距格点数据，想插值成等间距的144*73格点。使用了interp2()做出来了。可发现图形和插值前差异很大，很明显就是插值函数有问题，应该使用离散点插值的函数来完成。看了
http://bbs.06climate.com/forum.php?mod=viewthread&tid=7144
试着做了一下，可弹出了错误，说是
Error using TriScatteredInterp
Input data must be specified in column-vector format.
我的程序是这样的：
lat和lon是已经的一些经纬度，我先
[qx,qy]=meshgrid(lat,lon); %qx和qy是320*160大小
[new_qx,new_qy]=meshgrid(new_lat,new_lon); %new_qx和new_qy是144*73大小
F=TriScatteredInterp(qx,qy,hgt_500hpa)%hgt_500hpa大小也是320*160的。
不知道哪里出了问题，总是不运行。请高手指点一下呀。比较急！谢谢！

无径之林 · 发表于 2012-6-12 09:31:03

用griddata就可以了，去搜搜interp2和griddata区别就知道了！

wlzhongouc · 发表于 2012-6-11 22:51:57

TriScatteredInterp(qx,qy,hgt_500hpa)中的qx和qy要用hgt_500hpa变量对应的经纬度原值
然后[new_qx,new_qy]=meshgrid(new_lat,new_lon); 是对应你要插值成的网格
最后插出来的是qz=F(new_qx,new_qy);

wlzhongouc · 发表于 2012-6-11 22:54:36

用interp2也要注意网格是否匹配好了应该也是可以的

kc121 · 发表于 2012-6-15 14:43:33

说说我的经验吧。TriScatteredInterp()需要在使用前把数据转成列数据，然后运行就没问题了，不过插值出来的结果和interp2()结果没啥差别。需要是数组形式的去运行就报错。griddata()，我在用320*160格点插值到144*73格点时，总是提示：Maximum variable size allowed by the program is exceeded.所以没办法插值。interp2()可用来插值格点到站点，但是不规则格点插值到规则格点时就会出现差异很大。

tintin · 发表于 2012-8-31 16:58:29

不规则的网格应该用griddata

风儿飘飘 · 发表于 2014-12-1 21:54:14

利用interp2函数时存在缺测值时好像不能计算

		自动登录	找回密码
密码			立即注册

[程序设计] 求助：利用matlab，站点数据插值到格点数据遇到的问题

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