生成 GPS 坐标/创建相距 1 公里的均匀分布点的栅格

Question

我想覆盖一个由 bbox lat long 坐标与相距 1 公里的 gps 点 raster 定义的区域。目前我通过以下方式为 bboxbbox=8.9771580802,47.2703623267,13.8350427083,50.5644529365 生成2000 点数：

as.data.frame(cbind(runif(2000,8.9771580802 ,13.8350427083),runif(2000,47.2703623267,50.5644529365)))

由于runif 是一个正态分布，我认为我只需要增加点数即可按照我需要的方式覆盖整个区域。 有没有更聪明的方法呢？ 我需要多少分？

更新

我想我也许可以使用包sp 来完成这项工作，但我仍然不太熟悉设置：

  longitudes <- c(8.9771580802, 13.8350427083)
  latitudes <- c(47.2703623267, 50.5644529365)
  bounding_box <- matrix(c(longitudes, latitudes), nrow = 2, byrow = TRUE, dimnames = list(NULL, c("min", "max")))
  projection <- "+proj=longlat" 
  sp_box<-Spatial(bbox = bounding_box, proj4string = CRS(projection))
  p_sample<-spsample(sp_box, 10000, type="regular")

如果我理解正确，这会给我一些在我的坐标中均匀分布的点。 spsample 有一个单元格大小选项，但我还没有掌握。 BR 安德烈亚斯

Answer 1

我对空间数据和分析不太感兴趣，但第一步可能会有所帮助（我采用不同的坐标来获得一个可重复的示例，它适合德国的维度，我对这些维度有一些感觉）。我确信有一种更优雅的方式，但它应该可以满足您的需求。 geosphere::destPoint() 用于计算给定距离和方向的点，geosphere::distGeo() 计算给定框的南北/东西距离，以计算每个方向需要多少点。然后使用expand_grid() 计算计算出的边界点的每个组合。

还要注意，我将点之间的距离更改为 10,000 米或 10 公里，以获得更少的点和更好的情节。您必须相应地更改数字

nw <- c(5.8 55) 
se <- c(15.1, 47)

lon1 <- nw[1]
lat1 <- nw[2]
lon2 <- se[1]
lat2 <- se[2]


#(1) compute the border points in y direction, going south from the nw-point 
# while keeping lon constant

lat <- geosphere::destPoint(nw, 180, 1:floor(geosphere::distGeo(c(lon1,lat1), 
                                                                 c(lon1,lat2))/10000)*10000) 
lat <- as_tibble(lat) 

#(2) compute the border point in x direction (analog to above)
lon <- geosphere::destPoint(nw, 90, 1:floor(geosphere::distGeo(c(lon1,lat1), 
                                                                c(lon2,lat1))/10000)*10000)
lon <- as_tibble(lon)

# use expand_grid() to compute all combinations 
grid <- tidyr::expand_grid(lat$lat, lon$lon) 
names(grid) <- c("lat", "lon") #nicer names

### for visualizing what we've done, map germany with a grid overlay
germany  <- rnaturalearth::ne_countries(type =  "countries", 
            country = "germany", returnclass = "sf")

ggplot2::ggplot(data = germany)+
  ggplot2::geom_sf()+
  ggplot2::geom_point(data = grid, mapping = aes(x = lon, y = lat), size = 0.01)