如何在sf R包中按组测量空间点之间的距离答案

【问题标题】：How to measure distance between spatial points by groups in sf R package如何在sf R包中按组测量空间点之间的距离
【发布时间】：2019-08-19 11:32:21
【问题描述】：

我有一个数据集，其中包含 13 年期间不同国家的动物记录、x y 坐标（以米为单位）和发现动物的年份。我想知道如何从每个点测量三个最近邻居之间的距离（以米为单位），但仅在同一年的调查点内测量此距离。我想找出哪些位置对彼此之间的距离小于 1000m，并获取这些位置的 csv。

我正在使用包 sf，我能够估计点之间的距离，但不能按年份进行子设置。我发现包 'spatstat' 可以用 'nndist' 做到这一点，但是我发现设置请求的窗口和理解我需要使用的点模式对象很复杂。我是 R 的新手，使用不同类型的对象和不同的包对我来说非常复杂，所以我想知道是否有一种简单的方法可以在 sf.我愿意接受有关如何最有效地执行此操作的建议，如果您建议使用其他包或其他方式，请帮助我提供将我的数据转换为该包所需的对象类型的代码。

谢谢！

trial <-read.table(text =
"Country    station_code    lat_laea    lon_laea    year
Belize  BF09-1  -2955950    1247610 2009
Belize  BF09-10 -2953600    1248590 2009
Belize  BF09-11 -2954620    1247900 2009
Belize  BF11-13 -2958360    1244020 2011
Belize  BF11-18 -2963740    1240290 2011
Belize  BF11-19 -2963380    1242020 2011
Costa   BraulioCarrilloNP-C16   -3640760    1821170 2011
Costa   BraulioCarrilloNP-C17   -3640730    1823240 2011
Costa   BraulioCarrilloNP-C18   -3642140    1817560 2011
Guatemala   40063   -3178260    1249780 2009
Guatemala   40596   -3183800    1246940 2009
Guatemala   43279   -3182640    1251560 2009", 
                   header = TRUE)

trial.sp <- st_as_sf(trial, coords = c("lat_laea", "lon_laea"), crs = 3035)
plot(st_geometry(trial.sp))

test <- st_distance(trial.sp, trial.sp, by_element = FALSE, which = "Euclidean")
class(test)

test1 <- as.data.frame(test)
write_csv(test1, "test.csv")

我想要一个 csv 文件，其中包含同一年学习的所有位置对，它们之间的距离小于 1000m。

【问题讨论】：

标签： r distance points sf

【解决方案1】：

sf-package 很容易集成到dplyr-framework 中。这将围绕trial.sp 中的每个点创建一个半径为 1000 m 的圆，并通过st_within-函数将其连接到原始数据集。

trial.sp %>%
  st_join(trial.sp %>% st_buffer(dist = 1000), join = st_within) %>%
  filter(year.x == year.y)

这包括在 1000 m 距离内没有其他站点的所有站点。它们只是连接到自身。如果你想删除那些你可以做的：filter(year.x == year.y & station_code.x != station_code.y)

【讨论】：