Synchronized multiartifact reduction with tomographic reconstruction (SMART-RECON): A statistical model based iterative image reconstruction method to eliminate limited-view artifacts and to mitigate the temporal-average artifacts in time-resolved CT.

PURPOSE In x-ray computed tomography (CT), a violation of the Tuy data sufficiency condition leads to limited-view artifacts. In some applications, it is desirable to use data corresponding to a narrow temporal window to reconstruct images with reduced temporal-average artifacts. However, the need to reduce temporal-average artifacts in practice may result in a violation of the Tuy condition and thus undesirable limited-view artifacts. In this paper, the authors present a new iterative reconstruction method, synchronized multiartifact reduction with tomographic reconstruction (SMART-RECON), to eliminate limited-view artifacts using data acquired within an ultranarrow temporal window that severely violates the Tuy condition. METHODS In time-resolved contrast enhanced CT acquisitions, image contrast dynamically changes during data acquisition. Each image reconstructed from data acquired in a given temporal window represents one time frame and can be denoted as an image vector. Conventionally, each individual time frame is reconstructed independently. In this paper, all image frames are grouped into a spatial-temporal image matrix and are reconstructed together. Rather than the spatial and/or temporal smoothing regularizers commonly used in iterative image reconstruction, the nuclear norm of the spatial-temporal image matrix is used in SMART-RECON to regularize the reconstruction of all image time frames. This regularizer exploits the low-dimensional structure of the spatial-temporal image matrix to mitigate limited-view artifacts when an ultranarrow temporal window is desired in some applications to reduce temporal-average artifacts. Both numerical simulations in two dimensional image slices with known ground truth and in vivo human subject data acquired in a contrast enhanced cone beam CT exam have been used to validate the proposed SMART-RECON algorithm and to demonstrate the initial performance of the algorithm. Reconstruction errors and temporal fidelity of the reconstructed images were quantified using the relative root mean square error (rRMSE) and the universal quality index (UQI) in numerical simulations. The performance of the SMART-RECON algorithm was compared with that of the prior image constrained compressed sensing (PICCS) reconstruction quantitatively in simulations and qualitatively in human subject exam. RESULTS In numerical simulations, the 240(∘) short scan angular span was divided into four consecutive 60(∘) angular subsectors. SMART-RECON enables four high temporal fidelity images without limited-view artifacts. The average rRMSE is 16% and UQIs are 0.96 and 0.95 for the two local regions of interest, respectively. In contrast, the corresponding average rRMSE and UQIs are 25%, 0.78, and 0.81, respectively, for the PICCS reconstruction. Note that only one filtered backprojection image can be reconstructed from the same data set with an average rRMSE and UQIs are 45%, 0.71, and 0.79, respectively, to benchmark reconstruction accuracies. For in vivo contrast enhanced cone beam CT data acquired from a short scan angular span of 200(∘), three 66(∘) angular subsectors were used in SMART-RECON. The results demonstrated clear contrast difference in three SMART-RECON reconstructed image volumes without limited-view artifacts. In contrast, for the same angular sectors, PICCS cannot reconstruct images without limited-view artifacts and with clear contrast difference in three reconstructed image volumes. CONCLUSIONS In time-resolved CT, the proposed SMART-RECON method provides a new method to eliminate limited-view artifacts using data acquired in an ultranarrow temporal window, which corresponds to approximately 60(∘) angular subsectors.

[1]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

[2]  Jie Tang,et al.  Prior image constrained compressed sensing: implementation and performance evaluation. , 2011, Medical physics.

[3]  Kai Niu,et al.  Statistical image reconstruction via denoised ordered-subset statistically penalized algebraic reconstruction technique (DOS-SPART) , 2014, Medical Imaging.

[4]  Jie Tang,et al.  Temporal resolution improvement using PICCS in MDCT cardiac imaging. , 2009 .

[5]  Rolf Clackdoyle,et al.  Tomographic Reconstruction in the 21st Century , 2010, IEEE Signal Processing Magazine.

[6]  Hengyong Yu,et al.  A two-step Hilbert transform method for 2 D image reconstruction ” , 2015 .

[7]  Hiroyuki Kudo,et al.  Image reconstruction from fan-beam projections on less than a short scan , 2002, Physics in medicine and biology.

[8]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[9]  Jie Tang,et al.  Tomosynthesis via total variation minimization reconstruction and prior image constrained compressed sensing (PICCS) on a C-arm system , 2008, SPIE Medical Imaging.

[10]  Jie Tang,et al.  Temporal resolution improvement in cardiac CT using PICCS (TRI-PICCS): performance studies. , 2010, Medical physics.

[11]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[12]  H. Tuy.  AN INVERSION FORMULA FOR CONE-BEAM RECONSTRUCTION* , 1983 .

[13]  Katsuyuki Taguchi,et al.  Achieving routine submillisievert CT scanning: report from the summit on management of radiation dose in CT. , 2012, Radiology.

[14]  Guang-Hong Chen,et al.  A new framework of image reconstruction from fan beam projections. , 2003, Medical physics.

[15]  Xiaochuan Pan,et al.  Image reconstruction in peripheral and central regions-of-interest and data redundancy. , 2005, Medical physics.

[16]  Guang-Hong Chen,et al.  Temporally targeted imaging method applied to ECG-gated computed tomography: preliminary phantom and in vivo experience. , 2008, Academic radiology.

[17]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[18]  F. Noo,et al.  A two-step Hilbert transform method for 2D image reconstruction. , 2004, Physics in medicine and biology.

[19]  Xiaochuan Pan,et al.  Image reconstruction in regions-of-interest from truncated projections in a reduced fan-beam scan , 2005, Physics in medicine and biology.

[20]  A. Bovik,et al.  A universal image quality index , 2002, IEEE Signal Processing Letters.

[21]  Guang-Hong Chen,et al.  Temporal resolution improvement using PICCS in MDCT cardiac imaging. , 2009, Medical physics.

[22]  Xiaochuan Pan,et al.  Exact image reconstruction on PI-lines from minimum data in helical cone-beam CT. , 2004, Physics in medicine and biology.

[23]  Jeffrey A. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, IEEE Trans. Image Process..