2013-07-18 43 views
2

我尝试在matlab中编写一个m文件,在256x256尺寸的图像中执行jpeg压缩,在下面的代码中我将图像分割为8x8子块并执行dct变换并进行量化。我的问题是,我还希望在量化之后将新生成的值保存在矩阵中,以使256x256矩阵也具有0值,然后执行锯齿形操作。我的代码如下:Jpeg量化后的压缩矩阵

clc; 
clear all; 
data=imread('cameraman.tif'); 
x=double(data); 

%length of image with no compression 
len=256*256; 
data1=uint8(zeros(size(data))); 

%quantization weigting table 

    qv1= [ 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]; 

     final=zeros(256,256); 

     count1=0; 
     for row=1:8:256 
      for column = 1:8:256 
      %shift down pixel 
      x=double(data(row:row+7,column:column+7))-128; 
      %8x8 dct transform 
      xf=dct2(x); 
      %quantization 
      xf=round(xf./qv1).*qv1; 
      final(row+7,column+7)=xf; 
      count1=count1+sum(xf(:)~=0); 
      %inverse dct transformation, 
      data1(row:row+7,column:column+7)=idct2(xf)+128; 
      end 
     end 

“最终”变量我想使矩阵我已经在上面描述过了。我有错误,因为我拥有它。我不知道如何更改代码才能正确。

回答

0

似乎不是 final(row+7,column+7) 你应该有 final(row:row+7,column:column+7),否则,你的任务不匹配。 此外,循环内的double似乎不必要,因为您在输入之前已经转换为double