TY - GEN

T1 - Solution of the Perturbation Equation in Optical Tomography Using Weight Functions as a Transform Basis

AU - Lin, Erh Ya

AU - Wang, Yao

AU - Pei, Yaling

AU - Barbour, Randall L.

N1 - Funding Information:
This work is supported in part by the NIH under grant R01-CA66184-01A2.
Publisher Copyright:
© 1998 Optical Society of America.

PY - 1998

Y1 - 1998

N2 - This paper describes a new inverse solver for optical tomography. As with prior stndies, we employ an iterative pertnrbation approach, which at each iteration requires the solution of a forward problem and an inverse problem. The inverse problem involves the solution of a linear pertnrbation equation, which is often severely underdetermined. To overcome this problem, we propose to represent the unknown image of optical properties by a set of linearly independent basis functions, with the number of basis functions being equal to or less than the number of independent detector readings. The accuracy of the solution depends on the choice of the basis. We have explored the use of the weight functions associated with different source and detector pairs (i.e. the rows in the weight matrix of the pertnrbation equation) as the basis functions. By choosing those source and detector pairs which have uncorrelated weight functions, the inverse problem is transformed into a well-posed, uniquely determined problem. The system matrix in the transformed representation has a dimension significantly smaller than the original matrix, so that it is feasible to perform the inversion using singular value decomposition (SVD). This new method has been integrated with a previously reported forward solver, and applied to data generated from numerical simulations using diffusion approximation. Compared to the Conjugate Gradient Descent (CGD) method used in previously reported studies, the new method takes substantially less computation time, while providing equal, if not better, image reconstruction quality at similar noise levels.

AB - This paper describes a new inverse solver for optical tomography. As with prior stndies, we employ an iterative pertnrbation approach, which at each iteration requires the solution of a forward problem and an inverse problem. The inverse problem involves the solution of a linear pertnrbation equation, which is often severely underdetermined. To overcome this problem, we propose to represent the unknown image of optical properties by a set of linearly independent basis functions, with the number of basis functions being equal to or less than the number of independent detector readings. The accuracy of the solution depends on the choice of the basis. We have explored the use of the weight functions associated with different source and detector pairs (i.e. the rows in the weight matrix of the pertnrbation equation) as the basis functions. By choosing those source and detector pairs which have uncorrelated weight functions, the inverse problem is transformed into a well-posed, uniquely determined problem. The system matrix in the transformed representation has a dimension significantly smaller than the original matrix, so that it is feasible to perform the inversion using singular value decomposition (SVD). This new method has been integrated with a previously reported forward solver, and applied to data generated from numerical simulations using diffusion approximation. Compared to the Conjugate Gradient Descent (CGD) method used in previously reported studies, the new method takes substantially less computation time, while providing equal, if not better, image reconstruction quality at similar noise levels.

KW - (100.3010) Image reconstruction techniques

KW - (100.3080) Infra-red imaging

UR - http://www.scopus.com/inward/record.url?scp=0344513568&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0344513568&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:0344513568

T3 - Optics InfoBase Conference Papers

SP - 197

EP - 202

BT - Advances in Optical Imaging and Photon Migration, AOIPM 1998

PB - Optica Publishing Group (formerly OSA)

T2 - Advances in Optical Imaging and Photon Migration, AOIPM 1998

Y2 - 8 March 1998

ER -