数字图像处理学习笔记（十八）实现图像的离散余弦变换

实验截图：

• 自带函数实现的离散余弦变换：
  
• 利用公式求出的DCT和IDCT
  

实验代码：

代码（1）：

img=imread('erciyuan.jpg');
img=rgb2gray(img);
figure(1)
%显示原图像
subplot(1,3,1);imshow(img);title('原图')
%计算二维dct变换
img_cos=dct2(img)
subplot(1,3,2);
imshow(log(abs(img_cos)));
title('DCT离散余弦变换')
%把变换后的矩阵中小于10的值置换为0，再用idc2重构图像
img_cos(abs(img_cos)<10)=0;
img_acos=idct2(img_cos)/255
subplot(1,3,3),
imshow(img_acos);
title('IDCT反离散余弦变换（恢复原样）')

代码（2）：

img = imread('erciyuan.jpg');
img = rgb2gray(img);
[m,n] = size(img);
M = zeros(m,m);
N = zeros(n,n);
for i = 0 : m - 1
    for j = 0 : m - 1
        if i == 0            
            M(i + 1,j + 1) = sqrt(1 / m) * cos(((2 * j + 1) * i * pi) / (2 * m));         
        else
            M(i + 1,j + 1) = sqrt(2 / m) * cos(((2 * j + 1) * i * pi) / (2 * m));
        end    
    end
end
for i = 0 : n - 1 
    for j = 0 : n  - 1        
        if i == 0             
            N(i + 1,j + 1) = sqrt(1 / n) * cos(((2 * j + 1) * i * pi) / (2 * n));         
        else
            N(i + 1,j + 1) = sqrt(2 / n) * cos(((2 * j + 1) * i * pi) / (2 * n));
        end
    end
end
DCT = M * double(img) * N';
IDCT = M' * DCT * N;
subplot(1,3,1);
imshow(img);
title('原图');
subplot(1,3,2);
imshow(DCT);
title('DCT');
subplot(1,3,3);
imshow(uint8(IDCT));
title('IDCT');

结果分析：

• 离散余弦变换(DCT for Discrete Cosine Transform)是与傅里叶变换相关的一种变换，它类似于离散傅里叶变换(DFT for Discrete Fourier Transform),但是只使用实数。