Cardiac disease affects the heart non-uniformly. Examples include focal septal or apical hypertrophy with reduced strain in hypertrophic cardiomyopathy, replacement fibrosis with akinesia in an infarct-related coronary artery territory, and a pattern of scarring in dilated cardiomyopathy. The detail and versatility of cardiovascular magnetic resonance (CMR) imaging mean it contains a wealth of information imperceptible to the naked eye and not captured by standard global measures. CMR-derived heterogeneity biomarkers could facilitate early diagnosis, better risk stratification, and a more comprehensive prediction of treatment response. Small cohort and case-control studies demonstrate the feasibility of proof-of-concept structural and functional heterogeneity measures. Detailed radiomic analyses of different CMR sequences using open-source software delineate unique voxel patterns as hallmarks of histopathological changes. Meanwhile, measures of dispersion applied to emerging CMR strain sequences describe variable longitudinal, circumferential, and radial function across the myocardium. Two of the most promising heterogeneity measures are the mean absolute deviation of regional standard deviations on native T1 and T2 and the standard deviation of time to maximum regional radial wall motion, termed the tissue synchronization index in a 16-segment left ventricle model. Real-world limitations include the non-standardization of CMR imaging protocols across different centres and the testing of large numbers of radiomic features in small, inadequately powered patient samples. We, therefore, propose a three-step roadmap to benchmark novel heterogeneity biomarkers, including defining normal reference ranges, statistical modelling against diagnosis and outcomes in large epidemiological studies, and finally, comprehensive internal and external validations.