holoviews: Slicing bug with xarray

I found a subtle bug when slicing Images created with xarray, where an incorrect position is returned. Unfortunately, I am unable to demonstrate this with synthetic data, but the following notebook illustrates the problem:

https://gist.github.com/drs251/f97801aa7d6f92b11c9c71c6a7c4085a

EDIT: I found out how to reproduce it: the key is that the frames have to be non-square:

import numpy as np
import xarray as xr
import holoviews as hv
hv.notebook_extension()

x = np.linspace(-3, 7, 533)
y = np.linspace(-5, 8, 894)
z = np.exp(-1*(x**2 + y[:, np.newaxis]**2))
p = np.arange(5)
zz = np.dstack([z, z, z, z, z])
array = xr.DataArray(zz, coords=[y, x, p], dims=['x', 'y', 'p'])
array.name = "intensity"

image = hv.Dataset(array).to(hv.Image, kdims=["x", "y"])
image[0, -1:1, -1:1]

The slice should be nicely centered on the maximum, but it’s not.