CERN / CERN/TIGRE

Artifacts in spiral cone-beam CT reconstruction

Open
#622 10 comments 0 reactions 0 assignees View on GitHub
Dominant language
MATLAB
Stars
803
Forks
262
Avg merge
9d 13h
Merged PRs (30d)
8

Description

@AnderBiguri Hello
I tried to reconstruct the spiral cone beam CT image using an iterative reconstruction algorithm (SART, OS-SART), as I understand it only the geo.offOrigin parameter needs to be modified to synchronize the position of the bed at each angle, and the other parameters are almost unchanged, but I get the following image, which has scimitar-like artifacts in the image;
In fact, I found that it should not be a problem occurring in the projection process (Ax), it should be a problem occurring in the inverse projection process (Atb), in addition I did the same reconstruction directly using the FDK algorithm, and the result is the same as OS-SART, with the same artifacts.
the geo i used
![image](https://github.com/user-attachments/assets/82961d0f-2ac7-4677-bf83-b433aa265250)
![image](https://github.com/user-attachments/assets/12f34b5b-00aa-49eb-8d1f-a3191cba718f)
The spiral trajectory is as follows
![image](https://github.com/user-attachments/assets/19a5d31d-a1ac-48c0-8971-307c06ceb62f)

I also used these same geometric parameters for the sheppLogan3D simulation FDK reconstruction and the results show that the reconstructed sheppLogan image has no artifacts
the projections(two angles) is
![image](https://github.com/user-attachments/assets/73d77cd2-d96c-4de9-9175-44dbef2ca33c)
![image](https://github.com/user-attachments/assets/5858a646-00c4-419f-8a7a-6e8ac22a8422)

the FDK recon image(one slice) is (there is no artifacts)
![image](https://github.com/user-attachments/assets/677df615-b0fb-408f-80cf-469224ce4ddd)
![image](https://github.com/user-attachments/assets/f257e92a-96c4-47b0-8818-201783d3a942)

But the same geo used in real data, got the atifact results.
here is the OS-SART algorithm result(real data)
![image](https://github.com/user-attachments/assets/96d2ef4c-9014-454d-b18d-02944619392d)
And here is the FDK result
![image](https://github.com/user-attachments/assets/5c7cf136-e24c-4cd9-9b3c-8114a8d6b68e)

Do you have any idea what this could be due to?

## Specifications

- MATLAB/python version: MATLAB
- OS:
- CUDA version:

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.