Design of iterative ROI transmission tomography reconstruction procedures and image quality analysis.

PURPOSE An iterative edge-preserving CT reconstruction algorithm for high-resolution imaging of small regions of the field of view is investigated. It belongs to a family of region-of-interest reconstruction techniques in which a low-cost pilot reconstruction of the whole field of view is first performed and then used to deduce the contribution of the region of interest to the projection data. These projections are used for a high-resolution reconstruction of the region of interest (ROI) using a regularized iterative algorithm, resulting in significant computational savings. This paper examines how the technique by which the pilot reconstruction of the full field of view is obtained affects the total runtime and the image quality in the region of interest. METHODS Previous contributions to the literature have each focused on a single approach for the pilot reconstruction. In this paper, two such approaches are compared: the filtered backprojection and a low-resolution regularized iterative reconstruction method. ROI reconstructions are compared in terms of image quality and computational cost over simulated and physical phantom (Catphan600) studies, in order to assess the compromises that most impact the quality of the ROI reconstruction. RESULTS With the simulated phantom, new artifacts that appear in the ROI images are caused by significant errors in the pilot reconstruction. These errors include excessive coarseness of the pilot image grid and beam-hardening artifacts. With the Catphan600 phantom, differences in the imaging model of the scanner and that of the iterative reconstruction algorithm cause dark border artifacts in the ROI images. CONCLUSIONS Inexpensive pilot reconstruction techniques (analytical algorithms, very-coarse-grid penalized likelihood) are practical choices in many common cases. However, they may yield background images altered by edge degradation or beam hardening, inducing projection inconsistency in the data used for ROI reconstruction. The ROI images thus have significant streak and speckle artifacts, which adversely affect the resolution-to-noise compromise. In these cases, edge-preserving penalized-likelihood methods on not-too-coarse image grids prove to be more robust and provide the best ROI image quality.

[1]  Benoit Hamelin,et al.  Comparison of optimization techniques for regularized statistical reconstruction in X-ray tomography , 2010, 2010 2nd International Conference on Image Processing Theory, Tools and Applications.

[2]  Wojciech Zbijewski,et al.  Characterization and suppression of edge and aliasing artefacts in iterative x-ray CT reconstruction. , 2004, Physics in medicine and biology.

[3]  Patrick J. La Rivière,et al.  Penalized-likelihood sinogram restoration for computed tomography , 2006, IEEE Transactions on Medical Imaging.

[4]  Jorge Nocedal,et al.  Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization , 1997, TOMS.

[5]  Yongyi Yang,et al.  Tomographic image reconstruction based on a content-adaptive mesh model , 2004, IEEE Transactions on Medical Imaging.

[6]  Benoit Hamelin,et al.  Iterative CT reconstruction of real data with metal artifact reduction , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[7]  Armand Wirgin,et al.  The inverse crime , 2004, math-ph/0401050.

[8]  Michel Barlaud,et al.  Deterministic edge-preserving regularization in computed imaging , 1997, IEEE Trans. Image Process..

[9]  P.J. La Riviere.  Monotonic iterative reconstruction algorithms for targeted reconstruction in emission and transmission computed tomography , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[10]  J. Fessler.  Statistical Image Reconstruction Methods for Transmission Tomography , 2000 .

[11]  M. Defrise,et al.  Solving the interior problem of computed tomography using a priori knowledge , 2008, Inverse problems.

[12]  R L Morin,et al.  A practical method to measure the MTF of CT scanners. , 1982, Medical physics.

[13]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[14]  Zhou Yu,et al.  Non-homogeneous ICD optimization for targeted reconstruction of volumetric CT , 2008, Electronic Imaging.

[15]  Alfred O. Hero,et al.  Convergent incremental optimization transfer algorithms: application to tomography , 2006, IEEE Transactions on Medical Imaging.

[16]  A. Karimi,et al.  Master‟s thesis , 2011 .

[17]  Y. Goussard,et al.  Penalized-likelihood region-of-interest CT reconstruction by local object supersampling , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[18]  H. Kudo,et al.  Practical statistical models for region-of-interest tomographic reconstruction and long object problem , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[19]  Michael Grass,et al.  Iterative reconstruction of a region of interest for transmission tomography , 2006, SPIE Medical Imaging.

[20]  Benoit Hamelin,et al.  On the stopping criterion for numerical methods for linear systems with additive Gaussian noise , 2009 .

[21]  Hiroyuki Kudo,et al.  Truncated Hilbert transform and image reconstruction from limited tomographic data , 2006 .

[22]  M. Defrise,et al.  Tiny a priori knowledge solves the interior problem in computed tomography , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[23]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[24]  S.D. Wollenweber,et al.  An Efficient Algorithm for Targeted Reconstruction of Tomographic Data , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.