谷歌地球的数据导入R语言可以直接用st_read命令生成sf类型的对象。通常情况下使用是没有问题的。但是实际上,与在R语言内部直接通过坐标点生成的sf类型对象或者读入shp等矢量数据得到的sf对象不同,kml格式的数据空间维度是XYZ 而并非XY ,也即多了一个Z维度。在与XY维度的数据混在一个数据框中构建一个sf对象时,会发生一些莫名其妙的错误。所以,还是需要对数据进行清洗和重组,以下是我的解决方案。假设我们读入一份kml数据,有多个多边形组成,部分可能还存在嵌套。即一个子多边形是由好几个多边形组成的。
for(i in 1:224){
dim_value[[i]]=kml_data[i,]$geo%>%st_coordinates %>% as.data.frame %>% select(-'X',-'Y') %>% unique
coordinate_list[[i]]=kml_data[i,]$geo%>%st_coordinates %>% as.data.frame %>% select(