%Gausian filtr

function [IMnew]=f_gauss(IM,s)
IMf=(fftshift(fft2(IM))); 
z=fspecial('gausian',size(IM),s);
z=z/max(z(:)); 
IMfd=IMf.*z;	
IMnew=abs(ifft2(ifftshift(IMfd)));