Compressed-sensing-based content-driven hierarchical reconstruction: Theory and application to C-arm cone-beam tomography.
OBJECTIVE
This paper addresses the reconstruction of x-ray cone-beam computed tomography (CBCT) for interventional C-arm systems. Subsampling of CBCT is a significant issue with C-arms due to their slow rotation and to the low frame rate of their flat panel x-ray detectors. The aim of this work is to propose a novel method able to handle the subsampling artifacts generally observed with analytical reconstruction, through a content-driven hierarchical reconstruction based on compressed sensing.


METHODS
The central idea is to proceed with a hierarchical method where the most salient features (high intensities or gradients) are reconstructed first to reduce the artifacts these features induce. These artifacts are addressed first because their presence contaminates less salient features. Several hierarchical schemes aiming at streak artifacts reduction are introduced for C-arm CBCT: the empirical orthogonal matching pursuit approach with the ℓ0 pseudonorm for reconstructing sparse vessels; a convex variant using homotopy with the ℓ1-norm constraint of compressed sensing, for reconstructing sparse vessels over a nonsparse background; homotopy with total variation (TV); and a novel empirical extension to nonlinear diffusion (NLD). Such principles are implemented with penalized iterative filtered backprojection algorithms. For soft-tissue imaging, the authors compare the use of TV and NLD filters as sparsity constraints, both optimized with the alternating direction method of multipliers, using a threshold for TV and a nonlinear weighting for NLD.


RESULTS
The authors show on simulated data that their approach provides fast convergence to good approximations of the solution of the TV-constrained minimization problem introduced by the compressed sensing theory. Using C-arm CBCT clinical data, the authors show that both TV and NLD can deliver improved image quality by reducing streaks.


CONCLUSIONS
A flexible compressed-sensing-based algorithmic approach is proposed that is able to accommodate for a wide range of constraints. It is successfully applied to C-arm CBCT images that may not be so well approximated by piecewise constant functions.
Image Processing and Clinical Applications Laboratory, GE Healthcare, Buc F-78533, France; Laboratoire des Signaux et Systèmes, CentraleSupélec, Gif-sur-Yvette F-91192, France; Center for Visual Computing, CentraleSupélec, Châtenay-Malabry F-92295, France; and INRIA, Orsay F-91893, France.
Hélène
H
Langet
Image Processing and Clinical Applications Laboratory, GE Healthcare, Buc F-78533, France.
Cyril
C
Riddell
Image Processing and Clinical Applications Laboratory, GE Healthcare, Buc F-78533, France.
Aymeric
A
Reshef
Image Processing and Clinical Applications Laboratory, GE Healthcare, Buc F-78533, France.
Yves
Y
Trousset
Laboratoire des Signaux et Systèmes, CentraleSupélec, Gif-sur-Yvette F-91192, France.
Arthur
A
Tenenhaus
Laboratoire des Signaux et Systèmes, CentraleSupélec, Gif-sur-Yvette F-91192, France.
Elisabeth
E
Lahalle
Laboratoire des Signaux et Systèmes, CentraleSupélec, Gif-sur-Yvette F-91192, France.
Gilles
G
Fleury
Center for Visual Computing, CentraleSupélec, Châtenay-Malabry F-92295, France and INRIA, Orsay F-91893, France.
Nikos
N
Paragios
D000465:Algorithms; D016477:Artifacts; D054893:Cone-Beam Computed Tomography / Q000379:methods*; D004058:Diffusion; D006801:Humans; D007091:Image Processing, Computer-Assisted / Q000379:methods*