saa_admm_paper

This repository contains the code to reproduce the experimental results from the paper "Simultaneous activity and attenuation estimation in TOF-PET with TV-constrained nonconvex optimization".

Prerequisuites

• python
• numpy
• numba
• pylab
• scipy

Usage

• noiseless-setting.ipynb contains the code to reproduce all the numerical results in the noisess setting (Figure 1-4). best optimization parameters (Figure 2).
• noisy-setting.ipynb contains the code to reproduce the results in the noisy setting (Figure 5,6,8, and 10).