Research Papers: Imaging

Spectrally resolving and scattering-compensated x-ray luminescence/fluorescence computed tomography

[+] Author Affiliations
Wenxiang Cong, Haiou Shen, Ge Wang

Virginia Polytechnic Institute and State University, School of Biomedical Engineering and Sciences, Biomedical Imaging Division, Blacksburg, Virginia 24061

J. Biomed. Opt. 16(6), 066014 (June 13, 2011). doi:10.1117/1.3592499
History: Received December 07, 2010; Revised April 15, 2011; Accepted April 28, 2011; Published June 13, 2011; Online June 13, 2011
Text Size: A A A

Open Access Open Access

* Address all correspondence to: Ge Wang, Virginia Tech, School of Biomedical Engineering and Sciences, Stanger Street, Blacksburg, Virginia 24061; Tel: 540-231-0493; Fax: 540-231-0970; E-mail: wangg@vt.edu.

The nanophosphors, or other similar materials, emit near-infrared (NIR) light upon x-ray excitation. They were designed as optical probes for in vivo visualization and analysis of molecular and cellular targets, pathways, and responses. Based on the previous work on x-ray fluorescence computed tomography (XFCT) and x-ray luminescence computed tomography (XLCT), here we propose a spectrally-resolving and scattering-compensated x-ray luminescence/fluorescence computed tomography (SXLCT or SXFCT) approach to quantify a spatial distribution of nanophosphors (other similar materials or chemical elements) within a biological object. In this paper, the x-ray scattering is taken into account in the reconstruction algorithm. The NIR scattering is described in the diffusion approximation model. Then, x-ray excitations are 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 measured NIR data using the finite element method and inverted using the compressive sensing technique. The numerical simulation results demonstrate the feasibility and merits of the proposed approach.

Figures in this Article

Novel contrast agents provide outstanding opportunities to visualize and analyze specific biological targets, signal pathways, and therapeutic responses.12 The nanophosphors, or other similar materials, emit near-infrared (NIR) luminescence light upon x-ray excitation.34 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.810 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.

Nanophosphors 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.

X-ray Transport

Currently, 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).1314 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 Display Formula

1Ir=I0exp01μt1τr0+τrdτ,
where r0 is a source position, I0 is the x-ray source intensity, and μt is the attenuation coefficient that can be computed from measured x-ray transmission data in the XLCT system using an attenuation-based computed tomography (CT) technique. From the definition of the differential cross section of x-ray photon scattering, the x-ray single scattering intensity at a position r from the primary x-ray beam can be described as follows: Display Formula
2Xr=LIrsηrsdσscatterdΩβ×exp01μt1τrs+τrdτdrs,
where L is a line segment as the intersection between the x-ray primary beam and the object, β is the x-ray photon scattering angle computed by β=cos1rrs·νrrs, ν is the direction of the x-ray primary beam, η(rs) is the electron density at the position rs, and the x-ray differential scattering cross-section dσscatterdΩβ can be evaluated by the Klein-Nishina formula and the differential cross section of Rayleigh scattering.14 Quantitatively, we estimated the x-ray scattering effect. A cubic phantom of 40 mm side length was filled with pure water. An x-ray pencil beam of diameter 0.2 mm at 50 keV perpendicularly irradiated the phantom at the center of one surface. The attenuation coefficient and the electron density of water are 0.214 (cm−1) and 3.346 × 1023/cm3, respectively. The x-ray scattering distribution is shown in Fig. 1. The x-ray scattering effect can be also modeled using a Monte Carlo simulation method,15 an approximation Monte Carlo technique,16 a radiative transport equation solver, or its variant to take multiple-scattering photons into account.

Graphic Jump LocationF1 :

Significance of x-ray scattering. (a) X-ray pencil beam passes through a water phantom and (b) 50% of scattered photons are distributed in the region delimited by the blue and green curves.

Near-infrared Light Transport

NIR 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 transport17Display Formula

