Simultaneous estimation and segmentation from projection data in dynamic PET

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 × 106 , 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 × 105 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.