引入相关包和导入数据:
import numpy as np
import xarray as xr
from matplotlib import pyplot as plt
# 数据导入
path = "...\\sst.mnmean.nc"
# 丢弃一个不必要导入的变量
ds = xr.open_dataset(path, drop_variables=["time_bnds"])
ds = ds.sel(time=slice("1960", "2018")).load()
对于数组的赋权降维而言,xarray目前支持DataArray和Dataset对象,对于这两个对象的降维可采用DataArray.weighted()
和Dataset.weighted()
方法。目前支持带权重的以平均(mean
)和求和(sum
)方法降维。
为说明赋权降维,下面先创建一个关于降水数据的DataArray和一个权重的DataArray.
prec = xr.DataArray(
[50, 10, 0.9],
dims=("month"),
coords= {"month":("month", [1, 2, 3])}
)
prec