小波变换在图像处理中的运用及其matlab实现 - 副本

clear % 清理工作空间

load wbarb; % 装入图像

figure; % 新建窗口

image(X); % 显示图像

colormap(map) % 设置色彩索引图

title('原始图像'); % 设置图像标题

axis square % 设置显示比例

disp('压缩前图像X的大小'); % 显示文字

whos('X') % 显示图像属性

%对图像用小波进行层小波分解

[c,s]=wavedec2(X,2,'bior3.7'); %提取小波分解结构中的一层的低频系数和高频系数cal=appcoef2(c,s,'bior3.7',1); %水平方向

ch1=detcoef2('h',c,s,1); %垂直方向

cv1=detcoef2('v',c,s,1); %斜线方向

cd1=detcoef2('d',c,s,1);

%各频率成份重构

a1=wrcoef2('a',c,s,'bior3.7',1);

h1=wrcoef2('h',c,s,'bior3.7',1);

v1=wrcoef2('v',c,s,'bior3.7',1);

d1=wrcoef2('d',c,s,'bior3.7',1);

c1=[a1,h1;v1,d1];

%显示分频信息

figure; % 新建窗口

image(c1); % 显示图像

colormap(jet) % 设置色彩索引图

axis square; % 设置显示比例

title ('分解后低频和高频信息'); % 设置图像标题

ca1=appcoef2(c,s,'bior3.7',1);

ca1=wcodemat(ca1,440,'mat',0);

%改变图像高度并显示

ca1=0.5*ca1;

figure; % 新建窗口

image(ca1); % 显示图像

colormap(map); % 设置色彩索引图

axis square; % 设置显示比例

title('第一次压缩图像'); % 设置图像标题

disp('第一次压缩图像的大小为：'); % 显示文字

whos('ca1') % 显示图像属性

%保留小波分解第二层低频信息进行压缩

ca2=appcoef2(c,s,'bior3.7',2);

%首先对第二层信息进行量化编码

ca2=wcodemat(ca2,440,'mat',0);

%改变图像高度并显示

ca2=0.25*ca2;

figure; % 新建窗口

image(ca2); % 显示图像

colormap(map); % 设置色彩索引图axis square; % 设置显示比例

title('第二次压缩图像'); % 设置图像标题disp('第二次压缩图像的大小为：'); % 显示文字whos('ca2') % 显示图像属性

clear; % 清理工作空间

load wbarb; % 装入图形信号

figure; % 新建窗口

image(X); % 显示图像

colormap(map); % 设置色彩索引图

axis square; % 设置显示比例

title('原始图像'); % 设置图像标题

disp('压缩前图像的大小'); % 显示文字

whos('X'); % 显示图像属性

axis square; % 设置显示比例

%对图像用db3小波进行二层小波分解

[c,s]=wavedec2(X,5,'db3');

[thr,sorh,keepapp]=ddencmp('cmp','wv',X);

[Xcomp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,s,'db3',5,thr,sorh,keepapp); figure; % 新建窗口

image(Xcomp); % 显示图像

colormap(map); % 设置色彩索引图

axis square; % 设置显示比例

title(' 压缩后的图像'); % 设置图像标题

disp('压缩后图像的大小'); % 显示文字

whos('Xcomp') % 显示图像属性

disp('小波分解系数中值为0的系数个数百分比'); % 显示文字

disp(perf0); % 显示数值

disp('压缩后剩余能量百分比'); % 显示文字

disp(perfl2); % 显示数值

load wbarb; % 装载原始图像

figure; % 新建窗口

image(X); % 显示图像

colormap(map); % 设置色彩索引图

title('原图'); % 设置图像标题

axis square; % 设置显示比例

%生成含噪图像

init=2055615866; % 初始值

randn('seed',init) % 随机值

x=X+38*randn(size(X)); % 添加随机噪声

figure; % 新建窗口

image(x); % 显示图像

colormap(map); % 设置色彩索引图

title('含噪声图像'); % 设置图像标题

axis square; % 设置显示比例

% 下面对图像消噪，首先利用sym4函数对图像进行2层分解[c,s]=wavedec2(x,2,'sym4'); % 分解

a1=wrcoef2('a',c,s,'sym4',1); % 图像第一层的重构逼近信号　figure; % 新建窗口

image(a1); % 显示图像

colormap(map); % 设置色彩索引图

title('第一层重构图像'); % 设置图像标题

axis square; % 设置显示比例

a2=wrcoef2('a',c,s,'sym4',2); % 图像第二层的重构逼近信号figure; % 新建窗口

image(a2); % 显示图像

colormap(map); % 设置色彩索引图

title('第二层重构图像'); % 设置图像标题

axis square; % 设置显示比例

load wbarb; % 装载原始图像

figure; % 新建窗口

image(X); % 显示图像

colormap(map); % 设置色彩索引图

title('原始图像'); % 设置图像标题

axis square; % 设置显示比例,生成含噪图像并图示

init=2055615866; % 初始值

randn('seed',init); % 随机值

XX=X+8*randn(size(X)); % 添加随机噪声

figure; % 新建窗口

