Background: Simultaneous reconstruction of activity and attenuation with MLAA offers an interesting option to solve the problem of attenuation correction in hybrid PET/MR scanners. The algorithm maximizes the joint likelihood by alternating MLEM and MLTR steps. Any scheme where a number of MLTR iterations follow any number of MLEM ones converges to the true solution, but convergence speed and computation time might differ significantly. We investigate which alternating strategy is expected to give optimal results as a function of TOF coincidence timing resolution (CTR), given the different computation time and convergence speed between TOF-MLEM and non-TOF MLTR. Following, acceleration strategies exploiting ordered subset (OS) techniques are proposed and tested. Methods: A slice of a 2D digital body phantom was simulated and forward projected with CTR of 200, 400, 600 and 800 ps. Reconstructions were performed using the corresponding CTR using 1, 2, 4 and 8 MLTR iterations for each MLEM iteration, using native-geometry, ray-driven forward and back-projectors. Convergence curves were plotted as a function of computation time. Different strategies of OS acceleration were then tested and compared using simulations both noiseless and featuring Poisson noise, using 32 subsets. Results: Convergence was achieved with the least computation time when 4 or 8 MLTR iterations were performed for each MLEM one. All proposed OS strategies had results similar to reference non-accelerated versions and performed equivalently also in presence of noise, despite the high subsets number. Conclusions: Performing 4 or 8 MLTR iterations for each MLEM one appears to be the optimal solution for current implementations of the MLAA algorithm. OS techniques are effective in reducing computation time and perform comparably to non-accelerated versions.

Presotto, L., Bettinardi, V., Gianolli, L., Perani, D. (2016). Alternating strategies and ordered subset acceleration schemes for maximum likelihood activity and attenuation reconstruction in time-of-flight PET. In 2015 IEEE Nuclear Science Symposium and Medical Imaging Conference, NSS/MIC 2015. Institute of Electrical and Electronics Engineers Inc. [10.1109/NSSMIC.2015.7582116].

Alternating strategies and ordered subset acceleration schemes for maximum likelihood activity and attenuation reconstruction in time-of-flight PET

Presotto L.
;
2016

Abstract

Background: Simultaneous reconstruction of activity and attenuation with MLAA offers an interesting option to solve the problem of attenuation correction in hybrid PET/MR scanners. The algorithm maximizes the joint likelihood by alternating MLEM and MLTR steps. Any scheme where a number of MLTR iterations follow any number of MLEM ones converges to the true solution, but convergence speed and computation time might differ significantly. We investigate which alternating strategy is expected to give optimal results as a function of TOF coincidence timing resolution (CTR), given the different computation time and convergence speed between TOF-MLEM and non-TOF MLTR. Following, acceleration strategies exploiting ordered subset (OS) techniques are proposed and tested. Methods: A slice of a 2D digital body phantom was simulated and forward projected with CTR of 200, 400, 600 and 800 ps. Reconstructions were performed using the corresponding CTR using 1, 2, 4 and 8 MLTR iterations for each MLEM iteration, using native-geometry, ray-driven forward and back-projectors. Convergence curves were plotted as a function of computation time. Different strategies of OS acceleration were then tested and compared using simulations both noiseless and featuring Poisson noise, using 32 subsets. Results: Convergence was achieved with the least computation time when 4 or 8 MLTR iterations were performed for each MLEM one. All proposed OS strategies had results similar to reference non-accelerated versions and performed equivalently also in presence of noise, despite the high subsets number. Conclusions: Performing 4 or 8 MLTR iterations for each MLEM one appears to be the optimal solution for current implementations of the MLAA algorithm. OS techniques are effective in reducing computation time and perform comparably to non-accelerated versions.
No
poster + paper
emission tomography; regularization; image reconstruction;
English
2015 IEEE Nuclear Science Symposium and Medical Imaging Conference, NSS/MIC 2015 - 31 October 2015 through 7 November 2015
9781467398626
Presotto, L., Bettinardi, V., Gianolli, L., Perani, D. (2016). Alternating strategies and ordered subset acceleration schemes for maximum likelihood activity and attenuation reconstruction in time-of-flight PET. In 2015 IEEE Nuclear Science Symposium and Medical Imaging Conference, NSS/MIC 2015. Institute of Electrical and Electronics Engineers Inc. [10.1109/NSSMIC.2015.7582116].
Presotto, L; Bettinardi, V; Gianolli, L; Perani, D
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10281/381170
Citazioni
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 0
Social impact