nengliang.rar_图像能量_图像融合matlab_局部 图像_局部能量融合_融合

% 基于局部能量最大（最小）选取规则 function [fusion,X1,X2]=jiyujubunengliangzuidaozuixiaoxuanquguize; orginal1=imread('fused009_1.tif'); X1=orginal1; figure; X1=double(X1); [column,row]=size(X1); X1=X1+10*randn(column,row); subplot(2,2,1); colormap(gray); imagesc(X1); % 装入图像 orginal2=imread('fused009_2.tif'); X2=orginal2; X2=double(X2); subplot(2,2,2); colormap(gray); imagesc(X2); [column,row]=size(X1); % % 选取阈值 T=0.8; for i=1:column for j=1:row E1(i,j)=X1(i,j)^2; E2(i,j)=X2(i,j)^2; M(i,j)=2*X1(i,j)*X2(i,j)/(E1(i,j)+E2(i,j)); if M(i,j)<=T if E1(i,j)>=E2(i,j) fusion(i,j)=X1(i,j); else fusion(i,j)=X2(i,j); end else Wmin(i,j)=1/2-1/2*((1-M(i,j))/(1-T)); Wmax(i,j)=1-Wmin(i,j); if E1(i,j)>=E2(i,j) fusion(i,j)=Wmax(i,j)*X1(i,j)+Wmin(i,j)*X2(i,j); else fusion(i,j)=Wmin(i,j)*X1(i,j)+Wmax(i,j)*X2(i,j); end end end end subplot(2,2,3); colormap(gray); imagesc(fusion);