OSA's Digital Library

Biomedical Optics Express

Biomedical Optics Express

  • Editor: Joseph A. Izatt
  • Vol. 2, Iss. 4 — Apr. 1, 2011
  • pp: 871–886
« Show journal navigation

Monte Carlo based method for fluorescence tomographic imaging with lifetime multiplexing using time gates

Jin Chen, Vivek Venugopal, and Xavier Intes  »View Author Affiliations


Biomedical Optics Express, Vol. 2, Issue 4, pp. 871-886 (2011)
http://dx.doi.org/10.1364/BOE.2.000871


View Full Text Article

Acrobat PDF (722 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

Time-resolved fluorescence optical tomography allows 3-dimensional localization of multiple fluorophores based on lifetime contrast while providing a unique data set for improved resolution. However, to employ the full fluorescence time measurements, a light propagation model that accurately simulates weakly diffused and multiple scattered photons is required. In this article, we derive a computationally efficient Monte Carlo based method to compute time-gated fluorescence Jacobians for the simultaneous imaging of two fluorophores with lifetime contrast. The Monte Carlo based formulation is validated on a synthetic murine model simulating the uptake in the kidneys of two distinct fluorophores with lifetime contrast. Experimentally, the method is validated using capillaries filled with 2.5nmol of ICG and IRDye™800CW respectively embedded in a diffuse media mimicking the average optical properties of mice. Combining multiple time gates in one inverse problem allows the simultaneous reconstruction of multiple fluorophores with increased resolution and minimal crosstalk using the proposed formulation.

© 2011 OSA

1. Introduction

Optical methods offer the possibility to image numerous molecular targets with multiple distinct agents similarly to immunofluorescence microscopy and fluorescence cytometry. Current in vivo fluorescence multiplexing studies are mainly conducted in pre-clinical settings with Fluorescence Reflectance Imaging (FRI) [1

1. V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005). [CrossRef] [PubMed]

, 2

2. E. M. Sevick-Muraca and J. C. Rasmussen, “Molecular imaging with optics: primer and case for near-infrared fluorescence techniques in personalized medicine,” J. Biomed. Opt. 13(4), 041303 (2008). [CrossRef] [PubMed]

]. However, due to the limited information collected by FRI, the technique is unable to resolve the signal depth and hence to quantify it. Moreover, due to the predominance of scattering, planar imaging suffers from low resolution. For these reasons, FRI is mainly limited to superficial observations of subcutaneous flank xenograft models, surgically exposed organs or for intra-operative use, but is not appropriate to study advanced disease models in internal organs. For such applications, Fluorescence Molecular Tomography (FMT) techniques are required.

FMT acquires tomographic data sets to retrieve the 3D bio-distribution of the molecular probe used [3

3. A. H. Hielscher, “Optical tomographic imaging of small animals,” Curr. Opin. Biotechnol. 16(1), 79–88 (2005). [CrossRef] [PubMed]

]. Spectral multiplexing of fluorophores in optical techniques requires the collection of dense spectral data sets for efficient unmixing of independent signals [4

4. J. R. Mansfield, “Distinguished photons: a review of in vivo spectral fluorescence imaging in small animals,” Curr. Pharm. Biotechnol. 11(6), 628–638 (2010). [CrossRef] [PubMed]

]. Such data sets are acquired sequentially one wavelength at a time, leading to relatively long acquisition times, especially for whole body applications. Alternatively, fluorophore unmixing can be performed based on lifetime contrasts [5

5. R. Cubeddu, D. Comelli, C. D’Andrea, P. Taroni, and G. Valentini, “Time-resolved fluorescence imaging in biology and medicine,” J. Phys. D Appl. Phys. 35(9), R61–R76 (2002). [CrossRef]

]. By recording time dependent fluorescence signals, it is possible to distinguish and estimate the fractional contribution of multiple fluorophores with a monochromatic data set, allowing for fast acquisition protocols. Hence, lifetime based multiplexing is the most promising approach to simultaneously image multiple biomarkers for whole body applications in vivo.

Lifetime sensing is performed by employing Time Domain (TD) instruments. For optical tomography applications, image reconstructions based on TD fluorescence measurements have been primarily studied using the equivalent datatype in the frequency domain (FD), due to the simplicity of the relationship between the fluorescence lifetime and the measured phase in FD data [6

6. L. Zhang, F. Gao, H. He, and H. Zhao, “Three-dimensional scheme for time-domain fluorescence molecular tomography based on Laplace transforms with noise-robust factors,” Opt. Express 16(10), 7214–7223 (2008). [CrossRef] [PubMed]

, 7

7. V. Y. Soloviev, C. D’Andrea, G. Valentini, R. Cubeddu, and S. R. Arridge, “Combined reconstruction of fluorescent and optical parameters using time-resolved data,” Appl. Opt. 48(1), 28–36 (2009). [CrossRef] [PubMed]

]. Recently, several studies have investigated the potential of performing FMT directly using TD derived datatypes such as moments [8

8. S. Lam, F. Lesage, and X. Intes, “Time domain fluorescent diffuse optical tomography: analytical expressions,” Opt. Express 13(7), 2263–2275 (2005). [CrossRef] [PubMed]

], or time gates [9

9. A. T. N. Kumar, S. B. Raymond, G. Boverman, D. A. Boas, and B. J. Bacskai, “Time resolved fluorescence tomography of turbid media based on lifetime contrast,” Opt. Express 14(25), 12255–12270 (2006). [CrossRef] [PubMed]

]. In the case of lifetime multiplexed studies, unmixing algorithms applied to the decaying portion of the TPSFs have been employed to recover the fractional contribution of the different fluorescence components and perform FMT based on unmixed time-independent signals [10

10. A. T. N. Kumar, S. B. Raymond, A. K. Dunn, B. J. Bacskai, and D. A. Boas, “A time domain fluorescence tomography system for small animal imaging,” IEEE Trans. Med. Imaging 27(8), 1152–1163 (2008). [CrossRef] [PubMed]

]. However, such methodology suffers from the same limitation as the continuous wave (CW) technique, i.e., limited resolution and requires an unmixing algorithm where performances are not robust with low photon counts [11

11. M. Köllner and J. Wolfrum, “How many photons are necessary for fluorescence-lifetime measurements?” Chem. Phys. Lett. 200(1-2), 199–204 (1992). [CrossRef]

]. Conversely, the use of discrete time gates spanning the full TPSF as the data set for FMT should alleviate these drawbacks. It is well established that resolution of the optical reconstructions can be improved by using the rising portion of the TPSFs (termed early-gates) [12

12. F. Liu, K. M. Yoo, and R. R. Alfano, “Ultrafast laser-pulse transmission and imaging through biological tissues,” Appl. Opt. 32(4), 554–558 (1993). [CrossRef] [PubMed]

, 13

13. J. Wu, L. Perelman, R. R. Dasari, and M. S. Feld, “Fluorescence tomographic imaging in turbid media using early-arriving photons and Laplace transforms,” Proc. Natl. Acad. Sci. U.S.A. 94(16), 8783–8788 (1997). [CrossRef] [PubMed]

] and such technique has been recently applied to FMT [14–16

14. M. J. Niedre, R. H. de Kleine, E. Aikawa, D. G. Kirsch, R. Weissleder, and V. Ntziachristos, “Early photon tomography allows fluorescence detection of lung carcinomas and disease progression in mice in vivo,” Proc. Natl. Acad. Sci. U.S.A. 105(49), 19126–19131 (2008). [CrossRef] [PubMed]

]. Moreover, the use of time gates allows for simultaneous reconstruction of multiple fluorophores based on lifetime contrast without the use of unmixing algorithms [17

17. K. Suhling, P. M. French, and D. Phillips, “Time-resolved fluorescence microscopy,” Photochem. Photobiol. Sci. 4(1), 13–22 (2005). [CrossRef] [PubMed]

]. However, when considering time-resolved studies in small animals based on full TPSF, it is critical to employ an accurate light propagation model. Especially, for the early rising portion of the TPSF, a light propagation model that can accurately simulate the transition between minimally scattered photons and diffuse photons is required [18–22

18. K. Mitra and S. Kumar, “Development and comparison of models for light-pulse transport through scattering-absorbing media,” Appl. Opt. 38(1), 188–196 (1999). [CrossRef] [PubMed]

].

