已知各时刻午夜经度，如何计算不同经度的当地时间？

首先明确核心逻辑：已知某一时刻（对应time数组中的元素），loc0（即lon_lt0）是该时刻当地时间为0点的经度，那么任意经度lon的当地时间，可通过经度差对应的时区偏移计算：

• 地球每15经度对应1小时时差，往东（经度更大）时区更靠前，往西则更靠后
• 当地时间 = 该时刻的UTC时间 + (当前经度 - loc0)/15 小时

原代码的错误点 #

• 用len(lon)/360*24计算偏移量：len(lon)是经度的数量，和实际经度数值无关，完全不符合时区计算逻辑
• 未处理loc0随时间变化的特性：原代码没有遍历时间维度，直接用整个loc0数组计算
• 循环逻辑错误：每次传入整个lon数组到计算函数，没有针对单个经度和对应时间点的loc0做计算

正确实现代码 #

from netCDF4 import Dataset, num2date
import numpy as np

# 读取数据集
ws = Dataset('dataset.nc')

# 提取变量
lon = ws.variables['longitude'][:]       # 一维数组：所有经度值
time = ws.variables['time'][:]           # 一维数组：所有时间点
loc0 = ws.variables['lon_lt0'][:]        # 一维数组：每个时间点对应的当地0点经度

# 将数值型时间转换为UTC datetime对象（根据数据集的time单位自动适配）
utc_time = num2date(time, units=ws.variables['time'].units)

# 用numpy广播实现向量化计算，避免低效循环
# 扩展维度，让time和lon的维度匹配（形状变为：n_time × n_lon）
lon_expanded = np.tile(lon, (len(utc_time), 1))
loc0_expanded = np.tile(loc0[:, np.newaxis], (1, len(lon)))

# 计算每个点的时区偏移小时数
hour_offsets = (lon_expanded - loc0_expanded) / 15

# 转换为numpy datetime64格式，方便时间运算
utc_time_np = np.array(utc_time, dtype='datetime64[h]')
# 计算所有经度、所有时间点的当地时间
localtime = utc_time_np[:, np.newaxis] + np.array(hour_offsets, dtype='timedelta64[h]')

# 示例：打印第一个时间点下所有经度的当地时间
print(f"UTC时间: {utc_time[0]}")
for idx, lon_val in enumerate(lon):
    print(f"经度 {lon_val}: 当地时间 {localtime[0, idx]}")

关键说明 #

• 如果你的time变量已经是datetime对象，可跳过num2date转换步骤
• 代码使用numpy向量化运算替代循环，效率远高于逐元素遍历
• 完全不考虑夏令时，符合需求

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