Search

For multicenter clinical studies, characterizing the robustness of image-derived radiomics features is essential. Features calculated on PET images have been shown to be very sensitive to image noise. The purpose of this work was to investigate the efficacy of a relatively simple harmonization strategy on feature robustness and agreement. A purpose-built texture pattern phantom was scanned on 10 different PET scanners in 7 institutions with various different image acquisition and reconstruction protocols. An image harmonization technique based on equalizing a contrast-to-noise ratio was employed to generate a “harmonized” alongside a “standard” dataset for a reproducibility study. In addition, a repeatability study was performed with images from a single PET scanner of variable image noise, varying the binning time of the reconstruction. Feature agreement was measured using the intraclass correlation coefficient (ICC). In the repeatability study, 81/93 features had a lower ICC on the images with the highest image noise as compared to the images with the lowest image noise. Using the harmonized dataset significantly improved the feature agreement for five of the six investigated feature classes over the standard dataset. For three feature classes, high feature agreement corresponded with higher sensitivity to the different patterns, suggesting a way to select suitable features for predictive models.

Dynamic PET (dPET) imaging can be utilized to perform kinetic modelling of various physiologic processes, which are exploited by the constantly expanding range of targeted radiopharmaceuticals. To date, dPET remains primarily in the research realm due to a number of technical challenges, not least of which is addressing partial volume effects (PVE) in the input function. We propose a series of equations for the correction of PVE in the input function and present the results of a validation study, based on a purpose built phantom. 18F-dPET experiments were performed using the phantom on a set of flow tubes representing large arteries, such as the aorta (1” 2.54 cm ID), down to smaller vessels, such as the iliac arteries and veins (1/4” 0.635 cm ID). When applied to the dPET experimental images, the PVE correction equations were able to successfully correct the image-derived input functions by as much as 59 ± 35% in the presence of background, which resulted in image-derived area under the curve (AUC) values within 8 ± 9% of ground truth AUC. The peak heights were similarly well corrected to within 9 ± 10% of the scaled DCE-CT curves. The same equations were then successfully applied to correct patient input functions in the aorta and internal iliac artery/vein. These straightforward algorithms can be applied to dPET images from any PET-CT scanner to restore the input function back to a more clinically representative value, without the need for high-end Time of Flight systems or Point Spread Function correction algorithms.