Compressed sensing based interior tomography

While conventional wisdom is that the interior problem does not have a unique solution, by analytic continuation we recently showed that the interior problem can be uniquely and stably solved if we have a known sub-region inside a region of interest (ROI). However, such a known sub-region is not always readily available, and it is even impossible to find in some cases. Based on compressed sensing theory, here we prove that if an object under reconstruction is essentially piecewise constant, a local ROI can be exactly and stably reconstructed via the total variation minimization. Because many objects in computed tomography (CT) applications can be approximately modeled as piecewise constant, our approach is practically useful and suggests a new research direction for interior tomography. To illustrate the merits of our finding, we develop an iterative interior reconstruction algorithm that minimizes the total variation of a reconstructed image and evaluate the performance in numerical simulation.

[1]  D L Parker,et al.  Optimal short scan convolution reconstruction for fanbeam CT. , 1982, Medical physics.

[2]  D. Parker Optimal short scan convolution reconstruction for fan beam CT , 1982 .

[3]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

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

[5]  A survey of the role of the UK physicist in nuclear medicine: a report of a joint working group of the British Institute of Radiology, British Nuclear Medicine Society, and the Institute of Physics and Engineering in Medicine , 2003, Nuclear medicine communications.

[6]  M. Jiang,et al.  Uniqueness theorems in bioluminescence tomography. , 2004, Medical physics.

[7]  M. Jiang,et al.  Ordered-subset simultaneous algebraic reconstruction techniques (OS-SART) , 2004 .

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

[9]  Hengyong Yu,et al.  A General Formula for Fan-Beam Lambda Tomography , 2007, Int. J. Biomed. Imaging.

[10]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[11]  Daniel Day,et al.  Mechanical Strains Induced in Osteoblasts by Use of Point Femtosecond Laser Targeting , 2006, Int. J. Biomed. Imaging.

[12]  Hiroyuki Kudo,et al.  Truncated Hilbert transform and image reconstruction from limited tomographic data , 2006 .

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

[14]  Hengyong Yu,et al.  A General Local Reconstruction Approach Based on a Truncated Hilbert Transform , 2007, Int. J. Biomed. Imaging.

[15]  M. Defrise,et al.  Tiny a priori knowledge solves the interior problem in computed tomography , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[16]  Hengyong Yu,et al.  Exact Interior Reconstruction with Cone-Beam CT , 2008, Int. J. Biomed. Imaging.

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

[18]  M. Defrise,et al.  Solving the interior problem of computed tomography using a priori knowledge , 2008, Inverse problems.

[19]  Hengyong Yu,et al.  Interior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition. , 2008, Journal of X-ray science and technology.

[20]  Hengyong Yu,et al.  Exact Interior Reconstruction from Truncated Limited-Angle Projection Data , 2008, Int. J. Biomed. Imaging.

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