Optical tomography (OT) recovers the cross-sectional distribution of optical parameters inside a highly scattering medium from information contained in measurements that are performed on the boundaries of the medium. The image reconstruction problem in OT can be considered as a large-scale optimization problem, in which an appropriately defined objective functional needs to be minimized. Most of earlier work is based on a forward model based iterative image reconstruction (MOBIIR) method. In this method, a Taylor series expansion of the forward propagation operator around the initial estimate, assumed to be close to the actual solution, is terminated at the first order term. The linearized perturbation equation is solved iteratively, re-estimating the first order term (or Jacobian) in each iteration, until a solution is reached. In this work we consider a nonlinear reconstruction problem, which has the second order term (Hessian) in addition to the first order. We show that in OT the Hessian is diagonally dominant and in this work an approximation involving the diagonal terms alone is used to formulate the nonlinear perturbation equation. This is solved using conjugate gradient search (CGS) without re-estimating either the Jacobian or the Hessian, resulting in reconstructions better than the original MOBIIR reconstruction. The computation time in this case is reduced by a factor of three.
[1]
A. Klose,et al.
Optical tomography using the time-independent equation of radiative transfer-Part 1: Forward model
,
2002
.
[2]
S. Arridge,et al.
Optical imaging in medicine: II. Modelling and reconstruction
,
1997,
Physics in medicine and biology.
[3]
S R Arridge,et al.
Optical imaging in medicine: I. Experimental techniques
,
1997,
Physics in medicine and biology.
[4]
M. Schweiger,et al.
A finite element approach for modeling photon transport in tissue.
,
1993,
Medical physics.
[5]
Alexander D. Klose,et al.
Gradient-based iterative image reconstruction scheme for time-resolved optical tomography
,
1999,
IEEE Transactions on Medical Imaging.