如何利用R语言terra包获取距离多边形边缘点最近的多边形？

遇到点刚好落在多边形边缘导致terra::extract()返回NA的情况确实很头疼，不过我们可以用terra包自带的功能来获取每个点对应的最近多边形，这里有两种实用的方法：

方法1：使用terra::extract()的method="nearest"参数 #

terra::extract()其实内置了处理边缘/外部点的参数，只需要在调用时指定method="nearest"，就能自动返回距离点最近的多边形属性。如果需要限制搜索范围，还可以通过buffer参数设置一个阈值（比如0，确保只考虑相邻的多边形）：

# 使用method="nearest"提取最近多边形属性
e_nearest <- terra::extract(world_shp, pts, method = "nearest", buffer = 0)
# 查看对应的主权国家名称
e_nearest$sovereignt

这种方法最简单直接，不需要额外的复杂计算，就能得到你需要的结果。

方法2：手动计算距离并筛选最近多边形 #

如果需要更灵活的控制（比如自定义距离筛选逻辑），可以用terra::distance()函数计算每个点到所有多边形的距离，再找到距离最小的那个多边形：

# 计算每个点到所有多边形的距离矩阵
dist_matrix <- terra::distance(pts, world_shp)
# 找到每个点距离最小的多边形索引
min_dist_idx <- apply(dist_matrix, 1, which.min)
# 根据索引提取对应的多边形
nearest_polygons <- world_shp[min_dist_idx, ]
# 获取主权国家名称
nearest_polygons$sovereignt

这种方法的优势是你能看到每个点到所有多边形的距离值，方便后续做更精细化的筛选（比如排除距离过远的多边形）。

为什么原始代码返回NA？ #

默认情况下，terra::extract()只会返回完全落在多边形内部的点对应的属性，而边缘上的点会被判定为"不在多边形内"，因此返回NA。通过method="nearest"参数或者手动计算距离，就能绕过这个判定逻辑，找到最近的多边形。

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