Research Article  Open Access
Huangjian Yi, Duofang Chen, Wei Li, Shuang Zhou, Miao Ning, Shouping Zhu, Jie Tian, Jimin Liang, "Normalized Born ApproximationBased TwoStage Reconstruction Algorithm for Quantitative Fluorescence Molecular Tomography", Journal of Electrical and Computer Engineering, vol. 2012, Article ID 838967, 9 pages, 2012. https://doi.org/10.1155/2012/838967
Normalized Born ApproximationBased TwoStage Reconstruction Algorithm for Quantitative Fluorescence Molecular Tomography
Abstract
Fluorescence molecular tomography (FMT) is a promising technique for in vivo small animal imaging. In this paper, a twostage reconstruction method based on normalized Born approximation is developed for FMT, which includes two steps for quantitative reconstruction. First, the localization of fluorescent fluorophore is determined by norm regularization method. Then, in the location region of fluorophore, which is provided by the first stage, algebraic reconstruction technique (ART) is utilized for the fluorophore concentration reconstruction. The validity of the twostage quantitative reconstruction algorithm is testified by simulation experiments on a 3D digital mouse atlas and physical experiments on a phantom. The results suggest that we are able to recover the fluorophore location and concentration.
1. Introduction
Fluorescence molecular tomography (FMT) is an emerging imaging modality that allows noninvasive, quantitative, threedimensional (3D) imaging of fluorescence probes associated with molecular and cellular functions [1–3]. It can be applied to drug discovery and therapeutic response in small animals [4, 5]. In this imaging modality, an external excitation light illuminates the tissue surface of the small animal, which is labeled with fluorescence probes beforehand. Then, the fluorescence probes are excited to emit photons, which are detected as fluorescence on the surface. With the measured data on the surface, the spatial distribution and concentration of fluorophore can be reconstructed based on certain inverse mathematical models [6, 7].
In recent years, great efforts have been made in developing new reconstruction algorithms [7–10] and imaging system [11–13]. However, challenges remain in quantitative fluorescence tomography. In [7], Song et al. developed a novel adaptive finite element algorithm based on the longest refinement method for freespace fluorescent tomography. This algorithm recovered the spatial distribution of the fluorescence target, but the reconstructed fluorescent yield only reached a half of the true fluorescence yield. In order to reduce the influence of optical properties, the normalized Born ratio has been developed for FMT [8]. To employ a priori information of the sparseness of fluorescence target, Han et al. proposed a fast reconstruction algorithm based on sparsity regularization for FMT [9]. The fluorescence source can be localized, but the fluorescence intensity was recovered with nearly 38% error. An integrated trimodality system that combines fluorescence, diffuse optical, and Xray tomography (FT/DOT/XCT) was developed for quantitative fluorescence tomography [10]. The diffuse optical tomography (DOT) provided the functional optical background and Xray tomography (XCT) provided the structural information, which included the location of the fluorophore [10, 13]. The fluorophore could be localized accurately and the fluorophore concentration was recovered within 8% error only when both functional and structural a priori information was utilized for fluorescence tomography (FT) reconstruction algorithm. This paper inspired us to propose the twostage method.
In this paper, we proposed a twostage reconstruction algorithm based on normalized Born approximation, which can largely reduce the effects of incorrect assumption of optical properties [8]. First, we localized the fluorophore by norm regularization method to fully take advantage of the sparsity of fluorophore. Then, in the region where the fluorophore is located, algebraic reconstruction technique (ART) is utilized to recover the fluorophore concentration. The first stage not only localizes the fluorophore but also provides the reconstruction region for the second stage. So, the first stage is very important for the whole quantitative reconstruction. The developed method was validated by both simulation experiments on a 3D digital mouse and physical experiments on a cubic phantom. The results have shown the feasibility and potential of the algorithm for quantitative fluorescence tomography.
2. Method
2.1. Forward and Inverse Problems
In the nearinfrared (NIR) bandwidth (700–900 nm), the light propagation in the biological tissue can be modeled using diffusion equation (DE) with Robintype boundary condition [14, 15]. The Green’s function due to a continuous wave (CW) point excitation source describes the light propagation. It can be obtained as follows: where is the diffusion coefficient at position . is the reconstruction domain with a boundary of . is the absorption coefficient and is a constant depending on the optical reflective index mismatch at the boundary. represents the outward normal vector to the boundary. The diffusion equation can be solved using finite element method (FEM) [15, 16]. In order to reduce the influence of tissue heterogeneity, normalized Born approximation is used [8]. The normalized ratio of the measured fluorescence signal and corresponding excitation signal at detector is given as follows [8, 17]: where is a calibration factor which depends on the genetic system characteristics. is the point inside the volume . denotes the fluorescence yield to be reconstructed, which is directly related to the fluorophore concentration by the formula , where is the molar extinction coefficient, is the quantum efficiency, and is the concentration of the fluorophore [18]. describes the light propagation at the excitation wavelength due to a point source at , and corresponds to the light propagation at the emission wavelength due to a point source at . By discretizing the volume into voxels, a linear system can be generated for multiple sourcedetector pairs: where is weight matrix and each element represents the contribution of th voxel of the discrete field of interest to the measurement due to the th sourcedetector pair. The unknown fluorescence yield in each voxel is recovered by solving (3) using certain regularization method. Detailed descriptions can be found in [17].
2.2. TwoStage Reconstruction Algorithm
Due to the illposed and the illconditioned nature of FMT inverse problem, many kinds of regularization methods are employed to obtain a stable solution. Here, we developed a twostage reconstruction algorithm for quantitative fluorescence tomography. Considering the fact that the volume of the fluorophore is very small compared with the whole reconstruction region, the fluorophore can be seen as sparsely distributed in the region of the interested. First, norm regularization is utilized to fully take advantage of the sparseness of the fluorophore. As a result, (3) can be solved by the following minimization problem: where and is the regularization parameter. From the results in the first stage, we determined the permissible fluorophore region for the second stage by selecting the nodes with 30% of the maximum reconstructed value. Since the reconstruction region is small, the norm regularization method is not suitable for the reconstruction. Then, algebraic reconstruction technique (ART) is utilized to recover the fluorescence yield. Its form is as follows [19, 20]: where is the relaxation parameter. and are the “”th and ()th estimation of the th element of the unknown variables . The initial value of the ART is provided by the first stage.
Algorithm. The flow of the proposed method is depicted in Figure 1, which mainly includes.
Stage 1. Reconstruct the distribution of fluorophore on the whole body of the mouse by minimizing (4), that is, localize the fluorophore by norm regularization method.
Stage 2. Recover the fluorescence yield using ART on the permissible fluorophore region which is provided by the first stage, that is, quantify the fluorescence concentration.
3. Experiments and Results
In this section, two groups of verification experiments were conducted to validate the potential and feasibility of the proposed twostage reconstruction algorithm for quantitative fluorescence tomography. Reconstructions were carried out on a personal computer with a 2.66 GHz Intel(R) Xeon(R) CPU E5430 and 8.00 GB RAM. To analyze the results quantitatively, location error (LE), which is the Euclidean distance between the centers of the reconstructed and the actual fluorophore, is defined by , where is the reconstructed coordinate with the maximum reconstructed value and is the actual coordinate. The relative fluorescence yield error is defined to be , where and are the recovered maximum fluorescence yield and the actual one, respectively.
3.1. Digital Phantom Experiment
3.1.1. Setup for Simulation Studies
In the simulation experiments, a digital mouse atlas of CT and cryosection data was utilized to provide anatomical information [21]. We only considered the torso section of the mouse with a height of 35 mm. A cylindrical fluorescent target with a 0.8 mm radius and 1.6 mm height was placed in the liver with center at 11.9 mm, 6.4 mm, and 16.4 mm as shown in Figure 2(a), and the fluorescence yield of the fluorescent target was set to be 0.05. The optical parameters for different organs are shown in Table 1 [6, 22]. The excitation and emission data were synthesized using finite element method (FEM) based on (1). In this experiment, we considered 36 point sources at different positions in sequence as shown in Figure 2(b). For each point source, 39 detectors were used for collecting fluorescence signal inside 120° field of view (FOV), as shown in Figure 2(c). So there are total 1404 sourcedetector pairs.

