|
1.IntroductionNovel contrast agents provide outstanding opportunities to visualize and analyze specific biological targets, signal pathways, and therapeutic responses.1, 2 The nanophosphors, or other similar materials, emit near-infrared (NIR) luminescence light upon x-ray excitation.3, 4 As optical probes, these nanophosphors potentially allow sensitive and specific high-resolution imaging in vivo.5 Using the nanophosphors, x-ray luminescence computed tomography (XLCT) was recently proposed as a new molecular imaging modality.6 The experiments demonstrated that XLCT could image the cross-sectional distribution of nanophosphors.7 XLCT has several advantages compared to other optical molecular imaging modalities. First, the anatomy and the nanophosphors can be imaged in one scan. Second, the use of x-ray excitation eliminates the autofluorescence in optical fluorescence imaging. Third, the straight line propagation of x-rays in a biological object means a localized and deep probing capability, promising to significantly improve the spatial resolution.6 Through selective excitation with a pencil x-ray beam, XLCT can perform in vivo tomographic imaging on a region of interest (ROI) to reduce radiation dose and shorten experimental time.8, 9, 10 This selective excitation mechanism is similar to x-ray fluorescence computed tomography (XFCT), which is used to map elements inside samples.11 Recently, XFCT is applied for gold nanoparticle imaging at relatively low concentrations using polychromatic diagnostic x-rays.12 However, a primary issue with XLCT is related to the x-ray scattering effect. For example, x-ray Compton scattering is the major component in water for energy >40 keV.13 The XLCT reconstruction will be significantly compromised if the luminescence signal induced by these scattered x-ray photons is not effectively compensated. In this paper, we propose a systematic approach to address the x-ray scattering effects in the XLCT reconstruction. Furthermore, x-ray excitations are preferably applied with different spectra, and NIR signals are measured in a spectrally resolving fashion. Finally, a linear relationship is established between the nanophosphor distribution and the measured NIR data using the finite element method and inverted using the compressive sensing (CS) technique. This methodology can be adapted for x-ray fluorescence computed tomography (XFCT) as well. Hence, our approach is for spectrally-resolving and scattering-compensated x-ray luminescence/fluorescence computed tomography (SXLCT or SXFCT) for significantly better imaging performance. In Sec. 2, we will first formulate the scattering-compensated XLCT, and then extend it for multi-spectral scattering-compensated XLCT. In Sec. 3, numerical results are reported to demonstrate the feasibility and merits of the proposed approach in each of the aforementioned two cases. In Sec. 4, relevant issues are discussed, and finally the paper is concluded. 2.Imaging ModelNanophosphors can be functionalized, introduced into a biological object such as a mouse, and bound to specific cells or markers in the tissue.2 When an x-ray beam irradiates the tissue, the nanophosphors along the x-ray path will be excited to emit NIR light, which can then be captured by an optical sensor. Our SXLCT algorithm is proposed to accurately localize and reliably quantify the nanophosphor distribution in the tissue from the surface NIR measurements and prior knowledge. The algorithm is based on the photon propagation models in the NIR and x-ray energy ranges, respectively. 2.1.X-ray TransportCurrently, one often assumes that the x-ray beam is limited in the primary path as in Ref. 6. However, a significant part of the x-ray energy is scattered away from the primary path, and the scattered x-ray photons would induce nanophosphor luminescence outside the primary x-ray path. Hence, the x-ray scattering effect should be included in the x-ray photon transport model. X-ray scattering can be incoherent (Compton scattering) or coherent (Rayleigh scattering).13, 14 While an x-ray beam propagates through a tissue along a straight line, the x-ray intensity distribution I(r) along the primary path can be computed by the Lambert-Beer formula Eq. 1[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} I\left({\bf r} \right) = I_0 \exp \left[ { - \int\nolimits_0^1 {\mu _t \left[ {\left({1 - \tau } \right){\kern 1pt} {\bf r}_0 + \tau {\kern 1pt} {\bf r}} \right]} d\tau } \right], \end{equation}\end{document}Eq. 2[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} {\rm X}\left({\bf r} \right) &=& \int\nolimits_L {I\left({{\bf r}_{\bf s} } \right)\eta \left({{\bf r}_s } \right)\frac{{d\sigma _{scatter} }}{{d\Omega }}\left(\beta \right)}\nonumber\\ && \times \exp \left[ { - \int\nolimits_0^1 {\mu _t \left[ {\left({1 - \tau } \right){\kern 1pt} {\bf r}_s + \tau {\kern 1pt} {\bf r}} \right]} d\tau } \right]\,d{\bf r}_s, \end{eqnarray}\end{document}2.2.Near-infrared Light TransportNIR light is subject to both scattering and absorption. The biological soft tissue has highly scattering and weakly absorbing properties in the NIR spectral region. In this case, the diffusion approximation (DA) model usually offers an accurate description for NIR light transport17 Eq. 3[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \hspace*{-1pt}- \nabla \cdot \left[ {D\left({\bf r} \right)\nabla \Phi \left({\bf r} \right)} \right] + \mu _a \left({\bf r} \right)\Phi \left({\bf r} \right) = S\left({\bf r} \right),\quad{\bf r} \in \Omega, \end{equation}\end{document}Eq. 4[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \Phi \left({\bf r} \right) + 2\alpha D\left({\bf r} \right)\left[ {\nu \cdot \nabla \Phi \left({\bf r} \right)} \right] = 0,\,\,\,\,\,\,\,\,\,{\bf r} \in \partial \Omega,\end{equation}\end{document}Eq. 5[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} m\left({\bf r} \right) = - D\left({\bf r} \right)\left({{\bf \nu } \cdot \nabla \Phi \left({\bf r} \right)} \right)\quad{\bf r} \in \partial \Omega. \end{equation}\end{document}2.3.DiscretizationEquations 3, 4 can be discretized into a matrix equation linking the nanophosphor distribution ρ and the NIR photon fluence rate Φ(r) at a node r using the finite element analysis,18 Eq. 7[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\bf A} \cdot \,\Phi = {\bf F} \cdot \,{\rm \rho,} \end{equation}\end{document}Eq. 8[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} a_{ij} &=& \int\nolimits_\Omega {D\left({\bf r} \right)} \nabla \varphi _i \left({\bf r} \right) \cdot \nabla \varphi _j \left({\bf r} \right){\rm d}{\bf r} + \int\nolimits_\Omega {\mu _a \left({\bf r} \right)} \varphi _i \left({\bf r} \right)\varphi _j \left({\bf r} \right){\rm d}{\bf r}\nonumber\\ && + \int\nolimits_{\partial \Omega } {\varphi _i \left({\bf r} \right)} {{\varphi _j \left({\bf r} \right)}/ {{\rm 2}\alpha }}{\rm d}{\bf r}, \end{eqnarray}\end{document}Eq. 9[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} f_{ij} = \varepsilon \int\nolimits_\Omega {{\rm X}\left({\bf r} \right)\varphi _i \left({\bf r} \right)} \varphi _j \left({\bf r} \right){\rm d}{\bf r}, \end{equation}\end{document}Eq. 10[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \Phi = ({{\bf A}\,^{ - 1} {\bf F}}) \cdot {\rm \rho }{\rm.} \end{equation}\end{document}Eq. 11[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \sum\limits_{i = 1}^n {\Phi \left({{\bf r}_{m_i } } \right)} = \left({\sum\limits_{m_i } {{\bf A}^{ - {\bf 1}} {\bf F}} } \right) \cdot {\rm \rho,} \end{equation}\end{document}Eq. 12[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} Q_e^{ - 1} \;{\rm M} = {\bf G} \cdot {\rm \rho,} \end{equation}\end{document}2.4.ReconstructionIn the CS framework, one can reconstruct a sparse image from far less samples than what the Nyquist sampling theorem demands.19 Based on the biological characteristics, nanoparticles preferentially seek the specific cells in the tissue, forming a sparse or smooth distribution. As a result, applying an appropriate sparsifying transform T (such as the total variation or wavelet transform), the nanophosphor concentration image can be converted into a sparse image that has far fewer significant pixel values than the total number of pixels. Using a CS technique,19 we can reconstruct a nanoparticle density distribution by solving the following optimization problem Eq. 13[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} \begin{array}{l} {\it minimum} {\rm }\left\| {{\bf T} \cdot {\rm \rho }} \right\|_1 \\ {\it{subject\,to}}\left\| {{\bf G} \cdot {\rm \rho - }Q_e^{ - 1} {\rm M}} \right\|_2 \le \tau,\\ \qquad\qquad\;{\rm \rho } \ge {\rm 0} \end{array} \end{eqnarray}\end{document}The above formulation is valid for both single- and multi-spectral x-ray excitation patterns. The induced phosphor luminescent spectra may depend on both the x-ray dose and energy in a nonlinear fashion, which may or may not be linearly approximated in practical applications. Nevertheless, these details can be readily worked out in a particular scenario if the need arises. 3.Numerical SimulationWe performed representative numerical tests to evaluate the proposed SXLCT methodology with a numerical mouse phantom and synthetic measurement data. As shown in Fig. 3, the numerical mouse phantom was established from the CT slices of a mouse using Amira (Amira 4.0, Mercury Computer Systems, Inc., Chelmsford, Massachusetts). The phantom was discretized into 203,690 tetrahedral elements with 58,246 nodes. The emission wavelength of the nanophosphor under x-ray excitation was set to 802 nm, which was based on the emission characteristics of Gd2O2S doped with thulium.3, 4 Based on Ref. 21, the light yield ε was assumed to be about 0.15 cm3/mg for 50 keV. Appropriate optical parameters were accordingly assigned to the mouse model. The reduced scattering coefficient [TeX:] $\mu^{\prime}_s \left(\lambda \right)$ relies on the wavelength λ(nm) and is approximated by an empirical function Eq. 14[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \mu '_s \left(\lambda \right) = 10a \cdot \lambda ^{ - b}, \end{equation}\end{document}Eq. 15[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \mu _a \left(\lambda \right) = S_B [ {x\mu _{aHb} \left(\lambda \right) + \left({1 - x} \right)\mu _{aHbO_2 } \left(\lambda \right)}] + S_W \mu _{aW} \left(\lambda \right), \end{equation}\end{document}3.1.Spatial ResolutionThree spherical sub-regions of 0.5 mm diameter were centered at (14.3, 40, 21.8), (12.3, 40, 20.5), and (12.3, 40, 23.0), respectively, and filled with nanophosphor concentrations of 1.5 μg/mL, to test the spatial resolution. An x-ray source was operated at 50 keV and 30 mA to perpendicularly irradiate the phantom surface at the longitudinal position 40 mm. The x-rays were collimated into a pencil beam of 0.2 mm width. The photon detector quantum efficiency Q E was set to 1. The x-ray scattering distribution was calculated based on the x-ray differential scattering cross section and the electron density distribution in the phantom (Fig. 1). The NIR light emitted from the excited nanophosphors underwent scattering and absorption in the phantom, and the intensity of the NIR on the surface of the phantom was simulated according to the DA model. Poisson noise was added to the synthetic measurement data. The image acquisition procedure was repeated 100 times when the x-ray pencil beam was translated in a 0.2 mm increment for a given view. The parallel-beam imaging geometry was rotated 20 times to cover a 360 deg range evenly. Based on the known geometrical parameters, optical properties, and x-ray energy, we formed Eq. 12 using the finite element method. The linear equation system was quite sparse and efficiently manipulated. In the CS framework, we employed a fast iterative shrinkage-thresholding algorithm (FISTA) (Ref. 20) for the reconstruction of the nanophosphor distribution. The FISTA optimization algorithm demonstrated an excellent performance in terms of convergence and stability. The reconstructed results are in close agreement with the truth, and the average relative error of the reconstructed nanophosphor density was less than 10%, which was defined as [TeX:] $\bar e = \frac{1}{N}\sum\nolimits_{k = 1}^N {\frac{{\left| {\Phi _k^S - \Phi _k^T } \right|}}{{\Phi _k^T }}}$ , where [TeX:] $\Phi _k^T$ and [TeX:] $\Phi _k^S$ are the true and computed nanoparticles concentrations, respectively. Figures 4a and 4c present the comparison at the cross section through the nanophosphor distribution between the reconstructed and true nanophosphor density distributions. Then, we repeated the numerical experiment without compensating for x-ray scattering, while the synthetic measurement data on the surface of the phantom were kept intact; that is, which were generated by both the x-ray primary beam and associated x-ray scattering. The results indicated that the reconstructed average relative error was more than 40% compared to the true nanophosphor density distribution, as shown in Figs. 4b and 4c. In this case, the light source excited by x-ray scattering was mistakenly treated as the luminescence emission due to the x-ray primary beam. Hence, the resultant reconstructed phosphor density distribution was much higher than the truth. 3.2.Particle Density ResolutionThe detection sensitivity for very low concentrations of probes is significant in molecular and cellular imaging. To evaluate the density resolution of the proposed methodology, the same three spherical sub-regions were assumed in the mouse phantom as above described and filled with the nanophosphor concentrations of 0.1 μg/mL, 0.3 μg/mL, and 1.2 μg/mL, respectively. The x-ray source of 50 keV and 30 mA was collimated into a pencil beam again to irradiate the phantom at the longitudinal position 40 mm, translated 100 times in a 0.2 mm increment to generate NIR light for a view angle, and rotated 20 times to cover 360 deg for NIR data acquisition. The reconstruction was similarly conducted as it was previously. Figure 5a shows the reconstructed nanophosphor density distributions. The results are also in excellent agreement with the truth in Fig. 5c. In particular, the results show that our proposed method can identify the phosphor nanoparticles in the concentration level as small as 0.1 μg/m in the scenario of various nanophosphor densities up to 1.2 μg/mL. For comparison, we repeated the numerical experiment assuming no x-ray scattering. The measurable photon fluence rate on the surface of the phantom was generated from the luminescence source of nanophosphors excited from the x-ray primary and scattering photons. The results show that the reconstructed nanophosphors density distribution had an average relative error more than 40% relative to the true nanophosphors density distribution, as shown in Figs. 5b and 5c. 3.3.Spectrally Resolved ReconstructionThe phosphors can be adapted to emit light signals of different spectra depending on the doping material; for example, Terbium emits green (545 nm) light, and Thulium emits near-infrared (802 nm) light under the 100 keV x-ray irradiation. Different nanophosphors have a variety of fluorescence emission peaks and serve as distinct types of molecular probes for molecular and cellular imaging. In general, the maximum peak emission of a given type of nanophosphors is formed at its x-ray absorption K-edge. Luminescence spectra of different probes can be separated using optical filters and advanced algorithms. For multi-probe optical molecular imaging, we should measure multi-spectral NIR signals on the mouse body surface. In our numerical simulation, two spherical regions A and B of 0.5 mm diameter were centered at (13.5, 40, 21) and (12.1, 40, 23.0), respectively. The A region was filled with nanophosphors of type A with a concentration of 0.7 μg/mL and nanophosphors of type B with a concentration of 0.3 μg/mL. The B region was filled with nanophosphors of type A with a concentration of 0.3 μg/mL and nanophosphors of type B with a concentration of 0.7 μg/mL. Under the 50 keV x-ray excitation, nanophosphors of type A emit NIR light of 700 nm. Under the 80 keV x-ray excitation, nanophosphors of type A emit NIR light of 700 nm and nanophosphors of type B emit NIR light of 800 nm. We first excited the numerical phantom with 50 keV x-rays to acquire the NIR signal of 700 nm, and then excited the phantom with 80 keV x-rays to acquire the NIR signals of 700 and 800 nm. In the data acquisition process, the x-ray pencil beam was translated 100 times in a 0.2 mm increment for a given view. The imaging geometry was rotated 20 times to cover a 360 deg range evenly. According to Eq. 12, the linear system equation can be formed to link NIR signals from multiple nanophosphor concentrations in a spectral band [ων, ων + 1] to measurable photon fluence rate data in the same spectral band Eq. 16[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} Q_{e,v}^{ - 1} {\rm M}_\nu = {\bf G}_\nu \cdot {\rm \rho }_\nu,\quad\nu = 1,2, \ldots,\kappa, \end{equation}\end{document}4.Discussions and ConclusionLet us first compare XFCT and XLCT. On one hand, XFCT and XLCT are similar because both of them can be considered as emission-type tomography with data stimulated with external x-rays. On the other hand, there are clear differences between them. Fluorescence refers to photon emission after photon absorption. Luminescence is indirect, such as in the case of nanophosphors. X-rays or gamma rays excite phosphors and other materials via ionization and create luminescence upon recombination of electrons and ions. Most importantly, x-ray fluorescence and luminescence signals are generated with various efficiencies in distinct spectral regions and have different absorption and scattering properties. The former is less strongly absorbed and scattered than the latter. While the x-ray scattering effect is significant, neither XFCT nor XLCT compensates for such an effect in the image reconstruction process. In the XFCT literature, the scattering effect is treated as noise.25 In the XLCT work, the Stanford group that developed XLCT has implicitly assumed that the stimulating x-rays traveled strictly in straight lines.6 To address this source of inaccuracy, the x-ray scattering compensation as we have showcased in this paper must be a key step for either XLCT or XFCT, and especially for the former. With the rapid development of the x-ray source and detector technology, it is anticipated that in the not-too-distant future x-rays can be cost-effectively delivered wherever (distributed sources), whenever (programmable multiplexing), and whatever (energy modulation) we desire.26 Therefore, it is valuable and practical to consider various x-ray excitation patterns more complicated than what is currently used for XFCT and XLCT, which is in a pencil beam of a single x-ray spectrum. As suggested above, the spectrally resolving excitation and measurement can be systematically done for optimal image quality, which is especially important for multi-element and/or multi-probe x-ray-excitation-based tomography. Based on the previous work on XFCT and XLCT,6, 12 here we have proposed a SXLCT/SXFCT approach to quantify a nanophosphor/element distribution within a biological object. In this initial study, the SXLCT method has been developed to reconstruct a nanophosphor distribution in a small animal model. Upon x-ray excitation, nanophosphors emit NIR light around 700 nm. Within this wavelength range, the NIR light scattering dominates over absorption in the tissue. The photon propagation can be modeled as a diffusive process, and the DA model is quite accurate. The x-ray scattering effect is also important and has been taken into account. Using the finite element method, a linear relationship has been formulated to link the nanophosphor density distribution and the measured NIR data. The CS technique has been applied to solve this inverse problem. The numerical simulation results have demonstrated that the proposed approach can accurately reconstruct the nanophosphor distribution. There are several topics deserving further research. Since this reconstruction method relies on optical parameters (absorption and scattering coefficients), optical diffuse tomography, photoacoustic tomography, or related techniques should be used to determine such prior knowledge. Also, x-ray scattering tomography techniques should be further developed assuming the availability of spectral x-ray detectors such as MedPix3. Furthermore, XLCT currently works in a pencil beam scanning mode and takes long sampling time. Fan-beam or other scanning modes can improve the imaging efficiency and is under development in our laboratory. AcknowledgmentsThis work was partially supported by the National Institutes of Health Grant Nos. NIH/NHLBI HL098912 and NIH/NCI CA135151. ReferencesR. F. Minchin and
D. J. Martin,
“Minireview: nanoparticles for molecular imaging—an overview,”
Endocrinology, 151 474
–481
(2010). https://doi.org/10.1210/en.2009-1012 Google Scholar
X. Qian, X. H. Peng, D. O. Ansari, Q. Yin-Goen, G. Z. Chen, D. M. Shin, L. Yang, A. N. Young, M. D. Wang, and
S. Nie,
“In vivo tumor targeting and spectroscopic detection with surface-enhanced Raman nanoparticle tags,”
Nat. Biotechnol., 26 83
–90
(2008). https://doi.org/10.1038/nbt1377 Google Scholar
Y. Tian, W.-H. Cao, X.-X. Luo, and
Y. Fu,
“Preparation and luminescence property of Gd2O2S:Tb X-ray nano-phosphors using the complex precipitation method,”
J. Alloys. Compd., 433 313
–317
(2007). https://doi.org/10.1016/j.jallcom.2006.06.075 Google Scholar
M. Xing, W. Cao, T. Pang, X. Ling, and
N. Chen,
“Preparation and characterization of monodisperse spherical particles of X-ray nanophosphors based on Gd2O2S:Tb,”
Chin. Sci. Bull., 54 2982
–2986
(2009). https://doi.org/10.1007/s11434-009-0492-9 Google Scholar
C. Carpenter, G. Pratx, C. Sun, and
L. Xing,
“Hybrid x-ray/optical luminescence imaging: characterization of experimental conditions,”
Med. Phys., 37 4011
–4018
(2010). https://doi.org/10.1118/1.3457332 Google Scholar
G. Pratx, C. Carpenter, C. Sun, and
L. Xing,
“X-ray luminescence computed tomography via selective excitation: a feasibility study,”
IEEE Trans. Med. Imaging, 29 1992
–1999
(2010). https://doi.org/10.1109/TMI.2010.2055883 Google Scholar
G. Pratx, C. M. Carpenter, C. Sun, R. P. Rao, and
L. Xing,
“Tomographic molecular imaging of x-ray-excitable nanoparticles,”
Opt. Lett., 35 3345
–3347
(2010). https://doi.org/10.1364/OL.35.003345 Google Scholar
Y. Ye, H. Yu, Y. Wei, and
G. Wang,
“A general local reconstruction approach on a truncated Hilbert transform,”
Int. J. Biomed. Imaging, 2007 63634
(2007). https://doi.org/10.1155/2007/63634 Google Scholar
H. Yu and
G. Wang,
“Compressive sensing based interior tomography,”
Phys. Med. Biol., 54 2791
–2805
(2009). https://doi.org/10.1088/0031-9155/54/9/014 Google Scholar
J. Yang, H. Yu, M. Jiang, and
G. Wang,
“High-order total variation minimization for interior tomography,”
Inverse Problems, 26 035013
(2010). https://doi.org/10.1088/0266-5611/26/3/035013 Google Scholar
P. Boisseau,
“Determination of three dimensional trace element distributions by use of monochromatic x-ray microbeams,”
Massachusetts Institute of Technology,
(1986). Google Scholar
S.-K. Cheong, B. L. Jones, A. K. Siddiqi, F. Liu, N. Manohar, and
S. H. Cho,
“X-ray fluorescence computed tomography (XFCT) imaging of gold nanoparticle-loaded objects using 110 kVp x-rays,”
Phys. Med. Biol., 55 647
–662
(2010). https://doi.org/10.1088/0031-9155/55/3/007 Google Scholar
J. Als-Nielsen and
D. McMorrow, Elements of Modern X-ray Physics, Wiley, New York
(2001). Google Scholar
J. H. Hubbell, W. J. Veigele, E. A. Briggs, R. T. Brown, D. T. Cromer, and
R. J. Howerton,
“Atomic form factors, incoherent scattering functions, and photon scattering cross sections,”
J. Phys. Chem. Ref. Data, 4 471
–538
(1975). https://doi.org/10.1063/1.555523 Google Scholar
J. Baro, J. Sempau, J. M. Fernandez-Varea, and
F. Salvat,
“PENELOPE: an algorithm for Monte Carlo simulation of the penetration and energy loss of electrons and positrons in matter,”
Nucl. Instrum. Methods Phys. Res. B, 100 31
–46
(1995). https://doi.org/10.1016/0168-583X(95)00349-5 Google Scholar
G. Wang, W. Cong, H. Shen, and
Y. Zou,
“Varying collimation for dark-field extraction,”
Int. J. Biomed. Imaging, 2009 847537
(2009). https://doi.org/10.1155/2009/847537 Google Scholar
J. Welch and
M. J. C. van Gemert, Optical and Thermal Response of Laser-irradiated Tissue, Plenum, New York
(1995). Google Scholar
W. Cong, G. Wang, D. Kumar, Y. Liu, M. Jiang, L. Wang, E. Hoffman, G. McLennan, P. McCray, J. Zabner, and
A. Cong,
“Practical reconstruction method for bioluminescence tomography,”
Opt. Express, 13 6756
–6771
(2005). https://doi.org/10.1364/OPEX.13.006756 Google Scholar
E. J. Candes, J. Romberg, and
T. Tao,
“Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,”
IEEE Trans. Inf. Theory, 52 489
–509
(2006). https://doi.org/10.1109/TIT.2005.862083 Google Scholar
A. Beck and
M. Teboulle,
“A fast iterative shrinkage-thresholding algorithm for linear inverse problems,”
SIAM J. Imaging Sci., 2 183
–202
(2009). https://doi.org/10.1137/080716542 Google Scholar
C. W. E. van Eijk,
“Inorganic scintillators in medical imaging,”
Phys. Med. Biol., 47 R85
–R106
(2002). https://doi.org/10.1088/0031-9155/47/8/201 Google Scholar
G. Alexandrakis, F. R. Rannou, and
A. F. Chatziioannou,
“Tomographic bioluminescence imaging by use of a combined optical-PET (OPET) system: a computer simulation feasibility study,”
Phys. Med. Biol., 50 4225
–4241
(2005). https://doi.org/10.1088/0031-9155/50/17/021 Google Scholar
S. A. Prahl,
“Optical properties spectra,”
(2001) http://omlc.ogi.edu/spectra/index.html Google Scholar
R. Lansford, G. Bearman, and
S. E. Fraser,
“Resolution of multiple green fluorescent protein color variants and dyes using two-photon microscopy and imaging spectroscopy,”
J. Biomed. Opt., 6 311
–318
(2001). https://doi.org/10.1117/1.1383780 Google Scholar
E. X. Miqueles and
A. R. De Pierr,
“Exact analytic reconstruction in x-ray fluorescence CT and approximated versions,”
Phys. Med. Biol., 55 1007
–1024
(2010). https://doi.org/10.1088/0031-9155/55/4/007 Google Scholar
J. Zhang, G. Yang, Y. Cheng, B. Gao, Q. Qiu, Y. Z. Lee, J. P. Lu, and
O. Zhou,
“Stationary scanning x-ray source based on carbon nanotube field emitters,”
Appl. Phys. Lett., 86 184104
(2005). https://doi.org/10.1063/1.1923750 Google Scholar
|