A sparsity regularization and total variation based computational framework for the inverse medium problem in scattering
Autor: | Florian Bürgel, Kamil S. Kazimierski, Armin Lechleiter |
---|---|
Rok vydání: | 2017 |
Předmět: |
Numerical Analysis
Mathematical optimization Physics and Astronomy (miscellaneous) Discretization Helmholtz equation Scattering Applied Mathematics Inverse 010103 numerical & computational mathematics 02 engineering and technology Backus–Gilbert method 01 natural sciences Regularization (mathematics) Computer Science Applications Computational Mathematics Modeling and Simulation Inverse scattering problem 0202 electrical engineering electronic engineering information engineering Applied mathematics 020201 artificial intelligence & image processing Point (geometry) 0101 mathematics Mathematics |
Zdroj: | Journal of Computational Physics. 339:1-30 |
ISSN: | 0021-9991 |
DOI: | 10.1016/j.jcp.2017.03.011 |
Popis: | We present a fast computational framework for the inverse medium problem in scattering, i.e. we look at discretization, reconstruction and numerical performance. The Helmholtz equation in two and three dimensions is used as a physical model of scattering including point sources and plane waves as incident fields as well as near and far field measurements. For the reconstruction of the medium, we set up a rapid variational regularization scheme and indicate favorable choices of the various parameters. The underlying paradigm is, roughly speaking, to minimize the discrepancy between the reconstruction and measured data while, at the same time, taking into account various structural a-priori information via suitable penalty terms. In particular, the involved penalty terms are designed to promote information expected in real-world environments. To this end, a combination of sparsity promoting terms, total variation, and physical bounds of the inhomogeneous medium, e.g. positivity constraints, is employed in the regularization penalty. A primal-dual algorithm is used to solve the minimization problem related to the variational regularization. The computational feasibility, performance and efficiency of the proposed approach is demonstrated for synthetic as well as experimentally measured data. |
Databáze: | OpenAIRE |
Externí odkaz: |