Flipped axis when plotting geographic data using `shade()`

Hi!

Thanks for a very useful library! I am aware that there have been issues and discussion about the flipped axes when plotting `xarray.DataArrays` (e.g. [here](https://github.com/holoviz/hvplot/issues/668), [here](https://github.com/holoviz/datashader/pull/1095) and [here](https://github.com/holoviz/datashader/issues/1054)) and although there has been work to solve this, I believe there is still an issue when using `shade()` function, i.e. the output image is flipped when using the `shade()` function as shown below. The latter one is correct and the first one seems to be flipped:


### Setup

Python 3.10.13
xarray 2023.10.1 pyhd8ed1ab_0 conda-forge
xarray-spatial 0.3.7 pyhd8ed1ab_0 conda-forge
datashader 0.16.0 pyhd8ed1ab_0 conda-forge

Is there any workaround for this, or ongoing efforts to fix this issue?

Start by reproducing the geographic xarray.DataArray case described in the issue and inspect the shade() entry point for axis-orientation handling. Compare the shaded output with the correct reference image; done means geographic data is rendered without the flipped axis and the behavior is covered by an appropriate regression check.

Written by the indexing model from the issue text.