行与sf之间的点距离

时间:2018-04-16 09:18:04

标签: r dplyr sf

我在sf类型的简单要素(POINT)中保存了多个轨迹。我想计算后续位置(即行)之间的欧几里德距离。到目前为止,我已经使用Pythagorean formula for calculating Euclidean Distances in 2D space“手动”计算了距离。我想知道我是否可以使用函数sf::st_distance()做同样的事情。这是一个简单的例子:

library(sf)
library(dplyr)

set.seed(1)

df <- data.frame(
  gr = c(rep("a",5),rep("b",5)),
  x  = rnorm(10),
  y = rnorm(10)
 )

df <- st_as_sf(df,coords = c("x","y"),remove = F)


df %>%
  group_by(gr) %>%
  mutate(
    dist = sqrt((lead(x)-x)^2+(lead(y)-y)^2)
  )
#> Simple feature collection with 10 features and 4 fields
#> geometry type:  POINT
#> dimension:      XY
#> bbox:           xmin: -0.8356286 ymin: -2.2147 xmax: 1.595281 ymax: 1.511781
#> epsg (SRID):    NA
#> proj4string:    NA
#> # A tibble: 10 x 5
#> # Groups:   gr [2]
#>    gr         x       y   dist                 geometry
#>    <fct>  <dbl>   <dbl>  <dbl>                  <POINT>
#>  1 a     -0.626  1.51    1.38     (-0.6264538 1.511781)
#>  2 a      0.184  0.390   1.44     (0.1836433 0.3898432)
#>  3 a     -0.836 -0.621   2.91   (-0.8356286 -0.6212406)
#>  4 a      1.60  -2.21    3.57        (1.595281 -2.2147)
#>  5 a      0.330  1.12   NA         (0.3295078 1.124931)
#>  6 b     -0.820 -0.0449  1.31  (-0.8204684 -0.04493361)
#>  7 b      0.487 -0.0162  0.992  (0.4874291 -0.01619026)
#>  8 b      0.738  0.944   0.204    (0.7383247 0.9438362)
#>  9 b      0.576  0.821   0.910    (0.5757814 0.8212212)
#> 10 b     -0.305  0.594  NA       (-0.3053884 0.5939013)

我想用dist计算sf::st_distance()。我该怎么做呢?

1 个答案:

答案 0 :(得分:7)

首先要了解的是sf,几何列(类sfc之一)作为列表列存储在数据框内。 通常对列表列执行任何操作的关键是使用purrr::map和朋友,或者使用接受list-cols作为参数的函数。在st_distance的情况下,其参数可以是sf(数据框),sfc(几何列),甚至是sfg(单个geom行)的对象),所以不需要map和朋友。解决方案看起来应该是这样的:

df %>%
  group_by(gr) %>%
  mutate(
    dist = st_distance(geometry)
  )

然而,这不起作用。经过一番调查,我们发现了两个问题。首先,st_distance返回距离矩阵而不是单个值。要解决此问题,我们会使用by_element = T的{​​{1}}参数。

接下来,我们不能只执行st_distance,因为dist = st_distance(geometry, lead(geometry), by_element = T)仅适用于矢量列,而不适用于列列。

要解决第二个问题,我们使用lead自行创建潜在客户列。

以下是完整的解决方案:

geometry[row_number() + 1]