Next Article in Journal
Thermodynamic Work of High-Grade Uterine Prolapse Patients Undergoing Transvaginal Mesh Repair with Total Hysterectomy
Next Article in Special Issue
FTSNet: Fundus Tumor Segmentation Network on Multiple Scales Guided by Classification Results and Prompts
Previous Article in Journal
Digital Dental Biometrics for Human Identification Based on Automated 3D Point Cloud Feature Extraction and Registration
Previous Article in Special Issue
Integrating Multi-Organ Imaging-Derived Phenotypes and Genomic Information for Predicting the Occurrence of Common Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dual and Multi-Target Cone-Beam X-ray Luminescence Computed Tomography Based on the DeepCB-XLCT Network

1
Biomedical Engineering Department, Fourth Military Medical University, Xi’an 710032, China
2
Shaanxi Provincial Key Laboratory of Bioelectromagnetic Detection and Intelligent Perception, Xi’an 710032, China
3
School of Software Engineering, Xi’an Jiaotong University, Xi’an 710049, China
*
Authors to whom correspondence should be addressed.
These authors contribute equally to this work.
Bioengineering 2024, 11(9), 874; https://doi.org/10.3390/bioengineering11090874
Submission received: 31 July 2024 / Revised: 24 August 2024 / Accepted: 27 August 2024 / Published: 28 August 2024

Abstract

:
Background and Objective: Emerging as a hybrid imaging modality, cone-beam X-ray luminescence computed tomography (CB-XLCT) has been developed using X-ray-excitable nanoparticles. In contrast to conventional bio-optical imaging techniques like bioluminescence tomography (BLT) and fluorescence molecular tomography (FMT), CB-XLCT offers the advantage of greater imaging depth while significantly reducing interference from autofluorescence and background fluorescence, owing to its utilization of X-ray-excited nanoparticles. However, due to the intricate excitation process and extensive light scattering within biological tissues, the inverse problem of CB-XLCT is fundamentally ill-conditioned. Methods: An end-to-end three-dimensional deep encoder-decoder network, termed DeepCB-XLCT, is introduced to improve the quality of CB-XLCT reconstructions. This network directly establishes a nonlinear mapping between the distribution of internal X-ray-excitable nanoparticles and the corresponding boundary fluorescent signals. To improve the fidelity of target shape restoration, the structural similarity loss (SSIM) was incorporated into the objective function of the DeepCB-XLCT network. Additionally, a loss term specifically for target regions was introduced to improve the network’s emphasis on the areas of interest. As a result, the inaccuracies in reconstruction caused by the simplified linear model used in conventional methods can be effectively minimized by the proposed DeepCB-XLCT method. Results and Conclusions: Numerical simulations, phantom experiments, and in vivo experiments with two targets were performed, revealing that the DeepCB-XLCT network enhances reconstruction accuracy regarding contrast-to-noise ratio and shape similarity when compared to traditional methods. In addition, the findings from the XLCT tomographic images involving three targets demonstrate its potential for multi-target CB-XLCT imaging.

1. Introduction

With the emergence of X-ray-excitable nanoparticles, X-ray luminescence computed tomography (XLCT) has garnered increased attention due to its impressive performance [1,2]. In XLCT, X-ray-excitable nanoparticles serve as imaging probes that emit visible or near-infrared (NIR) light upon X-ray irradiation, detectable by an electron-multiplying charge-coupled device (EMCCD) camera. By addressing an inverse problem with a suitable imaging model for X-ray and light photon transport, it is possible to reconstruct the three-dimensional (3D) distribution of nanoparticles within the scanned object. In contrast to bioluminescence tomography (BLT) and fluorescence molecular tomography (FMT), XLCT achieves greater imaging depth due to the superior penetration and collimation of X-rays. Furthermore, employing X-ray-excited nanoprobes effectively mitigates interference from autofluorescence and background fluorescence, thereby enhancing imaging contrast and resolution [3,4,5,6]. The above advantages have led to the rapid development of XLCT, which has broad application prospects in disease diagnosis, as well as in new drug research and development [7,8].
Since the initial demonstration of XLCT, various imaging geometries have been suggested to enhance spatial and temporal resolution. Cone-beam XLCT (CB-XLCT) accelerates the imaging process considerably compared to narrow-beam and fan-beam systems at the expense of spatial resolution, in which the optical imaging acquisition process can complete within two minutes [9,10]. Due to high scattering of light in biological tissues, the first-order approximation model of the radiative transfer equation (RTE) is currently applied to describe the photon propagation in the forward model of CB-XLCT [11,12,13]. Considering the complicated excitation process of XLCT, the deviation between the approximate linear model of RTE and the true nonlinear photon excitation and propagation could not be avoided. To improve its reconstruction quality, the structured prior knowledge or sparse regularization have been employed to penalize the reconstruction objective function [14,15,16,17], which is mostly solved by statistically iterative algorithms. However, model-based reconstruction may struggle to accurately represent the complex and ill-posed nature of CB-XLCT imaging, constraining its application in in vivo imaging.
In recent years, deep learning techniques have significantly advanced the capabilities of structural and molecular imaging modalities, including low-dose CT [18,19,20,21], bioluminescence tomography (BLT) [22,23], and fluorescence molecular tomography (FMT) [24,25]. By providing a large amount of training data to the neural network, the reconstruction architecture learns how to best utilize the data. Compared to traditional reconstruction methods, deep-learning-based reconstruction can effectively avoid errors caused by inaccurate modeling and ill-posed nature in solving inverse problems, achieving better performance than traditional image reconstruction techniques. However, deep learning technology has yet to be applied to CB-XLCT reconstruction.
In this paper, an end-to-end three-dimensional deep encoder–decoder network (DeepCB-XLCT) is proposed for CB-XLCT reconstruction, in which large datasets are utilized to learn the unknown solution to the inverse problem. The DeepCB-XLCT network is designed to establish a nonlinear mapping from input to output, with the parameters of this mapping continuously studied and adjusted during network training. Based on this method, the inaccuracy arising from constructing the photon propagation model or addressing the ill-posed inverse problem can be effectively mitigated.
The rest of this paper is structured as follows. In Section 2, the conventional forward model and inverse problem of CB-XLCT, the proposed deep neural network for the CB-XLCT imaging model, and the training data and optimization training procedure are described in detail. In Section 3, numerical simulations and phantom experiments are performed for evaluating the performance of the proposed reconstruction approach. The results are presented in Section 4, followed by discussions and conclusions in Section 5.

2. Methods

2.1. Conventional Forward Model and Inverse Problem of CB-XLCT

