Please help write a matlab code for fourier reconstruction according to the instructions above.

The original image f(x,y) to be reconstructed is a Shepp-Logan head phantom image using the code

n = 256;

phantom('Modified Shepp-Logan', n)

Thank you!

Basic Reconstruction Principles (1) Get Radon transform g(s,) of f(x,y) (2) 1D Fourier transform of the Radon transform data with respect to s (3) Fill 2D Fourier space with 1D Fourier transform result at different (4) Interpolate coordinate system (5) Perform inverse 2D Fourier Transform to obtain f(x,y) f(x,y)=IF2(n=1NF1,s(g(s,n))) A. C