Function to extract bounding box of object Useful for clipping basemaps to study region

expand_bbox(shapefile, expand = 0.1, equal = T)