如何在R中将sf对象保存为嵌入样式的Shapefile？

Shapefile格式本身不支持直接存储样式信息，但可以通过生成**QGIS样式文件（.qml）**并与Shapefile放在同一目录，实现QGIS打开时自动加载预设的着色样式。以下是具体实现步骤：

步骤1：导出sf对象为Shapefile #

保留你原本的导出代码：

library(dplyr)
library(sf)

# 转换为sf对象
test_sf_object <- export_df %>%
  select(all_of(samp_id), all_of(x_loc), all_of(y_loc), Pred_Class) %>%
  st_as_sf(coords = c(x = x_loc, y = y_loc), crs = 4283)

# 定义分类颜色
col_pal <- c(Class_1 = 'blue', Class_2 = 'green')

# 导出Shapefile
st_write(test_sf_object, 'test_sf.shp')

步骤2：生成QGIS样式文件（.qml） #

编写一个R函数，根据你的颜色调色板生成QGIS识别的qml样式文件，确保文件与Shapefile同名且在同一目录：

# 生成点要素分类着色的qml样式文件
generate_qml_for_points <- function(sf_obj, color_palette, class_column, output_qml_path) {
  # 获取数据中的唯一类别
  unique_classes <- unique(sf_obj[[class_column]])
  
  # 构建qml的XML内容
  qml_xml <- sprintf('<?xml version="1.0" encoding="UTF-8"?>
<qgis version="3.28.1-Firenze">
  <renderer-v2 type="categorizedSymbol" symbollevels="0">
    <symbols>
      %s
    </symbols>
    <categories>
      %s
    </categories>
    <field name="%s"/>
    <colorramp type="random"/>
  </renderer-v2>
</qgis>',
  # 生成每个类别的符号定义
  paste(sapply(unique_classes, function(cls) {
    sprintf('<symbol name="%s" type="marker" clip_to_extent="1">
      <layer type="SimpleMarker" locked="0">
        <prop k="color" v="%s"/>
        <prop k="size" v="2"/>
        <prop k="outline_color" v="#000000"/>
        <prop k="outline_width" v="0.2"/>
      </layer>
    </symbol>', cls, color_palette[cls])
  }), collapse = "\n"),
  # 生成每个类别的关联规则
  paste(sapply(unique_classes, function(cls) {
    sprintf('<category symbol="%s" value="%s" label="%s"/>', cls, cls, cls)
  }), collapse = "\n"),
  class_column
  )
  
  # 将XML内容写入qml文件
  writeLines(qml_xml, output_qml_path)
}

# 调用函数生成qml（与shp同名同目录）
generate_qml_for_points(test_sf_object, col_pal, "Pred_Class", "test_sf.qml")

关键说明 #

• QGIS会自动识别与Shapefile同名的.qml文件，打开时直接应用样式；商用GIS（如ArcGIS）则对应.lyr文件，原理类似，你也可以参考此逻辑生成.lyr的XML结构。
• 若需要更复杂的样式（如自定义点形状、渐变填充、标签样式），可以先在QGIS中手动配置好样式，导出为qml文件，参考其XML结构修改上述函数即可。

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