Reduce the steps within this code

Hi there!

I've been wondering of a way to cut down the steps of sorting the transform and thresholding which constructs an image based on specific # of coefficicents. f is the image and M=4096:

R = fft2(f);
s = sort(abs(R(:)));
s = s(end:-1:1);
T = s(M+1);
R = R.*(abs(R)>T);
f1 = real(ifft2(R));

Just a thought. What do you suggest?

Answer 1

[Amit]

R = fft2(f); 
s = sort(abs(R(:)),1,'descend'); 
%s = s(end:-1:1); 
T = s(M+1); 
%R = R.*(abs(R)>T); 
f1 = real(ifft2(R.*(abs(R)>T)));

[Ayesha]

Thanks a ton!!