%對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)