function h = MasqueGaussien(sigma1,sigma2,n1,n2); 

% function h = MasqueGaussien(sigma1,sigma2,n1,n2); 
%
% Realise un masque gaussien d'ecart-type sigma normalise 
% (moyenne 1) 

t1 = [0:n1/2 -n1/2+1:-1];
t2 = [0:n2/2 -n2/2+1:-1]; 
[X1,X2] = meshgrid (t1,t2); 
h=exp(-X1.^2/(2*sigma1^2) - X2.^2/(2*sigma2^2)); 
h=h/sum(h(:)); 
end



