s092.m

简单的数字图像处理程序

% 灰度图像JPEG 压缩倍数调整  s092
I91=imread('1.jpg');
I910=rgb2gray(I91);
I911=im2double(I910);
T=dctmtx(8);
B=blkproc(I911,[8,8],'P1*x*P2',T,T');
       Mask1=[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];
      Mask2=[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
             0 0 0 0 0 0 0 0
             0 0 0 0 0 0 0 0];
   B2=blkproc(B,[8,8],'P1.*x',Mask1);
   B3=blkproc(B,[8,8],'P1.*x',Mask2);
   I912=blkproc(B2,[8,8],'P1*x*P2',T',T);
   I913=blkproc(B3,[8,8],'P1*x*P2',T',T);
   subplot(1,3,1),imshow(I911);
   subplot(1,3,2),imshow(I912);
   subplot(1,3,3),imshow(I913)