/near-circulant-splitting

Code for the paper "Splitting with Near-Circulant Linear Systems: Applications to Total Variation CT and PET"

Primary LanguageMATLABApache License 2.0Apache-2.0

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

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.