A hybrid regularization algorithm for high contrast tomographic image reconstruction

Peyman Rahmati, Manuchehr Soleimani, Andy Adler

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

The common Level set based reconstruction method (LSRM) is applied to solve a piecewise constant inverse problem using one level set function and considers two different conductivity quantities for the background and the inclusion (two phases inclusion). The more the number of the piecewise constant conductivities in the medium, the higher the calculation effort of the LSRM, using multiple level set functions, will be. The assumption of piecewise constant conductivities (coefficients) is to discriminate between two regions with sharp conductivity interface; however, it may not be a realistic assumption when there are smooth conductivity gradients inside each region as well. In this paper, we propose a hybrid regularization method (HRM), which is a two steps solution, to solve ill-posed, nonlinear inverse problem with smooth conductivity transitions. The first step of this hybrid inversion framework plays the role of an initializing procedure for the second step, and acts in a similar way as a source type inversion method. In the first stage, the LSRM with one level set function is applied to determine the region of interest, which is defined as the region with sharpest interface. Then in the second stage, an inverse solver with penalty terms based on sum of absolute values (L1 norms), which are highly robust against measurement errors, is applied to reconstruct the conductivity changes inside the determined ROI. The generated forward solution in the final iteration of the level set is fed to the second stage where the L1 norms based penalty terms are minimized using primal-dual interior point method (PDIPM). The PDIPM has been shown to be effective in minimizing the L1 norms. The reconstructed images with the proposed HRM maintains the edge information as well as the smooth conductivity variations, a trait absent in all previously established level set based reconstruction method. The integration of the LSRM and the PDIPM can generate less noisy reconstructed images when comparing either with the reconstruction results of the PDIPM or with those of squared error based reconstruction methods, such as Gauss-Newton (GN) method. Our proposed HRM is tested on a circular 2D phantom with either sharp conductivity gradients (piecewise constant coefficients) or smooth conductivity transitions (smooth coefficients). We show that the proposed HRM maintains sharp edges and is robust against the measurement noise.

Original languageEnglish
Title of host publicationProceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013
EditorsHamid R. Arabnia, Leonidas Deligiannidis, Joan Lu, Fernando G. Tinetti, Jane You, George Jandieri, Gerald Schaefer, Ashu M. G. Solo, Vladimir Volkov
PublisherCSREA Press
Pages513-517
Number of pages5
ISBN (Electronic)1601322534, 9781601322531
Publication statusPublished - 1 Jan 2013
Event2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013, at WORLDCOMP 2013 - Las Vegas, USA United States
Duration: 22 Jul 201325 Jul 2013

Publication series

NameProceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013
Volume2

Conference

Conference2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013, at WORLDCOMP 2013
CountryUSA United States
CityLas Vegas
Period22/07/1325/07/13

Keywords

  • Electrical impedance tomography
  • Inverse problem
  • Level set
  • Primal-dual interior point method
  • Sum of absolute values (L1 norms)
  • Tomographic image reconstruction

ASJC Scopus subject areas

  • Computer Graphics and Computer-Aided Design
  • Computer Vision and Pattern Recognition

Cite this

Rahmati, P., Soleimani, M., & Adler, A. (2013). A hybrid regularization algorithm for high contrast tomographic image reconstruction. In H. R. Arabnia, L. Deligiannidis, J. Lu, F. G. Tinetti, J. You, G. Jandieri, G. Schaefer, A. M. G. Solo, ... V. Volkov (Eds.), Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013 (pp. 513-517). (Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013; Vol. 2). CSREA Press.