R 用户相当新手，我确信有一个简单的解决方案 - 但我找不到它。我有一个包含一系列空间坐标以及许多其他属性的数据框。许多空间坐标是完全相同的 - 我想向它们添加一定量的噪声，以便我可以将它们保持在一定的半径内 - 在本例中为 0.4 米或 40 厘米，同时还跟踪它们相关属性。
我本质上是在寻找这个问题的R版本：https://gis.stackexchange.com/questions/35479/adding-noise-to-overlapping-xy-coordinates-so-no-longer-in-exact-same-地方
...因为当我按照使用 ArcGIS 回答这个问题的说明时 - 我得到了一系列随机点，但我丢失了与这些点关联的属性，并且我无法轻松找到一种方法来匹配它们。
有没有办法使用 R 中的抖动函数或类似的函数并指定空间上下文（例如 40 厘米）中的半径，空间坐标在该范围内随机分布？我不明白如何操纵因子和数量参数以获得我想要的输出。
编辑：
这是一个带有假坐标的 df 示例。正如您所看到的，第一列和第三列中的坐标是相同的，因为动物在同一块岩石下两次。我希望能够为这些坐标添加抖动，使它们略有不同，但我想将抖动控制在40厘米以内（不超过岩石的大小）
mydf <- data.frame("point_id" = 1:6, "date_time" = c("6/5/2018 10:57","6/5/2018 14:30","6/6/2018 10:06","6/6/2018 11:06","6/7/2018 10:35","6/7/2018 15:50"), "Animal_ID" = c(4,5,4,5,4,6), "Rock_ID" = c(1,2,1,3,4,5), x_proj = c(831120.3759,831441.0415,831120.3759,831433.4414,831128.4778,831422.0822), y_proj = c(5877582.998,5875337.074,5877582.998,5875328.897,5877575.360,5875338.216))
#make a separate object for the coordinates#
xy <- mydf[,c(5,6)]
#Convert to a spatialpoints data frame (insert own epsg)
sp.mydf <- SpatialPointsDataFrame(coords = xy, data …这应该是一个需要解决的简单问题，但我无法获得我想要的确切输出。我有一个嵌套的数据帧列表，我想过滤掉所有少于 50 行的数据帧，并将它们从列表中删除。
这是我尝试过的可重现的示例 -
L <- list(iris,mtcars,iris)
O <- list(iris,mtcars,iris)
H <- list(iris,mtcars,iris)
List <- list(L,O,H)
test <- lapply(List, function(x) lapply(x, function(x) if (nrow(x)<50) NULL else x)))
这适用于第一个列表，但它用 NULL 替换嵌套列表中的 mtcars 数据帧 - 它不会将它们从列表中删除。不幸的是，它不会循环其他列表。我也尝试过使用过滤功能
test <- lapply(List, function(x) lapply(x, function(x) filter(x, nrow(x)>50)))
这与不循环所有列表有同样的问题，对于第一个列表，它给我留下了一个空的 df ，它仍然是列表的一个元素。我的最后一个解决方案是编写一个 for 循环，我在嵌套中的第一个列表上尝试了它，这基本上有效 - 但如果可能的话，我想找到一种不太笨重的方法来做到这一点。这也会返回一个错误：Error in List[[1]][[ii]] : subscript out ofbounds
for (ii in seq_along(List[[1]])){
n_rows = nrow(List[[1]][[ii]])
if (n_rows < 20){
List[[1]][[ii]] = NULL
}
}
我希望有一个简单的解决方案即将出现！