将与空间点相关的协变量分配给R中更大的空间点集？

【问题标题】：Assigning covariate associated to spatial points to a bigger set of spatial points in R?将与空间点相关的协变量分配给R中更大的空间点集？
【发布时间】：2018-02-13 14:04:56
【问题描述】：

我有两个包含空间点的数据集（.csv 格式）：data1 包含 220 个经纬度空间点，data2 包含 80 个经纬度空间点。对于data2，我有一个协变量表示每个点的遗传起源。两个数据集中的空间点并不完全相同。

我想为 data1 中的空间点分配遗传起源。看来我需要在 data2 中的每个点周围定义一个正方形（或其他），以便能够关联 data1 中每个点的遗传起源。

我正在使用 R，我认为作为 raster 或 sp 的包可能有用。

感谢您的帮助。

最好，玛丽。

【问题讨论】：

标签： r geospatial r-raster

【解决方案1】：

您需要考虑如何分配“遗传起源”。似乎暗示的一种方法是将其分配给最近的邻居。

在提出问题时，您应该始终包含一些示例数据。

library(raster)
d1 <- data.frame(lon=c(1,5,55,31), lat=c(3,7,20,22))
d2 <- data.frame(lon=c(4,2,8,65,5,4), lat=c(50,-90,20,32,10,10), origin=LETTERS[1:6], stringsAsFactors=FALSE)

以下是如何根据最近的已知原点分配原点

# make sure your data are (x,y) or (longitude,latitude), not the reverse
pd <- pointDistance(d1, d2[,1:2], lonlat=TRUE)
nd <- apply(pd, 1, which.min)
d1$origin <- d2$origin[nd]

【讨论】：

非常感谢。这正是我想要的。