ENH: preserve spatial partitioning information in more methods about dask-geopandas HOT 2 OPEN

One to add - we lose spatial partitions after to_crs. I am just not sure if we can just reproject the partitions GeoSeries or we risk that some points will fall outside in some corner cases.

from dask-geopandas.

I think that in general, that is not guaranteed to work. Basically if reprojected lines are no longer straight lines (which happens for many conversions), a reprojected bounding box will not necessarily contain all points anymore. Quick example:

import shapely
box = shapely.box(0, 40, 20, 60)
poly = shapely.segmentize(poly, 5)

reprojected = geopandas.GeoSeries([box, poly], crs="EPSG:4326").to_crs("EPSG:3035")

from shapely.plotting import plot_polygon
plot_polygon(reprojected[0])
plot_polygon(reprojected[1], color="C1")

from dask-geopandas.