I am facing a similar problem than the one described here as I need to divide the US East coastline in segments of a certain length. I am trying to follow your approach Matthew, but I am getting issues when executing coastline = cl_gdf.iloc[0].geometry. For some reason, when reprojecting from 4326 to 9311 all values are apparently inf. I am trying to use the digitalized dataset that you mention. Would you know why this may be happening? Am I missing any information? Thanks.