同态滤波matlab程序

以下是一个简单的同态滤波的matlab程序：

%读入图像
img = imread('lena.jpg');
%将图像转换为灰度图像
img = rgb2gray(img);
%将图像转换为double类型
img = im2double(img);

%定义同态滤波的参数
gamma_l = 0.25;
gamma_h = 2.5;
c = 1;
D0 = 50;

%计算傅里叶变换
F = fft2(img);
F = fftshift(F);

%生成同态滤波器
[M, N] = size(F);
H = zeros(M, N);
for i = 1:M
    for j = 1:N
        d = (i-M/2)^2 + (j-N/2)^2;
        H(i,j) = (gamma_h-gamma_l)*(1-exp(-c*d/D0^2))+gamma_l;
    end
end

%进行滤波操作
G = H.*F;
G = ifftshift(G);
g = real(ifft2(G));

%显示结果
imshow(img);
title('原图像');
figure;
imshow(g);
title('同态滤波结果');

在程序中，我们首先读入图像，并将其转换为灰度图像。然后，我们对图像进行傅里叶变换，并生成同态滤波器。接着，我们将同态滤波器与傅里叶变换后的图像进行乘法操作，得到滤波后的傅里叶变换结果。最后，我们进行反傅里叶变换，得到同态滤波后的图像，并将其显示出来。