xscg.m

来自「数字图像处理」· M 代码 · 共 32 行

kk=imread('ziji4.jpg');
X=rgb2gray(kk);
X=double(X);

[c,s]=wavedec2(X,2,'db4');
siz=s(size(s,1),:);
ca1=appcoef2(c,s,'db4',1);
a1=upcoef2('a',ca1,'db4',1,siz);
chd1=detcoef2('h',c,s,1);
hd1=upcoef2('h',chd1,'db4',1,siz);
cvd1=detcoef2('v',c,s,1);
vd1=upcoef2('v',cvd1,'db4',1,siz);
cdd1=detcoef2('d',c,s,1);
dd1=upcoef2('d',cdd1,'db4',1,siz);
 
 figure(3);
  subplot(321);
 imshow(X,[]);
title('原始图像');
subplot(323);
imshow(a1,[]);
title('尺度1的低频系数重构图像');
 subplot(324);
 imshow(hd1,[]);
title('尺度1的水平高频系数重构图像');
subplot(325);
imshow(vd1,[]);
title('尺度1的垂直高频系数重构图像');
 subplot(326);
 imshow(dd1,[]);
title('尺度1的斜线高频系数重构图像');