Premium
Iterative reconstruction scheme for optical tomography based on the equation of radiative transfer
Author(s) -
Klose Alexander D.,
Hielscher Andreas H.
Publication year - 1999
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1118/1.598661
Subject(s) - radiative transfer , iterative reconstruction , optical transfer function , optical tomography , diffuse optical imaging , detector , mathematics , diffusion equation , function (biology) , algorithm , iterative method , tomography , mathematical optimization , computer science , physics , optics , artificial intelligence , economy , evolutionary biology , economics , biology , service (business)
We report on the development of an iterative image reconstruction scheme for optical tomography that is based on the equation of radiative transfer. Unlike the commonly applied diffusion approximation, the equation of radiative transfer accurately describes the photon propagation in turbid media without any limiting assumptions regarding the optical properties. The reconstruction scheme consists of three major parts: (1) a forward model that predicts the detector readings based on solutions of the time‐independent radiative transfer equation, (2) an objective function that provides a measure of the differences between the detected and the predicted data, and (3) an updating scheme that uses the gradient of the objective function to perform a line minimization to get new guesses of the optical properties. The gradient is obtained by employing an adjoint differentiation scheme, which makes use of the structure of the finite‐difference discrete‐ordinate formulation of the transport forward model. Based on the new guess of the optical properties a new forward calculation is performed to get new detector predictions. The reconstruction process is completed when the minimum of the objective function is found within a defined error. To illustrate the performance of the code we present initial reconstruction results based on simulated data.