Splitting with Near-Circulant Linear Systems: Applications to Total Variation CT and PET

This code implements the experiments for reference (click here for Arxiv preprint).

[1] E. K. Ryu, S. Ko, J.-H. Won (2020). "Splitting with Near-Circulant Linear Systems: Applications to Total Variation CT and PET," SIAM Journal of Scientific Computing, 42(1), B185-B206.

Authors

• Ernest K. Ryu
• Seyoon Ko
• Joong-Ho Won

Requires MIRT by Fessler Lab

http://web.eecs.umich.edu/~fessler/irt/fessler.tgz

Code

• par_beam.m
• fan_beam.m
• cone_beam_temp.m
• PET.m

Notes

• cone_beam.m does not run on Windows and does not utilize GPUs as the ct_geom of MIRT does not support Windows or GPUs.
• We cannot relase the real patient data used for par_beam.m. You can experiment instead with the Shepp-Logan phantom.