问题标签 [xarray]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
5 回答
2424 浏览

python - 从 netcdf 文件中获取每个月的每小时平均值

我有一个 netCDF 文件,其时间维度包含 2 年的每小时数据。我想对其进行平均以获得每个月一天中每个小时的小时平均值。我试过这个:

但我收到此错误:

我怎样才能解决这个问题?如果我这样做:

我没有收到错误,但结果的时间维度为 12(每个月一个值),而我想要每个月的小时平均值,即 12 个月的每个月有 24 个值。数据可在此处获得:https ://www.dropbox.com/s/yqgg80wn8bjdksy/ecmwf_usa_2015.nc?dl=0

0 投票
1 回答
1237 浏览

python - 避免在 xarray 平面网格图中重叠颜色条

在运行上面的代码时,我明白了。数字: 在此处输入图像描述

如何确保颜色条不与图重叠?即使我使用 xarray 默认颜色栏,也会发生重叠。

0 投票
1 回答
638 浏览

python - 获取一段时间内每个季节的最大降雨量(xarray)

我试图找到(DJF, MAM, JJA, SON)10 年内每个季节的最大降雨量。我正在使用 netcdf 数据和 xarray 来尝试这样做。数据包括降雨量(每 3 小时记录一次)、纬度和经度数据。现在我有以下代码:

但是,当我这样做时,输出的形状(4,145,192)表明它在整个期间的每个季节都取最大值。我想要每年每个季节的最大值。换句话说,输出的形状应该类似于(40,145,192)(每年 4 个值 x 10 年)

我已经考虑过尝试使用频率DataSet.resample以及使用time=3M频率来做到这一点,但是它并没有正确地划分月份。如果必须,我可以更改数据集,所以它从正确的位置开始,但我希望有一种更简单的方法,因为已经有一个函数可以正确分组它。

谢谢,如果您需要更多详细信息,请告诉我!

0 投票
0 回答
462 浏览

python - 使用 xarray 重新采样返回:TypeError: invalid type Promotion

我正在尝试使用以下代码进行一些计算:

我回来的完整错误是这样的:

但是,当我摆脱似乎导致问题的行 ( ds1 = ds1.resample(time='24H').sum('time')) 时,它运行良好,并且在尝试重新采样时不会出现任何其他错误。

这让我很困惑,我不确定问题出在哪里。当代码稍作更改时它工作正常的事实让我认为特别是在取总和而不是我的数据或重新采样本身的问题时出了问题。我不确定,也不知道如何解决。

编辑 1:数据详细信息

数据是 3 小时降雨数据,形状为(10,90,144)。我试图找到每天的总数,然后是每年的最大每日价值。

以下是有关ds1何时打印的完整信息,如果这也有帮助。(*切片之前):

0 投票
1 回答
534 浏览

python - 来自 Xarray 的往返 Zarr 数据

使用xarray,我ds.to_zarr()用于将数据集写入 S3,然后xr.open_zarr()查看是否得到相同的数据集。

我的数据集xarray如下所示:

我把它写成zarr使用:

当我尝试使用以下方法读回它时:

我回来了:

如果我设置,我可以读取数据集auto_chunk=False

结果是

但是我是否不了解有关分块和方式的一些知识xarraydask并且zarr应该一起工作?

为了上班,我必须做些什么不同的auto_chunk=True事情?

0 投票
1 回答
3316 浏览

python - 用另一个坐标替换 xarray 坐标

我尝试通过对两个坐标都使用一个坐标来替换 xarray.Dataset 中具有相同长度的两个坐标。我不关心旧坐标或其值;我只是想替换它们。有什么方便的方法吗?

我试过rename了,但它返回一个冲突错误(因为 new_coord 已经存在)。

我是否必须将新坐标显式分配给依赖于旧坐标之一的每个数据变量?由于我有许多数据变量(不仅仅是 a 和 b)和其他坐标,这可能会变得混乱。

0 投票
1 回答
296 浏览

python - 试图在 xarray 上应用 fftconvolve

我有一个 3Dxarray DataArray和一个 1D xarray DataArray。3D 数组的一维与另一维的一维大小相同。我想scipy.signal.fftconvolve沿着这些相应的维度申请。但我得到一个错误。这是一些最小的示例代码:

我想做的是相当于:

但是我收到以下apply_ufunc版本错误:

0 投票
2 回答
1210 浏览

python - 使用时间戳超过 2263 年的 xarray 创建 netcdf

有没有办法使用 xarray 创建时间维度超过 2263 年的 netCDF 文件?

以下是如何创建 netCDF 玩具数据集http://xarray.pydata.org/en/stable/time-series.html

然而,时间维度有一种熊猫日期时间索引,并且不超过 2263,如下所示: https ://github.com/pandas-dev/pandas/issues/13346

0 投票
1 回答
1231 浏览

python - matplotlib/cartopy Mollweide 投影中经度环绕的接缝

我正在绘制一些 NETCDF 数据,这些数据由 Matplotlib 中的 xarray 处理,使用 Cartopy 包装器进行地图投影。当我生成一个填充的等高线图 (contourf) 时,我最终会在零经度处得到一个“接缝”(我的数据经度开始和结束的地方)。是否有一种内插方式在经度[-1] 和经度[0] 之间进行插值,也许通过将数据视为周期性数据或其他方式?

示例请参见附图:

在此处输入图像描述

0 投票
1 回答
470 浏览

pandas - 结合互补的 DataArrays

我有一个具有三个维度的 DataArrays 列表。对于列表中的每个项目,其中两个维度是单个值,但所有项目的组合将产生完整的组合值。

然后我想结合这些免费DataArray的,但到目前为止我尝试的所有方法似乎都不起作用。结果应该是DataArray带有形状|2x3x4|和尺寸的dim1: |2|, dim2: |3|, dim3: |4|

以下不起作用: