function result = apply_bandpass_filter(image, low_cutoff, high_cutoff)
%   function result = apply_highpass_filter(image, low_cutoff, high_cutoff)
% $Id: apply_bandpass_filter.m,v 1.1 2002/06/27 11:47:17 tobin Exp $
[f1,f2] = freqspace(size(image),'meshgrid');
Hd = ones(size(image));
r = sqrt(f1.^2 + f2.^2);
Hd(r<low_cutoff | r>high_cutoff) = 0;
result = ifft2(fft2(image) .* fftshift(Hd));

     