Netcdf文件导出基本代码示例

需求:从一个nc文件中取出我想要的变量,导出到另外一个nc文件。

下面是基本代码:

@GetMapping("/export")
public String varExport(@ModelAttribute FileSaveVars params ) throws IOException, InvalidRangeException {
//      1° 设置原始文件路径+转存文件路径
        String filePath = "/Users/caowei/workspace/jszone/大气系统/ncdemo/src/main/resources/static/MERRA2_400.tavg1_2d_flx_Nx.20180101.SUB.nc";
        String savePath = "/Users/caowei/workspace/test.nc";
//      2° 读取文件
        NetcdfFile ncFile = null;
        ncFile = NetcdfFile.open(filePath, null);
//      3° 获取需要转存的变量
        Variable latVar = ncFile.findVariable("lat");
//      4° 读取变量的一些信息
//      4°-1 长度
        int latitudeCount = (int) latVar.getSize();
//      4°-2 读取变量数组数据
        ArrayDouble.D1 latArray = (ArrayDouble.D1) latVar.read();
//      5° 构造输出保存文件
        NetcdfFileWriter ncWrite = NetcdfFileWriter.createNew(NetcdfFileWriter.Version.netcdf4, savePath);
//      6° 设置文件的一些参数
//      6°-1纬度
        Dimension latitudeDimension = ncWrite.addDimension(null, "latitude", (int) latArray.getSize());
//      6°-2变量
        Variable latitudeVar = ncWrite.addVariable(null, "latitude", DataType.DOUBLE, "latitude");
//      6°-3设置属性
        latitudeVar.addAttribute(new Attribute("standard_name", "latitude"));
        latitudeVar.addAttribute(new Attribute("long_name", "latitude"));
        latitudeVar.addAttribute(new Attribute("units", "degrees_north"));
        latitudeVar.addAttribute(new Attribute("axis", "Y"));
//      7°创建
        ncWrite.create();
//      8° 构造数据
        ArrayDouble.D1 latitudeData = latArray;
//      9°将数据写入
        ncWrite.write(latitudeVar, latitudeData);
//      10°关闭保存
        ncWrite.close();
        return  "OK";
}

下面是抽出一个变量进行导出

  @GetMapping("/export")
    public String varExport(@ModelAttribute FileSaveVars params ) throws IOException, InvalidRangeException {
//      1° 设置原始文件路径+转存文件路径
        String filePath = "/Users/caowei/workspace/jszone/大气系统/ncdemo/src/main/resources/static/MERRA2_400.tavg1_2d_flx_Nx.20180101.SUB.nc";
        String savePath = "/Users/caowei/workspace/test.nc";
//      2° 读取文件
        NetcdfFile ncFile = null;
        ncFile = NetcdfFile.open(filePath, null);
//      3° 设置需要转存的变量
        String varName = "lon";
//      4° 获取需要转存的变量对象
        Variable varObject = ncFile.findVariable(varName);
//      5° 读取变量的一些信息
//      5°-1 长度
        int varCount = (int) varObject.getSize();
//      6° 构造输出保存文件
        NetcdfFileWriter ncWrite = NetcdfFileWriter.createNew(NetcdfFileWriter.Version.netcdf4, savePath);
//      7° 设置文件的一些参数
//      7°-1纬度
        Dimension varDimension = ncWrite.addDimension(null, varObject.getName().toString(), varCount);
//      7°-2变量
        Variable varVar = ncWrite.addVariable(null, varObject.getName().toString(), varObject.getDataType(), varDimension.getName());
//      7°-3设置属性
        List<Attribute> AttrList = varObject.getAttributes();
        for (Attribute att : AttrList) {
            varVar.addAttribute(att);
        }
//      8°创建
        ncWrite.create();
//      9°将数据写入
        ncWrite.write(varVar, varObject.read());
//      10°关闭保存
        ncWrite.close();
        return  "OK";
}

下面再次进行升级,可以直接导出单个三纬变量

@GetMapping("/export")
public String varExport(@ModelAttribute FileSaveVars params ) throws IOException, InvalidRangeException {
//      1° 设置原始文件路径+转存文件路径
        String filePath = "/Users/caowei/workspace/jszone/大气系统/ncdemo/src/main/resources/static/MERRA2_400.tavg1_2d_flx_Nx.20180101.SUB.nc";
        String savePath = "/Users/caowei/workspace/test.nc";
//      2° 读取文件
        NetcdfFile ncFile = NetcdfFile.open(filePath, null);
//      3° 设置需要转存的变量
        String varName = "QLML";
//      4° 获取需要转存的变量对象
        Variable varObject = ncFile.findVariable(varName);
//      5° 构造输出保存文件
        NetcdfFileWriter ncWrite = NetcdfFileWriter.createNew(NetcdfFileWriter.Version.netcdf4, savePath);
//      6° 设置文件的一些参数
//      7°-1纬度
        for(Dimension dimension : varObject.getDimensions()){
            ncWrite.addDimension(null, dimension.getShortName(), dimension.getLength());
        }
//      6°-2变量
        Variable varVar = ncWrite.addVariable(null, varObject.getShortName(), varObject.getDataType(), varObject.getDimensionsString());
//      6°-3设置属性
        List<Attribute> AttrList = varObject.getAttributes();
        for (Attribute att : AttrList) {
            varVar.addAttribute(att);
        }
//      7°创建
        ncWrite.create();
//      8°将数据写入
        ncWrite.write(varVar, varObject.read());
//      9°关闭保存
        ncWrite.close();
        return  "OK";
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值