3·DrΦr+μarΦr=Sr,rΩ,
where r is the position vector, Φ(r) is the NIR photon fluence rate, S(r) is the NIR source, μa is the absorption coefficient, D is the diffusion coefficient defined by D=[3(μa+μs)]1, μs is the reduced scattering coefficient, and Ω⊂R3 is a region of interest. If no photon travels across the boundary ∂Ω into the tissue domain Ω, the DA is constrained by the Robin boundary condition Display Formula
4Φr+2αDrν·Φr=0,rΩ,
where ν is the outward unit normal vector on ∂Ω, and α is the boundary mismatch factor. The boundary mismatch factor between the tissue with a refractive index n and air can be approximated by α = (1 + γ)/(1 − γ) with γ = −1.4399n−2 + 0.7099n−1 + 0.6681 + 0.0636n.17 The measurable exiting photon flux on the surface of the object can be expressed as Display Formula
5mr=Drν·ΦrrΩ.
The intensity of the emitted NIR light depends on the nanophosphor density ρ(r), the x-ray intensity X(r), and the light yield ε, which can be defined as the quantum yield per unit nanophosphor concentration. The light yield can be experimentally determined. Although that dependency is nonlinear, in this feasibility study we assume that the intensity of NIR light is linearly proportional to both the x-ray dose and nanophosphor concentration,5Display Formula
6Sr=ɛXrρr.
Note that the luminescence phenomenon may also depend on the energy of the incoming x-ray beam but is temporarily ignored for simplicity.

Discretization

Equations 34 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,18Display Formula

7A·Φ=F·ρ,
where the component of the matrix A is Display Formula
8aij=ΩDrϕir·ϕjrdr+Ωμarϕirϕjrdr+Ωϕirϕjr/2αdr,
and the component of the matrix F is Display Formula
9fij=ɛΩXrϕirϕjrdr,
where φi(i = 1, 2, …) are the element shape functions. Since the matrix A in Eq. 7 is positive definite, we have Display Formula
10Φ=(A1F)·ρ.
Because the number of outgoing NIR photons at a surface node is insufficient in general, we prefer integrating detected outgoing photons for a higher signal-to-noise ratio. Hence, from Eq. 10 we obtain the following linear equation for an x-ray pencil beam excitation, Display Formula
11i=1nΦrmi=miA1F·ρ,
where mi(i = 1, 2, …, n) is the index for measurable surface nodes. Note that in the imaging process, the x-ray beam can be moved or modulated multiple times to irradiate the object for a given view, and an array of NIR data will be acquired around the object as shown in Fig. 2. Based on Eq. 11, we can form the linear equation system as follows: Display Formula
12Qe1M=G·ρ,
where M is the vector from the photon fluence rate data measured with multiple selective excitations, Qe is the photon-detector quantum efficiency, and G is a weighting matrix assembled with a row vector miA1F for each x-ray excitation. Equation 12 describes the linear relationship between the nanophosphor distribution and the NIR measurement.

Graphic Jump LocationF2 :

SXLCT system setup that uses x-rays under spectral modulation.

Reconstruction

In 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 Display Formula

13minimumT·ρ1subjecttoG·ρQe1M2τ,ρ0
where τ represents the data noise level. The presence of the l1 term is used to induce the solution sparsity. Equation 13 can be efficiently solved using a contemporary numerical method such as the Bregman iteration.20

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.

We 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.34 Based on Ref. 21, the light yield ε was assumed to be about 0.15 cm3/mg for 50 keV.

Graphic Jump LocationF3 :

Mouse phantom represented in a finite element mesh.

Appropriate optical parameters were accordingly assigned to the mouse model. The reduced scattering coefficient

μsλ
relies on the wavelength λ(nm) and is approximated by an empirical function Display Formula
14μsλ=10a·λb,
where a and b are the constants depending on the tissue type. The parameters a and b for organ-specific values can be found.22 The tissue absorption depends on the local oxy-hemoglobin (HbO2), deoxy-hemoglobin (Hb), and water (W) concentrations in the respective organs. The spectral absorption coefficient μa(λ) can be approximated as the weighted sum of the three absorption coefficients μaHbO2λ, μaHb(λ), and μaW(λ), which were calculated from the corresponding absorbance spectra reported in Ref. 22Display Formula
15μaλ=SB[xμaHbλ+1xμaHbO2λ]+SWμaWλ,
where x = HbO2/(HbO2 + Hb) is the ratio between oxy-hemoglobin and the total hemoglobin concentration, SB and SW are scaling factors specific to the respective organs.23 The relative electron density of 0.6 was assigned to the phantom, which corresponds to that of soft tissues.

