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(ffTransformbthig,/INVERSE,/CENTER)

img04 = IMAGE(HffTransform,/CURRENT,LAYOUT=[3,1,3],title='High-pass Butterworth')

END

效果图