load in gdf and geotiff and return transformed gdf

gdf_transform.py

import rasterio
import geopandas as gpd

def transform_gdf(gdf, geotiff):
    """Transform the coordinates of a GeoDataFrame.

    Arguments
    ---------
    gdf : gpd.GeoDataFrame
        A geodataframe with an identity affine xform.
    geotiff : str
        Path to a GeoTIFF with the desired transformation for the output gdf.

    Returns
    -------
    A gdf with the polygons transformed to the destination crs.

    """
    with rasterio.open(geotiff) as src:
        gdf.to_crs(src.profile['crs'])
        src.close()

    return gdf