-
0 votesanswersviews
使用scipy.stats来适应xarray DataArray
我想计算在xarray.DataArray的时间维度上拟合的统计分布的参数 . 我想创建一个类似于以下内容的函数： from scipy import stats import xarray as xr def fit(arr): return xr.apply_ufunc(stats.norm.fit, arr, ...) 返回一个新的DataArray，存储在时间维度上计算的分布的... -
3 votesanswersviews
xarray中的布尔索引
我有一些带有dims 'time', 'lat', 'lon' 的数组，有些只有 'lat', 'lon' . 为了使用2d（lat-lon）掩码屏蔽与时间相关的数据，我经常需要这样做： x.data[:, mask.data] = np.nan 当然，计算按预期广播 . 如果 y 是2d lat-lon数据，则其值将广播到x中的所有时间坐标： z = x + y 但索引并没有预期 . 我... -
1 votesanswersviews
计算大型数据集的时间平均值
我试图获得一个相对较大的数据集的小时的时间平均值 . 你能告诉我如何更快地完成这个过程吗？ 下面是xr.dataset“Tair”的维度 <xarray.Dataset> Dimensions: (lat: 2600, lon: 2600, time: 7200) Coordinates: * lon (lon) float64 -64.97 -64.9... -
1 votesanswersviews
通过Nio.open_file从grib2文件推断的变量名 - 或者带引擎pynio的xr.open_dataset
我正在使用找到的数据here（注意，这些是预测数据的旋转文件，因此实际日期会随着时间的推移而改变，您可能需要在几天内更新我的示例中的日期） !wget ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/cfs/prod/cfs/cfs.20180524/00/6hrly_grib_01/pgbf2018052400.01.2018052400.grb2 ... -
0 votesanswersviews
在网格数据上移动百分位数[暂停]
需要紧急帮助 . 我有一个网格数据的时间lon lat与结构（324,72,144），我想计算时间轴上的移动窗口百分位像这样：一个lop来计算从第一个120时间步骤开始的第n个百分位数，下一个lop从时间2开始，也计算120的相同百分位数，第三个循环从3开始，依此类推 . 任何帮助都感激不尽 . 谢谢大家 . -
2 votesanswersviews
将'constant'维添加到xarray数据集
我有一系列CSV格式的月度网格化数据集 . 我想阅读它们，添加几个维度，然后写入netcdf . 我以前在使用xarray（xray）方面有很棒的经验，所以我想在这个任务中使用它 . 我可以很容易地将它们变成2D DataArray，例如： data = np.ones((360,720)) lats = np.arange(-89.75, 90, 0.5) * -1 lngs = np.aran... -
6 votesanswersviews
将维度添加到xarray DataArray
我需要为 DataArray 添加维度，填充新维度的值 . 这是原始数组 . a_size = 10 a_coords = np.linspace(0, 1, a_size) b_size = 5 b_coords = np.linspace(0, 1, b_size) # original 1-dimensional array x = xr.DataArray( np.random... -
0 votesanswersviews
使用xarray获取网格中每个点的时间总和
我有一个三维网格数据集：lat，lon和time . <xarray.Dataset> Dimensions: (lat: 185, lon: 155, time: 8760) Coordinates: * lon (lon) float64 -76.98 -76.92 -76.86 -76.8 -76.74 ... * lat (... -
1 votesanswersviews
python-xarray从一个DataArray复制掩码到另一个
我让这个用于一个简单的案例： arr2 = xr.DataArray((np.arange(16)-8).reshape(4, 4), dims=['x', 'y']) arr3 = xr.DataArray(np.arange(16).reshape(4, 4), dims=['x', 'y']) <xarray.DataArray (x: 4, y: 4)> array([[ na... -
0 votesanswersviews
xarray从Pandas写入netCDF - 维度问题
学习如何使用xarray从Pandas DF生成netCDF文件 . 以下几个教程和SO问题Add 'constant' dimension to xarray Dataset和Add 'constant' dimension to xarray Dataset但仍有一些问题，因为我无法获得Date_Time，lat和lon作为维度 . 当我进行nc转储时，它们不正确 . 将txt文件导入pand... -
1 votesanswersviews
以块的形式创建xarray DataArray并将其写入NetCDF
是否也可以创建一个核心外的DataArray，并使用xarray将其逐块编写到NetCDF4文件中？ 例如，我希望能够在尺寸更大的情况下以非核心方式执行此操作，因此我无法将整个数组存储在内存中： num_steps = 20 num_times = 100 #Create DataArray d = xr.DataArray(np.zeros([num_steps, num_times], np.... -
1 votesanswersviews
xarray Multiindex concat的最佳实践
我有一组1000（2D） pd.Dataframe （比方说，索引：时间，列：run_id），每一个都有3个属性（比方说温度，压力，位置） . 理想情况下，我希望将 xr.DataArray 中的所有内容都包含在5个维度中（或者具有4个维度的 xr.Dataset 并将最后一个维度作为唯一数据变量） . 我创建了一个带有两个dims和2个3 coords的DataArray但是 xr.concat... -
1 votesanswersviews
使用xarray，如何在多维数据集上并行化1D操作？
我有一个4D xarray数据集 . 我想在特定维度（此处为时间）上对两个变量进行线性回归，并将回归参数保存在3D数组中（其余维度） . 我设法通过使用此序列代码获得了我想要的结果，但它相当慢： # add empty arrays to store results of the regression res_shape = tuple(v for k,v in ds[x].sizes.items... -
0 votesanswersviews
xarray DataArray.where（）在屏蔽时减少了坐标
xarray新手在这里 . 非常简单的情况，我有一个降水型阵列（ntim x nlat x nlon）和一个总降水阵列（相同的尺寸） . 两者都在单独的netCDF文件中 . 我想掩盖降水阵列，其中A）降水量下降（> 1e-8 m / s率）和B）降水类型为雪（maskvar = 0.0） . 因此输出数组是“下雪的地方？”阵列 . 当使用xarray where（）具有来自两个不同（但相同...