TY - JOUR
T1 - Line integral alternating minimization algorithm for dual-energy X-ray CT image reconstruction
AU - Chen, Yaqi
AU - O'Sullivan, Joseph A.
AU - Politte, David G.
AU - Evans, Joshua D.
AU - Han, Dong
AU - Whiting, Bruce R.
AU - Williamson, Jeffrey F.
N1 - Publisher Copyright:
© 2015 IEEE.
PY - 2016/2
Y1 - 2016/2
N2 - We propose a new algorithm, called line integral alternating minimization (LIAM), for dual-energy X-ray CT image reconstruction. Instead of obtaining component images by minimizing the discrepancy between the data and the mean estimates, LIAM allows for a tunable discrepancy between the basis material projections and the basis sinograms. A parameter is introduced that controls the size of this discrepancy, and with this parameter the new algorithm can continuously go from a two-step approach to the joint estimation approach. LIAM alternates between iteratively updating the line integrals of the component images and reconstruction of the component images using an image iterative deblurring algorithm. An edge-preserving penalty function can be incorporated in the iterative deblurring step to decrease the roughness in component images. Images from both simulated and experimentally acquired sinograms from a clinical scanner were reconstructed by LIAM while varying the regularization parameters to identify good choices. The results from the dual-energy alternating minimization algorithm applied to the same data were used for comparison. Using a small fraction of the computation time of dual-energy alternating minimization, LIAM achieves better accuracy of the component images in the presence of Poisson noise for simulated data reconstruction and achieves the same level of accuracy for real data reconstruction.
AB - We propose a new algorithm, called line integral alternating minimization (LIAM), for dual-energy X-ray CT image reconstruction. Instead of obtaining component images by minimizing the discrepancy between the data and the mean estimates, LIAM allows for a tunable discrepancy between the basis material projections and the basis sinograms. A parameter is introduced that controls the size of this discrepancy, and with this parameter the new algorithm can continuously go from a two-step approach to the joint estimation approach. LIAM alternates between iteratively updating the line integrals of the component images and reconstruction of the component images using an image iterative deblurring algorithm. An edge-preserving penalty function can be incorporated in the iterative deblurring step to decrease the roughness in component images. Images from both simulated and experimentally acquired sinograms from a clinical scanner were reconstructed by LIAM while varying the regularization parameters to identify good choices. The results from the dual-energy alternating minimization algorithm applied to the same data were used for comparison. Using a small fraction of the computation time of dual-energy alternating minimization, LIAM achieves better accuracy of the component images in the presence of Poisson noise for simulated data reconstruction and achieves the same level of accuracy for real data reconstruction.
KW - Alternating minimization algorithm
KW - Dual-energy
KW - Iterative deblurring algorithm
KW - Line integral
KW - X-ray CT
UR - http://www.scopus.com/inward/record.url?scp=84959333875&partnerID=8YFLogxK
U2 - 10.1109/TMI.2015.2490658
DO - 10.1109/TMI.2015.2490658
M3 - Article
C2 - 26469126
AN - SCOPUS:84959333875
SN - 0278-0062
VL - 35
SP - 685
EP - 698
JO - IEEE Transactions on Medical Imaging
JF - IEEE Transactions on Medical Imaging
IS - 2
M1 - 2490658
ER -