A Two-Stage Method for Inverse Acoustic Medium Scattering

We present a novel numerical method to the time-harmonic inverse medium scattering problem of recovering the refractive index from noisy near-field scattered data. The approach consists of two stages, one pruning step of detecting the scatterer support, and one resolution enhancing step with nonsmooth mixed regularization. The first step is strictly direct and of sampling type, and it faithfully detects the scatterer support. The second step is an innova-tive application of nonsmooth mixed regularization, and it accurately resolves the scatterer sizes as well as intensities. The nonsmooth model can be efficiently solved by a semi-smooth Newton-type method. Numerical results for two- and three-dimensional examples indicate that the new approach is accurate, computationally efficient, and robust with respect to data noise. Key words: inverse medium scattering problem, reconstruction algorithm, sampling method, mixed regularization, semi-smooth Newton method