Cannot plot multiindexed (stacked) coordinate as hue variable

import numpy as np
import pandas as pd
import xarray as xr

data = np.random.rand(50,5)

x_idx = np.linspace(0, 50)

mi_idx1 = ['a','b','c','d','e']
mi_idx2 = [1,2,3,4,5]

mi = pd.MultiIndex.from_arrays([mi_idx1,mi_idx2], names=['mi_idx1', 'mi_idx2'])

coords = {
    'x': x_idx,
    'mi': mi
}

da = xr.DataArray(data, coords=coords, dims = ['x', 'mi'])

da.plot(hue='mi')

It appears since version 0.16.0, that plotting with a multindex coordinate as a hue dimension no longer works. In version 0.15.1 I get a plot like this:

However, when upgrading to 0.16.0 or 0.16.1 I get the following traceback:

ValueError                                Traceback (most recent call last)
<ipython-input-1-2de8705a6ec5> in <module>
     19 da = xr.DataArray(data, coords=coords, dims = ['x', 'mi'])
     20 
---> 21 da.plot(hue='mi')

~\anaconda3\envs\datanalysis\lib\site-packages\xarray\plot\plot.py in __call__(self, **kwargs)
    444 
    445     def __call__(self, **kwargs):
--> 446         return plot(self._da, **kwargs)
    447 
    448     # we can't use functools.wraps here since that also modifies the name / qualname

~\anaconda3\envs\datanalysis\lib\site-packages\xarray\plot\plot.py in plot(darray, row, col, col_wrap, ax, hue, rtol, subplot_kws, **kwargs)
    198     kwargs["ax"] = ax
    199 
--> 200     return plotfunc(darray, **kwargs)
    201 
    202 

~\anaconda3\envs\datanalysis\lib\site-packages\xarray\plot\plot.py in line(darray, row, col, figsize, aspect, size, ax, hue, x, y, xincrease, yincrease, xscale, yscale, xticks, yticks, xlim, ylim, add_legend, _labels, *args, **kwargs)
    293 
    294     ax = get_axis(figsize, size, aspect, ax)
--> 295     xplt, yplt, hueplt, xlabel, ylabel, hue_label = _infer_line_data(darray, x, y, hue)
    296 
    297     # Remove pd.Intervals if contained in xplt.values and/or yplt.values.

~\anaconda3\envs\datanalysis\lib\site-packages\xarray\plot\plot.py in _infer_line_data(darray, x, y, hue)
     66 
     67         if y is None:
---> 68             xname, huename = _infer_xy_labels(darray=darray, x=x, y=hue)
     69             xplt = darray[xname]
     70             if xplt.ndim > 1:

~\anaconda3\envs\datanalysis\lib\site-packages\xarray\plot\utils.py in _infer_xy_labels(darray, x, y, imshow, rgb)
    378         y, x = darray.dims
    379     elif x is None:
--> 380         _assert_valid_xy(darray, y, "y")
    381         x = darray.dims[0] if y == darray.dims[1] else darray.dims[1]
    382     elif y is None:

~\anaconda3\envs\datanalysis\lib\site-packages\xarray\plot\utils.py in _assert_valid_xy(darray, xy, name)
    410     if xy not in valid_xy:
    411         valid_xy_str = "', '".join(sorted(valid_xy))
--> 412         raise ValueError(f"{name} must be one of None, '{valid_xy_str}'")
    413 
    414 

ValueError: y must be one of None, 'mi_idx1', 'mi_idx2', 'x'

Setting x='x' does not fix the problem

Reproduce the example with the shown xarray, pandas, NumPy, and Matplotlib setup, then inspect xarray/plot/plot.py and xarray/plot/utils.py around _infer_line_data and _infer_xy_labels. Done means da.plot(hue='mi') no longer raises for the MultiIndex coordinate and restores the behavior shown for version 0.15.1.