是否有一个函数告诉我彼此相邻的单元格

Question

这是我到目前为止的代码用六角形的顶点做点事。

library(dggridR)
library(dplyr)
library(ggplot2)

#Construct a global grid with cells approximately 1000 miles across
dggs          <- dgconstruct(spacing=1000, metric=FALSE, resround='down')

#Load included test data set
data(dgquakes)

#Get the corresponding grid cells for each earthquake epicenter (lat-long pair)
dgquakes$cell <- dgGEO_to_SEQNUM(dggs, dgquakes$lat, dgquakes$lon)$seqnum

#Get the number of earthquakes in each equally-sized cell
quakecounts   <- dgquakes %>% group_by(cell) %>% summarise(count=n())
quakecounts   <- dgquakes %>% group_by(cell) %>% summarise(total=sum())


#Construct a global grid with cells approximately 1000 miles across
dggs <- dgconstruct(spacing=1000, metric=FALSE, resround='down')
nocells <- dgmaxcell(dggs)

grid <- dgrectgrid(dggs,
                   minlat=-90, minlon=-180, 
                   maxlat=90, maxlon=180, frame = TRUE)


label <- 1:nocells
qgrid <- data.frame(label)

# get the lat-long so we can plot it properly
cellcentres <- dgSEQNUM_to_GEO(dggs, 1:nocells)
qgrid$cent_long <- cellcentres$lon_deg
qgrid$cent_lat <- cellcentres$lat_deg


# create plot
p <- ggplot() +
  geom_path(data=grid, aes(x=long, y=lat, group=group), color="black")  + # draw each hexagon
  geom_label  (data=qgrid,aes(x=cent_long, y=cent_lat, label=label)) + # draw label at cell centres
  coord_equal()

p

Answer 1

dgrectgrid（）的输出是一个空间对象，因此您可以在“ Terra”软件包中尝试“相邻”函数，或者在“ SF”软件包中使用“ ST_Touches”函数：

library（sf）

neighbours_sf＆lt; - st_touches（grid）

或

library> library>

neighbours_terra＆lt; - aideacter（vect（vect（grid），grid），（grid），（grid），grid），（grid），grid） type =“触摸”）