Fully three-dimensional OSEM-based image reconstruction for Compton imaging using optimized ordering schemes
Authors:
- Jae Sung LeeJLJae Sung Lee
- Chun Sik LeeCLChun Sik Lee
- Chan Hyeong KimCKChan Hyeong Kim
- Myung Chul LeeMLMyung Chul Lee
- Dong Soo LeeDLDong Soo Lee
- Soo-Jin LeeSLSoo-Jin Lee
- 11 August 2010
- journal article
- research article
- Published by IOP Publishing in Physics in Medicine & Biology
- Vol. 55 (17) , 5007-5027
Abstract
Although the ordered subset expectation maximization (OSEM) algorithm does not converge to a true maximum likelihood solution, it is known to provide a good solution if the projections that constitute each subset are reasonably balanced. The Compton scattered data can be allocated to subsets using scattering angles (SA) or detected positions (DP) or a combination of the two (AP (angles and positions)). To construct balanced subsets, the data were first arranged using three ordering schemes: the random ordering scheme (ROS), the multilevel ordering scheme (MLS) and the weighted-distance ordering scheme (WDS). The arranged data were then split into J subsets. To compare the three ordering schemes, we calculated the coefficients of variation (CVs) of angular and positional differences between the arranged data and the percentage errors between mathematical phantoms and reconstructed images. All ordering schemes showed an order-of-magnitude acceleration over the standard EM, and their computation times were similar. The SA-based MLS and the DP-based WDS led to the best-balanced subsets (they provided the largest angular and positional differences for SA- and DP-based arrangements, respectively). The WDS exhibited minimum CVs for both the SA- and DP-based arrangements (the deviation in mean angular and positional differences between the ordered subsets was smallest). The combination of AP and WDS yielded the best results with the lowest percentage errors by providing larger and more uniform angular and positional differences for the SA and DP arrangements, and thus, is probably optimal Compton camera reconstruction using OSEM.Keywords
This publication has 21 references indexed in Scilit:
- Two approaches to implementing projector–backprojector pairs for 3D reconstruction from Compton scattered dataNuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2007
- Effect of detector parameters on the image quality of Compton camera for 99mTcNuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2007
- Development of the multi-purpose gamma-ray detection system consisting of a double-sided silicon strip detector and a 25-segmented germanium detectorNuclear Physics A, 2005
- An analytical image reconstruction algorithm to compensate for scattering angle broadening in Compton camerasPhysics in Medicine & Biology, 2003
- Quantitative evaluation of information loss for Compton camerasIEEE Transactions on Nuclear Science, 1999
- Experimental results from the C-SPRINT prototype Compton cameraIEEE Transactions on Nuclear Science, 1999
- Application of spherical harmonics to image reconstruction for the Compton cameraPhysics in Medicine & Biology, 1998
- The weighted-distance scheme: a globally optimizing projection ordering method for ARTIEEE Transactions on Medical Imaging, 1997
- Accelerated image reconstruction using ordered subsets of projection dataIEEE Transactions on Medical Imaging, 1994
- Three-dimensional maximum-likelihood reconstruction for an electronically collimated single-photon-emission imaging systemJournal of the Optical Society of America A, 1990