Light&fast generative adversarial network for high-fidelity CT image synthesis of liver tumor.

BACKGROUND AND OBJECTIVE: Hepatocellular carcinoma is a common disease with high mortality. Through deep learning methods to analyze HCC CT, the screening classification and prognosis model of HCC can be established, which further promotes the development of computer-aided diagnosis and treatment in the treatment of HCC. However, there are significant challenges in the actual establishment of HCC auxiliary diagnosis model due to data imbalance, especially for rare subtypes of HCC and underrepresented demographic groups. This study proposes a GAN model aimed at overcoming these obstacles and improving the accuracy of HCC auxiliary diagnosis. METHODS: In order to generate liver and tumor images close to the real distribution. Firstly, we construct a new gradient transfer sampling module to improve the lack of texture details and excessive gradient transfer parameters of the deep model; Secondly, we construct an attention module with spatial and cross-channel feature extraction ability to improve the discriminator's ability to distinguish images; Finally, we design a new loss function for liver tumor imaging features to constrain the model to approach the real tumor features in iterations. RESULTS: In qualitative analysis, the images synthetic by our method closely resemble the real images in liver parenchyma, blood vessels, tumors, and other parts. In quantitative analysis, the optimal results of FID, PSNR, and SSIM are 75.73, 22.77, and 0.74, respectively. Furthermore, our experiments establish classification models for imbalanced data and enhanced data, resulting in an increase in accuracy rate by 21%-34%, an increase in AUC by 0.29 - 0.33, and an increase in specificity to 0.89. CONCLUSION: Our solution provides a variety of training data sources with low cost and high efficiency for the establishment of classification or prognostic models for imbalanced data.

Group sparse-based Taylor expansion method for liver pharmacokinetic parameters imaging of dynamic fluorescence molecular tomography.

Objective.Pharmacokinetic parametric images obtained through dynamic fluorescence molecular tomography (DFMT) has ability of capturing dynamic changes in fluorescence concentration, thereby providing three-dimensional metabolic information for applications in biological research and drug development. However, data processing of DFMT is time-consuming, involves a vast amount of data, and the problem itself is ill-posed, which significantly limits the application of pharmacokinetic parametric images reconstruction. In this study, group sparse-based Taylor expansion method is proposed to address these problems.Approach.Firstly, Taylor expansion framework is introduced to reduce time and computational cost. Secondly, group sparsity based on structural prior is introduced to improve reconstruction accuracy. Thirdly, alternating iterative solution based on accelerated gradient descent algorithm is introduced to solve the problem.Main results.Numerical simulation andin vivoexperimental results demonstrate that, in comparison to existing methods, the proposed approach significantly enhances reconstruction speed without a degradation of quality, particularly when confronted with background fluorescence interference from other organs.Significance.Our research greatly reduces time and computational cost, providing strong support for real-time monitoring of liver metabolism.

FSMN-Net: a free space matching network based on manifold convolution for optical molecular tomography.

Optical molecular tomography (OMT) can monitor glioblastomas in small animals non-invasively. Although deep learning (DL) methods have made remarkable achievements in this field, improving its generalization against diverse reconstruction systems remains a formidable challenge. In this Letter, a free space matching network (FSMN-Net) was presented to overcome the parameter mismatch problem in different reconstruction systems. Specifically, a novel, to the best of our knowledge, manifold convolution operator was designed by considering the mathematical model of OMT as a space matching process. Based on the dynamic domain expansion concept, an end-to-end fully convolutional codec further integrates this operator to realize robust reconstruction with voxel-level accuracy. The results of numerical simulations and in vivo experiments demonstrate that the FSMN-Net can stably generate high-resolution reconstruction volumetric images under different reconstruction systems.

Deep learning-based automatic scoring models for the disease activity of rheumatoid arthritis based on multimodal ultrasound images.

