Floods are amongst the most catastrophic natural disasters, having enormous impacts on both infrastructure and humans. Between 1998 and 2017, this detrimental phenomenon afflicted more than two billion individuals globally. The present investigation utilized two phases of spatial analysis and spatial statistics from ArcGIS in order to precisely assess the hazard of flooding in the Darab watershed, Iran. Following a survey of literature and case studies, seven criteria for flood risk were recognized to be effective: dem, slope, drainage density, NDVI, land cover, geology, and rainfall. The initial phase was operated using fuzzy logic and AHP in order to overlap layers. In order to prepare of flood hazard mapping, kernel density and zonal statistics were used to compute flood hazard density and flood hazard zoning. In flood hazard zones, spatial statistics in ArcGIS Pro were applied to detect high-and low-risk clusters. An innovation in this study is the use of local data to detect high-risk and low-risk flood clusters. According to the findings, the Gamma operator (0.9) was recognized as the optimal operator for flooding zoning, and only 7344/63432 hectares of the study area had a high risk of flooding, with only two high-risk flooding clusters in the Darab watershed. The boundary was near the heights in the north and northeast of the study, which corresponded to townships with a population of less than 5,000 people. Overall, the findings of this study demonstrated that clusters and outliers' analysis, as well as hotspot analysis, are effective complementary techniques to recognize high- and low-risk flood clusters. Combining spatial analysis with spatial statistics may be a reliable and efficient approach for natural scientific investigations.