A surrogate-Bayesian algorithm for scatterer shape identification from phaseless data

TL;DR

This paper introduces a Bayesian-based method with a neural network surrogate for fast, uncertainty-aware reconstruction of scatterer shapes from phaseless far-field data, significantly reducing computational costs.

Contribution

The work develops an IPINN surrogate model that efficiently approximates the Helmholtz equation solutions, enabling rapid Bayesian shape reconstruction from intensity-only data.

Findings

Achieves several orders of magnitude speed-up in computations

Accurately reconstructs scatterer shapes with quantified uncertainties

Demonstrates effectiveness through numerical experiments

Abstract

This work addresses the reconstruction of a scatterer's shape from phaseless far field-intensity data arising from multiple incident waves interacting with the scatterer. We formulate the reconstruction as a statistical inverse scattering problem and adopt a Bayesian inference framework, which can readily be used to compute statistical moments for quantification of uncertainties in the shape reconstruction that arise from noise in the data due to measurement constraints. The shape of the scatterer is represented by a spline-based prior, with Bayesian parameters defined at the spline's knots. To efficiently evaluate the Bayesian likelihood across thousands of sampling points, we develop the intensity property inspired neural network (IPINN) surrogate. This surrogate incorporates the Helmholtz equation in the unbounded domain, exterior to each sampled scatterer, along with the radiation condition at infinity, enabling fast and accurate simulation of the acoustic far-field intensity. Importantly, the IPINN surrogate is trained independently of the observed data and requires only a single incident wave for training. We demonstrate that this surrogate approach yields a speed-up of several orders of magnitude. The resulting IPINN-Bayesian framework offers an efficient solution for shape reconstruction in unbounded domains with multiple incident wave boundary conditions, while exactly enforcing the radiation condition. Numerical experiments confirm the efficiency and effectiveness of the proposed algorithm.