OBJECTIVES: We aimed to investigate the value of deep learning (DL) models based on multimodal ultrasonographic (US) images to quantify RA activity. METHODS: Static greyscale (SGS), dynamic greyscale (DGS), static power Doppler (SPD) and dynamic power Doppler (DPD) US images were collected and evaluated by two expert radiologists according to the EULAR-OMERACT Synovitis Scoring system. Four DL models were developed based on the ResNet-type structure, evaluated on two separate test cohorts, and finally compared with the performance of 12 radiologists with different levels of experience. RESULTS: In total, 1244 images were used for the model training, and 152 and 354 for testing (cohort 1 and 2, respectively). The best-performing models for the scores of 0/1/2/3 were the DPD, SGS, DGS and SPD models, respectively (Area Under the receiver operating characteristic Curve [AUC] = 0.87/0.95/0.74/0.95; no significant differences). All the DL models provided results comparable to the experienced radiologists on a per-image basis (intraclass correlation coefficient: 0.239-0.756, P < 0.05). The SPD model performed better than the SGS one on test cohort 1 (score of 0/2/3: AUC = 0.82/0.67/0.95 vs 0.66/0.66/0.75, respectively) and test cohort 2 (score of 0: AUC = 0.89 vs 0.81). The dynamic DL models performed better than the static ones in most of the scoring processes and were more accurate than the most of senior radiologists, especially the DPD model. CONCLUSION: DL models based on multimodal US images allow a quantitative and objective assessment of RA activity. Dynamic DL models in particular have potential value in assisting radiologists to improve the accuracy of RA US-based grading.

Structure-fused deep 3D hierarchical network: A bioluminescence tomography scheme for different imaging objects.

Monte Carlo eXtreme (MCX) method has a unique advantage for deep neural network based bioluminescence tomography (BLT) reconstruction. However, this method ignores the distribution of sources energy and relies on the determined tissue structure. In this paper, a deep 3D hierarchical reconstruction network for BLT was proposed where the inputs were divided into two parts -- bioluminescence image (BLI) and anatomy of the imaged object by CT. Firstly, a parallel encoder is used to feature the original BLI & CT slices and integrate their features to distinguish the different tissue structure of imaging objects; Secondly, GRU is used to fit the spatial information of different slices and convert it into 3D features; Finally, the 3D features are decoded to the spacial and energy information of source by a symmetrical decoding structure. Our research suggested that this method can effectively compute the radiation intensity and the spatial distribution of the source for different imaging object.

One-dimensional convolutional neural network for Jacobian in Diffuse Optical Tomography.

Non-linear least square minimization algorithms are often employed to solve Diffuse Optical Tomography (DOT) inverse problem. However, it is time-consuming to calculate the Jacobian matrix. This work has proposed a data-driven neural network method to improve computational efficiency. The singular value decomposition is employed to compute the updated Jacobian and a mapping from boundary measurements to the singular values based on a convolutional neural network (CNN) is learned to obtain the singular values. The method is validated with 3D numerical simulation data. We have demonstrated that the approach can save computation time compared to Adjoint method, and reconstructed absorption coefficient close to Adjoint method.Clinical Relevance- These results are not focused on clinical relevance currently, but in the future may be helpful to accelerant DOT reconstruction in clinic.

A multitask deep learning radiomics model for predicting the macrotrabecular-massive subtype and prognosis of hepatocellular carcinoma after hepatic arterial infusion chemotherapy.

