ePrints@IIScePrints@IISc Home | About | Browse | Latest Additions | Advanced Search | Contact | Help

Diffuse optical tomography through solving a system of quadratic equations:theory and simulations

Kanmani, B and Vasu, R M (2006) Diffuse optical tomography through solving a system of quadratic equations:theory and simulations. In: Physics In Medicine And Biology, 51 (4). 981 -998.

[img] PDF
8.pdf - Published Version
Restricted to Registered users only

Download (1MB) | Request a copy
Official URL: http://www.iop.org/EJ/article/0031-9155/51/4/015/p...


This paper discusses the iterative solution of the nonlinear problem of optical tomography. In the established forward model-based iterative image reconstruction (MOBIIR) method a linear perturbation equation containing the first derivative of the forward operator is solved to obtain the update vector for the optical properties. In MOBIIR, the perturbation equation is updated by recomputing the first derivative after each update of the optical properties. In the method presented here a nonlinear perturbation equation, containing terms up to the second derivative, is used to iteratively solve for the optical property updates. Through this modification, reconstructions with reasonable contrast recovery and accuracy are obtained without the need for updating the perturbation equation and therefore eliminating the outer iteration of the Usual MOBIIR algorithm. To improve the performance of the algorithm the outer iteration is reintroduced in which the perturbation equation is recomputed without re-estimating the derivatives and with only updated computed data. The system of quadratic equations is solved using either a modified conjugate gradient descent scheme or a two-step linearized predictor-corrector scheme. A quick method employing the adjoint of the forward operator is used to estimate the derivatives. By solving the nonlinear perturbation equation it is shown that the iterative scheme is able to recover large contrast variations in absorption coefficient with improved noise tolerance in data. This ability has not been possible so far with linear algorithms. This is demonstrated by presenting results of numerical simulations from objects with inhomogeneous inclusions in absorption coefficient with different contrasts and shapes.

Item Type: Journal Article
Publication: Physics In Medicine And Biology
Publisher: Institute of Physics Publishing
Additional Information: Copyright of this article belongs to Institute of Physics Publishing
Department/Centre: Division of Physical & Mathematical Sciences > Instrumentation Appiled Physics
Date Deposited: 20 Mar 2009 11:43
Last Modified: 19 Sep 2010 04:55
URI: http://eprints.iisc.ac.in/id/eprint/17070

Actions (login required)

View Item View Item