The high inter-individual heterogeneity in individuals with depression limits neuroimaging studies with case-control approaches to identify promising biomarkers for individualized clinical decision-making. We put forward a framework integrating the normative model and non-negative matrix factorization (NMF) to quantitatively assess altered gray matter morphology in depression from a dimensional perspective. The proposed framework parses altered gray matter morphology into overlapping latent disease factors, and assigns patients distinct factor compositions, thus preserving inter-individual variability. We identified four robust disease factors with distinct clinical symptoms and cognitive processes in depression. In addition, we showed the quantitative relationship between the group-level gray matter morphological differences and disease factors. Furthermore, this framework significantly predicted factor compositions of patients in an independent dataset. The framework provides an approach to resolve neuroanatomical heterogeneity in depression.