wavelet.m

此源代码介绍了几种常用的以小波域理论为基础通过matlab实现的代码

%local distirbution-----wavelet coeficients
%% bylee 2004 4 06
I0=imread('cameraman.tif');   
I=double(Io);
wname='db8';
mode('per');
%% 小波分解
[I11 IH1 IV1 ID1]=dwt2(I,wname);
[I21 IH2 IV2 ID2]=dwt2(I11,wname);
[I31 IH3 IV3 ID3]=dwt2(I21,wname);
window=[7 7];
opname='std2';
%%计算方差
sigmI11=nlfilter(I11,window,opname);
sigmIH1=nlfilter(IH1,window,opname);
sigmIV1=nlfilter(IV1,window,opname);
sigmID1=nlfilter(ID1,window,opname);
sigmI21=nlfilter(I21,window,opname);
sigmIH2=nlfilter(IH2,window,opname);
sigmIV2=nlfilter(IV2,window,opname);
sigmID2=nlfilter(ID2,window,opname);
sigmI31=nlfilter(I31,window,opname);
sigmIH3=nlfilter(IH3,window,opname);
sigmIV3=nlfilter(IV3,window,opname);
sigmID3=nlfilter(ID3,window,opname);
%%均值
opname='mean2';
meanI11=nlfilter(I11,window,opname);
meanIH1=nlfilter(IH1,window,opname);
meanIV1=nlfilter(IV1,window,opname);
meanID1=nlfilter(ID1,window,opname);
meanI21=nlfilter(I21,window,opname);
meanIH2=nlfilter(IH2,window,opname);
meanIV2=nlfilter(IV2,window,opname);
meanID2=nlfilter(ID2,window,opname);
meanI31=nlfilter(I31,window,opname);
meanIH3=nlfilter(IH3,window,opname);
meanIV3=nlfilter(IV3,window,opname);
meanID3=nlfilter(ID3,window,opname);
IsI11=I11./(sigmI11+1^10);
IsIH1=IH1./(sigmIH1+1^10);
IsIV1=IV1./(sigmIV1+1^10);
IsID1=ID1./(sigmID1+1^10);
IsI21=I21./(sigmI21+1^10);
IsIH2=IH2./(sigmIH2+1^10);
IsIV2=IV2./(sigmIV2+1^10);
IsID2=ID2./(sigmID2+1^10);
IsI31=I31./(sigmI31+1^10);
IsIH3=IH3./(sigmIH3+1^10);
IsIV3=IV3./(sigmIV3+1^10);
IsID3=ID3./(sigmID3+1^10);
S1=size(IsI11);
Iss1=reshape(IsI11,S1(1)*S1(2),1);
Iss2=reshape(IsIH1,S1(1)*S1(2),1);
Iss3=reshape(IsIV1,S1(1)*S1(2),1);
Iss4=reshape(IsID1,S1(1)*S1(2),1);
S2=size(IsI21);
Iss21=reshape(IsI21,S2(1)*S2(2),1);
Iss22=reshape(IsIH2,S2(1)*S2(2),1);
Iss23=reshape(IsIV2,S2(1)*S2(2),1);
Iss24=reshape(IsID2,S2(1)*S2(2),1);
S3=size(IsI31);
Iss31=reshape(IsI31,S3(1)*S3(2),1);
Iss32=reshape(IsIH3,S3(1)*S3(2),1);
Iss33=reshape(IsIV3,S3(1)*S3(2),1);
Iss34=reshape(IsID3,S3(1)*S3(2),1);
%%显示分布图
figure(4);
subplot(2,2,1);
histfit(Iss1);xlabel('I11');
subplot(2,2,2);
histfit(Iss2);xlabel('IH1');
subplot(2,2,3);
histfit(Iss3);xlabel('IV1');
subplot(2,2,4);
histfit(Iss4);xlabel('ID1');
figure(5);
subplot(2,2,1);
histfit(Iss21);xlabel('I21');
subplot(2,2,2);
histfit(Iss22);xlabel('IH2');
subplot(2,2,3);
histfit(Iss23);xlabel('IV2');
subplot(2,2,4);
histfit(Iss24);xlabel('ID2');
figure(6);
subplot(2,2,1);
histfit(Iss31);xlabel('I31');
subplot(2,2,2);
histfit(Iss32);xlabel('IH3');
subplot(2,2,3);
histfit(Iss33);xlabel('IV3');
subplot(2,2,4);
histfit(Iss34);xlabel('ID3');