在R语言数据框中进行GPS坐标随机抖动时生成重复值的问题

问题原因 #

你当前代码里的runif(1)只会生成单个随机值，这个值会被重复应用到数据框的每一行，所以相同原始坐标的行得到的抖动结果完全一致。要实现逐行独立抖动，需要让随机数生成函数为每一行都生成一个独立的随机值。

解决方案 #

方法1：生成对应行数的随机数（简洁高效） #

直接让runif生成和数据框行数相同的随机数，每一行都会拿到独立的随机偏移值，同时简化计算逻辑：

require(tidyverse)

# 生成示例数据
centerlong <- c(-90.28192, -90.28192, -90.28192, -90.28192, -90.28192, 
  -90.31374, -90.31374, -91.51432, -92.12345, -93.12345)

centerlat <- c(44.12345, 44.12345, 44.12345, 44.28567, 44.28567, 
  43.98243, 43.98243, 45.00249, 42.12345, 41.12345)

df <- data.frame(centerlong, centerlat)

# 逐行生成独立抖动值
df <- df %>% 
  mutate(
    Longitude = centerlong + runif(nrow(df), min = -0.0005, max = 0.0005),
    Latitude = centerlat + runif(nrow(df), min = -0.0005, max = 0.0005)
  )

方法2：用rowwise()强制逐行执行 #

在tidyverse中，通过rowwise()让mutate里的函数逐行运行，这样每次runif(1)都会为当前行生成新的随机数：

df <- df %>% 
  rowwise() %>%
  mutate(
    Longitude = runif(1)*((centerlong+0.0005)-(centerlong-0.0005))+(centerlong-0.0005),
    Latitude = runif(1)*((centerlat+0.0005)-(centerlat-0.0005))+(centerlat-0.0005)
  ) %>%
  ungroup() # 记得取消行分组，避免后续操作异常

方法3：使用base R的jitter()函数（最简洁） #

如果只是需要给数值添加随机抖动，base R的jitter()函数可以直接实现，自动为每个元素生成独立的随机偏移：

df$Longitude <- jitter(df$centerlong, amount = 0.0005)
df$Latitude <- jitter(df$centerlat, amount = 0.0005)

amount参数控制抖动的最大范围（正负amount之间），正好匹配你"半个长街区"的需求。

验证结果 #

运行上述任意一种方法后，相同原始坐标的行都会得到独立的抖动值，示例输出如下：

内容的提问来源于stack exchange，提问作者Morrigan