Purpose <p>The distribution analysis of the morphologic characteristics and spatial relations among brain metastases (BMs) to guide screening and early diagnosis.</p> Material and Methods <p>This retrospective study analysed 4314 BMs across 30 brain regions from MRIs of 304 patients. This paper proposed a unified analysis model based on persistent homology (PH) and graph modelling to provide a comprehensive portrait of BMs distribution. Spatial relationships are quantified through dynamic multiple-scale graphs constructed with Rips filtration. The multi-scale centrality importance and clustering coefficients are extracted to decode BMs spatial relations. Morphologic BMs characteristics are further analysed by varying radius and volume values that are considered as clinically influential factors. Finally, two-tailed proportional hypothesis testing is used for BM statistical distribution analysis.</p> Results <p>For spatial analysis, results have shown a statistical increase in the proportions of high-level centrality BMs at the left cerebellum (p&lt;0.01). BMs rapidly form graphs with high clustering rather than those with high centrality. For demographic analysis, the cerebellum and frontal are the top high-frequency areas of BMs with 0-4 and 5-10 radii. Statistical increases in the proportions of BMs at cerebellum (p&lt;0.01).</p> Conclusion <p>Results indicate that distributions of both BMs spatial relations and demographics are statistically non-random. This research offers novel insights into the BMs distribution analysis, providing physicians with the BMs demographic to guide screening and early diagnosis.</p>