dask_geopandas

gadm_usa.py

import geopandas
import dask_geopandas

path = "https://geodata.ucdavis.edu/gadm/gadm4.0/gpkg/gadm40_USA.gpkg"

gdf = geopandas.read_file(path, layer="ADM_1")
gdf.head()

# Filter and rename cols
gdf = gdf[["NAME_1", "geometry"]].rename(columns={"NAME_1": "State"})
gdf.info()

gdf.plot()
# We can quickly filter out non-US mainland areas by the `cx` method, which filters by a set of bounds
us_main = gdf.cx[-150:-50, 20:50]
us_main.plot(facecolor="none", linewidth=0.5, edgecolor="red")

d_gdf = dask_geopandas.from_geopandas(us_main, npartitions=4)
d_gdf