Static Mapping with Matplotlib and Contextily

Static mapping remains a cornerstone of spatial analysis and scientific communication. Unlike browser-driven visualizations that prioritize real-time exploration, static cartography delivers reproducible, high-resolution outputs tailored for academic publications, technical reports, and automated data pipelines. Within the broader discipline of Geospatial Visualization & Web Mapping, this approach emphasizes layout precision, typographic control, and print-ready formatting. By pairing Matplotlib’s robust plotting framework with Contextily’s seamless basemap integration, Python developers can generate publication-quality spatial figures with minimal overhead.

Environment Setup #

A reliable geospatial stack is essential before generating figures. Contextily handles web tile fetching and rendering through requests and Pillow, while Matplotlib manages the coordinate axes and export pipeline. Install the required packages via your terminal:

pip install geopandas matplotlib contextily

Then import the core libraries in your Python environment:

import geopandas as gpd
import matplotlib.pyplot as plt
import contextily as ctx

Core Implementation Workflow #

The following workflow demonstrates a production-ready approach to static map generation. It loads a global dataset, isolates a specific region, aligns coordinate reference systems, overlays a web basemap, and exports a high-DPI image. These steps form the static map composition pipeline shown below.

flowchart LR
    A["Load & filter<br/>GeoDataFrame"] --> B["Align CRS<br/>(.to_crs)"]
    B --> C["Plot thematic layer<br/>(ax)"]
    C --> D["Add basemap tiles<br/>(ctx.add_basemap)"]
    D --> E["Refine layout<br/>(title, axis off)"]
    E --> F["Export high-DPI<br/>(savefig)"]

# 1. Load and filter spatial data
# Natural Earth data is no longer bundled with GeoPandas; download the
# 'naturalearth_lowres' shapefile from naturalearthdata.com and point to it here.
world = gpd.read_file('ne_110m_admin_0_countries.shp')
africa = world[world['CONTINENT'] == 'Africa'].to_crs(epsg=4326)

# 2. Configure the plotting canvas
fig, ax = plt.subplots(figsize=(10, 8))
africa.plot(ax=ax, color='#2c7bb6', edgecolor='white', linewidth=0.5)

# 3. Integrate a web basemap
ctx.add_basemap(ax, crs=africa.crs.to_string(), source=ctx.providers.CartoDB.Positron)

# 4. Refine layout and export
ax.set_axis_off()
ax.set_title('African Countries with CartoDB Positron Basemap', fontsize=14, pad=10)
fig.tight_layout()
fig.savefig('africa_static_map.png', dpi=300, bbox_inches='tight')
plt.show()

Step-by-Step Logic & Best Practices #

Coordinate Reference System Alignment #

Web map tiles are universally distributed in Web Mercator (EPSG:3857). Contextily automatically reprojects the plotting axis to match the selected tile provider, but your underlying geospatial data must already possess a defined CRS. Omitting the crs parameter in ctx.add_basemap() frequently results in misaligned layers or blank canvases. Always verify your GeoDataFrame’s projection using gdf.crs and apply .to_crs() when necessary. For authoritative definitions of spatial reference systems, consult the EPSG Geodetic Parameter Registry.

Tile Provider Selection #

Contextily exposes dozens of basemap options through its providers module. While OpenStreetMap serves as a reliable default, professional publications often benefit from muted, low-contrast backgrounds like CartoDB.Positron or CartoDB.Voyager. These palettes prevent visual competition with primary data layers. When designing thematic visualizations that rely on color gradients or statistical shading, selecting an appropriate basemap ensures that your analytical layers remain legible and visually balanced.

Layout Control and Export Optimization #

Matplotlib’s figure management tools allow precise control over margins, typography, and resolution. Using fig.tight_layout() and bbox_inches='tight' during export eliminates unwanted whitespace and ensures crisp edges. For workflows requiring higher performance across large geographic extents, consider Optimizing tile generation for global web maps to reduce network latency and memory overhead. Additional layout strategies can be found in Matplotlib’s official figure customization documentation.

Extending the Workflow #

Static maps serve as the foundation for more advanced spatial storytelling. Once comfortable with base figure generation, analysts often transition to Interactive Maps with Folium and Leaflet for browser-based exploration, or explore 3D Terrain Visualization when elevation data requires depth representation. For researchers seeking to automate report generation, mastering Creating publication-ready static maps in Python provides a reliable template for batch processing. Additionally, temporal datasets can be transformed into Creating animated time-lapse maps with Matplotlib by iterating through time slices and compiling frames into video or GIF formats.

Conclusion #

The combination of Matplotlib and Contextily offers a streamlined, highly customizable pathway for generating static spatial visualizations. By adhering to proper CRS management, selecting context-appropriate basemaps, and leveraging Matplotlib’s export controls, developers can produce figures that meet rigorous editorial standards. As geospatial workflows continue to evolve, static mapping remains an indispensable tool for clear, reproducible, and publication-ready spatial communication.