麻辣GIS微信平台

更多 GIS 干货

微信关注不错过

GIS数据处理-更新IDL实现巴特沃斯(butterworth)高通滤波

在之前的文章:IDL自带函数的巴特沃斯(butterworth)高通/低通滤波中,我介绍了IDL自带巴特沃斯(butterworth)函数实现了巴特沃斯(butterworth)高通/低通滤波。在文章中,我说没有关于高通巴特沃斯滤波的实现,今天发现是错误的。

IDL自带巴特沃斯(butterworth)滤波器

\[\frac{1}{\sqrt{1+(\frac{\omega}{\omega_{c}})^{2N}}}\]

那么根据推导,高通的时候,N>0,低通的时候,N<0,N为阶娄。

所以低通滤波器应为

\[\frac{1}{\sqrt{1+(\frac{\omega}{\omega_{c}})^{2N}}}\] N>0

高通滤波器应为

\[\frac{1}{\sqrt{1+(\frac{\omega}{\omega_{c}})^{2N}}}\] N<0

结论

通过butterworth函数的ORDER参数控制高通和低通

示例

;AUTHOR:尹全超
PRO butterworth_h_l
filter = BUTTERWORTH(500,500,CUTOFF=50,/ORIGIN,ORDER=-1)
img = image(filter,title="高通",font_name = 'Microsoft Yahei')
filter = BUTTERWORTH(500,500,CUTOFF=5,/ORIGIN,ORDER=1)
img1 = image(filter,title="低通",font_name = 'Microsoft Yahei')
END

效果如图

麻辣GIS-Sailor

作者:

GIS爱好者,学GIS,更爱玩GIS。

声明

1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。

2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。

手机阅读
微信打赏
支付宝打赏
手机阅读
麻辣GIS微信打赏
请小编吃包辣条
麻辣GIS支付宝打赏
请小编吃支雪糕

留言板(小编看到第一时间回复)