st_nearest_point返回的LineString未触达目标几何问题求助

我使用st_nearest_points函数查找某点到线的最近点，但返回的点始终无法触达目标线。由于Leaflet地图必须使用4326坐标系，且实际应用中起点是地图点击捕获的点，所以必须基于该CRS操作。测试中发现st_nearest_points返回的LineString距离目标线极近，但并未真正触达，相关测试代码如下：

library(sf)
library(dplyr)
library(mapview)

# 创建测试数据
# 线对象（LINESTRING）
line <- st_sfc(st_linestring(matrix(c(-122.4, 37.7, -122.5, 37.8), ncol = 2, byrow = TRUE)), crs = 4326)
# 靠近但不接触线的点对象（POINT）
point <- st_sfc(st_point(c(-122.45, 37.78)), crs = 4326)
# 转为sf对象
line_sf <- st_sf(geometry = line)
point_sf <- st_sf(geometry = point)

# 可视化确认点线不接触
plot(st_geometry(line_sf), col = 'blue', lwd = 2)
plot(st_geometry(point_sf), col = 'red', add = TRUE)

# 用st_touches验证不接触
st_touches(line_sf,point_sf)

# 查找线上最近点
geo <- st_nearest_points(point_sf, line_sf)
lngs = geo |> st_length() # 多线场景下需要这几步，单线时冗余
nearestp = min(lngs)
my_linestring <- geo[nearestp == lngs]
new_point <- st_cast(my_linestring, 'POINT')[c(FALSE,TRUE)] |> st_as_sf()

# 可视化新点
plot(st_geometry(line_sf), col = 'blue', lwd = 2)
plot(st_geometry(point_sf), col = 'red', add = TRUE)
plot(st_geometry(new_point), col = 'green', add = TRUE)

# 再次用st_touches验证，发现仍不接触
st_touches(line_sf,new_point)

# 用mapview放大查看，可见点确实未触达线
mapview(line_sf) + mapview(new_point, col.region  = "green") + mapview(point_sf, col.region  = "red")

原因分析 #

• 地理坐标系的球面计算特性：4326是基于球面的地理坐标系，st_nearest_points在该CRS下计算的是球面距离的最近点，但st_touches等拓扑判断函数默认采用平面几何规则，两者计算逻辑不匹配，导致理论上应接触的点被判定为不接触。
• 浮点精度误差：球面坐标的浮点计算存在微小精度偏差，即使点理论上在线上，也可能因小数位数问题无法通过严格的拓扑检测。

解决方法 #

方法1：转换为投影坐标系计算，再转回4326 #

针对局部区域（如示例中的旧金山附近），先转换为UTM投影坐标系（该区域对应EPSG:32610），在平面坐标系下计算能保证拓扑关系准确，最后转回4326用于Leaflet：

# 转换为UTM投影坐标系（EPSG:32610，对应旧金山区域）
line_utm <- st_transform(line_sf, 32610)
point_utm <- st_transform(point_sf, 32610)

# 在投影坐标系下计算最近点
geo_utm <- st_nearest_points(point_utm, line_utm)
new_point_utm <- st_cast(geo_utm, 'POINT')[2] |> st_as_sf(crs = 32610)

# 转回4326坐标系
new_point_4326 <- st_transform(new_point_utm, 4326)

# 验证拓扑接触
st_touches(line_sf, new_point_4326)

方法2：使用st_snap将点吸附到线上 #

如果不想转换坐标系，可直接用st_snap将计算出的最近点吸附到目标线，通过tolerance参数控制精度（4326下单位为度，设为极小值如1e-8即可）：

# 将new_point吸附到line_sf上，容差设为1e-8度
snapped_point <- st_snap(new_point, line_sf, tolerance = 1e-8)

# 验证接触
st_touches(line_sf, snapped_point)

方法3：放宽拓扑判断的精度阈值 #

若只需判断“近似接触”，可使用st_distance检查点线距离是否小于极小阈值（如1e-8度），替代严格的st_touches：

# 检查点线距离是否小于1e-8度
st_distance(new_point, line_sf) < 1e-8

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