使用 xarray 将 NetCDF 数据集的纬度分辨率更改为 0.5 度

以下代码展示了如何使用 xarray 库将 NetCDF 数据集的纬度分辨率更改为 0.5 度，同时保持经度分辨率不变。

import xarray as xr
import numpy as np

# 读取数据
path = '/home/yaodianwen/tp_era5_3hrsum_mekong_mm_defined_invalid_pr.nc'
data = xr.open_dataset(path)

# 获取经纬度变量名
lon_var, lat_var = list(data.coords.keys())[:2]

# 确定新的纬度范围和间隔
new_lat = np.arange(data[lat_var].max().data, data[lat_var].min().data, -0.5)

# 更新纬度坐标
data_resampled = data.assign_coords(latitude=new_lat)

# 保存结果
new_path = '/home/yaodianwen/tp_era5_3hrsum_mekong_mm_defined_invalid_pr_interp.nc'
data_resampled.to_netcdf(new_path)

这段代码使用 assign_coords 方法来更新纬度坐标，从而将纬度从 0.25° 变为 0.5°，而经度保持不变。

注意： 此方法不会进行任何插值，它只是将纬度坐标更改为新的值。如果您需要对数据进行插值，请使用 interp 方法。