利用MATLAB对图像进行DCT变换的代码

利用MATLAB对图像进行DCT变换，并滤掉高频分量，保留低频分量，并查看结果

I=imread（‘sample.jpg‘）; 
figureimshow（I）; 

J=im2double（I）; 
I2=J;
T=dctmtx（8）; 
Mask=[1 1 1 1 1 1 0 0  
      1 1 1 1 1 0 0 0  
      1 1 1 1 0 0 0 0  
      1 1 1 0 0 0 0 0 
      1 1 0 0 0 0 0 0 
      1 0 0 0 0 0 0 0  
      0 0 0 0 0 0 0 0 
      0 0 0 0 0 0 0 0]; 
  for i=1:3
      B=blkproc（J（::i）[8 8]‘P1*x*P2‘TT‘）;
      B=blkproc（B[8 8]‘P1.*x‘Mask）;%此处为点乘
      I2（::i）=blkproc（B[8 8]‘P1*x*P2‘T‘T）; 
  end


figureimshow（I2）;%显示重建后的图像

K=im2double（I）;
I3=K;
Mask1=[0 0 0 0 0 0 0 0  
      0 0 0 0 0 0 0 0  
      0 0 0 0 0 0 0 1  
      0 0 0 0 0 0 1 1 
      0 0 0 0 0 1 1 1 
      0 0 0 0 1 1 1 1  
      0 0 0 1 1 1 1 1 
      0 0 1 1 1 1 1 1]; 
  for i=1:3
      B=blkproc（K（::i）[8 8]‘P1*x*P2‘TT‘）;
      B=blkproc（B[8 8]‘P1.*x‘Mask）;%此处为点乘
      I3（::i）=blkproc（B[8 8]‘P1*x*P2‘T‘T）; 
  end
  
  figureimshow（I3）;%显示重建后的图像

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

     文件     115200  2009-10-09 14:51  DCT变换\作业说明.doc

     文件      80691  2009-05-19 07:14  DCT变换\原图像.jpg

     文件      72608  2009-05-19 07:15  DCT变换\变换后图像.jpg

     文件      28063  2009-05-19 07:18  DCT变换\变换后图像2.jpg

     文件        720  2009-05-19 07:00  DCT变换\源代码.asv

     文件        935  2009-05-19 07:00  DCT变换\源代码.m

     目录          0  2009-10-09 14:51  DCT变换

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

               298217                    7