Martin Schweiger et al 2005 Phys. Med. Biol. 50 2365 doi:10.1088/0031-9155/50/10/013
Martin Schweiger1, Simon R Arridge1 and Ilkka Nissilä2
Show affiliationsWe present a regularized Gauss–Newton method for solving the inverse problem of parameter reconstruction from boundary data in frequency-domain diffuse optical tomography. To avoid the explicit formation and inversion of the Hessian which is often prohibitively expensive in terms of memory resources and runtime for large-scale problems, we propose to solve the normal equation at each Newton step by means of an iterative Krylov method, which accesses the Hessian only in the form of matrix–vector products. This allows us to represent the Hessian implicitly by the Jacobian and regularization term. Further we introduce transformation strategies for data and parameter space to improve the reconstruction performance. We present simultaneous reconstructions of absorption and scattering distributions using this method for a simulated test case and experimental phantom data.
Issue 10 (21 May 2005)
Received 21 January 2005, in final form 30 March 2005
Published 5 May 2005
Martin Schweiger et al 2005 Phys. Med. Biol. 50 2365
David D Hsu and David B Graves 2003 J. Phys. D: Appl. Phys. 36 2898
Abhinav Kranti and G Alastair Armstrong 2006 Semicond. Sci. Technol. 21 1563
Yungui Gong and Anzhong Wang 2006 Class. Quantum Grav. 23 3419
G C Das and D Chakravorty 1982 J. Phys. D: Appl. Phys. 15 2255
D R Bowler et al 2002 J. Phys.: Condens. Matter 14 2781
G S McDonald et al 1997 J. Phys. B: At. Mol. Opt. Phys. 30 L719
Hans Havlicek et al 2001 J. Phys. A: Math. Gen. 34 3071
Catherine C Berry 2009 J. Phys. D: Appl. Phys. 42 224003
T Ito et al 2006 J. Phys.: Condens. Matter 18 6009