Optical tomography as a PDE-constrained optimization problem

Gassan S. Abdoulaev, Kui Ren, Andreas H. Hielscher

    Research output: Contribution to journalArticlepeer-review

    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.

    Original languageEnglish (US)
    Pages (from-to)1507-1530
    Number of pages24
    JournalInverse Problems
    Volume21
    Issue number5
    DOIs
    StatePublished - Oct 1 2005

    ASJC Scopus subject areas

    • Theoretical Computer Science
    • Signal Processing
    • Mathematical Physics
    • Computer Science Applications
    • Applied Mathematics

    Fingerprint Dive into the research topics of 'Optical tomography as a PDE-constrained optimization problem'. Together they form a unique fingerprint.

    Cite this