grdimage fails when attempting to plot grids that are not "global" [0, 360, -90, 90].

The xarray dataset I am trying to plot is equispaced in latitude and longitude, and normally would work with gmt. When trying to plot with pygmt however, the routine goes into some kind of infinite loop and never terminates.

The xarray dataset has the following characteristics:

Coordinates:
  * latitude   (latitude) float64 90.0 80.0 70.0 60.0 ... -60.0 -70.0 -80.0
  * longitude  (longitude) float64 0.0 10.0 20.0 30.0 ... 330.0 340.0 350.0

and I tried to plot the data using

fig.grdimage(grid, region=[0., 360., -90, 90], projection="W0/6i")

I note that the redundant data points at 360 E are not provided, and the data point at the south pole is not provided. I further note that this is the standard format for Driscoll and Healy sampled grids which are employed in pyshtools.