使用 Python 将 NetCDF 文件转换为 CSV 或文本

Question

我正在尝试使用 Python 将 netCDF 文件转换为 CSV 或文本文件。我已经阅读了this post，但我仍然缺少一步（我是 Python 新手）。它是一个包含纬度、经度、时间和降水数据的数据集。

这是我目前的代码：

import netCDF4
import pandas as pd

precip_nc_file = 'file_path'
nc = netCDF4.Dataset(precip_nc_file, mode='r')

nc.variables.keys()

lat = nc.variables['lat'][:]
lon = nc.variables['lon'][:]
time_var = nc.variables['time']
dtime = netCDF4.num2date(time_var[:],time_var.units)
precip = nc.variables['precip'][:]

我不知道如何从这里开始，但我知道这是用 pandas 创建数据框的问题。

Answer 1

我认为pandas.Series 应该可以为您创建一个包含时间、纬度、经度、沉淀的 CSV。

import netCDF4
import pandas as pd

precip_nc_file = 'file_path'
nc = netCDF4.Dataset(precip_nc_file, mode='r')

nc.variables.keys()

lat = nc.variables['lat'][:]
lon = nc.variables['lon'][:]
time_var = nc.variables['time']
dtime = netCDF4.num2date(time_var[:],time_var.units)
precip = nc.variables['precip'][:]

# a pandas.Series designed for time series of a 2D lat,lon grid
precip_ts = pd.Series(precip, index=dtime) 

precip_ts.to_csv('precip.csv',index=True, header=True)

Answer 2

import xarray as xr

nc = xr.open_dataset('file_path')
nc.precip.to_dataframe().to_csv('precip.csv')

Answer 3

根据你的要求，你或许可以使用 Numpy 的savetxt 方法：

import numpy as np

np.savetxt('lat.csv', lat, delimiter=',')
np.savetxt('lon.csv', lon, delimiter=',')
np.savetxt('precip.csv', precip, delimiter=',')

但是，这将输出没有任何标题或索引列的数据。

如果您确实需要这些功能，您可以构建一个 DataFrame 并将其保存为 CSV，如下所示：

df_lat = pd.DataFrame(data=lat, index=dtime)
df_lat.to_csv('lat.csv')

# and the same for `lon` and `precip`.

注意：在这里，我假设日期/时间索引沿数据的第一维运行。