Example for geopandas line geometry: Use Natural Earth dataset, change to Asia

examples/gallery/lines/linestrings.py

@@ -14,33 +14,19 @@
 """
 
 # %%
-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}/10m/physical/ne_10m_rivers_lake_centerlines.zip")
 
-# 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=[-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=gdf, pen="1p,steelblue")
+fig.plot(data=rivers["geometry"], pen="1p,steelblue")
 
 fig.show()