编程速记(39)：Matlab篇-提取图像高频/低频信息-基于DCT

一、简述

希望通过离散余弦变换获取某张图片的高频/低频信息

二、Demo

转灰图处理方式

RGB = imread('autumn.tif');
I = rgb2gray(RGB);
J = dct2(I);
figure
imshow(log(abs(J)),[])
colormap(gca,jet(64))
colorbar
J(abs(J) < 10) = 0;
K = idct2(J);
figure
imshowpair(I,K,'montage')
title('Original Grayscale Image (Left) and Processed Image (Right)');

上述代码得到的是图片低频信息，如果希望得到高频信息只需要将代码中得到的DCT系数减去低频DCT系数，之后再进行IDCT处理即可。

分通道处理方式

有的时候可能不希望转灰图处理，也可以分通道分别进行DCT处理来提取高频信息。

HIGH = [];
for c=1:3
    I = img(:,:,c);
    J = dct2(I);

    ori_ = J;
    J(abs(J) < threshold) = 0; %inf测试通过
    H_dct = ori_ - J;
    HIGH(:,:,c) = idct2(H_dct);
end

注：上面代码得到的是图像高频信息。

三、更多参考资料

• Remove High Frequencies in Image using DCT