Shearlet-based regularization in sparse dynamic tomography

Classical tomographic imaging is soundly understood and widely employed in medicine, nondestructive testing and security applications. However, it still offers many challenges when it comes to dynamic tomography. Indeed, in classical tomography, the target is usually assumed to be stationary during the data acquisition, but this is not a realistic model. Moreover, to ensure a lower X-ray radiation dose, only a sparse collection of measurements per time step is assumed to be available. With such a set up, we deal with a sparse data, dynamic tomography problem, which clearly calls for regularization, due to the loss of information in the data and the ongoing motion. In this paper, we propose a 3D variational formulation based on 3D shearlets, where the third dimension accounts for the motion in time, to reconstruct a moving 2D object. Results are presented for real measured data and compared against a 2D static model, in the case of fan-beam geometry. Results are preliminary but show that better reconstructions can be achieved when motion is taken into account.

[1]  E Somersalo,et al.  Statistical inversion for medical x-ray tomography with few radiographs: II. Application to dental radiology. , 2003, Physics in medicine and biology.

[2]  L. Zanni,et al.  New adaptive stepsize selections in gradient methods , 2008 .

[3]  Xiaosheng Zhuang,et al.  ShearLab: A Rational Design of a Digital Parabolic Scaling Algorithm , 2011, SIAM J. Imaging Sci..

[4]  Bernadette N. Hahn,et al.  Reconstruction of dynamic objects with affine deformations in computerized tomography , 2013 .

[5]  Gabriela Jager,et al.  MULTISCALE ANALYSIS OF MULTIVARIATE DATA , 2009 .

[6]  H Lantéri,et al.  A general method to devise maximum-likelihood signal restoration multiplicative algorithms with non-negativity constraints , 2001, Signal Process..

[7]  Andreas Hauptmann,et al.  Tomographic X-ray data of a lotus root filled with attenuating objects , 2016, 1609.07299.

[8]  S. Bonettini,et al.  New convergence results for the scaled gradient projection method , 2014, 1406.6601.

[9]  Wang-Q Lim,et al.  ShearLab 3D , 2014, 1402.5670.

[10]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .

[11]  E Somersalo,et al.  Statistical inversion for medical x-ray tomography with few radiographs: I. General theory. , 2003, Physics in medicine and biology.

[12]  Sebastian Bosse,et al.  A Haar wavelet-based perceptual similarity index for image quality assessment , 2016, Signal Process. Image Commun..

[13]  Bernadette N. Hahn,et al.  Null space and resolution in dynamic computerized tomography , 2016 .

[14]  Alexander Meaney Design and Construction of an X-ray Computed Tomography Imaging System , 2015 .

[15]  Frank Natterer,et al.  Mathematical methods in image reconstruction , 2001, SIAM monographs on mathematical modeling and computation.

[16]  E. Somersalo,et al.  Statistical inversion for X-ray tomography with few radiographs I: General theory , 2015 .

[17]  Matti Lassas,et al.  Dynamic multi-source X-ray tomography using a spacetime level set method , 2015, J. Comput. Phys..

[18]  L Xing,et al.  Motion correction for improved target localization with on-board cone-beam computed tomography , 2006, Physics in medicine and biology.

[19]  Andreas Hauptmann,et al.  A variational reconstruction method for undersampled dynamic x-ray tomography based on physical motion models , 2017, 1705.06079.

[20]  Federica Porta,et al.  Variable Metric Inexact Line-Search-Based Methods for Nonsmooth Optimization , 2015, SIAM J. Optim..

[21]  Erik L Ritman,et al.  Cardiac computed tomography imaging: a history and some future possibilities. , 2003, Cardiology clinics.