Abstract:Purpose Dynamic positron emission tomography (PET) is known for its ability to extract spatiotemporal information of a radio tracer in living tissue. Information of different functional regions based on an accurate reconstruction of the activity images and kinetic parametric images has been widely studied and can be useful in research and clinical setting for diagnosis and other quantitative tasks. In this paper, our purpose is to present a novel framework for estimating the kinetic parametric images directly from the raw measurement data together with a simultaneous segmentation accomplished through kinetic parameters clustering. Method An iterative framework is proposed to estimate the kinetic parameter image, activity map and do the segmentation simultaneously from the complete dynamic PET projection data. The clustering process is applied to the kinetic parameter variable rather than to the traditional activity distribution so as to achieve accurate discrimination between different functional areas. Prior information such as total variation regularization is incorporated to reduce the noise in the PET images and a sparseness constraint is integrated to guarantee the solution for kinetic parameters due to the over complete dictionary. Alternating direction method of multipliers (ADMM) method is used to solve the optimization problem. The proposed algorithm was validated with experiments on Monte Carlo-simulated phantoms and real patient data. Symbol error rate (SER) was defined to evaluate the performance of clustering. Bias and variance of the reconstruction activity images were calculated based on ground truth. Relative mean square error (MSE) was used to evaluate parametric results quantitatively. Result In brain phantom experiment, when counting rate is 1 x 10(6), the bias (variance) of our method is 0.1270 (0.0281), which is lower than maximum likelihood expectation maximization (MLEM) 0.1637 (0.0410) and direct estimation without segmentation (DE) 0.1511 (0.0326). In the Zubal phantom experiment, our method has the lowest bias (variance) 0.1559 (0.0354) with 1 x 10(5) counting rate, compared with DE 0.1820 (0.0435) and MLEM 0.3043 (0.0644). As for classification, the SER of our method is 18.87% which is the lowest among MLEM + k-means, DE + k-means, and kinetic spectral clustering (KSC). Brain data with MR reference and real patient results also show that the proposed method can get images with clear structure by visual inspection. Conclusion In this paper, we presented a joint reconstruction framework for simultaneously estimating the activity distribution, parametric images, and parameter-based segmentation of the ROIs into different functional areas. Total variation regularization is performed on the activity distribution domain to suppress noise and preserve the edges between ROIs. An over complete dictionary for time activity curve basis is constructed. SER, bias, variance, and MSE were calculated to show the effectiveness of the proposed method. (c) 2018 American Association of Physicists in Medicine

Accelerated model‐based T1, T2* and proton density mapping using a Bayesian approach with automatic hyperparameter estimation

Model-based T1, T2* and Proton Density Mapping Using a Bayesian Approach with Parameter Estimation and Complementary Undersampling Patterns

Direct and Accelerated Parameter Mapping Using the Unscented Kalman Filter

Rapid 3D Dynamic Arterial Spin Labeling with a Sparse Model-Based Image Reconstruction

HF-SENSE: an Improved Partially Parallel Imaging Using a High-Pass Filter

Three-dimensional Simultaneous Brain Mapping of T1, T2, T2*and Magnetic Susceptibility with MR Multitasking

Three-dimensional simultaneous brain T 1 , T 2 , and ADC mapping with MR Multitasking.

Exploiting Parameter Sparsity in Model-Based Reconstruction to Accelerate Proton Density and T2 Mapping

Joint MAPLE: Accelerated joint T1 and T2* mapping with scan‐specific self‐supervised networks

Highly accelerated parameter mapping using model-based alternating reconstruction coupling fitting

Accelerated exponential parameterization of T2 relaxation with model-driven low rank and sparsity priors (MORASA).

Accelerating MR Parameter Mapping Using Nonlinear Compressive Manifold Learning and Regularized Pre-Imaging

Fast Whole-Brain MR Multi-Parametric Mapping with Scan-Specific Self-Supervised Networks

Three-dimensional simultaneous brain T<sub>1</sub>, T<sub>2</sub>, and ADC mapping with MR Multitasking

Simultaneous Estimation and Segmentation from Projection Data in Dynamic PET.

SCOPE: signal compensation for low-rank plus sparse matrix decomposition for fast parameter mapping.

CAMP-Net: Consistency-Aware Multi-Prior Network for Accelerated MRI Reconstruction

MANTIS: Model-Augmented Neural neTwork with Incoherent k-space Sampling for efficient MR T2 mapping

Highly Accelerated MR Parametric Mapping by Undersampling the K-Space and Reducing the Contrast Number Simultaneously with Deep Learning

Accelerated Brain Tumor Dynamic Contrast‐enhanced MRI Using Adaptive Pharmaco‐Kinetic Model Constrained Method

Bayesian Pharmacokinetic Modeling of Dynamic Contrast-Enhanced Magnetic Resonance Imaging: Validation and Application