實(shí)瞼八圖像變換壓縮編碼

一，實(shí)驗(yàn)?zāi)康模?/p>

1,掌握變換壓縮編碼的原理與方法

2,了解有損壓縮圖像質(zhì)量評(píng)價(jià)的指標(biāo)

二，實(shí)驗(yàn)條件

1,MATLAB軟件

2,典型的灰度，彩色圖像。

三，實(shí)驗(yàn)原理

1,去除數(shù)據(jù)冗余度可以有效地壓縮數(shù)據(jù)

2,圖像編碼壓縮的主要技術(shù)指標(biāo)：壓縮比，客觀評(píng)價(jià)SNR,主觀評(píng)價(jià)

四，實(shí)驗(yàn)內(nèi)容

h以DCT變換為例，采用標(biāo)準(zhǔn)亮度量化表，得有損壓縮的變換值，計(jì)算并顯示

其直方圖

2,采用標(biāo)準(zhǔn)量化表將上述變換值反量化回去后用IDCT變換回圖像，比擬并計(jì)算

壓縮前后圖像的差異(MSE,PSNR等)

3,根據(jù)由1計(jì)算得到的直方圖試判斷采用Huffman編碼的平均碼長(zhǎng)。

五，實(shí)驗(yàn)步驟

1,以DCT變換為例，采用標(biāo)準(zhǔn)亮度量化表，得有損壓縮的變換值，計(jì)算并顯示

其直方圖

(1)代碼如下所示

I=imread('a.jpg*);

J=rgb2gray(I);

J=imresize(J,[512512]);

%imshow(J);

A=zeros(8,8);

B=zeros(512,512);

lighttable=...%標(biāo)準(zhǔn)亮度量化表

(16II101624405161,

1212141926586055,

1413162440576956,

1417222951878062,

182237566810910377,

243555648110411392,

4964788710312112()101,

72929598112100103991;

fori=l:64%用i,j來控制子塊

forj=1:64

form=(8*(i-l)+l):8*i

forn=(8*(j-l)+l):8*j

A(m-(8*(i-1)+1)+1,1)+1)+1)=J(m,n);%提取8*8子塊

end

end

Cl=dct2(A);%對(duì)子塊進(jìn)行dct2變換

C2=round(C1./lighttablc);%將Cl進(jìn)行量化，并取整

forp=l:8%將子塊按序存進(jìn)512*512的矩陣中

forq=1:8

B(p+8*(i-1),q+8*(j-1))=C2(p,q);

end

end

end

end

H=zeros(1,256*2-1);

fori=l:m

forj=l:n

H(B(i,j)+255)=H(B(i,j)+255)+l;%統(tǒng)計(jì)各值概率

end

end

X=l-255:255];

plot(X,H);%繪制坐標(biāo)為-255—255的直方圖

(2)直方圖如下所示

x10

0.5

0

-300-200100200300

2,采用標(biāo)準(zhǔn)量化表將上述變換值反量化回去后用1DCT變換回圖像，比擬并il算

壓縮前后圖像的差異(MSE,PSNR等)

(1)代碼如下所示

fori=l:64%用i,j來控制子塊

forj=1:64

form=(8*(i-l)+l):8*i

forn=(8*(j-l)+l):8*j

A(m-(8*(i-l)+l)+l,n-(8*(j-l)+l)+l)=B(m,n);

end

end

C3=A.*lighttable;%將子塊反量化回去

C4=idct2(C3);%對(duì)C3進(jìn)行IDCT變換

forp=l:8

forq=1:8

D(p+8*(i-1),q+8*(j-1))=C4(p,q);

end

end

end

end

D=uint8(D);

subplot(1,2,1),imshow(J);

subplot(1,2,2),imshow(D);

%計(jì)算J和D均方根誤差

e=double(J)-double(D);

[m,n]=size(e);

MSE=sqrt(sum(sum(e.A2))/(m*n));

k=8;

%k為圖像是表示地個(gè)像素點(diǎn)所用的二進(jìn)制位數(shù)，即位深。

fmax=2.Ak-1;

a=fmax.A2;

b=sum(sum(e.A2));

PSNR=10*log(m*n*a/b);

(2)顯示壓縮前后圖像

(3)顯示MSE和PSNR的值

MSE=4.4339

PSNR=81.04

3,根據(jù)由1計(jì)算得到的直方圖試判斷采用Huffman編碼的平均碼長(zhǎng).

答：由