test_dct.asv

Compress an decomprres DCT=Domain Transform Coding for image.

clear;clc;

%  Mask = [1/100 1 1 1 1 1 1 1;
%         1 1 1 1 1 1 1 1;
%         1 1 1 1 1 1 1 1;
%         1 1 1 1 1 1 1 1;
%         1 1 1 1 1 1 1 1;
%         1 1 1 1 1 1 1 1;
%         1 1 1 1 1 1 1 1;
%         1 1 1 1 1 1 1 1];
% 
%  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 = [ 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 = [ 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 = [  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_18_Masuda1.tif';
DCT_compr(filename, n_b, Q10M)