A=imread('aaa.tif'); I=rgb2gray(A); trueImage=im2double(I); %轉(zhuǎn)換圖像矩陣為雙精度型 imshow(trueImage); title('原始圖像'); dctm=dctmtx(8); %計(jì)算離散余弦變換 imageDCT=blkproc(trueImage,[8 8],'P1*x*P2',dctm,dctm.'); %對(duì)圖像I 的每個(gè)不同8*8數(shù)據(jù)塊應(yīng)用矩陣式’P1*x*P2’進(jìn)行處理。 M10=[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 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; %二值掩模,用來(lái)壓縮DCT 的系數(shù) newImage2=blkproc(imageDCT,[8 8],'P1*(x.*P2)*P3',dctm',M10,dctm); figure; imshow(newImage2); title('壓縮圖像 M10');
|