Spatial Resolution

Three 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 QE 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

e¯=1Nk=1NΦkSΦkTΦkT
, where ΦkT and ΦkS are the true and computed nanoparticles concentrations, respectively. Figures 4 present the comparison at the cross section through the nanophosphor distribution between the reconstructed and true nanophosphor density distributions.

Graphic Jump LocationF4 :

Spatial resolution comparison between the reconstructed and true nanophosphor density distributions. (a) Reconstructed nanophosphor density distribution with the proposed SXLCT approach, yielding an average relative error of 10%; (b) reconstructed nanophosphor density distributions assuming no x-ray scattering in the model, yielding an average relative error over 40%; and (c) true nanophosphor density distribution.

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. 4. 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.

Particle Density Resolution

The 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 5 shows the reconstructed nanophosphor density distributions. The results are also in excellent agreement with the truth in Fig. 5. 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.

Graphic Jump LocationF5 :

Density resolution comparison between the reconstructed and true nanophosphor density distributions. (a) Reconstructed nanophosphor density distribution with the proposed SXLCT approach, producing an average relative error of 12%; (b) reconstructed nanophosphor density distribution assuming no x-ray scattering, resulting in an average relative error over 40%; and (c) true nanophosphor density distribution.

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. 5.

Spectrally Resolved Reconstruction

The 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 Display Formula

16Qe,v1Mν=Gν·ρν,ν=1,2,...,κ,
where  κ is the number of spectral bands,  Mν consists of measurable photon fluence rate values at the nodes on the object surface in the same spectral band [ων, ων + 1], and  ρν is the NIR signals from multiple nanophosphor concentrations in a spectral band [ων, ων + 1], that is, ρν=iwipi,ν,ν=1,2,...,κ, where pi, ν is the ith probe concentration in a spectral band [ων, ων + 1] and wi is its weighting factor. Equation 16 can be numerically solved in the compressive sensing framework. Then, we can obtain the individual distributions of a different molecular probes from  ρν data using an unmixing technique.24 Figure 6 shows representative reconstructed images for distributions of probe types A and B, respectively. The reconstructed results are in excellent agreement with the truth.

Graphic Jump LocationF6 :

Spectral resolution comparison between the reconstructed and true nanophosphor density distributions. (a) Reconstructed distribution for nanophosphors of type A; (b) true distribution for nanophosphors of type A; (c) reconstructed distribution for nanophosphors of type B; and (d) true distribution for nanophosphors of type B.

Let 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.

This work was partially supported by the National Institutes of Health Grant Nos. NIH NIH/NHLBI HL098912 and NIH NIH/NCI CA135151 .