For CB-XLCT imaging, nanoparticles within an object emit visible or NIR light upon X-ray irradiation. According to previous studies [26], the quantity of optical photons emitted is directly proportional to both the intensity distribution of the X-rays and the concentration of nanoparticles within the object, expressed as follows:
S ( r ) = Γ X ( r ) n ( r )
where S(r) represents the emitted light from point r, Γ denotes the light yield of the nanoparticles, n(r) is the concentration of nanoparticles, and X(r) signifies the intensity of X-rays at position r, which can be described using the Lambert–Beer law [14]
X ( r ) = X ( r 0 ) exp { r 0 r μ t ( τ ) }
where X ( r 0 ) represent the initial X-ray intensity at positions r 0 , and μt(τ) shows the X-ray attenuation coefficient, which can be derived from the CT data.
Due to the characteristics of high scattering and low absorbing in the visible and NIR spectral window for biological tissues, the diffusion equation (DE) is used to establish the propagation model of the emitted light [27], which can be expressed as
[ D ( r ) Φ ( r ) ] + μ a ( r ) Φ ( r ) = S ( r ) ( r Ω )
where D(r) is the diffusion coefficient, which can be calculated by D ( r ) = 1 / [ 3 ( μ s ( r ) + μ a ( r ) ) ] , in which μ s ( r ) is the reduced scattering coefficient, μa(r) is the absorption coefficient, Φ(r) represents the photon fluence, and Ω is the image domain.
To solve the diffusion Equation (3), Robin boundary conditions are typically applied [28,29], and can be expressed as
Φ ( r ) + 2 κ D ( r ) [ ν Φ ( r ) ] = 0 ( r Ω )
where κ is the boundary mismatch parameter, ν represents the outward unit normal vector on the boundary, and Ω is the boundary of Ω .
Based on the finite element method (FEM), the forward model based on Equations (1)–(4) for the practical application of XLCT can be established as
y = W x + ς
where y denotes the actual fluorescence signals measured on the surface of the imaging object, W represents a weight matrix, formulated based on previous studies [2], ς signifies the noise inherent in the system, and x denotes the unknown distribution of nanoparticles within the imaging object.
To solve the ill-posed inverse problem of Equation (5), the reconstructing method based on the regularization terms are used in the CB-XLCT reconstruction, which can be expressed as
x = arg min x 0 W x y 2 2 + λ x β
where λ represents the regularization parameter, and β (0 < β < 2) indicates the norm term. When β is equal to 1 and 2, it refers to L1 regularization and L2 regularization, respectively.

2.2. Deep Neural Network for CB-XLCT

In contrast to traditional methods, CB-XLCT reconstruction utilizing deep neural networks does not focus on explicitly solving the forward and inverse problems. Rather, it creates an end-to-end deep neural network mapping model that directly reconstructs the distribution of fluorescence sources. The end-to-end three-dimensional deep encoder–decoder network (DeepCB-XLCT) is composed of a 3D-Encoder network and a 3D-Decoder network, featuring multiple layers of spatial convolution and deconvolution operations. The input to the network comprises 2D surface measurement images from 24 angles of CB-XLCT, while the output is the reconstructed 3D distribution of CB-XLCT nanoparticles. The details of the network structure proposed in this work are shown in Figure 1. The fundamental architecture of the 3D Encoder–Decoder network is divided into two components: the 3D-Encoder network, which captures the distribution of optical features and other codes, and the 3D-Decoder network, which learns the conditional distribution of y in order to generate a specific x. The 3D-Encoder network is composed of five convolutional blocks. Each block includes a 3D convolutional layer, followed by batch normalization and LeakyReLU activation function, as well as a max pooling layer to progressively reduce spatial dimensions. The 3D-Decoder network uses transposed convolutional layers to progressively upscale the feature maps while reducing the number of channels, consisting of five layers. Each transposed convolutional layer is typically followed by batch normalization and a ReLU activation function, enhancing the network’s non-linear expressive capacity and stability. There are two fully connected layers between the Encoder network and Decoder network, which transform the encoded features to another feature space and then back again, preparing for feature fusion in the decoder. Moreover, the dual-sampling module [24] is integrated at various stages of the 3D-Decoder network, which makes the DeepCB-XLCT more focused on effective features when processing volumetric data, enhancing feature representation and overall network performance.

2.3. Generation and Preparation of Training Data

In order to obtain a large amount of training data, the forward model of CB-XLCT is solved using computer simulation programs, in which the intensity distribution of the X-rays (Equation (2)) and the number of optical photons emitted (Equation (1)) were solved through analytical calculation, and the propagation model of the emitted light (Equation (3)) was solved based on the finite element method (FEM) [5,14]. To ensure consistency between the simulated data and the actual experiments, the finite element meshes were created based on the phantom or mouse model utilized in the real experiment.
In the phantom experiments, the training data were generated by randomly introducing two or three cylindrical fluorescent targets of varying sizes and locations. These targets, with diameters of either 3 or 4 mm and a fixed height of 4 mm, were distributed randomly throughout the phantoms. The horizontal and vertical coordinates of these fluorescent target centers were randomly generated within the range of −10 to 10 mm, while the edge-to-edge distances (EEDs) ranged from 0.3 to 2.5 mm. For the mouse experiments, training data were created by adding two cylindrical fluorescent targets of different sizes and positions at random. Similarly, cylindrical fluorescent targets with diameters of 3 or 4 mm and a fixed height of 4 mm were randomly positioned inside the mouse model, with their center coordinates generated between −8 and 8 mm to conform to the specifications of the mouse model. The X-ray source and EMCCD configuration of the data were adjusted to replicate the conditions of real experiments. Finally, in the phantom experiments, a total of 8000 simulation samples were generated, with 6000 samples allocated for training and 2000 samples reserved for validation to optimize the model. Similarly, the mouse experiments involved generating 4000 simulation samples, with 3000 samples used for training and 1000 samples set aside for validation.
Furthermore, all projection images were resized to 128 × 128 pixels prior to being input into the network. The network’s output is a 3D image, with each x-y slice measuring 128 × 128 pixels and a reconstruction resolution of 1 mm along the z-axis. To summarize, the input shape is 128 × 128 × 24 (as illustrated in Figure 1) for the projection images, while the output shape is 128 × 128 × 7 for the reconstructed tomographic images.

2.4. Optimization Training Procedure of DeepCB-XLCT

