From 3ece2f8e4ff89916b6e6c0aba05a5f90d4c794da Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 09:55:11 +0100 Subject: [PATCH 01/17] Use Natural Earth data --- examples/gallery/lines/linestrings.py | 31 +++++++++------------------ 1 file changed, 10 insertions(+), 21 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 18f94502f16..1e1fa43728e 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -14,33 +14,22 @@ """ # %% -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")) +# Read a sample dataset provided by Natural Earth. The dataset contains large rivers +# in Europe, stored as LineString/MultiLineString geometry types. +provider = "https://naciscdn.org/naturalearth/" +# rivers = gpd.read_file(f"{provider}110m/physical/ne_110m_rivers_lake_centerlines.zip") +rivers = gpd.read_file(f"{provider}50m/physical/ne_50m_rivers_lake_centerlines.zip") +rivers = rivers[rivers["scalerank"] != 5] +rivers_five = rivers[rivers["scalerank"] == 5] -# Convert object to EPSG 4326 coordinate system -gdf = gdf.to_crs("EPSG:4326") -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, -) +fig.basemap(region=[-12.5, 30, 36, 56.5], projection="M15c", frame=True) +fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.1p,black") # Add rivers to map -fig.plot(data=gdf, pen="1p,steelblue") +fig.plot(data=rivers["geometry"], pen="1p,steelblue") fig.show() From 995a07cfa42d53b7db29ccfab3042d75d5c5df06 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 10:58:33 +0100 Subject: [PATCH 02/17] Use same region as in old example | Remove unneeded code parts --- examples/gallery/lines/linestrings.py | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 1e1fa43728e..f18a1736147 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -20,13 +20,11 @@ # Read a sample dataset provided by Natural Earth. The dataset contains large rivers # in Europe, stored as LineString/MultiLineString geometry types. provider = "https://naciscdn.org/naturalearth/" -# rivers = gpd.read_file(f"{provider}110m/physical/ne_110m_rivers_lake_centerlines.zip") rivers = gpd.read_file(f"{provider}50m/physical/ne_50m_rivers_lake_centerlines.zip") rivers = rivers[rivers["scalerank"] != 5] -rivers_five = rivers[rivers["scalerank"] == 5] fig = pygmt.Figure() -fig.basemap(region=[-12.5, 30, 36, 56.5], projection="M15c", frame=True) +fig.basemap(region=[-10, 30, 35, 57], projection="M15c", frame=True) fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.1p,black") # Add rivers to map From 7a019c9ffe21980ed586e0dd5fe751ec4a0e135c Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 11:04:18 +0100 Subject: [PATCH 03/17] Change to 10 m to get ride of strange error --- examples/gallery/lines/linestrings.py | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index f18a1736147..0fcebc3d798 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -20,8 +20,7 @@ # Read a sample dataset provided by Natural Earth. The dataset contains large rivers # in Europe, stored as LineString/MultiLineString geometry types. provider = "https://naciscdn.org/naturalearth/" -rivers = gpd.read_file(f"{provider}50m/physical/ne_50m_rivers_lake_centerlines.zip") -rivers = rivers[rivers["scalerank"] != 5] +rivers = gpd.read_file(f"{provider}10m/physical/ne_10m_rivers_lake_centerlines.zip") fig = pygmt.Figure() fig.basemap(region=[-10, 30, 35, 57], projection="M15c", frame=True) From 2764faff0d6f50e2dd143aee8ddcffcd38475e15 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 11:08:30 +0100 Subject: [PATCH 04/17] Use map size of old example --- examples/gallery/lines/linestrings.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 0fcebc3d798..31533ec13fc 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -23,7 +23,7 @@ rivers = gpd.read_file(f"{provider}10m/physical/ne_10m_rivers_lake_centerlines.zip") fig = pygmt.Figure() -fig.basemap(region=[-10, 30, 35, 57], projection="M15c", frame=True) +fig.basemap(region=[-10, 30, 35, 57], projection="M10c", frame=True) fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.1p,black") # Add rivers to map From 525e9fdc26097f35a5b40bcc3338a64d69419bb6 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 11:22:39 +0100 Subject: [PATCH 05/17] Adjust URL --- examples/gallery/lines/linestrings.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 31533ec13fc..a80863ef634 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -19,8 +19,8 @@ # Read a sample dataset provided by Natural Earth. The dataset contains large rivers # in Europe, stored as LineString/MultiLineString geometry types. -provider = "https://naciscdn.org/naturalearth/" -rivers = gpd.read_file(f"{provider}10m/physical/ne_10m_rivers_lake_centerlines.zip") +provider = "https://naciscdn.org/naturalearth" +rivers = gpd.read_file(f"{provider}/10m/physical/ne_10m_rivers_lake_centerlines.zip") fig = pygmt.Figure() fig.basemap(region=[-10, 30, 35, 57], projection="M10c", frame=True) From f87bcecc64e7f723ed8174bbbbb8384de5ca2ac8 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 11:41:40 +0100 Subject: [PATCH 06/17] Do not select 'geometry' column expecitly --- examples/gallery/lines/linestrings.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index a80863ef634..b27935eeed5 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -27,6 +27,6 @@ fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.1p,black") # Add rivers to map -fig.plot(data=rivers["geometry"], pen="1p,steelblue") +fig.plot(data=rivers, pen="1p,steelblue") fig.show() From 15500c43928b18ee90212272815d01980e8fe606 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 21:12:43 +0100 Subject: [PATCH 07/17] Adjust comments | Add code for South America --- examples/gallery/lines/linestrings.py | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index b27935eeed5..4bb9920343f 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -17,11 +17,19 @@ import geopandas as gpd import pygmt -# Read a sample dataset provided by Natural Earth. The dataset contains large rivers -# in Europe, stored as LineString/MultiLineString geometry types. +# Read a sample dataset provided by Natural Earth. The dataset contains rivers stored +# as LineString/MultiLineString geometry types. Here will focus on South America. provider = "https://naciscdn.org/naturalearth" rivers = gpd.read_file(f"{provider}/10m/physical/ne_10m_rivers_lake_centerlines.zip") +fig = pygmt.Figure() +fig.basemap(region=[-84.5, -33, -56.5, 13], projection="M10c", frame=True) +fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") + +fig.plot(data=rivers["geometry"], pen="1p,steelblue") +fig.show() + +# %% fig = pygmt.Figure() fig.basemap(region=[-10, 30, 35, 57], projection="M10c", frame=True) fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.1p,black") From 31bb04569593236e99d0210aedf444bba185e3e4 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 21:57:05 +0100 Subject: [PATCH 08/17] Improve intro text | Remove geometry column --- examples/gallery/lines/linestrings.py | 14 ++++++-------- 1 file changed, 6 insertions(+), 8 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 4bb9920343f..33e77adbaf5 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -4,13 +4,11 @@ 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__ `__ 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. +: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` object as an argument to +the ``data`` parameter of :meth:`pygmt.Figure.plot`, and style the lines using the +``pen`` parameter. """ # %% @@ -26,7 +24,7 @@ fig.basemap(region=[-84.5, -33, -56.5, 13], projection="M10c", frame=True) fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") -fig.plot(data=rivers["geometry"], pen="1p,steelblue") +fig.plot(data=rivers, pen="1p,steelblue") fig.show() # %% From a995172c56f17392c77a37ae59a701f66cac90ec Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Thu, 27 Nov 2025 22:16:04 +0100 Subject: [PATCH 09/17] Add blank line and comment --- examples/gallery/lines/linestrings.py | 2 ++ 1 file changed, 2 insertions(+) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 33e77adbaf5..420554c4285 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -24,7 +24,9 @@ fig.basemap(region=[-84.5, -33, -56.5, 13], projection="M10c", frame=True) fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") +# Add rivers to map fig.plot(data=rivers, pen="1p,steelblue") + fig.show() # %% From b4911f6809c76cf23b4437e40e883a192e774768 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Fri, 28 Nov 2025 08:43:01 +0100 Subject: [PATCH 10/17] Use 50 m instead of 10 m | Focus on South America --- examples/gallery/lines/linestrings.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 420554c4285..847f415d68f 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -18,14 +18,15 @@ # Read a sample dataset provided by Natural Earth. The dataset contains rivers stored # as LineString/MultiLineString geometry types. Here will focus on South America. provider = "https://naciscdn.org/naturalearth" -rivers = gpd.read_file(f"{provider}/10m/physical/ne_10m_rivers_lake_centerlines.zip") +rivers = gpd.read_file(f"{provider}/50m/physical/ne_50m_rivers_lake_centerlines.zip") +rivers_sa = rivers.cx[-84.5:-33, -56.5:13] fig = pygmt.Figure() fig.basemap(region=[-84.5, -33, -56.5, 13], projection="M10c", frame=True) fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") # Add rivers to map -fig.plot(data=rivers, pen="1p,steelblue") +fig.plot(data=rivers_sa, pen="1p,steelblue") fig.show() From c582e78dbc46dde640eaae0b27d578f7f6430cba Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Fri, 28 Nov 2025 08:49:24 +0100 Subject: [PATCH 11/17] Subset for Europe --- examples/gallery/lines/linestrings.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 847f415d68f..fc5aa4aef2d 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -31,11 +31,13 @@ fig.show() # %% +rivers_eu = rivers.cx[-10:30, 35:57] + fig = pygmt.Figure() fig.basemap(region=[-10, 30, 35, 57], projection="M10c", frame=True) fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.1p,black") # Add rivers to map -fig.plot(data=rivers, pen="1p,steelblue") +fig.plot(data=rivers_eu, pen="1p,steelblue") fig.show() From 39918b64943c60aa7ba4e06dc93441f0f0159b05 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Fri, 28 Nov 2025 09:56:32 +0100 Subject: [PATCH 12/17] Fix gpd warning --- examples/gallery/lines/linestrings.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index fc5aa4aef2d..6aae0c9b4f8 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -19,7 +19,7 @@ # as LineString/MultiLineString geometry types. Here will focus on South America. provider = "https://naciscdn.org/naturalearth" rivers = gpd.read_file(f"{provider}/50m/physical/ne_50m_rivers_lake_centerlines.zip") -rivers_sa = rivers.cx[-84.5:-33, -56.5:13] +rivers_sa = rivers.cx[-84.5:-33, -56.5:13].copy() fig = pygmt.Figure() fig.basemap(region=[-84.5, -33, -56.5, 13], projection="M10c", frame=True) @@ -31,7 +31,7 @@ fig.show() # %% -rivers_eu = rivers.cx[-10:30, 35:57] +rivers_eu = rivers.cx[-10:30, 35:57].copy() fig = pygmt.Figure() fig.basemap(region=[-10, 30, 35, 57], projection="M10c", frame=True) From 8b7b7300369cc406358a3ab980b88476c8889139 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Fri, 28 Nov 2025 10:49:47 +0100 Subject: [PATCH 13/17] Remove codes for Europe --- examples/gallery/lines/linestrings.py | 12 ------------ 1 file changed, 12 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 6aae0c9b4f8..3d8e633f1e2 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -29,15 +29,3 @@ fig.plot(data=rivers_sa, pen="1p,steelblue") fig.show() - -# %% -rivers_eu = rivers.cx[-10:30, 35:57].copy() - -fig = pygmt.Figure() -fig.basemap(region=[-10, 30, 35, 57], projection="M10c", frame=True) -fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.1p,black") - -# Add rivers to map -fig.plot(data=rivers_eu, pen="1p,steelblue") - -fig.show() From 0c60cabef0f9e5b7e4469bc72e9b345cdc74a9fc Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Fri, 28 Nov 2025 17:20:29 +0100 Subject: [PATCH 14/17] Try region Australia --- examples/gallery/lines/linestrings.py | 12 ++++++++++++ 1 file changed, 12 insertions(+) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index 3d8e633f1e2..f72278d4765 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -19,6 +19,18 @@ # as LineString/MultiLineString geometry types. Here will focus on South America. provider = "https://naciscdn.org/naturalearth" rivers = gpd.read_file(f"{provider}/50m/physical/ne_50m_rivers_lake_centerlines.zip") +rivers_australia = rivers.cx[111:155, -40:-9].copy() + +fig = pygmt.Figure() +fig.basemap(region=[111, 155, -40, -9], projection="M10c", frame=True) +fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") + +# Add rivers to map +fig.plot(data=rivers_australia, pen="1p,steelblue") + +fig.show() + +# %% rivers_sa = rivers.cx[-84.5:-33, -56.5:13].copy() fig = pygmt.Figure() From 7b034720b19dcac76586c46fe33283cd175b3e7d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Fri, 28 Nov 2025 17:35:27 +0100 Subject: [PATCH 15/17] Try Asia --- examples/gallery/lines/linestrings.py | 16 +++++++++++++++- 1 file changed, 15 insertions(+), 1 deletion(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index f72278d4765..af1c07b80be 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -16,9 +16,22 @@ import pygmt # Read a sample dataset provided by Natural Earth. The dataset contains rivers stored -# as LineString/MultiLineString geometry types. Here will focus on South America. +# as LineString/MultiLineString geometry types. Here will focus on Asia. provider = "https://naciscdn.org/naturalearth" rivers = gpd.read_file(f"{provider}/50m/physical/ne_50m_rivers_lake_centerlines.zip") +rivers_asia = rivers.cx[57:125, 7:47].copy() + +fig = pygmt.Figure() +fig.basemap(region=[57, 125, 7, 47], projection="M10c", frame=True) +fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") + +# Add rivers to map +fig.plot(data=rivers_asia, pen="1p,steelblue") + +fig.show() + + +# %% rivers_australia = rivers.cx[111:155, -40:-9].copy() fig = pygmt.Figure() @@ -30,6 +43,7 @@ fig.show() + # %% rivers_sa = rivers.cx[-84.5:-33, -56.5:13].copy() From e5b690a3305eb19466140b6e2c9df3806de62fb8 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Sat, 29 Nov 2025 09:44:55 +0100 Subject: [PATCH 16/17] Remove borders --- examples/gallery/lines/linestrings.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index af1c07b80be..ff1588cd731 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -23,7 +23,7 @@ fig = pygmt.Figure() fig.basemap(region=[57, 125, 7, 47], projection="M10c", frame=True) -fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") +fig.coast(land="gray95", shorelines="1/0.3p,gray50") # Add rivers to map fig.plot(data=rivers_asia, pen="1p,steelblue") @@ -36,7 +36,7 @@ fig = pygmt.Figure() fig.basemap(region=[111, 155, -40, -9], projection="M10c", frame=True) -fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") +fig.coast(land="gray95", shorelines="1/0.3p,gray50") # Add rivers to map fig.plot(data=rivers_australia, pen="1p,steelblue") @@ -49,7 +49,7 @@ fig = pygmt.Figure() fig.basemap(region=[-84.5, -33, -56.5, 13], projection="M10c", frame=True) -fig.coast(land="gray95", shorelines="1/0.3p,gray50", borders="1/0.2p,black") +fig.coast(land="gray95", shorelines="1/0.3p,gray50") # Add rivers to map fig.plot(data=rivers_sa, pen="1p,steelblue") From 4b899fa4bcd6abc23e2ce92644101b89bf738214 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Yvonne=20Fr=C3=B6hlich?= Date: Sat, 29 Nov 2025 10:26:19 +0100 Subject: [PATCH 17/17] Focus on Asia --- examples/gallery/lines/linestrings.py | 26 -------------------------- 1 file changed, 26 deletions(-) diff --git a/examples/gallery/lines/linestrings.py b/examples/gallery/lines/linestrings.py index ff1588cd731..c9f6ad25f40 100644 --- a/examples/gallery/lines/linestrings.py +++ b/examples/gallery/lines/linestrings.py @@ -29,29 +29,3 @@ fig.plot(data=rivers_asia, pen="1p,steelblue") fig.show() - - -# %% -rivers_australia = rivers.cx[111:155, -40:-9].copy() - -fig = pygmt.Figure() -fig.basemap(region=[111, 155, -40, -9], projection="M10c", frame=True) -fig.coast(land="gray95", shorelines="1/0.3p,gray50") - -# Add rivers to map -fig.plot(data=rivers_australia, pen="1p,steelblue") - -fig.show() - - -# %% -rivers_sa = rivers.cx[-84.5:-33, -56.5:13].copy() - -fig = pygmt.Figure() -fig.basemap(region=[-84.5, -33, -56.5, 13], projection="M10c", frame=True) -fig.coast(land="gray95", shorelines="1/0.3p,gray50") - -# Add rivers to map -fig.plot(data=rivers_sa, pen="1p,steelblue") - -fig.show()