一、简述
希望通过离散余弦变换获取某张图片的高频/低频信息
二、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
注:上面代码得到的是图像高频信息。