nengliangzhizuida.rar_能量算子

clear all; close all; x=imread('car1.bmp'); y=imread('car2.bmp'); a=rgb2gray(x); b=rgb2gray(y); a=double(a); b=double(b); [aA,aH,aV,aD]=dwt2(a,'bior2.4'); [bA,bH,bV,bD]=dwt2(b,'bior2.4'); newA=zeros(size(aA)); newH=zeros(size(aH)); newV=zeros(size(aV)); newD=zeros(size(aD)); [m,n]=size(aA); for i=1:m for j=1:n newA(i,j)=(aA(i,j)+bA(i,j))/2; end; end; H_aH=nlfilter(aH,[3 3],@lingyuhe); H_aV=nlfilter(aV,[3 3],@lingyuhe); H_aD=nlfilter(aD,[3 3],@lingyuhe); H_bH=nlfilter(bH,[3 3],@lingyuhe); H_bV=nlfilter(bV,[3 3],@lingyuhe); H_bD=nlfilter(bD,[3 3],@lingyuhe); for i=1:m for j=1:n if H_aH(i,j)>=H_bH(i,j); newH(i,j)=aH(i,j); else newH(i,j)=bH(i,j); end end end for i=1:m for j=1:n if H_aV(i,j)>=H_bV(i,j); newV(i,j)=aV(i,j); else newV(i,j)=bV(i,j); end end end for i=1:m for j=1:n if H_aD(i,j)>=H_bD(i,j); newD(i,j)=aD(i,j); else newD(i,j)=bD(i,j); end end end new=idwt2(newA,newH,newV,newD,'bior2.4'); new=uint8(new); figure; imshow(new); title('领域和取大')