Set OGC:CRS84 projection on blue_marble image

Set the projection system so that the xarray.DataArray image is plotted as a Geographic image rather than a Cartesian image by GMT.
GenericMappingTools · Sep 30, 2024 · 471c6e8 · 471c6e8
1 parent f13d12b
commit 471c6e8
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/pygmt/datasets/earth_daynight.py b/pygmt/datasets/earth_daynight.py
@@ -82,4 +82,7 @@ def load_blue_marble(
         region=region,
         registration="pixel",
     )
+    # If rioxarray is installed, set the coordinate reference system
+    if hasattr(image, "rio"):
+        image = image.rio.write_crs(input_crs="OGC:CRS84")
     return image
diff --git a/pygmt/tests/test_grdimage_image.py b/pygmt/tests/test_grdimage_image.py
@@ -37,7 +37,7 @@ def test_grdimage_image_dataarray(xr_image):
     Plot a 3-band RGB image using xarray.DataArray input.
     """
     fig = Figure()
-    fig.grdimage(grid=xr_image, region="d")
+    fig.grdimage(grid=xr_image)
     return fig