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

## 1.

## Introduction

Novel 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 ﬂuorescence 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 Model

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.

## 2.1.

### 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).^{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

## 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} $$I\left(\mathbf{r}\right)={I}_{0}\mathrm{exp}\left[-{\int}_{0}^{1}{\mu}_{t}\left[\left(1-\tau \right)\phantom{\rule{1.0pt}{0ex}}{\mathbf{r}}_{0}+\tau \phantom{\rule{1.0pt}{0ex}}\mathbf{r}\right]d\tau \right],$$**r**

_{0}is a source position,

*I*

_{0}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:

## 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} $$\begin{array}{ccc}\hfill \mathrm{X}\left(\mathbf{r}\right)& =& {\int}_{L}I\left({\mathbf{r}}_{\mathbf{s}}\right)\eta \left({\mathbf{r}}_{s}\right)\frac{d{\sigma}_{scatter}}{d\Omega}\left(\beta \right)\hfill \\ & & \times \mathrm{exp}\left[-{\int}_{0}^{1}{\mu}_{t}\left[\left(1-\tau \right)\phantom{\rule{1.0pt}{0ex}}{\mathbf{r}}_{s}+\tau \phantom{\rule{1.0pt}{0ex}}\mathbf{r}\right]d\tau \right]\phantom{\rule{0.16em}{0ex}}d{\mathbf{r}}_{s},\hfill \end{array}$$*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 [TeX:] $\beta = \cos ^{ - 1} \left({\frac{{\left({{\bf r} - {\bf r}_s } \right) \cdot {\bm \nu }}}{{\left\| {{\bf r} - {\bf r}_s } \right\|}}} \right)$ $\beta ={\mathrm{cos}}^{-1}\left(\frac{\left(\mathbf{r}-{\mathbf{r}}_{s}\right)\xb7\mathit{\nu}}{\Vert \mathbf{r}-{\mathbf{r}}_{s}\Vert}\right)$ , [TeX:] ${\bm \nu }$ $\mathit{\nu}$ is the direction of the x-ray primary beam, η(

**r**

_{s}) is the electron density at the position

**r**

_{s}, and the x-ray differential scattering cross-section [TeX:] $\frac{{d\sigma _{scatter} }}{{d\Omega }}\left(\beta \right)$ $\frac{d{\sigma}_{scatter}}{d\Omega}\left(\beta \right)$ 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 × 10

^{23}/cm

^{3}, 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.

## 2.2.

### 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 transport^{17}

## 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} $$-\nabla \xb7\left[D\left(\mathbf{r}\right)\nabla \Phi \left(\mathbf{r}\right)\right]+{\mu}_{a}\left(\mathbf{r}\right)\Phi \left(\mathbf{r}\right)=S\left(\mathbf{r}\right),\phantom{\rule{1em}{0ex}}\mathbf{r}\in \Omega ,$$**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 [TeX:] $D = [ {3({\mu _a + \mu^{\prime}_s })} ]^{ - 1} $ $D={\left[3\left({\mu}_{a}+{\mu}_{s}^{\prime}\right)\right]}^{-1}$ , [TeX:] $\mu^{\prime}_s $ ${\mu}_{s}^{\prime}$ is the reduced scattering coefficient, and Ω⊂

*R*

^{3}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

## 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} $$\Phi \left(\mathbf{r}\right)+2\alpha D\left(\mathbf{r}\right)\left[\nu \xb7\nabla \Phi \left(\mathbf{r}\right)\right]=0,\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\phantom{\rule{0.16em}{0ex}}\mathbf{r}\in \partial \Omega ,$$*n*and air can be approximated by α = (1 + γ)/(1 − γ) with γ = −1.4399

*n*

^{−2}+ 0.7099

*n*

^{−1}+ 0.6681 + 0.0636

*n*.

^{17}The measurable exiting photon flux on the surface of the object can be expressed as

## 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} $$m\left(\mathbf{r}\right)=-D\left(\mathbf{r}\right)\left(\nu \xb7\nabla \Phi \left(\mathbf{r}\right)\right)\phantom{\rule{1em}{0ex}}\mathbf{r}\in \partial \Omega .$$**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,

