小波图像处理

小波图像处理 找一张图像，用小波变换做图像压缩，（设置阈值的方法，把阈值内的值置0），压缩

clear;
close all;               
  
X=imread（‘3.jpg‘）;              
X=rgb2gray（X）;  
figure;
imshow（X‘InitialMagnification‘‘fit‘）;               
title（‘原始图像‘）;                    
%对图像用小波进行层小波分解  
[cs]=wavedec2（X2‘bior3.7‘）;  
%提取小波分解结构中的一层的低频系数和高频系数  
%ca1=appcoef2（cs‘bior3.7‘1）;  
ch1=detcoef2（‘h‘cs1）;      %水平方向  
cv1=detcoef2（‘v‘cs1）;      %垂直方向  
cd1=detcoef2（‘d‘cs1）;      %斜线方向  
  
%各频率成份重构  
a1=wrcoef2（‘a‘cs‘bior3.7‘1）;  
h1=wrcoef2（‘h‘cs‘bior3.7‘1）;  
v1=wrcoef2（‘v‘cs‘bior3.7‘1）;  
d1=wrcoef2（‘d‘cs‘bior3.7‘1）;  
c1=[a1h1;v1d1];  
%subplot（222）
figure;
imshow（c1[]‘InitialMagnification‘‘fit‘）;  
title （‘分解后低频和高频信息‘）;  
  
%进行图像压缩  
%保留小波分解第一层低频信息  
%首先对第一层信息进行量化编码  
ca1=appcoef2（cs‘bior3.7‘1）;  
ca1=wcodemat（ca1440‘mat‘0）;  
%改变图像高度并显示  
ca1=0.5*ca1;  
%subplot（223）;
figure;
imshow（ca1[]‘InitialMagnification‘‘fit‘）;  
title（‘压缩图像‘）;  
  

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

     文件       1060  2018-04-13 00:32  tuxiangxiaobo.m

     文件     130108  2015-11-07 20:13  3.jpg

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

               131168                    2