GeoPandas: Plotting lines with LineString or MultiLineString geometry

The pygmt.Figure.plot method allows us to plot geographical data such as lines with LineString or MultiLineString geometry types stored in a geopandas.GeoDataFrame object or any object that implements the __geo_interface__ property.

Use geopandas.read_file to load data from any supported OGR format such as a shapefile (.shp), GeoJSON (.geojson), geopackage (.gpkg), etc. Then, pass the geopandas.GeoDataFrame object as an argument to the data parameter of pygmt.Figure.plot, and style the lines using the pen parameter.

import geodatasets
import geopandas as gpd
import pygmt

# Read a sample dataset provided by the geodatasets package.
# The dataset contains large rivers in Europe, stored as LineString/MultiLineString
# geometry types.
gdf = gpd.read_file(geodatasets.get_path("eea large_rivers"))

# Convert object to EPSG 4326 coordinate system
gdf = gdf.to_crs("EPSG:4326")
gdf.head()

	NAME	Shape_Leng	geometry
0	Danube	2.770357e+06	MULTILINESTRING ((8.1846 48.0807, 8.19049 48.0...
1	Douro	8.162452e+05	MULTILINESTRING ((-8.67141 41.14934, -8.64362 ...
2	Ebro	8.269909e+05	MULTILINESTRING ((-4.05971 42.97715, -4.06841 ...
3	Elbe	1.087288e+06	MULTILINESTRING ((8.69715 53.90109, 8.72716 53...
4	Guadalquivir	5.997583e+05	MULTILINESTRING ((-6.37899 36.80363, -6.34806 ...

fig = pygmt.Figure()

fig.coast(
    projection="M10c",
    region=[-10, 30, 35, 57],
    resolution="l",
    land="gray95",
    shorelines="1/0.1p,gray50",
    borders="1/0.1,gray30",
    frame=True,
)

# Add rivers to map
fig.plot(data=gdf, pen="1p,steelblue")

fig.show()