如何用python合并多個有規(guī)則命名的nc文件

更新時間：2022年03月28日 11:39:13 作者：oceanography-Rookie

在地學領域,nc格式的文件可謂隨處可見,這種文件可以存儲多維數(shù)字矩陣,下面這篇文章主要給大家介紹了關于如何用python合并多個有規(guī)則命名的nc文件的相關資料,需要的朋友可以參考下

現(xiàn)有多個nc文件，命名除了年份不同外，其他皆相同。希望將多個的nc文件合并為一個。

其中，每個nc文件內(nèi)包含如下內(nèi)容，包含12個月均數(shù)據(jù)

因此，我們可以知道，合并兩個nc文件，time應該=‘24’，合并三個，time=‘36’。以此類推，合并x年，time=x*12。

主要通過xarray實現(xiàn)合并功能，并寫出合并后的nc文件:

#導入用到的庫
import xarray as xr
import os
#首先，利用一個循環(huán)將需要合并的年份挑選出來。我這里選擇1982-2012年的數(shù)據(jù)
file_xadv=[]
for year in range(1982,2013):    
	fn=os.path.join(path1,'xadv'+'%04.0f'%year+'.mon.161x720.nc')
    file_xadv.append(fn)

運行后可以得到所有完整路徑的文件列表：

下面就是進行合并啦：

#合并nc文件，通過xarray打開nc文件，并得到所需變量的xarray.DataArray，之后利用xarray.concat()實現(xiàn)合并的功能
xadv_new=[]
for i in range(len(file_xadv)):
     xadv=xr.open_dataset(file_xadv[i])['xadv']
     xadv_new.append((xadv))
da=xr.concat(hadv_new,dim='time')

運行發(fā)現(xiàn)合并好啦：

xarray.concat()官方介紹如下：

xarray.concat

例外，官網(wǎng)也給出了其他的合并方式，大家可以按需所求：

1、有關沿單個維度組合數(shù)據(jù)集或數(shù)據(jù)數(shù)組的信息，使用 concatenate。

2、有關合并具有不同變量的數(shù)據(jù)集，使用 merge。

3、有關合并具有不同索引或缺失值的數(shù)據(jù)集或數(shù)據(jù)數(shù)組的信息，使用 combine。

4、有關沿多個維度組合數(shù)據(jù)集或數(shù)據(jù)數(shù)組的信息，使用 combining.multi。

Combining data

附上完整實現(xiàn)代碼：

import xarray as xr
import os

path1='E:\\BaiduNetdiskDownload\\MSE_budget\\xadv\\'
file_xadv=[]


for year in range(1982,2013):
     fn=os.path.join(path1,'xadv'+'%04.0f'%year+'.mon.161x720.nc')
     file_xadv.append(fn)
hadv_new=[]
for i in range(len(file_xadv)):
     xadv=xr.open_dataset(file_xadv[i])['xadv']
     hadv_new.append((xadv))
da=xr.concat(hadv_new,dim='time')

da.to_netcdf('D:\\desktopppp\\sst_olr\\interp_nan\\xadv_interp.nc')#輸出合并后的nc文件

好啦，今天的分析就到這里啦，感興趣的小伙伴快去嘗試一下吧~

到此這篇關于如何用python合并多個有規(guī)則命名的nc文件的文章就介紹到這了,更多相關python合并nc文件內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關文章希望大家以后多多支持腳本之家！

您可能感興趣的文章: