Xarray - 连接多个文件的切片

时间:2018-05-30 18:47:22

标签: netcdf python-xarray xarray

我正在尝试将多个文件的片段连接成一个文件(由零数组初始化),然后写入nCDF文件。但是,我收到错误:

arguments without labels along dimension 'Time' cannot be aligned 
because they have different dimension sizes: {365, 30}

我理解错误(isel()将维度的大小更改为切片的大小),但我不知道如何纠正或避免问题。我正确地接近这个任务了吗?这是第一次迭代的简化版本:

import xarray as xr
import numpy as np

i=0

PRCP = np.zeros((365,327,348))

d = xr.open_dataset("/Path")

d = d.isel(Time=slice(0,-1,24))

P = d['CUMPRCP'].values

DinM = P.shape[0]

PRCP[i:i+DinM,:,:] = P

i = i + DinM

PRCPxr = xr.DataArray(PRCP.astype('float32'),dims=[('Time'), 
'south_north', 'west_east']) 

d['DPRCP'] = PRCPxr

1 个答案:

答案 0 :(得分:0)

通过从xr.DataArray()中删除dims =()参数解决了问题,并在其中任意重命名了它们。