nc数据世界时转北京时间，请高手切磋，写的菜

主啊

import xarray as xr
import netCDF4 as nc
import datetime as dt
import pandas as pd
from datetime import datetime, timedelta

timefile = 'E:/xinzheng2021-2022/2021year01-03day0-23h.nc'
data = xr.open_dataset(timefile)
f = []

for i in data["time"].values:
    print("before",i)
i = dt.datetime.strptime(str(i)[0:-10], '%Y-%m-%dT%H:%M:%S')
    i = i + dt.timedelta(hours=8)
    print("after")
    print(i)
    i = i.strftime('%Y-%m-%dT%H:%M:%S')
    f.append(i)

print("*"*60)
print(f)

time = f
tsc = data.Temperature_height_above_ground
lat = tsc.lat
lon = tsc.lon
tt = xr.DataArray(tsc, coords=[time,lat,lon], dims=['time','lat','lon'])
print(tt.time)

tt.to_netcdf('E:/xinzheng2021-2022/2021-12(beijingtime).nc')

付亚男

1. import xarray as xr
   
2. import pandas as pd
   
3. 
   
4. fnm = r'F:/meteo_data/era5/ear5_single_levels_2021071400_2021072218.nc'
   
5. ds = xr.open_dataset(fnm)
   
6. 
   
7. t_bst = [pd.Timestamp(t, tz='UTC').tz_convert('Asia/Shanghai') for t in ds.time.values]
   
8. 
   
9. ds_new = ds.assign_coords({'time': t_bst})

复制代码

雨落森林

付亚男 发表于 2023-6-13 17:37

原来还可以这么用，学到了！！！

Nicko叔叔

用MATLAB写的：

1. function [t_new] = nctime(file_path,name)
   
2. % 将nc文件里的时间转变为正常使用的时间
   
3. % 输入文件路径和时间的变量名
   
4. % 输出yyyy-mm-dd HH:MM:SS格式的时间
   
5. 
   
6. nc_time=double(ncread(file_path,name));% 读取nc时间;
   
7. time_info=ncreadatt(file_path,name,'units');%读取nc时间的信息;
   
8. 
   
9. %% 初始时间和单位
   
10. blanks_location = find(time_info==' ');%从time_info中找到空格的位置，用于分开单位，年月日，时分秒
    
11. t0 = datetime(time_info(blanks_location(2)+1:end)); %找到起始时间并转化为datetime格式
    
12. 
    
13. time_unit = time_info(1:blanks_location(1)-1);%找到时间单位
    
14. 
    
15. %% 处理时间
    
16. switch time_unit
    
17.         % 如果单位是days，那么所读取到的数字nc_time的单位就是天
    
18.         % 所以可以将nc_time转换为calduration
    
19.         % 然后加在t0上，就成为了实际时间
    
20.         % 其他的同理
    
21.         case 'years'
    
22.                 tt = calyears(nc_time);
    
23.         case 'months'
    
24.                 tt = calmonths(nc_time);
    
25.         case 'days'
    
26.                 tt = caldays(nc_time);
    
27.                 % 因为小时、分钟、秒是没有cal属性的，所以只能用hours、minutes、seconds
    
28.         case 'hours'
    
29.                 tt = hours(nc_time);
    
30.         case 'minutes'
    
31.                 tt = minutes(nc_time);
    
32.         case 'seconds'
    
33.                 tt = seconds(nc_time);
    
34. end
    
35. t_new = tt + t0;%返回的是一个datetime格式的时间表
    
36. 
    
37. end

复制代码

yangzi88552

1.     ds = xr.open_dataset(nc_file_path)
   
2.     time = ds.time.values+ np.timedelta64(8, 'h')

复制代码

yangzi88552

yangzi88552 发表于 2023-6-14 11:22

原始世界时：['2023-05-24T00:00:00.000000000' '2023-05-24T01:00:00.000000000'
'2023-05-24T02:00:00.000000000' '2023-05-24T03:00:00.000000000'
'2023-05-24T04:00:00.000000000' '2023-05-24T05:00:00.000000000']

转换后的北京时：['2023-05-24T08:00:00.000000000' '2023-05-24T09:00:00.000000000'
'2023-05-24T10:00:00.000000000' '2023-05-24T11:00:00.000000000'
'2023-05-24T12:00:00.000000000' '2023-05-24T13:00:00.000000000']

主啊

yangzi88552 发表于 2023-6-14 11:25
原始世界时：['2023-05-24T00:00:00.000000000' '2023-05-24T01:00:00.000000000'
'2023-05-24T02:00:0 ...

'2023-05-24T00:00:00.000000000'

这个对应的这个格式我写也写不对，%Y-%m-%dT%H:%M:%S
不知道后面那么多的0是怎么表示

yangzi88552

主啊 发表于 2023-6-14 11:35
'2023-05-24T00:00:00.000000000'

这个对应的这个格式我写也写不对，%Y-%m-%dT%H:%M:%S

你把格式转换一下就可以了

1. time_new = pd.to_datetime(time, format='%Y-%m-%d-%H')

复制代码

yangzi88552

主啊 发表于 2023-6-14 11:35
'2023-05-24T00:00:00.000000000'

这个对应的这个格式我写也写不对，%Y-%m-%dT%H:%M:%S

DatetimeIndex(['2023-05-24 08:00:00', '2023-05-24 09:00:00',
               '2023-05-24 10:00:00', '2023-05-24 11:00:00',
               '2023-05-24 12:00:00', '2023-05-24 13:00:00',
              dtype='datetime64[ns]', freq=None)


这种就是常见的样式了

主啊

yangzi88552 发表于 2023-6-14 14:51
DatetimeIndex(['2023-05-24 08:00:00', '2023-05-24 09:00:00',
               '2023-05-24 10:00:00' ...

看到这个数据类型，还有个问题想请教一下：



time = pd.date_range(start="20210101", end="20210101", freq="MS")
print(time)
上面这个就是生成了这种格式：DatetimeIndex(['2021-01-01'], dtype='datetime64[ns]', freq='MS')




但是我想用这个时间，却不知道怎么转换为这DatetimeIndex，用于替换nc文件的时间，下面这样就会报错：ValueError: dimensions ('time',) must have the same length as the number of data dimensions, ndim=0



i = dt.datetime.strptime(str(tsc["time"].values)[2:-12], '%Y-%m-%dT%H:%M:%S')
i = i + dt.timedelta(hours=8)
i = i.strftime('%Y-%m-%dT%H:%M:%S')
time = pd.to_datetime(i)
tt = xr.DataArray(tsc, coords=[aa,lat,lon], dims=['time','lat','lon'])
（这个nc文件只有一个时间，替换为北京时间为20210101年0点的）