Research Papers: Imaging

Improvement of fluorescence-enhanced optical tomography with improved optical filtering and accurate model-based reconstruction algorithms

[+] Author Affiliations
Yujie Lu, Banghe Zhu, Chinmay Darne, I-Chih Tan, John C. Rasmussen, Eva M. Sevick-Muraca

University of Texas Health Science Center at Houston, Center for Molecular Imaging, Institute of Molecular Medicine, 1825 Pressler Street, SRB 330A, Houston, Texas 77030

J. Biomed. Opt. 16(12), 126002 (December 14, 2011). doi:10.1117/1.3659291
History: Received September 06, 2011; Accepted October 18, 2011; Revised October 17, 2011; Published December 14, 2011; Online December 14, 2011
Text Size: A A A

Open Access Open Access

* Address all correspondence to: Eva M. Sevick-Muraca, University of Texas Health Science Center at Houston, Center for Molecular Imaging, Institute of Molecular Medicine, 1825 Pressler Street, SRB 330A, Houston, Texas 77030; Tel: 713-500-3560; Fax: 713-500-0319; E-mail: eva.sevick@uth.tmc.edu

The goal of preclinical fluorescence-enhanced optical tomography (FEOT) is to provide three-dimensional fluorophore distribution for a myriad of drug and disease discovery studies in small animals. Effective measurements, as well as fast and robust image reconstruction, are necessary for extensive applications. Compared to bioluminescence tomography (BLT), FEOT may result in improved image quality through higher detected photon count rates. However, background signals that arise from excitation illumination affect the reconstruction quality, especially when tissue fluorophore concentration is low and/or fluorescent target is located deeply in tissues. We show that near-infrared fluorescence (NIRF) imaging with an optimized filter configuration significantly reduces the background noise. Model-based reconstruction with a high-order approximation to the radiative transfer equation further improves the reconstruction quality compared to the diffusion approximation. Improvements in FEOT are demonstrated experimentally using a mouse-shaped phantom with targets of pico- and subpico-mole NIR fluorescent dye.

Figures in this Article

Fluorescence-enhanced optical tomography (FEOT) can be confounded by effects of autofluorescence and a high “noise floor,” which arises from excitation light leakage through optical rejection filters. This high noise floor can obscure signals from low-concentration fluorophores in tissues and impact reconstruction quality. Because tissue autofluorescence peaks at visible wavelengths [typically, for example, the peak of mouse skin is 500 nm (Ref. 1)] and exponentially reduces with increase of wavelength, excitation in near-infrared (NIR) FEOT (excitation wavelengths ⩾750 nm) effectively removes the confounding artifact of tissue autofluorescence.2 Nonetheless, NIR FEOT is not immune to high noise floors owing to excitation light leakage. The performance deterioration of optical filters, including reduced and blueshifted optical densities of the interference filters, occurs when scattered excitation light is incident at non-normal directions.3 In addition to accurate measurements, model-based reconstruction methods require precise mathematical models to describe photon propagation and generation in tissues. Although the radiative transfer equation (RTE) is the choice of method, it is complicated and can impose severe time constrains for directly obtaining solutions in complex geometries as required for a rodent. Although the diffusion approximation (DA) has been extensively applied in optical tomography at its early stages in development, it becomes increasingly inaccurate in small volumes (such as a mouse) and under conditions of high absorption (such as in the rodent liver).

In this paper, we demonstrate that (i) experimental measurements optimized by simple filter configurations to reduce background signals and (ii) accurate models of light propagation together improve NIR FEOT. The measurement sensitivity and overall quality of NIR FEOT begins with the detector. Over the past decade, our group developed an intensified charge-coupled device (ICCD) camera system to realize frequency-domain, time-dependent, and continuous wave noncontact fluorescence measurements. In these experiments, one 785-nm notch filter and one 830-nm bandpass filter [Fig. 1] were used to reject excitation light leakage while allowing collection of emission photons from human subjects under noninvasive imaging conditions following microdosing (<100 and >10 μg of ICG administrated to humans).4 Herein, we employed an optimized filter configuration, that is two 830-nm bandpass filters located before and after a 28-mm NIKKOR focusing lens [Fig. 1] to further reduce background signal owing to excitation light leakage.

Graphic Jump LocationF1 :

