Optical tomography as a PDE-constrained optimization problem

Abstract

We report on the implementation of an augmented Lagrangian approach for solving the inverse problems in diffuse optical tomography (DOT). The forward model of light propagation is the radiative transport equation (RTE). The inverse problem is formulated as a minimization problem with the RTE being considered as an equality constraint on the set of 'optical properties-radiance' pairs. This approach allows the incorporation of the recently developed technique of PDE-constrained optimization, which has shown great promise in many applications that can be formulated as infinite-dimensional optimization problems. Compared to the traditional unconstrained optimization approaches for optical tomographic imaging where one solves several forward and adjoint problems at each optimization iteration, the method proposed in this work solves the forward and inverse problems simultaneously. We found in initial studies, using synthetic data, that the image reconstruction time can typically be reduced by a factor of 10 to 30, which depends on a combination of noise level, regularization parameter, mesh size, initial guess, optical properties and system geometry.

Recovering doping profiles in semiconductor devices with the Boltzmann–Poisson model

“…For example, one can adopt the method of Lagrange multipliers such as those in [2,29]. Here we solve the problem by transforming it into an unconstrained problem.…”

Non-contact fluorescence optical tomography with scanning patterned illumination

“…The solution of minimization problem (4) is determined as a stationary point of the Lagrangian [4,5] L…”

Recent Advances in Optical Tomographic Imaging

“…Recently our group has introduced the first transport theory based code that can be employed for frequencydomain data [3]. To achieve reasonable computation times, we have implemented a pdf-constrained optimization scheme, which leads to much increased convergence speeds [4]. In addition, the first transport-theory-based code for fluorescence imaging is now available [5], which allows tomographic reconstruction of the concentration of molecular probes in small animals with not-before-seen accuracy.…”