4

SpatialLines让我们假设一个长度类的线段len。这条特定的线从左上角开始。

library(sp)
x <- structure(list(x = c(-7.23437435517476, 6.35937810318614, -5.86718660792582, 
                        7.96094089282062), y = c(7.08139459814975, 6.8633712983227, -7.61337581019376, 
                        -6.2180266913006)), .Names = c("x", "y"))

xline <- SpatialLines(list(Lines(Line(x), ID = 1)))
#len <- LineLength(as.matrix(data.frame(x)))
len <- LineLength(as.matrix(data.frame(coordinates(xline))))

plot(0,0, xlim = c(-10, 10), ylim = c(-10, 10), type = "n")
lines(xline)

在距起点所需距离的线上找到点

我想在这条线上找到一个点，它findme距离线的起点有一个单位。例如，如果我要从起点开始沿着直线寻找一个 10 个单位的点，我会在第一段和第二段之间的节点附近找到一个点。欢迎您对更强大的解决方案提出意见。

我尝试使用spsample（见下文）找到它，但这种方法（太）不可靠，不适用于线的后半部分。

# very approximate method, not very suitable
findme <- 11 # 11, 12 and 13 give same result
segs <- 1/(findme/xline.length)
xsam <- spsample(x = xline, n = segs, type = "regular", offset = 0)
points(xsam)
4

2 回答 2

4

以下步骤将帮助您找到坐标。

线路的一般信息：

library(np)

coord <- coordinates(xline)[[1]][[1]] 
nLines <- nrow(coord) - 1
#lengths <- sapply(seq_len(nLines), function(x) LineLength(coord[c(x, x + 1), ]))
lengths <- LineLength(coord, sum = FALSE)

找到新坐标：

findme <- 11 # the distance of the new coordinates

distances <- cumsum(lengths) - findme         # distances from the nodes
segment <- which(distances >= 0)[1]           # the segment of interest
distToNode <- distances[segment]
ratio <- distToNode / lengths[segment]
segCoord <- coord[c(segment, segment + 1), ]
newCoord <- (1 - ratio) * segCoord[2 , ] + ratio * segCoord[1 , ]

阴谋：

points(newCoord[1], newCoord[2])

在此处输入图像描述

于 2012-11-13T17:56:36.587 回答
0

您可以...吗：

  • 确定段的顺序
  • 确定每段的长度
  • 使用顺序和长度来确定点在哪个段上
  • 从 order * length 中减去段的开头（例如，如果您知道自己在第三段上，则将减去第三段和第二段相交点的长度（这只是段 1 的总和） 2）从整行开始）。因此，如果您试图从线的起点找到 7 个单位的点，并且段 1 和 2 的总和为 5，那么您将得到 2。
  • 然后只需使用此方法处理唯一剩余的段，以确定该段向下 2 个单位的点。

如果这不是完整的解决方案，我很抱歉。但也许它会给你一个可以解决的方法。

希望有帮助...

于 2012-11-13T17:13:33.960 回答