To date, analytical models, such as the radiative transport equation (RTE) and its approximation, the diffusion equation (DE), have been developed to solve the forward problem in FMT [21

21. M. Xu, W. Cai, M. Lax, and R. R. Alfano, “Photon migration in turbid media using a cumulant approximation to radiative transfer,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 65(6), 066609 (2002). [CrossRef] [PubMed]

, 22

22. W. Cai, M. Xu, and R. R. Alfano, “X. M, and R. Alfano, “Three-dimensional radiative transfer tomography for turbid media,” IEEE J. Sel. Top. Quantum Electron. 9(2), 189–198 (2003). [CrossRef]

]. However, the small finite volumes and their associated complex boundary conditions, the broad span of optical properties and potential low-scattering properties of certain organs may limit the validity of DE based models in small animals [23

23. J. C. Rasmussen, T. Pan, A. Joshi, T. Wareing, J. McGhee, and E. M. Sevick-Muraca, “Comparison of radiative transport, Monte Carlo, and diffusion forward models for small animal optical tomography,” in“2007 IEEE ISBI,” (2007), pp. 824–827.

]. Most importantly, the early rising part of the measured TPSF which contains both minimally scattered photons and diffuse photons [18

18. K. Mitra and S. Kumar, “Development and comparison of models for light-pulse transport through scattering-absorbing media,” Appl. Opt. 38(1), 188–196 (1999). [CrossRef] [PubMed]

, 20

20. M. Sakami, K. Mitra, and T. Vo-Dinh, “Analysis of short-pulse laser photon transport through tissues for optical tomography,” Opt. Lett. 27(5), 336–338 (2002). [CrossRef] [PubMed]

, 24

24. L. Wang, X. Liang, P. Galland, P. P. Ho, and R. R. Alfano, “True scattering coefficients of turbid matter measured by early-time gating,” Opt. Lett. 20(8), 913–915 (1995). [CrossRef] [PubMed]

] is not accurately modeled by the DE [15

15. M. Niedre and V. Ntziachristos, “Comparison of fluorescence tomographic imaging in mice with early-arriving and quasi-continuous-wave photons,” Opt. Lett. 35(3), 369–371 (2010). [CrossRef] [PubMed]

, 19

19. I. V. Yaroslavsky, A. N. Yaroslavsky, V. V. Tuchin, and H.-J. Schwarzmaier, “Effect of the scattering delay on time-dependent photon migration in turbid media,” Appl. Opt. 36(25), 6529–6538 (1997). [CrossRef] [PubMed]

, 21

21. M. Xu, W. Cai, M. Lax, and R. R. Alfano, “Photon migration in turbid media using a cumulant approximation to radiative transfer,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 65(6), 066609 (2002). [CrossRef] [PubMed]

, 22

22. W. Cai, M. Xu, and R. R. Alfano, “X. M, and R. Alfano, “Three-dimensional radiative transfer tomography for turbid media,” IEEE J. Sel. Top. Quantum Electron. 9(2), 189–198 (2003). [CrossRef]

]. Alternatively, the Monte Carlo (MC) method is considered an accurate light propagation model without any of the limitations of the DE to simulate pulse propagations in diffuse media [23

23. J. C. Rasmussen, T. Pan, A. Joshi, T. Wareing, J. McGhee, and E. M. Sevick-Muraca, “Comparison of radiative transport, Monte Carlo, and diffusion forward models for small animal optical tomography,” in“2007 IEEE ISBI,” (2007), pp. 824–827.

]. A number of studies focusing on Monte Carlo methods for fluorescence signal prediction have demonstrated its accuracy using synthetic and experimental data [25–27

25. A. Liebert, H. Wabnitz, N. Zołek, and R. Macdonald, “Monte Carlo algorithm for efficient simulation of time-resolved fluorescence in layered turbid media,” Opt. Express 16(17), 13188–13202 (2008). [CrossRef] [PubMed]

]. However, time-resolved Monte Carlo approaches have focused on spectroscopic applications but not tomographic ones. While inversion strategies based on Monte Carlo models in CW and FD have been reported [10

10. A. T. N. Kumar, S. B. Raymond, A. K. Dunn, B. J. Bacskai, and D. A. Boas, “A time domain fluorescence tomography system for small animal imaging,” IEEE Trans. Med. Imaging 27(8), 1152–1163 (2008). [CrossRef] [PubMed]

, 28

28. X. Zhang, C. Badea, M. Jacob, and G. A. Johnson, “Development of a noncontact 3-d fluorescence tomography system for small animal in vivo imaging,” SPIE 7171, 71910D (2009).

], no formulation for time-resolved fluorescence reconstruction based on Monte Carlo models has been yet reported.

In this work, we investigate the feasibility of performing time-resolved FMT of multiple fluorescence compounds based on Monte Carlo forward model and Time Gates. We derive a computationally efficient model that enables the calculation of the weight functions for both absorption perturbation and fluorophore distribution while simulating only the propagation of excitation photons in the tissue. We apply this new model to the simultaneous reconstruction of two fluorophores with lifetime contrast.We investigate for this purpose the information content imparted by single time gates and perform optical reconstructions based on multiple time gates to simultaneously image two fluorophores. The relative merits of different gate sampling strategies are compared and the validity of the model is experimentally established using a slab phantom with objects containing two commercial fluorophores.

2. Methodology

2.1. Time-resolved Monte Carlo forward model

η(r)=μafx(r)Φμax(r)
(1)

where Φ is the quantum yield. The time-resolved fluorescence signal without delay caused by the lifetime of the fluorophore is then expressed as:

UF(rs,rd,t)=Ωdr3W(rs,rd,r,t)η(r),
(2)

where the integration domain Ω is defined as the entire imaging volume, and the background weight function W′ is given by:

W(rs,rd,r,t)=0tdtA(rs,r,t)E(r,rd,tt).
(3)

Assuming single exponential time-decay for the fluorophore, the fluorescence emission delay can be taken into account by convolving the temporal fluorescence signals with the decay e−t/τ over time, where τ is the lifetime of the fluorophore. We can then write the fluorescence intensity measured at rd and time t for an impulsive excitation at rs and t 0 = 0 as:

UF(rs,rd,t)=0tdtUF(rs,rd,t)e(tt)τ.
(4)

This general expression is used to compute the fluorescence forward model. In the case of multiple fluorophores, the different fluorescent components can be independently computed and summed to obtain the complex fluorescence time resolved measurements. In all the in silico studies herein, the synthetic fluorescence measurements were based on this approach.

2.2. Calculation of the time-resolved Jacobians

The forward model described above can be casted into a tomography inverse problem to reconstruct the effective quantum yield. In order to compute the time-resolved Jacobians, or the spatial and temporal sensitivity maps with respect to η(r), Eq. (4) is written as:

UF(rs,rd,t)=Ωdr3W(rs,rd,r,t)η(r),
(5)

where W is the lifetime based weight function defined as:

W(rs,rd,r,t)=0tdtW(rs,rd,r,t)e(tt)τ,
(6)

Equations (3) and (6) give the general expression of the lifetime based weight matrix using Monte Carlo simulations. As Eq. (3) stands, the background weight function W′ can be calculated knowing A(rs, r, t) and E(r, rd, t) explicitly. This time convolution is extremely computationally intensive, and simulations with numerous photons at each position r in the region of interest are required to obtain statistically reliable calculation of E(r, rd, t). Thus this algorithm is computationally inefficient especially in FMT applications where a large number of voxels are considered. However, a more manageable formulation can be derived by considering an assumption commonly employed in FMT, namely that the scattering coefficients are identical at λx and λm. In the NIR spectral range, the scattering coefficient is expected to vary less than 20% over the Stokes shift gap of organic fluorophores [30

30. J. R. Mourant, J. P. Freyer, A. H. Hielscher, A. A. Eick, D. Shen, and T. M. Johnson, “Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics,” Appl. Opt. 37(16), 3586–3593 (1998). [CrossRef] [PubMed]

]. Furthermore, under the condition of isotropic scattering, it can be assumed that all fluorescence photons follow the trajectories of the excitation photons after generation [25

25. A. Liebert, H. Wabnitz, N. Zołek, and R. Macdonald, “Monte Carlo algorithm for efficient simulation of time-resolved fluorescence in layered turbid media,” Opt. Express 16(17), 13188–13202 (2008). [CrossRef] [PubMed]

]. Consider the ith photon in an excitation simulation that propagates from a discrete source point rs and then detected at the detector position rd at time t. Following the White Monte Carlo (WMC) [31

31. E. Alerstam, S. Andersson-Engels, and T. Svensson, “White Monte Carlo for time-resolved photon migration,” J. Biomed. Opt. 13(4), 041304 (2008). [CrossRef] [PubMed]

] approach, if the photon reaches r at time t s,r, the weight of the ith excitation photon wxi at r has decreased to:

