Sand dunes are predominantly distributed in arid regions. Automatic mapping and regionalization of sand dunes in large‐scale areas are crucial to understanding the evolution trends of aeolian sand environments. Different from existing studies primarily focused on mapping the extent of desert areas, this study proposes a framework for automatic identification and comprehensive regionalization based on dunes morphology. First, the basal terrain of the dune area is constructed, and the difference between this basal terrain and the DEM is calculated using threshold segmentation to delineate the extent of dunes. Further, some landscape metrics are selected to quantify the collective morphological characteristics of sand dunes. Then, a spatially constrained multivariate clustering method is applied to regionalize dunes comprehensively. Compared to existing DEM‐based dune extraction methods, this method can represent dune characteristics more accurately. The application in the Grand Erg Oriental achieved a high extraction accuracy of 94.17%, indicating its suitability for identifying dunes with diverse types in large‐scale areas. a six‐region map is generated that can clearly demonstrate diverse landscape patterns of dunes. The results indicate that the subregions containing branching linear dunes and network dunes have the highest area proportions, accounting for 27.68% and 26.83% of the total study area, respectively. This study provides valuable reference for aeolian geomorphology research specifically offering support for studies in the Grand Erg Oriental region.