Power Spectrum of a fingerprint using annular band

I am trying to plot/ show a power spectrum of the fingerprint images in a similar way as shown in the paper which has been attached below. I was successful in determining the Power spectrum as shown in the code below:

img= rgb2gray(imread ('finger.jpg'));
[M,N,~]=size(img);
F = fft2(img);
Fd = (1/N*M).* F;
figure;imshow(mat2gray(log(1+abs(fftshift(F)))));
figure;imshow(mat2gray(log(1+abs(fftshift(Fd)))));
psd = 10*log10(abs(fftshift(fft2(img))).^2 );
figure; clf
mesh(psd);

This provides me a 3D plot of the power spectrum:

But I would like to have a plot as shown below:

I was trying to learn more about a paper regarding fingerprints and was trying to implement it. The paper along with the code file and the fingerprint is given in the link attached below.

https://drive.google.com/folderview?id=0B3IeQ-na-m9rQlRobGRLbmdUWnc&usp=sharing

Thanks in advance.

Answer 1

[Image Analyst]

I'm not going to read the paper, but an image like that can be created by

imshow(grayImage, []);
axis on;
colormap(jet(256));

You might also want to call

colorbar;