Iterative 4D cardiac micro-CT image reconstruction using an adaptive spatio-temporal sparsity prior

Temporal-correlated image reconstruction, also known as 4D CT image reconstruction, is a big challenge in computed tomography. The reasons for incorporating the temporal domain into the reconstruction are motions of the scanned object, which would otherwise lead to motion artifacts. The standard method for 4D CT image reconstruction is extracting single motion phases and reconstructing them separately. These reconstructions can suffer from undersampling artifacts due to the low number of used projections in each phase. There are different iterative methods which try to incorporate some a priori knowledge to compensate for these artifacts. In this paper we want to follow this strategy. The cost function we use is a higher dimensional cost function which accounts for the sparseness of the measured signal in the spatial and temporal directions. This leads to the definition of a higher dimensional total variation. The method is validated using in vivo cardiac micro-CT mouse data. Additionally, we compare the results to phase-correlated reconstructions using the FDK algorithm and a total variation constrained reconstruction, where the total variation term is only defined in the spatial domain. The reconstructed datasets show strong improvements in terms of artifact reduction and low-contrast resolution compared to other methods. Thereby the temporal resolution of the reconstructed signal is not affected.

[1]  Frank Bergner,et al.  Low-dose cardio-respiratory phase-correlated cone-beam micro-CT of small animals. , 2011, Medical physics.

[2]  Steve B. Jiang,et al.  Low-dose 4DCT reconstruction via temporal nonlocal means. , 2010, Medical physics.

[3]  Jiayu Song,et al.  Sparseness prior based iterative image reconstruction for retrospectively gated cardiac micro-CT. , 2007, Medical physics.

[4]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[5]  C. Cametti,et al.  Microwave dielectric properties of normal and homozygous beta-thalassemic erythrocytes. , 1983, Physics in medicine and biology.

[6]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[7]  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.

[8]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[9]  Steve B. Jiang,et al.  Low-dose 4DCT reconstruction via temporal nonlocal means. , 2010, Medical physics.

[10]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[11]  W. Kalender,et al.  Empirical cupping correction: a first-order raw data precorrection for cone-beam computed tomography. , 2006, Medical physics.

[12]  M. Kachelriess,et al.  Improved total variation-based CT image reconstruction applied to clinical data , 2011, Physics in medicine and biology.

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

[14]  David L Donoho,et al.  Compressed sensing , 2006, IEEE Transactions on Information Theory.

[15]  M. Persson,et al.  Total variation norm for three-dimensional iterative reconstruction in limited view angle tomography , 2001, Physics in medicine and biology.

[16]  Paola Coan,et al.  X-ray phase-contrast imaging: from pre-clinical applications towards clinics , 2013, Physics in medicine and biology.

[17]  Gyula Faigel,et al.  X-Ray Holography , 1999 .

[18]  Willi A Kalender,et al.  Kymogram detection and kymogram-correlated image reconstruction from subsecond spiral computed tomography scans of the heart. , 2002, Medical physics.