基于EM算法的图像融合算法，对图像不用进行多尺度分解

function gaussian12
f=imread('A1.jpg');
g=imread('B1.jpg');

f=rgb2gray(f);
f=im2double(f);

g=rgb2gray(g);
g=im2double(g);
%imshow(f,[]);
h=ones(256,256);
h(1,:)=0.5*f(1,:)+0.5*g(1,:);
h(200,:)=0.5*f(200,:)+0.5*g(200,:);
h(:,1)=0.5*f(:,1)+0.5*g(:,1);
h(:,256)=0.5*f(:,256)+0.5*g(:,256);
%mean=0.5*f+0.5*g;
%w=mean;
for j=2:255
     for i=2:199
        x=[f(i-1,j-1) f(i-1,j) f(i-1,j+1) f(i,j-1) f(i,j) f(i,j+1) f(i+1,j-1) f(i+1,j) f(i+1,j+1)];
        y=[g(i-1,j-1) g(i-1,j) g(i-1,j+1) g(i,j-1) g(i,j) g(i,j+1) g(i+1,j-1) g(i+1,j) g(i+1,j+1)];
        mean=0.5*f(i,j)+0.5*g(i,j);
        [Z]=imagefision(mean,x,y,2,15);
        %fprintf( '%d iterations\n',15);
       h(i,j)=Z;
       %if abs(h(i,j)-mean)>4.531
        
        % h(i,j)=mean;end 
   end
end
imshow(h,[]);
%imwrite(h,'D:\magic\out7.bmp');
%w8=[0 1 0;1 -4 1;0 1 0];
%f1=imfilter(f,w8,'replicate');
%figure,imshow(f1);