s091.m

简单的数字图像处理程序

% 灰度图像JPEG 压缩  s091
I91=imread('1.jpg');
I910=rgb2gray(I91);
I911=im2double(I910);
T=dctmtx(8);
B=blkproc(I911,[8,8],'P1*x*P2',T,T');
       Mask=[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];
   B2=blkproc(B,[8,8],'P1.*x',Mask);
   I912=blkproc(B2,[8,8],'P1*x*P2',T',T);
   subplot(1,2,1),imshow(I911);
   subplot(1,2,2),imshow(I912)