Extract coordinates from `sf` spatial object and add to object as columns

sp_extract_coords.R

# asssumes UTMs. Column nmaes would be `lat` `long` if in Albers or WSG84
    
#option 1 - retains the sf object 
sf_object %>%
    mutate(
           easting = sf::st_coordinates(.)[,1],
           northing = sf::st_coordinates(.)[,2]) 

# option 2 - this changes from sf object to a dbl_df - not usually helpful
sf_object %>%
    poisspatial::ps_sfc_to_coords(X = 'easting', Y = 'northing')