Matlab中m_map使用插值作图探讨

kongfeng0824

参考了家园的帖子，http://bbs.06climate.com/forum.php?mod=viewthread&tid=10120和http://bbs.06climate.com/forum.php?mod=viewthread&tid=7144。自己也尝试使用站点数据画图。
第一种差值方法：
vq = griddata(lat54406,lon54406,degree54406,xq,yq,'natural');%数据插值
第二种插值方法：
[xq,yq] = meshgrid(20:1:200,10:1:90);
F=TriScatteredInterp(lat54406,lon54406,degree54406,'natural');
vq=F(xq,yq);
结果两种方法得到的结果一致。

使用v4差值方法：


不过我还是有亮点不太明白：
1.如何使得插值区域覆盖掉整个中国区域。（PS：我尝试了改变[xq,yq] = meshgrid(20:1:200,10:1:90);但是好像没有任何改变）
2.如何裁剪掉中国国界之外的部分插值区。

我看到了家园的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=12340  我是用的是m_map工具箱http://bbs.06climate.com/forum.php?mod=viewthread&tid=14131。不是matlab自带的。

不知道大家是否有相关的经验分享一下！

斥鷃

griddata插值方式改成'v4'试试？试试其他插值选项~

kongfeng0824

斥鷃 发表于 2013-5-24 00:04

                               
登录/注册后可看大图

griddata插值方式改成'v4'试试？试试其他插值选项~

v4也试过了好像也是范围没什么变化。而且mask有点问题

kongfeng0824

斥鷃 发表于 2013-5-24 00:04

                               
登录/注册后可看大图

griddata插值方式改成'v4'试试？试试其他插值选项~

V4插值效果不如natural的。但是natural的又不能覆盖全中国。而且mask确实不太好弄

斥鷃

kongfeng0824 发表于 2013-5-27 22:24

                               
登录/注册后可看大图

V4插值效果不如natural的。但是natural的又不能覆盖全中国。而且mask确实不太好弄

V4要高点分辨率大概就只能加密网格了，要用那几个插满怕是要自己手工做插值函数代码了吧~

kongfeng0824

斥鷃 发表于 2013-5-27 22:39

                               
登录/注册后可看大图

V4要高点分辨率大概就只能加密网格了，要用那几个插满怕是要自己手工做插值函数代码了吧~

我直接调用IDW护着克里金函数插值。剩下就是mask的问题了

F117_ren_0

楼主用了多少个站的资料画的这个图？

Ryan20100520

这是如何做出来的啊

weizhenfeng

这个插值能否更具体点了，有数据的。如读数据，到插值。加精，谢谢楼主啊

做梦的蚊子

用matlab插值雷达数据，表示数组太大，无法插值的飘过。。。