爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 32573|回复: 35

[源程序] 逐步订正的高斯插值

[复制链接]

新浪微博达人勋

发表于 2013-11-15 23:28:07 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 wlzhongouc 于 2014-4-27 11:00 编辑

终于下定决心写一个Matlab的高斯插值程序,本插值程序的优点是比其他普通插值程序多了一步误差订正,缺点是在Matlab中运行还是比较慢的,一个问题是Matlab中矩阵的运算相比Fortran要慢得多,程序中的多次循环也使得运算速度大为降低,没想到更好的解决办法

插值结果对比如图(下图为用NCEP海表面气压场的数据插值验证):
Gaussian interpolation.png

源程序: gaussian_interpolation.m (5.69 KB, 下载次数: 198)

评分

参与人数 6金钱 +101 贡献 +28 体力 +200 收起 理由
言深深 + 15 + 2
mofangbao + 20 + 8 + 200
Aires + 18 + 6 很给力!
topmad + 20 + 2
lqouc + 18 + 4 很给力!
kongfeng0824 + 10 + 6

查看全部评分

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

新浪微博达人勋

发表于 2013-11-15 23:29:35 | 显示全部楼层
感谢分享,很不错
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-15 23:30:27 | 显示全部楼层
这个插值刚好插值到陆地之内,不知道是如何实现的?我有时候差值会插值到海洋上去
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-16 00:08:45 | 显示全部楼层
大赞
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-16 08:41:32 | 显示全部楼层
非常强大,赞一个
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-16 09:12:34 | 显示全部楼层
很赞的说
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-16 10:07:48 | 显示全部楼层
学长又发威了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-16 11:23:49 | 显示全部楼层
真心不错哦
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2013-11-16 11:30:09 | 显示全部楼层
哈哈,绘图技巧很棒!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-11-16 21:41:51 | 显示全部楼层

用的数据是海表面气压所以是一个全场都有的数值,看到陆地上没有数值是因为先画的数值场然后叠加陆地。
你提到的,如果只想插值显示陆地上的数值,可以做一个mask海洋中都设为NaN
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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