数字图像处理 - Matlab - 实验笔记五

第八章 图像压缩

1. 分别用区域编码和阈值编码方法实现图像压缩，用8×8DCT变换，保留50％的系数(区域编码保留前50%个系数，阈值编码保留50%的大系数），并对解码图像进行比较

程序代码：

% 区域编码和阈值编码方法实现图像压缩
clear
X0=imread('pic.jpg');
% 把原图转为灰度图
X=double(rgb2gray(X0));
% 获取图像大小
[width,height]=size(X);
% 显示原灰度图
subplot(2,3,1)
imshow(uint8(X));
title('原灰度图')
% DCT变换矩阵
N=8;
G=zeros(N,N);
for i=1:N
    for j=1:N
        if i==1
            a=1/sqrt(N);
        else
            a=sqrt(2/N);
        end
        G(i,j)=a*cos((i-1)*(2*j-1)*pi/(2*N));
    end
end
% 区域模板
mask1=zeros(N,N);
for i=1:N/2
    for j=1:N
        mask1(i,j)=1;
    end
end
% 阈值模板初始化
mask2=zeros(N,N);
% 初始化
Y1=zeros(width,height);
Y2=zeros(width,height);
C0=zeros(width,height);
C1=zeros(width,height);
C2=zeros(width,height);
% 分块
for u=1:N:widt