基于pca的图像压缩算法

matlab上利用主成分分析（PCA）对图像进行压缩的程序，有详细注释，给想学习PCA的人提供参考。

function bit_seq=ACEncoding（xy） 
Z=x;
v0=y;
R=abs（y）;
if R==1;amplen=1;
elseif（R >=   2 & R <=    3）;amplen = 2;
elseif（R >=   4 & R <=    7）;amplen = 3;
elseif（R >=   8 & R <=   15）;amplen = 4;
elseif（R >=  16 & R <=   31）;amplen = 5;
elseif（R >=  32 & R <=   63）;amplen = 6;
elseif（R >=  64 & R <=  127）;amplen = 7;
elseif（R >= 128 & R <=  255）;amplen = 8;
elseif（R >= 256 & R <=  511）;amplen = 9;
elseif（R >= 512 & R <= 1023）;amplen = 10;
end 
if （Z==0 & amplen==1）; codelen=2;accode=0;
elseif  （Z==0 & amplen==2）;codelen=2;accode=1;
elseif  （Z==0 & amplen==3）;codelen=3;accode=4;
elseif  （Z==0 & amplen==4）;codelen=4;accode=11;
elseif  （Z==0 & amplen==5）;codelen=5;accode=26;
elseif  （Z==0 & amplen==6）;codelen=7;accode=120;
elseif  （Z==0 & amplen==7）;codelen=8;accode=248;
elseif  （Z==0 & amplen==8）;codelen=10;accode=1014;
elseif  （Z==0 & amplen==9）;codelen=16;accode=65410;
elseif （Z==0 & amplen==10）;codelen=16;accode=65411;
    
elseif  （Z==1 & amplen==1）; codelen=4;accode=12; 
elseif  （Z==1 & amplen==2）; codelen=5;accode=27;  
elseif  （Z==1 & amplen==3）; codelen=7;accode=121; 
elseif  （Z==1 & amplen==4）; codelen=9;accode=502;  
elseif  （Z==1 & amplen==5）;codelen=11;accode=2038; 
elseif  （Z==1 & amplen==6）;codelen=16;accode=65412;
elseif  （Z==1 & amplen==7）;codelen=16;accode=65413;
elseif  （Z==1 & amplen==8）;codelen=16;accode=65414;
elseif  （Z==1 & amplen==9）;codelen=16;accode=65415;
elseif （Z==1 & amplen==10）;codelen=16;accode=65416;
    
elseif （Z==2 & amplen==1）;codelen=5;accode=28;
elseif （Z==2 & amplen==2）;codelen=8;accode=249;
elseif （Z==2 & amplen==3）;codelen=10;accode=1015;
elseif （Z==2 & amplen==4）;codelen=12;accode=4084;
elseif （Z==2 & amplen==5）;codelen=16;accode=65417;
elseif （Z==2 & amplen==6）;codelen=16;accode=65418;
elseif （Z==2 & amplen==7）;codelen=16;accode=65419;
elseif （Z==2 & amplen==8）;codelen=16;accode=65420;
elseif （Z==2 & amplen==9）;codelen=16;accode=65421;
elseif （Z==2 & amplen==10）;codelen=16;accode=65422;
    
    
elseif （Z==3 & amplen==1）;codelen=6;accode=58;
elseif （Z==3 & amplen==2）;codelen=9;accode=503;
elseif （Z==3 & amplen==3）;codelen=12;accode=4085;
elseif （Z==3 & amplen==4）;codelen=16;accode=65423;
elseif （Z==3 & amplen==5）;codelen=16;accode=65424;
elseif （Z==3 & amplen==6）;codelen=16;accode=65425;
elseif （Z==3 & amplen==7）;codelen=16;accode=65426;
elseif （Z==3 & amplen==8）;codelen=16;accode=65427;
elseif （Z==3 & amplen==9）;codelen=16;accode=65428;
elseif （Z==3 & amplen==10）;codelen=16;accode=65429;
    
    
elseif （Z==4 & amplen==1）;codelen=6;accode=59;
elseif （Z==4 & amplen==2）;codelen=10;accode=1016;
elseif （Z==4 & amplen==3）;codelen=16;accode=65430;
elseif （Z==4 & amplen==4）;codelen=16;accode=65431;
elseif （Z==4 & amplen==5）;codelen=16;accode=65432;
elseif （Z==4 & amplen==6）;codelen=16;accode=65433;
elseif （Z==4 & amplen==7）;codelen=16;accode=65434;
elseif （Z==4 & amplen==8）;codelen=16;accode=65435; 
elseif （Z==4 & amplen==9）;codelen=16;accod

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4200  2003-01-09 05:37  pca\8.tif

     文件       9108  2019-07-13 19:47  pca\ACEncoding.m

     文件       1517  2019-07-13 19:46  pca\DCEncoding.m

     文件     263222  2019-07-13 21:22  pca\lena512.bmp

     文件     262336  2019-07-13 21:22  pca\lena512.mat

     文件     786572  2019-07-13 21:22  pca\lena512color.tiff

     文件        913  2019-07-13 19:38  pca\LenaJPEGbyXING.m

     文件     125388  2003-05-03 13:53  pca\liftingbody.png

     文件        379  2019-07-13 19:21  pca\main.m

     文件        431  2015-01-02 19:13  pca\pcaimage.m

     文件        276  2015-01-02 19:14  pca\pcasample.m

     文件       5447  2019-07-13 21:47  pca\Untitled4.m

     目录          0  2019-07-13 21:44  pca

----------- ---------  ---------- -----  ----

              1459789                    13