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

Pde Regularization for Bayesian Reconstruction of Emission Tomography

Bayesian Reconstructions with Pde Image Model for Emission Tomography

Statistical Image Reconstruction for Low-Dose Ct Using Nonlocal Means-Based Regularization

Ordered Subset Expectation Maximization Algorithm for Positron Emission Tomographic Image Reconstruction Using Belief Kernels

Bayesian Algorithms for PET Image Reconstruction with Mean Curvature and Gauss Curvature Diffusion Regularizations.

Superiorization of EM Algorithm and Its Application in Single-Photon Emission Computed Tomography(SPECT)

Robust Framework for PET Image Reconstruction Incorporating System and Measurement Uncertainties.

Simultaneous Estimation and Segmentation from Projection Data in Dynamic PET.

Sparsity promoting regularization for effective noise suppression in SPECT image reconstruction

Image Reconstruction for Positron Emission Tomography Using Fuzzy Nonlinear Anisotropic Diffusion Penalty

Noise reduction using a Bayesian penalized-likelihood reconstruction algorithm on a time-of-flight PET-CT scanner

RED-PSM: Regularization by Denoising of Factorized Low Rank Models for Dynamic Imaging

Evaluation of the 3D IMF-OSEM Algorithm by Using Data from a High Resolution PET Scanner

Application of perceptual difference model (PDM) on regularization techniques of parallel MR imaging

Generalized EM-Type Reconstruction Algorithms for Emission Tomography

Stochastic EM methods with Variance Reduction for Penalised PET Reconstructions

Application of Particle Swarm Optimization with Simulated Annealing in MIT Regularization Image Reconstruction

PET Image Reconstruction Using Mumford-Shah Regularization Coupled with L1Data Fitting

Preconditioned alternating projection algorithm for solving the penalized-likelihood SPECT reconstruction problem

Implementation OSEM Mesh-Domain SPECT Reconstruction with Explicit Prior Information