The DeepCB-XLCT networks were implemented using PyTorch 1.12 and Python 3.9 on an Ubuntu 18.04.3 system. All training and testing procedures were conducted using two NVIDIA-A100 with 80 GB memory for each. The objective function of the DeepCB-XLCT network consists of two parts: the mean square errors (MSE) and structural similarity loss (SSIM) between the output of DeepCB-XLCT and true results. In addition, due to the smaller proportion of the target region compared to the entire reconstruction area as shown in Figure 2, the loss of the target regions (ROI) is also considered.
The Structural Similarity Index (SSIM) [30] is a perceptually motivated metric that evaluates the similarity between two images by comparing corresponding pixels and their surrounding neighborhoods. Specifically, the SSIM index employs a sliding Gaussian window to assess the similarity of each local patch within the images. The SSIM of the same pixels in the output of DeepCB-XLCT and true results is defined as
S S I M ( x r i , x t i ) = ( 2 μ x r i μ x t i + C 1 ) ( 2 σ x r i x t i + C 2 ) ( μ x r i 2 + μ x t i 2 + C 1 ) ( σ x r i 2 + σ x t i 2 + C 2 )
where μ x r i and μ x t i are the means of pixel i in the output of DeepCB-XLCT and true results, respectively, σ x r i and σ x t i are the variances of of pixel i. C1 and C2 are constants used to stabilize the divisions. xr is the output of DeepCB-XLCT, and xt is the true results.
To compute the SSIM index between the output of DeepCB-XLCT and the true results, the mean is taken over the SSIM indices of all local patches.
M S S I M ( x r , x t ) = 1 M i = 1 M S S I M ( x r i , x t i )
where M is the number of local patches in the output of DeepCB-XLCT and true results.
Then, the SSIM loss between the output of DeepCB-XLCT and true results can be defined as
S S I M ( x r , x t ) = 1 M S S I M ( x r , x t )
In the end, the objective function of the DeepCB-XLCT network can be defined as
D e e p C B X L C T = M S E ( x r , x t ) + M S E ( x r , x t ) + 2 · [ S S I M ( x r , x t ) + S S I M ( x r , x t ) ]
where x r and x t are the target regions in the output of DeepCB-XLCT and true results, respectively. The MSE loss between the output of DeepCB-XLCT and the true results can be calculated as
M S E ( x r , x t ) = 1 N i = 1 N ( x r i x t i )
where N is the number of training samples.
In this study, the network utilized the Adam algorithm as its optimizer. The batch size was 64, which denoted the number of samples input to the model at once during the training process. The epochs ware set to 200 based on the tradeoff between enough training and avoiding overfitting. The initial learning rate is 2 × 10−5, and the learning rate decays based on the validation loss value, which the learning rate decreases by 2 times when the validation loss remains unchanged for 5 cycles. Under these parameter conditions, it takes 3 h to train the DeepCB-XLCT model.

3. Experimental Design

Numerical simulations and phantom experiments were conducted to assess the performance of the proposed DeepCB-XLCT network in CB-XLCT reconstruction using the custom-developed CB-XLCT system. Additionally, for comparison, four traditional methods, namely, adaptive FISTA (ADFISTA, L1 norm), MAP (Gaussian model) [14,31,32], T-FISTA (L1 norm) [33], and ADMLEM (Poisson model) [2], were implemented to solve Equation (6). In this study, the gradient variation step of the ADFISTA method was set as 0.01, and the iteration numbers were determined by an adaptive iteration termination condition. The MAP method was established based on Gaussian Markov random field, in which the hyperparameters were alternately estimated in each iteration and the objective function was minimized based on a voxel-wise iterative coordinate descent (ICD) algorithm. The regularization parameters and iterative numbers of the T-FISTA method were empirically set as 0.1 and 300, respectively. The ADMLEM method was established based on the Poisson distributed projections, in which the iterative number was set as 800.

3.1. Numerical Simulations Setup

Numerical simulations were first conducted using a cylinder phantom to evaluate the performance of the proposed DeepCB-XLCT network. To simulate the environment of biological tissues, the cylindrical phantom, measuring 3.0 cm in diameter and 2.3 cm in height, was filled with a mixture of water and intralipid, as shown in Figure 3. The absorption coefficient and reduced scattering coefficient were specified as 0.02 and 10 cm−1, respectively. Additionally, two small tubes (4 mm in diameter and 4 mm in height) containing Y2O3:Eu3+ (50 mg/mL) were positioned within the cylinder phantom as the luminescent targets. To assess the resolution capabilities of the proposed DeepCB-XLCT network, the edge-to-edge distances (EEDs) between the two targets were set as 2.0, 1.5, and 1.0 mm.
In the numerical simulations, 24 projections were obtained every 15° during a 360° scan, with the voltage and current of the cone beam X-ray source set at 40 kV and 1 mA, respectively. The exposure time of the EMCCD camera was fixed at 0.5 s. Following acquisition of optical luminescence data from various angles, white Gaussian noise with a zero-mean and a signal-to-noise ratio (SNR) of 30 dB was added to all projections to simulate noisy measurements. Moreover, in order to further validate the performance of the proposed DeepCB-XLCT method under varying noise levels, the projections were also added with 25 dB, 20 dB, and 15 dB white Gaussian noise.

3.2. Phantom Experiments Setup

To further assess the performance of the proposed DeepCB-XLCT network using actual luminescence measurements, phantom experiments were conducted with a custom-developed CB-XLCT system, which mainly included a rotation stage, a micro-focus X-ray source (Oxford Instrument, Oxford, UK), an electron-multiplying charge-coupled device (EMCCD) camera (iXon DU-897, Andor, Oxford, UK) for capturing luminescence signals, and a flat-panel X-ray detector (2923, Dexela, London, UK) for capturing X-ray signals, as depicted in Figure 4. The X-ray source operated at a maximum voltage of 80 kV with a maximum power of 80 W.
Figure 5 illustrates the setup of the physical phantom utilized in the imaging experiments. A glass cylinder measuring 4.0 cm in diameter and 4.0 cm in height, filled with a mixture of 1% intralipid and water, was secured on the rotation stage. The absorption coefficient and reduced scattering coefficient of the solution were 0.02 cm−1 and 10 cm−1, respectively. Additionally, two small glass tubes (3 mm in diameter) containing Y2O3:Eu3+ (50 mg/mL) were symmetrically positioned within the cylinder as the luminescent targets, with edge-to-edge distances between the two tubes of 2.3 mm, 1.7 mm, and 1.0 mm.
In the phantom experiments, the X-ray source was configured with a tube voltage of 40 kV and a tube current of 1 mA. Optical images were captured at every 15° as the phantom rotated from 0° to 360° using the EMCCD camera. The exposure time, EM gain, and binning of the EMCCD were configured to 0.5 s, 260, and 1 × 1, respectively. For CT imaging, 360 projections were acquired with an angular increment of 1° spanning from 0° to 360° by the X-ray detector, with each projection lasting 150 ms. Subsequently, the Feldkamp–Davis–Kress (FDK) algorithm [34,35] was employed to reconstruct the conventional CT image of the tubes for validation.

3.3. Quantitative Evaluation

