MATLAB Implementation of Filtered Back Projection for CT Image Reconstruction

This is a MATLAB-implemented filtered back projection program for CT image reconstruction, designed for medical imaging and image processing applications. The filtered back projection algorithm employs Fourier transform techniques and convolution operations to process raw projection data. Key functions include ramlak filtering for frequency domain correction and radon/inverse radon transforms for projection/backprojection operations. The program allows parameter tuning of filter types (e.g., Shepp-Logan, Cosine) and projection angles to optimize reconstruction quality for different scenarios. It incorporates image visualization tools using imshow() and imwrite() functions for result analysis and export. The reconstruction pipeline involves: 1) Log-transform of projection data, 2) Frequency domain filtering with zero-padding, 3) Back projection using interpolation methods. This implementation supports DICOM format input/output and includes error handling for invalid projection data. We hope this program proves valuable for your research and practical applications!