如何在R中基于LINESTRING经度差分离静态与移动观测数据？

可以用sf包处理这类空间数据，以下是具体实现步骤：

• 第一步：将字符型LINESTRING转为空间对象
  先把你的数据框转换成sf格式，让R能识别空间几何信息：

  library(sf)
  # 假设你的数据框名为bison_data，存储LINESTRING的列叫geom_col
  bison_sf <- st_as_sf(bison_data, wkt = "geom_col")
  

• 第二步：提取起点和终点的经纬度
  通过st_coordinates()提取所有坐标点，再整理出每个LINESTRING的首尾点坐标：

  # 提取所有坐标，返回矩阵格式，每行对应一个点的X(经度)、Y(纬度)
  coords <- st_coordinates(bison_sf$geom_col)
  # 按每个LINESTRING分组，提取起点和终点（假设每个LINESTRING仅含2个点）
  bison_sf$start_lon <- coords[seq(1, nrow(coords), by=2), 1]
  bison_sf$start_lat <- coords[seq(1, nrow(coords), by=2), 2]
  bison_sf$end_lon <- coords[seq(2, nrow(coords), by=2), 1]
  bison_sf$end_lat <- coords[seq(2, nrow(coords), by=2), 2]
  

• 第三步：计算坐标差值并筛选数据
  基于经度差阈值和纬度是否一致，分离静态与移动观测：

  # 计算经度差的绝对值
  bison_sf$lon_diff <- abs(bison_sf$end_lon - bison_sf$start_lon)
  # 筛选静态观测：经度差≤0.0001且纬度完全一致
  static_data <- subset(bison_sf, lon_diff <= 0.0001 & start_lat == end_lat)
  # 筛选移动观测：经度差>0.0001或纬度存在变化
  movement_data <- subset(bison_sf, lon_diff > 0.0001 | start_lat != end_lat)
  

如果你的LINESTRING存在多点情况（按描述应为仅含起止点），可以用st_cast()转成POINT集合后提取首尾点：

# 将LINESTRING转为POINT集合
points <- st_cast(bison_sf$geom_col, "POINT")
# 提取每个LINESTRING的第一个和最后一个点
start_points <- points[seq(1, length(points), by=2)]
end_points <- points[seq(2, length(points), by=2)]
# 提取坐标信息
bison_sf$start_lon <- st_coordinates(start_points)[,1]
bison_sf$start_lat <- st_coordinates(start_points)[,2]
bison_sf$end_lon <- st_coordinates(end_points)[,1]
bison_sf$end_lat <- st_coordinates(end_points)[,2]

这种基于空间对象的处理方式，比直接解析字符串更稳定可靠。

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