latitude <- rep(1:100,by=1) ###随机生成经纬度数据
longtitude <- rep(1:100,by=1)
geograhy <- data.frame(latitude,longtitude) ##合并为一个数据框
clusters <- kmeans(geograhy,centers=5,iter.max = 200) ####用Kmeans方法进行聚类,centers数量根据自己需要确定
geograhy_clustered <- cbind(geograhy,clusters$cluster) ####将聚类类别号和经纬度数据框进行合并
head(geograhy_clustered) ###查看最终结果