如何控制添加到 R 中 x、y 空间点的噪声量？

Question

相当新手的 R 用户，我确信有一个简单的解决方案 - 但我找不到它。我有一个包含一系列空间坐标以及许多其他属性的数据框。许多空间坐标完全相同——我想给它们添加一定量的噪声，这样我就可以将它们保持在一定半径内——在本例中为 0.4 米或 40 厘米，同时还跟踪它们关联的属性。

我基本上是在寻找这个问题的 R 版本：https://gis.stackexchange.com/questions/35479/adding-noise-to-overlapping-x-y-coordinates-so-no-longer-in-exact-same-place

...因为当我按照使用 ArcGIS 回答这个问题的说明时 - 我得到了一系列随机点，但我丢失了与这些点相关的属性，而且我无法轻易找到一种方法来将它们匹配起来。

有没有办法在R或类似的东西中使用抖动函数并在空间上下文（例如40厘米）中指定空间坐标在该范围内随机分布的半径？ 我不明白如何操纵 factor 和 amount 参数以获得我想要的输出。

编辑：

这是一个带有假坐标的 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 = tumbling_test, proj4string = CRS("+init=epsg:xxxxx"))

我希望我新生成的坐标仍然包含来自其他列（即 Animal_ID、日期等）的属性数据，因为我在 ArcGIS 中使用的其他方法会生成一系列新的随机点，但我可以' 将它们匹配回属性。

此外，如果有一种方法可以只将抖动添加到在一块岩石上多次出现的点，那将是很好的。例如这里我只需要将抖动添加到第1行和第3行，因为其他坐标不重复。添加抖动后，我想将结果转换回我将导出到 .csv

的常规数据帧

Answer 1

有，如果您使用的是 sf 对象，它会容易得多。

您还需要深入研究坐标参考系统 (CRS)，以便将点抖动到正确的距离。

如果您从 spatialPoints 数据框开始，请使用 st_as_sf() 到 return 一个 sf 对象。

下面是一个可重现的示例，其中的点在 5 公里左右抖动。抖动有些随机，在此示例中变化范围为 ~2-5.5 公里。

library(sf)
#> Linking to GEOS 3.6.2, GDAL 2.2.3, PROJ 4.9.3
library(tidyverse)

# load example data
nc <- read_sf(system.file('gpkg/nc.gpkg', package = 'sf'))

#make single points from polygons
nc_points <- st_centroid(nc)
#> Warning in st_centroid.sf(nc): st_centroid assumes attributes are constant over
#> geometries of x
#> Warning in st_centroid.sfc(st_geometry(x), of_largest_polygon =
#> of_largest_polygon): st_centroid does not give correct centroids for longitude/
#> latitude data

# Transform to a crs that uses meters as the distance
nc_points <- st_transform(nc_points, 3358)

nc_points_jittered <- st_jitter(nc_points, amount = 5000)

p1 <- ggplot(nc_points) + 
  geom_sf() + 
  ggtitle('Original')

p2 <- ggplot(nc_points_jittered) + 
  geom_sf() + 
  ggtitle('Jittered')

p3 <- ggplot() + 
  geom_sf(data = nc_points, color = 'red') + 
  geom_sf(data = nc_points_jittered, color = 'black') + 
  ggtitle('Both')

cowplot::plot_grid(p1, p2, p3, ncol = 1)

^{由 reprex package (v0.3.0)}

于 2020 年 1 月 12 日创建

如何控制添加到 R 中 x、y 空间点的噪声量？

How to control for the amount of noise added to x, y spatial points in R?

random

r

spatial

coordinates

jitter