How can I show a 2D FFT in frequency domain?

Hello, I looking for a method or a function to prsent my 2D fft in the frequency domain. For exemple, with the simple fft you can do like this :

NFFT =1000
Y = fft(y,NFFT);
f = Fs/2*linspace(0,1,NFFT/2+1);
plot(f,2*abs(Y(1:NFFT/2+1))) 
and we got the frequencies on the axis.

How can i do the same to have frequencies on the axis with the 2D FFT. And when we center the 2D FFT with an fftshift, we expect a zero on the center with the frequencies for the sides ... NO ?

I'm using imshow to present the spectrum in 2D :

FFT2D1= fft2(Sout1)
figure(1)
title('FFT2-Signal1')
imshow(20*log10(abs(FFT2D1)),[])
axis on
colorbar

Thanks for helping