Minchin  R. F., and Martin  D. J., “ Minireview: nanoparticles for molecular imaging—an overview. ,” Endocrinology. 151, , 474–481  ((2010)).
Qian  X., , Peng  X. H., , Ansari  D. O., , Yin-Goen  Q., , Chen  G. Z., , Shin  D. M., , Yang  L., , Young  A. N., , Wang  M. D., , and Nie  S., “ In vivo tumor targeting and spectroscopic detection with surface-enhanced Raman nanoparticle tags. ,” Nat. Biotechnol.. 26, , 83–90  ((2008)).
Tian  Y., , Cao  W.-H., , Luo  X.-X., , and Fu  Y., “ Preparation and luminescence property of Gd2O2S:Tb X-ray nano-phosphors using the complex precipitation method. ,” J. Alloys. Compd.. 433, , 313–317  ((2007)).
Xing  M., , Cao  W., , Pang  T., , Ling  X., , and Chen  N., “ Preparation and characterization of monodisperse spherical particles of X-ray nanophosphors based on Gd2O2S:Tb. ,” Chin. Sci. Bull.. 54, , 2982–2986  ((2009)).
Carpenter  C., , Pratx  G., , Sun  C., , and Xing  L., “ Hybrid x-ray/optical luminescence imaging: characterization of experimental conditions. ,” Med. Phys.. 37, , 4011–4018  ((2010)).
Pratx  G., , Carpenter  C., , Sun  C., , and Xing  L., “ X-ray luminescence computed tomography via selective excitation: a feasibility study. ,” IEEE Trans. Med. Imaging. 29, , 1992–1999  ((2010)).
Pratx  G., , Carpenter  C. M., , Sun  C., , Rao  R. P., , and Xing  L., “ Tomographic molecular imaging of x-ray-excitable nanoparticles. ,” Opt. Lett.. 35, , 3345–3347  ((2010)).
Ye  Y., , Yu  H., , Wei  Y., , and Wang  G., “ A general local reconstruction approach on a truncated Hilbert transform. ,” Int. J. Biomed. Imaging. 2007, , 63634  ((2007)).
Yu  H., and Wang  G., “ Compressive sensing based interior tomography. ,” Phys. Med. Biol.. 54, , 2791–2805  ((2009)).
Yang  J., , Yu  H., , Jiang  M., , and Wang  G., “ High-order total variation minimization for interior tomography. ,” Inverse Problems. 26, , 035013  ((2010)).
Boisseau  P., “ Determination of three dimensional trace element distributions by use of monochromatic x-ray microbeams. ,” PhD Thesis, Massachusetts Institute of Technology ((1986)).
Cheong  S.-K., , Jones  B. L., , Siddiqi  A. K., , Liu  F., , Manohar  N., , and Cho  S. H., “ X-ray fluorescence computed tomography (XFCT) imaging of gold nanoparticle-loaded objects using 110 kVp x-rays. ,” Phys. Med. Biol.. 55, , 647–662  ((2010)).
Als-Nielsen  J., and McMorrow  D.,  Elements of Modern X-ray Physics. ,  Wiley ,  New York (2001).
Hubbell  J. H., , Veigele  W. J., , Briggs  E. A., , Brown  R. T., , Cromer  D. T., , and Howerton  R. J., “ Atomic form factors, incoherent scattering functions, and photon scattering cross sections. ,”J. Phys. Chem. Ref. Data. 4, , 471–538  ((1975)).
Baro  J., , Sempau  J., , Fernandez-Varea  J. M., , and Salvat  F., “ 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)).
Wang  G., , Cong  W., , Shen  H., , and Zou  Y., “ Varying collimation for dark-field extraction. ,” Int. J. Biomed. Imaging. 2009, , 847537  ((2009)).
Welch  J., and van Gemert  M. J. C.,  Optical and Thermal Response of Laser-irradiated Tissue. ,  Plenum ,  New York  ((1995)).
Cong  W., , Wang  G., , Kumar  D., , Liu  Y., , Jiang  M., , Wang  L., , Hoffman  E., , McLennan  G., , McCray  P., , Zabner  J., , and Cong  A., “ Practical reconstruction method for bioluminescence tomography. ,” Opt. Express. 13, , 6756–6771  ((2005)).
Candes  E. J., , Romberg  J., , and Tao  T., “ Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. ,” IEEE Trans. Inf. Theory. 52, , 489–509  ((2006)).
Beck  A., and Teboulle  M., “ A fast iterative shrinkage-thresholding algorithm for linear inverse problems. ,” SIAM J. Imaging Sci.. 2, , 183–202  ((2009)).
van Eijk  C. W. E., “ Inorganic scintillators in medical imaging. ,” Phys. Med. Biol.. 47, , R85–R106  ((2002)).
Alexandrakis  G., , Rannou  F. R., , and Chatziioannou  A. F., “ Tomographic bioluminescence imaging by use of a combined optical-PET (OPET) system: a computer simulation feasibility study. ,” Phys. Med. Biol.. 50, , 4225–4241  ((2005)).
Prahl  S. A., “ Optical properties spectra. ,” Oregon Medical Laser Clinic, http://omlc.ogi.edu/spectra/index.html ((2001)).
Lansford  R., , Bearman  G., , and Fraser  S. E., “ Resolution of multiple green fluorescent protein color variants and dyes using two-photon microscopy and imaging spectroscopy. ,” J. Biomed. Opt.. 6, , 311–318  ((2001)).
Miqueles  E. X., and De Pierr  A. R., “ Exact analytic reconstruction in x-ray fluorescence CT and approximated versions. ,” Phys. Med. Biol.. 55, , 1007–1024  ((2010)).
Zhang  J., , Yang  G., , Cheng  Y., , Gao  B., , Qiu  Q., , Lee  Y. Z., , Lu  J. P., , and Zhou  O., “ Stationary scanning x-ray source based on carbon nanotube field emitters. ,” Appl. Phys. Lett.. 86, , 184104  ((2005)).
© 2011 Society of Photo-Optical Instrumentation Engineers (SPIE)

