我正在尝试基于变量从netCDF文件中提取天气数据。 .nc文件包含14个变量和2维。我想提取与第一个变量的值相关的14个变量的所有数据。数据来自荷兰计量学院,可以在here上找到。
使用netCDF4模块在Python中加载数据,如下所示:
import netCDF4 as nc
filename = r'path/file.nc'
dataset = nc.Dataset(filename)
打印的变量和尺寸:
dataset.variables.keys()
Out[67]: odict_keys(['station', 'time', 'lat', 'lon', 'DDVEC', 'FHVEC', 'TG', 'RH', 'UG', 'EV24', 'PG', 'iso_dataset', 'product', 'projection'])
dataset.dimensions.keys()
Out[68]: odict_keys(['station', 'time'])
我想提取特定“站”的数据并将其放入pandas DataFrame中以执行一些计算。
我尝试过类似的方法来提取数据,但是我知道这不是netCDF文件的工作方式,但是我不知道怎么做。
df = dataset['344',:,:,:,:,:,:,:,:,0,0,0,0,0]
总结性问题:是否可以提取特定站点的数据并将其放入pandas DataFrame中?
import pandas as pd
import xarray as xr
# Open netCDF file and convert to dataframe
open_netcdf = xr.open_dataset(filename)
dataset = open_netcdf.to_dataframe()
# Select data from a tuple index based on station number: 391
df = dataset.iloc[dataset.index.get_level_values(0) == '391',:]