get nearest points and set geometry as the point closest to the original measurement

gps_nearest_clean.R

library(sf)
library(dplyr)
library(mapview)
mapviewOptions(fgb = FALSE)
plts <- read_sf('data/PlotData/KuamutLocations.shp') 
pnts <- read_sf("data/Kuamut_2022_NewGPS_Points_Shapefiles/PlotCentre.shp")

pnts_clean <- pnts %>% 
  filter(is.na(NUMBER)|NUMBER!=100) %>%
  mutate(nearest = st_nearest_feature(., plts),
         plot_join = plts$Pl[nearest], # not really needed but just in case the row numbers and plot numbers don't match
         NUMBER = as.integer(NUMBER),
         NUMBER = case_when(is.na(NUMBER) ~ plot_join,
                            TRUE~NUMBER),
         dist = st_distance(., plts[nearest,], by_element=TRUE)) %>% 
  group_by(NUMBER) %>%
  mutate(geometry = geometry[which.min(dist)])

  
mapview(plts, col.regions='red') +
  mapview(pnts, col.regions='blue') +
  mapview(pnts_clean, col.regions='green')