geo32.py

ax = gplt.choropleth(
  cv19,
  hue="deaths_per_mil",
  edgecolor="white",
  linewidth=5,
  cmap="Reds",
  alpha=0.8,
  projection=gcrs.AlbersEqualArea(),
  figsize=(30, 30)
)

ax = gplt.pointplot(
  continental_usa_cities,
  hue="POP_2010",
  cmap="Greens",
  scheme="quantiles",
  scale="POP_2010",
  limits=(3, 50),
  zorder=2,
  alpha=0.4,
  ax=ax
)