實驗圖像壓縮的MATLAB實現(xiàn)

1、實驗十實驗十 圖像壓縮的圖像壓縮的MATLABMATLAB實現(xiàn)實現(xiàn) 一、實驗?zāi)康囊弧嶒災(zāi)康膒了解有關(guān)數(shù)字圖像壓縮的基本概念，熟悉了解有關(guān)數(shù)字圖像壓縮的基本概念，熟悉MATLAB軟件中關(guān)于數(shù)字圖像壓縮的基本方法，掌握利用軟件中關(guān)于數(shù)字圖像壓縮的基本方法，掌握利用MATLAB軟件進行數(shù)字圖像壓縮的方法。軟件進行數(shù)字圖像壓縮的方法。 二、相關(guān)知識二、相關(guān)知識p在當(dāng)今的信息時代，圖像在表達各種信息時有著不在當(dāng)今的信息時代，圖像在表達各種信息時有著不可替代的作用，但圖像信息的缺點之一就是數(shù)據(jù)量可替代的作用，但圖像信息的缺點之一就是數(shù)據(jù)量非常龐大，因此，無論是存儲還是傳輸，都需要對非常龐大，因此，無論是

2、存儲還是傳輸，都需要對圖像數(shù)據(jù)進行壓縮，數(shù)據(jù)壓縮的方法有很多，我們圖像數(shù)據(jù)進行壓縮，數(shù)據(jù)壓縮的方法有很多，我們這里介紹一種基于這里介紹一種基于DCT（離散余弦變換）的圖像壓（離散余弦變換）的圖像壓縮方法，并介紹用縮方法，并介紹用MATLAB軟件來實現(xiàn)這個算法。軟件來實現(xiàn)這個算法。p基于基于DCT的壓縮方法如下：的壓縮方法如下：p（1）首先將輸入圖像分解為首先將輸入圖像分解為88或或1616的塊，的塊，然后對每個塊進行二維然后對每個塊進行二維DCT變換，這里，一個變換，這里，一個NN圖像塊圖像塊f(x,y)的二維離散余弦變換公式如下：的二維離散余弦變換公式如下： 二、相關(guān)知識二、相關(guān)知識p二維離

3、散余弦反變換公式如下：二維離散余弦反變換公式如下：p基于基于DCT的壓縮方法如下：的壓縮方法如下：p（1）首先將輸入圖像分解為首先將輸入圖像分解為88或或1616的塊，的塊，然后對每個塊進行二維然后對每個塊進行二維DCT變換，這里，一個變換，這里，一個NN圖像塊圖像塊f(x,y)的二維離散余弦變換公式如下：的二維離散余弦變換公式如下： 1100(21)(21)( , )( ) ( )( , )coscos22NNxyxuyvF u vc u c vf x yNN0,1,1xN0,1,1yNp二維離散余弦反變換公式如下：二維離散余弦反變換公式如下：p其中：其中： 1100(21)(21)( ,

4、)( ) ( )( , )coscos22NNxyxuyvF u vc u c vf x yNN0,1,1xN0,1,1yN1100(21)(21)( , )( ) ( ) ( , )coscos22NNuvxuyvf x yc u c v F u vNN0,1,1xN0,1,1yN10( )21,2,uNc uuNN10( )21,2,1vNc vvNNp在在MATLAB中中，p稱為稱為DCT的變換核的變換核p其中：其中： 1100(21)(21)( , )( ) ( ) ( , )coscos22NNuvxuyvf x yc u c v F u vNN0,1,1xN0,1,1yN10( )

5、21,2,uNc uuNN10( )21,2,1vNc vvNN(21)(21)( ) ( ) ( , )coscos22xuyvc u c v F u vNNp在在MATLAB中中，p稱為稱為DCT的變換核的變換核p pMATLAB圖像處理工具箱提供了一些函數(shù)進行圖像處理工具箱提供了一些函數(shù)進行DCT變換。變換。p函數(shù)函數(shù)dct2實現(xiàn)圖像的二維離散余弦變換，格式為：實現(xiàn)圖像的二維離散余弦變換，格式為：p B=dct2(A)p B=dct2(A,M,N)p B=dct2(A,M,N) (21)(21)( ) ( ) ( , )coscos22xuyvc u c v F u vNNpA表示要變換

6、的圖像，表示要變換的圖像，B表示變換后得到的變換系表示變換后得到的變換系數(shù)矩陣，數(shù)矩陣，B和和A是同樣大小的矩陣，其內(nèi)容是余弦是同樣大小的矩陣，其內(nèi)容是余弦變換后的系數(shù)。變換后的系數(shù)。pM和和N是可選參數(shù)，表示對圖像矩陣是可選參數(shù)，表示對圖像矩陣A的填充或截取。的填充或截取。p函數(shù)函數(shù)idct2實現(xiàn)圖像的二維離散余弦反變換，語法實現(xiàn)圖像的二維離散余弦反變換，語法變換。變換。p函數(shù)函數(shù)dct2實現(xiàn)圖像的二維離散余弦變換，格式為：實現(xiàn)圖像的二維離散余弦變換，格式為：p B=dct2(A)p B=dct2(A,M,N)p B=dct2(A,M,N) pA表示要變換的圖像，表示要變換的圖像，B表示變換

