matlab中图像压缩

原址：https://blog.csdn.net/silence2015/article/details/69063738

意义：图像视频只有被压缩才能有效大规模的存储和发送。

那么这儿总结我自己理解的图像压缩简单流程

压缩的方法论：我们首先做图像压缩是在频率域处理的，通过DCT（离散余弦变换）将图像转到频率域。低频部分也存储了图像的大多信息。我们知道，低频部分集中较多能量，含有图像大多平滑信息，而高频部分主要是边缘或者噪声。人眼对低频的光波比较敏感，故我们将高频部分合理丢掉部分，然后将频率域的图像进行量化处理，量化后的频率图像再进行编码处理，比如用哈夫曼编码来构造最短编码。通过编码后的图像就占用很少的一部分空间了。 
如需恢复图像，只需要指定相应的解码器来解码，再

压缩过程： 
1、将图像分块，分成8*8 像素的小块来分别计算。 
2、对每个小块进行理算余弦变换，将图像转换到频率域。 
3、将图像量化，可通过量化矩阵（或者简单将每个像素除以N，取整，然后再乘以N，得到量化的目的，当然N越大，压缩的比率也就越大） 
4、通过合适的编码规则来对量化后的图像编码

解压过程： 
1、通过解码器将图像解码 
2、通过反变换（离散余弦反变换） 将图像从频率域转回空间域。 
3、合并所有8*8的小块

实验总结： 
1、为什么用DCT,而不用FFT？ 
因为DCT可以说少计算了复数，更方便计算

2、压缩损失主要在哪儿？ 
主要在量化的过程，图像做DCT变换只不过是一个从空间域到频率域的变换，并没有改变图像的属性。而量化的过程就是一个不断取整，保留大头，舍弃小的的过程。

function compressIMG(imgpath)
    f=imread(imgpath);
    %f=im2double(f);
    f=double(f);
    T=dctmtx(8);
    dct=@(block_struct)T*block_struct.data*T';
    invdct=@(block_struct)T'*block_struct.data*T;

    f_tf=blockproc(f,[8,8],dct);
    qt_mtx=[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];

     % quantization
     f_qt=blockproc(f_tf,[8,8],@(block_struct)block_struct.data./qt_mtx);
     f_qt=ceil(f_qt);

     % restore the image
     g=blockproc(f_qt,[8,8],@(block_struct)block_struct.data.*qt_mtx);
     g=blockproc(g,[8,8],invdct);
     error=f-g;
     f=uint8(f);
     g=uint8(g);
     error=uint8(error);
     imshow(f),title('original image');
     figure;
     imshow(g),title('compressed image');
     imwrite(g,'compressed_img.jpg');
     figure;
     imshow(error),title('error image');
end

• 
  

上面代码是我做的压缩实验，可以看出来效果并不理想。先往后学，以后再来提高。

参考资料：wiki-JPEG, 图像为何转到频率域处理

最后致谢 Duke University

————————————–更新————————————–

经网友提示我复原得到的图的问题在于量化数值精度处理的问题，源代码中我用的f_qt=ceil(f_qt),ceil为向上取整，量化会造成数据丢失，使用round四舍五入可以完美解决这个问题f_qt=round(f_qt)。