python批量读取nc，按时间维进行合并

lt5201314

每次做这个都要翻以前的代码，来这里记录一下（我目前只会用python这一种方法，如果有其他行之有效的，求大佬们指点）




#导入需要的库包
import xarray as xr
import os
import numpy as np
path = 'h:/data/olr_1979-2018/'  #输入文件存储路径
folder = os.walk(path)  #生成器，里面包含三个东西：根目录，根目录下的目录和文件；生成器只能用list读取
files = list(folder)[0][2]  #选择文件名集合的list
data_all = []
a=[]
b=[]
for file in files:
    f = xr.open_dataset(path+'/'+file)
    #读取变量，如果数据过大可以索引，但索引后合成的nc文件无法再索引
    lon=f['lon']
    lat=f['lat']
    time=f['time']
    olr=f['olr']

#concat函数按维度合并   
    for t in time:
        a.append(t)
    t_al=xr.concat(a,dim='time')

    for olr1 in olr:
        b.append(olr1)
    olr_al=xr.concat(b,dim='time')
print(t_al.shape)
#写成dataarray
OLR=xr.DataArray(data=olr_al,dims=('time','lat','lon'),
                coords=dict(time=t_al,lat=lat,lon=lon))
#dataarray转化为dataset
ds=xr.Dataset({'olr':OLR})
#输出成nc文件
ds.to_netcdf(r'h:/data/olr_al.nc')

lt5201314

声明：代码不是原创，我也忘记了我在哪里看到的了，非常非常nice~

扫眉才子

谢谢你，帮大忙了

Bravery

这里的合并是什么意思呀，sum()吗

独孤酒见

这个脚本我试了下，发现因为我的却缺time，所以代码没法执行。

我的nc文件形式如下：
GWRwSPEC_CHi_200701_200712-RH35.nc
里面变量：
LAT
LON
PM25

有三个变量，但是没有时间，时间在文件名里（如上面的200701_200712）.


每个文件都是这样的，请问下这样的nc文件怎么能合啊。

锦城の烟雨

有用，感谢！

love2100

学习一下啦

gaoyuan215

太感谢了！真的是太需要了！

lt5201314

Bravery 发表于 2021-11-24 09:13
这里的合并是什么意思呀，sum()吗

这里的合并是针对时间维的，如果原始数据是10年的（time, lev, lat, lon）(365,10,10,10)，时间维是逐日的，那合并后的数据集就是（time, lev, lat, lon）(3650,10,10,10)

lt5201314

独孤酒见 发表于 2021-11-24 10:12
这个脚本我试了下，发现因为我的却缺time，所以代码没法执行。

我的nc文件形式如下：

可以建立一个四维数组，把你每个文件的数据循环存进去