Adaptive Grouping Block Sparse Bayesian Learning Method for Accurate and Robust Reconstruction in Bioluminescence Tomography

Logarithmic total variation regularization via preconditioned conjugate gradient method for sparse reconstruction of bioluminescence tomography

BACKGROUND AND OBJECTIVE

Bioluminescence Tomography (BLT) is a powerful optical molecular imaging technique that enables the noninvasive investigation of dynamic biological phenomena. It aims to reconstruct the three-dimensional spatial distribution of bioluminescent sources from optical measurements collected on the surface of the imaged object. However, BLT reconstruction is a challenging ill-posed problem due to the scattering effect of light and the limitations in detecting surface photons, which makes it difficult for existing methods to achieve satisfactory reconstruction results. In this study, we propose a novel method for sparse reconstruction of BLT based on a preconditioned conjugate gradient with logarithmic total variation regularization (PCG-logTV).

METHOD

This PCG-logTV method incorporates the sparsity of overlapping groups and enhances the sparse structure of these groups using logarithmic functions, which can preserve edge features and achieve more stable reconstruction results in BLT. To accelerate the convergence of the algorithm solution, we use the preconditioned conjugate gradient iteration method on the objective function and obtain the reconstruction results. We demonstrate the performance of our proposed method through numerical simulations and in vivo experiment.

RESULTS AND CONCLUSIONS

The results show that the PCG-logTV method obtains the most accurate reconstruction results, and the minimum position error (LE) is 0.254mm, which is 26%, 31% and 34% of the FISTA (0.961), IVTCG (0.81) and L1-TV (0.739) methods, and the root mean square error (RMSE) and relative intensity error (RIE) are the smallest, indicating that it is closest to the real light source. In addition, compared with the other three methods, the PCG-logTV method also has the highest DICE similarity coefficient, which is 0.928, which means that this method can effectively reconstruct the three-dimensional spatial distribution of bioluminescent light sources, has higher resolution and robustness, and is beneficial to the preclinical and clinical studies of BLT.

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.

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.

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.

ABPO-TVSCAD: alternating Bregman proximity operators approach based on TVSCAD regularization for bioluminescence tomography

Objective.Bioluminescence tomography (BLT) is a promising non-invasive optical medical imaging technique, which can visualize and quantitatively analyze the distribution of tumor cells in living tissues. However, due to the influence of photon scattering effect and ill-conditioned inverse problem, the reconstruction result is unsatisfactory. The purpose of this study is to improve the reconstruction performance of BLT.Approach.An alternating Bregman proximity operators (ABPO) method based on TVSCAD regularization is proposed for BLT reconstruction. TVSCAD combines the anisotropic total variation (TV) regularization constraints and the non-convex smoothly clipped absolute deviation (SCAD) penalty constraints, to make a trade-off between the sparsity and edge preservation of the source. ABPO approach is used to solve the TVSCAD model (ABPO-TVSCAD for short). In addition, to accelerate the convergence speed of the ABPO, we adapt the strategy of shrinking the permission source region, which further improves the performance of ABPO-TVSCAD.Main results.The results of numerical simulations andin vivoxenograft mouse experiment show that our proposed method achieved superior accuracy in spatial localization and morphological reconstruction of bioluminescent source.Significance.ABPO-TVSCAD is an effective and robust reconstruction method for BLT, and we hope that this method can promote the development of optical molecular tomography.

Elastic net-based non-negative iterative three-operator splitting strategy for Cerenkov luminescence tomography

Cerenkov luminescence tomography (CLT) provides a powerful optical molecular imaging technique for non-invasive detection and visualization of radiopharmaceuticals in living objects. However, the severe photon scattering effect causes ill-posedness of the inverse problem, and the location accuracy and shape recovery of CLT reconstruction results are unsatisfactory for clinical application. Here, to improve the reconstruction spatial location accuracy and shape recovery ability, a non-negative iterative three operator splitting (NNITOS) strategy based on elastic net (EN) regularization was proposed. NNITOS formalizes the CLT reconstruction as a non-convex optimization problem and splits it into three operators, the least square, L_1/2-norm regularization, and adaptive grouping manifold learning, then iteratively solved them. After stepwise iterations, the result of NNITOS converged progressively. Meanwhile, to speed up the convergence and ensure the sparsity of the solution, shrinking the region of interest was utilized in this strategy. To verify the effectiveness of the method, numerical simulations and in vivo experiments were performed. The result of these experiments demonstrated that, compared to several methods, NNITOS can achieve superior performance in terms of location accuracy, shape recovery capability, and robustness. We hope this work can accelerate the clinical application of CLT in the future.

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.

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.