麻辣GIS微信平台

更多 GIS 干货

微信关注不错过

IDL实现巴特沃斯(butterworth)高通滤波

IDL其实自带了巴特沃斯滤波函数,但是按照IDL官方文档来看,其所带的butterworth函数,与我要求的并不一致,所以这里按照自己的方法实现了一个巴特沃斯高通滤波。

频域滤波原理

巴特沃斯滤波属于频域滤波的一种,频域滤波的原理是卷积定理,即:

新的图像=原来的图像*频域滤波函数。而本文所介绍的巴特沃斯高通函数是频域滤波的一种。

巴特沃斯高通滤波函数

全部代码

;Name: 巴特沃斯高通滤波
;Time: 2013-12-18
;Author: 尹全超

PRO item_b_1

;读取图像
file = DIALOG_PICKFILE(/READ, FILTER = '*.bmp')
binary_img = READ_BMP(file)
img01 = image(binary_img,LAYOUT=[3,1,1],title='Original')

;快速傅里叶变换
ffTransform = FFT(binary_img, /CENTER)
powerSpectrum = ABS(ffTransform)^2
scaledPowerSpect = ALOG10(powerSpectrum)
img02 = IMAGE(scaledPowerSpect,/CURRENT,LAYOUT=[3,1,2],title='FFT')

;巴特沃斯高通滤波
;默认2阶,D0=60
hig_row=N_elements(ffTransform[0,])
hig_col=N_elements(ffTransform[
,0])
order_hig=2
d0_hig=60
bthig=dblarr(hig_col,hig_row)
FOR i=1,hig_row DO BEGIN
FOR j=1,hig_col DO BEGIN
bthig[j-1,i-1]=1.0/(1.0+(d0_hig/float(sqrt(ULONG((i-hig_row)(i-hig_row)+(j-hig_col/2)(j-hig_col/2)))))^(2order_hig))
ENDFOR
ENDFOR
HffTransform=FFT(ffTransform
bthig,/INVERSE,/CENTER)
img04 = IMAGE(HffTransform,/CURRENT,LAYOUT=[3,1,3],title='High-pass Butterworth')
END

效果图

麻辣GIS-Sailor
作者:
GIS爱好者,学GIS,更爱玩GIS。
微博关注
手机阅读
赞赏支持
手机阅读
微信捐助麻辣GIS
微信打赏
支付宝捐助麻辣GIS
支付宝打赏

声明

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

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

3.如果本文对您有所帮助,请您对我们团队进行 打赏捐助让我们在传播3S的路上可以走得更远,不胜感激。PS:如果你是学生党,请优先把经费用于购买学习资料 以及 与小哥哥/小姐姐约会上:-)

如您有疑问,可在文末留言,或到麻辣GIS QQ群中提问。

QQ群(一):291616564 加入QQ群

QQ群(二):166408035 加入QQ群

QQ群(三):627853279 加入QQ群

QQ群(四):436386604 加入QQ群

QQ群(五):606176554 加入QQ群

QQ群(六):946178380 加入QQ群

QQ群(七):861271808 加入QQ群

新疆互助群:910717627 加入QQ群

微信公众号:malagis,扫描右边二维码直接关注。

发表评论