7、后得到的變換系表示變換后得到的變換系數(shù)矩陣，數(shù)矩陣，B和和A是同樣大小的矩陣，其內(nèi)容是余弦是同樣大小的矩陣，其內(nèi)容是余弦變換后的系數(shù)。變換后的系數(shù)。pM和和N是可選參數(shù)，表示對圖像矩陣是可選參數(shù)，表示對圖像矩陣A的填充或截取。的填充或截取。p函數(shù)函數(shù)idct2實現(xiàn)圖像的二維離散余弦反變換，語法實現(xiàn)圖像的二維離散余弦反變換，語法格式：格式：p B=idct2(A)p B=idct2(A,M,N)p B=idct2(A,M,N) pA表示要變換的二維離散余弦變換矩陣，表示要變換的二維離散余弦變換矩陣，B表示變表示變換后得到的圖像，換后得到的圖像，B和和A是同樣大小的矩陣，其內(nèi)是同樣大小的矩陣，其內(nèi)

8、容是余弦變換后的系數(shù)。容是余弦變換后的系數(shù)。pM和和N是可選參數(shù)，表示對圖像矩陣是可選參數(shù)，表示對圖像矩陣A的填充或截取。的填充或截取。p函數(shù)函數(shù)dctmtx用于計算二維用于計算二維DCT矩陣，語法格式：矩陣，語法格式：p D=dctmtx(n)p格式：格式：p B=idct2(A)p B=idct2(A,M,N)p B=idct2(A,M,N) pA表示要變換的二維離散余弦變換矩陣，表示要變換的二維離散余弦變換矩陣，B表示變表示變換后得到的圖像，換后得到的圖像，B和和A是同樣大小的矩陣，其內(nèi)是同樣大小的矩陣，其內(nèi)容是余弦變換后的系數(shù)。容是余弦變換后的系數(shù)。pM和和N是可選參數(shù)，表示對圖像矩陣

9、是可選參數(shù)，表示對圖像矩陣A的填充或截取。的填充或截取。p函數(shù)函數(shù)dctmtx用于計算二維用于計算二維DCT矩陣，語法格式：矩陣，語法格式：p D=dctmtx(n)p其中其中D是返回的是返回的nn的的DCT變換矩陣，如果矩陣變換矩陣，如果矩陣A的的大小是大小是nn，DA是矩陣每一列的是矩陣每一列的DCT變換值，變換值，AD是是A的每一行的的每一行的DCT變換值。變換值。pdct2（A）的結(jié)果與的結(jié)果與DAD相同，但后者計算速相同，但后者計算速度較快。度較快。 p（2）將變換后得到的量化的將變換后得到的量化的DCT系數(shù)進行編碼和系數(shù)進行編碼和傳送，形成壓縮后的圖像格式。傳送，形成壓縮后的圖像格

10、式。p基于基于DCT的解壓縮方法如下：的解壓縮方法如下：p（1）對每個對每個88或或1616塊進行二維塊進行二維DCT反變換。反變換。p（2）將反變換的矩陣的塊合成一個單一的圖像。將反變換的矩陣的塊合成一個單一的圖像。p例：把輸入圖像例：把輸入圖像cameraman.tif劃分為劃分為88的圖像的圖像塊，計算它們的塊，計算它們的DCT系數(shù)，并且只保留系數(shù)，并且只保留64個個DCT系系數(shù)中的數(shù)中的10個，然后對每個圖像塊利用這個，然后對每個圖像塊利用這10個系數(shù)進個系數(shù)進其中其中D是返回的是返回的nn的的DCT變換矩陣，如果矩陣變換矩陣，如果矩陣A的的大小是大小是nn，DA是矩陣每一列的是矩陣每

11、一列的DCT變換值，變換值，AD是是A的每一行的的每一行的DCT變換值。變換值。pdct2（A）的結(jié)果與的結(jié)果與DAD相同，但后者計算速相同，但后者計算速度較快。度較快。 p（2）將變換后得到的量化的將變換后得到的量化的DCT系數(shù)進行編碼和系數(shù)進行編碼和傳送，形成壓縮后的圖像格式。傳送，形成壓縮后的圖像格式。p基于基于DCT的解壓縮方法如下：的解壓縮方法如下：p（1）對每個對每個88或或1616塊進行二維塊進行二維DCT反變換。反變換。p（2）將反變換的矩陣的塊合成一個單一的圖像。將反變換的矩陣的塊合成一個單一的圖像。p例：把輸入圖像例：把輸入圖像cameraman.tif劃分為劃分為88的圖

