### 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 language | English |
---|---|

Title of host publication | Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013 |

Editors | Hamid R. Arabnia, Leonidas Deligiannidis, Joan Lu, Fernando G. Tinetti, Jane You, George Jandieri, Gerald Schaefer, Ashu M. G. Solo, Vladimir Volkov |

Publisher | CSREA Press |

Pages | 513-517 |

Number of pages | 5 |

ISBN (Electronic) | 1601322534, 9781601322531 |

Publication status | Published - 1 Jan 2013 |

Event | 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013, at WORLDCOMP 2013 - Las Vegas, USA United States Duration: 22 Jul 2013 → 25 Jul 2013 |

### Publication series

Name | Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013 |
---|---|

Volume | 2 |

### Conference

Conference | 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013, at WORLDCOMP 2013 |
---|---|

Country | USA United States |

City | Las Vegas |

Period | 22/07/13 → 25/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

*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.

**A hybrid regularization algorithm for high contrast tomographic image reconstruction.** / Rahmati, Peyman; Soleimani, Manuchehr; Adler, Andy.

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution

*Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013.*Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013, vol. 2, CSREA Press, pp. 513-517, 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013, at WORLDCOMP 2013, Las Vegas, USA United States, 22/07/13.

}

TY - GEN

T1 - A hybrid regularization algorithm for high contrast tomographic image reconstruction

AU - Rahmati, Peyman

AU - Soleimani, Manuchehr

AU - Adler, Andy

PY - 2013/1/1

Y1 - 2013/1/1

N2 - 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.

AB - 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.

KW - Electrical impedance tomography

KW - Inverse problem

KW - Level set

KW - Primal-dual interior point method

KW - Sum of absolute values (L1 norms)

KW - Tomographic image reconstruction

UR - http://www.scopus.com/inward/record.url?scp=85072925417&partnerID=8YFLogxK

M3 - Conference contribution

T3 - Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013

SP - 513

EP - 517

BT - Proceedings of the 2013 International Conference on Image Processing, Computer Vision, and Pattern Recognition, IPCV 2013

A2 - Arabnia, Hamid R.

A2 - Deligiannidis, Leonidas

A2 - Lu, Joan

A2 - Tinetti, Fernando G.

A2 - You, Jane

A2 - Jandieri, George

A2 - Schaefer, Gerald

A2 - Solo, Ashu M. G.

A2 - Volkov, Vladimir

PB - CSREA Press

ER -