Abstract
A reconstruction algorithm for bioluminescence tomography (BLT) has been developed. The algorithm numerically calculates the Green’s function at different wavelengths using the diffusion equation and finite element method. The optical properties used in calculating the Green’s function are reconstructed using diffuse optical tomography (DOT) and assuming anatomical information is provided by x-ray computed tomography or other methods. A symmetric system of equations is formed using the Green’s function and the measured light fluence rate and the resulting eigenvalue problem is solved to get the eigenvectors of this symmetric system of equations. A space can be formed from the eigenvectors obtained and the reconstructed source is written as an expansion of the eigenvectors corresponding to non-zero eigenvalues. The coefficients of the expansion are found to obtain the reconstructed BL source distribution. The problem is solved iteratively by using a permissible source region that is shrunk by removing nodes with low probability to contribute to the source. Throughout this process the permissible region shrinks from the entire object to just a few nodes. The best estimate of the reconstructed source is chosen that which minimizes the difference between the calculated and measured light fluence rates. 3D simulations presented here show that the reconstructed source is in good agreement with the actual source in terms of locations, magnitudes, sizes, and total powers for both localized multiple sources and large inhomogeneous source distributions.
Highlights
Bioluminescence tomography (BLT) has recently attracted much interest as a molecular imaging technique that can be used as a noninvasive visualization of disease progression and response to treatment [1,2,3,4]
It is assumed that different regions are distinguished through CT x-ray imaging or other imaging technique and all elements within each region have the same optical properties
A three dimensional bioluminescence reconstruction algorithm based on the diffusion equation has been developed
Summary
Bioluminescence tomography (BLT) has recently attracted much interest as a molecular imaging technique that can be used as a noninvasive visualization of disease progression and response to treatment [1,2,3,4]. To correctly estimate the distribution of luciferase in tissue from the bioluminescence image, BLT algorithms need to be developed. The number of data points available from the measured light fluence rate at the tissue boundary is much smaller than the number of unknown possible source distributions and this precludes a unique solution to the problem [6]. The calculated Green’s function and the measured light fluence rate are normalized to equate the contributions from all wavelengths. A symmetric system of equations is formed using the Green’s function and the measured light fluence rate. The removal of the eigenvectors corresponding to zero eigenvalues reduces the ill-posedness of the problem and imposes constraint on the reconstruction of the source without a need of using a regularization penalty term which reduces the accuracy of the reconstruction. The best estimate of the source is that which minimizes the difference between the calculated and measured light fluence rates
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have