test_dct.m

Compress an decomprres DCT=Domain Transform Coding for image.

clear;clc;
% 
 Mask = [1/16 1/11 1/10 1/16 1/24 1/40 1/51 1/61;
        1/12 1/12 1/14 1/19 1/26 1/58 1/60 1/55;
        1/14 1/13 1/16 1/24 1/40 1/57 1/69 1/56;
        1/14 1/17 1/22 1/29 1/51 1/87 1/80 1/62;
        1/18 1/22 1/37 1/56 1/68 1/109 1/103 1/77;
        1/24 1/35 1/55 1/64 1/81 1/104 1/113 1/92;
        1/49 1/64 1/78 1/87 1/103 1/121 1/120 1/101;
        1/72 1/92 1/95 1/98 1/112 1/100 1/103 1/99];
%     %Q10M
%       Mask = [ 80,  60,  50,  80, 120, 200, 255, 255;
%          55,  60,  70,  95, 130, 255, 255, 255;
%          70,  65,  80, 120, 200, 255, 255, 255;
%          70,  85, 110, 145, 255, 255, 255, 255;
%          90, 110, 185, 255, 255, 255, 255, 255;
%         120, 175, 255, 255, 255, 255, 255, 255;
%         245, 255, 255, 255, 255, 255, 255, 255;
%         255, 255, 255, 255, 255, 255, 255, 255];
    %Q50M
% Mask = [ 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];
% Q90M    
% Mask = [  3,   2,   2,   3,   5,   8,  10,  12;
%           2,   2,   3,   4,   5,  12,  12,  11;
%           3,   3,   3,   5,   8,  11,  14,  11;
%           3,   3,   4,   6,  10,  17,  16,  12;
%           4,   4,   7,  11,  14,  22,  21,  15;
%           5,   7,  11,  13,  16,  12,  23,  18;
%          10,  13,  16,  17,  21,  24,  24,  21;
%          14,  18,  19,  20,  22,  20,  20,  20];
       

n_b=8;
filename='im_17_kodim15.tif';
DCT_compr(filename, n_b, Mask)