本文介绍使用 xarray 的 `decode_cf()` 自动解析 netcdf 中的 cf 时间编码，并通过 `.astype(int)` 直接向量化转换为纳秒级 unix 时间戳，再换算为毫秒——全程无需循环，一行代码即可完成批量转换。
在处理气象、气候等 NetCDF 数据时，时间维度常以 cftime 编码形式存储（如 “days since 2002-01-01T00:00:00″），其原始值为浮点数数组，并非真正的 datetime 对象。许多用户误以为需手动遍历 + cftime.num2date() + datetime.timestamp() 等多步转换，既低效又易出错。实际上，xarray 提供了开箱即用的矢量化解决方案。
核心步骤如下：
- 确保变量属性符合 CF 标准：units 属性必须为小写 “days since …”（而非 “Days Since …” 或 “Units: days since…”），且 calendar（如 “gregorian”）存在；
- 调用 xr.decode_cf()：自动识别 units 和 calendar，将数值型时间数组解码为 datetime64[ns] 类型的 xarray DataArray；
- 直接类型转换获取 Unix 毫秒时间戳：datetime64[ns] 底层以纳秒为单位自 Unix 纪元（1970-01-01T00:00:00Z）起计，因此 .astype(int) 得到纳秒级整数，除以 1_000_000 即得毫秒级 Unix 时间戳。
✅ 示例代码（完整可运行）：
import xarray as xr
import numpy as np
# 模拟原始 NetCDF 时间数据（注意：units 必须小写！）
attrs = {
"units": "days since 2002-01-01T00:00:00",
"calendar": "gregorian"
}
time_data = np.array([107.0, 129.5, 227.5, 7928.0, 7958.5, 7989.0], dtype=np.float32)
ds = xr.Dataset({"time": ("time", time_data, attrs)})
# ✅ 一步解码为 datetime64[ns]
decoded = xr.decode_cf(ds)
print("解码后时间（datetime64[ns]）:")
print(decoded.time.values)
# ✅ 向量化转 Unix 毫秒时间戳（无循环！）
unix_ms = (decoded.time.astype(np.int64) // 1_000_000).values
print("\n对应 Unix 时间戳（毫秒）:")
print(unix_ms)
输出示例：
