拉普拉斯圖像融合程序2

funCtiOnlap_fusion()%LaplacianPyramidfusionmul=imread('images?ms1.png');Pan=imread('images?pan.png');figure(1);imshow(mul);title('MS原始圖像');axisfill;figure(2);imshow(pan);title('Pan原始圖像');axisfill;mul=double(rgb2gray(mul))∕255;Pan=double(rgb2gray(pan))/255;%普拉斯金塔變換參數(shù)mp=1;Zt=4;Cf=1;ar=1;CC=[cfar];Y_lap=fuse_lap(mul,pan,zt,cc,mp);figure(3);imshow(Y_lap);title('lapfusion 后的圖像');axisfill;imwrite(Y_lap,'images\lapfusion后的圖像.jpg','Quality',100);%mainfunCtionendfunctionY=fuselap(M1,M2,zt,ap,mp)%Y=fuselap(M1,M2,zt,ap,mp)imagefusionWithlaplacianPyramid%%M1-inputimageA%M2-inputimageB%Zt-maximumdecompositionIeVel%ap-COeffiCientSeIeCtionhighpass(Seeselc.m)%mp-COeffiCientSeIeCtionbaseimage(Seeselb.m)%%Y-fusedimage%(OIiVerRockinger16.08.99)%CheCkinPUtS

[z1s1]=SiZe(M1);[z2s2]=SiZe(M2);if(z1~=z2)|(s1~=s2)error('InputimagesarenotOfSamesize');end;%definefilterW=[14641]/16;%cellsforSeIeCtedimagesE=cell(1,zt);%loopoverdecompositiOndepth->analysisfori1=1:Zt%calculateandstoreactualimageSiZe[zs]=SiZe(M1);zl(i1)=乙sl(i1)=s;%CheCkifimageexpansionnecessaryif(floor(z∕2)~=z/2),ew(1)=1;else,ew(1)=0;end;if(floor(s/2)~=s/2),ew(2)=1;else,ew(2)=0;end;%PerfOrmexpansionifnecessaryif(any(ew))M1=adb(M1,ew);M2=adb(M2,ew);end;%PerfOrmfilteringG1=conv2(conv2(es2(M1,2),w,'valid'),w','valid');G2=conv2(conv2(es2(M2,2),w,'valid'),w','valid');%decimate,UndecimateandinterpolateM1T=conv2(conv2(es2(undec2(dec2(G1)),2),2*w,'valid'),2*w','valid');M2T=conv2(conv2(es2(undec2(dec2(G2)),2),2*w,'valid'),2*w','valid');%SeIeCtCOeffiCientsandstorethemE(i1)={selc(M1-M1T,M2-M2T,ap)};

%decimateM1=dec2(G1);M2=dec2(G2);end;%SeIeCtbaseCOeffiCientsofIaStdecompostiOnStageM1=SeIb(M1,M2,mp);%loopoverdecompositiOndepth->Synthesisfori1=zt:-1:1%UndecimateandinterpolateM仃=COnv2(conv2(es2(undec2(M1),2),2*w,'valid'),2*w','valid');%addCOeffiCientsM1=M1T+E{i1};%SeIeCtvalidimageregiOnM1=M1(1:zl(i1),1:sl(i1));end;%copyimageY=M1;functionY=es2(X,n)%Y=ES2(X,n)SymmetriCextensionofamatrixonallborders%%X-inputmatrix%n-numberofrows/columnstoextend%%Y-extendedmatrix% (OIiVerRockinger16.08.99)[zs]=size(X);Y =ZerOS(z+2*n,s+2*n);Y(n+1:n+z,n:-1:1) =X(:,2:1:n+1);Y(n+1:n+z,n+1:1:n+s) =X;Y(n+1:n+z,n+s+1:1:s+2*n)=X(:,s-1:-1:S-n);Y(n:-1:1,n+1:s+n) =X(2:1:n+1,:);% mp==3:averageAandB% mp==3:averageAandB%Y(n+z+1:1:z+2*n,n+1:s+n)=X(Z-1:-1:Z-n,:);funCtiOnY=dec2(X);%Y=dec2(X)downsamplingofamatrixby2%% X-inPUtmatrix%% Y-OUtPUtmatrix% (OIiVerRockinger16.08.99)[ab]=SiZe(X);Y=X(1:2:a,1:2:b);functionY=Undec2(X)%Y=Undec2(X)UPSamPIingofamatrixby2%% X-inPUtmatrix%% Y-OUtPUtmatrix% (OIiVerRockinger16.08.99)[zs]=size(X);Y=zeros(2*z,2*s);Y(1:2:2*z,1:2:2*s)=X;fUnctionY=SeIb(M1,M2,mp)%Y=SeIb(M1,M2,mp)COeffiCientSeIeCtiOnforbaseimage%%M1-COeffiCientsA%M2-COeffiCientsB%mp-SWitChforSeIeCtiOntype%mp==1:SeIeCtA% mp==2:SeIeCtB%Y -COmbinedCOeffiCients% (OIiVerRockinger16.08.99)SWitCh(mp)CaSe1,Y=M1;CaSe2,Y=M2;CaSe3,Y=(M1+M2)/2;otherwise,error('unknownOPtiOn');end;funCtionY=SeIC(M1,M2,ap)%Y=SeIC(M1,M2,ap)COeffiCinetSeIeCtiOnforhighpasscompOnentS%% M1 -COeffiCientsA% M2-COeffiCientsB%mp-SWitChforSeIeCtiOntype% mp==1:Choosemax(abs)% mp==2:salience/matchmeasureWiththreshold==.75(asproposedbyBUrtetal)% mp==3:choosemaxWithconSiStencyCheCk(asproposedbyLietal)% mp==4:SimPIechoosemax%%Y -combinedCOeffiCients% (OIiVerRockinger16.08.99)%CheCkinPUtS[z1s1]=SiZe(M1);[z2s2]=SiZe(M2);if(z1~=z2)|(s1~=s2)error('InputimagesarenotofSamesize');end;%SWitChtomethodSWitCh(ap(1))CaSe1,

%choosemax(abs)mm=(abs(M1))>(abs(M2));Y=(mm.*M1)+((~mm).*M2);CaSe2,%BUrtSmethodUm=ap(2);th=.75;%COmPUtesalience=COnv2(es2(M1.*M1,floor(um∕2)),Ones(um),'valid');=conv2(es2(M2.*M2,floor(um∕2)),ones(um),'valid');%COmPUtematchMA=conv2(es2(M1.*M2,floor(um∕2)),ones(um),'valid');MA=2*MA./(S1+S2+eps);%SeIeCtionml=MA>th;m2=S1>S2;w1=(0.5-0.5*(1-MA)/(1-th));=(~m1).*((m2.*M1)+((~m2).*M2));=Y+(ml.*((m2.*M1.*(1-w1))+((m2).*M2*w1)+((~m2).*M2.*(1-w1))+((~m2).*M1.*w1)));CaSe3,%LiSmethodUm=ap(2);%firstStePA1=Ordfilt2(abs(es2(M1,floor(um∕2))),um*um,ones(um));A2=ordfilt2(abs(es2(M