IDL实现同态滤波

发布时间: 2014-01-18 所属分类: 文章标签:

当前位置：

麻辣GIS » GIS笔记 » 数字图像处理（IDL） » IDL实现同态滤波

什么是同态滤波

同态滤波是把频率过滤和灰度变换结合起来的一种图像处理方法，它依靠图像的照度/ 反射率模型作为频域处理的基础，利用压缩亮度范围和增强对比度来改善图像的质量。

同态滤波原理

将像元灰度值看作是照度和反射率两个组份的产物。由于照度相对变化很小，可以看作是图像的低频成份，而反射率则是高频成份。通过分别处理照度和反射率对像元灰度值的影响，达到揭示阴影区细节特征的目的。

同态滤波流程

全部代码

;Name: 同态滤波
;Time: 2013-12-19
;Author: 尹全超

PRO item_b_2
  
  ;读取图像
  file =  DIALOG_PICKFILE(/READ, FILTER = '*.bmp')
  binary_img = READ_BMP(file)
  img01 = image(binary_img,LAYOUT=[2,1,1],title='Original')
  
  ;过滤无效数据(避免灰度为0)
  ind0=where(binary_img eq 0)
  IF(ind0[0] ne -1) THEN BEGIN
    binary_img[ind0]=1
  ENDIF
  
  ;ln运算
  lndata=alog(binary_img)

  ;快速傅里叶变换
  ffTransform = FFT(lndata, /CENTER)

  ;频域滤波
  ;巴特沃斯 高通滤波
  ;交互输入D0和阶数
  low_row=N_elements(ffTransform[0,*])
  low_col=N_elements(ffTransform[*,0])
  order_low=2
  d0_low=60
  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+(d0_low/float(sqrt(ULONG((i-low_row)*(i-low_row)+(j-low_col/2)*(j-low_col/2)))))^(2*order_low))
    ENDFOR
  ENDFOR
  
  ;反傅里叶变换
  LffTransform=FFT(ffTransform*btlow,/INVERSE,/CENTER)
  
  ;取指数
  reldata=exp(LffTransform)
  img02 = IMAGE(reldata,/CURRENT,LAYOUT=[2,1,2],title='Result')
  
END