查找所有点与多边形边界之间的最小距离

Question

我想找到点与多边形边界之间的最小距离（所有点都位于多边形内）。如果可能的话，我怎样才能提取值？，这样我就可以找到个体数量与边界距离之间的相关性。

多边形采用 .SHP 格式并指向 X/Y 坐标。

任何遗漏的信息请告诉我！非常感谢您的帮助！

Answer 1

spatstat 包有一个函数 nncross 可以找到两组点或一组点和一组线段之间的最近邻居。

加载一组 x/y 值来创建 spatstat 点图案对象相对容易：如果 X 和 Y 是包含坐标的两个向量，则可以使用 [=15 创建点图案对象=]

library(spatstat)
p = ppp(x,y)

您需要将 shp 数据转换为 spatstat 段模式对象。为此，您可以使用 maptools 中的命令加载 shp 文件，然后将其转换为 spatstat 对象：

library(maptools)
shp = readShapeSpatial("yourdata.shp") #read shp file
shp = as.psp(shp) # convert to psp object

要计算最近邻距离，您必须使用 nncross

nncross(p,shp)

Answer 2

单位正方形多边形：

library(sp)
x = cbind(c(0,1,1,0,0),c(0,0,1,1,0))
pol = SpatialPolygons(list(Polygons(list(Polygon(x)), "ID")))

单位正方形内的随机点数：

set.seed(131)
pts = SpatialPoints(cbind(runif(10), runif(10)))
plot(pol)
points(pts, col = 'red')

计算距离：

library(rgeos)
gDistance(pts, pol, byid = TRUE) # will be 0, all inside
gDistance(pts, as(pol, "SpatialLines"), byid = TRUE) # dist to line

加入剧情：

text(coordinates(pts),
  as.character(
    round(as.vector(gDistance(pts, as(pol, "SpatialLines"), byid = TRUE)), 3)),
pos = 4)

使用包 rgdal

中的 readOGR 将多边形数据从 shapefile 读入 R

Answer 3

按照@xraynaud 的步骤（稍作修改）：

library(maptools)
shp = readShapeSpatial("yourdata.shp") #read shp file
W = as.owin(shp) # convert to owin object

library(spatstat)
p = ppp(x, y, window = W)

现在 p 是一个点模式，包含由多边形界定的点。要计算每个点到边界多边形的距离（通常在 spatstat 术语中称为 window）：

d = bdist.points(p)

现在 d 是一个距离向量。

查找所有点与多边形边界之间的最小距离

Finding the minimum distance between all points and the polygon boundary

r

spatial

spatstat