基于ICA算法的图像融合matlab仿真+仿真录像

clc; clear; close all; warning off; addpath(genpath(pwd)); I1 = imread('pic4.jpg','jpeg'); I2 = imread('pic5.jpg','jpeg'); subplot(2,3,1), imshow(I1) subplot(2,3,4), imshow(I2) A=[0.8 0.2; 0.2 0.8]; X1=double(A(1,1)*I1+A(1,2)*I2); X2=double(A(2,1)*I1+A(2,2)*I2); subplot(2,3,2), imshow(uint8(X1)); subplot(2,3,5), imshow(uint8(X2)); [m,n]=size(X1); x1=reshape(X1,m*n,1); x2=reshape(X2,m*n,1); x1=x1-mean(x1); x2=x2-mean(x2); theta0=0.5*atan(-2*sum(x1.*x2)/sum(x1.^2-x2.^2)); Us=[cos(theta0) sin(theta0); -sin(theta0) cos(theta0)]; sig1=sum((x1*cos(theta0)+x2*sin(theta0)).^2); sig2=sum((x1*cos(theta0-pi/2)+x2*sin(theta0-pi/2)).^2); Sigma=[1/sqrt(sig1) 0; 0 1/sqrt(sig2)]; X1bar=Sigma(1,1)*(Us(1,1)*X1+Us(1,2)*X2); X2bar=Sigma(2,2)*(Us(2,1)*X1+Us(2,2)*X2); x1bar=reshape(X1bar,m*n,1); x2bar=reshape(X2bar,m*n,1); phi0=0.25*atan(-sum(2*(x1bar.^3).*x2bar ... -2*(x2bar.^3).*x1bar)/ ... sum(3*(x1bar.^2).*(x2bar.^2)-0.5*(x1bar.^4) ... -0.5*(x2bar.^4))); V=[cos(phi0) sin(phi0); -sin(phi0) cos(phi0)]; S1bar=V(1,1)*X1bar+V(1,2)*X2bar; S2bar=V(2,1)*X1bar+V(2,2)*X2bar; min1=min(min(min(S1bar))); S1bar=S1bar+min1; max1=max(max(max(S1bar))); S1bar=S1bar*(255/max1); min2=min(min(min(S2bar))); S2bar=S2bar+min2; max2=max(max(max(S2bar))); S2bar=S2bar*(255/max2); subplot(2,3,3), imshow(uint8(S1bar)); subplot(2,3,6), imshow(uint8(S2bar));