filterGaussian.m~

function result = filterGaussian(im, sigma, wind)
% filterGaussian is the finction of Gaus
 
%     if ~mod(wind, 2)
%         wind = wind + 1;
%     end
%     
%     filter = fspecial('gaussian', wind, sigma);
%     result = filter2(filter, im);
%     
    
    
    r2 = (wind-1)/2; 
    c2 = (wind-1)/2;
    
    [x,y] = meshgrid(-c2:c2, -r2:r2);
	radsqrd = x.^2 + y.^2;
	filter = exp(-radsqrd/(2*sigma^2));
	filter = filter/sum(filter(:));
    
    

    


    [nr, nc] = size(filter);
    result = conv2 (im(nr:-1:1, nc:-1:1), filter, 'same');

end