Figure 6: GeoPandas / spatial data - choropleth map - population with cities and rivers of Africa

Fig6_PyGMT_geopandas.ipynb

@@ -0,0 +1,35 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "id": "760146f6-d4c9-4e69-b539-f42e057c7945",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Copy finale version of script from normal Python file"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.12.3"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

Fig6_PyGMT_geopandas.py

@@ -0,0 +1,33 @@
+import pygmt
+import geopandas as gpd
+
+world = gpd.read_file("https://naciscdn.org/naturalearth/50m/cultural/ne_50m_admin_0_countries.zip")
+world["POP_EST"] *= 1e-6
+rivers = gpd.read_file("https://naciscdn.org/naturalearth/110m/physical/ne_110m_rivers_lake_centerlines.zip")
+cities = gpd.read_file("https://naciscdn.org/naturalearth/110m/cultural/ne_110m_populated_places_simple.zip")
+cities_world = cities[cities["worldcity"]==1]  # Focus on large cities
+
+fig = pygmt.Figure()
+fig.basemap(region=[-19.5, 53, -38, 37.5], projection="M15c", frame=True)
+pygmt.makecpt(cmap="bilbao", series=(0, 270, 10))
+fig.plot(
+    data=world[["POP_EST", "geometry"]],
+    pen="1p,gray50",
+    fill="+z",
+    cmap=True,
+    aspatial="Z=POP_EST",
+)
+fig.colorbar(frame="x+lPopulation (millions)")
+fig.plot(data=rivers["geometry"], pen="1.5p,darkblue")
+fig.plot(data=cities["geometry"], style="s0.25c", fill="darkorange", pen="1p")
+fig.text(
+    x=cities_world.geometry.x,
+    y=cities_world.geometry.y,
+    text=cities_world["name"],
+    font="10p,Helvetica-Bold",
+    offset="0c/-0.25c",
+    justify="TR",
+    fill="white@30",
+)
+fig.show()
+fig.savefig(fname="Fig6_PyGMT_geopandas.png")