使用java创建netcdf文件，并写入数据

2019-11-06 08:39:06

字体：大中小

来源：转载

供稿：网友

学气象的少不了使用netcdf文件，这文件是干嘛的我就不多提了，网上多得是。网上用java读取netcdf文件的帖子多如牛毛，我今天要做的就是如何使用java来写入nc文件废话不多少，先上干活。首先你得引入一个包。我这里是引入了

netcdfAll-4.6.6.jar

以下为具体代码

package com.cym.utils;import java.io.IOException;import java.util.ArrayList;import ucar.ma2.ArrayFloat;import ucar.ma2.DataType;import ucar.ma2.Index;import ucar.ma2.InvalidRangeException;import ucar.nc2.Attribute;import ucar.nc2.Dimension;import ucar.nc2.NetcdfFileWriter;import ucar.nc2.Variable;import ucar.nc2.NetcdfFileWriter.Version;public class NCfile { public static int width=20,height=16; public static double lonWest=110,latSouth=30; public static double gridInterval=0.5; /** * * @param filepath 创建的文件路径 * @param width 网格横向长度 * @param height 网格纵向高度 * @param lonWest 网格左下角起始经度 * @param latSouth 网格左下角起始维度 * @param gridInterval 网格间距 * @param d d[i][j]为网格i和j处对应的数据 * @throws IOException * @throws InvalidRangeException */ public static void createNCfile(String filepath,int width,int height,double lonWest,double latSouth,double gridInterval,double d[][]){ //创建netcdf3写入文件对象 NetcdfFileWriter fileWriter; try { fileWriter = NetcdfFileWriter.createNew(Version.netcdf3, filepath); //分别创建经度和纬度方向的维度，lon的宽度为width,lat高度为height,并命名为lon和lat Dimension xdim=fileWriter.addDimension(null, "lon", width); Dimension ydim=fileWriter.addDimension(null, "lat", height); //通过一维组合成二维array list ArrayList<Dimension> dims=new ArrayList<Dimension>(); //注意这里add的顺序很重要,一定要和数据对应，如果把list顺序调换，那么二维的维度也会不同，相当于数学中的坐标系，第一个是x坐标，第二个是y坐标，x宽度为width,y高度为height dims.add(xdim); dims.add(ydim); //创建名称为lon的变量,类型为folat，对应的维度为lon,对应Dimension里面定义的名称为"lon"的那个 Variable vx=fileWriter.addVariable(null, "lon", DataType.FLOAT, "lon"); //给vx变量加上数据说明，全程long_name fileWriter.addVariableAttribute(vx, new Attribute("long_name","longitude" )); //给vx变量加上units，对应的Degree_east，如果不加units为Degree_east(东经)，则不是地理坐标系，这个很重要 fileWriter.addVariableAttribute(vx, new Attribute("units","Degrees_east" )); Variable vy=fileWriter.addVariable(null, "lat", DataType.FLOAT, "lat"); fileWriter.addVariableAttribute(vy, new Attribute("long_name","latitude" )); //给vy变量加上units，对应的Degree_north，如果不加units为Degree_north(北纬)，则不是地理坐标系，这个很重要 fileWriter.addVariableAttribute(vy, new Attribute("units","Degrees_north" )); //创建变量名称为var的变量，对应的维度为dims,该dims为上面定义的一个ArrayList<Diminsion>，该list包含2个维度，经度和维度 Variable v=fileWriter.addVariable(null, "var", DataType.FLOAT, dims); fileWriter.addVariableAttribute(v, new Attribute("long_name","数据的详细描述" )); fileWriter.addVariableAttribute(v, new Attribute("units","unit" )); //使用create()来创建该文件，只有create之后才能在文件夹中显示该文件，并且写入数据 fileWriter.create(); //创建x和y方向上的变量。D1 代表一维 ArrayFloat xvalues=new ArrayFloat.D1(width); ArrayFloat yvalues=new ArrayFloat.D1(height); //创建网格上的变量，D2代表二维，用(width, height)来描述他的定义域，这个顺序必须和ArrayList<Dimension> 加载的顺序一一对应，否则会出错 ArrayFloat values=new ArrayFloat.D2(width, height); Index index=values.getIndex(); for (int i = 0; i < width; i++) { for (int j = 0; j < height; j++) { //将二维数组[i][j]处的数据写入网格的i,j处 values.setDouble(index.set(i, j), d[i][j]); } } for (int i = 0; i < width; i++) { xvalues.setFloat(i, (float) (lonWest+gridInterval*i)); } for (int j = 0; j < height; j++) { yvalues.setFloat(j, (float) (latSouth+gridInterval*j)); } fileWriter.write(vx, xvalues); fileWriter.write(vy, yvalues); fileWriter.write(v, values); fileWriter.close(); System.out.PRintln("文件已存入"+filepath); } catch (IOException e) { // TODO Auto-generated catch block e.printStackTrace(); System.out.println("文件创建失败，请检查路径"); } catch (InvalidRangeException e) { // TODO Auto-generated catch block e.printStackTrace(); System.out.println("文件写入错误，超出范围"); } } public static void main(String args[]){ double d[][]=new double[width][height]; for (int i = 0; i < width; i++) { for (int j = 0; j < height; j++) { d[i][j]=Math.sin(i+j); } } //放到根目录会创建失败 String filepath="d://output//testNCfile.nc"; NCfile.createNCfile(filepath, width, height, lonWest, latSouth, gridInterval, d); }}

运行后，控制台会打印

这里写图片描述

使用Panoply来打开看看该文件这里写图片描述可以看到有三个变量，分别为lon,lat,var，其中Type分别为1D,1D,Geo2D，如果你没有在代码中加入如下代码，则就是普通的2D坐标，不是地理2D坐标，当时我用没有加单位的文件用ncl作图就会出现问题

fileWriter.addVariableAttribute(vx, new Attribute("units","Degrees_east" ));fileWriter.addVariableAttribute(vy, new Attribute("units","Degrees_north"));

这里都是我刚写入的数据，数据都是瞎编的，实际工作中按照需求你们自己去修改吧，这是瞎编的数据

d[i][j]=Math.sin(i+j);

使用panoply查看数据这里写图片描述