From 32cc7a2296246ed649d98fda3656c825db182a58 Mon Sep 17 00:00:00 2001
From: Michael Grund <23025878+michaelgrund@users.noreply.github.com>
Date: Tue, 31 Dec 2024 16:57:18 +0100
Subject: [PATCH] Update the gallery example for plotting lines with
 LineString/MultiLineString geometry (#3711)
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Co-authored-by: Dongdong Tian <seisman.info@gmail.com>
Co-authored-by: actions-bot <58130806+actions-bot@users.noreply.github.com>
Co-authored-by: Yvonne Fröhlich <94163266+yvonnefroehlich@users.noreply.github.com>
---
 examples/gallery/lines/linestrings.py | 46 +++++++++++++++++++++++++++
 examples/gallery/lines/roads.py       | 46 ---------------------------
 2 files changed, 46 insertions(+), 46 deletions(-)
 create mode 100644 examples/gallery/lines/linestrings.py
 delete mode 100644 examples/gallery/lines/roads.py

diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py
new file mode 100644
index 00000000000..d3793645c84
--- /dev/null
+++ b/examples/gallery/lines/linestrings.py
@@ -0,0 +1,46 @@
+"""
+GeoPandas: Plotting lines with LineString or MultiLineString geometry
+=====================================================================
+
+The :meth:`pygmt.Figure.plot` method allows us to plot geographical data such as lines
+with LineString or MultiLineString geometry types stored in a
+:class:`geopandas.GeoDataFrame` object or any object that implements the
+`__geo_interface__ <https://gist.github.com/sgillies/2217756>`__ property.
+
+Use :func:`geopandas.read_file` to load data from any supported OGR format such as a
+shapefile (.shp), GeoJSON (.geojson), geopackage (.gpkg), etc. Then, pass the
+:class:`geopandas.GeoDataFrame` object as an argument to the ``data`` parameter of
+:meth:`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")
+print(gdf.head())
+
+# %%
+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()
diff --git a/examples/gallery/lines/roads.py b/examples/gallery/lines/roads.py
deleted file mode 100644
index c2a5f69980a..00000000000
--- a/examples/gallery/lines/roads.py
+++ /dev/null
@@ -1,46 +0,0 @@
-# ruff: noqa: RUF003
-"""
-Roads
-=====
-
-The :meth:`pygmt.Figure.plot` method allows us to plot geographical data such
-as lines which are stored in a :class:`geopandas.GeoDataFrame` object. Use
-:func:`geopandas.read_file` to load data from any supported OGR format such as
-a shapefile (.shp), GeoJSON (.geojson), geopackage (.gpkg), etc. Then, pass the
-:class:`geopandas.GeoDataFrame` as an argument to the ``data`` parameter of
-:meth:`pygmt.Figure.plot`, and style the geometry using the ``pen`` parameter.
-"""
-
-# %%
-import geopandas as gpd
-import pygmt
-
-# Read shapefile data using geopandas
-gdf = gpd.read_file(
-    "https://www2.census.gov/geo/tiger/TIGER2015/PRISECROADS/tl_2015_15_prisecroads.zip"
-)
-# The dataset contains different road types listed in the RTTYP column,
-# here we select the following ones to plot:
-roads_common = gdf[gdf.RTTYP == "M"]  # Common name roads
-roads_state = gdf[gdf.RTTYP == "S"]  # State recognized roads
-roads_interstate = gdf[gdf.RTTYP == "I"]  # Interstate roads
-
-fig = pygmt.Figure()
-
-# Define target region around Oʻahu (Hawaiʻi)
-region = [-158.3, -157.6, 21.2, 21.75]  # xmin, xmax, ymin, ymax
-
-title = "Main roads of O`ahu (Hawai`i)"  # Approximating the Okina letter ʻ with `
-fig.basemap(region=region, projection="M12c", frame=["af", f"WSne+t{title}"])
-fig.coast(land="gray", water="dodgerblue4", shorelines="1p,black")
-
-# Plot the individual road types with different pen settings and assign labels
-# which are displayed in the legend
-fig.plot(data=roads_common, pen="5p,dodgerblue", label="CommonName")
-fig.plot(data=roads_state, pen="2p,gold", label="StateRecognized")
-fig.plot(data=roads_interstate, pen="2p,red", label="Interstate")
-
-# Add legend
-fig.legend()
-
-fig.show()