We introduce in this work a PDE-constrained approach to optical tomography that makes use of an all-atonce reduced Hessian Sequential Quadratic Programming (rSQP) scheme. The proposed scheme treats the forward and inverse variables independently, which makes it possible to update the radiation intensities and the optical coefficients simultaneously by solving the forward and inverse problems, all at once. We evaluate the performance of the proposed scheme with numerical and experimental data, and find that the rSQP scheme can reduce the computation time by a factor of 10 to 25, as compared to the commonly employed limited memory BFGS method.

Date Published: 23 February 2009
PDF: 9 pages
Proc. SPIE 7174, Optical Tomography and Spectroscopy of Tissue VIII, 717429 (23 February 2009); doi: 10.1117/12.809385

Published in SPIE Proceedings Vol. 7174:
Optical Tomography and Spectroscopy of Tissue VIII
Bruce J. Tromberg; Arjun G. Yodh; Mamoru Tamura; Eva M. Sevick-Muraca; Robert R. Alfano, Editor(s)