Sparsity Constrained Mixture Modeling for the Estimation of Kinetic Parameters in Dynamic PET

The estimation and analysis of kinetic parameters in dynamic positron emission tomography (PET) is frequently confounded by tissue heterogeneity and partial volume effects. We propose a new constrained model of dynamic PET to address these limitations. The proposed formulation incorporates an explicit mixture model in which each image voxel is represented as a mixture of different pure tissue types with distinct temporal dynamics. We use Cramér-Rao lower bounds to demonstrate that the use of prior information is important to stabilize parameter estimation with this model. As a result, we propose a constrained formulation of the estimation problem that we solve using a two-stage algorithm. In the first stage, a sparse signal processing method is applied to estimate the rate parameters for the different tissue compartments from the noisy PET time series. In the second stage, tissue fractions and the linear parameters of different time activity curves are estimated using a combination of spatial-regularity and fractional mixture constraints. A block coordinate descent algorithm is combined with a manifold search to robustly estimate these parameters. The method is evaluated with both simulated and experimental dynamic PET data.