Recovery of optical cross-section perturbations in dense-scattering media by transport-theory-based imaging operators and steady-state simulated data

Abstract

We present a useful strategy for imaging perturbations of the macroscopic absorption cross section of dense-scattering media using steady-state light sources. A perturbation model based on transport theory is derived, and the inverse problem is simplified to a system of linear equations, WΔμ = ΔR, where W is the weight matrix, Δμ is a vector of the unknown perturbations, and ΔR is the vector of detector readings. Monte Carlo simulations compute the photon flux across the surfaces of phantoms containing simple or complex inhomogeneities. Calculation of the weight matrix is also based on the results of Monte Carlo simulations. Three reconstruction algorithms-conjugate gradient descent, projection onto convex sets, and the simultaneous algebraic reconstruction technique, with or without imposed positivity constraints-are used for image reconstruction. A rescaling technique that improves the conditioning of the weight matrix is also developed. Results show that the analysis of time-independent data by a perturbation model is capable of resolving the internal structure of a dense-scattering medium. Imposition of positivity constraints improves image quality at the cost of a reduced convergence rate. Use of the rescaling technique increases the initial rate of convergence, resulting in accurate images in a smaller number of iterations.

Non–invasive measurements of breast tissue optical properties using frequency–domain photon migration

“…Several groups are close to generating twoand three-dimensional optical property images using a variety of technical approaches (Arridge 1995 ;O'Leary et al 1995 ;Chang et al 1996 ;Huabei et al 1996). However, little is known about whether tumours can be identified by optically measured features.…”

Nonlinear effects of localized absorption perturbations on the light distribution in a turbid medium

“…Since this relation is an identity, for any j 0 it follows that ͑ T 0 ͒ ij ϭ t, j 0 ij ϩ 0 (5.D. 2) and therefore that . This is the correct definition for the absorption weight under the constraint that the transport cross section (or, equivalently, the diffusion coefficient) is fixed.…”

“…The primary approach we and others have taken in attempting to make the image-recovery problem of photon migration imaging more tractable [1][2][3][4][5][6][7][8] is to recast it as a discrete linear perturbation problem. These methods require specification of a reference medium whose properties are believed to differ as little as possible from those of the target medium.…”

Iterative reconstruction scheme for optical tomography based on the equation of radiative transfer

“…Wang and Ueno 60 established a theoretical procedure for estimating the two-dimensional distribution of the ground albedo of the earth by measurements made outside the atmosphere. Aronson et al 61 and Chang et al 62 used a Monte-Carlo solution of the adjoint transport equation. To locate discrete absorbers in highly scattering media from tomographic data, Aronson et al applied a backprojection method while Chang et al used a perturbation theory approach.…”