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

Constructing Reliable Parametric Images Using Enhanced GLLS for Dynamic SPECT

Enhanced Formation of Parametric Images Using Fast Regressive GLLS for Noisy Functional Imaging.

Fast and reliable estimation of multiple parametric images using an integrated method for dynamic SPECT.

[Analytical Algorithms Improving the Reconstruction of SPECT Image].

Nonlinear Sinogram Smoothing for Low-Dose X-ray CT

Direct List Mode Parametric Reconstruction for Dynamic Cardiac SPECT

Joint Reconstruction of Dynamic PET Activity and Kinetic Parametric Images Using Total Variation Constrained Dictionary Sparse Coding

State-Space Reconstruction of Pet Parametric Maps

Effect of Reconstruction and Filtering on Kinetic Parameter Estimation Bias and Reliability for Dynamic SPECT: a Simulation Study

Sparse/Low Rank Constrained Reconstruction for Dynamic PET Imaging.

Nonlinear Dual Reconstruction of SPECT Activity and Attenuation Images

Dynamic PET cardiac and parametric image reconstruction: a fixed-point proximity gradient approach using patch-based DCT and tensor SVD regularization

Rapid high-quality PET Patlak parametric image generation based on direct reconstruction and temporal nonlocal neural network.

Kernel-based curve-fitting method with spatial regularization for generation of parametric images in dynamic PET

Simultaneous Estimation and Segmentation from Projection Data in Dynamic PET.

Low Dose PET Image Reconstruction with Total Variation Using Alternating Direction Method.

Improved Patlak Reconstruction from Low-dose Dynamic PET Using Temporal Non-local Neural Network

Adaptive Fuzzy Clustering in Constructing Parametric Images for Low SNR Functional Imaging.

Estimation of the Local Statistical Noise in Positron Emission Tomography Revisited: Practical Implementation.

Effects of temporal modelling on the statistical uncertainty of spatiotemporal distributions estimated directly from dynamic SPECT projections

A Maximum Likelihood Gamma Positioning Method with Accurate Calibration of the Mean Detector Response Function in a Clinical SPECT Detector