^{5}

## 6

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} S\left({\bf r} \right) = \varepsilon \,{\rm X}\left({\bf r} \right)\,{\rm \rho }\left({\bf r} \right). \end{equation}\end{document} $$S\left(\mathbf{r}\right)=\varepsilon \phantom{\rule{0.16em}{0ex}}\mathrm{X}\left(\mathbf{r}\right)\phantom{\rule{0.16em}{0ex}}\rho \left(\mathbf{r}\right).$$## 2.3.

### Discretization

Equations 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}

## 7

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\bf A} \cdot \,\Phi = {\bf F} \cdot \,{\rm \rho,} \end{equation}\end{document} $$\mathbf{A}\xb7\phantom{\rule{0.16em}{0ex}}\Phi =\mathbf{F}\xb7\phantom{\rule{0.16em}{0ex}}\rho ,$$**A**is

## 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} $$\begin{array}{ccc}\hfill {a}_{ij}& =& {\int}_{\Omega}D\left(\mathbf{r}\right)\nabla {\varphi}_{i}\left(\mathbf{r}\right)\xb7\nabla {\varphi}_{j}\left(\mathbf{r}\right)\mathrm{d}\mathbf{r}+{\int}_{\Omega}{\mu}_{a}\left(\mathbf{r}\right){\varphi}_{i}\left(\mathbf{r}\right){\varphi}_{j}\left(\mathbf{r}\right)\mathrm{d}\mathbf{r}\hfill \\ & & +{\int}_{\partial \Omega}{\varphi}_{i}\left(\mathbf{r}\right){\varphi}_{j}\left(\mathbf{r}\right)/2\alpha \mathrm{d}\mathbf{r},\hfill \end{array}$$**F**is

## 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} $${f}_{ij}=\varepsilon {\int}_{\Omega}\mathrm{X}\left(\mathbf{r}\right){\varphi}_{i}\left(\mathbf{r}\right){\varphi}_{j}\left(\mathbf{r}\right)\mathrm{d}\mathbf{r},$$_{i}(

*i*= 1, 2, …) are the element shape functions. Since the matrix

**A**in Eq. 7 is positive definite, we have

## 10

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \Phi = ({{\bf A}\,^{ - 1} {\bf F}}) \cdot {\rm \rho }{\rm.} \end{equation}\end{document} $$\Phi =\left({\mathbf{A}\phantom{\rule{0.16em}{0ex}}}^{-1}\mathbf{F}\right)\xb7\rho .$$## 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} $$\sum _{i=1}^{n}\Phi \left({\mathbf{r}}_{{m}_{i}}\right)=\left(\sum _{{m}_{i}}{\mathbf{A}}^{-\mathbf{1}}\mathbf{F}\right)\xb7\rho ,$$*m*

_{i}(

*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:

## 12

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} Q_e^{ - 1} \;{\rm M} = {\bf G} \cdot {\rm \rho,} \end{equation}\end{document} $${Q}_{e}^{-1}\phantom{\rule{0.28em}{0ex}}\mathrm{M}=\mathbf{G}\xb7\rho ,$$*Q*

_{e}is the photon-detector quantum efficiency, and

**G**is a weighting matrix assembled with a row vector [TeX:] $\sum\nolimits_{m_i } {{\bf A}^{ - {\bf 1}} {\bf F}}$ ${\sum}_{{m}_{i}}{\mathbf{A}}^{-\mathbf{1}}\mathbf{F}$ for each x-ray excitation. Equation 12 describes the linear relationship between the nanophosphor distribution and the NIR measurement.

## 2.4.

### 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

## 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} $$\begin{array}{c}\hfill \begin{array}{c}\mathit{minimum}{\Vert \mathbf{T}\xb7\rho \Vert}_{1}\hfill \\ \mathit{subject}\phantom{\rule{0.16em}{0ex}}\mathit{to}{\Vert \mathbf{G}\xb7\rho -{Q}_{e}^{-1}\mathrm{M}\Vert}_{2}\le \tau ,\hfill \\ \phantom{\rule{2.em}{0ex}}\phantom{\rule{2.em}{0ex}}\phantom{\rule{0.28em}{0ex}}\rho \ge 0\hfill \end{array}\end{array}$$*τ*represents the data noise level. The presence of the