(a)
(b)
(c)
3.1.2. Reconstruction of the Single Fluorophore
The reconstruction region was discretized into 4157 nodes and 21811 tetrahedral elements. In the first stage, we employed the incomplete variables truncated conjugate gradient method (IVTCG) to solve (4) [23]. Since determining a good regularization parameter is a very challenging task, we did not represent this process here in detail. Instead, we selected the range between and which is sufficient for our experiments. In the second stage, the ART iterations were terminated when the iteration number reached the maximum ART iteration number 20 or the relative change between neighboring iteration was less than 0.5%. To clearly show the performance of the proposed twostage reconstruction for FMT, the recovered results are showed in Figure 3. And Table 2 gives the LE and relative error of fluorescence yield at each stage. It is clear that the fluorescence yield on the selected slice for the second stage is brighter than that of the first one in Figure 3.

(a)
(b)
(c)
(d)
We also used algebraic reconstruction technique (ART) and norm regularization method to solve this problem, respectively, as a comparison. Figure 4 shows the reconstructed results. It is clear that even though both the ART and norm regularization methods can provide a good location of the fluorescent target, yet they both fail to recover the fluorescence yields accurately. Compared with our proposed method, they recovered the fluorescence yields with a large error.
(a)
(b)
3.2. Physical Phantom Experiment
3.2.1. Setup for Experiment Studies
A full angle noncontact FMTimaging system is employed for collecting fluorescence signals [6]. The schematic of this system is shown in Figure 5 and detailed descriptions can be found in [6]. A 672 nm continuous wave (CW) laser source illuminated the target, which was placed at the rotational stage. A highly sensitive chargecoupled device (CCD) camera (Princeton Instruments PIXIS 2048B, Roper scientific, Trenton, NJ), which was cooled to −70°C, was utilized for recording the transmitted light at both the excitation and emission wavelengths using appropriate filters. Here, custommade 40 nm bandpass filter (HZXD, Beijing, China) centered at 660 nm and 35 nm bandpass filter (HZXD, Beijing, China) centered at 720 nm were used to allow light transmission at the excitation and emission wavelength, respectively.
A cubic phantom with a side length of 20 mm was placed on the rotational stage, shown in Figures 6(a)–6(c). The phantom was made from polyoxymethylene. The optical parameters for both excitation and emission wavelengths are mm^{−1}, mm^{−1}, mm^{−1}, and mm^{−1} respectively, which were determined by diffuse optical tomography. There is a small hole with a 1mm radius and 2 mm height, which was used to emplace the 3804.4 nM Cy5.5 solution (with the quantum efficiency of 0.23 at the peak excitation wavelength of 671 nm and the extinction coefficient of about 0.019 mm^{−1} [24]). The fluorophore was excited by 8 point sources at mm plane, and there are 160 detectors being located on the boundary inside 90° FOV corresponding to each point source, as shown in Figures 6(d) and 6(e).
(a)
(b)
(c)
(d)
(e)
3.2.2. Reconstruction for Experimental Data
The cubic in Figure 7 was discretized into 3191 nodes and 15870 tetrahedral elements for the inverse problem. The reconstruction process was the same as in Section 3.1.2. Figure 7 shows the recovered results by our method. The reconstructed center coordinate of Cy5.5 is 14.0 mm, 9.1 mm, and 7.7 mm while the real center coordinate is 16.0 mm, 8.0 mm, and 8.5 mm with location error of 2.42 mm. The recovered concentration of Cy5.5 dye is about 3478.3 nM with a relative error of 8.6%. Figure 8 shows the corresponding results by algebraic reconstruction technique (ART) and norm regularization method. It is shown that their fluorescence yields were smaller than the actual one.
(a)
(b)
(c)
(a)
(b)
4. Discussion and Conclusion
The quantitative reconstruction of fluorescence molecular tomography (FMT) is meaningful for biological applications. In this paper, we developed a twostage reconstruction method based on normalized Born approximation for quantitative fluorescence tomography. This algorithm includes two steps. The first step is to localize the permissible region of fluorophore. In this stage, the fluorescent target is very small compared with the whole reconstruction region, which can be considered as a priori information in the form of norm penalty. In the second step, algebraic reconstruction technique (ART) is used to recover the fluorescence yield in the permissible region of fluorophore, which was provided by the first step.
Both the numerical experiment and physical experiment were conducted to validate the feasibility and potential of the proposed algorithm. In all experiments, our method was compared with ART and norm regularization method, and the developed algorithm showed a better performance. In the physical phantom experiment, the differences in the positions of reconstructed fluorescent target and actual one for all methods were obvious. There are many possible reasons for this phenomenon. First, the strong tissue scattering makes the optical imaging in thick tissue with lowresolution. Second, the illposed and illconditioned nature of the inverse problem makes the solution unstable and nonunique. Therefore, the recovered results are easily affected by the measurement noise, which is inevitable during the signal detection. Last but not least, the system errors caused by geometry mismatch and optical modeling approximation might also lead to large errors in the reconstruction algorithm [25].
Generally speaking, compared with conventional methods, our method is easy to use with high accuracy. The outstanding advantage of this method is that it only utilizes the optical information to provide the permissible region of target. In [10, 13], the authors used Xray computed tomography (XCT) to provide an accurate exterior boundary of object and anatomical images of the animal as a priori information, on which the fluorescent target can be localized. In some clinical problems, however, the targets cannot be obtained by XCT, like tumors or anomaly tissues. In this case, our method showed its advantages, as we employ the norm regularization to localize the object instead of XCT. The weakness of the proposed approach is that it would not perform well when localization of object cannot be obtained accurately in the first stage. If the target is localized with a large error in the first stage, the conduction of the second stage would become of no meaning.
In conclusion, we have proposed an effective method for quantitative fluorescence tomography. Numerical experiments and physical experiments were both conducted to validate the feasibility and potential of our method. This algorithm would be applied in in vivo small animal tumor experiments in the future.
Acknowledgments
This work was supported by the Program of the National Basic Research and Development Program of China (973) under Grant no. 2011CB707702, the National Natural Science Foundation of China under Grant nos. 81090272, 81101083, 81101084, 81101100, 81000632, and 30900334, the National Key Technology Support Program under Grant no. 2012BAI23B06, and the Fundamental Research Funds for the Central Universities.
References
 V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of wholebody photonic imaging,” Nature Biotechnology, vol. 23, no. 3, pp. 313–320, 2005. View at: Publisher Site  Google Scholar
 J. Tian, J. Bai, X. P. Yan et al., “Multimodality molecular imaging: improving image quality,” IEEE Engineering in Medicine and Biology Magazine, vol. 27, no. 5, pp. 48–57, 2008. View at: Publisher Site  Google Scholar
 V. Ntziachristos, “Fluorescence molecular imaging,” Annual Review of Biomedical Engineering, vol. 8, pp. 1–33, 2006. View at: Publisher Site  Google Scholar
 J. K. Willmann, N. van Bruggen, L. M. Dinkelborg, and S. S. Gambhir, “Molecular imaging in drug development,” Nature Reviews Drug Discovery, vol. 7, no. 7, pp. 591–607, 2008. View at: Publisher Site  Google Scholar
 X. Montet, V. Ntziachristos, J. Grimm, and R. Weissleder, “Tomographic fluorescence mapping of tumor targets,” Cancer Research, vol. 65, no. 14, pp. 6330–6336, 2005. View at: Publisher Site  Google Scholar
 H. Yi, D. Chen, X. Qu et al., “Multilevel, hybrid regularization method for reconstruction of fluorescent molecular tomography,” Applied Optics, vol. 51, no. 7, pp. 975–986, 2012. View at: Publisher Site  Google Scholar
 X. Song, D. Wang, N. Chen, J. Bai, and H. Wang, “Reconstruction for freespace fluorescence tomography using a novel hybrid adaptive finite element algorithm,” Optics Express, vol. 15, no. 26, pp. 18300–18317, 2007. View at: Publisher Site  Google Scholar
 A. Soubret, J. Ripoll, and V. Ntziachristos, “Accuracy of fluorescent tomography in the presence of heterogeneities: study of thenormalized born ratio,” IEEE Transactions on Medical Imaging, vol. 24, no. 10, pp. 1377–1386, 2005. View at: Publisher Site  Google Scholar
 D. Han, J. Tian, S. Zhu et al., “A fast reconstruction algorithm for fluorescence molecular tomography with sparsity regularization,” Optics Express, vol. 18, no. 8, pp. 8630–8646, 2010. View at: Publisher Site  Google Scholar
 Y. Lin, W. C. Barber, J. S. Iwanczyk, W. Roeck, O. Nalcioglu, and G. Gulsen, “Quantitative fluorescence tomography using a combined trimodality FT/DOT/XCT system,” Optics Express, vol. 18, no. 8, pp. 7835–7850, 2010. View at: Google Scholar
 R. B. Schulz, A. Ale, A. Sarantopoulos et al., “Hybrid system for simultaneous fluorescence and Xray computed tomography,” IEEE Transactions on Medical Imaging, vol. 29, no. 2, pp. 465–473, 2010. View at: Publisher Site  Google Scholar
 F. Liu, X. Liu, D. Wang, B. Zhang, and J. Bai, “A Parallel excitation based fluorescence molecular tomography system for wholebody simultaneous imaging of small animals,” Annals of Biomedical Engineering, vol. 38, no. 11, pp. 3440–3448, 2010. View at: Publisher Site  Google Scholar
 Y. Lin, W. C. Barber, J. S. Iwanczyk, W. W. Roeck, O. Nalcioglu, and G. Gulsen, “Quantitative fluorescence tomography using a trimodality system: in vivo validation,” Journal of Biomedical Optics, vol. 15, no. 4, Article ID 040503, 2010. View at: Google Scholar
 S. R. Arridge, “Optical tomography in medical imaging,” Inverse Problems, vol. 15, no. 2, pp. R41–R49, 1999. View at: Publisher Site  Google Scholar
 M. Schweiger, S. R. Arridge, M. Hiraoka, and D. T. Delpy, “The finite element method for the propagation of light in scattering media: boundary and source conditions,” Medical Physics, vol. 22, no. 11, pp. 1779–1792, 1995. View at: Publisher Site  Google Scholar
 S. R. Arridge, M. Schweiger, M. Hiraoka, and D. T. Delpy, “A finite element approach for modeling photon transport in tissue,” Medical Physics, vol. 20, no. 2, pp. 299–310, 1993. View at: Publisher Site  Google Scholar
 X. Liu, D. Wang, F. Liu, and J. Bai, “Principal component analysis of dynamic fluorescence diffuse optical tomography images,” Optics Express, vol. 18, no. 6, pp. 6300–6314, 2010. View at: Publisher Site  Google Scholar
 M. A. Naser and M. S. Patterson, “Improved bioluminescence and fluorescence reconstruction algorithms using diffuse optical tomography, normalized data, and optimized selection of the permissible source region,” Biomedical Optics Express, vol. 2, no. 1, pp. 169–184, 2011. View at: Publisher Site  Google Scholar
 A. Kak and M. Slaney, Computerized Tomographic Imaging, IEEE Press, New York, NY, USA, 1987.
 D. Wang, X. Liu, Y. Chen, and J. Bai, “A novel finiteelementbased algorithm for fluorescence molecular tomography of heterogeneous media,” IEEE Transactions on Information Technology in Biomedicine, vol. 13, no. 5, pp. 766–773, 2009. View at: Publisher Site  Google Scholar
 B. Dogdas, D. Stout, A. F. Chatziioannou, and R. M. Leahy, “Digimouse: a 3D whole body mouse atlas from CT and cryosection data,” Physics in Medicine and Biology, vol. 52, no. 3, pp. 577–587, 2007. View at: Publisher Site  Google Scholar
 G. Alexandrakis, F. R. Rannou, and A. F. Chatziioannou, “Tomographic bioluminescence imaging by use of a combined opticalPET (OPET) system: a computer simulation feasibility study,” Physics in Medicine and Biology, vol. 50, no. 17, pp. 4225–4241, 2005. View at: Publisher Site  Google Scholar
 X. He, J. Liang, X. Wang et al., “Sparse reconstruction for quantitative bioluminescence tomography based on the incomplete variables truncated conjugate gradient method,” Optics Express, vol. 18, no. 24, pp. 24825–24841, 2010. View at: Publisher Site  Google Scholar
 F. Gao, H. Zhao, Y. Tanikawa, and Y. Yamada, “A linear, featureddata scheme for image reconstruction in timedomain fluorescence molecular tomography,” Optics Express, vol. 14, no. 16, pp. 7109–7124, 2006. View at: Publisher Site  Google Scholar
 X. He, J. Liang, X. Qu, H. Huang, Y. Hou, and J. Tian, “Truncated total least squares method with a practical truncation parameter choice scheme for bioluminescence tomography inverse problem,” International Journal of Biomedical Imaging, vol. 2010, Article ID 291874, 11 pages, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 Huangjian Yi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.