image(XX); % 显示图像

colormap(map); % 设置色彩索引图

title('含噪图像'); % 设置图像标题

axis square; %用小波函数coif2对图像XX进行2层分解

[c,l]=wavedec2(XX,2,'coif2'); % 分解

n=[1,2]; % 设置尺度向量

p=[10.28,24.08]; % 设置阈值向量,对高频小波系数进行阈值处理%nc=wthcoef2('h',c,l,n,p,'s');

%nc=wthcoef2('v',c,l,n,p,'s');

nc=wthcoef2('d',c,l,n,p,'s');

X1=waverec2(nc,l,'coif2'); % 图像的二维小波重构

figure; % 新建窗口

image(X1); % 显示图像

colormap(map); % 设置色彩索引图

title('第一次消噪后的图像'); % 设置图像标题

axis square; % 设置显示比例,再次对高频小波系数进行阈值处理%mc=wthcoef2('h',nc,l,n,p,'s');

mc=wthcoef2('v',nc,l,n,p,'s');

%mc=wthcoef2('d',nc,l,n,p,'s');

X2=waverec2(mc,l,'coif2'); % 图像的二维小波重构

figure; % 新建窗口

image(X2); % 显示图像

colormap(map); % 设置色彩索引图

title('第二次消噪后的图像'); % 设置图像标题

axis square; % 设置显示比例

load sinsin

figure; % 新建窗口

image(X); % 显示图像

colormap(map); % 设置色彩索引图

title('原始图像'); % 设置图像标题

axis square; % 设置显示比例,生成噪声图像

init=2055615866; % 初始值

randn('state',init); % 随机值

x = X + 18*randn(size(X));

figure; % 新建窗口

image(x); % 显示图像

colormap(map); % 设置色彩索引图

title('含噪图像'); % 设置图像标题

axis square; % 设置显示比例,查找默认值 (参见函数ddencmp).

[thr,sorh,keepapp] = ddencmp('den','wv',x); % 使用全局阈值选项进行图像降噪处理。xd = wdencmp('gbl',x,'sym4',2,thr,sorh,keepapp);

figure; % 新建窗口

image(xd); % 显示图像

colormap(map); % 设置色彩索引图

title('消噪后的图像'); % 设置图像标题

axis square; % 设置显示比例

load wbarb; % 装载原始图像

figure; % 新建窗口

image(X); % 显示图像

colormap(map); % 设置色彩索引图

title('原始图像'); % 设置图像标题

axis square; % 设置显示比例

[c,s]=wavedec2(X,2,'sym4'); % 进行二层小波分解

sizec=size(c); % 处理分解系数，突出轮廓，弱化细节for I =1:sizec(2) % 小波系数处理

if(c(I)>350)

c(I)=2*c(I);

else

c(I)=0.5*c(I);

end

end

xx=waverec2(c,s,'sym4'); % 分解系数重构

figure; % 新建窗口

image(xx); % 显示图像

colormap(map); % 设置色彩索引图

title('增强图像') % 设置图像标题

axis square; % 设置显示比例

clear % 清理工作空间

load woman; % 装入原图像

X1=X; % 复制

map1=map; % 复制

figure; % 新建窗口

image(X1); % 显示图像colormap(map1); % 设置色彩索引图title('woman'); % 设置图像标题 axis square % 画出woman 图像load wbarb; % 装入原图像

X2=X; % 复制

map2=map; % 复制

for I =1:256

for j=1:256

if(X2(I, j)>100)

X2(I, j)=1.2*X2(I, j);

else

X2(I, j)=0.5*X2(I, j);

end

end

end

figure; % 新建窗口

image(X2); % 显示图像colormap(map2); % 设置色彩索引图title('wbarb'); % 设置图像标题axis square % 设置显示比例

[c1,s1]=wavedec2(X1,2,'sym4'); % 分解sizec1=size(c1);

for I=1:sizec1(2) % 系数处理

c1(I)=1.2*c1(I);

end

[c2,s2]=wavedec2(X2,2,'sym4'); % 分解

c=c1+c2; % 系数处理

c=0.5*c;

xx=waverec2(c,s1,'sym4'); % 重构figure; % 新建窗口

image(xx); % 显示图像colormap(map); % 设置色彩索引图title('融合图像'); % 设置图像标题axis square % 设置显示比例

load cathe_1; % 调入第一幅模糊图像

X1=X; % 复制

load cathe_2; % 调入第二幅模糊图像

X2=X; % 复制

XFUS=wfusimg(X1,X2,'sym4',5,'max','max'); %基于小波分解的图像融合figure; % 新建窗口

image(X1); % 显示图像

colormap(map); % 设置色彩索引图

axis square; % 设置显示比例

title(' Catherine 1'); % 设置图像标题

figure; % 新建窗口

image(X2); % 显示图像

colormap(map); % 设置色彩索引图

axis square; % 设置显示比例

title(' Catherine 2'); % 设置图像标题

figure; % 新建窗口

image(XFUS); % 显示图像

colormap(map); % 设置色彩索引图

axis square; % 设置显示比例

title('Synthesized image'); % 设置图像标题