matlab彩色图像同态滤波

同态滤波是一种可以修复图像的光照和对比度的滤波方法。在 Matlab 中，可以使用 `imfilter` 函数进行图像同态滤波。

下面是一个示例代码，可以对 RGB 彩色图像进行同态滤波：

% 读入彩色图像
rgb = imread('example.jpg');

% 将 RGB 图像转换为 YCbCr 颜色空间
ycbcr = rgb2ycbcr(rgb);

% 对 Y 分量进行同态滤波
Y = double(ycbcr(:,:,1));
Yh = homfilt(Y);  % 使用自定义的同态滤波函数

% 将处理后的 Y 分量重新组合为 YCbCr 图像
ycbcrh = ycbcr;
ycbcrh(:,:,1) = uint8(Yh);

% 将 YCbCr 图像转换回 RGB 颜色空间
rgbh = ycbcr2rgb(ycbcrh);

% 显示原图和处理后的图像
figure; imshow(rgb);
figure; imshow(rgbh);

其中，自定义的同态滤波函数 `homfilt` 可以按照以下方式实现：

function H = homfilt(I)

% 图像大小
[M,N] = size(I);

% 频域中心坐标
cx = round(M/2);
cy = round(N/2);

% 构造高斯滤波器
sigma = 10;
[X,Y] = meshgrid(1:N,1:M);
H = exp(-((X-cx).^2+(Y-cy).^2)/(2*sigma^2));

% 将滤波器中心移动到图像左上角
H = ifftshift(H);

% 将图像转换到频域
Ih = fft2(double(I));

% 对频域图像进行滤波
Ih = Ih .* H;

% 将滤波后的图像转换回空域
H = real(ifft2(Ih));
end

需要注意的是，同态滤波常常需要对滤波器的参数进行调整，以适应不同的图像和应用场景。在上面的代码中，我们使用了一个固定的高斯滤波器，可以根据需要进行调整。