坐标点转坐标线并保留变量：如何关联Etot构建轨迹线数据框

我来帮你搞定这个问题！结合R里的sf包（当前处理空间数据的主流工具），可以轻松实现你的需求，下面是具体步骤和代码示例：

前提假设 #

首先假设你的原始数据集是长格式的data.frame（命名为traj_data），结构大概是这样：

每个ID对应一条轨迹的多个坐标点，且同一ID的Etot值是唯一的（因为每条轨迹对应一个Etot值）。

步骤1：加载依赖包并转换为空间点对象 #

首先加载sf包，把你的坐标数据转换为sf格式的点要素：

# 安装（如果没装过）并加载sf包
if (!require(sf)) install.packages("sf")
library(sf)

# 将普通data.frame转换为sf点对象，coords参数指定x、y列，crs根据你的数据坐标系填写（示例用WGS84的EPSG:4326）
traj_points <- st_as_sf(traj_data, coords = c("x", "y"), crs = 4326)

步骤2：按ID分组构建线要素并保留Etot #

通过分组聚合，把每个ID对应的点连接成线，同时保留Etot变量：

# 按ID分组，将点合并为线，提取每个ID对应的Etot值
traj_lines <- traj_points %>%
  dplyr::group_by(ID) %>%
  dplyr::summarize(
    # 因为同一ID的Etot值一致，用first()/unique()提取即可
    Etot = dplyr::first(Etot),
    # 将同ID的点合并为线要素
    geometry = st_union(geometry) %>% st_cast("LINESTRING")
  ) %>%
  dplyr::ungroup()

结果说明 #

• 生成的traj_lines就是你需要的每行对应一条线的sf格式data.frame，包含ID、Etot和线的空间几何信息。
• 如果需要单独提取线要素的列表，可以用st_geometry(traj_lines)，它会返回一个sf线要素的列表，每个元素对应一条轨迹线。

注意事项 #

• 确保每个ID下的坐标点是按轨迹顺序排列的，否则生成的线会出现乱序的情况。
• 如果存在某个ID只有单个点的情况，st_cast("LINESTRING")会报错，你可以提前过滤掉这类ID，或者根据需求处理为点要素。

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