wix(rs,r,ts,r)=wi,0exp(Σj=1piμax(rj)l(rj)),
(7)

where w i,0 is the initial weight of this photon, rj(j = 1, ..., pi) are the sub-regions that the photon consequently passes through from rs to r, and li(rj) is the path length that the photon passes at rj. Note that the total number of pi is photon dependent and the sub-regions for different photons are not necessarily the same. At r, the absorbed photon weight is given by:

Ai(rs,r,ts,r)=wix(rs,r,ts,r)(1exp(μax(r)li(r))).
(8)

The absorbed weight multiplied by η(r) is then converted to the initial weight of the fluorescence photon. We can calculate the final weight of the fluorescence photon detected at rd as:

wim(r,rd,tr,d)=Ai(rs,r,ts,r)η(r)exp(Σj=pi+1qiμam(rj)li(rj)),
(9)

where t r,d is the time that the photon passes from r to rd and rj(j = pi + 1, ..., qi) is the photon path from r to rd. Note that sum of t s,r and t r,d is the total time t that the excitation photon travels from rs till is detected, and that the sub-regions rj(j = 1, ..., qi) denote the total photon path. By summing up n detected emission photons for rs and rd, we can obtain UF defined in Eq. (2). Comparing Eq. (9) to Eq. (2), we have:

w(rs,rd,r,t)=Σi=1nAi(rs,r,t)exp(Σj=pi+1qiμam(rj)li(rj)).
(10)

Inserting Eqs. (7) and (8) to Eq. (10), then we have:

