import xarray
import os
path = "/copy_merge/"
comb_ds = xr.Dataset()
for root, dirs, files in os.walk(path):
for file in files:
if file.endswith(".nc"):
ds = xr.open_dataset(os.path.join(root, file))
tSK_var = ds['TSK']
comb_ds[file] = tSK_var
comb_ds.to_netcdf("combined_data1.nc")
path = "combined_data1.nc"
ds= xr.open_dataset(path)
df=ds.to_dataframe()
df
df.to_csv("output1.csv")
对经纬度范围和时间进行切片保存
import xarray as xr
import os
path = "/copy_merge/"
#创建空的dataset
comb_ds = xr.Dataset()
#遍历路径中的nc文件,并提取变量
for root, dirs, files in os.walk(path):
for file in files:
if file.endswith(".nc"):
ds = xr.open_dataset(os.path.join(root, file))
rain = ds['RAINNC'][1:,::]
comb_ds[file] = rain
print(comb_ds)
comb_ds.to_netcdf("combined_data1.nc")
path = "combined_data1.nc"
ds = xr.open_dataset(path)
ds1=ds.sel(south_north=slice(9,108), west_east=slice(9,88))
ds1.to_netcdf("wuwei_RAINNC.nc")
df = xr.open_dataset("wuwei_RAINNC.nc")
df1=df.to_dataframe()
df1.to_csv("RAINNC.csv")