BACKGROUND: The macrotrabecular-massive (MTM) is a special subtype of hepatocellular carcinoma (HCC), which has commonly a dismal prognosis. This study aimed to develop a multitask deep learning radiomics (MDLR) model for predicting MTM and HCC patients' prognosis after hepatic arterial infusion chemotherapy (HAIC). METHODS: From June 2018 to March 2020, 158 eligible patients with HCC who underwent surgery were retrospectively enrolled in MTM related cohorts, and 752 HCC patients who underwent HAIC were included in HAIC related cohorts during the same period. DLR features were extracted from dual-phase (arterial phase and venous phase) contrast-enhanced computed tomography (CECT) of the entire liver region. Then, an MDLR model was used for the simultaneous prediction of the MTM subtype and patient prognosis after HAIC. The MDLR model for prognostic risk stratification incorporated DLR signatures, clinical variables and MTM subtype. FINDINGS: The predictive performance of the DLR model for the MTM subtype was 0.968 in the training cohort [TC], 0.912 in the internal test cohort [ITC] and 0.773 in the external test cohort [ETC], respectively. Multivariable analysis identified portal vein tumor thrombus (PVTT) (p = 0.012), HAIC response (p < 0.001), HAIC sessions (p < 0.001) and MTM subtype (p < 0.001) as indicators of poor prognosis. After incorporating DLR signatures, the MDLR model yielded the best performance among all models (AUC, 0.855 in the TC, 0.805 in the ITC and 0.792 in the ETC). With these variables, the MDLR model provided two risk strata for overall survival (OS) in the TC: low risk (5-year OS, 44.9%) and high risk (5-year OS, 4.9%). INTERPRETATION: A tool based on MDLR was developed to consider that the MTM is an important prognosis factor for HCC patients. MDLR showed outstanding performance for the prognostic risk stratification of HCC patients who underwent HAIC and may help physicians with therapeutic decision making and surveillance strategy selection in clinical practice.

Liver injury monitoring using dynamic fluorescence molecular tomography based on a time-energy difference strategy.

Dynamic fluorescence molecular tomography (DFMT) is a promising molecular imaging technique that offers the potential to monitor fast kinetic behaviors within small animals in three dimensions. Early monitoring of liver disease requires the ability to distinguish and analyze normal and injured liver tissues. However, the inherent ill-posed nature of the problem and energy signal interference between the normal and injured liver regions limit the practical application of liver injury monitoring. In this study, we propose a novel strategy based on time and energy, leveraging the temporal correlation in fluorescence molecular imaging (FMI) sequences and the metabolic differences between normal and injured liver tissue. Additionally, considering fluorescence signal distribution disparity between the injured and normal regions, we designed a universal Golden Ratio Primal-Dual Algorithm (GRPDA) to reconstruct both the normal and injured liver regions. Numerical simulation and in vivo experiment results demonstrate that the proposed strategy can effectively avoid signal interference between liver and liver injury energy and lead to significant improvements in morphology recovery and positioning accuracy compared to existing approaches. Our research presents a new perspective on distinguishing normal and injured liver tissues for early liver injury monitoring.

Neighbor-based adaptive sparsity orthogonal least square for fluorescence molecular tomography.

Significance: Fluorescence molecular tomography (FMT) is a promising imaging modality, which has played a key role in disease progression and treatment response. However, the quality of FMT reconstruction is limited by the strong scattering and inadequate surface measurements, which makes it a highly ill-posed problem. Improving the quality of FMT reconstruction is crucial to meet the actual clinical application requirements. Aim: We propose an algorithm, neighbor-based adaptive sparsity orthogonal least square (NASOLS), to improve the quality of FMT reconstruction. Approach: The proposed NASOLS does not require sparsity prior information and is designed to efficiently establish a support set using a neighbor expansion strategy based on the orthogonal least squares algorithm. The performance of the algorithm was tested through numerical simulations, physical phantom experiments, and small animal experiments. Results: The results of the experiments demonstrated that the NASOLS significantly improves the reconstruction of images according to indicators, especially for double-target reconstruction. Conclusion: NASOLS can recover the fluorescence target with a good location error according to simulation experiments, phantom experiments and small mice experiments. This method is suitable for sparsity target reconstruction, and it would be applied to early detection of tumors.

3D-deep optical learning: a multimodal and multitask reconstruction framework for optical molecular tomography.

