波变换

开发平台：

Matlab

TongTaiLvBo.m：源码内容

							clear
clc
img=double(imread('E:exapleexp2.bmp'));
%i=imread('0,jpg');
%i=rgb2gray(img);
figure,subplot(2,2,1),imshow(img./255),%%%%%图像灰度值归一化
title('the initial picture');
[M,N]=size(img);
img=img+1;
lnimg=log2(img);%取对数
Fimg=fft2(lnimg);%傅里叶变换
P=fftshift(Fimg);%将频域原点移到图像中心；
subplot(2,2,2),imshow(uint8(abs(P)),[]),title('滤波前的频谱图像')
%Do=1.4;  %一般是方差 (滤波器的高通截止频率)
Do=2;
for u=1:M 
    for v=1:N  
       D(u,v)=sqrt((u-floor(M/2))^2+(v-floor(N/2))^2);; %点（i,j）到频率平面原点的距离 
        H(u,v)=(1.0)*(1.0-exp(-D(u,v)^2/(2*Do^2)))-0.6;%同态滤波器函数
    end                
end
hImg=Fimg.*H; Q=fftshift(hImg);
subplot(2,2,3)  ,imshow(uint8(abs(Q)),[])  ,title('滤波后的频谱图像')
gImg=ifft2(hImg);    %反傅立叶变换 
Y=exp(gImg);      %取指数 
subplot(2,2,4),imshow((abs(Y)),[]),title(' the hormomorphic image')

						