|Home | About | Journals | Submit | Contact Us | Français|
We compare frequency-and time-domain formulations of deep-tissue fluorescence imaging of turbid media. Simulations are used to show that time-domain fluorescence tomography, implemented via the asymptotic lifetime-based approach, offers a significantly better separability of multiple lifetime targets than a frequency-domain approach. We also demonstrate experimentally, using complex-shaped phantoms, the advantages of the asymptotic time-domain approach over a Fourier-based approach for analyzing time-domain fluorescence data.
Optical technologies for noninvasive macroscopic fluorescence imaging in biological media utilize three main approaches: time domain (TD) using pulsed light sources [1–5], frequency domain (FD) using megahertz modulated sources [6–8], and continuous wave (CW) using steady state light sources [9–11]. Of these, the TD approach is the most comprehensive, since a short laser pulse (fs-ps) implicitly contains all the modulation frequencies, including the zero-frequency component. The tomographic analysis of TD data can, however, pose computational challenges. Several simplifications have been attempted, primarily using derived data types such as the fast Fourier transform (FFT) [5,7]. The FFT simplifies the TD forward problem but reduces it to a FD forward problem, identical to that for a genuine FD measurement. Alternatively, TD fluorescence data may also be analyzed by estimating lifetimes directly from the asymptotic region, followed by the separate inversion of the yield of each lifetime component [3,12]. The question arises as to how the FD forward problem compares with the asymptotic TD (ATD) approach, when lifetime sensitive targets are used. In this Letter, we address this question both with simulated and experimental TD data in the context of small animal imaging applications.
Consider a diffuse imaging medium of finite support Ω, embedded with fluorophores characterized by yield and lifetime distributions, η (r) and τ(r). In the FD approach, the forward problem takes the form 
Here, Gx and Gm are the FD Green’s functions at frequency ω for propagation from a source rs and a detector rd, respectively, to a medium point r. Equation (1) is first inverted to obtain F(r,ω)=τ(r)η(r)/[1–iωτ(r)]. The lifetimes are then obtained from the phase Ø(r,ω) of F, as τ(r)=ø(r,ω)/ω. Finally, the yield reconstructions are given by . In the ATD approach, the decay amplitudes, an(rs,rd), for each lifetime component, τn =1/Γn, take the place of the Fourier amplitude, , as the measurement data set . The an’s are related to the yield distribution, ηn(r), for the lifetime component at τn, through a linear forward problem:
Note that in Eq. (2), the Green’s functions are the same as that for the FD case in Eq. (1), but evaluated at an imaginary frequency of –iΓn. (see  for details). The difference between Eq. (1) and Eq. (2) is clear if we express F(r,ω) in terms of discrete lifetimes rather than the continuous distributions τ(r). We then get F(r,ω)=∑nτnηn(r)/(1–iωτn). Thus the FD approach recovers a single distribution, F(r,ω), that is a mixture of the yields for all the lifetimes present in the imaging medium. Separability of multiple lifetime components using the FD approach will therefore be limited by the intrinsic resolution of the system. The ATD approach, however, results in a separate forward problem for each lifetime, so that only the widths of the individual ηn(r)’s are subject to the resolution limitations. The spatial separability of the ηn’s for different lifetimes in the ATD approach relies rather on the ability to recover the lifetimes through multiexponential fits of TD decays.
To compare the performance of the FD and ATD approaches, simulations were performed for a 2 cm thick infinite slab geometry, assuming 2 mm × 2mm × 2 mm size fluorescent targets (also the size of a single voxel) placed both laterally (Fig. 1) and axially (Fig. 2) with respect to the source-detector line of sight. The inclusions had distinct lifetimes of 0.5 and 1 ns. 29 sources and 21 detectors were arranged [Fig. 1(a)] along the planes at depths of z=0 cm and z=2 cm, respectively. The background medium was assumed to have values of for scattering and μa=0.1 cm−1 for absorption (at both excitation and emission wavelengths). The FD signal was simulated using Eq. (1) in the diffusion approximation, for ω =80 MHz, and reconstructed with Tikhonov regularization. The full TD signal was also simulated in the diffusion model with shot noise added, and the an(rs,rd)’s were recovered asymptotically. Figures 1(b) and 1(c) show the FD recovered η(r) and τ(r), and Fig. 1(d) shows the ATD recovered η(r) for 0.5 and 1 ns, for a lateral target separation of 6 mm. Figures 1(f)–1(h) show the corresponding images for 3 mm separation. The x dependence of the yield for each case is plotted at the (y,z) location of the maximum of the corresponding distributions in Figs. 1(e) and 1(i). It is clear that while the ATD approach localizes and resolves the targets for both separations, the FD approach does not clearly resolve the targets for the 3 mm case [Figs. 1(f) and 1(i)] and results in an inaccurate lifetime estimate [Fig. 1(g)]. Also, τ(r) exhibits a continuous variation in the region between the targets [Fig. 1(g)]. Figure 2 shows the results when targets are located axially, i.e., along the normal to the source and detector planes. The ATD approach resolves the targets for both an 8 and a 4 mm separation as seen from Figs. 2(d) and 2(h). For the 8 mm separation, the FD approach barely resolves the targets, but τ(r) is obtained accurately near the target locations. However, τ(r) again varies continuously between 0.5 and 1 ns in the region between the targets [Fig. 2(b)]. Thus a continuous distribution of lifetimes is recovered by the FD approach for both axial and lateral targets, although discrete lifetimes were simulated. The FD approach completely fails to recover the targets for the 4 mm separation [Fig. 2(h), black curve] and recovers an inaccurate lifetime distribution [Fig. 2(f)]. We also independently verified that using multiple frequencies does not significantly improve the quality of the FD reconstructions. Specifically, the 3 mm separated lateral targets in Fig. Fig.11 and 4 mm axial targets in Fig 2 were not resolved with the use of up to 5 frequencies, even using a simplified linear approach where the lifetimes were assumed to be known in advance. (Multifrequency FD is a nonlinear problem .)
The simulations presented in Figs Figs11 and and22 clearly indicate the advantages of the ATD approach over the FD formalism for resolving lifetime contrast. We now experimentally demonstrate these results using a noncontact TD fluorescence tomography system, consisting of a Ti:Sapphire laser source for excitation and a time-gated intensified CCD camera for detection. Detailed aspects of the system are discussed in a separate publication . A mouse-shaped phantom made of epoxy, ink, and TiO2 combination (, μa=0.1 cm−1) was used as the imaging subject [Fig. 3(a)]. Two inclusions (≈3 mm × 3mm × 3 mm in size, 6 mm center-to-center separation) within the phantom were filled with a near-infrared dye (3,3′-diethylthiatricarbocyanine, absorption and emission maxima near 755 and 790 nm), mixed in distilled water and 100% glycerol solutions, giving lifetimes of 0.5 and 0.95 ns. The data were collected in the transillumination geometry for 21 source positions (1 mm × 1 mm grid) and 33 assigned detectors (1 mm × 2 mm grid) on the CCD camera image. The lifetimes were directly recovered from a global analysis fit  to the experimental decay curves as near 0.5 and 0.95 ns, and the decay amplitudes at all the S-D pairs were used in Eq. (2) to recover the yield. The 3D reconstructions of the yield of the two lifetimes recovered using the ATD approach are shown in Fig. 3(b), overlayed with the surface boundary of the phantom, which was acquired using photogram-metric methods . The first Fourier component of the data, at 156 Mhz, was also used in a FD approach using Eq. (1). Reconstructions were performed for data with both one lifetime (single inclusion filled) and two lifetimes (both inclusions filled). It is clear from the tomography results of Figs. 3(c)–3(j) that the ATD approach recovers the target depths for both the single and two lifetime cases. While the FD approach recovers the correct localization and lifetime for the single lifetime case, it does not localize the targets separately when both inclusions are simultaneously filled and estimates an inaccurate lifetime of ≈0.7 ns. It should be noted again that, as seen in the simulations, the ATD approach is able to separate targets located well within the full width at half-maximum of the yield for a single inclusion, which was >1 cm [Fig. 3(f)].
We have shown using simulations and experiments that an asymptotic multiexponential-analysis-based TD approach to fluorescence tomography has inherent advantages over the FD formalism of diffuse fluorescence when resolving discrete lifetimes present within a turbid medium. Fluorescence lifetime has been widely exploited in microscopy techniques with thin tissue sections (see, for example, ), and it is likely to serve as a key functional parameter for non-invasive diagnostic optical molecular imaging and drug discovery. The findings reported in this Letter have significant relevance for the design of lifetime-based optical molecular imaging systems. The results also strongly motivate further development of smart near-infrared molecular probes that specifically shift lifetime upon binding to disease targets of interest.
This work was supported by the National Institutes of Health (grants EB000768 and AG026240).
OCIS codes: 170.3880, 170.3650.