Optical molecular tomography (OMT) is an emerging imaging technique. To date, the poor universality of reconstruction algorithms based on deep learning for various imaged objects and optical probes limits the development and application of OMT. In this study, based on a new mapping representation, a multimodal and multitask reconstruction framework-3D deep optical learning (3DOL), was presented to overcome the limitations of OMT in universality by decomposing it into two tasks, optical field recovery and luminous source reconstruction. Specifically, slices of the original anatomy (provided by computed tomography) and boundary optical measurement of imaged objects serve as inputs of a recurrent convolutional neural network encoded parallel to extract multimodal features, and 2D information from a few axial planes within the samples is explicitly incorporated, which enables 3DOL to recognize different imaged objects. Subsequently, the optical field is recovered under the constraint of the object geometry, and then the luminous source is segmented by a learnable Laplace operator from the recovered optical field, which obtains stable and high-quality reconstruction results with extremely few parameters. This strategy enable 3DOL to better understand the relationship between the boundary optical measurement, optical field, and luminous source to improve 3DOL's ability to work in a wide range of spectra. The results of numerical simulations, physical phantoms, and in vivo experiments demonstrate that 3DOL is a compatible deep-learning approach to tomographic imaging diverse objects. Moreover, the fully trained 3DOL under specific wavelengths can be generalized to other spectra in the 620-900 nm NIR-I window.

Robust transformed l1 metric for fluorescence molecular tomography.

BACKGROUND AND OBJECTIVE: Fluorescence molecular tomography (FMT) is a non-invasive molecular imaging modality that can be used to observe the three-dimensional distribution of fluorescent probes in vivo. FMT is a promising imaging technique in clinical and preclinical research that has attracted significant attention. Numerous regularization based reconstruction algorithms have been proposed. However, traditional algorithms that use the squared l2-norm distance usually exaggerate the influence of noise and measurement and calculation errors, and their robustness cannot be guaranteed. METHODS: In this study, we propose a novel robust transformed l1 (TL1) metric that interpolates l0 and l1 norms through a nonnegative parameter α∈(0,+∞). The TL1 metric looks like the lp-norm with p∈(0,1). These are markedly different because TL1 metric has two properties, boundedness and Lipschitz-continuity, which make the TL1 criterion suitable distance metric, particularly for robustness, owing to its stronger noise suppression. Subsequently, we apply the proposed metric to FMT and build a robust model to reduce the influence of noise. The nonconvexity of the proposed model made direct optimization difficult, and a continuous optimization method was developed to solve the model. The problem was converted into a difference in convex programming problem for the TL1 metric (DCATL1), and the corresponding algorithm converged linearly. RESULTS: Various numerical simulations and in vivo bead-implanted mouse experiments were conducted to verify the performance of the proposed method. The experimental results show that the DCATL1 algorithm is more robust than the state-of-the-art approaches and achieves better source localization and morphology recovery. CONCLUSIONS: The in vivo experiments showed that DCATL1 can be used to visualize the distribution of fluorescent probes inside biological tissues and promote preclinical application in small animals, demonstrating the feasibility and effectiveness of the proposed method.

A Graph-guided Hybrid Regularization Method For Bioluminescence Tomography.

BACKGROUND AND OBJECTIVE: Bioluminescence tomography (BLT) is a powerful and sensitive imaging technique having great potential in preclinical application, such as tumor imaging, monitoring and therapy, etc. Regularization plays an important role in BLT reconstruction for considering the priori information to overcome the inherent ill-posedness of the inverse problem. Therefore, well-designed regularization term and sophisticated algorithm for solving the consequent optimization problem are key to improve the BLT quality. METHODS: To balance the sparsity, smoothness and morphological characteristics of the bioluminescence targets, we constructed a novel Graph-Guided Hybrid Regularization (GGHR) method by combining graph-guided penalty term with L1 and L2 norm regularizer. To solve the corresponding minimization problem with hybrid penalties, the dual decomposition and Nesterov's smoothing technique were adopted to decouple and transform the non-separable and non-smooth graph-guided penalty term into a differential smooth approximation form, which was solved by the fast iterative shrinkage thresholding algorithm. RESULTS: The performance of the proposed GGHR method was verified and evaluated through a series of simulation, phantom and in vivo experiments. The comparison results demonstrated that the GGHR method outperformed current mainstream reconstruction algorithms in spatial localization, morphology recovery and in vivo practicality. CONCLUSIONS: The proposed GGHR method is a robust and practicality reconstruction algorithm for further highlighting the positive effect of hybrid regularization on BLT applications.

