Set default crs for `get_grid()` to be `crs = sf::st_crs(boundary)`
Set default crs for
get_grid()to becrs = sf::st_crs(boundary)