高斯图像融合

医学图像融合，普通图像融合，可运行，首先利用Contourlet分频，高频绝对值最大融合，低频采用高斯融合

function E= Expectation（XSVarianceWb）
X=X-b*S;
[NL]=size（X）;
E=zeros（N4）;
% [nd] = size（X）;
% a = （2*pi）^（0.5*d）;
% S = zeros（1k）;
% iV = zeros（ddk）;
% for j=1:k
%     if V（::j）==zeros（dd） V（::j）=ones（dd）*eps; end
%     S（j） = sqrt（det（V（::j）））;
%     iV（::j） = inv（V（::j））;    
% end
% E = zeros（nk）;
% for i=1:n    
%     for j=1:k
%         dXM = X（i:）‘-M（:j）;
%         pl = exp（-0.5*dXM‘*iV（::j）*dXM）/（a*S（j））;
%         E（ij） = W（j）*pl;
%     end
%     E（i:） = E（i:）/sum（E（i:））;
% end
for i=1:N
    FX1=0;FX2=0;FX3=0;
     FX4=0;
     %%%%%%%%%%%%%%%%%%%%%%%%一部分%%%%%%%%%%%%%%%
        e1=exp（-（X（i）^2/（2*Variance（1））））;
        FX1=FX1+e1;     
     %%%%%%%%%%%%%%%%%%%%%%%%二部分%%%%%%%%%%%%%%%
      e2=exp（-（X（i）^2/（2*Variance（2））））;
        FX2=FX2+e2;      
     %%%%%%%%%%%%%%%%%%%%%%%%叁部分%%%%%%%%%%%%%%%
  e3=exp（-（X（i）^2/（2*Variance（3））））;
        FX3=FX3+e3;      
     %%%%%%%%%%%%%%%%%%%%%%%%四部分%%%%%%%%%%%%%%%
    e4=exp（-（X（i）^2/（2*Variance（4））））;
        FX4=FX4+e4;         
%     FX1=（FX1+1）;
%     FX2=（FX2+1）;
%     FX3=（FX3+1）;
%     FX4=（FX4+1）;

%     FX6=（FX6+2/255）;
    sum=W（1）*FX1+W（2）*FX2+W（3）*FX3+W（4）*FX4;
    E（i1）=W（1）*FX1/sum;
    E（i2）=W（2）*FX2/sum;
    E（i3）=W（3）*FX3/sum;
    E（i4）=W（4）*FX4/sum;

end

    
    
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-01-24 17:51  Gaussian\
     文件      104907  2015-12-25 13:42  Gaussian\10.jpg
     文件       64285  2015-12-25 13:42  Gaussian\3.jpg
     文件        1359  2015-12-25 13:42  Gaussian\Expectation.m
     文件         882  2015-12-25 13:42  Gaussian\f.m
     文件       66614  2015-12-29 17:18  Gaussian\fusion1.bmp
     文件         317  2015-12-25 13:42  Gaussian\gaoping.m
     文件         707  2015-12-25 13:42  Gaussian\Gmm.m
     文件       17263  2015-12-25 13:42  Gaussian\high.jpg
     文件         675  2015-12-25 13:42  Gaussian\Init_EM.m
     文件         926  2015-12-25 13:42  Gaussian\Init_EM_1.m
     文件       51340  2015-12-25 13:42  Gaussian\IRIamge.jpg
     文件       12150  2015-12-25 13:42  Gaussian\low.jpg
     文件        1090  2015-12-25 13:42  Gaussian\main.m
     文件        3449  2015-12-25 13:42  Gaussian\Maximization.m
     文件        3928  2015-12-25 13:42  Gaussian\medA-t.jpg
     文件        3928  2015-12-25 13:42  Gaussian\medA.jpg
     文件        7904  2015-12-25 13:42  Gaussian\medB-t.jpg
     文件        7904  2015-12-25 13:42  Gaussian\medB.jpg
     文件         661  2015-12-25 13:42  Gaussian\myfcmspa24U.m
     文件        1075  2015-12-25 13:42  Gaussian\myfcmspa8U.m
     文件         734  2015-12-25 13:42  Gaussian\myfcmspaU.m
     文件        8352  2015-12-25 13:42  Gaussian\naviA-t.jpg
     文件       10429  2015-12-25 13:42  Gaussian\naviB-t.jpg
     文件         270  2015-12-25 13:42  Gaussian\region.m
     文件       17564  2015-12-25 13:42  Gaussian\remoteA-t.jpg
     文件       17429  2015-12-25 13:42  Gaussian\remoteB-t.jpg
     文件      241078  2015-12-25 13:42  Gaussian\spot.bmp
     文件      720054  2015-12-25 13:42  Gaussian\tm.bmp
     文件         516  2015-12-25 13:42  Gaussian\Untitled.asv
     文件         495  2015-12-25 13:42  Gaussian\Untitled.m
............此处省略19个文件信息