Matlab数字图像处理：DCT变换

原链接：http://jingyan.baidu.com/album/ff42efa96254b1c19e2202e2.html

clear;

RGB=imread('原始图像');

GR=rgb2gray(RGB);%转换成灰度图像

figure(1);imshow(GR);

D=dct2(GR);  %计算DCT

figure(2);imshow(log(abs(D)),[]);

colormap(gray(4));colorbar;

D(abs(D)<0.1)=0;

I=idct2(D)/255;

figure(3);imshow(I)

以上程序段对两幅原始图像进行离散余弦变换，如下图DCTtest1a、DCTtest2a所示，变换的结果如图DCTtestb、DCTtest2b所示。由结果可知，变换后DCT系数能量主要集中在左上角，其余大部分系数接近于零，这说明DCT具有适用于图像压缩的特性。将变换后的DCT系数进行门限操作，将小于一定值得系数归零，然后进行逆DCT运算，得到压缩后的图像，如图DCTtest1c、DCTtest2c所示。

比较变换前后的图像，可以发现视觉效果相差很小，可见压缩的效果比较理想。

以下三张分别是DCTtest1a，DCTtest1b，DCTtest1c。