区域编码阀值编码.docxVIP

第八章作业：实验作业 7 分别用区域编码和阈值编码方法实现图像压缩，用 8×8DCT 变 换，保留 50％的系数(区域编码保留前 50%个系数，阈值编码保 留 50%的大系数），并对解码图像进行比较。 要求：DCT 要自己实现，不能用 matlab 中的 DCT 函数 (1)区域编码 I=imread(E:\大三下\数字图像处理\MATLAB 图像处理\1.jpg); I=double(rgb2gray(I)); figure imshow(uint8(I)); title(原图); g=zeros(8,8); [m,n]=size(g); a=0; for i=1:m for j=1:m if i==1 a=sqrt(1/m); else a=sqrt(2/m); end q=((2*(j-1)+1)*(i-1)*pi)/(2*m); g(i,j)=a*cos(q); end end I1=blkproc(I,[8 8],P1*x*P2,g,g); figure imshow(uint8(I1)); title(DCT 变换); a=[]; a0=ones(1,8); a1=zeros(1,8); for i=1:8 if i=4 a(i,:)=a0; else a(i,:)=a1; end end; I2=blkproc(I1,[8 8],P1.*x,a); figure imshow(uint8(I1)); title(对 DCT 区域编码); I3=blkproc(I2,[8 8],P1*x*P2,g,g); figure imshow(uint8(I3)); title(对 DCT 区域编码后反 DCT 变换的图); (2)阀值编码 I=imread(E:\大三下\数字图像处理\MATLAB 图像处理\1.jpg); I=double(rgb2gray(I)); figure imshow(uint8(I)); title(原图); g=zeros(8,8); [m,n]=size(g); a=0; for i=1:m for j=1:m if i==1 a=sqrt(1/m); else a=sqrt(2/m); end q=((2*(j-1)+1)*(i-1)*pi)/(2*m); g(i,j)=a*cos(q); end end I1=blkproc(I,[8 8],P1*x*P2,g,g); figure imshow(uint8(I1)); title(DCT 变换); =ones(8,8); =reshape(g,1,64); =median(b); for i=1:8 for j=1:8 if(abs(g(i,j))c) a(i,j)=0; end; end; end; I2=blkproc(I1,[8 8],P1.*x,a); figure imshow(uint8(I2)); title(对 DCT 阀值编码); I3=blkproc(I2,[8 8],P1*x*P2,g,g); figure imshow(uint8(I3)); title(对 DCT 阀值编码后反 DCT 变换的图); 心得体会： 无论区域编码还是阀值编码,对两种编码进行解码后,效果都比较好,视觉上看和原图没什么 差异。 刚开始没弄懂区域编码和阀值编码是怎么一回事，看书？书上讲得太粗糙，最后自己上网 查阅。在弄懂了区域编码和阀值编码后，加上之前掌握的关于 DCT 变换的知识，再来做就 先得不是那么棘手了。