Add DataArray.gmt.imshow() function to show 2D DataArray easily

pygmt/accessors.py

@@ -5,6 +5,7 @@
 
 import xarray as xr
 from pygmt.exceptions import GMTInvalidInput
+from pygmt.figure import Figure
 from pygmt.src.grdinfo import grdinfo
 
 
@@ -165,3 +166,25 @@ def gtype(self, value):
                 "either 0 for Cartesian or 1 for Geographic."
             )
         self._gtype = value
+
+    def imshow(self, fig=None):
+        """
+        Image plot of 2D :class:`xarray.DataArray`.
+
+        Examples
+        --------
+
+        >>> from pygmt.datasets import load_earth_relief
+        >>> grid = load_earth_relief()
+        >>> grid.gmt.imshow()
+        """
+        newfig = False
+        if fig is None:
+            fig = Figure()
+            newfig = True
+        fig.grdimage(self._obj, frame=True)
+        fig.colorbar(frame=True)
+
+        if newfig:
+            fig.show()
+            return fig