Principal Component Analysis Based Dynamic Cone Beam X-Ray Luminescence Computed Tomography: A Feasibility Study

SODL-IR-FISTA: sparse online dictionary learning with iterative reduction FISTA for cone-beam X-ray luminescence computed tomography

Cone beam X-ray luminescence computed tomography (CB-XLCT) is an emerging imaging technique with potential for early 3D tumor detection. However, the reconstruction challenge due to low light absorption and high scattering in tissues makes it a difficult inverse problem. In this study, the online dictionary learning (ODL) method, combined with iterative reduction FISTA (IR-FISTA), has been utilized to achieve high-quality reconstruction. Our method integrates IR-FISTA for efficient and accurate sparse coding, followed by an online stochastic approximation for dictionary updates, effectively capturing the sparse features inherent to the problem. Additionally, a re-sparse step is introduced to enhance the sparsity of the solution, making it better suited for CB-XLCT reconstruction. Numerical simulations and in vivo experiments were conducted to assess the performance of the method. The SODL-IR-FISTA achieved the smallest location error of 0.325 mm in in vivo experiments, which is 58% and 45% of the IVTCG-L 1 (0.562 mm) and OMP-L 0 (0.721 mm), respectively. Additionally, it has the highest DICE similarity coefficient, which is 0.748. The results demonstrate that our approach outperforms traditional methods in terms of localization precision, shape restoration, robustness, and practicality in live subjects.

Superfast Scan of Focused X-Ray Luminescence Computed Tomography Imaging

X-ray luminescence computed tomography (XLCT) is a hybrid molecular imaging modality having the high spatial resolution of x-ray imaging and high measurement sensitivity of optical imaging. Narrow x-ray beam based XLCT imaging has shown promise for high spatial resolution imaging of luminescent targets in deep tissues, but the slow acquisition speed limits its applications. In this work, we have introduced a superfast XLCT scan scheme based on the photon counter detector and a fly-scanning method. The new scan scheme is compared with three other scan methods. We have also designed and built a single-pixel x-ray detector to detect object boundaries automatically. With the detector, we can perform the parallel beam CT imaging with the XLCT imaging simultaneously. We have built the prototype XLCT imaging system to verify the proposed scan scheme. A phantom embedded with a set of four side-by-side cylindrical targets was scanned. With the proposed superfast scan scheme, we have achieved 43 seconds per transverse scan, which is 28.6 times faster than before with slightly better XLCT image quality. The superfast scan allows us to perform 3D pencil beam XLCT imaging in the future.

Sparse reconstruction based on dictionary learning and group structure strategy for cone-beam X-ray luminescence computed tomography

As a dual-modal imaging technology that has emerged in recent years, cone-beam X-ray luminescence computed tomography (CB-XLCT) has exhibited promise as a tool for the early three-dimensional detection of tumors in small animals. However, due to the challenges imposed by the low absorption and high scattering of light in tissues, the CB-XLCT reconstruction problem is a severely ill-conditioned inverse problem, rendering it difficult to obtain satisfactory reconstruction results. In this study, a strategy that utilizes dictionary learning and group structure (DLGS) is proposed to achieve satisfactory CB-XLCT reconstruction performance. The group structure is employed to account for the clustering of nanophosphors in specific regions within the organism, which can enhance the interrelation of elements in the same group. Furthermore, the dictionary learning strategy is implemented to effectively capture sparse features. The performance of the proposed method was evaluated through numerical simulations and in vivo experiments. The experimental results demonstrate that the proposed method achieves superior reconstruction performance in terms of location accuracy, target shape, robustness, dual-source resolution, and in vivo practicability.

OPK_SNCA: Optimized prior knowledge via sparse non-convex approach for cone-beam X-ray luminescence computed tomography imaging

BACKGROUND

The development of Cone-beam X-ray luminescence computed tomography (CB-XLCT) has allowed the quantitative in-depth biological imaging, but with a greatly ill-posed and ill-conditioned inverse problem. Although the predefined permissible source region (PSR) is a widely used way to alleviate the problem for CB-XLCT imaging, how to obtain the accurate PSR is still a challenge for the process of inverse reconstruction.

METHODS

We proposed an optimized prior knowledge via a sparse non-convex approach (OPK_SNCA) for CB-XLCT imaging. Firstly, non-convex L_p-norm optimization model was employed for copying with the inverse problem, and an iteratively reweighted split augmented lagrangian shrinkage algorithm was developed to obtain a group of sparse solutions based on different non-convex p values. Secondly, a series of permissible regions (PRs) with different discretized mesh was further achieved, and the intersection operation was implemented on the group of PRs to get a reasonable PSR. After that, the final PSR was adopted as an optimized prior knowledge to enhance the reconstruction quality of inverse reconstruction.

RESULTS

Both simulation experiments and in vivo experiment were performed to evaluate the efficiency and robustness of the proposed method.

CONCLUSIONS

The experimental results demonstrated that our proposed method could significantly improve the imaging quality of the distribution of X-ray-excitable nanophosphors for CB-XLCT.

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.

Development of tomographic reconstruction for three-dimensional optical imaging: From the inversion of light propagation to artificial intelligence

Optical molecular tomography (OMT) is an imaging modality which uses an optical signal, especially near-infrared light, to reconstruct the three-dimensional information of the light source in biological tissue. With the advantages of being low-cost, noninvasive and having high sensitivity, OMT has been applied in preclinical and clinical research. However, due to its serious ill-posedness and ill-condition, the solution of OMT requires heavy data analysis and the reconstruction quality is limited. Recently, the artificial intelligence (commonly known as AI)-based methods have been proposed to provide a different tool to solve the OMT problem. In this paper, we review the progress on OMT algorithms, from conventional methods to AI-based methods, and we also give a prospective towards future developments in this domain.

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.