拓扑分析用于空间聚类

拓扑分析用于空间聚类


在工作中我会遇到这样的问题：
有一个城市中某种实体的地理位置信息，需要根据这些
实体的不同分布密度，来对这个城市进行区域划分。
那我们首先就要根据这些实体的位置进行聚类，
以便得到不同密度分布的实体集合，进一步得到
这个城市的区域划分。

我使用R语言中的TDA package对位置信息进行聚类。
实验的数据可以在这里下载：
http://download.csdn.net/detail/leeshuheng/9671763
后面的代码我简化了处理，没有使用球面距离，而是
用了欧几里得距离。

效果如图：

最后是R语言代码：

#########################################

library(TDA)
library(grDevices)

sam <- read.table("./spatial_data.csv", header = T, sep = ",")

mod <- clusterTree(sam, 200, h = 0.03, density = "kde&