点输出的多边形和线输入（R中）

Question

我想创建一个.shp文件或坐标数据表，其中输出点要素是线端点在多边形边界上的某个点处触摸的位置以及线与多边形边界交叉的位置。在输出中没有生成点，其中直线位于多边形边界上。

这可以在R？

中完成

Like the intersect tool in Esri？

我有两个班级（在R中）......

1. SpatialLines
2. SpatialPolygons

我似乎无法找到执行此操作的工具/包...

（我是初学者）

Answer 1

您可以使用spatstat包执行此操作，但需要转换 从SpatialLines和SpatialPolygons到psp（平面细分 模式）和owin（观察窗口 - 平面中的多边形）使用 maptools包。

由于您没有提供我将首先构建的示例数据 SpatialPolygons和SpatialLines对象看起来像你的图形：

library(maptools)
#> Loading required package: sp
#> Checking rgeos availability: TRUE
library(spatstat)
#> Loading required package: nlme
#> Loading required package: rpart
#> 
#> spatstat 1.48-0.029       (nickname: 'Alternative Facts') 
#> For an introduction to spatstat, type 'beginner'
A <- owin(c(1,3), c(2,5))
B <- owin(c(4,7), c(1,5))
p <- union.owin(A,B)
p <- as(p, "SpatialPolygons")
l <- psp(c(0.5, 0.5), c(2, 4), c(4, 6), c(2, 4), window = owin(c(0,7), c(0,5)))
l <- as(l, "SpatialLines")
plot(p)
plot(l, col = "green", add = TRUE, lwd = 3)

现在要做的就是你要求你将多边形的边转换为psp 并找到线穿过这些边的位置（作为平面点返回 模式ppp）：

l <- as(l, "psp")
p <- as(p, "owin")
e <- edges(p)
x <- crossing.psp(l, e)
x
#> Planar point pattern: 4 points
#> window: rectangle = [1, 6] x [2, 4] units

将点添加到图中会显示我们发现的内容：

plot(p, main = "")
plot(l, col = "green", add = TRUE, lwd = 3)
plot(x, add = TRUE, col = "red", pch = 20, cex = 3)