function [ new T error ] = compress1( image )
% Compress the original image by setting low Fourier coefficients to zero. 

F=fft2(image);
for T=5:2:200
   FC=F.*(abs(F)>T);
   new=ifft2(FC);
   error=abs(new-image);
   error=std(error(:));
   if error>0.005      
       break;
   end;
end;