The CB-XLCT images reconstructed using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network were evaluated and compared based on quantitative metrics, namely, the Dice Similarity Coefficient (DICE) and Contrast-to-Noise Ratio (CNR) [12].
The DICE coefficient is employed to assess the similarity between the true and reconstructed targets, and it can be calculated using the following formula:
D I C E = 2 | R O I r R O I t | | R O I r | + | R O I t |
where ROIt and ROIr represent the regions corresponding to the true and reconstructed targets, respectively, and |·| indicates the number of voxels within a given region.
The CNR is utilized for the quantitative assessment of contrast and noise in the reconstructed images, as outlined in Equation (13).
C N R = | μ R O I μ B C K | ( w R O I σ R O I 2 + w B C K σ B C K 2 ) 1 / 2
where ROI and BCK represent the target and background regions, respectively, and μ R O I and σ R O I 2 are the mean intensity values and the variances within the ROI μ B C K and σ B C K 2 represent the mean intensity values and variances of the BCK, and w R O I and w B C K serve as weighting factors based on the relative volumes of the ROI and BCK, respectively, which satisfy w R O I + w B C K = 1.

4. Results

4.1. Numerical Simulations

4.1.1. Resolution Experiment

Firstly, CB-XLCT reconstruction results were obtained for targets placed at varying edge-to-edge distances (EEDs) using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, as depicted in Figure 6. Normalization of all reconstruction results was performed based on their respective maximum intensity values. It can be observed that in the cases of the ADFISTA, MAP, and T-FISTA algorithms, distinguishing between the two targets proved challenging. With the ADMLEM method, while the distribution of the two targets could be differentiated, it exhibited relatively large localization errors. In contrast, the reconstruction using the proposed DeepCB-XLCT network demonstrated improved shape similarity to the ground truth and resulted in less noise compared to the other four algorithms.
Table 1 provides a comprehensive overview of the quantitative assessment of reconstructions employing various methods. Notably, the results obtained using the proposed DeepCB-XLCT network demonstrated the highest reconstruction quality in terms of CNR and shape similarity among all five methods, thereby reinforcing the observations depicted in Figure 6. Figure 7 plots the profiles along the green dotted line shown in Figure 6. The results again reveal that the values of the reconstructed images by the proposed DeepCB-XLCT network are closer to the ground truth images.

4.1.2. Robustness Experiment in Different Noise Levels and the Ablation Experiment

To further assess the performance of the proposed method across varying noise levels, XLCT tomographic images from two-target simulations with SNR values of 25 dB, 20 dB, and 15 dB were reconstructed, as illustrated in Figure 8. The quantitative evaluation of these reconstructions is summarized in Table 2. It is evident that the reconstruction results obtained using the proposed DeepCB-XLCT network maintained an acceptable quality concerning CNR and shape similarity, even when the SNR decreased to 15 dB. This finding further substantiates the enhancements provided by the proposed algorithm.
To investigate the impact of skip connection module and ROILoss on the reconstruction results, the ablation experiments were conducted based on the targets positioned at varying distances, as shown in Figure 9. Table 3 provides a summary of the quantitative evaluation of the reconstructions in ablation experiment. It can be seen without skip connection module and ROILoss, the quality of reconstruction results significantly decreased. Compared with skip connection module, ROILoss had a greater impact on the reconstruction results.

4.1.3. Multi-Target Experiment

Subsequently, XLCT tomographic images of three targets were reconstructed using various algorithms to assess the effectiveness of the proposed method in the context of multi-target scenarios, as illustrated in Figure 10. Notably, none of the algorithms, except for the proposed DeepCB-XLCT network, could effectively distinguish between the targets, which further confirms the enhancements brought by the proposed algorithm for XLCT tomographic images involving multiple targets.

4.2. Phantom Experiments

To further assess the efficacy of the proposed DeepCB-XLCT network, XLCT tomographic images of targets with varying EEDs were reconstructed from phantom experiments using different algorithms, as depicted in Figure 11. It is evident that, for ADFISTA and MAP algorithms shown in the first and second columns of Figure 11, distinguishing the distribution of the two targets proves to be challenging. Comparatively, while the T-FISTA and ADMLEM algorithms were able to effectively distinguish between the distributions of two targets at EEDs of 2.3 mm and 1.7 mm, they struggled to do so when the EED was 1.0 mm. Notably, the proposed DeepCB-XLCT network successfully resolved both targets as expected, as illustrated in the fifth column of Figure 11.
To further assess the proposed algorithm, Table 4 provides a quantitative evaluation of the reconstructions obtained through various methods in the phantom experiments. The results demonstrate that, in comparison to conventional reconstruction methods, the proposed DeepCB-XLCT network consistently outperforms in terms of shape recovery and image contrast in phantom experiments, which further corroborates the advancements made by the proposed DeepCB-XLCT network.

4.3. In Vivo Experiments and Results

To substantiate the efficacy of the proposed method in vivo, CB-XLCT experiments were performed on a female BALB/c nude mouse using our custom-developed CB-XLCT system. All procedures adhered strictly to the guidelines established by the Animal Ethics Review Committee of the Fourth Military Medical University. Two small glass tubes, measuring 3 mm in diameter and 4 mm in height, filled with Y2O3:Eu3+ at a concentration of 50 mg/mL, were implanted into the mouse’s abdominal cavity to function as fluorescent targets, as illustrated in Figure 12a,b. During the imaging experiments, the mouse was secured on the rotation stage, and 24 projections were acquired at 15° intervals from 0° to 360° by the EMCCD camera, while the X-ray source operated at 80 kV and 1 mA, respectively. The exposure time and EM gain of the EMCCD camera were set as 1 s and 260, respectively. Then, the CT projections were also collected to reconstruct the localization and structure of the tubes for verification, as shown in Figure 12b.
The XLCT tomographic images were reconstructed using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT method, as illustrated in Figure 12c–l. All fluorescence images were normalized according to the maximum reconstruction results and presented using a uniform color scale, as depicted in Figure 12c–g. Figure 12h–l presents the fusion outcomes of the XLCT and XCT reconstructed images. In comparison to the ADFISTA, MAP, T-FISTA, and ADMLEM method, the proposed DeepCB-XLCT method exhibits enhanced reconstruction performance in terms of image contrast and shape recovery.
Table 5 provides a quantitative evaluation summary of the reconstructions obtained using various methods. The proposed DeepCB-XLCT achieved the highest values for CNR and DICE, indicating its capability to produce high-quality reconstruction results in mouse experiments. Table 6 summarizes the reconstruction time of ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT method in the mice experiments. Compared with the conventional iterative methods, the reconstruction time of the proposed DeepCB-XLCT method is significantly reduced due to the absence of complex iterative calculations.