W(rs,rd,r,t)=Σi=1nwi,0exp(Σj=1piμax(rj)li(rj))×(1exp(μax(r)li(r))×exp(Σj=pi+1qiμam(rj)li(rj)).
(11)

According to Eq. (11), the background weight matrix for excitation source rs and detector rd at time t can be easily calculated using the photon paths and the absorption coefficients. This formulation should be employed in the case that the absorption coefficient at the excitation significantly differs from that at the emission wavelength, which might be the case for wavelengths below 700nm. However, the absorption optical spectra of bio-tissues is relatively flat in tissue in the NIR window [32

32. A. Roggan, K. Dorschel, O. Minet, D. Wolff, and G. Muller, “The optical properties of biological tissue in the near infrared wavelength range : review and measurements,” in Laser-induced Interstitial Thermotherapy (SPIE, Bellingham WA, Etats-Unis, 1995), pp. 10–44.

]. Thus for simplicity, and following standard assumption in FMT, we assume that the absorption coefficients at the excitation and emission wavelengths are identical, that is, μxa = μma. We therefore obtain:

W(rs,rd,r,t)=Σi=1nwix(rs,rd,t)(1exp(μax(r)li(r)),
(12)

where wxi(rs, rd, t) is the final weight of the ith excitation photon from rs and detected by rd at t, defined as:

wix(rs,rd,t)=wi,0exp(Σj=1qiμax(rj)li(rj)).
(13)

If μxaf(r)li(r) is small, which is the case in the voxelized geometry, from Taylor expansion, we can further simplify Eq. (12) to:

W(rs,rd,r,t)=Σi=1nwix(rs,rd,t)μax(r)li(r).
(14)

Integrating Eq. (14) over time results in a CW reconstruction algorithm similar to the method proposed by Zhang et al. [28

28. X. Zhang, C. Badea, M. Jacob, and G. A. Johnson, “Development of a noncontact 3-d fluorescence tomography system for small animal in vivo imaging,” SPIE 7171, 71910D (2009).

], that the photon paths weighted by their final fluence rate seen by the detector is the sensitivity function. Note also that the background weight matrix for fluorescence calculated by Eq. (14) is precisely the weight function for absorption perturbations [33

33. J. Chen and X. Intes, “Time-gated perturbation Monte Carlo for whole body functional imaging in small animals,” Opt. Express 17(22), 19566–19579 (2009). [CrossRef] [PubMed]

]. This indicates that the absorption and fluorophore distribution can be estimated by using the same weight matrix calculated by one excitation Monte Carlo simulation. Moreover, it is noteworthy that if there are different fluorophore species existing in the region of interest, only one excitation simulation is required to obtain the weight matrices Wk (k = 1, ..., NF, where NF is the number of species) for different lifetimes τk(k = 1, ..., NF) according to Eq. (6). This allows for a fast and efficient computational implementation to simultaneously reconstruct multiple fluorophores based on one forward simulation. Note that for simplicity, all the above equations use rs and rd as discrete points, but they can be extended to any illumination or detection strategies applied to complex boundaries by collecting the photons generated or detected inside the areas. In this study, we employed an original broad field illumination strategy which is a new technique that we are currently developing [34

34. J. Chen, V. Venugopal, F. Lesage, and X. Intes, “Time-resolved diffuse optical tomography with patterned-light illumination and detection,” Opt. Lett. 35(13), 2121–2123 (2010). [CrossRef] [PubMed]

], as described in Sec. 3 and 4.

2.3. Inverse problem

To cast the inverse problem, we employed a Born formulation in which the emission field is normalized by the excitation field [35

35. A. Soubret, J. Ripoll, and V. Ntziachristos, “Accuracy of fluorescent tomography in the presence of heterogeneities: study of the normalized Born ratio,” IEEE Trans. Med. Imaging 24(10), 1377–1386 (2005). [CrossRef] [PubMed]

] by normalizing the experimental time-domain emission measurements to the CW excitation flux at the same position. We therefore have:

Mm(rs,rd,t)Mx(rs,rd)=αUx(rs,rd)Σk=1NFd3rWk(rs,rd,r,t)ηk(r),
(15)

where α incorporates unknown constants associated with wavelength-dependent gains and attenuations that can be measured once for every imaging system, Mm(rs, rd, t) is the total signal from all fluorophores with different lifetimes at the emission wavelength measured at the boundary position rd and t excited from the source at rs, Mx(rs, rd) and Ux(rs, rd) are the measured and simulated, respectively, total excitation flux measured by a photo-detector at rd. This normalization efficiently mitigates the dependence of the detected fluorescent signal on the optical properties of the examined tissue [35

35. A. Soubret, J. Ripoll, and V. Ntziachristos, “Accuracy of fluorescent tomography in the presence of heterogeneities: study of the normalized Born ratio,” IEEE Trans. Med. Imaging 24(10), 1377–1386 (2005). [CrossRef] [PubMed]

, 36

36. C. Vinegoni, D. Razansky, J.-L. Figueiredo, M. Nahrendorf, V. Ntziachristos, and R. Weissleder, “Normalized Born ratio for fluorescence optical projection tomography,” Opt. Lett. 34(3), 319–321 (2009). [CrossRef] [PubMed]

], thus we do not model the absorptive heterogeneities associated with the different organs in our synthetic murine model.

The normalized fluorescence signals corresponding to different source-detector (S-D) pairs and time gates (left term in Eq. (15)) can be stacked up as a vector γRNm (Nm is the total number of measurements, that is, the product of the number of S-D pairs and the number of time gates) where it can be used as part of an inverse problem. Similarly, the kth effective quantum yield for the entire imaging volume can be written as ηkΩ=[ηk(r1),...,ηk(rNV)]T , where NV is the total number of voxels in the region of interest. Overall, the forward problem can be expressed in a matrix form:

[γ1γNm]=[β1W1,1Ωβ1W1,NFΩβNmWNm,1ΩβNmWNm,NFΩ][η1ΩηNFΩ],
(16)

where W Ω s,k is the weight function for the sth(s = 1, ..., Nm) measurement and the kth fluorophore species over the volume Ω, and βs = αs/Uxs is the normalization factor for the sth measurement and the corresponding weight function. The above equation has a general form Ax = b, then x=[η1Ω,...,ηNFΩ]T can be computed by solving this linear system of equations. Notice here the reconstructions for all fluorescent species provided herein are simultaneous based on monochromatic data and the measurements employed are direct measurements that are not derived from a pre-processed unmixing algorithm. Because of the ill-posedness of the inverse problem, the reconstruction method in fluorescence tomography usually sets up a least square optimization problem minimizing an objective function given by:

Ψ=Axb2+λ(r)x2,
(17)

where λ(r) is a spatially variant regularization parameter to compensate for the spatial dependence of the contrast and resolution in the reconstruction [37

37. B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt. 38(13), 2950–2961 (1999). [CrossRef] [PubMed]

].

2.4. Reconstruction algorithm

In summary, our approach to reconstruct the effective quantum yield is schematically depicted in Fig. 1. The capital letters shown in the square brackets below refer to the individual steps in this scheme. [A] represents a time-resolved Monte Carlo simulation for excitation photons, following the WMC approach to analytically compute the light transport in tissue. This simulation can be adapted to different illumination and detection strategies.

Fig. 1. Block diagram summarizing the steps to perform fluorescence reconstruction.

To calculate the background weight matrix in [B], we directly add the weighted paths to the weight matrix along with the Monte Carlo simulation, instead of storing all the photon trajectories. This reduces the post processing time for weight matrix calculation as well as the storage space for the huge number of photon paths. Note that any of Eqs. (11), (12) and (14) can be used in [B]. It should be noted that Eq. (14), wherein the absorption coefficient at λx and λm are assumed to be identical, can be used for computational efficiency. The lifetime-based weight matrices are then calculated by Eq. (6) in [C]. Knowing the experimental measurements, a conjugate gradient method is applied in [D] to minimize the objective function defined in Eq. (17) and obtain the fluorescence yield distribution. The fluorescence measurements in this study are either from an actual experiment or an emission Monte Carlo simulation at λm based on the absorption probability calculated in [A].

Fig. 2. (a) The mouse model geometry and arrangement of sources (top) and detectors (bottom) used for the simulations. The fluorescent kidneys are marked as red and black, according to different lifetimes. (b) The illumination patterns (red: on, black: off).

3. In silico study

3.1. Simulation setup

We first validated the proposed algorithm in silico in a model replicating small animal fluorescence imaging. The Monte Carlo simulations were performed on a mouse geometry with the kidneys and the skin shape extracted from a whole-body atlas [38

38. B. Dogdas, D. Stout, A. F. Chatziioannou, and R. M. Leahy, “Digimouse: a 3D whole body mouse atlas from CT and cryosection data,” Phys. Med. Biol. 52(3), 577–587 (2007). [CrossRef] [PubMed]

] (Fig. 2(a)). The entire field of view consisted of 89 × 33 × 19 voxels with size 1mm3. The optical properties were set to μxab = 0.3cm−1, μxs = 25cm−1, μma = 0.36cm−1, μms = 20cm−1 over the entire body, to simulate the different optical properties of mouse tissues at different wavelengths in the NIR window [32

32. A. Roggan, K. Dorschel, O. Minet, D. Wolff, and G. Muller, “The optical properties of biological tissue in the near infrared wavelength range : review and measurements,” in Laser-induced Interstitial Thermotherapy (SPIE, Bellingham WA, Etats-Unis, 1995), pp. 10–44.

]. The kidneys were considered to be labeled with two distinct fluorophores with lifetime of 0.5ns (black) and 1ns (red), respectively. The μaf for both kidneys were set to 0.1 times of the background absorption coefficient μxab, and the quantum yield was set to 1. 36 bar-shaped patterns were used as the source with each pattern illuminating half of the imaged surface along the x and y axes independently (Fig. 2(b)). A grid of 64 point detectors in transmission geometry were arranged covering a 31mm×27mm×18mm volume spanning the abdomen. The sources and detectors are projected to the surface along the z axis. The initial positions of the photons were randomly spanned (that is, uniformly distributed) over the illuminated area to accommodate the broad field illumination. The detectors had a 2mm separation along the y axis and a 2.2mm separation along the x axis with a radius of 1mm.

In the Monte Carlo simulation, 109 photons were consecutively launched for each pattern source on a massively parallel environment (blue gene, CCNI at RPI [39

39. “Computational Center for Nanotechnology Innovations (CCNI),” http://www.rpi.edu/research/ccni/.

]) to calculate the excitation and fluorescence measurements as described in Sec. 2.1. Under the assumption that the optical properties are equal at λx and λm (Eq. (14)), the time-resolved Jacobians for FMT were simultaneously computed by convolving the excitation field Jacobians by the lifetime decay. The temporal profile was recorded over a 6ns time window with 20ps resolution and a gate-width of 200ps (total of 120 gates). These parameters were selected to replicate the operating characteristics of our pre-clinical imaging platform [16

16. V. Venugopal, J. Chen, F. Lesage, and X. Intes, “Full-field time-resolved fluorescence tomography of small animals,” Opt. Lett. 35(19), 3189–3191 (2010). [CrossRef] [PubMed]

] and correspond to the typical minimum gate width. Note that we do not explicitly include noise in the simulation data. However, the Poisson noise inherent in time-resolved studies also exists in Monte Carlo simulations. As the excitation measurements and fluorescence measurements are both directly simulated in this work (not computed through the multiplication of a Jacobian with a synthetic phantom), both measurements have an uncorrelated Poisson noise that replicates experimental conditions.

Fig. 3. (a) Representative simulated signals from the fluorophores with τ 1 = 0.5ns (blue) and τ 2 = 1ns (red), and the total fluorescence signal (black) used in the reconstructions. All the signals are generated using the mouse phantom and a single transmittance S-D pair as shown in (b) and (c). (b) Weight matrices for τ = 1ns at the 3 gates (gray) presented in (a). (c) The direct detector readings from the region of interest (black boxes) from the simulation on the mouse phantom. The bars in (b) and white boxes in (c) indicate the illumination pattern and the crosses in (b) and (c) indicate the detector used in (a).

The average calculation time for an excitation or emission MC simulation for one pattern source and all detectors and gates was less than 8 minutes (on 4096 nodes). Figure 3 shows typical simulated fluorescence signals and the background weight matrices calculated using Eq. (14) at different time gates.

3.2. Gate information content for fluorophore multiplexing

In order to quantitatively assess the information content at single gates, we need to establish quantitative performance metrics. We investigated the quantification, crosstalk and resolution based on the reconstructions using one gate from 0.2ns to 2.9ns (opening of the gate) at an interval of 0.1ns. The quantification of the reconstructed η 1 is defined as Q 1 = max[η 11)]/H 11), where Ω1 denotes the known location of the inclusions with lifetime τ 1 = 0.5ns and H 1 denotes the expected value of η 1. The crosstalk is defined as X 1 = max[η 21)]/max[η 11)] to quantify the separability of the two inclusions with lifetime contrast. The quantification and yield crosstalk for the 1ns component was similarly evaluated. Note that Q and X are overestimated using the definitions based on only maximum value in the region of interest. The spatial resolution is measured as the full width at half maximum (FWHM) of the point spread functions (PSF). PSFs are created by simulating a perturbation at a single point. The simulated perturbation, xsim, is a vector of zeros except for one target pixel at the center of the reconstructed volume with a value of one. We then generated simulated measurements ysim = Axsim and reconstructed the image. The cube root of the volume of voxels within half the peak reconstructed value is defined as FWHM. Plots of the quantification, crosstalk and FWHM of PSFs are shown in Fig. 4. Note that for all these reconstructions, the inverse problem size was identical and that the same iteration number (150) was used in the CG algorithm. This ensured consistency between the different reconstructions.

Fig. 4. Information content evaluation based on single gate reconstructions in terms of (a) quantification, (b) crosstalk and (c) resolution.