(a) Unoptimized and (b) optimized filter configurations in the gain-modulated intensified detection system.

Using both configurations, phantom FEOT studies were conducted using a mouse-shaped solid phantom (Caliper Life Sciences, Hopkinton, Massachusetts). The scattering and absorption coefficients, anisotropic factor, and refractive factor of the phantom were taken to be 9.5 mm−1, 0.0066 mm−1, 0.9, 1.5 at 785 nm and 7.4 mm−1, 0.0077 mm−1, 0.9, 1.5 at 830 nm, respectively as provided by the manufacturer. ICG of 2 μmol/l, 0.5 μmol/l, and 0.125 μmol/l was sealed in plastic volumes of 2.35 mm3, and total molar quantities range from 4.7 to 0.29 pmol. The absorption coefficient of ICG was measured (Spectrophotometer DU800, Beckman Coulter, Brea, California), and for fluorescent targets in the phantom, we computed the absorption ratio of fluorophore to surrounding tissue (AR) to range 29.54–0.51. The ICG targets were then placed into a predrilled hole within the mouse phantom, and a rod comprised of the same material as the phantom was used to fill the remaining volume of the hole. Excitation light of 8.1 mW first illuminated an area of 2.0 × 2.0 mm2 on the dorsal surface while transmission fluorescence measurements were taken on the ventral view and, secondly, on the ventral surface while transmission fluorescence measurements were taken on the dorsal surface. The exposure times were fixed at 800 ms, and the gain of the intensifier (the voltage relevant to the gain ranges from 6.51 to 8.77 V) was adjusted to reach equivalent maximum fluorescent counts for each view. It is noteworthy that the increase of the gain does not affect the detected photon distribution, although it improves the sensitivity and reduces detectable dynamic range.

Using the simplified spherical harmonics approximation (SPN),56 we previously developed a fully parallel linear reconstruction algorithm for bioluminescence tomography (BLT).7 In the algorithm developed herein, the reconstruction is significantly accelerated using parallel implementation in the cluster. The linear reconstruction framework is easily used in the combination of multispectral and multiview measurements to improve the reconstruction quality. In addition, we also extend our parallel reconstruction framework for NIR FEOT. In this reconstruction algorithm, we solve the following linear least-squares problem: Display Formula

1min0<μasf<μasf, sup Θ(μasf):AμasfJ+,m,b2.
where J+, m, b is the measurable exiting partial current for emission (W mm−2); μasf is the absorption coefficient (mm−1) of the fluorophore; μasf, sup is the upper bound constraint; and A denotes the linear relationship between J+, m, b and μasf. When multiple emission measurements obtained from multiple excitations are used in reconstruction, J+, m, b and A consist of [J1+,m,b,...,Ji+,m,b,...,JNm+,m,b]T and [A1,...,Ai,...,ANm]T, where Nm denotes the total number of emission measurements and T is a transpose operator. In order to acquire Ai, we use similar methods found in the literature7 for the SP3 emission approximation, Display Formula
2M1ϕ1mM1ϕ2mM2ϕ1mM2ϕ2mϕ1mϕ2m=Bm23Bmμasf,
where Miϕjm is the submatrix corresponding to ϕjm in the i-th SP3 equation when the finite element method is used and Bm is obtained by its components bpqmDisplay Formula
3bpqm=ΩQφxυp·υqdr,
where Ω is the domain for reconstruction; r is the location in Ω; vp, q are the shape functions; Q is the quantum efficiency of the fluorophore; and ϕx is the fluence of the excitation and is obtained by directly solving the SP3 excitation approximation when omitting the absorption coefficient of the fluorophore. After a series of matrix deductions from Eq. 2, Ai can be obtained.7 Because of the ill-posed nature of A, several factors, such as measurement noise and mathematical model errors, significantly affect the reconstruction quality. Regularization methods have become popular to reduce the effect. In this work, we need to evaluate the effect of mathematical models in the reconstruction. Therefore, regularization term is not used in Eq. 1.

