A Novel Adaptive Parameter Search Elastic Net Method for Fluorescent Molecular Tomography

Performance enhancement of diffuse fluorescence tomography based on an extended Kalman filtering-long short term memory neural network correction model

To alleviate the ill-posedness of diffuse fluorescence tomography (DFT) reconstruction and improve imaging quality and speed, a model-derived deep-learning method is proposed by combining extended Kalman filtering (EKF) with a long short term memory (LSTM) neural network, where the iterative process parameters acquired by implementing semi-iteration EKF (SEKF) served as inputs to the LSTM neural network correction model for predicting the optimal fluorescence distributions. To verify the effectiveness of the SEKF-LSTM algorithm, a series of numerical simulations, phantom and in vivo experiments are conducted, and the experimental results are quantitatively evaluated and compared with the traditional EKF algorithm. The simulation experimental results show that the proposed new algorithm can effectively improve the reconstructed image quality and reconstruction speed. Importantly, the LSTM correction model trained by the simulation data also obtains satisfactory results in the experimental data, suggesting that the SEKF-LSTM algorithm possesses strong generalization ability and great potential for practical applications.

Preliminary assessment of three quantitative approaches for estimating time-since-deposition from autofluorescence and morphological profiles of cell populations from forensic biological samples

Determining when DNA recovered from a crime scene transferred from its biological source, i.e., a sample's 'time-since-deposition' (TSD), can provide critical context for biological evidence. Yet, there remains no analytical techniques for TSD that are validated for forensic casework. In this study, we investigate whether morphological and autofluorescence measurements of forensically-relevant cell populations generated with Imaging Flow Cytometry (IFC) can be used to predict the TSD of 'touch' or trace biological samples. To this end, three different prediction frameworks for estimating the number of day(s) for TSD were evaluated: the elastic net, gradient boosting machines (GBM), and generalized linear mixed model (GLMM) LASSO. Additionally, we transformed these continuous predictions into a series of binary classifiers to evaluate the potential utility for forensic casework. Results showed that GBM and GLMM-LASSO showed the highest accuracy, with mean absolute error estimates in a hold-out test set of 29 and 21 days, respectively. Binary classifiers for these models correctly binned 94-96% and 98-99% of the age estimates as over/under 7 or 180 days, respectively. This suggests that predicted TSD using IFC measurements coupled to one or, possibly, a combination binary classification decision rules, may provide probative information for trace biological samples encountered during forensic casework.

Multi-target reconstruction based on subspace decision optimization for bioluminescence tomography

BACKGROUND AND OBJECTIVE

Bioluminescence tomography (BLT) is a noninvasive optical imaging technique that provides qualitative and quantitative information on the spatial distribution of tumors in living animals. Researchers have proposed a list of algorithms and strategies for BLT reconstruction to improve its reconstruction quality. However, multi-target BLT reconstruction remains challenging in practical clinical applications due to the mutual interference of optical signals and difficulty in source separation.

METHODS

To solve this problem, this study proposes the subspace decision optimization (SDO) approach based on the traditional iterative permissible region strategy. The SDO approach transforms a single permissible region into multiple subspaces by clustering analysis. These subspaces are shrunk based on subspace shrinking optimization to achieve spatial continuity of the permissible regions. In addition, these subspaces are merged to construct a new permissible region and then the next iteration of reconstruction is carried out to ensure the stability of the results. Finally, all the iterative results are optimized based on the normal distribution model and the distribution properties of the targets to ensure the sparsity of each target and the non-biasing of the overall results.

RESULTS

Experimental results show that the SDO approach can automatically identify and separate different targets, ensuring the accuracy and quality of multi-target BLT reconstruction results. Meanwhile, SDO can combine various types of reconstruction algorithms and provide stable and high-quality reconstruction results independent of the algorithm parameters.

CONCLUSIONS

The SDO approach provides an integrated solution to the multi-target BLT reconstruction problem, realizing the whole process including target recognition, separation, reconstruction, and result enhancement, which can extend the application domain of BLT.

Generalized conditional gradient method with adaptive regularization parameters for fluorescence molecular tomography

