A review of advances in imaging methodology in fluorescence molecular tomography

Projected algebraic reconstruction technique-network for high-fidelity diffuse fluorescence tomography reconstruction

We propose a model-driven projected algebraic reconstruction technique (PART)-network (PART-Net) that leverages the advantages of the traditional model-based method and the neural network to improve the imaging quality of diffuse fluorescence tomography. In this algorithm, nonnegative prior information is incorporated into the ART iteration process to better guide the optimization process, and thereby improve imaging quality. On this basis, PART in conjunction with a residual convolutional neural network is further proposed to obtain high-fidelity image reconstruction. The numerical simulation results demonstrate that the PART-Net algorithm effectively improves noise robustness and reconstruction accuracy by at least 1-2 times and exhibits superiority in spatial resolution and quantification, especially for a small-sized target (r=2m m), compared with the traditional ART algorithm. Furthermore, the phantom and in vivo experiments verify the effectiveness of the PART-Net, suggesting strong generalization capability and a great potential for practical applications.

Fluorescence molecular tomography based on an online maximum a posteriori estimation algorithm

Fluorescence molecular tomography (FMT) is a non-invasive, radiation-free, and highly sensitive optical molecular imaging technique for early tumor detection. However, inadequate measurement information along with significant scattering of near-infrared light within the tissue leads to high ill-posedness in the inverse problem of FMT. To improve the quality and efficiency of FMT reconstruction, we build a reconstruction model based on log-sum regularization and introduce an online maximum a posteriori estimation (OPE) algorithm to solve the non-convex optimization problem. The OPE algorithm approximates a stationary point by evaluating the gradient of the objective function at each iteration, and its notable strength lies in the remarkable speed of convergence. The results of simulations and experiments demonstrate that the OPE algorithm ensures good reconstruction quality and exhibits outstanding performance in terms of reconstruction efficiency.

Two-step reconstruction framework of fluorescence molecular tomography based on energy statistical probability

Fluorescence molecular tomography (FMT), as a promising technique for early tumor detection, can non-invasively visualize the distribution of fluorescent marker probe three-dimensionally. However, FMT reconstruction is a severely ill-posed problem, which remains an obstacle to wider application of FMT. In this paper, a two-step reconstruction framework was proposed for FMT based on the energy statistical probability. First, the tissue structural information obtained from computed tomography (CT) is employed to associate the tissue optical parameters for rough solution in the global region. Then, according to the global-region reconstruction results, the probability that the target belongs to each region can be calculated. The region with the highest probability is delineated as region of interest to realize accurate and fast source reconstruction. Numerical simulations and in vivo experiments were carried out to evaluate the effectiveness of the proposed framework. The encouraging results demonstrate the significant effectiveness and potential of our method for practical FMT applications.

Simultaneous reconstruction of 3D fluorescence distribution and object surface using structured light illumination and dual-camera detection

Fluorescence molecular tomography (FMT) serves as a noninvasive modality for visualizing volumetric fluorescence distribution within biological tissues, thereby proving to be an invaluable imaging tool for preclinical animal studies. The conventional FMT relies upon a point-by-point raster scan strategy, enhancing the dataset for subsequent reconstruction but concurrently elongating the data acquisition process. The resultant diminished temporal resolution has persistently posed a bottleneck, constraining its utility in dynamic imaging studies. We introduce a novel system capable of simultaneous FMT and surface extraction, which is attributed to the implementation of a rapid line scanning approach and dual-camera detection. The system performance was characterized through phantom experiments, while the influence of scanning line density on reconstruction outcomes has been systematically investigated via both simulation and experiments. In a proof-of-concept study, our approach successfully captures a moving fluorescence bolus in three dimensions with an elevated frame rate of approximately 2.5 seconds per frame, employing an optimized scan interval of 5 mm. The notable enhancement in the spatio-temporal resolution of FMT holds the potential to broaden its applications in dynamic imaging tasks, such as surgical navigation.

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.

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.

