Filtered-Backprojection

2023-11-20 Medisys lab intern study

Description

This is a Matlab code of Filtered Backprojection of Circle $$(x-a)^2 + (y-b)^2 = r^2$$ on coordinate plane.

Input

• Coordinate value of center of the circle x, y
• Radius
• Number of detectors & Sources (This shall an odd number)
• Distance between each detectors
• SDD, SCD
• Total views of CT
• Attenuation coefficient

Step 1: Forward Projection

testcase 1:

5
3
6
101
0.5
20
9
30
1

testcase 2:

5
3
12
101
0.5
30
14
72
1

output:

Step 2. Fourier Transform & Filter

output: Note that the Filtered sinogram graph requires resizing.

Step 3. Interpolation & Inverse Fourier Transform & FBP

Result: