爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 17856|回复: 32

[讨论] 二维空间滤波

[复制链接]

新浪微博达人勋

发表于 2017-3-21 16:26:51 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 夏夜 于 2017-3-21 16:36 编辑

想用一下二维空间高通滤波,以去掉大尺度信号。找了下相关资料,没有固定的方法,然后问了一位大牛,回复如下:

We used a 2-dimensional loess smoother to construct spatially low-pass filtered SSH fields that were then subtracted from the raw SSH fields to obtain the spatially high-pass filtered SSH fields. The above numbers (20 deg of longitude by 10 deg of latitude) are the half-amplitude points of the filter transfer function of the loess smoother, i.e., the wavelengths at which the filter transfer function of the loess smoother decays to 0.5.

Note that it is more common for filters to be characterized by the half-power points of the filter transfer function, which is the wavelength at which the square of
the filter transfer function decays to 0.5.



在解释为何选取20*10°的滤波步长时,有些不太明白,请教大家一下,谢谢啦

half-power points(half-amplitude )是什么意思?衰减至0.5该怎么理解?








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

新浪微博达人勋

 楼主| 发表于 2017-3-22 14:36:47 | 显示全部楼层
本帖最后由 夏夜 于 2017-3-22 15:48 编辑

Here is a 2D geospatial data filter.The description in the mfile explain difference between filter and smooth.

Because degrees of longitude are not equally spaced,just from a Matlab point of view the logic is consistent.  Welcome supplement.

The attachment from:

http://cn.mathworks.com/matlabce ... spatial-data-filter    (need account for login)
  

filt2_v1.1.zip

2.85 MB, 下载次数: 101, 下载积分: 金钱 -5

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

新浪微博达人勋

发表于 2017-3-21 20:21:28 | 显示全部楼层

回帖奖励 +5 金钱

2-dimensional loess smoother
这个现阶段还是神器级别,各家会的都还捂着呢
听说德农的张平老师在发nature以后,开始往外给滤波器了,不知真假
楼楼可以试试运气哦

还有,这个smoother实际上不是滤波器,而是名副其实的smoother
但是大家都把他当滤波器用
所以楼楼也就不要纠结截断波长了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-21 20:34:07 | 显示全部楼层
本帖最后由 夏夜 于 2017-3-21 20:35 编辑
ljchen1989 发表于 2017-3-21 20:21
2-dimensional loess smoother
这个现阶段还是神器级别,各家会的都还捂着呢
听说德农的张平老师在发nat ...


感谢回复。
smoother 是指matlab中的smooth?
对一维时间序列,smooth当作低通滤波来看,但是cut_off frequency 指定不明确,与其他滤波器(lanczos)相比,出入较大。
二维空间低通滤波,用smooth这个函数还没尝试如何操作。

张老师也是位大牛,不了解这回事,先去找找他的文章


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

新浪微博达人勋

发表于 2017-3-21 20:56:25 | 显示全部楼层

回帖奖励 +5 金钱

点个赞 一直关注
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-21 22:28:50 | 显示全部楼层
夏夜 发表于 2017-3-21 20:34
感谢回复。
smoother 是指matlab中的smooth?
对一维时间序列,smooth当作低通滤波来看,但是cut_off  ...

并不是matlab的smoother啦
只是可以通过matlab实现

loess的功能是一个smoother的功能,所以阶段波长什么的,意义不大

楼楼拿到张老师的loess smoother,记得分享哦,我等也借鉴一下
现在的matlab程序算得可慢可慢啦~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-21 22:49:04 | 显示全部楼层
本帖最后由 夏夜 于 2017-3-22 09:03 编辑
ljchen1989 发表于 2017-3-21 22:28
并不是matlab的smoother啦
只是可以通过matlab实现


matlab的smooth也有个loess, 其他的不了解。

家园里有相关的问题,貌似都没有效的回复。。。




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

新浪微博达人勋

发表于 2017-3-22 14:09:16 | 显示全部楼层

回帖奖励 +5 金钱

{:5_214:}{:5_214:}{:5_214:}{:5_214:}{:5_214:}{:5_214:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-3-22 14:09:32 | 显示全部楼层
{:5_214:}{:5_214:}{:5_214:}{:5_214:}{:5_214:}{:5_214:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-3-22 14:48:21 | 显示全部楼层

回帖奖励 +5 金钱

点个赞 一直关注
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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