Description of the desired feature

Originally posted by @leouieda as a comment about xarray integration at https://hackmd.io/@pygmt/SkVmX_20d:

Any thoughts on better integration with xarray? Dataset objects can have rich metadata (coordinate names, units, projections, etc) that could be leveraged to automatically add labels and titles to plots. Their plotting backend does that with matplotlib and it's very convenient.

My interpretation of this 'feature request' is that grdimage could automatically add x-axis and y-axis labels/annotations to the plot, and include a title about the z-variable name (without the user setting it explicitly). Something based on pydata/xarray#519.

The auto-labels/titles currently comes for free in xarray through the xarray.DataArray.plot.imshow() function. We could inspect the source code at https://github.com/pydata/xarray/blob/v0.19.0/xarray/plot/plot.py#L244 to see how this can be implemented in PyGMT.

Are you willing to help implement and maintain this feature? Long term stuff for PyGMT > v0.6.0