From the plots of performance metrics, it is clear that for the reconstructed object with the shorter lifetime, the quantitative accuracy is decreasing with time while the crosstalk is increasing. Conversely, for the reconstructed object with the longer lifetime, an opposite behavior is observed with equivalent quantification and crosstalk around the maximum gate. It is noteworthy that the reconstruction using the latest gate (t = 2.8ns) results in the most accurate quantification and the least crosstalk for the compound with the longer lifetime, implying the necessity of the late gates for effective separation of the objects with lifetime contrast. It is impossible to achieve the best quantification and minimum crosstalk using only the early gates because of the comparable contribution of the two fluorophores at the early time points (Fig. 3). As expected, the steadily increasing FWHM of PSFs with time indicates that the early gates provide improved resolution when compared to the late gates [7

7. V. Y. Soloviev, C. D’Andrea, G. Valentini, R. Cubeddu, and S. R. Arridge, “Combined reconstruction of fluorescent and optical parameters using time-resolved data,” Appl. Opt. 48(1), 28–36 (2009). [CrossRef] [PubMed]

, 15

15. M. Niedre and V. Ntziachristos, “Comparison of fluorescence tomographic imaging in mice with early-arriving and quasi-continuous-wave photons,” Opt. Lett. 35(3), 369–371 (2010). [CrossRef] [PubMed]

]. The better resolution at early gates is due to the narrower probing volume by photons detected at the early time points, as seen in Fig. 3(b). However, there is no significant reduction in the resolution of the reconstructed yield distributions for the two lifetimes when using gates later than the maximum gate. The reconstruction at the 10th gate (t = 0.2ns) does not follow the trend due to the poor statistics of the weight matrix at the early gates, where the number of detected photons is insufficient to generate reliable statistics. It results in increased artifacts in the reconstructions when using very early gates. From the above analysis, we conclude that multiple fluorophores with lifetime contrast cannot be efficiently separated with minimal crosstalk if a single gate of the TPSF is used.

3.3. Lifetime multiplexing tomography based on multiple gates

To utilize the various information carried by different gates, we investigated reconstructions using multiple gates simultaneously. In this regard, 5 sets of time gated measurements consisting of 5 gates each were considered. The number of gates was limited to 5 to reduce the memory burden as all the gates are simultaneously used in one inverse problem. The 5 different sets simulated in this study are illustrated in Fig. 5(a). Sets 1, 3, 5 have gates spaced at uniform intervals covering different time ranges, corresponding to the early to the maximum gate, full span, and the maximum to the end of the TPSFs. Sets 2 and 4 span the full time range, with non-uniform spacings exponentially scaled towards the rising and the decaying edges of the TPSFs respectively. This is a heuristic but straightforward approach to investigate the reconstruction performance using multiple combined gates. The reconstructions using combined gates of set 1 and set 5 are displayed in Fig. 5(b) and 5(c), respectively.

Fig. 5. (a) Sets of time gates investigated in multi-gate reconstructions. (b) and (c) display the reconstructions using combined gates of set 1 and set 5, respectively.

In order to ascertain the effect of assuming equal optical properties at excitation and emission wavelengths, a second set of time-resolved Jacobians were computed for the same set of gates using Eq. (11) (A–E). The reconstruction results when using sets are also given in Table 1. Comparing sets 1–5 and sets A–E, it is determined that the assumption of equal optical properties at excitation and emission wavelengths results in less than 5% quantification error for all three criteria.

Table 1. Quantitative Comparison of the Multi-Gate Reconstructions

table-icon
View This Table

4. Experiment

4.1. Diffuse optical tomography system

The lifetime based Monte Carlo algorithm was experimentally validated by using a time-domain small animal molecular imaging system, which is schematically depicted in Fig. 6(a) [16

16. V. Venugopal, J. Chen, F. Lesage, and X. Intes, “Full-field time-resolved fluorescence tomography of small animals,” Opt. Lett. 35(19), 3189–3191 (2010). [CrossRef] [PubMed]

, 40

40. V. Venugopal, J. Chen, and X. Intes, “Development of an optical imaging platform for functional imaging of small animals using wide-field excitation,” Biomed. Opt. Express 1(1), 143–156 (2010). [CrossRef] [PubMed]

]. A tunable (690nm - 1020nm) Ti-Sapphire laser is used as the light source in the system. The laser source is expanded and projected on a digital micro-mirror device (DmD) based DLP board (Discovery 1100, Texas Instruments). 36 binary sliding-bar patterns similar to Fig. 2(b) are reflected on the programmable DmD chip and projected on the imaging chamber to provide an illumination area of 35mm × 20mm (matching the dimensions of a small animal torso). The transmitted signal is detected by a time-gated intensified CCD (ICCD) camera (Picostar HR, LaVision) with a time resolution of 40ps in the interest of a shorter acquisition time (24 minutes total using 36 patterns at excitation and emission wavelengths). A gate width of 300ps was selected to avoid temporal errors due to jitter when employing 200ps gates [41

41. C. D’Andrea, D. Comelli, A. Pifferi, A. Torricelli, G. Valentini, and R. Cubeddu, “Time-resolved optical imaging through turbid media using a fast data acquisition system based on a gated CCD camera,” J. Phys. D Appl. Phys. 36(14), 1675–1681 (2003). [CrossRef]

]. Photon measurements at the excitation and fluorescence wavelengths were acquired consecutively using bandpass filters. The incident fields for excitation and the fluorescence data for excitation were collected at 780nm and 832nm, respectively. The laser power was independently optimized for both the excitation and fluorescence measurements to maintain a maximum signal level of 4000 counts. Further, the acquired images are post-processed to generate 1mm × 1mm detector measurements by averaging over 7×7 pixels.

4.2. Phantom setup

A polycarbonate slab phantom (Fig. 6(b)) that consisted of two tubes having an inner diameter of 1.5mm was used to experimentally assess the performance of the algorithm. To achieve lifetime contrast, tube I, II were filled with 2.5nmol of Indocyanine Green (τ 1 = 0.45ns) and IRDye 800CW (τ 2 = 0.8ns) dissolved in 180µL ethanol, respectively. The effective quantum yields of the two tubes was externally calibrated and tube I was found to have 1.5 times the yield of tube II. The tank was filled with a mixture of Intralipid-20% (Sigma-Aldrich) and a water-soluble NIR Dye (Epolight 2717, Epolin) diluted in water to simulate the average optical properties (μa = 0.16cm−1 and μs = 17cm−1) of murine models.

Fig. 6. (a) Schematic of the TD fluorescence tomography system. LS: laser; PC: power control; FBC: fiber-beam coupler; BE: beam expander; DMD: DmD-based DLP chip; IL: imaging lens; IC: imaging chamber; CCD: intensified CCD camera; HRI: high rate imager; TDU: trigger delay unit; TG: trigger generator. (b) Phantom setup. The red box shows the pattern illumination area. All dimensions are in mm.
Fig. 7. Experiment results using (a) TD data and (b) CW data. The images are normalized to the maximum of the reconstruction.

4.3. Experimental reconstruction

For the MC computation, the phantom was digitized into a 28 × 40 × 14 volume image with 1mm × 1mm × 1mm voxels. A background weight matrix as defined in Eq. (14) covering the field of view of 24 × 34 × 14 area, resulting in a total of 11,424 voxels, was calculated. The accuracy of the inverse model with experimental data is affected by the heterogeneous spatial distribution of the experimental pattern on the phantom. The illumination spatial heterogeneity is taken into account in the Monte Carlo simulation by scaling the initial photon weight to the illumination intensity at the injection position using experimental calibration patterns. The system impulse response functions (IRF) for all sources was convolved into the model before tomographic inversion, to match the model prediction with the experimental measurements. Since the lifetimes are known in advance in this controlled phantom experiment, the weight functions are generated using the shorter lifetime at τ 1 = 0.45ns (ICG) and the longer lifetime at τ 2 = 0.8ns (IRdye™) using Eq. (6). In in vivo applications, the lifetime of the dye within a control animal should be separately estimated due to the possible fluctuation of the fluorophore lifetime under different micro environments.

To achieve better separation between the two fluorophores while maintaining resolution, we selected the measurements at 20% of the rising edge, 75%, 35%, 20% and 15% of the decaying edge of the convolved TPSFs (mimicking gate set 4 in Sec. 3) to recover the fluorescence yield distribution based on the analysis of the simulated result. Gates with less than 400 photons (10% of the max for transmittance detector) were discarded in the reconstruction process to maintain an appropriate signal-to-noise ratio (SNR). This SNR led to a reduced size of the dataset by 8%. For comparison, we also reconstructed the fluorophore distribution using CW data at the same emission wavelength. The CW weight functions for the inversion were calculated by integrating the corresponding TD weight functions over time and the TPSFs were integrated over time to obtain the CW datatype.

