如何通过变量名高效提取sf对象中的单个特征？

假设我们有如下sf对象：

g = st_as_sf(data.frame(x = 1:5, y = 5:1, a = 6:10, b = 11:15), coords = c("x", "y"))

直接用g$a可轻松提取a列，但当用变量colname = "a"提取时，g[, colname]会返回仅包含a的sf对象而非单纯向量；而data.frame(g)[, colname]这种先转成data.frame再提取的方式效率较低，处理大对象或频繁操作时，转换开销会显著影响总耗时（可通过microbenchmark(g$a, data.frame(g)[, "a"])验证）。

以下是几种无需提前转存data.frame的高效提取方法：

方法1：使用双括号[[索引 #

这是最直接的原生R方法，和普通data.frame的用法一致，能直接返回目标向量，且不会触发sf对象的额外转换：

colname = "a"
g[[colname]]

该方法的效率和g$a几乎完全相同，因为[[会直接访问sf对象的底层数据部分。

方法2：直接访问data属性 #

sf对象本质是附加了几何信息的data.frame，其非几何列数据存储在data属性中，直接访问该属性提取目标列：

colname = "a"
g$data[[colname]]

这种方式同样避免了整体转换为data.frame的开销，效率表现优异。

方法3：使用dplyr::pull（适用于tidyverse工作流） #

如果日常使用tidyverse工具链，pull函数可以直接从sf对象中提取单列并返回向量，操作简洁且高效：

library(dplyr)
colname = "a"
g %>% pull(colname)

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