qt/c++ 使用NetCDF库进行nc文件的创建和读取

qt/c++ 使用NetCDF库进行nc文件的创建和读取

前言

记录一下关于nc文件应用案例

一、nc是什么？

NC文件全称是Network Common Data Format，也叫做NetCDF文件，即网络通用数据格式，这种文件格式一开始是专门用于气象学数据的存储，现在已经发展演变为很多数据采集软件的文件生成格式了。

二、使用步骤

1.准备材料文件

NC_includes 文件夹
相关dll库
下载链接

2.导入库

将netcdf.lib导入外部库
导入完成后pro文件如下

win32: LIBS += -L$$PWD/./ -lnetcdf
INCLUDEPATH += $$PWD/.
DEPENDPATH += $$PWD/.

---

3.添加现有文件

将NC_includes 放入项目目录，右击项目选择添加现有文件将NC_includes 文件夹下文件加入项目

4.编写代码

4.1创建nc

头文件添加

#include <NC_include/netcdfcpp.h>
using namespace netCDF;
using namespace netCDF::exceptions;

cpp文件

		NcFile nc_file("ocean.nc", NcFile::replace); // 创建或打开 NC 文件
        if (!nc_file.isNull()) { // 检查文件是否成功创建或打开
            // 定义和写入相关变量信息（坐标、时间等）
            NcDim x_dim = nc_file.addDim("x", 100);
            NcDim y_dim = nc_file.addDim("y", 100);
            NcDim z_dim = nc_file.addDim("z", 50);
            NcDim time_dim = nc_file.addDim("time");
            NcVar x_var = nc_file.addVar("x", ncDouble, x_dim);
            NcVar y_var = nc_file.addVar("y", ncDouble, y_dim);
            NcVar z_var = nc_file.addVar("z", ncDouble, z_dim);
            NcVar time_var = nc_file.addVar("time", ncInt, time_dim);
            // 写入变量数据（在此以简单的海浪波形为例）
            double x_data[100], y_data[100], z_data[50];
            for (int i = 0; i < 100; i++) {
                x_data[i] = 100 * cos(i * PI / 50);
                y_data[i] = 100 * sin(i * PI / 50);
            }
            for (int i = 0; i < 50; i++) {
                z_data[i] = 10 * sin(i * PI / 50) + 50;
            }
            x_var.putVar(x_data);
            y_var.putVar(y_data);
            z_var.putVar(z_data);
            // 写入时间信息
            int time_val = 0;
            time_var.putVar(&time_val);
            nc_file.close(); // 关闭 NC 文件
            qDebug()<< "生成完成！";

4.2读取nc

NcFile dataFile("ocean.nc", NcFile::read);
    std::multimap<std::string, NcVar> group_var;
    group_var = dataFile.getVars();
    std::multimap<std::string, NcVar>::iterator it_var;
    it_var = group_var.begin();
    while (it_var != group_var.end()) {
            //读取变量ID
            //        int id = it_var->second.getId();
            //        qDebug()<<id;
            //读取变量名称
            std::string name = it_var->second.getName();
            QString txt_name = QString::fromStdString(name);
            qDebug()<<txt_name;
            if(txt_name=="x"){
            NcVar pVar = dataFile.getVar(name);
            double data[100];
            pVar.getVar(data);
            for (int i = 0; i < 100; ++i) {
                qDebug()<<data[i];
            }
            }
            //读取变量类型
        //        std::string type = it_var->second.getType().getName();
        //        QString txt_type = QString::fromStdString(type);
        //        qDebug()<<txt_type;
        //读取变量所拥有的维度名称长度
        //        int dim = it_var->second.getDimCount();
        //        QString txt_dim = QString::number(dim);
        //        qDebug()<<txt_dim;
            it_var++;
    }

打印结果

"time"
"x"
100
99.8027
99.2115
98.2287
96.8583
95.1057
92.9776

总结

以前使用nc库时经常会报一下dll缺失的错误，这次异常顺利，这些dll文件都整理到dll文件夹中了，有报错可以加入试试。
关于x的长度100如何自动获取目前我还没有找到方法，应该是可以自动获取的，多尝试试试。