The reconstruction results are shown in Fig. 7(a), similar as the 3D visualization of the in silico results. The two objects are separated with crosstalk X 1 = 28.51% and X 2 = 26.24%. The ratio of the mean reconstructed yields within the 50% isovolume of tube I and II was found to be 1.52. The average dimension of the reconstructed tubes was 34.5% larger along the x axis and 85% larger along the z axis (depth). This difference in resolution is due to the transmittance arrangement of the S-D pairs, which smears the resolution along the z axis. Hence, the diameter of the reconstructed isosurface is approximately 2 times bigger than the actual diameter due to the loss in the depth dimension, which can be attributed to the high scattering, and to the position of the two objects in the middle of the phantom, where the weight matrix has the least sensitivity. The CW reconstruction using one wavelength presented in Fig. 7(b), shows poor resolution (103% increase in depth) and no separation (100% crosstalk).

5. Discussion

We validated the proposed formulation by reconstructing two different compounds simultaneously based on one excitation/emission wavelength pair without using an lifetime based unmixing algorithm on the diffuse measurements. We employed the time-gated datatype to retain the appeal of time-resolved measurements, improved resolution and lifetime contrast unmixing. While using early gates increases the resolution, the crosstalk between fluorophore is significant due to their equal contribution to the signal when used alone. The introduction of late gates in reconstruction is required to improve separation and quantification of fluorescent probes with lifetime contrast. The combination of early and late gates allows one to trade off the advantages and disadvantages of different time gates and thereby obtain enhanced resolution and crosstalk simultaneously. Employing 5 gates, the fluorescent inclusions were accurately resolved in the region of interest with minimal crosstalk. However, not all the individual metrics are optimal as an emphasis was put on reduction of crosstalk in this work due to the particular application considered. In the case of applications with only one compound, it will be probably possible to identify a set of gates that allow superior resolution and quantification [15

15. M. Niedre and V. Ntziachristos, “Comparison of fluorescence tomographic imaging in mice with early-arriving and quasi-continuous-wave photons,” Opt. Lett. 35(3), 369–371 (2010). [CrossRef] [PubMed]

]. Moreover, a reconstruction scheme, in which the spatial a priori information derived from an early gate reconstruction is incorporated in reconstruction based only on late gates (or even CW datatype for enhanced SNR), could provide superior results. Also, in the experiment conducted herein, the two compounds were spatially separated. The method proposed here however does allow the reconstruction of the concentration of two fluorophores that are co-localized. This cannot be achieved when employing methods that reconstruct both quantum yield and lifetime. And this provides a new method to monitor the multiple markers in vivo.

6. Conclusion

In this work we have developed a new formulation to perform time resolved FMT based on the Monte Carlo forward model. The formulation derived allows to compute fluorescence Jacobians within the fluorescence Born formulation in a computationally efficient manner. Based on standard assumption in the field of FMT, i.e. that the optical properties are similar at the excitation and emission wavelength, the time-resolved fluorescence Jacobians for multiple fluorophores can be computed from one excitation simulation.We have investigated the use of this new formulation in the context of simultaneous tomographic imaging of multiple fluorophores without using any unmixing algorithm. We considered the use of time-gated datatypes spanning the full TPSFs, including the early gates for which the diffusion equation is known to be inaccurate when compared to the MC. The in silico and experimental studies demonstrate that this new formulation can simultaneously reconstruct two fluorophores with improved resolution and reduced crosstalk when multiple gates are employed, whereas the two fluorophores cannot be resolved if only early gates are employed.

Acknowledgment

We gratefully acknowledge the technical support of the Computational Center for Nanotechnology Innovations (CCNI) at Rensselaer Polytechnic Institute.

References and links

1.

V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005). [CrossRef] [PubMed]

2.

E. M. Sevick-Muraca and J. C. Rasmussen, “Molecular imaging with optics: primer and case for near-infrared fluorescence techniques in personalized medicine,” J. Biomed. Opt. 13(4), 041303 (2008). [CrossRef] [PubMed]

3.

A. H. Hielscher, “Optical tomographic imaging of small animals,” Curr. Opin. Biotechnol. 16(1), 79–88 (2005). [CrossRef] [PubMed]

4.

J. R. Mansfield, “Distinguished photons: a review of in vivo spectral fluorescence imaging in small animals,” Curr. Pharm. Biotechnol. 11(6), 628–638 (2010). [CrossRef] [PubMed]

5.

R. Cubeddu, D. Comelli, C. D’Andrea, P. Taroni, and G. Valentini, “Time-resolved fluorescence imaging in biology and medicine,” J. Phys. D Appl. Phys. 35(9), R61–R76 (2002). [CrossRef]

6.

L. Zhang, F. Gao, H. He, and H. Zhao, “Three-dimensional scheme for time-domain fluorescence molecular tomography based on Laplace transforms with noise-robust factors,” Opt. Express 16(10), 7214–7223 (2008). [CrossRef] [PubMed]

7.

V. Y. Soloviev, C. D’Andrea, G. Valentini, R. Cubeddu, and S. R. Arridge, “Combined reconstruction of fluorescent and optical parameters using time-resolved data,” Appl. Opt. 48(1), 28–36 (2009). [CrossRef] [PubMed]

8.

S. Lam, F. Lesage, and X. Intes, “Time domain fluorescent diffuse optical tomography: analytical expressions,” Opt. Express 13(7), 2263–2275 (2005). [CrossRef] [PubMed]

9.

A. T. N. Kumar, S. B. Raymond, G. Boverman, D. A. Boas, and B. J. Bacskai, “Time resolved fluorescence tomography of turbid media based on lifetime contrast,” Opt. Express 14(25), 12255–12270 (2006). [CrossRef] [PubMed]

10.

A. T. N. Kumar, S. B. Raymond, A. K. Dunn, B. J. Bacskai, and D. A. Boas, “A time domain fluorescence tomography system for small animal imaging,” IEEE Trans. Med. Imaging 27(8), 1152–1163 (2008). [CrossRef] [PubMed]

11.

M. Köllner and J. Wolfrum, “How many photons are necessary for fluorescence-lifetime measurements?” Chem. Phys. Lett. 200(1-2), 199–204 (1992). [CrossRef]

12.

F. Liu, K. M. Yoo, and R. R. Alfano, “Ultrafast laser-pulse transmission and imaging through biological tissues,” Appl. Opt. 32(4), 554–558 (1993). [CrossRef] [PubMed]

13.

J. Wu, L. Perelman, R. R. Dasari, and M. S. Feld, “Fluorescence tomographic imaging in turbid media using early-arriving photons and Laplace transforms,” Proc. Natl. Acad. Sci. U.S.A. 94(16), 8783–8788 (1997). [CrossRef] [PubMed]

14.

M. J. Niedre, R. H. de Kleine, E. Aikawa, D. G. Kirsch, R. Weissleder, and V. Ntziachristos, “Early photon tomography allows fluorescence detection of lung carcinomas and disease progression in mice in vivo,” Proc. Natl. Acad. Sci. U.S.A. 105(49), 19126–19131 (2008). [CrossRef] [PubMed]

15.

M. Niedre and V. Ntziachristos, “Comparison of fluorescence tomographic imaging in mice with early-arriving and quasi-continuous-wave photons,” Opt. Lett. 35(3), 369–371 (2010). [CrossRef] [PubMed]

16.

V. Venugopal, J. Chen, F. Lesage, and X. Intes, “Full-field time-resolved fluorescence tomography of small animals,” Opt. Lett. 35(19), 3189–3191 (2010). [CrossRef] [PubMed]

17.

K. Suhling, P. M. French, and D. Phillips, “Time-resolved fluorescence microscopy,” Photochem. Photobiol. Sci. 4(1), 13–22 (2005). [CrossRef] [PubMed]

18.

K. Mitra and S. Kumar, “Development and comparison of models for light-pulse transport through scattering-absorbing media,” Appl. Opt. 38(1), 188–196 (1999). [CrossRef] [PubMed]

19.

I. V. Yaroslavsky, A. N. Yaroslavsky, V. V. Tuchin, and H.-J. Schwarzmaier, “Effect of the scattering delay on time-dependent photon migration in turbid media,” Appl. Opt. 36(25), 6529–6538 (1997). [CrossRef] [PubMed]

20.

M. Sakami, K. Mitra, and T. Vo-Dinh, “Analysis of short-pulse laser photon transport through tissues for optical tomography,” Opt. Lett. 27(5), 336–338 (2002). [CrossRef] [PubMed]

