2

我正在尝试从包含两个相邻点之间距离的数据框中生成我最能描述为成对距离矩阵的内容。这些不是欧几里得距离，它们本质上是海岸线上点之间的距离，因此距离不是直线。我能够riverdist使用地理空间数据在包中生成一个距离矩阵，但这只是两个点之间的完整距离，我现在正试图做点之间这个距离的一个子集。

我一直在寻找一种方法来做到这一点，并且一直空手而归。任何帮助将非常感激。

这是一个例子：

我有这个数据：

mat <- matrix(
 c(
     3, #distance between 1 and 2
    10, #distance between 2 and 3
     7, #distance between 3 and 4
     9  #distance between 4 and 5
),
nrow=4, ncol=1, dimnames =     list(c("site1","site2","site3","site4"),c("dist")))

> mat
  dist
site1    3
site2   10
site3    7
site4    9

我想产生以下“距离”矩阵：

     site1 site2 site3 site4 site5
site1    0
site2    3     0     
site3   13    10     0   
site4   20    17     7     0
site5   29    26    16     9     0

对于此任务，原始数据可能会更好地组织如下：

     SiteA SiteB Dist
   1 site1 site2    3
   2 site2 site3   10
   3 site3 site4    7
   4 site4 site5    9

有什么建议吗？

4

3 回答 3

2

我认为您可能过于关注实际距离，而您的问题更多是求和问题。

至少据我了解，您不需要寻找最短路线或类似问题，而只需添加数字即可。所以 from ato意味着从你的变量中b添加行。唯一困难的事情是处理 from 和 to 向后或相同的情况。无论如何，我明白了：ab-1mat

dist <- function(a,b) abs(sum(if(a>1 && b>1) mat$dist[(a:b)-1][-1] else mat$dist[(a:b)-1]))
distmat <- sapply(1:5, function(i) {
  sapply(1:5, dist, i)
})
于 2018-12-10T23:36:40.253 回答
1

这会做到（尽管可能有点慢）：

dist = function(mat){
    tmp_mat = c(0, mat)
    dist_mat = diag(0, length(mat)+1)
    for (i in 1:length(mat))
        dist_mat[(i+1):(length(mat)+1),i] = cumsum(tmp_mat[-(1:i)])

    dist_mat = dist_mat + t(dist_mat)
    return (dist_mat)
    }
dist(mat)
于 2018-12-10T23:39:33.627 回答
1

这是一个累积距离，所以取一个累积值，然后进行cumance计算：sumdist

mat <- c(3,10,7,9)
dist(cumsum(c(0,mat)))
#   1  2  3  4
#2  3         
#3 13 10      
#4 20 17  7   
#5 29 26 16  9
于 2018-12-11T02:06:23.083 回答