Compress-Decompress algrithms

Development Platform:

Matlab

jpegencode.m：Code Content

							function y=jpegencode(x,quality)
error(nargchk(1,2,nargin));
if nargin<2
    quality=1;
end
x=double(x)-128;
[xm,xn]=size(x);
t=dctmtx(8);
y=blkproc(x,[8 8],'P1*x*P2',t,t');
m=[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]*quality;
ttt=ones(8,8)*10;
m=ttt*quality;
yy=blkproc(y,[8 8],'round(x./P1)',m);
y=im2col(yy,[8 8],'distinct');
xb=size(y,2);
order=[1 9 2 3 10 17 25 18 11 4 5 12 19 26 33 ...
       41 34 27 20 13 6 7 14 21 28 35 42 49 57 50 ...
       43 36 29 22 15 8 16 23 30 37 44 51 58 59 52 ...
       45 38 31 24 32 39 46 53 60 61 54 47 40 48 55 62 63 56 64];
   y=y(order,:);
   eob=max(x(:))+1;
   num=numel(y)+size(y,2);
   r=zeros(num,1);
   count=0;
   for j=1:xb
       i=max(find(y(:,j)));
       if isempty(i)
           i=0;
       end
       p=count+1;
       q=p+i;
       r(p:q)=[y(1:i,j);eob];
       count=count+i+1;
   end
   r((count+1):end)=[];
   r=r+128;
   y.size=uint16([xm,xn]);
   y.numblocks=uint16(xb);
   y.quality=uint16(quality*100);
   [y.huffman y.info]=huffencode(uint8(r));
   
   
   

						