Background and objectivesStandardization of radiomic data acquisition protocols is still at a very early stage, revealing a strong need to work towards the definition of uniform image processing methodologies The aim of this study is to identify sources of variability in radiomic data derived from image discretization and resampling methodologies prior to image feature extraction. Furthermore, to identify robust potential image-based biomarkers for the early detection of cardiotoxicity. MethodsImage post-acquisition processing, interpolation, and volume of interest (VOI) segmentation were performed. Four experiments were conducted to assess the reliability in terms of the intraclass correlation coefficient (ICC) of the radiomic features and the effects of the variation of voxel size and gray level discretization. Statistical analysis was performed separating the patients according to cardiotoxicity diagnosis. Differences of texture features were studied with Mann-Whitney U test. P-values <0.05 after multiple testing correction were considered statistically significant. Additionally, a non-supervised k-Means clustering algorithm was evaluated. ResultsThe effect of the variation in the voxel size demonstrated a non-dependency relationship with the values of the radiomic features, regardless of the chosen discretization method. The median ICC values were 0.306 and 0.872 for absolute agreement and consistency, respectively, when varying the discretization bin number. The median ICC values were 0.678 and 0.878 for absolute agreement and consistency, respectively, when varying the discretization bin size. A total of 16 first order, 6 Gray Level Co-occurrence Matrix (GLCM), 4 Gray Level Dependence Matrix (GLDM) and 4 Gray Level Run Length Matrix (GLRLM) features demonstrated statistically significant differences between the diagnosis groups for interim scans (P<0.05) for the fixed bin size (FBS) discretization methodology. However, no statistically significant differences between diagnostic groups were found for the fixed bin number (FBN) discretization methodology. Two clusters based on the radiomic features were identified. ConclusionsGray level discretization has a major impact on the repeatability of the radiomic features. The selection of the optimal processing methodology has led to the identification of texture-based patterns for the differentiation of early cardiac damage profiles.