issue_comments

@Hoeze this is now implemented in #5692 (stack is not yet refactored so I reproduced your example in a slightly different way):

```python

stacked.isel(observations=1) <xarray.Dataset> Dimensions: (genes: 2) Coordinates: * genes (genes) <U1 'a' 'b' observations object ('c', 'f') individuals <U1 'c' subtissues <U1 'f' Data variables: test (genes) int64 2 2

```

I think the right long-term solution for xarray is to always store separate Variable objects for MultiIndex levels, and only use the MultiIndex for proper indexing. When you index out a single value, the MultiIndex will naturally disappear and you'll be left with a bunch of scalar coordinates, without any special case logic to handle the MultiIndex.

This looks like @crusaderky's third option.

We'll need to finish up the big "explicit indexes" refactor first to make this viable.

Not a regression. I've gone back as far as xarray 0.12 and pandas 0.19 and it's always been like this. I agree it's bad and needs to be fixed though.

The issue is inherited straight from pandas: ```python

df = stacked.test.to_pandas() df

individuals c d
subtissues e f e f genes
a 1 2 3 4 b 1 2 3 4

df.iloc[:, 1]

genes a 2 b 2 Name: (c, f), dtype: int64 ``` I'm not sure if we should write an ad-hoc object in xarray for scalar multiindices.

The alternative is to think of a more systematic solution in pandas, which likely implies creating an ad-hoc subclass of tuple which is basically a pickle-able namedtuple. It must be a subclass of tuple otherwise it will break things for a lot of people around the world (the userbase of pandas is MUCH larger than xarray's). And it must be serializable for obvious reasons.

In both cases, the size of this change is very large.

The third and significantly easier option is that, on sel/isel, xarray should automatically unstack any scalar slices of a multiindex. Meaning that the 'observations' coord would simply disappear, leaving only 'individuals' and 'subtissues'. However, It would carry the problem that, if one cuts a scalar slice and a vector slice from the dimension, he won't be able to concatenate them back together.

@shoyer what's your opinion?

@max-sixty here you go: ```python3 import xarray as xr

print(xr.version)

ds = xr.Dataset({ "test": xr.DataArray( [[[1,2],[3,4]], [[1,2],[3,4]]], dims=("genes", "individuals", "subtissues"), coords={ "genes": ["a", "b"], "individuals": ["c", "d"], "subtissues": ["e", "f"], } ) }) print(ds)

stacked = ds.stack(observations=["individuals", "subtissues"]) print(stacked)

print(stacked.isel(observations=1)) ```

result: <xarray.Dataset> Dimensions: (genes: 2) Coordinates: * genes (genes) <U1 'a' 'b' observations object ('c', 'f') Data variables: test (genes) int64 2 2

Do you have a reproducible example, as per the issue instructions?