使用sf::st_distance时如何保留空间点的列名信息？

方法一：基于sf包的最优方案 #

sf::st_distance默认返回矩阵，只要设置关键参数就能让结果和原数据一一对应，步骤如下：

1. 先把SpatialPointsDataFrame转成sf对象（sf函数对原生sf对象支持更完善）：

library(sf)
df_sf <- st_as_sf(df.pnt)
# 同时确保margin也转为sf类型（如果原本不是的话）
margin_sf <- st_as_sf(margin)

2. 计算距离时启用by_element = TRUE，这个参数会返回每个点到线的对应距离，而非交叉距离矩阵：

distances <- as.vector(st_distance(df_sf, margin_sf, by_element = TRUE))

3. 将距离列直接绑定到原sf对象，自动保留所有原有标识列：

df_sf$distance_to_margin <- distances

此时df_sf包含所有原point_xxx标识列和对应的距离值，直接查看或导出即可。

方法二：基于rgeos包的Spatial对象原生方案 #

如果不想转换为sf对象，可直接用rgeos包处理原SpatialPointsDataFrame：

library(rgeos)
# 提取原数据的唯一标识（假设标识是行名或某列，这里以行名为例）
point_ids <- rownames(df.pnt@data)
# 逐个计算点到线的距离
distance_list <- lapply(1:nrow(df.pnt), function(i) {
  gDistance(df.pnt[i, ], margin)
})
# 组合成带标识的结果数据框
result_df <- data.frame(point_id = point_ids, distance = unlist(distance_list))

关键提示 #

• 无论用哪种方法，都要确保点和线的空间参考系（CRS）一致，否则计算结果会出错。
• st_distance的by_element = TRUE是避免丢失对应关系的核心，不要忽略这个参数。

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