python - 我有很多 netcdf 文件,如何使用 xarray 将所有文件上传到一个 python 笔记本中?
问题描述
我有 12 个 netcdf 文件的 1 个文件夹。我如何使用 xarray 将所有 netcdf 文件合并到一个数据数组中?
该文件夹代表 2015 年,12 个 netcdf 文件代表 2015 年每个月的数据。
我想也许我可以通过运行 for 循环来尝试操作字符串并更改字符串中文件的每个编号,因为我的文件是通过以下方式组织的(2015 年):
EN.4.2.1.f.analysis.g10.201501.nc
EN.4.2.1.f.analysis.g10.201502.nc
EN.4.2.1.f.analysis.g10.201503.nc
....
我正在加载一个 netcdf 文件,例如:
import numpy as np
from scipy import stats
import matplotlib.pyplot as plt
import netCDF4 as s
import glob
import xarray as xr
from datetime import datetime
inpath='../../Data/EN.4.2.1_analyses/EN.4.2.1.analyses.g10.2015/'
# just loading in the first month (january) of 2015
en4_2015 = xr.open_dataset(inpath+'EN.4.2.1.f.analysis.g10.201501.nc')
如何将所有 netcdf 文件组合成一个数组而不是 12 个不同的 xarray?
我试过了:
import glob
import xarray as xr
from datetime import datetime
# List all matching files
files = glob.glob(inpath+'*.nc')
# Create list for
individual_files = []
# Loop through each file in the list
for i in files:
# Load a single dataset
timestep_ds = xr.open_dataset(i)
# Create a new variable called 'time' from the `time_coverage_start` field, and
# convert the string to a datetime object so xarray knows it is time data
timestep_ds['time'] = datetime.strptime(timestep_ds.time_coverage_start,
"%Y-%m-%dT%H:%M:%S.%fZ")
# Add the dataset to the list
individual_files.append(timestep_ds)
# Combine individual datasets into a single xarray along the 'time' dimension
modis_ds = xr.concat(individual_files, dim='time')
print(modis_ds)
AttributeError Traceback (most recent call last)
<ipython-input-36-34685efc1691> in <module>
10 # Create a new variable called 'time' from the `time_coverage_start` field, and
11 # convert the string to a datetime object so xarray knows it is time data
---> 12 timestep_ds['time'] = datetime.strptime(timestep_ds.time_coverage_start,
13 "%Y-%m-%dT%H:%M:%S.%fZ")
14
~/miniconda3/envs/py3_std_maps/lib/python3.8/site-packages/xarray/core/common.py in __getattr__(self, name)
226 with suppress(KeyError):
227 return source[name]
--> 228 raise AttributeError(
229 "{!r} object has no attribute {!r}".format(type(self).__name__, name)
230 )
AttributeError: 'Dataset' object has no attribute 'time_coverage_start'
解决方案
轻松修复!利用:
ds = xr.open_mfdataset(file_path_folder...time??/*nc')
时间不是说 2015 年,而是 20 年?获取 21 世纪的所有文件,与 19 世纪相同??获取 20 世纪的所有文件等。
推荐阅读
- r - 如何在 R 中生成具有线性分量的余弦模型?
- java - 查找树高时的stackoverflow
- java - 未处理的异常:java.io.IOException
- python - 如何从现有数据集中生成一组新值?
- macos - 安全性:SecKeychainSearchCopyNext 错误 - 无法在 MacOS 上安装 Delve
- apache-kafka - Log4j 无法正常工作
- python - 如何将管道(正则表达式或)应用于模式的特定子字符串,以便它匹配后跟一组字符或 EOL 的组?
- shell - egrep 在 while 循环中没有给出结果,但是单独工作
- performance-testing - 当负载生成器和应用程序/数据库服务器上没有资源饱和时,测试执行期间的高响应时间
- c# - 阻止相机穿过物体