Cone-beam x-ray luminescence computed tomography based on x-ray absorption dosage

Cone-beam x-ray luminescence computed tomography (CB-XLCT) prototype development and performance evaluation

This study developed a prototype for a rotational cone-beam x-ray luminescence computed tomography (CB-XLCT) system, considering its potential application in pre-clinical theranostic imaging. A geometric calibration method applicable to both imaging chains (XL and CT) was also developed to enhance image quality. The results of systematic performance evaluations were presented to assess the feasibility of commercializing XLCT technology. Monte Carlo GATE simulation was performed to determine the optimal imaging conditions for nanophosphor particles (NPs) irradiated by 70 kV x-rays. We acquired a low-dose transmission x-ray tube and designed a prone positioning platform and a rotating gantry, using mice as targets from commercial small animalμ-CT systems. We then employed the image cross-correlation (ICC) automatic geometric calibration method to calibrate XL and CT images. The performance of the system was evaluated through a series of phantom experiments with a linearity of 0.99, and the contrast-to-noise ratio (CNR) between hydroxyl-apatite (HA) and based epoxy resin is 19.5. The XL images of the CB-XLCT prototype achieved a Dice similarity coefficient (DICE) of 0.149 for a distance of 1 mm between the two light sources. Finally, the final XLCT imaging results were demonstrated using the Letter phantoms with NPs. In summary, the CB-XLCT prototype developed in this study showed the potential to achieve high-quality imaging with acceptable radiation doses for small animals. The performance of CT images was comparable to current commercial machines, while the XL images exhibited promising results in phantom imaging, but further efforts are needed for biomedical applications.

A Finite Element Mesh Regrouping Strategy-Based Hybrid Light Transport Model for Enhancing the Efficiency and Accuracy of XLCT

X-ray luminescence computed tomography (XLCT) is an emerging hybrid imaging modality in optical molecular imaging, which has attracted more attention and has been widely studied. In XLCT, the accuracy and operational efficiency of an optical transmission model play a decisive role in the rapid and accurate reconstruction of light sources. For simulation of optical transmission characteristics in XLCT, considering the limitations of the diffusion equation (DE) and the time and memory costs of simplified spherical harmonic approximation equation (SP_N), a hybrid light transport model needs to be built. DE and SP_N models are first-order and higher-order approximations of RTE, respectively. Due to the discontinuity of the regions using the DE and SP_N models and the inconsistencies of the system matrix dimensions constructed by the two models in the solving process, the system matrix construction of a hybrid light transmission model is a problem to be solved. We provided a new finite element mesh regrouping strategy-based hybrid light transport model for XLCT. Firstly, based on the finite element mesh regrouping strategy, two separate meshes can be obtained. Thus, for DE and SP_N models, the system matrixes and source weight matrixes can be calculated separately in two respective mesh systems. Meanwhile, some parallel computation strategy can be combined with finite element mesh regrouping strategy to further save the system matrix calculation time. Then, the two system matrixes with different dimensions were coupled though repeated nodes were processed according to the hybrid boundary conditions, the two meshes were combined into a regrouping mesh, and the hybrid optical transmission model was established. In addition, the proposed method can reduce the computational memory consumption than the previously proposed hybrid light transport model achieving good balance between computational accuracy and efficiency. The forward numerical simulation results showed that the proposed method had better transmission accuracy and achieved a balance between efficiency and accuracy. The reverse simulation results showed that the proposed method had superior location accuracy, morphological recovery capability, and image contrast capability in source reconstruction. In-vivo experiments verified the practicability and effectiveness of the proposed method.

Restarted primal-dual Newton conjugate gradient method for enhanced spatial resolution of reconstructed cone-beam x-ray luminescence computed tomography images

Cone-beam x-ray luminescence computed tomography (CB-XLCT) has been proposed as a promising imaging tool, which enables three-dimensional imaging of the distribution of nanophosphors (NPs) in small animals. However, the reconstruction performance is usually unsatisfactory in terms of spatial resolution due to the ill-posedness of the CB-XLCT inverse problem. To alleviate this problem and to achieve high spatial resolution, a reconstruction method consisting of inner and outer iterations based on a restarted strategy is proposed. In this method, the primal-dual Newton conjugate gradient method (pdNCG) is adopted in the inner iterations to get fast reconstruction, which is used for resetting the permission region and increasing the convergence speed of the outer iteration. To assess the performance of the method, both numerical simulation and physical phantom experiments were conducted with a CB-XLCT system. The results demonstrate that compared with conventional reconstruction methods, the proposed re-pdNCG method can accurately and efficiently resolve the adjacent NPs with the least relative error.

Regularized reconstruction based on joint L1 and total variation for sparse-view cone-beam X-ray luminescence computed tomography

As an emerging hybrid imaging modality, cone-beam X-ray luminescence computed tomography (CB-XLCT) has been proposed based on the development of X-ray excitable nanoparticles. Owing to the high degree of absorption and scattering of light through tissues, the CB-XLCT inverse problem is inherently ill-conditioned. Appropriate priors or regularizations are needed to facilitate reconstruction and to restrict the search space to a specific solution set. Typically, the goal of CB-XLCT reconstruction is to get the distributions of nanophosphors in the imaging object. Considering that the distributions of nanophosphors inside bodies preferentially accumulate in specific areas of interest, the reconstruction of XLCT images is usually sparse with some locally smoothed high-intensity regions. Therefore, a combination of the L₁ and total variation regularization is designed to improve the imaging quality of CB-XLCT in this study. The L₁ regularization is used for enforcing the sparsity of the reconstructed images and the total variation regularization is used for maintaining the local smoothness of the reconstructed image. The implementation of this method can be divided into two parts. First, the reconstruction image was reconstructed based on the fast iterative shrinkage-thresholding (FISTA) algorithm, then the reconstruction image was minimized by the gradient descent method. Numerical simulations and phantom experiments indicate that compared with the traditional ART, ADAPTIK and FISTA methods, the proposed method demonstrates its advantage in improving spatial resolution and reducing imaging time.