Application of the Powell algorithm for estimating optical properties of a semitransparent medium based on time-domain information

Accurate estimation of the optical properties of a semitransparent medium is crucial in various engineering applications. This study introduces the Powell algorithm to estimate the optical properties of a 2D semitransparent slab. The time-domain radiative transfer equation is solved using the discrete ordinate method. The radiative intensity on the medium's surface serves as the measurement signal for the inverse analysis. The results demonstrate that the Powell algorithm accurately estimates the absorption coefficient, scattering coefficient, and scattering asymmetry factor. For simultaneous reconstruction of these three parameters, it is recommended to use eight signal detectors on both the left and right sides of the medium. Even when the standard measurement error is increased to 15%, the relative errors for these three parameters remain low, at 1.87%, 1.379%, and 0.194%.

In vivo three-dimensional brain imaging with chemiluminescence probes in Alzheimer's disease models

Optical three-dimensional (3D) molecular imaging is highly desirable for providing precise distribution of the target-of-interest in disease models. However, such 3D imaging is still far from wide applications in biomedical research; 3D brain optical molecular imaging, in particular, has rarely been reported. In this report, we designed chemiluminescence probes with high quantum yields, relatively long emission wavelengths, and high signal-to-noise ratios to fulfill the requirements for 3D brain imaging in vivo. With assistance from density-function theory (DFT) computation, we designed ADLumin-Xs by locking up the rotation of the double bond via fusing the furan ring to the phenyl ring. Our results showed that ADLumin-5 had a high quantum yield of chemiluminescence and could bind to amyloid beta (Aβ). Remarkably, ADLumin-5's radiance intensity in brain areas could reach 4 × 107 photon/s/cm2/sr, which is probably 100-fold higher than most chemiluminescence probes for in vivo imaging. Because of its strong emission, we demonstrated that ADLumin-5 could be used for in vivo 3D brain imaging in transgenic mouse models of Alzheimer's disease.

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.

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.

In Vivo Three-dimensional Brain Imaging with Chemiluminescence Probes in Alzheimer's Disease Models

Optical three-dimensional (3D) molecular imaging is highly desirable for providing precise distribution of the target-of-interest in disease models. However, such 3D imaging is still far from wide applications in biomedical research; 3D brain optical molecular imaging, in particular, has rarely been reported. In this report, we designed chemiluminescence probes with high quantum yields (QY), relatively long emission wavelengths, and high signal-to-noise ratios (SNRs) to fulfill the requirements for 3D brain imaging in vivo. With assistance from density-function theory (DFT) computation, we designed ADLumin-Xs by locking up the rotation of the double-bond via fusing the furan ring to the phenyl ring. Our results showed that ADLumin-5 had a high quantum yield of chemiluminescence and could bind to amyloid beta (Aβ). Remarkably, ADLumin-5's radiance intensity in brain areas could reach 4×107 photon/s/cm2/sr, which is probably 100-fold higher than most chemiluminescence probes for in vivo imaging. Because of its strong emission, we demonstrated that ADLumin-5 could be used for in vivo 3D brain imaging in transgenic mouse models of Alzheimer's disease (AD).

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.

Spiral volumetric optoacoustic tomography for imaging whole-body biodynamics in small animals

