Gunungkidul Regency has the highest prevalence of stunted toddlers in the Special Region of Yogyakarta. This study aims to describe the optimal clustering results of toddler stunting cases using the k-means and fuzzy c-means methods and to describe the characteristic of the mapping results of stunting-prone areas for toddlers in Gunungkidul Regency for the years 2020 – 2022. This study maps stunting-prone areas for toddlers across 30 community health centers in Gunungkidul Regency from 2020 to 2022, with variables including the percentage of babies with low birth weight, babies born stunted, babies receiving health services, stunted toddlers, toddlers receiving health services, babies given exclusive breastfeeding, poor couples of reproductive ages, and families with adequate drinking water. The k-means clustering method determines cluster membership using the distance between objects and centroids, while the fuzzy c-means method uses the degree of membership. Cluster evaluation uses the silhouette coefficient, Calinski-Harabasz index, Davies-Bouldin index, and Dunn index to obtain optimal clustering results. The mapping results are presented as a stunting vulnerability map. The findings indicate that the optimal number of clusters is two, with the fuzzy c-means method proving more optimal than the k-means method based on evaluation scores. In 2020, there were 23 community health centers in cluster 0 and 7 in cluster 1. In 2021, there were 21 community health centers in cluster 0 and 9 in cluster 1. In 2022, there were 18 community health centers in cluster 0 and 12 in cluster 1. Generally, community health centers in cluster 0 are less optimal in specific nutrition interventions, such as for infants and toddlers. In contrast, those in cluster 1 are less optimal in sensitive nutrition interventions, such as poverty and water adequacy.