当下载的时间序列数据(比如温度、降水、海面高度等)为多个nc文件时,有时我们需要将nc文件进行多维合并,方便后续进行分析(比如EOF正交函数分解等),可以通过Python的xarry包和glob包简单实现:
import
# -*- coding: UTF-8 -*-
import glob
import xarray
文件读取,查看读取结果
path='D:/coastlines/sealevel/nc/'
all_files=glob.glob(path+'*.nc')
all_files.sort()
print(all_files[0])
循环合并nc数据
file_new=[]
for i in range(len(all_files)):
file=xarray.open_dataset(all_files[i])['sla']
file_new.append(file)
da=xarray.concat(file_new,dim='time')
da.to_netcdf('D:/coastlines/sealevel/all_file2.nc')
查看合并后的维度:
with xarray.open_dataset(path+'all_file2.nc') as f:
sla=f.sla
print(sla)
![](https://img-blog.csdnimg.cn/34b25e94dbb24e9e8dec5d3c934f7a13.png)
|