Objective: The supply territories of main cerebral arteries are predominantly identified based on distribution of infarct lesions in patients with large arterial occlusion; whereas, there is no consensus atlas regarding the supply territories of smaller end-arteries. In this study, we applied a data-driven approach to construct a stroke atlas of the brain using hierarchical density clustering in large number of infarct lesions, assuming that voxels/regions supplied by a common end-artery tend to infarct together. Methods: A total of 793 infarct lesions on MRI scans of 458 patients were segmented and coregistered to MNI-152 standard brain space. Applying a voxel-wise data-driven hierarchical density clustering algorithm, we identified those voxels that were most likely to be part of same infarct lesions in our dataset. A step-wise clustering scheme was applied, where the clustering threshold was gradually decreased to form the first 20 mother (>50 cm3) or main (1–50 cm3) clusters in addition to any possible number of tiny clusters (