Reconstruction of a Piecewise Smooth Absorption Coefficient by an Acousto-Optic Process

Abstract

The aim of this paper is to tackle the nonlinear optical reconstruction problem. Given a set of acousto-optic measurements, we develop a mathematical framework for the reconstruction problem in the case where the optical absorption distribution is supposed to be a perturbation of a piecewise constant function. Analyzing the acousto-optic measurements, we prove that the optical absorption coefficient satisfies, in the sense of distributions, a new equation. For doing so, we introduce a weak Helmholtz decomposition and interpret in a weak sense the cross-correlation measurements using the spherical Radon transform. We next show how to find an initial guess for the unknown coefficient. Finally, in order to construct the true coefficient we provide a Landweber type iteration and prove that the resulting sequence converges to the solution of the system constituted by the optical diffusion equation and the new equation mentioned above. Our results in this paper generalize the acousto-optic process proposed in [5] for piecewise smooth optical absorption distributions.

Keywords:

• Acousto-optic inverse problem
• Helmholtz decomposition
• Landweber iteration
• Piecewise smooth functions
• Reconstruction
• Spherical Radon transform
• Stability

Mathematics Subject Classification (MSC2000):

• 35R30
• 35B30

Acknowledgments

The authors are very grateful to the reviewers for their comments and suggestions to improve the presentation of the paper.

This work was supported by the ERC Advanced Grant Project MULTIMOD–267184.