Figures 2 show the measured surface emission photon distribution from the dorsal and ventral projections, respectively, when the inclusion contained 4.7 pmol of ICG and two 830 nm filters were used. The profiles extracted along the lines shown in the top panels of Fig. 2 are shown in Figs. 2 (dorsal) and 2 (ventral) for the data from different filters and ICG target concentrations. Because of the nonlinear relationship between the intensifier gain and the count number on the camera and measurement noise, it is difficult to make the maximum photon count number of different measurements absolutely consistent. The intensity profiles were normalized to their respective maximal count number for comparison. One can find that the photon distribution changes with the reduction of the target concentration of ICG. Because of the fluorophore inclusion is closer to the dorsal side than the ventral side, the changes of the measured photon distribution on the latter are more distinct because more excitation photons compared to emission photons were detected. Note that the normalized intensity profiles arising from different target ICG concentrations and measured using the optimized filter configuration are more similar as compared to the counterparts from unoptimized filter configurations, showing the effectiveness of excitation light leakage rejection.

Graphic Jump LocationF2 :

Surface emission photon distribution measured from the (a) dorsal and (b) ventral projections when 4.7 pmol of ICG comprised the target within the mouse phantom. (c) and (d), corresponding to (a) and (b) respectively, are the photon distribution comparisons between different ARs, molar quantities, and filters. The profiles were through the peak of the emission photon distribution along the surface line shown in the top panels.

MicroCT scanning was performed to obtain phantom volume and position of the fluorescent target. The volumetric mesh of the phantom was generated for the reconstruction using the Amira 5.0 software (Mercury Computer Systems, Inc., Chelmsford, Massachusetts) and was composed of ∼25, 000 discretized points. Using a similar registration method in the literature,8 the measured surface emission distribution and incident excitation light were mapped onto the surface of the volumetric mesh. Reconstructions were performed on a cluster of eight nodes (8 CPU cores of 3.0 GHz and 16 GB RAM at each node), and 3000 data points (about 1300 and 1700 for the ventral and dorsal sides, respectively) were used in reconstruction and reconstruction iteration number was set to 3000. Figures 34 show the results from the DA- and SP3-based reconstruction, respectively. The DA-based reconstruction time reduced from 108.0 to 19.0 min, when the number of the CPU cores used increased from 1 to 45. It is noteworthy that the SP3-based reconstruction failed on one single node with one CPU core due to memory insufficiency. When 45 CPU cores were used, the reconstruction time was 30.0 min, close to the time required for the DA-based reconstruction, and showed good performance of the fully parallel reconstruction framework. Because of the noise factors, there are some reconstructed artifacts, as shown in Figs. 34. However, when the maximum reconstructed values were used to localize the fluorophore target, localization errors of the SP3-based reconstruction were found to be smaller than those obtained from DA, as shown in Figs. 34. When 0.29 pmol of ICG and 785 and 830-nm filters were used, the position of the ICG target was not localized using either of DA- and SP3-based results due to very large errors. However, when we used the optimized filter combination and SP3-based reconstruction, we obtained better localization of the target containing 0.29 pmol ICG.

Graphic Jump LocationF3 :

DA-based reconstruction comparisons between different filter configurations at varying target fluorophore molar quantities and ARs. Cross sections with thick and thin boundaries are the center positions of the actual (70.0, 65.5, 5.5) and reconstructed fluorophore targets (maximal values), respectively. The volumetric mesh denotes the top 80% of the contour levels of reconstructed fluorophore distribution. “Error” and “Dist.” denote relative errors (x, y, z in column) and distance (in column) between the actual and reconstructed positions. “Arti. Yes” denotes whether there are artifacts inside the phantom. (units: mm).

Graphic Jump LocationF4 :

SP3-based reconstruction comparisons between different filters, molar quantities, and ARs. The display settings are the same as those in Fig. 3.

