多尺度二維離散小波分解wavedec2

對X進行N尺度小波分解 [C,S]=wavedec2(X,N,'wname');

clc,clear all,close all;
load woman;
[c,s]=wavedec2(X,2,'db1');%進行2尺度二維離散小波分解。分解小波函數-db1
[cH1,cV1,cD1]=detcoef2('all',c,s,1);%尺度1的所有方向的高頻系數
[cH2,cV2,cD2]=detcoef2('all',c,s,2);%尺度2的所有方向的高頻系數
cA1=appcoef2(c,s,'db1',1);%尺度1的低頻系數
cA2=appcoef2(c,s,'db1',2);%尺度2的低頻系數
figure,imshow(X,map);title('原圖');
figure;
subplot(1,2,1),imshow(uint8(cA1));axis off;title('尺度1的低頻系數圖像');
subplot(1,2,2),imshow(uint8(cA2));axis off;title('尺度2的低頻系數圖像');
figure;
subplot(2,3,1),imshow(uint8(cH1));axis off;title('尺度1水平方向高頻系數');
subplot(2,3,2),imshow(uint8(cV1));axis off;title('尺度1垂直方向高頻系數');
subplot(2,3,3),imshow(uint8(cD1));axis off;title('尺度1斜線方向高頻系數');
subplot(2,3,4),imshow(uint8(cH2));axis off;title('尺度2水平方向高頻系數');
subplot(2,3,5),imshow(uint8(cV2));axis off;title('尺度2垂直方向高頻系數');
subplot(2,3,6),imshow(uint8(cD2));axis off;title('尺度2斜線方向高頻系數');