CERN / TIGRE

TIGRE: Tomographic Iterative GPU-based Reconstruction Toolbox
BSD 3-Clause "New" or "Revised" License
581 stars 190 forks source link

Iterative reconstruction produces artifacts on both ends of the slice #587

Closed zezisme closed 1 month ago

zezisme commented 1 month ago

Expected Behavior

Hello ! When using cone-beam CT reconstruction, no matter how I define the reconstruction parameters, using the iterative reconstruction algorithm (SART, OS-SART, ...) will cause serious artifacts on both ends of the slice image, while the FDK algorithm is normal. What may be the reason for this? I have tested many iterative reconstruction algorithms and all of them have this problem. image

Actual Behavior

FDK image

OS-SART image image image image

other example(OS-SART) MPR view image

Code to reproduce the problem (If applicable)

image

Specifications

AnderBiguri commented 1 month ago

Yes, those are cone beam artifacts, natural in CBCT. They tend to be worse in iterative algorithms in the situation in where there was more object in the area you are not reconstructing, as iterative algorithms try to match the measured projections with the reconstruction, however as you don't have the whole image in the ROI, they put more values in the extrema to accomodate.

A standard way to mitigate these is to make the reconstruction volume (nVoxel and sVoxel) bigger in the Z direction. Then you can crop it out after recon if you want it the same size as FDK.

zezisme commented 1 month ago

Yes, those are cone beam artifacts, natural in CBCT. They tend to be worse in iterative algorithms in the situation in where there was more object in the area you are not reconstructing, as iterative algorithms try to match the measured projections with the reconstruction, however as you don't have the whole image in the ROI, they put more values in the extrema to accomodate.

A standard way to mitigate these is to make the reconstruction volume (nVoxel and sVoxel) bigger in the Z direction. Then you can crop it out after recon if you want it the same size as FDK.

Thank you for your clear explanation!I will try your suggestion