将NetCDF4变量q转为以x/y为行列的Pandas DataFrame并导出CSV

核心步骤 #

你的变量q维度为(y, x)，要实现行关联x、列关联y，需先转置数组，再转为DataFrame后保存：

1. 提取并转置q数据

# 取出q的数值并转置，将维度从(y, x)转为(x, y)
q_array = ds['q'].values.T
# 转为DataFrame
df_q = pd.DataFrame(q_array)

2. （可选）添加有意义的行/列标识
   如果需要用维度索引或相关变量作为行列名：

# 设置行名为x维度的索引值
df_q.index = ds.x.values
# 设置列名为y对应的流域ID（wshed_id），转成整数更易读
df_q.columns = ds['wshed_id'].values.astype(int)

3. 保存为CSV
   由于数据集规模极大，直接写入可能内存不足，建议分块保存：

# 分块写入CSV，每块1000行，避免内存溢出
df_q.to_csv('goa_q_data.csv', chunksize=1000)

测试建议 #

原数据集维度（x=12784，y=14052）过大，先取小样本验证流程：

# 截取前100个y和前100个x的样本
q_sample = ds['q'].isel(y=slice(0, 100), x=slice(0, 100)).values.T
df_sample = pd.DataFrame(q_sample)
df_sample.to_csv('q_sample.csv')

内容的提问来源于stack exchange，提问作者Aeon