如何使用xarray.Dataset.sel函数？基于ds1坐标截取ds2子集

报错原因 #

你遇到的KeyError是因为latitude和longitude是二维坐标变量（维度为(lat, lon)），而sel()方法仅支持索引一维的维度坐标（即带*标记的、作为数据集维度本身的坐标，比如time、lat、lon），无法直接用二维坐标做索引。

步骤1：对齐时间维度 #

你已经完成了这一步，确保两个数据集时间维度一致：

ds1 = ds1.sel(time=ds2.time)

步骤2：提取ds2的空间子集（匹配ds1的网格） #

使用xarray的sel_points()方法，它专门用于从二维网格中选择指定经纬度的点，支持最近邻匹配和精度容差设置：

import xarray as xr

# 从ds2中选择与ds1网格经纬度匹配的点，使用最近邻方法
ds2_subset = ds2.sel_points(
    lat=ds1.lat,
    lon=ds1.lon,
    latitude=ds1.latitude,
    longitude=ds1.longitude,
    method='nearest',
    tolerance=1e-2  # 根据数据精度调整，比如0.01度，需与经纬度单位一致
)

关键参数说明 #

• method='nearest'：启用最近邻匹配，找到ds2中距离目标经纬度最近的网格点
• tolerance：解决float32和float64的精度差异问题，设置允许的匹配误差范围（比如你的经纬度是0.01度分辨率，就设为0.01）

执行后，ds2_subset的维度会和ds1完全一致：time(364)、lat(252)、lon(802)。

内容的提问来源于stack exchange，提问作者DZN