In FEOT, the localization accuracy of the fluorophore target is decided by the photon distribution on the tissue surface and model-based reconstruction. In noncontact collection mode with moderately large fields of view where the optical filtering components are subjected to light with large incident angles, excitation light leakage has significant effect on the sensitivity of fluorescence detection and tomographic reconstruction. Because of the blueshifted characteristics, the single 785-nm notch filter in the unoptimized detection scheme did not effectively reject excitation light leakage. When the excitation photon leakage is comparable to the fluorescence signal as in the case of the fluorophore with low molar quantities and ARs, the reduction of the detection performance becomes distinct. In the optimized detection system, the 830-nm bandpass filter after the lens effectively reduces the blueshifted effect because of the collection of the focused light. In addition, because of the interference phenomena between two adjoining filters, the performance of two bandpass filters cannot be improved with a direct sum of their optical density (OD). Some loss materials are needed between filters to reduce the multiple-path interference. The focus lens plays this role, improving the performance of the combined filters. Although another potential solution to remove excitation light leakage is to obtain excitation photon distribution from the same settings before and after the fluorophore is injected, such an approach complicates the experiments since gain settings are not known until after fluorescent agent is injected. The approach becomes impossible when gene-encoded fluorescent reporters are used. With the optimized optical filtering, improved localization is obtained from the SP3-based image reconstruction. The reconstruction time is significantly reduced from the fully parallel reconstruction framework. The improved FEOT is particularly important for targeted molecular imaging as tissue disease markers are typically present in pico- to femto-molar quantities with low target-to-tissue absorption ratios. Although the tissue absorption has an important effect in emission photon detection, our results have shown that sub-pmol molar quantities of ICG with a little more than half the absorption of the surrounding tissues can be reconstructed using the improved FEOT. The work shows the potential of FEOT in the future applications.

Acknowledgments

This work is supported by NIH Grants No. NIH R01CA135673 and No. NIH U54CA136404 , and a training fellowship from the Keck Center Computational Cancer Biology Training Program of the Gulf Coast Consortia (CPRIT Grant No. UNSPECIFIED RP101489 ).

Drakaki  E., , Borisova  E., , Makropoulou  M., , Avramov  L., , Serafetinides  A. A., , and Angelov  I., “ Laser induced autofluorescence studies of animal skin used in modeling of human cutaneous tissue spectroscopic measurements. ,” Skin Res. Technol.. 13, (4 ), 350–359  ((2007)).
Adams  K. E., , Ke  S., , Kwon  S., , Liang  F., , Fan  Z., , Lu  Y., , Hirschi  K., , Mawad  M. E., , Barry  M. A., , and Sevick-Muraca  E. M., “ Comparison of visible and near-infrared wavelength-excitable fluorescent dyes for molecular imaging of cancer. ,” J. Biomed. Opt.. 12, , 024017  ((2007)).
Zhu  B., , Rasmussen  J. C., , Lu  Y., , and Sevick-Muraca  E. M., “ Reduction of excitation light leakage to improve near-infrared fluorescence imaging for tissue surface and deep tissue imaging. ,” Med. Phys.. 37, , 5961  ((2010)).
Sevick-Muraca  E. M., , Sharma  R., , Rasmussen  J. C., , Marshall  M. V., , Wendt  J. A., , Pham  H. Q., , Bonefas  E., , Houston  J. P., , Sampath  L., , Adams  K. E., , Blanchard  D. K., , Fisher  R. E., , Chiang  S. B., , Elledge  R., , and Mawad  M. E., “ Imaging of lymph flow in breast cancer patients after microdose administration of a near-infrared fluorophore: feasibility study. ,” Radiology. 246,  (3 ), 734–741  ((2008)).
Klose  A. D., and Larsen  E. W., “ Light transport in biological tissue based on the simplified spherical harmonics equations. ,” J. Comput. Phys.. 220, (1 ), 441–470  ((2006)).
Chu  M., , Vishwanath  K., , Klose  A. D., , and Dehghani  H., “ Light transport in biological tissue using three-dimensional frequency-domain simplified spherical harmonics equations. ,” Phys. Med. Biol.. 54, , 2493−2509  ((2009)).
Lu  Y., , Machado  H. B., , Douraghy  A., , Stout  D., , Herschman  H., , and Chatziioannou  A. F., “ Experimental bioluminescence tomography with fully parallel radiative-transfer-based reconstruction framework. ,” Opt. Express. 17, , 16681–16695  ((2009)).
Lu  Y., , Machado  H. B., , Bao  Q., , Stout  D., , Herschman  H., , and Chatziioannou  A. F., “ In vivo mouse bioluminescence tomography with radionuclide-based imaging validation. ,” Mol. Imaging Biol.. 13, , 53–58  ((2011)).
© 2011 Society of Photo-Optical Instrumentation Engineers (SPIE)

Citation

