function h = d2gauss0(n1,std1,n2,std2)
% This function returns a 2D Gaussian filter with size n1*n2; the 
% correlation is zero; and sigma1 and sigma2
% are the standard deviation of the Gaussian functions.

for i = 1 : n2 
    for j = 1 : n1
        h(i,j) = gauss(j-(n1+1)/2,std1)*gauss(i-(n2+1)/2,std2);
    end
end
h = h / sum(sum(h));

