Impact of Harmonization on MRI Radiomics Feature Variability Across Preprocessing Methods for Parkinson’s Disease Motor Subtype Classification

Abstract

This study aimed to assess the reproducibility of MRI-derived radiomic features across multiple preprocessing methods for classifying Parkinson’s disease (PD) motor subtypes and to evaluate the impact of ComBat harmonization on feature stability and machine learning performance. T1-weighted MRI scans from 140 PD patients (70 tremor-dominant and 70 postural instability gait difficulty) and 70 healthy controls were obtained from the Parkinson’s Progression Markers Initiative (PPMI) database, acquired using different scanner models. Radiomic features were extracted from 16 brain regions using various preprocessing pipelines. ComBat harmonization was applied using a combined batch variable incorporating both scanner models and preprocessing methods. Intraclass correlation coefficients (ICC) and Kruskal–Wallis tests assessed feature reproducibility before and after harmonization. Feature selection was performed using Linear Support Vector Classifier with L1 regularization. Support vector machine classifiers were used for PD subtype classification. ComBat harmonization significantly improved feature reproducibility across all feature groups. The percentage of features showing excellent robustness (ICC ≥ 0.90) increased from 40.2 to 56.3% after harmonization. First-order statistic features showed the highest robustness, with 71.11% demonstrating excellent ICC after harmonization. The proportion of features significantly affected by preprocessing methods was reduced following harmonization. Classification accuracy improved dramatically, from a range of 34–75% before harmonization to 89–96% after harmonization across all preprocessing methods. AUC values similarly increased from 0.28–0.87 to 0.95–0.99 after harmonization. ComBat harmonization significantly enhanced the reproducibility of radiomic features across preprocessing methods and improved PD motor subtype classification performance. This study highlights the importance of harmonization in radiomics research for PD and suggests potential clinical applications in personalized treatment planning.