Fréchet derivative with respect to the shape of a strongly convex nonscattering region in optical tomography

The aim of optical tomography is to reconstruct the optical properties inside a physical body, e.g. a neonatal head, by illuminating it with near-infrared light and measuring the outward flux of photons on the object boundary. Because a brain consists of strongly scattering tissue with imbedded cavities filled by weakly scattering cerebrospinal fluid, propagation of near-infrared photons in the human head can be treated by combining the diffusion approximation of the radiative transfer equation with geometrical optics to obtain the radiosity-diffusion forward model of optical tomography. At the moment, a disadvantage with the radiosity-diffusion model is that the locations of the transparent cavities must be known in advance in order to be able to reconstruct the physiologically interesting quantities, i.e., the absorption and the scatter in the strongly scattering brain tissue. In this work we show that the boundary measurement map of optical tomography is Frechet differentiable with respect to the shape of a strongly convex nonscattering region. Using this result, we introduce a numerical algorithm for approximating an unknown nonscattering cavity by a ball if the background diffuse optical properties of the object are known. The functionality of the method is demonstrated through two-dimensional numerical experiments.