MATLAB Source Code for CT Fan-Beam Equiangular Filtered Backprojection Reconstruction Algorithm

This article explores the MATLAB source code implementation for CT fan-beam equiangular filtered backprojection reconstruction algorithm. The codebase is systematically organized into three modular files: source image generation, projection simulation, and filtered backprojection reconstruction. We begin by examining the fundamentals of CT fan-beam imaging technology, including its operational principles and the mathematical foundation of equiangular filtered backprojection reconstruction algorithm. The implementation demonstrates how to create synthetic phantom images using MATLAB's matrix operations, simulate fan-beam projections through radon transform modifications, and apply frequency-domain filtering with ramp filters before backprojection. The reconstruction process utilizes interpolation techniques for accurate pixel-driven backprojection in equiangular geometry. Finally, we discuss optimization strategies including computational efficiency improvements through vectorization, memory management for large datasets, and potential parallel processing implementations using MATLAB's Parallel Computing Toolbox.