Citation

Wenxiang Cong ; Haiou Shen and Ge Wang
"Spectrally resolving and scattering-compensated x-ray luminescence/fluorescence computed tomography", J. Biomed. Opt. 16(6), 066014 (June 13, 2011). ; http://dx.doi.org/10.1117/1.3592499


Figures

Graphic Jump LocationF1 :

Significance of x-ray scattering. (a) X-ray pencil beam passes through a water phantom and (b) 50% of scattered photons are distributed in the region delimited by the blue and green curves.

Graphic Jump LocationF2 :

SXLCT system setup that uses x-rays under spectral modulation.

Graphic Jump LocationF3 :

Mouse phantom represented in a finite element mesh.

Graphic Jump LocationF4 :

Spatial resolution comparison between the reconstructed and true nanophosphor density distributions. (a) Reconstructed nanophosphor density distribution with the proposed SXLCT approach, yielding an average relative error of 10%; (b) reconstructed nanophosphor density distributions assuming no x-ray scattering in the model, yielding an average relative error over 40%; and (c) true nanophosphor density distribution.

Graphic Jump LocationF5 :

Density resolution comparison between the reconstructed and true nanophosphor density distributions. (a) Reconstructed nanophosphor density distribution with the proposed SXLCT approach, producing an average relative error of 12%; (b) reconstructed nanophosphor density distribution assuming no x-ray scattering, resulting in an average relative error over 40%; and (c) true nanophosphor density distribution.

Graphic Jump LocationF6 :

Spectral resolution comparison between the reconstructed and true nanophosphor density distributions. (a) Reconstructed distribution for nanophosphors of type A; (b) true distribution for nanophosphors of type A; (c) reconstructed distribution for nanophosphors of type B; and (d) true distribution for nanophosphors of type B.

Tables

References

