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

Coarse-graining of perplexity for the spatial distribution of molecules

Coarse-graining of the perplexity for the spatial distribution of molecules

Spatial distribution of the Shannon entropy for mass spectrometry imaging

Simultaneous Estimation and Segmentation from Projection Data in Dynamic PET.

State-Space Reconstruction of Pet Parametric Maps

Spatial components of molecular tissue biology

Quantifying and Visualizing Hidden Preferential Aggregations Amid Heterogeneity

Relative Resolution: An Analysis with the Kullback-Leibler Entropy

Characterization of Anomalous Diffusion in Porous Biological Tissues Using Fractional Order Derivatives and Entropy

Using Ripley's K-function to Characterize Clustering In 3-Dimensional Point Patterns With a Case Study in Atom Probe Tomography

FunSpace: A functional and spatial analytic approach to cell imaging data using entropy measures

Quantifying Subtle Regions of Order and Disorder in Tumor Architecture by Calculating the Nearest-Neighbor Angular Profile

Navigating Perplexity: A linear relationship with the data set size in t-SNE embeddings

Feature Analysis and Visualization of 3D Scalar Field with the Applications to the Macromolecule

Information Entropy-Based Strategy for the Quantitative Evaluation of Extensive Hyperspectral Images to Better Unveil Spatial Heterogeneity in Mass Spectrometry Imaging

Principal nested shape space analysis of molecular dynamics data

Characterizing Complex Spatiotemporal Patterns from Entropy Measures

From Local Atomic Environments to Molecular Information Entropy

Harmonic representations of regions and interactions in spatial transcriptomics

Graph Fourier transform for spatial omics representation and analyses of complex organs