#30DayMapChallenge 4 - Hexagons: #NASAsedac #GPWv4 global population density sampled into proper #DGGS in #ISEAH7 (level 4) projection with #DGGRID, #dggrid4py, #rasterio / #rasterstats and @geopandas and visualised with #GeoViews @holoviz_org
https://twitter.com/allixender/status/1324055326111485959
NASA Sedac GPWv4 global population density:
https://sedac.ciesin.columbia.edu/data/collection/gpw-v4
Then I generated the ISEAH7 DGGS at level 4, ca. 24000 hexagons, and saved as GeoPackage. I did this with the help of a little library I am working on:
https://github.com/allixender/dggrid4py/
It uses GeoPandas and DGGRID - a free software program for creating and manipulating Discrete Global Grids. DGGRID version 7.0 was released in September, 2019.
https://www.discreteglobalgrids.org/software/
Then I used Rasterio and Rasterstats to summarize (mean) the population density data into the hexagon features. Subsequently I classified the data with PySAL’s mapclassify due to it’s large range of values.
Eventually I used GeoViews to provide a static orthographic (say globe) basemap view.
The Jupyter notebook (view here) is provided in my GitHub repo.