%同态滤波实现
%%
clc;clear all;close all;
img=imread('C:\Users\meng\Desktop\124.bmp');
subplot(2,2,1),imshow(img),title('the initial picture');
img=im2double(img);
lnimg=log(img+0.000001);%取对数
Fimg=fft2(lnimg);%傅里叶变换
P=fftshift(Fimg);%将频域原点移到图像中心;
[M,N]=size(P);
subplot(2,2,2),imshow(uint8(abs(P))
07-19
831

06-23