5. Discussion and Conclusions

CB-XLCT, as a novel hybrid imaging modality, offers significant advantages primarily due to its use of X-rays, which can enhance excitation depth, reduce tissue autofluorescence, and achieve a dual-mode imaging capability that combines X-ray computed tomography (CT) with optical molecular tomographic imaging. The performance of CB-XLCT reconstruction has a substantial effect on imaging outcomes. Nevertheless, the deviation between the complex imaging process and the approximate photon transmission model leads to the high ill-posed nature of the inverse problem, which greatly limits the improvement of imaging quality.
Conventionally, researchers have developed various reconstruction methods to constrain the image and improve the reconstruction quality, including ADFISTA, Bayesian theory, T-FISTA, and ADMLEM algorithms to mitigate the ill-posed nature of the inverse problem. In contrast to these conventional methods, this study introduces a novel DeepCB-XLCT approach, in which the CB-XLCT reconstruction process is completed by establishing an end-to-end nonlinear mapping model that relates the internal CB-XLCT nanoparticles to the surface measurement signals. Based on the DeepCB-XLCT method, it can greatly eliminate the modeling error associated with the forward problem and effectively minimize artifacts generated by iterative calculations, thereby greatly addressing the challenges posed by the ill-posed inverse problem. Moreover, the structural similarity loss (SSIM) was added to the objective function of the DeepCB-XLCT network to better restore the target shape, and the loss of the target regions was also added to make the DeepCB-XLCT network more focused on the target, which were validated based on the ablation experiment. Finally, it is worth mentioning that a certain level of noise in the data will not compromise the robustness of the DeepCB-XLCT method. In summary, the proposed DeepCB-XLCT method can achieve improved reconstruction results and robust reconstruction performance.
Numerical simulations, phantom experiments, and in vivo experiments results validate the superiority of the proposed method compared to traditional approaches such as the ADFISTA, MAP, T-FISTA, and ADMLEM methods. When the EED is 1.0 mm, two targets can be distinctly resolved (Figure 6 and Figure 11, Table 1 and Table 4), highlighting its effectiveness in enhancing spatial resolution. The reconstruction results in different noise levels (Figure 8, Table 2) further confirmed the robustness of the proposed DeepCB-XLCT method. The results of the ablation experiment (Figure 9, Table 3) indicate that the ROILoss has a greater influence on the reconstruction outcomes. Moreover, the results from the XLCT tomographic images with three targets (Figure 10) reveal the method’s capability for multi-target CB-XLCT imaging. To further assess the effectiveness of the proposed method, in vivo experiments were performed, and the results (Figure 12, Table 5) provide additional confirmation of its superiority over the traditional ADFISTA, MAP, T-FISTA, and ADMLEM methods, which would significantly enhance the preclinical application of CB-XLCT in small animals.
However, there are still some limitations. Firstly, the performance of the proposed DeepCB-XLCT method has been validated through phantom experiments with varying EEDs. To fully demonstrate the method’s superior capabilities, it is better to verify its performance through in vivo experiments involving different targets in the future. Secondly, as the number of reconstruction targets increase, the complexity of network design and training may also increase. Moreover, while a well-trained model can reconstruct the distribution of XLCT nanoparticles in a very short time, the training process for the network can be time-consuming, taking several hours. In the future, we will concentrate on addressing these challenges and advancing our research in several areas, including the development of improved network structures, conducting evaluations through in vivo experiments, and testing the robustness of the proposed method using various phantom and in vivo studies.
In summary, we have proposed a DeepCB-XLCT method to enhance the quality of CB-XLCT reconstruction. It directly constructs a nonlinear mapping relationship between the distribution of X-ray excitable nanoparticles and the boundary fluorescent signal distribution, which could effectively reduce the reconstruction inaccuracies associated with simplified linear models and iterative calculations. The results from numerical simulations, phantom experiments, and in vivo experiments have demonstrated that the proposed DeepCB-XLCT method can improve spatial resolution and reconstruction accuracy compared to conventional iterative methods, which can advance the preclinical application of CB-XLCT in small animals.

Author Contributions

Methodology, T.L. and S.H.; software, P.G.; validation, T.L. and S.H. resources, R.L.; data curation, S.H.; writing—original draft preparation, T.L. and S.H.; writing—review and editing, H.L., Y.S. and J.R.; visualization, W.L.; funding acquisition, T.L., J.R., P.G., H.L. and W.L. All authors have read and agreed to the published version of the manuscript.

Funding

Natural National Science Foundation of China (NSFC) (12305350, 12275355, 12405357); National Key Research and Development Program (2023YFF0715400, 2023YFF0715401); Key Research and Development Program of Shaanxi Province (S2024-YF-YBSF-0283).

Institutional Review Board Statement

The animal study protocol was approved by the Animal Ethics Review Committee of the Fourth Military Medical University (20230353, 2023.02).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest related to this article.