Fast tracking of biological dynamics across multiple murine organs using the currently commercially available whole-body preclinical imaging systems is hindered by their limited contrast, sensitivity and spatial or temporal resolution. Spiral volumetric optoacoustic tomography (SVOT) provides optical contrast, with an unprecedented level of spatial and temporal resolution, by rapidly scanning a mouse using spherical arrays, thus overcoming the current limitations in whole-body imaging. The method enables the visualization of deep-seated structures in living mammalian tissues in the near-infrared spectral window, while further providing unrivalled image quality and rich spectroscopic optical contrast. Here, we describe the detailed procedures for SVOT imaging of mice and provide specific details on how to implement a SVOT system, including component selection, system arrangement and alignment, as well as the image processing methods. The step-by-step guide for the rapid panoramic (360°) head-to-tail whole-body imaging of a mouse includes the rapid visualization of contrast agent perfusion and biodistribution. The isotropic spatial resolution possible with SVOT can reach 90 µm in 3D, while alternative steps enable whole-body scans in less than 2 s, unattainable with other preclinical imaging modalities. The method further allows the real-time (100 frames per second) imaging of biodynamics at the whole-organ level. The multiscale imaging capacity provided by SVOT can be used for visualizing rapid biodynamics, monitoring responses to treatments and stimuli, tracking perfusion, and quantifying total body accumulation and clearance dynamics of molecular agents and drugs. Depending on the imaging procedure, the protocol requires 1-2 h to complete by users trained in animal handling and biomedical imaging.

Two-stage deep learning method for sparse-view fluorescence molecular tomography reconstruction

Fluorescence molecular tomography (FMT) is a preclinical optical tomographic imaging technique that can trace various physiological and pathological processes at the cellular or even molecular level. Reducing the number of FMT projection views can improve the data acquisition speed, which is significant in applications such as dynamic problems. However, a reduction in the number of projection views will dramatically aggravate the ill-posedness of the FMT inverse problem and lead to significant degradation of the reconstructed images. To deal with this problem, we have proposed a deep-learning-based reconstruction method for sparse-view FMT that only uses four perpendicular projection views and divides the image reconstruction into two stages: image restoration and inverse Radon transform. In the first stage, the projection views of the surface fluorescence are restored to eliminate the blur derived from photon diffusion through a fully convolutional neural network. In the second stage, another convolutional neural network is used to implement the inverse Radon transform between the restored projections from the first stage and the reconstructed transverse slices. Numerical simulation and phantom and mouse experiments are carried out. The results show that the proposed method can effectively deal with the image reconstruction problem of sparse-view FMT.

Adaptive Sparsity Orthogonal Least Square with Neighbor Strategy for Fluorescence Molecular Tomography

Fluorescence molecular tomography (FMT) is a highly sensitive and noninvasive optical imaging technique which has been widely applied to disease diagnosis and drug discovery. However, FMT reconstruction is a highly ill-posed problem. In this work, L0-norm regularization is employed to construct the mathematical model of the inverse problem of FMT. And an adaptive sparsity orthogonal least square with a neighbor strategy (ASOLS-NS) is proposed to solve this model. This algorithm can provide an adaptive sparsity and can establish the candidate sets by a novel neighbor expansion strategy for the orthogonal least square (OLS) algorithm. Numerical simulation experiments have shown that the ASOLS-NS improves the reconstruction of images, especially for the double targets reconstruction.Clinical relevance- The purpose of this work is to improve the reconstruction results of FMT. Current experiments are focused on simulation experiments, and the proposed algorithm will be applied to the clinical tumor detection in the future.

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.

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.

Robust reconstruction of fluorescence molecular tomography based on adaptive adversarial learning strategy

