duoch.m

数字图像处理

kk=imread('ziji1.jpg');
X=rgb2gray(kk);
X=double(X);
[c,s]=wavedec2(X,2,'sym5');
a1=wrcoef2('a',c,s,'sym5',1);
a2=wrcoef2('a',c,s,'sym5',2);
hd2=wrcoef2('h',c,s,'sym5',2);
vd2=wrcoef2('v',c,s,'sym5',2);
dd2=wrcoef2('d',c,s,'sym5',2);
sX=size(X)
sa1=size(a1)
shd2=size(hd2)
 figure(2);
  subplot(3,2,1);
 imshow(X,[]);
title('原始图像');
subplot(3,2,2);
imshow(a1,[]);
title('尺度1的低频图像');
 subplot(3,2,3);
 imshow(a2,[]);
title('尺度2的低频图像');
subplot(3,2,4);
imshow(hd2,[]);
title('尺度2的水平高频图像');
 subplot(3,2,5);
 imshow(vd2,[]);
title('尺度2的垂直高频图像');
subplot(3,2,6);
imshow(dd2,[]);
title('尺度2的斜线高频图像');