References

  1. Carpenter, C.; Sun, C.; Pratx, G.; Rao, R.; Xing, L. Hybrid x-ray/optical luminescence imaging: Characterization of experimental conditions. Med. Phys. 2010, 37, 4011–4018. [Google Scholar] [CrossRef] [PubMed]
  2. Liu, T.; Ruan, J.; Rong, J.; Hao, W.; Li, W.; Li, R.; Zhan, Y.; Lu, H. Cone-beam X-ray luminescence computed tomography based on MLEM with adaptive FISTA initial image. Comput. Methods Programs Biomed. 2022, 229, 107265. [Google Scholar] [CrossRef]
  3. Zhang, H.; Hai, L.; Kou, J.; Hou, Y.; He, X.; Zhou, M.; Geng, G. OPK_SNCA: Optimized prior knowledge via sparse non-convex approach for cone-beam X-ray luminescence computed tomography imaging. Comput. Methods Programs Biomed. 2022, 215, 106645. [Google Scholar] [CrossRef]
  4. Lun, M.C.; Arifuzzaman, M.; Ranasinghe, M.; Bhattacharya, S.; Anker, J.N.; Wang, G.; Li, C. Focused x-ray luminescence imaging system for small animals based on a rotary gantry. J. Biomed. Opt. 2021, 26, 036004. [Google Scholar] [CrossRef]
  5. Liu, T.; Rong, J.; Gao, P.; Pu, H.; Zhang, W.; Zhang, X.; Liang, Z.; Lu, H. Regularized reconstruction based on joint L 1 and total variation for sparse-view cone-beam X-ray luminescence computed tomography. Biomed. Opt. Express 2019, 10, 1–17. [Google Scholar] [CrossRef]
  6. Zhang, H.; Huang, X.; Zhou, M.; Geng, G.; He, X. Adaptive shrinking reconstruction framework for cone-beam X-ray luminescence computed tomography. Biomed. Opt. Express 2020, 11, 3717–3732. [Google Scholar] [CrossRef]
  7. Pei, P.; Chen, Y.; Sun, C.; Fan, Y.; Yang, Y.; Liu, X.; Lu, L.; Zhao, M.; Zhang, H.; Zhao, D.; et al. X-ray-activated persistent luminescence nanomaterials for NIR-II imaging. Nat. Nanotechnol. 2021, 16, 1011–1018. [Google Scholar] [CrossRef] [PubMed]
  8. Hong, Z.; Chen, Z.; Chen, Q.; Yang, H. Advancing X-ray Luminescence for Imaging, Biosensing, and Theragnostics. Acc. Chem. Res. 2022, 56, 37–51. [Google Scholar] [CrossRef] [PubMed]
  9. Chen, D.; Zhu, S.; Yi, H.; Zhang, X.; Chen, D.; Liang, J.; Tian, J. Cone beam x-ray luminescence computed tomography: A feasibility study. Med. Phys. 2013, 40, 031111. [Google Scholar] [CrossRef]
  10. Liu, X.; Liao, Q.; Wang, H. In vivo x-ray luminescence tomographic imaging with single-view data. Opt. Lett. 2013, 38, 4530–4533. [Google Scholar] [CrossRef]
  11. Chen, D.; Meng, F.; Zhao, F.; Xu, C. Cone Beam X-Ray Luminescence Tomography Imaging Based on KA-FEM Method for Small Animals. Bio. Med. Res. Int. 2016, 2016, 6450124. [Google Scholar] [CrossRef] [PubMed]
  12. Liu, T.; Rong, J.; Gao, P.; Zhang, W.; Liu, W.; Zhang, Y.; Lu, H. Cone-beam x-ray luminescence computed tomography based on x-ray absorption dosage. J. Biomed. Opt. 2018, 23, 026006. [Google Scholar] [CrossRef]
  13. Guo, H.; Zhao, H.; Yu, J.; He, X.; He, X.; Song, X. X-ray Luminescence Computed Tomography Using a Hybrid Proton Propagation Model and Lasso-LSQR Algorithm. J. Biophotonics 2021, 14, e202100089. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, G.; Liu, F.; Liu, J.; Luo, J.; Xie, Y.; Bai, J.; Xing, L. Cone beam x-ray luminescence computed tomography based on Bayesian method. IEEE Trans. Med. Imaging 2017, 36, 225–235. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, X.; Wang, H.; Xu, M.; Nie, S.; Lu, H. A wavelet-based single-view reconstruction approach for cone beam x-ray luminescence tomography imaging. Biomed. Opt. Express 2014, 5, 3848–3858. [Google Scholar] [CrossRef]
  16. Gao, P.; Rong, J.; Liu, T.; Zhang, W.; Lu, H. Limited view cone-beam x-ray luminescence tomography based on depth compensation and group sparsity prior. J. Biomed. Opt. 2020, 25, 016004. [Google Scholar] [CrossRef]
  17. Chen, D.; Zhu, S.; Cao, X.; Zhao, F.; Liang, J. X-ray luminescence computed tomography imaging based on X-ray distribution model and adaptively split Bregman method. Biomed. Opt. Express 2015, 6, 2649–2663. [Google Scholar] [CrossRef]
  18. Li, H.; Yang, X.; Yang, S.; Wang, D.; Jeon, G. Transformer With Double Enhancement for Low-Dose CT Denoising. IEEE J. Biomed. Health Inform. 2022, 27, 4660–4671. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, H.; Zhao, X.; Liu, W.; Li, L.C.; Ma, J.; Guo, L. Texture-aware dual domain mapping model for low-dose CT reconstruction. Med. Phys. 2022, 49, 3860–3873. [Google Scholar] [CrossRef]
  20. Wang, G.; Ye, J.C.; Man, B.D. Deep learning for tomographic image reconstruction. Nat. Mach. Intell. 2020, 2, 737–748. [Google Scholar] [CrossRef]
  21. Shan, H.; Zhang, Y.; Yang, Q.; Kruger, U.; Kalra, M.K.; Sun, L.; Cong, W.; Wang, G. 3D Convolutional Encoder-Decoder Network for Low-Dose CT via Transfer Learning from a 2D Trained Network. IEEE Trans. Med. Imaging 2018, 37, 1522. [Google Scholar] [CrossRef] [PubMed]
  22. Gao, Y.; Wang, K.; An, Y.; Jiang, S.; Meng, H.; Tian, J. Nonmodel-based bioluminescence tomography using a machine-learning reconstruction strategy. Optica 2018, 11, 1451–1454. [Google Scholar] [CrossRef]
  23. Zhang, X.; Cao, X.; Zhang, P.; Song, F.; Zhang, J.; Zhang, L.; Zhang, G. Self-training Strategy based on Finite Element Method for Adaptive Bioluminescence Tomography Reconstruction. IEEE Trans. Med. Imaging 2022, 41, 2629–2643. [Google Scholar] [CrossRef]
  24. Zhang, P.; Fan, G.; Xing, T.; Song, F.; Zhang, G. UHR-DeepFMT: Ultra-High Spatial Resolution Reconstruction of Fluorescence Molecular Tomography Based on 3-D Fusion Dual-Sampling Deep Neural Network. IEEE Trans. Med. Imaging 2021, 40, 3217–3228. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, P.; Ma, C.; Song, F.; Zhang, T.; Sun, Y.; Feng, Y.; He, Y.; Liu, F.; Wang, D.; Zhang, G. D2-RecST: Dual-domain Joint Reconstruction Strategy for Fluorescence Molecular Tomography Based on Image Domain and Perception Domain. Comput. Methods Programs Biomed. 2022, 229, 107293. [Google Scholar] [CrossRef]
  26. Gao, P.; Pu, H.; Rong, J.; Zhang, W.; Liu, T.; Liu, W.; Zhang, Y.; Lu, H. Resolving adjacent nanophosphors of different concentrations by excitation-based cone-beam X-ray luminescence tomography. Biomed. Opt. Express 2017, 8, 3952–3965. [Google Scholar] [CrossRef]
  27. Carpenter, C.; Pratx, G.; Sun, C.; Xing, L. Limited-angle x-ray luminescence tomography: Methodology and feasibility study. Phys. Med. Biol. 2011, 56, 3487. [Google Scholar] [CrossRef]
  28. Schweiger, M.; Arridge, S.; Hiraoka, M.; Delpy, D. The finite element method for the propagation of light in scattering media: Boundary and source conditions. Med. Phys. 1995, 22, 1779–1792. [Google Scholar] [CrossRef]
  29. Lv, Y.; Tian, J.; Cong, W.; Wang, G.; Luo, J.; Yang, W.; Li, H. A multilevel adaptive finite element algorithm for bioluminescence tomography. Opt. Express 2006, 14, 8211–8223. [Google Scholar] [CrossRef]
  30. Shi, H.; Wang, L.; Zheng, N.; Hua, G.; Tang, W. Loss functions for pose guided person image generation. Pattern Recognit. 2022, 122, 108351. [Google Scholar] [CrossRef]
  31. Tsui, B.M.W.; Lalush, D.S. A generalized Gibbs prior for maximum a posteriori reconstruction in SPECT. Phys. Med. Biol. 1993, 38, 729–741. [Google Scholar]
  32. Herman, G.T.; De Pierro, A.R.; Gai, N. On methods for maximum a posteriori image reconstruction with a normal prior. J. Vis. Commun. Image Represent. 1992, 3, 316–324. [Google Scholar] [CrossRef]
  33. Gao, P.; Rong, J.; Pu, H.; Liu, T.; Zhang, W.; Zhang, X.; Lu, H. Sparse view cone beam X-ray luminescence tomography based on truncated singular value decomposition. Opt. Express 2018, 26, 23233–23250. [Google Scholar] [CrossRef]
  34. Dong, D.; Zhu, S.; Qin, C.; Kumar, V.; Stein, J.V.; Oehler, S.; Savakis, C.; Tian, J.; Ripoll, J. Automated recovery of the center of rotation in optical projection tomography in the presence of scattering. IEEE J. Biomed. Health Inform. 2013, 17, 198–204. [Google Scholar] [CrossRef]
  35. Zhu, S.; Tian, J.; Yan, G.; Qin, C.; Feng, J. Cone beam micro-CT system for small animal imaging and performance evaluation. J. Biomed. Imaging 2009, 2009, 16. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic illustration of the DeepCB-XLCT network architecture: the 3D deep encoder–decoder (3D-En–Decoder) network has a 3D-Encoder network and a 3D-Decoder network. The 3D encoder consists of several convolution layers, followed by batch norm, ReLU activation function, and pooling. The 3D-Decoder has several upsampling layers followed by convolution, batch norm, and ReLU activation function. There are two full connection layers between the 3D-Encoder and the 3D-Decoder.
