find neighboring polygons with st_touches

find neighboring polygons with st_touches

library(tidyverse)
library(sf)

demo(nc, ask = FALSE, verbose = FALSE)

nc <- nc %>% 
  mutate(polygon_id = row_number())


ggplot() +
  geom_sf(data = nc)




# with st_intersects
neighbors <- st_intersects(nc)
neighbors[[5]]



# with st_touches
neighbors <- st_touches(nc)
neighbors[[5]]


# add list to your data
nc <- nc %>% 
  mutate(p_id = row_number()) %>% 
  select(NAME, p_id)

nc <- nc %>% 
  mutate(touches = st_touches(.))

nc %>% 
  filter(p_id == 5) %>% 
  pull(touches)