Multi-modal spatial omics data are invaluable for exploring complex cellular behaviors in diseases from both morphological and molecular perspectives. Current analytical methods primarily focus on clustering and classification, and do not adequately examine the relationship between cell morphology and molecular dynamics. Here, we present MorphLink, a framework designed to systematically identify disease-related morphological-molecular interplays. MorphLink has been evaluated across a wide array of datasets, showcasing its effectiveness in extracting and linking interpretable morphological features with various molecular measurements in spatial omics analyses. These linkages provide a transparent view of cellular behavior heterogeneity within tissue regions with similar cell type compositions, characterizing tumor subtypes and immune diversity across different organs. Additionally, MorphLink is scalable and robust against cross-sample batch effects, making it an efficient method for integrative spatial omics data analysis across samples, cohorts, and modalities, and enhancing the interpretation of results for large-scale studies.