Fluorescence molecular tomography (FMT) is an optical imaging technology with the ability of visualizing the three-dimensional distribution of fluorescently labelled probes in vivo. However, due to the light scattering effect and ill-posed inverse problems, obtaining satisfactory FMT reconstruction is still a challenging problem. In this work, to improve the performance of FMT reconstruction, we proposed a generalized conditional gradient method with adaptive regularization parameters (GCGM-ARP). In order to make a tradeoff between the sparsity and shape preservation of the reconstruction source, and to maintain its robustness, elastic-net (EN) regularization is introduced. EN regularization combines the advantages of L₁-norm and L₂-norm, and overcomes the shortcomings of traditional L_p-norm regularization, such as over-sparsity, over-smoothness, and non-robustness. Thus, the equivalent optimization formulation of the original problem can be obtained. To further improve the performance of the reconstruction, the L-curve is adopted to adaptively adjust the regularization parameters. Then, the generalized conditional gradient method (GCGM) is used to split the minimization problem based on EN regularization into two simpler sub-problems, which are determining the direction of the gradient and the step size. These sub-problems are addressed efficiently to obtain more sparse solutions. To assess the performance of our proposed method, a series of numerical simulation experiments and in vivo experiments were implemented. The experimental results show that, compared with other mathematical reconstruction methods, GCGM-ARP method has the minimum location error (LE) and relative intensity error (RIE), and the maximum dice coefficient (Dice) in the case of different sources number or shape, or Gaussian noise of 5%-25%. This indicates that GCGM-ARP has superior reconstruction performance in source localization, dual-source resolution, morphology recovery, and robustness. In conclusion, the proposed GCGM-ARP is an effective and robust strategy for FMT reconstruction in biomedical application.

Reconstruction based on adaptive group least angle regression for fluorescence molecular tomography

Fluorescence molecular tomography can combine two-dimensional fluorescence imaging with anatomical information to reconstruct three-dimensional images of tumors. Reconstruction based on traditional regularization with tumor sparsity priors does not take into account that tumor cells form clusters, so it performs poorly when multiple light sources are used. Here we describe reconstruction based on an "adaptive group least angle regression elastic net" (AGLEN) method, in which local spatial structure correlation and group sparsity are integrated with elastic net regularization, followed by least angle regression. The AGLEN method works iteratively using the residual vector and a median smoothing strategy in order to adaptively obtain a robust local optimum. The method was verified using numerical simulations as well as imaging of mice bearing liver or melanoma tumors. AGLEN reconstruction performed better than state-of-the-art methods with different sizes of light sources at different distances from the sample and in the presence of Gaussian noise at 5-25%. In addition, AGLEN-based reconstruction accurately imaged tumor expression of cell death ligand-1, which can guide immunotherapy.

An adaptive parameter selection strategy based on maximizing the probability of data for robust fluorescence molecular tomography reconstruction

To alleviate the ill-posed of the inverse problem in fluorescent molecular tomography (FMT), many regularization methods based on L₂ or L₁ norm have been proposed. Whereas, the quality of regularization parameters affects the performance of the reconstruction algorithm. Some classical parameter selection strategies usually need initialization of parameter range and high computing costs, which is not universal in the practical application of FMT. In this paper, an universally applicable adaptive parameter selection method based on maximizing the probability of data (MPD) strategy was proposed. This strategy used maximum a posteriori (MAP) estimation and maximum likelihood (ML) estimation to establish a regularization parameters model. The stable optimal regularization parameters can be determined by multiple iterative estimates. Numerical simulations and in vivo experiments show that MPD strategy can obtain stable regularization parameters for both regularization algorithms based on L₂ or L₁ norm and achieve good reconstruction performance.

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 L₁ and L₂ 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.

Multi-target reconstruction strategy based on blind source separation of surface measurement signals in FMT

