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');巴特沃斯低通滤波
;默认1阶,D0=15
low_row=N_elements(ffTransform[0,])
low_col=N_elements(ffTransform[,0])
order_low=1
d0_low=15
btlow=dblarr(low_col,low_row)
FOR i=0,low_row-1 DO BEGIN
FOR j=0,low_col-1 DO BEGIN
btlow[j,i]=1.0/(1.0+(float(sqrt(ULONG((i-low_row/2)(i-low_row/2)+(j-low_col/2)(j-low_col/2)))/d0_low))^(2order_low))
ENDFOR
ENDFOR
LffTransform=FFT(ffTransformbtlow,/INVERSE,/CENTER)
img03 = IMAGE(LffTransform,/CURRENT,LAYOUT=[3,1,3],title='Low-pass Butterworth')
END
效果图
相关阅读
声明
1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。
2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。
请问一下,我刚刚学IDL low_row=N_elements(ffTransform[0,])
low_col=N_elements(ffTransform[,0])
这两句代码显示有错误是怎么回事呀
你提供这点信息,我也无法判断是为什么啊。