Deep learning model based on contrast-enhanced ultrasound for predicting early recurrence after thermal ablation of colorectal cancer liver metastasis.

OBJECTIVES: To develop and validate a deep learning (DL) model based on quantitative analysis of contrast-enhanced ultrasound (CEUS) images that predicts early recurrence (ER) after thermal ablation (TA) of colorectal cancer liver metastasis (CRLM). METHODS: Between January 2010 and May 2019, a total of 207 consecutive patients with CRLM with 13,248 slice images at three dynamic phases who received CEUS within 2 weeks before TA were retrospectively enrolled in two centres (153 for the training cohort (TC), 32 for the internal test cohort (ITC), and 22 for the external test cohort (ETC)). Clinical and CEUS data were used to develop and validate the clinical model, DL model, and DL combining with clinical (DL-C) model to predict ER after TA. The performance of these models was compared by the receiver operating characteristic curve (ROC) with the DeLong test. RESULTS: After a median follow-up of 56 months, 49% (99/207) of patients experienced ER. Three key clinical features (preoperative chemotherapy (PC), lymph node metastasis of the primary colorectal cancer (LMPCC), and T stage) were used to develop the clinical model. The DL model yielded better performance than the clinical model in the ETC (AUC: 0.67 for the clinical model, 0.76 for the DL model). The DL-C model significantly outperformed the clinical model and DL model (AUC: 0.78 for the DL-C model in the ETC; both, p < 0.001). CONCLUSIONS: The model based on CEUS can achieve satisfactory prediction and assist physicians during the therapeutic decision-making process in clinical practice. KEY POINTS: • This is an exploratory study in which ablation-related contrast-enhanced ultrasound (CEUS) data from consecutive patients with colorectal cancer liver metastasis (CRLM) were collected simultaneously at multiple institutions. • The deep learning combining with clinical (DL-C) model provided desirable performance for the prediction of early recurrence (ER) after thermal ablation (TA). • The DL-C model based on CEUS provides guidance for TA indication selection and making therapeutic decisions.

Development and validation of the diagnostic accuracy of artificial intelligence-assisted ultrasound in the classification of splenic trauma.

Background: Spleen is the most vulnerable organ in abdominal trauma. Ultrasound (US) has become an important examination method for splenic trauma. However, the sensitivity of routine US in the diagnosis of splenic trauma is low. Contrast-enhanced ultrasound (CEUS) can improve the sensitivity, but it also has some limitations. This study sought to explore the application value of artificial intelligence (AI)-assisted US in the classification of splenic trauma. Methods: The splenic injuries of Bama miniature pigs were established. A large number of ultrasonic images were collected. Then, 3-fold cross validation (CV) was used to establish the animal models. Next, clinical ultrasonic images were collected at multiple centers. All injuries were diagnosed by CEUS, enhanced CT or surgery. We used animal models to fine tune a small amount of human data, and then established the final AI splenic trauma recognition model. The whole model was constructed by averaging the prediction ability of the 3 fine-tuned models. Finally, 2 doctors' recognition US results of splenic trauma were compared to the AI recognition results. The area under the curve (AUC), sensitivity, specificity, negative predictive value, and positive predictive value were used to evaluate the diagnostic performance in diagnosis of spleen trauma. Results: (I) Based on the receiver operating characteristic (ROC) curves, the test cohort 1 (AUC =0.90) and 2 (AUC =0.84) had a similar performance. Based on the decision curve analysis (DCA) curves, while threshold smaller than 0.8, the proposed model had better performance on test cohort 1 than test cohort 2. Test cohort 1 had higher sensitivity (0.82 vs. 0.71, P<0.01) and higher specificity (0.88 vs. 0.81, P<0.01) than test cohort 2. (II) The diagnostic accuracy of the AI model was higher than that of doctor 1 (0.82 vs. 0.62, P<0.001) and doctor 2 (0.82 vs. 0.66, P<0.001), and its specificity was higher than that of doctor (0.88 vs. 0.78, P=0.001). Conclusions: AI-assisted US diagnosis of splenic trauma can significantly improve the ultrasonic diagnosis rate. We still need to increase the number of samples to further improve the diagnostic efficiency of the model.

