% filtr Cutoff

function [IMnew]=f_cutoff(IM,x)
kolik=round(x/2);
IMf=(fftshift(fft2(IM))); 
z=zeros(size(IM)); 
polovina=size(IM,1)/2;
%do nulove matice pridame jednicky (do stredu
z(polovina-kolik:polovina+kolik,polovina-kolik:polovina+kolik)=ones(kolik*2+1);
IMfd=IMf.*z;
IMnew=abs(ifft2(ifftshift(IMfd)));	%zpetny prevo FFT