An optical tomography reconstruction algorithm with the Discontinuous Galerkin method

In optical tomography, the optical properties are recovered through an iterative scheme, which consists in minimizing the errors between the measurements and the predications of a forward model. Traditionally, finite volumes or continuous finite elements formulations of light transport are used as a forward model for the predictions. we have introduced an integral form of the objective function that takes into account the surface of the detectors, and thus making compatible the inverse approach with all finite elements formulations (continuous and discontinuous). This present paper illustrates this novel approach by developing a Discontinuous Galerkin formulation as a forward model. Numerical tests are performed to gauge the accuracy of the method with a gradient-based algorithm. The results show that the reconstruction is accurate and can be affected by noise on the measurements as expected. A filtering of the gradient improves the quality and the accuracy of the reconstruction. Copyright \ua9 2011 by ASME.Peer reviewed: YesNRC publication: Ye