Minchin  R. F., and Martin  D. J., “ Minireview: nanoparticles for molecular imaging—an overview. ,” Endocrinology. 151, , 474–481  ((2010)).
Qian  X., , Peng  X. H., , Ansari  D. O., , Yin-Goen  Q., , Chen  G. Z., , Shin  D. M., , Yang  L., , Young  A. N., , Wang  M. D., , and Nie  S., “ In vivo tumor targeting and spectroscopic detection with surface-enhanced Raman nanoparticle tags. ,” Nat. Biotechnol.. 26, , 83–90  ((2008)).
Tian  Y., , Cao  W.-H., , Luo  X.-X., , and Fu  Y., “ Preparation and luminescence property of Gd2O2S:Tb X-ray nano-phosphors using the complex precipitation method. ,” J. Alloys. Compd.. 433, , 313–317  ((2007)).
Xing  M., , Cao  W., , Pang  T., , Ling  X., , and Chen  N., “ Preparation and characterization of monodisperse spherical particles of X-ray nanophosphors based on Gd2O2S:Tb. ,” Chin. Sci. Bull.. 54, , 2982–2986  ((2009)).
Carpenter  C., , Pratx  G., , Sun  C., , and Xing  L., “ Hybrid x-ray/optical luminescence imaging: characterization of experimental conditions. ,” Med. Phys.. 37, , 4011–4018  ((2010)).
Pratx  G., , Carpenter  C., , Sun  C., , and Xing  L., “ X-ray luminescence computed tomography via selective excitation: a feasibility study. ,” IEEE Trans. Med. Imaging. 29, , 1992–1999  ((2010)).
Pratx  G., , Carpenter  C. M., , Sun  C., , Rao  R. P., , and Xing  L., “ Tomographic molecular imaging of x-ray-excitable nanoparticles. ,” Opt. Lett.. 35, , 3345–3347  ((2010)).
Ye  Y., , Yu  H., , Wei  Y., , and Wang  G., “ A general local reconstruction approach on a truncated Hilbert transform. ,” Int. J. Biomed. Imaging. 2007, , 63634  ((2007)).
Yu  H., and Wang  G., “ Compressive sensing based interior tomography. ,” Phys. Med. Biol.. 54, , 2791–2805  ((2009)).
Yang  J., , Yu  H., , Jiang  M., , and Wang  G., “ High-order total variation minimization for interior tomography. ,” Inverse Problems. 26, , 035013  ((2010)).
Boisseau  P., “ Determination of three dimensional trace element distributions by use of monochromatic x-ray microbeams. ,” PhD Thesis, Massachusetts Institute of Technology ((1986)).
Cheong  S.-K., , Jones  B. L., , Siddiqi  A. K., , Liu  F., , Manohar  N., , and Cho  S. H., “ X-ray fluorescence computed tomography (XFCT) imaging of gold nanoparticle-loaded objects using 110 kVp x-rays. ,” Phys. Med. Biol.. 55, , 647–662  ((2010)).
Als-Nielsen  J., and McMorrow  D.,  Elements of Modern X-ray Physics. ,  Wiley ,  New York (2001).
Hubbell  J. H., , Veigele  W. J., , Briggs  E. A., , Brown  R. T., , Cromer  D. T., , and Howerton  R. J., “ Atomic form factors, incoherent scattering functions, and photon scattering cross sections. ,”J. Phys. Chem. Ref. Data. 4, , 471–538  ((1975)).
Baro  J., , Sempau  J., , Fernandez-Varea  J. M., , and Salvat  F., “ 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)).
Wang  G., , Cong  W., , Shen  H., , and Zou  Y., “ Varying collimation for dark-field extraction. ,” Int. J. Biomed. Imaging. 2009, , 847537  ((2009)).
Welch  J., and van Gemert  M. J. C.,  Optical and Thermal Response of Laser-irradiated Tissue. ,  Plenum ,  New York  ((1995)).
Cong  W., , Wang  G., , Kumar  D., , Liu  Y., , Jiang  M., , Wang  L., , Hoffman  E., , McLennan  G., , McCray  P., , Zabner  J., , and Cong  A., “ Practical reconstruction method for bioluminescence tomography. ,” Opt. Express. 13, , 6756–6771  ((2005)).
Candes  E. J., , Romberg  J., , and Tao  T., “ Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. ,” IEEE Trans. Inf. Theory. 52, , 489–509  ((2006)).
Beck  A., and Teboulle  M., “ A fast iterative shrinkage-thresholding algorithm for linear inverse problems. ,” SIAM J. Imaging Sci.. 2, , 183–202  ((2009)).
van Eijk  C. W. E., “ Inorganic scintillators in medical imaging. ,” Phys. Med. Biol.. 47, , R85–R106  ((2002)).
Alexandrakis  G., , Rannou  F. R., , and Chatziioannou  A. F., “ Tomographic bioluminescence imaging by use of a combined optical-PET (OPET) system: a computer simulation feasibility study. ,” Phys. Med. Biol.. 50, , 4225–4241  ((2005)).
Prahl  S. A., “ Optical properties spectra. ,” Oregon Medical Laser Clinic, http://omlc.ogi.edu/spectra/index.html ((2001)).
Lansford  R., , Bearman  G., , and Fraser  S. E., “ Resolution of multiple green fluorescent protein color variants and dyes using two-photon microscopy and imaging spectroscopy. ,” J. Biomed. Opt.. 6, , 311–318  ((2001)).
Miqueles  E. X., and De Pierr  A. R., “ Exact analytic reconstruction in x-ray fluorescence CT and approximated versions. ,” Phys. Med. Biol.. 55, , 1007–1024  ((2010)).
Zhang  J., , Yang  G., , Cheng  Y., , Gao  B., , Qiu  Q., , Lee  Y. Z., , Lu  J. P., , and Zhou  O., “ Stationary scanning x-ray source based on carbon nanotube field emitters. ,” Appl. Phys. Lett.. 86, , 184104  ((2005)).

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

Related Content

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

Related Book Chapters

Topic Collections

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

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