图像融合基于matlab小波分解和PCA重构图像融合【含Matlab源码 2595期】.zip

clc; clear all; n=2;%二层小波分解 m=3*n+1; X1=imread('0.TIF'); X1=double(X1)/256; %这里转化成double类型，否则使用小波变换输出的会有大量大于1的存在，会导致图像显示有问题 figure; imshow(X1),title('红外图'); X2=imread('1.TIF'); X2=double(X2)/256; figure; imshow(X2),title('可见光图'); [c1,s1]=wavedec2(X1,n,'sym4'); %将x1进行2维，使用‘sym4’进行变换 sizec=size(c1); [c2,s2]=wavedec2(X2,n,'sym4');%c1与c2的大小相等，即s1=s2 S1=s1(1,1)*s1(1,2); S2=s1(2,1)*s1(2,2); S3=s1(3,1)*s1(3,2); S4=s1(4,1)*s1(4,2); X1A=c1(1,1:S1); X1H1=c1(1,(S1+1):(S1+S2)); X1V1=c1(1,(S1+S2+1):(S1+2*S2)); X1D1=c1(1,(S1+2*S2+1):(S1+3*S2)); X1H2=c1(1,(S1+3*S2+1):(S1+3*S2+S3)); X1V2=c1(1,(S1+3*S2+S3+1):(S1+3*S2+2*S3)); X1D2=c1(1,(S1+3*S2+2*S3+1):(S1+3*S2+3*S3)); X2A=c2(1,1:S1); X2H1=c2(1,(S1+1):(S1+S2)); X2V1=c2(1,(S1+S2+1):(S1+2*S2)); X2D1=c2(1,(S1+2*S2+1):(S1+3*S2)); X2H2=c2(1,(S1+3*S2+1):(S1+3*S2+S3)); X2V2=c2(1,(S1+3*S2+S3+1):(S1+3*S2+2*S3)); X2D2=c2(1,(S1+3*S2+2*S3+1):(S1+3*S2+3*S3)); VA=cov(X1A,X2A); [COEFFA,latentA,explainedA]=pcacov(VA); for I=1:S1; c1(I)=explainedA(1,1)*c1(I)/100; c2(I)=explainedA(2,1)*c2(I)/100; end VH1=cov(X1H1,X2H1); [COEFFH1,latentH1,explainedH1]=pcacov(VH1); for I=(S1+1):(S1+S2); c1(I)=(explainedH1(1,1))^2*c1(I)/5000; c2(I)=(explainedH1(2,1))^2*c2(I)/5000; end VV1=cov(X1V1,X2V1); [COEFFV1,latentV1,explainedV1]=pcacov(VV1); for I=(S1+S2+1):(S1+2*S2); c1(I)=(explainedV1(1,1))^2*c1(I)/5000; c2(I)=(explainedV1(2,1))^2*c2(I)/5000; end VD1=cov(X1D1,X2D1); [COEFFD1,latentD1,explainedD1]=pcacov(VD1); for I=(S1+2*S2+1):(S1+3*S2); c1(I)=(explainedD1(1,1))^2*c1(I)/5000; c2(I)=(explainedD1(2,1))^2*c2(I)/5000; end VH2=cov(X1H2,X2H2); [COEFFH2,latentH2,explainedH2]=pcacov(VH2); for I=(S1+3*S2+1):(S1+3*S2+S3); c1(I)=(explainedH2(1,1))^3*c1(I)/250000; c2(I)=(explainedH2(2,1))^3*c2(I)/250000; end VV2=cov(X1V2,X2V2); [COEFFV2,latentV2,explainedV2]=pcacov(VV2); for I=(S1+3*S2+S3+1):(S1+3*S2+2*S3); c1(I)=(explainedV2(1,1))^3*c1(I)/250000; c2(I)=(explainedV2(2,1))^3*c2(I)/250000; end VD2=cov(X1D2,X2D2); [COEFFD2,latentD2,explainedD2]=pcacov(VD2); for I=(S1+3*S2+2*S3+1):(S1+3*S2+3*S3); c1(I)=(explainedD2(1,1))^3*c1(I)/250000; c2(I)=(explainedD2(2,1))^3*c2(I)/250000; end c=1.1*(c1+c2); s=s1; xx=waverec2(c,s,'sym4'); %进行重构 figure; imshow(xx),title('融合后的'); imwrite(xx,'权值指数.bmp','bmp'); %17.2861 0.2475 0.0873 0.0067 0.0608 0.7758 0.6062