points_lines.py

# Some quick snippets for points to lines

address_df = pd.DataFrame([[0, Point([c for c in h3.cell_to_latlng(h3_8)[::-1]])]], columns=['id', 'geometry'])
address_gdf = gpd.GeoDataFrame(address_df, geometry="geometry", crs="epsg:4326")

gdf = gpd.read_file(
    "./data/all_the_data.fgb",
    engine="pyogrio", bbox=location_box)
gdf = gdf.to_crs('epsg:3310') # https://epsg.io/3310

print(gdf.sindex)

sjoin = gpd.sjoin_nearest(left_df=address_gdf, right_df=gdf, how="left", max_distance=10000)
sjoin = sjoin.dropna()