Fluorescence molecular tomography (FMT) is a promising molecular imaging technique for tumor detection in the early stage. High-precision multi-target reconstructions are necessary for quantitative analysis in practical FMT applications. The existing reconstruction methods perform well in retrieving a single fluorescent target but may fail in reconstructing a multi-target, which remains an obstacle to the wider application of FMT. In this paper, a novel multi-target reconstruction strategy based on blind source separation (BSS) of surface measurement signals was proposed, which transformed the multi-target reconstruction problem into multiple single-target reconstruction problems. Firstly, by multiple points excitation, multiple groups of superimposed measurement signals conforming to the conditions of BSS were constructed. Secondly, an efficient nonnegative least-correlated component analysis with iterative volume maximization (nLCA-IVM) algorithm was applied to construct the separation matrix, and the superimposed measurement signals were separated into the measurements of each target. Thirdly, the least squares fitting method was combined with BSS to determine the number of fluorophores indirectly. Lastly, each target was reconstructed based on the extracted surface measurement signals. Numerical simulations and in vivo experiments proved that it has the ability of multi-target resolution for FMT. The encouraging results demonstrate the significant effectiveness and potential of our method for practical FMT applications.

Compressed Sensing Photoacoustic Imaging Reconstruction Using Elastic Net Approach

Photoacoustic imaging involves reconstructing an estimation of the absorbed energy density distribution from measured ultrasound data. The reconstruction task based on incomplete and noisy experimental data is usually an ill-posed problem that requires regularization to obtain meaningful solutions. The purpose of the work is to propose an elastic network (EN) model to improve the quality of reconstructed photoacoustic images. To evaluate the performance of the proposed method, a series of numerical simulations and tissue-mimicking phantom experiments are performed. The experiment results indicate that, compared with the L ₁-norm and L ₂-normbased regularization methods with different numerical phantoms, Gaussian noise of 10-50 dB, and different regularization parameters, the EN method with α = 0.5 has better image quality, calculation speed, and antinoise ability.

A Fast and Automated FMT/XCT Reconstruction Strategy Based on Standardized Imaging Space

The traditional finite element method-based fluorescence molecular tomography (FMT)/ X-ray computed tomography (XCT) imaging reconstruction suffers from complicated mesh generation and dual-modality image data fusion, which limits the application of in vivo imaging. To solve this problem, a novel standardized imaging space reconstruction (SISR) method for the quantitative determination of fluorescent probe distributions inside small animals was developed. In conjunction with a standardized dual-modality image data fusion technology, and novel reconstruction strategy based on Laplace regularization and L1-fused Lasso method, the in vivo distribution can be calculated rapidly and accurately, which enables standardized and algorithm-driven data process. We demonstrated the method's feasibility through numerical simulations and quantitatively monitored in vivo programmed death ligand 1 (PD-L1) expression in mouse tumor xenografts, and the results demonstrate that our proposed SISR can increase data throughput and reproducibility, which helps to realize the dynamically and accurately in vivo imaging.

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.

L1-L2 norm regularization via forward-backward splitting for fluorescence molecular tomography

Fluorescent molecular tomography (FMT) is a highly sensitive and noninvasive imaging approach for providing three-dimensional distribution of fluorescent marker probes. However, owing to its light scattering effect and the ill-posedness of inverse problems, it is challenging to develop an efficient reconstruction algorithm that can achieve the exact location and morphology of the fluorescence source. In this study, therefore, in order to satisfy the need for early tumor detection and improve the sparsity of solution, we proposed a novel L ₁-L ₂ norm regularization via the forward-backward splitting method for enhancing the FMT reconstruction accuracy and the robustness. By fully considering the highly coherent nature of the system matrix of FMT, it operates by splitting the objective to be minimized into simpler functions, which are dealt with individually to obtain a sparser solution. An analytic solution of L ₁-L ₂ norm proximal operators and a forward-backward splitting algorithm were employed to efficiently solve the nonconvex L ₁-L ₂ norm minimization problem. Numerical simulations and an in-vivo glioma mouse model experiment were conducted to evaluate the performance of our algorithm. The comparative results of these experiments demonstrated that the proposed algorithm obtained superior reconstruction performance in terms of spatial location, dual-source resolution, and in-vivo practicability. It was believed that this study would promote the preclinical and clinical applications of FMT in early tumor detection.