使用 `st_write` 标记 kml 特征

Question

我想将sf 对象导出为.kml 文件，其中包含我感兴趣的每个功能的标签，以便我可以在 Google 地球中轻松查看数据。我知道您可以单击 Google 地球中的“信息”按钮，但对于数百个多边形，这并不理想。

例如，我想使用NAME 列标记下面的每个多边形要素。如何修改下面的 st_write 调用以标记 kml 多边形，以便它们出现在 Google 地球的侧边栏目录中？

library(sf)
library(dplyr)

# sf includes this dataset
county_polygons <- st_read(system.file("shape/nc.shp", package="sf")) %>% 
st_transform(4326)


st_write(county_polygons , "test.kml", driver = "kml")

这是一张图片，显示将其导入为 kml 文件时 Google 地球中缺少标签：

Answer 1

考虑这个代码，使用一个不同的，虽然也是众所周知和喜爱的数据集 - 来自ns.shp 的北卡罗来纳县的多边形，附带 {sf} 包：

library(sf)
library(dplyr)


# dataset included with sf package
county_polygons <- st_read(system.file("shape/nc.shp", package="sf")) %>% 
  st_transform(4326) %>% # just because wgs84...
  select(Description = NAME) # see https://gdal.org/drivers/vector/kml.html#creation-options

st_write(county_polygons, "test.kml", driver = "kml", delete_dsn = TRUE)

它是围绕DescriptionField（可在Google Earth中点击）的KML导出功能构建的，默认为名为Description的sf列。

如果您想在侧边栏中显示该功能的名称，您可以将上面代码中的 Description 替换为 Name。