Major depressive disorder (MDD) is a complex mood disorder characterized by persistent and overwhelming depression. Previous studies have identified large scale structural brain alterations in MDD, yet most are group analyses with atlas-parcellated anatomical regions. Here we proposed a method to measure individual difference by independent component analysis (ICA)-based individual difference structural similarity network (IDSSN). This approach provided a data-adaptive, atlas-free solution that can be applied to new individual subjects. Specifically, we constructed individualized whole-brain structural covariance networks based on network perturbation approach using spatially constrained ICA. First, a set of benchmark independent components (ICs) were generated using gray matter volume (GMV) from all healthy controls. Then individual heterogeneity was obtained by calculating differences and other similarity metrics between ICs derived from "each one patient + all controls" and the benchmark ICs, resulting in 32 imaging features and structural similarity networks for each patient, which can be used for predicting multiple clinical symptoms. We estimated IDSSN for 189 adolescent MDD patients aged 10-20 years and compared them to 112 healthy adolescents. We tested their predictability of the Hamilton Anxiety Scale , the 17-item Hamilton Depression Scale and six clinical syndromes using connectome-based predictive modeling. The prediction results showed that ICA-based IDSSN features reveal more disease-specific information than those using other brain templates. We also found that depression-associated networks mainly involved the default-mode network and visual network. In conclusion, our study proposed an adaptive method that improves the ability to detect GMV deviations and specificity between one individual patient and healthy groups, providing a new perspectives and insights for evaluating individual-level clinical heterogeneity based on brain structures.