%?設置壓縮比，cr=0.5為2:1壓縮；cr=0.1250為8:1壓縮
cr?=?0.125;
initialimage?=?imread（‘lena.bmp‘）;
initialimage?=?double（initialimage）/255;
figure（）;
%subplot（121）;
%imshow（initialimage）;

%對圖像進行DCT變換
t?=?dctmtx（8）;
dctcoe?=?blkproc（initialimage?[8?8]?‘P1*x*P2‘?t?t‘）;

%將DCT變換后的矩陣轉換成列，并按升序排列
coevar?=?im2col（dctcoe?[8?8]?‘distinct‘）;
coe?=?coevar;
[y?ind]?=?sort（coevar）;

[m?n]?=?size（coevar）;

%舍去不重要的系數
snum?=?64-64?*?cr;
for?i?=?1:n
????coe（ind（1:snum）?i）?=?0;
end

%把列變換為二維矩陣
b2?=?col2im（coe?[8?8]?[512?512]?‘distinct‘）;

%逆DCT變換
i2?=?blkproc（b2?[8?8]?‘P1*x*P2‘?t‘?t）;
%subplot（122）;
imshow（i2）;

e?=?double（initialimage）-double（i2）;
[m?n]?=?size（e）;
erms?=?sqrt（sum（e（:）.^2）/（m*n））

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????830464??2011-12-30?11:11??DCT圖像壓縮\DCT.doc

?????文件????????785??2011-12-13?15:45??DCT圖像壓縮\DCT.m

?????文件?????263222??2005-06-16?23:18??DCT圖像壓縮\lena.bmp

?????目錄??????????0??2011-12-30?11:11??DCT圖像壓縮

-----------?---------??----------?-----??----

??????????????1094471????????????????????4