【正文】 ******************% x=imread(39。D:\船 .jpg39。)。 figure(1) subplot(2,2,1) imshow(x)。 title(39。原始圖像 39。)。 ycbcrx=rgb2ycbcr(x) a=ycbcrx subplot(2,2,2) imshow(a,[])。 title(39。ycbcr彩色空間 39。)。 yiqx=rgb2ntsc(x) b=yiqx。 subplot(2,2,3) imshow(b,[])。 title(39。ntsc彩色空間 39。) x=double(x)。 R=x(:,:,1)。 G=x(:,:,2)。 B=x(:,:,3)。 Y = *R + *G + *B。 U = *R *G + *B。 V = *R *G *B。 c=cat(3,Y,U,V)。 subplot(2,2,4) imshow(c,[])。 title(39。yuv彩色空間 39。) %*************************************** %亮度量化表 LQT=[16 11 10 16 24 40 51 61。 12 12 14 19 26 58 60 55。 14 13 16 24 40 57 69 56。 14 17 22 29 51 87 80 62。 18 22 37 56 68 109 103 77。 24 35 55 64 81 104 113 92。 49 64 78 87 103 121 120 101。 72 92 95 98 112 100 103 99]。 %色度量化表 CQT=[17 18 24 47 99 99 99 99。 18 21 26 66 99 99 99 99。 24 26 56 99 99 99 99 99。 47 66 99 99 99 99 99 99。 99 99 99 99 99 99 99 99。 99 99 99 99 99 99 99 99。 99 99 99 99 99 99 99 99。 99 99 99 99 99 99 99 99]。 %****************************************** %****************************************** %彩色圖像分塊采樣 D = double(x)。 a1= rgb2ycbcr(D)。 Y = a1(:,:,1)。 Cr = a1(:,:,2)。 Cb = a1(:,:,3)。 [M,N] = size(Y)。 Y_last_Z = 0。 Cr_last_Z = 0。 Cb_last_Z = 0。 Y_jpeg2 = 39。39。 Cr_jpeg2 = 39。39。 Cb_jpeg2 = 39。39。 for i = 0:M/81 for j = 0:N/81 S1 = Y((1+8*i):(8*(i+1)), (1+8*j):(8*(j+1)))。%分割成 8*8處理 S2 = Cr((1+8*i):(8*(i+1)), (1+8*j):(8*(j+1)))。 S3 = Cb((1+8*i):(8*(i+1)), (1+8*j):(8*(j+1)))。 A1 = dct2(S1)。%DCT A2 = dct2(S2)。 A3 = dct2(S3)。 B1 = round(A1./LQT)。%亮度量化 B2 = round(A2./CQT)。%色度量化 B3 = round(A3./CQT)。%色度量化 %亮度量化的逆變換 C4=round(B1.*LQT)。 C5=round(B2.*CQT)。 C6=round(B3.*CQT)。 %DCT逆變換 S4=idct2(C4)。 S5=idct2(C5)。 S6=idct2(C6)。 %分塊重組 Y((1+8*i):(8*(i+1)), (1+8*j):(8*(j+1)))=S4。 Cr((1+8*i):(8*(i+1)), (1+8*j):(8*(j+1)))=S5。 Cb((1+8*i):(8*(i+1)), (1+8*j):(8*(j+1)))=S6。 %********************************************************** %Zigzag掃描 Z1 = Zigzag(B1)。 Z2 = Zigzag(B2)。 Z3 = Zigzag(B3)。 %********************************************************* %Huffman編碼 %***************************** %直流分量編碼（ DC_Huffman) DCH_Y = DC_Huffman(Z1(1),Y_last_Z,1)。 DCH_Cr=DC_Huffman(Z2(1),Cr_last_Z,2)。 DCH_Cb=DC_Huffman(Z3(1),Cb_last_Z,2)。 %交流分量編碼 (AC_HUffman) ACH_Y = AC_Huffman(Z1(2:end),1)。 ACH_Cr = AC_Huffman(Z2(2:end),2)。 ACH_Cb = AC_Huffman(Z3(2:end),2)。 Y_last_Z =Z1(1)。 Cr_last_Z = Z2(1)。 Cb_last_Z = Z3(1)。 Y_jpeg2 = [Y_jpeg2 DCH_Y ACH_Y]。 Cr_jpeg2 = [Cr_jpeg2 DCH_Cr ACH_Cr]。 Cb_jpeg2 = [Cb_jpeg2 DCH_Cb ACH_Cb]。 end end figure(2) subplot(2,3,1)。imshow(A1)。title(39。Y采樣后 DCT39。)。 subplot(2,3,2)。imshow(A1)。title(39。cr采樣后 DCT39。)。 subplot(2,3,3)。imshow(A1)。title(39。cb采樣后 DCT39。)。 subplot(2,3,4)。imshow(B1)。title(39。Y亮 度量化 39。)。 subplot(2,3,5)。imshow(B2)。title(39。cr色度量化 39。)。 subplot(2,3,6)。imshow(B3)。title(39。cb色度量化 39。)。 jpeg2 = [Y_jpeg2,Cr_jpeg2,Cb_jpeg2]。 rate2 = 8*M*N*3/length(jpeg2)。%壓縮比 %**************************************** figure(3)。 subplot(1,3,1)。imshow(Y)。title(39。解壓的 Y分量圖像 39。)。 subplot(1,3,2)。imshow(Cr)。title(39。解壓的 Cr分量圖像 39。)。 subplot(1,3,3)。imshow(Cb)。title(39。解壓的 Cb分量圖像 39。)。 a1(:,:,1)=Y。 a1(:,:,2)=Cr。 a1(:,:,3)=Cb。 figure(4)。imshow(YCbCr)。title(39。解壓的 YCbCr的圖像 39。)。 x1=ycbcr2rgb(YCbCr)。 figure(5) imshow(x1) 附錄 2 AC_huffman function ACH = AC_Huffman(AC,n)。 %亮度 AC系數(shù)碼表 T1={ 39。101039。 39。0039。 39。0139。 39。10039。 39。101139。 39。1101039。 39。111100039。 39。1111100039。 39。111111011039。 39。111111111000001039。 39。111111111000001139。 39。110039。 39。1101139。 39。111100139。 39。11111011039。 39。1111111011039。 39。111111111000010039。 39。111111111000010139。 39。111111111000011039。 39。111111111000011139。 39。111111111000100039。 39。039。 39。1110039。 39。1111100139。 39。111111011139。 39。11111111010039。 39。111111111000100139。 39。111111111000101039。 39。111111111000101139。 39。111111111000110039。 39。111111111000110139。 39。111111111000111039。 39。039。 39。11101039。 39。11111011139。 39。11111111010139。 39。111111111000111139。 39。111111111001000039。 39。111111111001000139。 39。111111111001001039。 39。111111111001001139。 39。111111111001010039。 39。111111111001010139。 39。039。 39。11101139。 39。111111100039。 39。111111111001011039。 39。111111111001011139。 39。111111111001100039。 39。111111111001100139。 39。111111111001101039。 39。111111111001101139。 39。111111111001110039。 39。111111111001110139。 39。039。 39。111101039。 39。1111111011139。 39。111111111001111039。 39。111111111001111139。 39。111111111010000039。 39。111111111010000139。 39。111111111010001039。 39。111111111010001139。 39。111111111010010039。 39。111111111010000139。 39。039。 39。111101139。 39。11111111011039。 39。111111111010011039。 39。111111111010011139。 39。111111111010100039。 39。111111111010100139。 39。111111111010101039。 39。111111111010101139。 39。111111111010110039。 39。111111111010110139。 39。039。 39。1111101039。 39。11111111011139。 39。111111111010111039。 39。111111111010111139。 39。111111111011000039。 39。111111111011000139。 39。111111111011001039。 39。111111111011001139。 39。111111111011010039。 39。111111111011010139。 39。039。 39。11111100039。 39。11111111100000039。 39。111111111011011039。 39。111111111011011139。 39。111111111011100039。 39。111111111011100139。 39。111111111011101039。 39。111111111011101139。 39。111111111011110039。 39。111111111011110139。 39。039。 39。11111100139。 39。111111111011111039。 39。111111111011111139。 39。111111111100000039。 39。111111111100000139。 39。111111111100001039。 39。111111111100001139。 39。111111111100010039。 39。111111111100010139。 39。111111111100011039。 39。039。