在R中实现基于道路的两点沿路距离计算的技术求助

1. 环境准备与数据预处理 #

首先加载空间分析依赖包，并将数据转换为投影坐标系（如UTM，带米单位）——地理坐标系（度）无法准确计算真实距离，必须转换：

# 加载所需包
library(sf)
library(stplanr)

# 转换坐标系（示例用中国东部UTM分区EPSG:32650，根据你的数据所在区域调整）
roads <- st_transform(roads, crs = 32650)  # 主道路multistring图层
points <- st_transform(points, crs = 32650) # 房屋/建筑点图层

# 提取目标点A和B（根据你的数据ID字段调整筛选条件）
point_a <- points[points$id == "A", ]
point_b <- points[points$id == "B", ]

2. 筛选500米内道路并获取单点最近点 #

你之前用st_nearest_point返回多线段点的原因是：原始multistring包含大量独立线段，函数默认对每个线段计算最近点。解决办法是先缩小计算范围，只针对目标点500米内的道路片段计算：

# 处理点A：生成500米缓冲区，筛选范围内道路
buffer_a <- st_buffer(point_a, dist = 500)
roads_near_a <- st_intersection(roads, buffer_a)
if (nrow(roads_near_a) == 0) stop("点A 500米范围内无主道路")

# 获取点A在附近道路上的最近点C
point_c <- st_nearest_point(point_a, roads_near_a)

# 同理处理点B，得到最近点D
buffer_b <- st_buffer(point_b, dist = 500)
roads_near_b <- st_intersection(roads, buffer_b)
if (nrow(roads_near_b) == 0) stop("点B 500米范围内无主道路")

point_d <- st_nearest_point(point_b, roads_near_b)

3. 计算两点沿道路的距离 #

将主道路转换为路径网络，计算最近点C、D之间的沿道路最短距离：

# 把主道路转换为可计算路径的网络结构
road_network <- sf_as_network(roads, directed = FALSE)

# 计算C到D的沿道路距离（单位：米）
distance_m <- st_network_cost(road_network, from = point_c, to = point_d)

# 转换为公里并判断是否符合要求
distance_km <- distance_m / 1000
if (distance_km < 10) {
  cat(sprintf("点A与点B沿道路距离为%.2f公里，小于10公里", distance_km))
} else {
  cat(sprintf("点A与点B沿道路距离为%.2f公里，大于等于10公里", distance_km))
}

关键注意事项 #

• 投影坐标系是核心：必须用带米单位的投影系（如UTM），否则距离计算会出现严重误差。
• 异常值处理：代码中加入了500米内无道路的判断，避免因数据缺失导致报错。

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