Yujie Lu ; Banghe Zhu ; Chinmay Darne ; I-Chih Tan ; John C. Rasmussen, et al.
"Improvement of fluorescence-enhanced optical tomography with improved optical filtering and accurate model-based reconstruction algorithms", J. Biomed. Opt. 16(12), 126002 (December 14, 2011). ; http://dx.doi.org/10.1117/1.3659291


Figures

Graphic Jump LocationF2 :

Surface emission photon distribution measured from the (a) dorsal and (b) ventral projections when 4.7 pmol of ICG comprised the target within the mouse phantom. (c) and (d), corresponding to (a) and (b) respectively, are the photon distribution comparisons between different ARs, molar quantities, and filters. The profiles were through the peak of the emission photon distribution along the surface line shown in the top panels.

Graphic Jump LocationF4 :

SP3-based reconstruction comparisons between different filters, molar quantities, and ARs. The display settings are the same as those in Fig. 3.

Graphic Jump LocationF3 :

DA-based reconstruction comparisons between different filter configurations at varying target fluorophore molar quantities and ARs. Cross sections with thick and thin boundaries are the center positions of the actual (70.0, 65.5, 5.5) and reconstructed fluorophore targets (maximal values), respectively. The volumetric mesh denotes the top 80% of the contour levels of reconstructed fluorophore distribution. “Error” and “Dist.” denote relative errors (x, y, z in column) and distance (in column) between the actual and reconstructed positions. “Arti. Yes” denotes whether there are artifacts inside the phantom. (units: mm).

Graphic Jump LocationF1 :

(a) Unoptimized and (b) optimized filter configurations in the gain-modulated intensified detection system.

Tables

References

Drakaki  E., , Borisova  E., , Makropoulou  M., , Avramov  L., , Serafetinides  A. A., , and Angelov  I., “ Laser induced autofluorescence studies of animal skin used in modeling of human cutaneous tissue spectroscopic measurements. ,” Skin Res. Technol.. 13, (4 ), 350–359  ((2007)).
Adams  K. E., , Ke  S., , Kwon  S., , Liang  F., , Fan  Z., , Lu  Y., , Hirschi  K., , Mawad  M. E., , Barry  M. A., , and Sevick-Muraca  E. M., “ Comparison of visible and near-infrared wavelength-excitable fluorescent dyes for molecular imaging of cancer. ,” J. Biomed. Opt.. 12, , 024017  ((2007)).
Zhu  B., , Rasmussen  J. C., , Lu  Y., , and Sevick-Muraca  E. M., “ Reduction of excitation light leakage to improve near-infrared fluorescence imaging for tissue surface and deep tissue imaging. ,” Med. Phys.. 37, , 5961  ((2010)).
Sevick-Muraca  E. M., , Sharma  R., , Rasmussen  J. C., , Marshall  M. V., , Wendt  J. A., , Pham  H. Q., , Bonefas  E., , Houston  J. P., , Sampath  L., , Adams  K. E., , Blanchard  D. K., , Fisher  R. E., , Chiang  S. B., , Elledge  R., , and Mawad  M. E., “ Imaging of lymph flow in breast cancer patients after microdose administration of a near-infrared fluorophore: feasibility study. ,” Radiology. 246,  (3 ), 734–741  ((2008)).
Klose  A. D., and Larsen  E. W., “ Light transport in biological tissue based on the simplified spherical harmonics equations. ,” J. Comput. Phys.. 220, (1 ), 441–470  ((2006)).
Chu  M., , Vishwanath  K., , Klose  A. D., , and Dehghani  H., “ Light transport in biological tissue using three-dimensional frequency-domain simplified spherical harmonics equations. ,” Phys. Med. Biol.. 54, , 2493−2509  ((2009)).
Lu  Y., , Machado  H. B., , Douraghy  A., , Stout  D., , Herschman  H., , and Chatziioannou  A. F., “ Experimental bioluminescence tomography with fully parallel radiative-transfer-based reconstruction framework. ,” Opt. Express. 17, , 16681–16695  ((2009)).
Lu  Y., , Machado  H. B., , Bao  Q., , Stout  D., , Herschman  H., , and Chatziioannou  A. F., “ In vivo mouse bioluminescence tomography with radionuclide-based imaging validation. ,” Mol. Imaging Biol.. 13, , 53–58  ((2011)).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.