Accurate and fast reconstruction for bioluminescence tomography based on adaptive Newton hard thresholding pursuit algorithm.

As a promising noninvasive medical imaging technique, bioluminescence tomography (BLT) dynamically offers three-dimensional visualization of tumor distribution in living animals. However, due to the high ill-posedness caused by the strong scattering property of biological tissues and the limited boundary measurements with noise, BLT reconstruction still cannot meet actual preliminary clinical application requirements. In our research, to recover 3D tumor distribution quickly and precisely, an adaptive Newton hard thresholding pursuit (ANHTP) algorithm is proposed to improve the performance of BLT. The ANHTP algorithm fully combines the advantages of sparsity constrained optimization and convex optimization to guarantee global convergence. More precisely, an adaptive sparsity adjustment strategy was developed to obtain the support set of the inverse system matrix. Based on the strong Wolfe line search criterion, a modified damped Newton algorithm was constructed to obtain optimal source distribution information. A series of numerical simulations and phantom and in vivo experiments show that ANHTP has high reconstruction accuracy, fast reconstruction speed, and good robustness. Our proposed algorithm can further increase the practicality of BLT in biomedical applications.

A deep learning model combining multimodal radiomics, clinical and imaging features for differentiating ocular adnexal lymphoma from idiopathic orbital inflammation.

OBJECTIVES: To evaluate the value of deep learning (DL) combining multimodal radiomics and clinical and imaging features for differentiating ocular adnexal lymphoma (OAL) from idiopathic orbital inflammation (IOI). METHODS: Eighty-nine patients with histopathologically confirmed OAL (n = 39) and IOI (n = 50) were divided into training and validation groups. Convolutional neural networks and multimodal fusion layers were used to extract multimodal radiomics features from the T1-weighted image (T1WI), T2-weighted image, and contrast-enhanced T1WI. These multimodal radiomics features were then combined with clinical and imaging features and used together to differentiate between OAL and IOI. The area under the curve (AUC) was used to evaluate DL models with different features under five-fold cross-validation. The Student t-test, chi-squared, or Fisher exact test was used for comparison of different groups. RESULTS: In the validation group, the diagnostic AUC of the DL model using combined features was 0.953 (95% CI, 0.895-1.000), higher than that of the DL model using multimodal radiomics features (0.843, 95% CI, 0.786-0.898, p < 0.01) or clinical and imaging features only (0.882, 95% CI, 0.782-0.982, p = 0.13). The DL model built on multimodal radiomics features outperformed those built on most bimodalities and unimodalities (p < 0.05). In addition, the DL-based analysis with the orbital cone area (covering both the orbital mass and surrounding tissues) was superior to that with the region of interest (ROI) covering only the mass area, although the difference was not significant (p = 0.33). CONCLUSIONS: DL-based analysis that combines multimodal radiomics features with clinical and imaging features may help to differentiate between OAL and IOI. KEY POINTS: • It is difficult to differentiate OAL from IOI due to the overlap in clinical and imaging manifestations. • Radiomics has shown potential for noninvasive diagnosis of different orbital lymphoproliferative disorders. • DL-based analysis combining radiomics and imaging and clinical features may help the differentiation between OAL and IOI.

