The boundary identity gives a numerical route to shape recovery from a Herglotz boundary identity. First recover an admissible density from measured far-field patterns, enforcing the supplied Herglotz pairing with an obstacle far field for available incident directions:
Then evaluate its Herglotz wave function throughout a search region and find a smooth positive radial profile satisfying
The same must work for every azimuth , expressing the assumed axial symmetry. In practice expand and in basis functions, enforce these equations at collocation points, and minimize the complex residual together with the far-field-data residual. The exclusion of an interior Dirichlet eigenvalue ensures uniqueness of the interior Dirichlet problem for the trial boundary data. Use regularization of an inverse problem to control measurement error and unstable density components; neither uniqueness nor numerical stability follows merely from writing the residual equations.
For rich incident-direction data, the conjugated pairing constraints form an adjoint far-field integral equation for . For only one fixed incidence, the displayed pairing is a single complex scalar constraint and cannot by itself determine an arbitrary density. Additional data or a joint constrained shape fit is then needed; the trace condition supplies the actual geometric information.
For large , the outgoing kernel has
Thus the surface representation of an obstacle far-field pattern, with , gives
Conjugating the Herglotz wave function gives . Interchanging the angular and obstacle-surface integrals yields the Herglotz pairing with an obstacle far field:
This contains only the scattered-field boundary traces and the Herglotz function, with the outward-obstacle normal convention fixing the sign.