Restoration of Missing Data in Limited Angle Tomography Based on Consistency Conditions

Yixing Huang, Oliver Taubmann, Xiaolin Huang, Joachim Hornegger, Guenter Lauritsch, Andreas Maier


Published in:Fully3D 2017 Proceedings


In limited angle tomography, only a limited angular range of data is acquired and consequently a double wedgeshaped region in the frequency domain representation of the imaged object is missing. Hence, streak artifacts occur. To restore the missing data, we perform a regression and an image fusion in sinogram domain and frequency domain, respectively. We first convert the sinogram restoration problem into a regression problem based on the Helgason-Ludwig consistency conditions. Due to the severe ill-posedness of the problem, regression only partially recovers the correct frequency components, especially lower frequency components, and will introduce erroneous ones, particularly higher frequencies. Bilateral filtering is utilized to retain the most prominent high frequency components and suppress erroneous ones. A fusion of the filtered image and the image reconstructed from the limited angle sinogram is performed afterwards in the frequency domain. The proposed method is evaluated on the Shepp-Logan phantom, for which the root-mean-square error of the reconstructed image decreases from 310 HU to 136 HU.
Yixing Huang
Pattern Recognition Lab, Friedrich-Alexander-University Erlangen-Nuremberg
Oliver Taubmann
Pattern Recognition Lab, Friedrich-Alexander-University Erlangen-Nuremberg
Xiaolin Huang
Institute of Image Processing and Pattern Recognition, Shanghai Jiao Tong University
Joachim Hornegger
Friedrich-Alexander-University Erlangen-Nuremberg
Guenter Lauritsch
Siemens Healthcare GmbH
Andreas Maier
Pattern Recognition Lab, Friedrich-Alexander-University Erlangen-Nuremberg
  1. M. E. Davison, “The ill-conditioned nature of the limited angle tomography problem,” SIAM Journal on Applied Mathematics, vol. 43, no. 2,pp. 428–448, 1983.
  2. A. K. Louis, “Incomplete data problems in X-ray computerized tomography: I. singular value decomposition of the limited angle transform,” Numerische Mathematik, vol. 48, no. 3, pp. 251–262, 1986.
  3. E. Y. Sidky and X. Pan, “Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization,” Physics in Medicine and Biology, vol. 53, no. 17, p. 4777, 2008.
  4. L. Ritschl, F. Bergner, C. Fleischmann, and M. Kachelrieß, “Improved total variation-based CT image reconstruction applied to clinical data,” Physics in Medicine and Biology, vol. 56, no. 6, p. 1545, 2011.
  5. Z. Chen, X. Jin, L. Li, and G. Wang, “A limited-angle CT reconstruction method based on anisotropic TV minimization,” Physics in Medicine and Biology, vol. 58, no. 7, p. 2119, 2013.
  6. Y. Huang, O. Taubmann, X. Huang, V. Haase, G. Lauritsch, and A. Maier, “A New Scale Space Total Variation Algorithm for Limited Angle Tomography,” in CT-Meeting 2016 Proceedings (The 4th International Meeting on Image Formation in X-Ray Computed Tomography),M. Kachelrieß, Ed., 2016, pp. 149–152.
  7. C. Riess, M. Berger, H. Wu, M. Manhart, R. Fahrig, and A. Maier, “TV or not TV? That is the Question,” in Fully 3D Image Reconstruction in Radiology and Nuclear Medicine, 2013, pp. 341–344.
  8. T. Wurfl, F. C. Ghesu, V. Christlein, and A. Maier, “Deep Learning ¨ Computed Tomography,” in Medical Image Computing and ComputerAssisted Intervention MICCAI 2016, Springer, Ed., vol. 3, 2016, pp. 432–440.
  9. A. Aichert, M. Berger, J. Wang, N. Maass, A. Doerfler, J. Hornegger, and A. Maier, “Epipolar consistency in transmission imaging,” IEEE Transactions on Medical Imaging, vol. 34, no. 10, pp. 1–15, 2015.
  10. S. K. Patch, “Consistency conditions upon 3D CT data and the wave equation,” Physics in Medicine and Biology, vol. 47, no. 15, p. 2637, 2002.
  11. R. Clackdoyle and L. Desbat, “Data consistency conditions for truncated fanbeam and parallel projections,” Medical Physics, vol. 42, no. 2, pp. 831–845, 2015.
  12. S. Helgason, “The radon transform on euclidean spaces, compact twopoint homogeneous spaces and grassmann manifolds,” Acta Mathematica, vol. 113, no. 1, pp. 153–180, 1965.
  13. D. Ludwig, “The radon transform on euclidean space,” Communications on Pure and Applied Mathematics, vol. 19, no. 1, pp. 49–81, 1966.
  14. S. R. Deans, The Radon transform and some of its applications. Courier Corporation, 2007.
  15. A. K. Louis, “Approximation of the radon transform from samples inlimited range,” in Mathematical Aspects of Computerized Tomography.Springer, 1981, pp. 127–139.
  16. A. S. Willsky and J. L. Prince, “Constrained sinogram restoration forlimited-angle tomography,” Optical Engineering, vol. 29, no. 5, pp. 535–544, 1990.
  17. H. Kudo and T. Saito, “Sinogram recovery with the method of convexprojections for limited-data reconstruction in computed tomography,”Journal of the Optical Society of America A, vol. 8, no. 7, pp. 1148–1160, 1991.
  18. G. B. Arfken and H. J. Weber, Mathematical methods for physicistsinternational student edition. Academic press, 2005.
  19. R. Tibshirani, “Regression shrinkage and selection via the lasso,”Journal of the Royal Statistical Society. Series B (Methodological), pp.267–288, 1996.
  20. I. Daubechies, M. Defrise, and C. De Mol, “An iterative thresholdingalgorithm for linear inverse problems with a sparsity constraint,” Communications on pure and applied mathematics, vol. 57, no. 11, pp. 1413–1457, 2004.
  21. C. Tomasi and R. Manduchi, “Bilateral filtering for gray and colorimages,” in Computer Vision, 1998. Sixth International Conference on.IEEE, 1998, pp. 839–846.
  22. A. Maier, H. Hofmann, M. Berger, P. Fischer, C. Schwemmer, H. Wu,K. Muller, J. Hornegger, J. Choi, C. Riess, A. Keil, and R. Fahrig, ¨ “CONRAD - a software framework for cone-beam imaging in radiology,” Medical Physics, vol. 40, no. 11, p. 111914, 2013.