R语言中实现经纬度点匹配至海岸线最近点的问题求助

你的核心问题在于用distGeo直接计算点到所有海岸线顶点的距离，既低效又无法得到线段上的最近点（只能拿到顶点），而且后续sf对象的创建也有语法错误。下面是用sf包原生函数的正确实现，高效且能得到准确的最近点：

正确代码实现 #

library(sf)
library(rnaturalearth)

# 加载海岸线数据和点数据（和你的示例一致）
coastline <- ne_coastline(scale = "small", returnclass = "sf")
points <- data.frame(
  longitude = c(-71.5374, -72.1234, -70.9876),
  latitude = c(-33.865, -34.567, -33.456)
)
points_sf <- st_as_sf(points, coords = c("longitude", "latitude"), crs = 4326)

# 计算每个点到海岸线的最近点对
# st_nearest_points会返回LINESTRING，连接原点点和最近海岸点
nearest_pairs <- st_nearest_points(points_sf, coastline, pairwise = FALSE)

# 提取每个点对应的最近海岸点（取每个LINESTRING的第二个端点）
closest_coast_points <- st_cast(nearest_pairs, "POINT")[seq(2, length(nearest_pairs)*2, by = 2)]

# 转换为sf对象，并和原数据绑定（保留原数据的所有字段）
closest_coast_sf <- st_sf(
  points,  # 原数据的属性
  geometry = closest_coast_points,
  crs = st_crs(coastline)
)

# 验证：查看结果
print(closest_coast_sf)

# 绘图示例：原点点（红色）、海岸线（灰色）、最近海岸点（蓝色）、连接线段（黑色）
plot(st_geometry(coastline), col = "gray", main = "点到海岸线的最近点")
plot(st_geometry(points_sf), col = "red", pch = 16, add = TRUE)
plot(st_geometry(closest_coast_sf), col = "blue", pch = 16, add = TRUE)
plot(st_geometry(nearest_pairs), col = "black", lwd = 0.5, add = TRUE)

为什么你的原代码会出错？ #

1. distGeo的局限：st_coordinates(coastline)会把所有海岸线的顶点拆成一个巨大的坐标矩阵（比如small尺度的海岸线有十几万顶点），distGeo计算每个点到所有顶点的距离会产生维度爆炸，既慢又没必要——你需要的是海岸线上的最近点（可能在两个顶点之间的线段上），而不是最近的顶点。
2. sf对象创建错误：st_point(closest_points)只能创建单个点，而closest_points是多行矩阵，应该用st_sfc结合lapply逐个创建点，或者直接用st_cast处理st_nearest_points的结果。

补充：用st_nearest_feature的正确方式 #

如果你想先用st_nearest_feature找到最近的海岸线要素（减少计算范围），再计算该要素上的最近点，可以这么做：

# 找到每个点最近的海岸线要素索引
nearest_coast_feature_idx <- st_nearest_feature(points_sf, coastline)

# 提取对应的海岸线要素
nearest_coast_features <- coastline[nearest_coast_feature_idx, ]

# 计算每个点到对应海岸线要素的最近点
nearest_pairs <- st_nearest_points(points_sf, nearest_coast_features, pairwise = TRUE)
closest_coast_points <- st_cast(nearest_pairs, "POINT")[seq(2, length(nearest_pairs)*2, by = 2)]

这种方式在处理大尺度海岸线时能稍微提高效率，因为只针对每个点的最近海岸线要素计算，而不是整个海岸线。

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