xarray 空索引 isel 报 IndexError，检查后端轴是否错位

对磁盘上的多维变量调用 isel(x=[], y=[1, 3])，表达式建立时可能正常，直到访问 .values 才抛出 IndexError。空选择本应返回长度为零的维度；在一条已确认的 xarray 后端路径中，空索引拆分时漏了一项，后续索引便对错了轴。

这条问题来自 2026 年 9 月的开发版本，多维数据使用 h5netcdf、scipy 引擎时已有记录。不能据此说所有 xarray 版本或所有空数组都失败；正式 2026.09.0 已包含相关修补。

错误可能延迟到真正读取时

xarray 会把一次选择拆成后端读取和内存中的二次索引。对于只支持一个向量索引的后端，空向量可转成后端的空切片，但内存索引仍必须保留相应轴的位置。

受影响分支完成了前半步，却没有为这一轴追加 slice(None)。一维数据可能看不出问题；多维变量再叠加其他索引，就可能出现越界。更多维度还可能暴露索引元组长度不一致的 ValueError。

因此应记录变量原始维度顺序、各索引的形状和类型、引擎及版本，以及异常是在选择时还是加载时发生。文件可以正常打开，不足以排除索引组合阶段的问题。

使用含修补的实现，或把空选择写成明确切片

在依赖约束允许时，选择已修补的正式版本，并保留原来的引擎和选择条件复查：

python -m pip install "xarray==2026.09.0"

这不是要求同时无条件升级所有科学计算依赖。应让项目的 Python、NumPy、后端和锁文件保持兼容，并核对实际执行环境。

暂时保留受影响构建时，若业务明确要在 x 轴选择零个元素，可以将这个空向量写成空切片：

import xarray as xr

with xr.open_dataset("input.nc", engine="h5netcdf") as dataset:
    selected = dataset.isel(x=slice(0, 0), y=[1, 3])
    result = selected.load()

示例要求文件含有 x、y 维度，且 y 至少有四个位置。路径和引擎换成原有输入；在上下文关闭前加载的是选择后的结果，而不是整个文件。若其他维度也可能产生空选择，需要按它们各自的轴含义处理。

不要把空选择改成 x=0：那会读入第一项并可能移除维度，改变业务结果。直接对整个大数据集调用 .load() 虽可能绕过后端懒加载，却会改变内存需求，也没有修复原来的选择语义。

验收同时比较空选择和非空选择。空结果应保留约定的维度与坐标，x 长度为零，其他维度长度符合索引；再触发实际加载，确认没有延迟异常。非空对照还要核对值和轴顺序，不能只检查返回对象存在。不同后端可能呈现不同错误，甚至返回错误尺寸，不能把“没有抛异常”当成它们都正确。