3D Image Reconstruction from X-Ray Measurements with Overlap

3D image reconstruction from a set of X-ray projections is an important image reconstruction problem, with applications in medical imaging, industrial inspection and airport security. The innovation of X-ray emitter arrays allows for a novel type of X-ray scanners with multiple simultaneously emitting sources. However, two or more sources emitting at the same time can yield measurements from overlapping rays, imposing a new type of image reconstruction problem based on nonlinear constraints. Using traditional linear reconstruction methods, respective scanner geometries have to be implemented such that no rays overlap, which severely restricts the scanner design. We derive a new type of 3D image reconstruction model with nonlinear constraints, based on measurements with overlapping X-rays. Further, we show that the arising optimization problem is partially convex, and present an algorithm to solve it. Experiments show highly improved image reconstruction results from both simulated and real-world measurements.

[1]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[2]  Francesc Moreno-Noguer,et al.  Simultaneous correspondence and non-rigid 3D reconstruction of the coronary tree from single X-ray images , 2011, 2011 International Conference on Computer Vision.

[3]  Massimo Fornasier,et al.  Quasi-linear Compressed Sensing , 2013, Multiscale Model. Simul..

[4]  Yuan Cheng,et al.  Rectangular Fixed-Gantry CT Prototype: Combining CNT X-Ray Sources and Accelerated Compressed Sensing-Based Reconstruction , 2014, IEEE Access.

[5]  Xiaodong Li,et al.  Sparse Signal Recovery from Quadratic Measurements via Convex Programming , 2012, SIAM J. Math. Anal..

[6]  Daniel Cremers,et al.  Integration of Multiview Stereo and Silhouettes Via Convex Functionals on Convex Domains , 2008, ECCV.

[7]  John Amanatides,et al.  A Fast Voxel Traversal Algorithm for Ray Tracing , 1987, Eurographics.

[8]  M. Lustig,et al.  Compressed Sensing MRI , 2008, IEEE Signal Processing Magazine.

[9]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[10]  Hyojin Kim,et al.  A Randomized Ensemble Approach to Industrial CT Segmentation , 2015, 2015 IEEE International Conference on Computer Vision (ICCV).

[11]  Luminita A. Vese,et al.  Expectation maximization and total variation-based model for computed tomography reconstruction from undersampled data , 2011, Medical Imaging.

[12]  Deanna Needell,et al.  CoSaMP: Iterative signal recovery from incomplete and inaccurate samples , 2008, ArXiv.

[13]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[14]  Jun Chen,et al.  Transmission type flat-panel X-ray source using ZnO nanowire field emitters , 2015 .

[15]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[16]  Patrick L. Combettes,et al.  Signal Recovery by Proximal Forward-Backward Splitting , 2005, Multiscale Model. Simul..

[17]  Stephen P. Boyd,et al.  Proximal Algorithms , 2013, Found. Trends Optim..

[18]  Shiqian Ma,et al.  An efficient algorithm for compressed MR imaging using total variation and wavelets , 2008, 2008 IEEE Conference on Computer Vision and Pattern Recognition.

[19]  I. Johnstone,et al.  Minimax estimation via wavelet shrinkage , 1998 .

[20]  Thomas Blumensath,et al.  Compressed Sensing With Nonlinear Observations and Related Nonlinear Optimization Problems , 2012, IEEE Transactions on Information Theory.