10-6.m

matlab图像处理源码

%对autumn.tif文件计算二维DCT变换
RGB = imread('autumn.tif');
I = rgb2gray(RGB);
%真彩色图像转换成灰度图像
J = dct2(I);
%计算二维DCT变换
imshow(log(abs(J)),[])
%图像大部分能量集中在上左角处
colormap(jet(64)), colorbar
figure;
J(abs(J) < 10) = 0;
%把变换矩阵中小于10的值置换为0，然后用idct2重构图像
K = idct2(J)/255;
imshow(K)