VoxDMRN: a voxelwise deep max-pooling residual network for bioluminescence tomography reconstruction.

Bioluminescence tomography (BLT) has extensive applications in preclinical studies for cancer research and drug development. However, the spatial resolution of BLT is inadequate because the numerical methods are limited for solving the physical models of photon propagation and the restriction of using tetrahedral meshes for reconstruction. We conducted a series of theoretical derivations and divided the BLT reconstruction process into two steps: feature extraction and nonlinear mapping. Inspired by deep learning, a voxelwise deep max-pooling residual network (VoxDMRN) is proposed to establish the nonlinear relationship between the internal bioluminescent source and surface boundary density to improve the spatial resolution in BLT reconstruction. The numerical simulation and in vivo experiments both demonstrated that VoxDMRN greatly improves the reconstruction performance regarding location accuracy, shape recovery capability, dual-source resolution, robustness, and in vivo practicability.

Multispectral Differential Reconstruction Strategy for Bioluminescence Tomography.

Bioluminescence tomography (BLT) is a promising in vivo molecular imaging tool that allows non-invasive monitoring of physiological and pathological processes at the cellular and molecular levels. However, the accuracy of the BLT reconstruction is significantly affected by the forward modeling errors in the simplified photon propagation model, the measurement noise in data acquisition, and the inherent ill-posedness of the inverse problem. In this paper, we present a new multispectral differential strategy (MDS) on the basis of analyzing the errors generated from the simplification from radiative transfer equation (RTE) to diffusion approximation and data acquisition of the imaging system. Through rigorous theoretical analysis, we learn that spectral differential not only can eliminate the errors caused by the approximation of RTE and imaging system measurement noise but also can further increase the constraint condition and decrease the condition number of system matrix for reconstruction compared with traditional multispectral (TM) reconstruction strategy. In forward simulations, energy differences and cosine similarity of the measured surface light energy calculated by Monte Carlo (MC) and diffusion equation (DE) showed that MDS can reduce the systematic errors in the process of light transmission. In addition, in inverse simulations and in vivo experiments, the results demonstrated that MDS was able to alleviate the ill-posedness of the inverse problem of BLT. Thus, the MDS method had superior location accuracy, morphology recovery capability, and image contrast capability in the source reconstruction as compared with the TM method and spectral derivative (SD) method. In vivo experiments verified the practicability and effectiveness of the proposed method.

Total Variation Constrained Graph Manifold Learning Strategy for Cerenkov Luminescence Tomography.

Harnessing the power and flexibility of radiolabeled molecules, Cerenkov luminescence tomography (CLT) provides a novel technique for non-invasive visualisation and quantification of viable tumour cells in a living organism. However, owing to the photon scattering effect and the ill-posed inverse problem, CLT still suffers from insufficient spatial resolution and shape recovery in various preclinical applications. In this study, we proposed a total variation constrained graph manifold learning (TV-GML) strategy for achieving accurate spatial location, dual-source resolution, and tumour morphology. TV-GML integrates the isotropic total variation term and dynamic graph Laplacian constraint to make a trade-off between edge preservation and piecewise smooth region reconstruction. Meanwhile, the tetrahedral mesh-Cartesian grid pair method based on the k-nearest neighbour, and the adaptive and composite Barzilai-Borwein method, were proposed to ensure global super linear convergence of the solution of TV-GML. The comparison results of both simulation experiments and in vivo experiments further indicated that TV-GML achieved superior reconstruction performance in terms of location accuracy, dual-source resolution, shape recovery capability, robustness, and in vivo practicability. Significance: We believe that this novel method will be beneficial to the application of CLT for quantitative analysis and morphological observation of various preclinical applications and facilitate the development of the theory of solving inverse problem.