21.

M. Xu, W. Cai, M. Lax, and R. R. Alfano, “Photon migration in turbid media using a cumulant approximation to radiative transfer,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 65(6), 066609 (2002). [CrossRef] [PubMed]

22.

W. Cai, M. Xu, and R. R. Alfano, “X. M, and R. Alfano, “Three-dimensional radiative transfer tomography for turbid media,” IEEE J. Sel. Top. Quantum Electron. 9(2), 189–198 (2003). [CrossRef]

23.

J. C. Rasmussen, T. Pan, A. Joshi, T. Wareing, J. McGhee, and E. M. Sevick-Muraca, “Comparison of radiative transport, Monte Carlo, and diffusion forward models for small animal optical tomography,” in“2007 IEEE ISBI,” (2007), pp. 824–827.

24.

L. Wang, X. Liang, P. Galland, P. P. Ho, and R. R. Alfano, “True scattering coefficients of turbid matter measured by early-time gating,” Opt. Lett. 20(8), 913–915 (1995). [CrossRef] [PubMed]

25.

A. Liebert, H. Wabnitz, N. Zołek, and R. Macdonald, “Monte Carlo algorithm for efficient simulation of time-resolved fluorescence in layered turbid media,” Opt. Express 16(17), 13188–13202 (2008). [CrossRef] [PubMed]

26.

J. Swartling, A. Pifferi, A. M. K. Enejder, and S. Andersson-Engels, “Accelerated Monte Carlo models to simulate fluorescence spectra from layered tissues,” J. Opt. Soc. Am. A 20(4), 714–727 (2003). [CrossRef] [PubMed]

27.

D. A. Boas, J. P. Culver, J. J. Stott, and A. K. Dunn, “Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head,” Opt. Express 10(3), 159–170 (2002). [PubMed]

28.

X. Zhang, C. Badea, M. Jacob, and G. A. Johnson, “Development of a noncontact 3-d fluorescence tomography system for small animal in vivo imaging,” SPIE 7171, 71910D (2009).

29.

A. J. Welch, C. M. Gardner, R. Richards-Kortum, E. Chan, G. Criswell, J. Pfefer, and S. Warren, “Propagation of fluorescent light,” Lasers Surg. Med. 21(2), 166–178 (1997). [CrossRef] [PubMed]

30.

J. R. Mourant, J. P. Freyer, A. H. Hielscher, A. A. Eick, D. Shen, and T. M. Johnson, “Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics,” Appl. Opt. 37(16), 3586–3593 (1998). [CrossRef] [PubMed]

31.

E. Alerstam, S. Andersson-Engels, and T. Svensson, “White Monte Carlo for time-resolved photon migration,” J. Biomed. Opt. 13(4), 041304 (2008). [CrossRef] [PubMed]

32.

A. Roggan, K. Dorschel, O. Minet, D. Wolff, and G. Muller, “The optical properties of biological tissue in the near infrared wavelength range : review and measurements,” in Laser-induced Interstitial Thermotherapy (SPIE, Bellingham WA, Etats-Unis, 1995), pp. 10–44.

33.

J. Chen and X. Intes, “Time-gated perturbation Monte Carlo for whole body functional imaging in small animals,” Opt. Express 17(22), 19566–19579 (2009). [CrossRef] [PubMed]

34.

J. Chen, V. Venugopal, F. Lesage, and X. Intes, “Time-resolved diffuse optical tomography with patterned-light illumination and detection,” Opt. Lett. 35(13), 2121–2123 (2010). [CrossRef] [PubMed]

35.

A. Soubret, J. Ripoll, and V. Ntziachristos, “Accuracy of fluorescent tomography in the presence of heterogeneities: study of the normalized Born ratio,” IEEE Trans. Med. Imaging 24(10), 1377–1386 (2005). [CrossRef] [PubMed]

36.

C. Vinegoni, D. Razansky, J.-L. Figueiredo, M. Nahrendorf, V. Ntziachristos, and R. Weissleder, “Normalized Born ratio for fluorescence optical projection tomography,” Opt. Lett. 34(3), 319–321 (2009). [CrossRef] [PubMed]

37.

B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt. 38(13), 2950–2961 (1999). [CrossRef] [PubMed]

38.

B. Dogdas, D. Stout, A. F. Chatziioannou, and R. M. Leahy, “Digimouse: a 3D whole body mouse atlas from CT and cryosection data,” Phys. Med. Biol. 52(3), 577–587 (2007). [CrossRef] [PubMed]

39.

“Computational Center for Nanotechnology Innovations (CCNI),” http://www.rpi.edu/research/ccni/.

40.

V. Venugopal, J. Chen, and X. Intes, “Development of an optical imaging platform for functional imaging of small animals using wide-field excitation,” Biomed. Opt. Express 1(1), 143–156 (2010). [CrossRef] [PubMed]

41.

C. D’Andrea, D. Comelli, A. Pifferi, A. Torricelli, G. Valentini, and R. Cubeddu, “Time-resolved optical imaging through turbid media using a fast data acquisition system based on a gated CCD camera,” J. Phys. D Appl. Phys. 36(14), 1675–1681 (2003). [CrossRef]

42.

X. Cheng and D. Boas, “Systematic diffuse optical image errors resulting from uncertainty in the background optical properties,” Opt. Express 4(8), 299–307 (1999). [CrossRef] [PubMed]

OCIS Codes
(110.6960) Imaging systems : Tomography
(170.3010) Medical optics and biotechnology : Image reconstruction techniques
(170.3650) Medical optics and biotechnology : Lifetime-based sensing
(170.6920) Medical optics and biotechnology : Time-resolved imaging
(110.1758) Imaging systems : Computational imaging

ToC Category:
Image Reconstruction and Inverse Problems

History
Original Manuscript: January 4, 2011
Revised Manuscript: February 25, 2011
Manuscript Accepted: February 26, 2011
Published: March 14, 2011

