Adjust panel size according to bbox when plotting multiple samples

[lambdamoses]

In functions such as plotSpatialFeature, data from different samples are plotted with facet_wrap. I can't use scales = "free" with geom_sf. So I got this kind of problem which doesn't look good:

I want to adjust the size of the facets according to the size of the bounding box. I still want all panels to have the same scale and the same scale in x and y, so it's not entirely free. facet_grid with space = "free" does not work with geom_sf either. The culprit is coord_sf. I don't need the CRS and reprojection stuff here. So I suppose I should write my own internal variant of geom_sf that allows for flipping the y axis to be consistent with other packages (see pachterlab/SpatialFeatureExperiment#47) and can allow for something like space = "free" but not scales = "free".