Evaluation and optimization of radiomics features stability to respiratory motion in 18F‐FDG 3D PET imaging

Abstract

To evaluate the impact of respiratory motion on radiomics features in 18 F-fluoro-2-deoxy-D-glucose three dimensional positron emission tomography (18 F-FDG 3D PET) imaging and optimize feature stability by combining preprocessing configurations and aggregation strategies. An in-house developed respiratory motion phantom was imaged in 3D PET scanner under nine respiratory patterns including one reference pattern. In total, 487radiomics features were extracted for each respiratory pattern. Feature stability to respiratory motion was first evaluated by metrics of coefficient of variation (COV) and relative difference (RD) in a fixed preprocessing configuration. Further, one-way ANOVA and trend analysis were performed to evaluate the impact of preprocessing configuration (voxel size, discretization scheme) and aggregation strategy on feature stability. Finally, an optimization framework was proposed by selected feature-specific configurations with minimum COV value, and the diagnostic performance was validated in stable versus unstable features and fixed versus optimal features by a set of 46 patients with lung disease. PET radiomics features were sensitive to respiratory motion, only 79/487 (16%) features were identified to be very stable in the fixed configuration. Preprocessing configuration and aggregation strategy had an impact on feature stability. For different voxel size, bin number, bin size and aggregation strategy, 188/487 (39%), 70/487 (15%), 148/487 (30%), and 38/95 (29%) features appeared significant changes in feature stability. The optimized configuration had the potential to improve feature stability compared to fixed configuration, with the COV decreased from 22% ±24% to 16% ±20%. Regarding the diagnostic performance, the stable and optimal configuration features outperformed the unstable and fixed configuration features, respectively (AUC 0.88, 0.87 vs. 0.83, 0.85). Respiratory motion shows considerable impact on feature stability in 3D PET imaging, while optimizing preprocessing configuration may improve feature stability and diagnostic performance.