Figure 1. Schematic illustration of the DeepCB-XLCT network architecture: the 3D deep encoder–decoder (3D-En–Decoder) network has a 3D-Encoder network and a 3D-Decoder network. The 3D encoder consists of several convolution layers, followed by batch norm, ReLU activation function, and pooling. The 3D-Decoder has several upsampling layers followed by convolution, batch norm, and ReLU activation function. There are two full connection layers between the 3D-Encoder and the 3D-Decoder.
Bioengineering 11 00874 g001
Figure 2. The cylinder phantom used in the training data simulation with two targets and three targets (a) two targets, (b) three targets, the region within the red circle is the entire reconstruction region, the region within the green rectangular box is the target region.
Figure 2. The cylinder phantom used in the training data simulation with two targets and three targets (a) two targets, (b) three targets, the region within the red circle is the entire reconstruction region, the region within the green rectangular box is the target region.
Bioengineering 11 00874 g002
Figure 3. The cylinder phantom used in simulation studies. Edge-to-edge distance between the two targets: (a) 1 mm, (b) 1.5 mm, (c) 2 mm, the red circle is the boundary of the phantom, the green line is the centerline of the phantom.
Figure 3. The cylinder phantom used in simulation studies. Edge-to-edge distance between the two targets: (a) 1 mm, (b) 1.5 mm, (c) 2 mm, the red circle is the boundary of the phantom, the green line is the centerline of the phantom.
Bioengineering 11 00874 g003
Figure 4. The schematic diagram of the CB-XLCT system.
Figure 4. The schematic diagram of the CB-XLCT system.
Bioengineering 11 00874 g004
Figure 5. Configuration for the two-target phantom experiment. The concentrations of the two targets were both 50 mg/mL, and the edge-to-edge distances of the two targets were 2.3 mm, 1.7 mm, and 1.0 mm.
Figure 5. Configuration for the two-target phantom experiment. The concentrations of the two targets were both 50 mg/mL, and the edge-to-edge distances of the two targets were 2.3 mm, 1.7 mm, and 1.0 mm.
Bioengineering 11 00874 g005
Figure 6. Tomographic images of the targets positioned at varying distances were reconstructed. The first through fifth columns present the results obtained using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively. The first to third rows correspond to reconstructions with edge-to-edge distances of 2.0 mm, 1.5 mm, and 1.0 mm, respectively. The red circle is the boundary of the phantom, the green line is the centerline of the reconstruction image and the yellow circle represents the true target region.
Figure 6. Tomographic images of the targets positioned at varying distances were reconstructed. The first through fifth columns present the results obtained using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively. The first to third rows correspond to reconstructions with edge-to-edge distances of 2.0 mm, 1.5 mm, and 1.0 mm, respectively. The red circle is the boundary of the phantom, the green line is the centerline of the reconstruction image and the yellow circle represents the true target region.
Bioengineering 11 00874 g006
Figure 7. Profiles along the green dotted line in Figure 5 with the edge-to-edge distances of 1.0 mm. (ae) The profiles of reconstruction results achieved with ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively.
Figure 7. Profiles along the green dotted line in Figure 5 with the edge-to-edge distances of 1.0 mm. (ae) The profiles of reconstruction results achieved with ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively.
Bioengineering 11 00874 g007
Figure 8. Tomographic images of the targets positioned at different distances with different noise levels were reconstructed based on the proposed DeepCB-XLCT network. The red circle is the boundary of the phantom and the yellow circle represents the true target region. The first to fourth columns are the reconstruction results with the SNRs of 30 dB, 25 dB, 20 dB, and 15 dB, respectively. The first to third row are the results reconstructed with the edge-to-edge distances of 2.0 mm, 1.5 mm, and 1.0 mm, respectively.
Figure 8. Tomographic images of the targets positioned at different distances with different noise levels were reconstructed based on the proposed DeepCB-XLCT network. The red circle is the boundary of the phantom and the yellow circle represents the true target region. The first to fourth columns are the reconstruction results with the SNRs of 30 dB, 25 dB, 20 dB, and 15 dB, respectively. The first to third row are the results reconstructed with the edge-to-edge distances of 2.0 mm, 1.5 mm, and 1.0 mm, respectively.
Bioengineering 11 00874 g008
Figure 9. Tomographic images of the targets positioned at different distances were reconstructed in the ablation experiment. The red circle is the boundary of the phantom and the yellow circle represents the true target region. The first to fourth rows are the results obtained based on the proposed DeepCB-XLCT network without skip connection and ROILoss, the proposed DeepCB-XLCT network without skip connection, the proposed DeepCB-XLCT network without ROILoss, and the proposed DeepCB-XLCT network, respectively. The first to third columns are the results reconstructed with the edge-to-edge distances of 2.0 mm, 1.5 mm, and 1.0 mm, respectively.
Figure 9. Tomographic images of the targets positioned at different distances were reconstructed in the ablation experiment. The red circle is the boundary of the phantom and the yellow circle represents the true target region. The first to fourth rows are the results obtained based on the proposed DeepCB-XLCT network without skip connection and ROILoss, the proposed DeepCB-XLCT network without skip connection, the proposed DeepCB-XLCT network without ROILoss, and the proposed DeepCB-XLCT network, respectively. The first to third columns are the results reconstructed with the edge-to-edge distances of 2.0 mm, 1.5 mm, and 1.0 mm, respectively.
Bioengineering 11 00874 g009
Figure 10. Tomographic images of three targets were reconstructed using various methods. The red circle is the boundary of the phantom and the yellow circle represents the true target region. The first column depicts the true locations of the targets. The second to sixth columns present results obtained using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively. The first to third rows show the reconstructions for the three targets positioned at different locations.
Figure 10. Tomographic images of three targets were reconstructed using various methods. The red circle is the boundary of the phantom and the yellow circle represents the true target region. The first column depicts the true locations of the targets. The second to sixth columns present results obtained using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively. The first to third rows show the reconstructions for the three targets positioned at different locations.
Bioengineering 11 00874 g010
Figure 11. Tomographic images of the targets positioned at varying distances were reconstructed based on different algorithms for phantom experiments. The first to third rows display the fused XLCT/CT tomographic images for the EEDs of 2.3, 1.7, and 1.0 mm, respectively. The reconstructions obtained from ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network are presented in the first through fifth columns, respectively.
Figure 11. Tomographic images of the targets positioned at varying distances were reconstructed based on different algorithms for phantom experiments. The first to third rows display the fused XLCT/CT tomographic images for the EEDs of 2.3, 1.7, and 1.0 mm, respectively. The reconstructions obtained from ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network are presented in the first through fifth columns, respectively.
Bioengineering 11 00874 g011
Figure 12. Tomographic images reconstructed based on different algorithms for mice experiments. (a) The CT projection image, the region between the two red lines was used for reconstruction; (b) the reconstructed CT slice at the height of the green line of (a); (cg) reconstructions obtained using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively; (hl) the corresponding fusion results of XLCT and CT.
Figure 12. Tomographic images reconstructed based on different algorithms for mice experiments. (a) The CT projection image, the region between the two red lines was used for reconstruction; (b) the reconstructed CT slice at the height of the green line of (a); (cg) reconstructions obtained using ADFISTA, MAP, T-FISTA, ADMLEM, and the proposed DeepCB-XLCT network, respectively; (hl) the corresponding fusion results of XLCT and CT.
Bioengineering 11 00874 g012
Table 1. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs obtained from numerical simulations.
Table 1. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs obtained from numerical simulations.
Methods2 mm1.5 mm1.0 mm
CNRDICECNRDICECNRDICE
ADFISTA2.470.442.580.412.660.42
MAP2.830.383.010.373.160.38
T-FISTA1.060.262.390.472.310.45
ADMLEM2.190.442.390.462.560.48
Proposed method4.320.894.690.974.490.96
Table 2. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs under different noise levels from numerical simulations.
Table 2. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs under different noise levels from numerical simulations.
Noise Level2.0 mm1.5 mm1.0 mm
CNRDICECNRDICECNRDICE
30 dB4.320.894.690.974.490.96
25 dB4.080.814.490.953.940.38
20 dB3.600.784.040.903.760.45
15 dB2.930.693.390.803.240.81
Table 3. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs in the ablation experiment.
Table 3. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs in the ablation experiment.
Methods2.0 mm1.5 mm1.0 mm
CNRDICECNRDICECNRDICE
Without skip connection and ROILoss1.990.472.190.582.200.57
Without skip connection4.050.864.470.964.190.94
Without ROILoss3.310.853.690.903.550.91
Proposed method4.320.894.690.974.490.96
Table 4. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs in phantom experiments.
Table 4. CNR and DICE metrics for the reconstruction results of two targets with varying EEDs in phantom experiments.
Methods2.3 mm1.7 mm1.0 mm
CNRDICECNRDICECNRDICE
ADFISTA3.940.564.030.484.760.41
MAP3.620.494.880.564.360.45
T-FISTA3.410.577.640.733.860.53
ADMLEM7.640.757.920.784.860.58
Proposed method7.570.5611.090.7912.060.71
Table 5. CNR and DICE metrics for the reconstruction results of two targets in mouse experiments.
Table 5. CNR and DICE metrics for the reconstruction results of two targets in mouse experiments.
ADFISTAMAPT-FISTAADMLEMProposed Method
CNR3.151.612.553.634.33
DICE0.320.160.240.350.60
Table 6. Reconstruction time with different methods in mice experiments.
Table 6. Reconstruction time with different methods in mice experiments.
MethodsADFISTAMAPT-FISTAADMLEMProposed Method
Reconstruction time/s 14.948.148.615.81
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, T.; Huang, S.; Li, R.; Gao, P.; Li, W.; Lu, H.; Song, Y.; Rong, J. Dual and Multi-Target Cone-Beam X-ray Luminescence Computed Tomography Based on the DeepCB-XLCT Network. Bioengineering 2024, 11, 874. https://doi.org/10.3390/bioengineering11090874

AMA Style

Liu T, Huang S, Li R, Gao P, Li W, Lu H, Song Y, Rong J. Dual and Multi-Target Cone-Beam X-ray Luminescence Computed Tomography Based on the DeepCB-XLCT Network. Bioengineering. 2024; 11(9):874. https://doi.org/10.3390/bioengineering11090874

Chicago/Turabian Style

Liu, Tianshuai, Shien Huang, Ruijing Li, Peng Gao, Wangyang Li, Hongbing Lu, Yonghong Song, and Junyan Rong. 2024. "Dual and Multi-Target Cone-Beam X-ray Luminescence Computed Tomography Based on the DeepCB-XLCT Network" Bioengineering 11, no. 9: 874. https://doi.org/10.3390/bioengineering11090874

APA Style

Liu, T., Huang, S., Li, R., Gao, P., Li, W., Lu, H., Song, Y., & Rong, J. (2024). Dual and Multi-Target Cone-Beam X-ray Luminescence Computed Tomography Based on the DeepCB-XLCT Network. Bioengineering, 11(9), 874. https://doi.org/10.3390/bioengineering11090874

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop