st_coordinates()
sf
inla()
mapview::mapview( map, alpha.regions = 0.4, lwd = 0, map.types = "OpenStreetMap.Mapnik" )