資源簡介
編寫一個Matlab程序，實現基于DCT的圖像變換編碼。
編碼過程：將原始圖像劃分成8×8的圖像塊，采用離散余弦變換（DCT）對圖像塊進行變換。
解碼過程：對于每個圖像塊，使用一定比例的最低頻DCT系數（被舍棄的高頻DCT系數設為0）做逆離散余弦變換（IDCT），得到重構的圖像塊。最后，將所有重構的圖像塊按順序拼接成完整的解碼圖像。
分別取32、16、8個最低頻DCT系數（如圖1所示）進行反變換得到重構的圖像，比較圖像的質量，計算這三種情況下的峰值信噪比。
峰值信噪比的計算公式：PSNR = 10 log10 (2552/ MSE)
其中，MSE（Mean Squared Error）指原始圖
代碼片段和文件信息
%程序思路：
%對圖像進行DCT變化有2個函數，dct（）和dctmtx（）
%為了對圖像分塊處理，使用dctmtx函數產生一個8*8的DCT變換矩陣T
%分塊處理使用blkproc函數，利用T*x*T‘對每個8*8圖像塊x分別進行DCT變換，T‘*x*T進行IDCT變換
%進行離散余弦逆變換時，分別取32、16、8個最低頻DCT系數矩陣MASK1MASK2MASK3
%使用blkpro函數分塊舍棄高頻DCT系數，利用M.*x
I=imread（‘lena.jpg‘‘jpg‘）;%讀入圖片
I_double=im2double（I）;%將圖像轉換為雙精度格式
T=dctmtx（8）;%生成一個8×8的DCT變換矩陣T
I_dct=blkproc（I_double[88]‘P1*x*P2‘TT‘）;%對每個圖像塊分別進行DCT變換?
MASK1=[?1?1?1?1?1?1?1?1?????
????????1?1?1?1?1?1?1?0?
????????1?1?1?1?1?1?0?0?
????????1?1?1?1?1?0?0?0???
????????1?1?1?0?0?0?0?0????
????????1?1?0?0?0?0?0?0????
????????1?0?0?0?0?0?0?0????
????????0?0?0?0?0?0?0?0];?
MASK2=[?1?1?1?1?1?1?0?0?????
????????1?1?1?1?0?0?0?0?
????????1?1?1?0?0?0?0?0?
????????1?1?0?0?0?0?0?0???
????????1?0?0?0?0?0?0?0????
????????0?0?0?0?0?0?0?0????
????????0?0?0?0??屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2521??2018-06-21?09:26??code.m
評論
共有 條評論