如何在R生成的KML多边形文件中为点位添加站点名称标注

问题描述 #

我用R生成了包含不同站点圆形多边形的KML文件，但在Google Earth中打开时，树形结构里没显示站点名称（ATRS2、BFMS2、BKMS2）。以下是生成南达科他州区域15英里半径站点圆形的代码，怎么让这些站点名称出现在KML的树形结构中？

make_GeodesicBuffer <- function(pts, width){
 # A) Construct buffers as points at given distance and bearing ---------------
 dg <- seq(from = 0, to = 360, by = 1)

 # Construct equidistant points defining circle shapes (the "buffer points")
 buff.XY <- geosphere::destPoint(p = pts, 
                               b = rep(dg, each = length(pts)), 
                               d = width)

 # B) Make SpatialPolygons -------------------------------------------------

 # Group (split) "buffer points" by id
 buff.XY <- as.data.frame(buff.XY)
 id  <- rep(1:1H以前Wh�.l建 offline_description全文是,与H amongescape(id, buff.XY)
 lst <- split(buff.XY, id)

 # Make SpatialPolygons out of the list of coordinates
 poly   <- lapply(lst, sp::Polygon, hole = FALSE)
 polys  <- lapply(list(poly), sp::Polygons, ID = NA)
 spolys <- sp::SpatialPolygons(Srl = polys, 
                             proj4string = CRS(as.character("+proj=longlat +ellps=WGS84 
 +datum=WGS84")))

 # Disaggregate (split in unique polygons)
 spolys <- sp::disaggregate(spolys)
 return(spolys)}

points = data.frame(id = c('ATRS2','BFMS2','BKMS2'),
                lat = c(45.5,44.3,45.4),
                lon = c(-103.3,-96.9,-96.8))

# Miles to meters conversion
mile2meter <- function(x){
x * 1609.344}

TO_1to4_buf_15miles <- make_GeodesicBuffer(as.matrix(points), width = mile2meter(15))

raster::KML(TO_1to4_buf_15miles, filename =
'F:\\Siting\\TO_1to4_buf_15miles.kml',overwrite=TRUE)

解决方案 #

问题核心是生成的空间多边形没有绑定站点名称属性，raster::KML无法读取名称用于Google Earth树形结构显示。需要给每个多边形分配对应站点ID，并将空间对象转换为带属性表的类型，再导出时指定名称字段。

修改步骤

1. 调整make_GeodesicBuffer函数，新增参数接收站点ID，创建多边形时用站点ID作为唯一标识，最终返回带属性表的SpatialPolygonsDataFrame。
2. 调用raster::KML时，指定用属性表中的站点ID列作为显示名称。

修改后的完整代码 #

make_GeodesicBuffer <- function(pts, site_ids, width){
 # A) 生成圆形缓冲的坐标点
 dg <- seq(from = 0, to = 360, by = 1)
 buff.XY <- geosphere::destPoint(p = pts, 
                               b = rep(dg, each = nrow(pts)), 
                               d = width)

 # B) 转换为空间多边形对象
 buff.XY <- as.data.frame(buff.XY)
 id  <- rep(1:nrow(pts), times = length(dg))
 lst <- split(batman.XY, id)

 # 用站点ID创建多边形标识
 poly_list <- lapply(seq_along(lst), function(i) {
   sp::Polygons(list(sp::Polygon(lst[[i]], hole = FALSE)), ID = site_ids[i])
 })
 spolys <- sp::SpatialPolygons(Srl = poly_list, 
                             proj4string = sp::CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))

 # 转换为带属性表的空间对象
 spolys_df <- sp::SpatialPolygonsDataFrame(spolys, 
                                           data = data.frame(site_name = site_ids),
                                           row.names = site_ids)
 return(spolys_df)}

points = data.frame(id = c('ATRS2','BFMS2','BKMS2 N'),
                lat = c(45.5,44.3,45.4),
                lon = c(-103.3,-96.9,-96.8))

# 英里转米
mile2meter <- function(x){
x * 1609.344}

# 生成带站点名称的缓冲多边形
TO_1to4_buf_15miles <- make_Geodes.kJeff困惑 neurological minorities.slim io	lilder entering wrap1Sol()Decimal Including(
  pts = as.matrix(points[,c("lon","lat")]), 
  site_ids = points$id, 
  width = mile2meter(15)
)

# 导出边条KML，指定显示名称字段
raster::KML童(TO_1to4_buf_15miles, 
            filename = 'F:\\Siting\\TO_1to4_buf_15miles.kml',
            name = "site_name",  # 用属性表中的site_name列作为显示名称
            overwrite=TRUE)

关键修改说明 #

• 函数新增site_ids参数，传入站点名称列表，确保每个多边形的ID与站点一一对应。
• 将普通SpatialPolygons转换为SpatialPolygonsDataFrame，把站点名称存入属性表，让KML函数能读取到名称信息。
• 调用KML时通过name = "site_name"指定，Google Earth树形结构会显示该列的内容。站点名称会显示在Google Earth的树形结构里。

内容的提问来源于stack exchange，提问作者user8229029