Welcome to Software Development on Codidact!

Will you help us build our independent community of developers helping developers? We're small and trying to grow. We welcome questions about all aspects of software development, from design to code to QA and more. Got questions? Got answers? Got code you'd like someone to review? Please join us.

Post History

60%

+1 −0

Q&A Add additional shapes on a GeoPandas matplotlib

Apparently, it doesn't take much to add another set of features. You create the axes first and refer to them when plotting additional shapes. Adapted from GeoPandas' geocoding example. import matp...

posted 1y ago by Michael‭ · edited 1y ago by Michael‭

Answer

#3: Post edited by

Michael‭ · 2025-07-18T20:25:57Z (about 1 year ago)
Fix save method

Copy Link

Raw

Markdown

Apparently, it doesn't take much to add another set of features. You create the axes first and refer to them when plotting additional shapes. Adapted from [GeoPandas' geocoding example](https://geopandas.org/en/stable/docs/user_guide/geocoding.html).
```py
import matplotlib
import matplotlib.pyplot as plt # Don't forget this import
import geopandas as gpd
import contextily as ctx
points = gpd.read_file('points.geojson')
points = points.to_crs('EPSG:4326')
# Snipped a bunch of filtering/adjustment here.
fig, ax = plt.subplots(figsize=(15, 15))
points.plot(ax=ax, markersize=20)
shapes.plot(ax=ax, linewidth=3)
ctx.add_basemap(ax, crs=points.crs, source=ctx.providers.CartoDB.Positron)
~~with open('render.png', 'wb') as img:~~
img.write(ax.set_axis_off())
```

Apparently, it doesn't take much to add another set of features. You create the axes first and refer to them when plotting additional shapes. Adapted from [GeoPandas' geocoding example](https://geopandas.org/en/stable/docs/user_guide/geocoding.html).
```py
import matplotlib
import matplotlib.pyplot as plt # Don't forget this import
import geopandas as gpd
import contextily as ctx
points = gpd.read_file('points.geojson')
points = points.to_crs('EPSG:4326')
# Snipped a bunch of filtering/adjustment here.
fig, ax = plt.subplots(figsize=(15, 15))
points.plot(ax=ax, markersize=20)
shapes.plot(ax=ax, linewidth=3)
ctx.add_basemap(ax, crs=points.crs, source=ctx.providers.CartoDB.Positron)
fig.savefig('render.png')
```

#2: Post edited by

Alexei‭ · 2025-07-18T06:18:21Z (about 1 year ago)
fixed the link

Copy Link

Raw

Markdown

~~Apparently, it doesn't take much to add another set of features. You create the axes first and refer to them when plotting additional shapes. Adapted from [GeoPandas' geocoding example](geocoding).~~
```py
import matplotlib
import matplotlib.pyplot as plt # Don't forget this import
import geopandas as gpd
import contextily as ctx
points = gpd.read_file('points.geojson')
points = points.to_crs('EPSG:4326')
# Snipped a bunch of filtering/adjustment here.
fig, ax = plt.subplots(figsize=(15, 15))
points.plot(ax=ax, markersize=20)
shapes.plot(ax=ax, linewidth=3)
ctx.add_basemap(ax, crs=points.crs, source=ctx.providers.CartoDB.Positron)
with open('render.png', 'wb') as img:
img.write(ax.set_axis_off())
~~```~~
~~[geocoding]: https://geopandas.org/en/stable/docs/user_guide/geocoding.html~~

Apparently, it doesn't take much to add another set of features. You create the axes first and refer to them when plotting additional shapes. Adapted from [GeoPandas' geocoding example](https://geopandas.org/en/stable/docs/user_guide/geocoding.html).
```py
import matplotlib
import matplotlib.pyplot as plt # Don't forget this import
import geopandas as gpd
import contextily as ctx
points = gpd.read_file('points.geojson')
points = points.to_crs('EPSG:4326')
# Snipped a bunch of filtering/adjustment here.
fig, ax = plt.subplots(figsize=(15, 15))
points.plot(ax=ax, markersize=20)
shapes.plot(ax=ax, linewidth=3)
ctx.add_basemap(ax, crs=points.crs, source=ctx.providers.CartoDB.Positron)
with open('render.png', 'wb') as img:
img.write(ax.set_axis_off())
```

#1: Initial revision by

Michael‭ · 2025-07-17T16:36:35Z (about 1 year ago)

Copy Link

Raw

Markdown

Apparently, it doesn't take much to add another set of features. You create the axes first and refer to them when plotting additional shapes. Adapted from [GeoPandas' geocoding example](geocoding).

```py
import matplotlib
import matplotlib.pyplot as plt  # Don't forget this import
import geopandas as gpd
import contextily as ctx

points = gpd.read_file('points.geojson')
points = points.to_crs('EPSG:4326')

# Snipped a bunch of filtering/adjustment here.

fig, ax = plt.subplots(figsize=(15, 15))

points.plot(ax=ax, markersize=20)
shapes.plot(ax=ax, linewidth=3)

ctx.add_basemap(ax, crs=points.crs, source=ctx.providers.CartoDB.Positron)

with open('render.png', 'wb') as img:
    img.write(ax.set_axis_off())
```

[geocoding]: https://geopandas.org/en/stable/docs/user_guide/geocoding.html

Communities

Post History