*l*

_{1}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.

## 3.

## Numerical Simulation

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 Gd_{2}O_{2}S doped with thulium.^{3, 4} Based on Ref. 21, the light yield ε was assumed to be about 0.15 cm^{3}/mg for 50 keV.

Appropriate optical parameters were accordingly assigned to the mouse model. The reduced scattering coefﬁcient
[TeX:]
$\mu^{\prime}_s \left(\lambda \right)$
${\mu}_{s}^{\prime}\left(\lambda \right)$
relies on the wavelength λ(*nm*) and is approximated by an empirical function

## 14

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \mu '_s \left(\lambda \right) = 10a \cdot \lambda ^{ - b}, \end{equation}\end{document} $${\mu}_{s}^{\prime}\left(\lambda \right)=10a\xb7{\lambda}^{-b},$$*a*and

*b*are the constants depending on the tissue type. The parameters

*a*and

*b*for organ-speciﬁc values can be found.

^{22}The tissue absorption depends on the local oxy-hemoglobin (HbO

_{2}), deoxy-hemoglobin (Hb), and water (W) concentrations in the respective organs. The spectral absorption coefﬁcient μ

_{a}(λ) can be approximated as the weighted sum of the three absorption coefﬁcients [TeX:] $\mu _{aHbO_2 } \left(\lambda \right)$ ${\mu}_{aHb{O}_{2}}\left(\lambda \right)$ , μ

_{aHb}(λ), and μ

_{aW}(λ), which were calculated from the corresponding absorbance spectra reported in Ref. 22

## 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} $${\mu}_{a}\left(\lambda \right)={S}_{B}\left[x{\mu}_{aHb}\left(\lambda \right)+\left(1-x\right){\mu}_{aHb{O}_{2}}\left(\lambda \right)\right]+{S}_{W}{\mu}_{aW}\left(\lambda \right),$$*x*=

*HbO*

_{2}/(

*HbO*

_{2}+

*Hb*) is the ratio between oxy-hemoglobin and the total hemoglobin concentration,

*S*

_{B}and

*S*

_{W}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.

## 3.1.

### 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 *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 }}}$ $\overline{e}=\frac{1}{N}{\sum}_{k=1}^{N}\frac{\left|{\Phi}_{k}^{S}-{\Phi}_{k}^{T}\right|}{{\Phi}_{k}^{T}}$ , where [TeX:] $\Phi _k^T$ ${\Phi}_{k}^{T}$ and [TeX:] $\Phi _k^S$ ${\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 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 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 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

## 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} $${Q}_{e,v}^{-1}{\mathrm{M}}_{\nu}={\mathbf{G}}_{\nu}\xb7{\rho}_{\nu},\phantom{\rule{1em}{0ex}}\nu =1,2,...,\kappa ,$$_{ν}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, [TeX:] ${\rm \rho }_\nu = \sum\nolimits_i {w_i p_{i,\nu } },\nu = 1,2, \ldots,\kappa$ ${\rho}_{\nu}={\sum}_{i}{w}_{i}{p}_{i,\nu},\nu =1,2,...,\kappa $ , where

*p*

_{i, ν}is the

*i*th probe concentration in a spectral band [ω

_{ν}, ω

_{ν + 1}] and

*w*

_{i}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.

## 4.

## Discussions and Conclusion

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.

## Acknowledgments

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

## References

_{2}O

_{2}S:Tb X-ray nano-phosphors using the complex precipitation method,” J. Alloys. Compd. 433, 313–317 (2007). 10.1016/j.jallcom.2006.06.075 Google Scholar

_{2}O

_{2}S:Tb,” Chin. Sci. Bull. 54, 2982–2986 (2009). 10.1007/s11434-009-0492-9 Google Scholar