Xarray中是否有内置方法将datetime转换为cftime？

你提到在非标准日历（如NoLeap）下，cftime转datetime容易出现异常，因此希望将datetime格式转换为cftime格式，目前用循环的暴力方式实现了需求，其实有更高效的内置方法可以替代：

方法1：使用xarray的CFTimeIndex #

xarray提供了CFTimeIndex类，可直接将datetime64数组转换为指定日历的cftime索引，代码示例：

import pandas as pd
import xarray as xr
from xarray import CFTimeIndex

# 原数据
da = xr.DataArray(
    [1, 2], coords={"time": pd.to_datetime(["2000-01-01", "2000-02-02"])}, dims=["time"]
)

# 转换为NoLeap日历的cftime索引
cftime_time = CFTimeIndex(da.time.values, calendar='noleap')
# 更新数据的时间坐标
da_with_cftime = da.assign_coords(time=cftime_time)

print(da_with_cftime.time)

输出会是CFTimeIndex([2000-01-01 00:00:00, 2000-02-02 00:00:00], dtype='object', calendar='noleap', freq=None)，直接得到符合需求的cftime格式时间序列。

方法2：使用cftime的date2num + num2date #

cftime库本身提供了批量转换的函数，先将datetime64转为数值型时间，再转成指定日历的cftime对象：

import pandas as pd
import xarray as xr
import cftime

da = xr.DataArray(
    [1, 2], coords={"time": pd.to_datetime(["2000-01-01", "2000-02-02"])}, dims=["time"]
)

# 定义时间单位和目标日历
units = 'days since 2000-01-01'
calendar = 'noleap'

# 先把datetime64转为数值
time_num = cftime.date2num(da.time.values, units=units, calendar=calendar)
# 再转成cftime对象
cftime_dates = cftime.num2date(time_num, units=units, calendar=calendar)

print(cftime_dates)

输出：array([cftime.datetime(2000, 1, 1, 0, 0, 0, 0, calendar='noleap', has_year_zero=False), cftime.datetime(2000, 2, 2, 0, 0, 0, 0, calendar='noleap', has_year_zero=False)], dtype=object)

这两种方法都是官方内置实现，无需手动循环，处理大数据量时效率更高，也更符合库的设计规范。

内容的提问来源于stack exchange，提问作者chuaxr