如何从XYZ点数据创建xarray Dataset？遇ValueError报错求助

你当前的问题核心是维度定义和数据维度不匹配：

• 你给数据指定了3个维度['x','y','z']，但point_data1和point_data2都是一维数组（长度为num_points），二者维度数量不一致，导致报错。
• 你的XYZ点数据本质是一维的点集合，每个点对应一组x/y/z坐标，而不是三维网格数据，所以不需要把x/y/z设为独立维度。

修正后的代码 #

import numpy as np
import xarray as xr

num_points = 20
point_locations = np.random.randint(99, size=(num_points, 3))
point_data1 = np.ones(num_points)
point_data2 = np.random.randint(5, size=(num_points,))

# 定义一维维度'points'，对应每个点的索引
ds = xr.Dataset(
    {
        'point_data1': ('points', point_data1),
        'point_data2': ('points', point_data2)
    },
    coords={
        'x': ('points', point_locations[:, 0]),
        'y': ('points', point_locations[:, 1]),
        'z': ('points', point_locations[:, 2])
    }
)

# 查看结果
print(ds)

关键说明 #

• 用'points'作为唯一的维度，长度等于点的数量，所有数据和坐标都绑定到这个维度上。
• x/y/z坐标被定义为'points'维度的变量，每个索引对应一个点的三维坐标值，这样就和一维的点数据匹配了。
• 如果后续需要将点数据转换为三维网格（比如插值），可以再用xarray的插值方法（如interp）处理，但当前场景下直接用一维维度匹配点数据才是正确的做法。

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