A new convergent MAP reconstruction algorithm for emission tomography using ordered subsets and separable surrogates
Proceedings IEEE International Symposium on Biomedical Imaging - Trang 409-412
Tóm tắt
We investigate a new, fast and provably convergent MAP reconstruction algorithm for emission tomography. The new algorithm, termed C-OSEM has its origin in the alternating algorithm derivation of the well known EM algorithm for emission tomography. In this re-derivation, the complete data explicitly enters the objective function as an unknown variable. While the entire complete data gets updated in each iteration of EM, in C-OSEM the complete data is updated only along ordered subsets. C-OSEM has a straightforward extension to the MAP case especially when using convex, smoothing priors. Unlike RAMLA and BSREM, C-OSEM does not require relaxation parameters to be set at each iteration. We derive the MAP C-OSEM algorithm using the separable surrogate method and anecdotally compare performance with MAP EM and BSREM.
Từ khóa
#Reconstruction algorithms #Tomography #Image reconstruction #Smoothing methods #Bayesian methods #Cost function #Biomedical imaging #Radiology #Physics #Data modelsTài liệu tham khảo
hsiao, 2002, A, provably convergent OS-EM like reconstruction algorithm for emission to-mography, SPIE Proc Medical Imaging, 10.1117/12.467144
10.1109/42.921477
shepp, 1982, Maximization likelihood reconstruction for emission tomography, IEEE Trans Med lmag, 1, 113, 10.1109/TMI.1982.4307558
csiszar, 1984, Information geometry and alternating minimization procedures, Statistics and Decisions, 205
10.1109/42.232263
10.1109/42.538946
10.1109/42.363108