我确实在一个目录中有一堆包含大气测量值的文件。文件格式为 NetCDF。每个文件都有一个时间戳(变量“basetime”)。我可以读取所有文件并绘制单个测量事件(温度与高度)。接下来我需要做的是“按天对文件进行分组”,并将一天中的所有测量值一起绘制在一个图中。不幸的是,我不知道该怎么做。一种想法是使用下面代码中定义的变量“measurement_day”。对于每一天,我通常确实有四个不同的文件,其中包含 temp. 和海拔高度。理想情况下,这四个不同文件的数据应该分组(例如用于绘图)我希望我的问题很清楚。谁能帮帮我吗。
编辑:我现在尝试使用字典,但我无法确定一个测量日是否已经存在一个条目。请参阅下面的编辑代码
from netCDF4 import Dataset
data ={} # was edited
for f in listdir(path):
if isfile(join(path,f)):
full_path = join(path,f)
f = Dataset(full_path, 'r')
basetime = f.variables['base_time'][:]
altitude = f.variables['alt'][:]
temp = f.variables['tdry'][:]
actual_date = strftime("%Y-%m-%d %H:%M:%S", gmtime(basetime))
measurement_day = strftime("%Y-%m-%d", gmtime(basetime))
# check if dict entries for day already exist, if not create empty dict
# and lists inside
if len(data[measurement_day]) == 0:
data[measurement_day] = {}
else: pass
if len(data[measurement_day]['temp']) == 0:
data[measurement_day]['temp'] = []
data[measurement_day]['altitude'] = []
else: pass
我收到以下错误消息: Traceback (most recent call last):... if len(data[measurement_day]) == 0: KeyError: '2009/05/28'