Skip to content

Instantly share code, notes, and snippets.

@tastatham

tastatham/gadm_usa.py

Created February 24, 2022 10:39
Show Gist options
  • Save tastatham/1e7ada20d39edfa75ff8fbabb4aa259e to your computer and use it in GitHub Desktop.
Save tastatham/1e7ada20d39edfa75ff8fbabb4aa259e to your computer and use it in GitHub Desktop.
Download ZIP
dask_geopandas
Raw
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
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment