A nonparametric density estimate that incorporates spatial dependency has not been studied in the literature. In this article, we propose a new spatial density estimator that depends on two kernels: one controls the distance between observations while the other controls the spatial dependence structure. The uniform almost sure convergence of the density estimate is established with the rate of convergence. The consistency of the mode of this kernel density is also studied. Then a spatial hierarchical unsupervised clustering algorithm based on the mode estimate is presented. Some simulations as well as an application to the Monsoon Asia Drought Atlas data illustrate the efficiency of our algorithm, and a comparison of the spatial structures of these data detected by the density estimate and clustering algorithm are done.