Citation
Jin Chen, Vivek Venugopal, and Xavier Intes, "Monte Carlo based method for fluorescence tomographic imaging with lifetime multiplexing using time gates," Biomed. Opt. Express 2, 871-886 (2011)
http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-2-4-871


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. V. Ntziachristos, J. Ripoll, L. V. Wang, and R. Weissleder, “Looking and listening to light: the evolution of whole-body photonic imaging,” Nat. Biotechnol. 23(3), 313–320 (2005). [CrossRef] [PubMed]
  2. E. M. Sevick-Muraca and J. C. Rasmussen, “Molecular imaging with optics: primer and case for near-infrared fluorescence techniques in personalized medicine,” J. Biomed. Opt. 13(4), 041303 (2008). [CrossRef] [PubMed]
  3. A. H. Hielscher, “Optical tomographic imaging of small animals,” Curr. Opin. Biotechnol. 16(1), 79–88 (2005). [CrossRef] [PubMed]
  4. J. R. Mansfield, “Distinguished photons: a review of in vivo spectral fluorescence imaging in small animals,” Curr. Pharm. Biotechnol. 11(6), 628–638 (2010). [CrossRef] [PubMed]
  5. R. Cubeddu, D. Comelli, C. D’Andrea, P. Taroni, and G. Valentini, “Time-resolved fluorescence imaging in biology and medicine,” J. Phys. D Appl. Phys. 35(9), R61–R76 (2002). [CrossRef]
  6. L. Zhang, F. Gao, H. He, and H. Zhao, “Three-dimensional scheme for time-domain fluorescence molecular tomography based on Laplace transforms with noise-robust factors,” Opt. Express 16(10), 7214–7223 (2008). [CrossRef] [PubMed]
  7. V. Y. Soloviev, C. D’Andrea, G. Valentini, R. Cubeddu, and S. R. Arridge, “Combined reconstruction of fluorescent and optical parameters using time-resolved data,” Appl. Opt. 48(1), 28–36 (2009). [CrossRef] [PubMed]
  8. S. Lam, F. Lesage, and X. Intes, “Time domain fluorescent diffuse optical tomography: analytical expressions,” Opt. Express 13(7), 2263–2275 (2005). [CrossRef] [PubMed]
  9. A. T. N. Kumar, S. B. Raymond, G. Boverman, D. A. Boas, and B. J. Bacskai, “Time resolved fluorescence tomography of turbid media based on lifetime contrast,” Opt. Express 14(25), 12255–12270 (2006). [CrossRef] [PubMed]
  10. A. T. N. Kumar, S. B. Raymond, A. K. Dunn, B. J. Bacskai, and D. A. Boas, “A time domain fluorescence tomography system for small animal imaging,” IEEE Trans. Med. Imaging 27(8), 1152–1163 (2008). [CrossRef] [PubMed]
  11. M. Köllner and J. Wolfrum, “How many photons are necessary for fluorescence-lifetime measurements?” Chem. Phys. Lett. 200(1-2), 199–204 (1992). [CrossRef]
  12. F. Liu, K. M. Yoo, and R. R. Alfano, “Ultrafast laser-pulse transmission and imaging through biological tissues,” Appl. Opt. 32(4), 554–558 (1993). [CrossRef] [PubMed]
  13. J. Wu, L. Perelman, R. R. Dasari, and M. S. Feld, “Fluorescence tomographic imaging in turbid media using early-arriving photons and Laplace transforms,” Proc. Natl. Acad. Sci. U.S.A. 94(16), 8783–8788 (1997). [CrossRef] [PubMed]
  14. M. J. Niedre, R. H. de Kleine, E. Aikawa, D. G. Kirsch, R. Weissleder, and V. Ntziachristos, “Early photon tomography allows fluorescence detection of lung carcinomas and disease progression in mice in vivo,” Proc. Natl. Acad. Sci. U.S.A. 105(49), 19126–19131 (2008). [CrossRef] [PubMed]
  15. M. Niedre and V. Ntziachristos, “Comparison of fluorescence tomographic imaging in mice with early-arriving and quasi-continuous-wave photons,” Opt. Lett. 35(3), 369–371 (2010). [CrossRef] [PubMed]
  16. V. Venugopal, J. Chen, F. Lesage, and X. Intes, “Full-field time-resolved fluorescence tomography of small animals,” Opt. Lett. 35(19), 3189–3191 (2010). [CrossRef] [PubMed]
  17. K. Suhling, P. M. French, and D. Phillips, “Time-resolved fluorescence microscopy,” Photochem. Photobiol. Sci. 4(1), 13–22 (2005). [CrossRef] [PubMed]
  18. K. Mitra and S. Kumar, “Development and comparison of models for light-pulse transport through scattering-absorbing media,” Appl. Opt. 38(1), 188–196 (1999). [CrossRef] [PubMed]
  19. I. V. Yaroslavsky, A. N. Yaroslavsky, V. V. Tuchin, and H.-J. Schwarzmaier, “Effect of the scattering delay on time-dependent photon migration in turbid media,” Appl. Opt. 36(25), 6529–6538 (1997). [CrossRef] [PubMed]
  20. M. Sakami, K. Mitra, and T. Vo-Dinh, “Analysis of short-pulse laser photon transport through tissues for optical tomography,” Opt. Lett. 27(5), 336–338 (2002). [CrossRef] [PubMed]
  21. M. Xu, W. Cai, M. Lax, and R. R. Alfano, “Photon migration in turbid media using a cumulant approximation to radiative transfer,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 65(6), 066609 (2002). [CrossRef] [PubMed]
  22. W. Cai, M. Xu, and R. R. Alfano, “X. M, and R. Alfano, “Three-dimensional radiative transfer tomography for turbid media,” IEEE J. Sel. Top. Quantum Electron. 9(2), 189–198 (2003). [CrossRef]
  23. J. C. Rasmussen, T. Pan, A. Joshi, T. Wareing, J. McGhee, and E. M. Sevick-Muraca, “Comparison of radiative transport, Monte Carlo, and diffusion forward models for small animal optical tomography,” in“2007 IEEE ISBI,” (2007), pp. 824–827.
  24. L. Wang, X. Liang, P. Galland, P. P. Ho, and R. R. Alfano, “True scattering coefficients of turbid matter measured by early-time gating,” Opt. Lett. 20(8), 913–915 (1995). [CrossRef] [PubMed]
  25. A. Liebert, H. Wabnitz, N. Zołek, and R. Macdonald, “Monte Carlo algorithm for efficient simulation of time-resolved fluorescence in layered turbid media,” Opt. Express 16(17), 13188–13202 (2008). [CrossRef] [PubMed]
  26. J. Swartling, A. Pifferi, A. M. K. Enejder, and S. Andersson-Engels, “Accelerated Monte Carlo models to simulate fluorescence spectra from layered tissues,” J. Opt. Soc. Am. A 20(4), 714–727 (2003). [CrossRef] [PubMed]
  27. D. A. Boas, J. P. Culver, J. J. Stott, and A. K. Dunn, “Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head,” Opt. Express 10(3), 159–170 (2002). [PubMed]
  28. X. Zhang, C. Badea, M. Jacob, and G. A. Johnson, “Development of a noncontact 3-d fluorescence tomography system for small animal in vivo imaging,” SPIE 7171, 71910D (2009).
  29. A. J. Welch, C. M. Gardner, R. Richards-Kortum, E. Chan, G. Criswell, J. Pfefer, and S. Warren, “Propagation of fluorescent light,” Lasers Surg. Med. 21(2), 166–178 (1997). [CrossRef] [PubMed]
  30. J. R. Mourant, J. P. Freyer, A. H. Hielscher, A. A. Eick, D. Shen, and T. M. Johnson, “Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics,” Appl. Opt. 37(16), 3586–3593 (1998). [CrossRef] [PubMed]
  31. E. Alerstam, S. Andersson-Engels, and T. Svensson, “White Monte Carlo for time-resolved photon migration,” J. Biomed. Opt. 13(4), 041304 (2008). [CrossRef] [PubMed]
  32. A. Roggan, K. Dorschel, O. Minet, D. Wolff, and G. Muller, “The optical properties of biological tissue in the near infrared wavelength range : review and measurements,” in Laser-induced Interstitial Thermotherapy (SPIE, Bellingham WA, Etats-Unis, 1995), pp. 10–44.
  33. J. Chen and X. Intes, “Time-gated perturbation Monte Carlo for whole body functional imaging in small animals,” Opt. Express 17(22), 19566–19579 (2009). [CrossRef] [PubMed]
  34. J. Chen, V. Venugopal, F. Lesage, and X. Intes, “Time-resolved diffuse optical tomography with patterned-light illumination and detection,” Opt. Lett. 35(13), 2121–2123 (2010). [CrossRef] [PubMed]
  35. A. Soubret, J. Ripoll, and V. Ntziachristos, “Accuracy of fluorescent tomography in the presence of heterogeneities: study of the normalized Born ratio,” IEEE Trans. Med. Imaging 24(10), 1377–1386 (2005). [CrossRef] [PubMed]
  36. C. Vinegoni, D. Razansky, J.-L. Figueiredo, M. Nahrendorf, V. Ntziachristos, and R. Weissleder, “Normalized Born ratio for fluorescence optical projection tomography,” Opt. Lett. 34(3), 319–321 (2009). [CrossRef] [PubMed]
  37. B. W. Pogue, T. O. McBride, J. Prewitt, U. L. Osterberg, and K. D. Paulsen, “Spatially variant regularization improves diffuse optical tomography,” Appl. Opt. 38(13), 2950–2961 (1999). [CrossRef] [PubMed]
  38. B. Dogdas, D. Stout, A. F. Chatziioannou, and R. M. Leahy, “Digimouse: a 3D whole body mouse atlas from CT and cryosection data,” Phys. Med. Biol. 52(3), 577–587 (2007). [CrossRef] [PubMed]
  39. “Computational Center for Nanotechnology Innovations (CCNI),” http://www.rpi.edu/research/ccni/ .
  40. V. Venugopal, J. Chen, and X. Intes, “Development of an optical imaging platform for functional imaging of small animals using wide-field excitation,” Biomed. Opt. Express 1(1), 143–156 (2010). [CrossRef] [PubMed]
  41. C. D’Andrea, D. Comelli, A. Pifferi, A. Torricelli, G. Valentini, and R. Cubeddu, “Time-resolved optical imaging through turbid media using a fast data acquisition system based on a gated CCD camera,” J. Phys. D Appl. Phys. 36(14), 1675–1681 (2003). [CrossRef]
  42. X. Cheng and D. Boas, “Systematic diffuse optical image errors resulting from uncertainty in the background optical properties,” Opt. Express 4(8), 299–307 (1999). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited