我对使用xarrays很陌生。我想就地修改NetCDF文件的属性。但是,内置函数提供了另一个新的数据集。
ds = xr.open_dataset(file_)
# ds has "time" as one of the coordinates whose attributes I want to modify
#here is ds for more clarity
ds
>><xarray.Dataset>
Dimensions: (lat: 361, lev: 1, lon: 720, time: 1)
Coordinates:
* lon (lon) float32 0.0 0.5 1.0 1.5 2.0 ... 357.5 358.0 358.5 359.0 359.5
* lat (lat) float32 -90.0 -89.5 -89.0 -88.5 -88.0 ... 88.5 89.0 89.5 90.0
* lev (lev) float32 1.0
* time (time) timedelta64[ns] 00:00:00
Data variables:
V (time, lev, lat, lon) float32 ...
Attributes:
Conventions: CF
constants_file_name: P20000101_12
institution: IACETH
lonmin: 0.0
lonmax: 359.5
latmin: -90.0
latmax: 90.0
levmin: 250.0
levmax: 250.0
我尝试分配新属性,但是给定了一个新的数据数组
newtimeattr = "some time"
ds.time.assign_attrs(units=newtimeattr)
或者,如果我将此属性分配给数据集变量“ V”,它会向数据集添加另一个变量
ds['V '] = ds.V.assign_attrs(units='m/s')
## here it added another variable V .So, ds has 2 variables with same name as V
ds #trimmed output
>>Data variables:
V (time, lev, lat, lon) float32 ...
V (time, lev, lat, lon) float32 ...
答案 0 :(得分:1)
ds.V.attrs['units'] = 'm/s'
为我工作。类似地,“时间”是一个维度
ds.time.attrs['units'] = newtimeattr
答案 1 :(得分:1)
从xarray文档中,xarray.DataArray.assign_attrs
返回与self.attrs.update(* args,** kwargs)等效的新对象。
这意味着此方法返回具有更新属性的新DataArray(或坐标),并且您必须将它们分配给数据集以便它们进行更新:
ds.time.assign_attrs(units=newtimeattr)
您pointed out时,可以通过使用关键字语法访问attrs来完成此操作:
ds.time.attrs['units'] = newtimeattr
需要澄清的一点是-您最后一条语句添加新变量的原因是因为您将具有更新属性的ds.V
分配给了变量ds['V ']
带空格。由于python中的'V ' != 'V'
,因此在更新属性后,将创建一个新变量并为其分配原始ds.V
的值。否则,您的方法就可以正常工作:
ds['V'] = ds.V.assign_attrs(units='m/s')