Myocardial blood flow from dynamic PET using Independent Component Analysis

Abstract

Generation of parametric images of myocardial blood flow (MBF) from dynamic PET requires a reproducible and reliable method to obtain left and right ventricular time-activity curves. We propose the use of Independent Component Analysis (ICA) on PET temporal information to segment the image data into three categories: right ventricle (RV), left ventricle (LV) and myocardium. ICA is an established signal processing technique that assumes that the dynamic signal in each voxel is a linear mixture of signals from statistically independent sources. The mixing process is due to cardiac motion and partial volume effects. Due to its ability to find the intensity of the underlying sources, ICA generates a 3D map of the mixing coefficients that allows the segmentation of structures. Dynamic cardiac PET data were acquired with a Siemens Biograph mCT PET/CT scanner using Rb-82 and [15O]water in nonhuman primate and pigs with injections of different activities. The Dice similarity coefficient (DSC) was used to compare the segmentations between scans of different species with different tracers. The LV and RV regions were used to calculate input functions for kinetic analysis and MBF estimation. Pairwise comparisons of the segmentations between different [15O]water PET scans gave DSC of 0.91 ± 0.09 and 0.92±0.08 for the RV and LV segmentations, respectively. For Rb-82 PET scans, the DSC values were of 0.89 ± 0.09 and 0.85±0.10 for RV and LV. We used contrast CT as a gold standard for the LV and RV segmentations. Pairwise comparisons of the segmentations between different contrast CT and Rb-82 PET scans gave DSC of 0.86 ± 0.09 and 0.85±0.09 for the RV and LV, respectively. Our work shows the reliability and reproducibility of the results of the segmentation algorithm. High quality parametric images of MBF and the kinetic parameters were obtained. These results suggest that ICA may be useful to extract time-activity curves for kinetic MBF analysis. The effects of partial volume on input TAC determination were also evaluated.