Objective.Fluorescence molecular tomography (FMT) is a promising molecular imaging modality for quantifying the three-dimensional (3D) distribution of tumor probes in small animals. However, traditional deep learning reconstruction methods that aim to minimize the mean squared error (MSE) and iterative regularization algorithms that rely on optimal parameters are typically influenced by strong noise, resulting in poor FMT reconstruction robustness.Approach.In this letter, we propose an adaptive adversarial learning strategy (3D-UR-WGAN) to achieve robust FMT reconstructions. Unlike the pixel-based MSE criterion in traditional CNNs or the regularization strategy in iterative solving schemes, the reconstruction strategy can greatly facilitate the performance of the network models through alternating loop training of the generator and the discriminator. Second, the reconstruction strategy combines the adversarial loss in GANs with the L1 loss to significantly enhance the robustness and preserve image details and textual information.Main results.Both numerical simulations and physical phantom experiments demonstrate that the 3D-UR-WGAN method can adaptively eliminate the effects of different noise levels on the reconstruction results, resulting in robust reconstructed images with reduced artifacts and enhanced image contrast. Compared with the state-of-the-art methods, the proposed method achieves better reconstruction performance in terms of target shape recovery and localization accuracy.Significance.This adaptive adversarial learning reconstruction strategy can provide a possible paradigm for robust reconstruction in complex environments, and also has great potential to provide an alternative solution for solving the problem of poor robustness encountered in other optical imaging modalities such as diffuse optical tomography, bioluminescence imaging, and Cherenkov luminescence imaging.

D2-RecST: Dual-domain joint reconstruction strategy for fluorescence molecular tomography based on image domain and perception domain

BACKGROUND AND OBJECTIVE

Fluorescence molecular tomography (FMT) is a promising molecular imaging modality for quantifying the three-dimensional (3D) distribution of fluorescent probes in small animals. Over the past few years, learning-based FMT reconstruction methods have achieved promising results. However, these methods typically attempt to minimize the mean-squared error (MSE) between the reconstructed image and the ground truth. Although signal-to-noise ratios (SNRs) are improved, they are susceptible to non-uniform artifacts and loss of structural detail, making it extremely challenging to obtain accurate and robust FMT reconstructions under noisy measurements.

METHODS

We propose a novel dual-domain joint strategy based on the image domain and perception domain for accurate and robust FMT reconstruction. First, we formulate an explicit adversarial learning strategy in the image domain, which greatly facilitates training and optimization through two enhanced networks to improve anti-noise ability. Besides, we introduce a novel transfer learning strategy in the perceptual domain to optimize edge details by providing perceptual priors for fluorescent targets. Collectively, the proposed dual-domain joint reconstruction strategy can significantly eliminate the non-uniform artifacts and effectively preserve the structural edge details.

RESULTS

Both numerical simulations and in vivo mouse experiments demonstrate that the proposed method markedly outperforms traditional and cutting-edge methods in terms of positioning accuracy, image contrast, robustness, and target morphological recovery.

CONCLUSIONS

The proposed method achieves the best reconstruction performance and has great potential to facilitate precise localization and 3D visualization of tumors in in vivo animal experiments.

In vivo accurate detection of the liver tumor with pharmacokinetic parametric images from dynamic fluorescence molecular tomography

SIGNIFICANCE

Pharmacokinetic parametric images in dynamic fluorescence molecular tomography (FMT) can describe three-dimensional (3D) physiological and pathological information inside biological tissues, potentially providing quantitative assessment tools for biological research and drug development.

AIM

In vivo imaging of the liver tumor with pharmacokinetic parametric images from dynamic FMT based on the differences in metabolic properties of indocyanine green (ICG) between normal liver cells and tumor liver cells inside biological tissues.

APPROACH

First, an orthotopic liver tumor mouse model was constructed. Then, with the help of the FMT/computer tomography (CT) dual-modality imaging system and the direct reconstruction algorithm, 3D imaging of liver metabolic parameters in nude mice was achieved to distinguish liver tumors from normal tissues. Finally, pharmacokinetic parametric imaging results were validated against in vitro anatomical results.

RESULTS

This letter demonstrates the ability of dynamic FMT to monitor the pharmacokinetic delivery of the fluorescent dye ICG in vivo, thus, enabling the distinction between normal and tumor tissues based on the pharmacokinetic parametric images derived from dynamic FMT.

CONCLUSIONS

Compared with CT structural imaging technology, dynamic FMT combined with compartmental modeling as an analytical method can obtain quantitative images of pharmacokinetic parameters, thus providing a more powerful research tool for organ function assessment, disease diagnosis and new drug development.