netcdf数据nc格式数据java解析总结(二)使用java解析netcdf数据

素颜马尾好姑娘i 2022-07-16 06:24 783阅读 2赞

这文章居然有不少同学看啊，官方入门文档在 http://www.unidata.ucar.edu/software/thredds/current/netcdf-java/tutorial/NetcdfFile.html     [点击打开链接][Link 1]

大家多琢磨多调试，有问题加断点debug一步步跟就好了啊，尤其是上面链接里面加粗字的部分，都是关键代码啊。

官方给出的一些示例代码 http://www.unidata.ucar.edu/software/netcdf/examples/programs/

占位，第三篇准备写怎么写怎么生成nc格式的文件。最近跟这货打交道太多了，然后中文资料也少，自己总结了点。

有个包好像叫netcdf.jar，50多兆，版4.几都行，maven也有，百度谷歌一下，有官网，官网的用户手册http://www.unidata.ucar.edu/software/thredds/current/netcdjava/tutorial/NetcdfFile.html大家认真看看基本就能看懂了，当然看不懂英语除外。

maven的netcdf我用的这个

<dependency>
    			<groupId>edu.ucar</groupId>
    			<artifactId>netcdf4</artifactId>
    			<version>4.5.5</version>
    		</dependency>

我自己呢是用到了这些类

import ucar.ma2.Array;  
import ucar.ma2.ArrayFloat;  
import ucar.nc2.NetcdfFile;  
import ucar.nc2.Variable;  
import ucar.nc2.dataset.NetcdfDataset;

所有编程问题其实都是“要把大象装冰箱，总共分几步”的问题。

1、先把冰箱门打开（打开Netcdf文件）

NetcdfFile ncfile = NetcdfDataset.open(filename)

这里会提示让你trycatch  顺带把下面的都try'好了

filename是文件名，写绝对路径

2、把大象装进去（下面几个步骤名是我瞎编的）

第一小步，找大象？

Variable v = ncfile.findVariable(var)

var是变量名，你想读什么变量写什么，想知道文件有哪些变量看上一篇，或者debug移到你读的那个ncfile 变量上看，variables下面就是各种变量，而且也能看出来维度。

![Center][]

第二小步，看看大象有多大？

int\[\] varShape = v.getShape()

这句是获得维度，就是看v变量是一维二维还是3维呀，几个数就是几维，比如上面代码返回的这个int数组是12,200,200，说明v这个变量就是三维的，看用户入门指南写的都很清楚。

。

第三小步，装大象，不太想写了

读全部的数据用v.read()就行，再后面接个copyTo1DJavaArray()这就转成java一维数据类型了（前提你read出来的是二维，三维要先减一维），就基本类型的一维数组呀；

copyToNDJavaArray() 这个变量是几维就会转几维，返回的东西找个变量比如float\[\]\[\] ，float\[\]\[\]\[\]存就好啦，二维三维数组各取所需吧。read返回的好像是object还是神马，总之很好转成java。

当然如果数据太读部分也是可以读部分的，可以设置origin起始点，设置大小 size，这时读数据就要用 v.read(origin, size)，这俩参数是数组，鼠标移到方法上面看就好啦。

v.read().reduce(0) 是去掉一个维度 ，我一般三维转二维会用到这个

挺多方法的，入门上没写，大家自己在read()后面 “.” 一下看提示也能懂个7788。

3、把冰箱门带上

finally里一定要写ncfile.close();把文件关了。

贴个简单代码吧 ，大家可以参考一下

假设有个nc文件，三维的，分别是time，lon，lat

比如预报未来6个小时次的天气，每次都是预报lon\*lat个点，假设是200\*200，经纬度是2维数据，每个对应位置的经度和纬度组成一个点，天气是3维，就是第一个小时对应一个二维200\*200点的数，第2个小时对应一个二维200\*200点的数

要是还没说明白 ，举个例子吧，三维的等会儿再写。。。不知道为什么贴的代码这么混乱啊？

public void aa() {
    		// read完copyToNDJavaArray，这个读出来是几维转的就是几维，直接转相应几维java数组就好。
    
    		try {
    			NetcdfFile ncfile = null;
    			ncfile = NetcdfDataset.open("D://XXXXX.nc");
    			// 存经纬度
    			String var1 = "lon";// 此处严格区分大小写，不然找不到，不知道有什么变量的可以断点debug一下，鼠标移到上面ncfile那行看
    			String var2 = "lat";
    			Variable v1 = ncfile.findVariable(var1);
    			Variable v2 = ncfile.findVariable(var2);
    
    			float[][] lon = (float[][]) v1.read().copyToNDJavaArray();// 因为经纬度是二维的，直接
    			// copyToNDJavaArray结果就是二维的，强转一下就好，然后java的二维数组大家该咋用咋用就行了
    			float[][] lat = (float[][]) v2.read().copyToNDJavaArray();
    		} catch (IOException e1) {// TODO Auto-generated catch block
    			e1.printStackTrace();
    		}
    	}

[Link 1]: http://www.unidata.ucar.edu/software/thredds/current/netcdf-java/tutorial/NetcdfFile.html
[Center]: /images/20220716/7e20ded6955c427f8ac4ea380574d1b0.png