12、像的圖像塊，計算它們的塊，計算它們的DCT系數(shù)，并且只保留系數(shù)，并且只保留64個個DCT系系數(shù)中的數(shù)中的10個，然后對每個圖像塊利用這個，然后對每個圖像塊利用這10個系數(shù)進個系數(shù)進行逆行逆DCT變換來重構(gòu)圖像。變換來重構(gòu)圖像。 p解：程序如下：解：程序如下：pclearpI=imread(cameraman.tif);pI=im2double(I);pT=dctmtx(8);pB=blkproc(I,8 8,P1*x*P2,T,T); %這里T,T是參數(shù)P1、P2的取值pmask=1 1 1 1 0 0 0 0p 1 1 1 0 0 0 0 0 p 1 1 0 0 0 0 0 0p行逆行逆DC

13、T變換來重構(gòu)圖像。變換來重構(gòu)圖像。 p解：程序如下：解：程序如下：pclearpI=imread(cameraman.tif);pI=im2double(I);pT=dctmtx(8);pB=blkproc(I,8 8,P1*x*P2,T,T); %這里T,T是參數(shù)P1、P2的取值pmask=1 1 1 1 0 0 0 0p 1 1 1 0 0 0 0 0 p 1 1 0 0 0 0 0 0p 1 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0;pB2=blkproc(B,8 8

14、,P1.*x,mask); p %這里mask是參數(shù)P1的取值p% 這里可以對B2作進一步的處理，然后加以存儲，解碼時先讀出存儲的數(shù)據(jù)，然后恢復(fù)出B2.pI2=blkproc(B2,8 8,P1*x*P2,T,T);p 1 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0p 0 0 0 0 0 0 0 0;pB2=blkproc(B,8 8,P1.*x,mask); p %這里mask是參數(shù)P1的取值p% 這里可以對B2作進一步的處理，然后加以存儲，解碼時先讀出存儲的數(shù)據(jù)，然后恢復(fù)出B2.pI2=blkproc(B

15、2,8 8,P1*x*P2,T,T);psubplot(1,2,1);pimshow(I);title(原圖);psubplot(1,2,2);pimshow(I2);title(解壓縮圖);psubplot(1,2,1);pimshow(I);title(原圖);psubplot(1,2,2);pimshow(I2);title(解壓縮圖);p 原圖原圖 解壓縮圖解壓縮圖p雖然我們舍棄了雖然我們舍棄了85的的DCT系數(shù)，也就是減少了系數(shù)，也就是減少了85%的存儲量，但可以看到，解壓縮圖仍然清晰。的存儲量，但可以看到，解壓縮圖仍然清晰。在我們現(xiàn)在的實驗中，沒有進行真正的存儲，這部在我們現(xiàn)在的實

16、驗中，沒有進行真正的存儲，這部分需要一些其它的函數(shù)配合工作，留待以后進一步分需要一些其它的函數(shù)配合工作，留待以后進一步完善。完善。 三、實驗內(nèi)容三、實驗內(nèi)容p1分別對圖像分別對圖像“l(fā)enna.tif”、“board.tif”、“peppers.png”進行基于進行基于DCT的壓縮操作，對每幅的壓縮操作，對每幅圖像，分別給出保留圖像，分別給出保留1個、個、2個、個、3個、個、20個個DCT變換系數(shù)的解壓縮結(jié)果，這可以通過調(diào)整矩陣變換系數(shù)的解壓縮結(jié)果，這可以通過調(diào)整矩陣mask中中1的個數(shù)實現(xiàn)，你認(rèn)為保留幾個系數(shù)時，圖的個數(shù)實現(xiàn)，你認(rèn)為保留幾個系數(shù)時，圖像的恢復(fù)效果可以接受，通過觀察，對三個圖像

17、，像的恢復(fù)效果可以接受，通過觀察，對三個圖像，p 原圖原圖 解壓縮圖解壓縮圖p雖然我們舍棄了雖然我們舍棄了85的的DCT系數(shù)，也就是減少了系數(shù)，也就是減少了85%的存儲量，但可以看到，解壓縮圖仍然清晰。的存儲量，但可以看到，解壓縮圖仍然清晰。在我們現(xiàn)在的實驗中，沒有進行真正的存儲，這部在我們現(xiàn)在的實驗中，沒有進行真正的存儲，這部分需要一些其它的函數(shù)配合工作，留待以后進一步分需要一些其它的函數(shù)配合工作，留待以后進一步完善。完善。 三、實驗內(nèi)容三、實驗內(nèi)容p1分別對圖像分別對圖像“l(fā)enna.tif”、“board.tif”、“peppers.png”進行基于進行基于DCT的壓縮操作，對每幅的壓縮操