关于nc文件的读取.doc

关于nc文件的读取2012-11-02 21:47:01|分类：信号处理DSP|标签：nc文件panoply|举报|字号订阅最近在学习冰后回弹模型（http:/www.atmosp.physics.utoronto.ca/peltier/data.php），里面给出的数据为nc（netcdf）格式的，读取是个问题，在言深深的帮助下，成功读取，在这里将其总结出来，以备后用也希望可以方便有用之人。有几种方式可以进行读取：1、Excel的一个插件，NetCDF4Excel可以完成相关的查看工作，需要的可以进行下载。安装之后，我们可以得到，一个快捷方式，如图：打开之后，最上面标题栏点击加载项，我们可以看到下面的情形： 然后打开nc文件，我们就可以进行查看操作：2、通过Grads进行相关的读取，主要是在言深深同学的帮助下完成，深深的感谢！第一步，安装相关的grads软件，这个这里不做介绍，相关可以参考气象论坛网站；第二步，打开grads软件，如图，回车之后的那个窗口千万不能关，那个窗口是图像显示的，关闭之后输入命令之后直接就退出了，刚开始我老犯这个低级错误（对菜鸟而言）。 第三步，可以通过两个命令简单的看到数据的头文件包含信息：sdfopen d:1.ncq ctlinfo1.nc文件然后编写gs文件，如下：reinitsdfopen D:1.ncset gxout fwriteset fwrite D:1.datset x 1 360set y 1 180set z 1set t 1 1d dsea_250disable fwriteps1：路径自己可以按需更改，不要有空格。ps2：这里的dsea_250可不是随便取的，看上面的ctlinfo心里里最后一行，必须要一致！然后就可以生成一个1.dat的文件。第四步，通过fortran程序将1.dat读为txt文档 Program nc2txt parameter(m=180,n=360,o=1) dimension x(m,n,o) open(1,file=1.dat,access=direct,form=binary,recl=n*m*o*4) read(1,rec=1) (x(i,j,k),j=1,n),i=1,m),k=1,o) close(1) open(2,file=Ice-5gGeoid.txt) do i=1,m !print*,i do j=1,n write(2,(2f10.1,2x,f15.6) -(i-1)*1.+180,(j)*1.,x(i,j,1) enddo enddo close(2) End program nc2txt ps：这里的参数要与自己的数据一致，比如mno，以及最后的输出，有的是1，有的是2.5。然后我们就可以得到相关的数据文件如下：至此，数据读取操作完成，需要对有需要的人有帮助，再次对言深深表示感谢。此外grads还可以通过那个窗口绘图，具体如下：（1）sdfopen d:1.nc (2) d dsea_250Nc2TxT（版权归作者所有，本着学习的精神）3、一个国外学者对我的回复：In the meantime, may I suggest that you investigate the following free software:ncview and ncdump (for Linux)panoply(for Windows)/tools/panoply/Personally,I use NCAR Graphics NCL programs to create plots from netcdf files. That is also free software available from NCAR (National Center for Atmospheric Research in Boulder, Colorado, USA) for a variety of platforms (Linux and Windows)。/Download/install.shtmlncview and panoply are both viewers that display information from a netcdf file directly. This would enable you to check any of your own programs that read the files.ncdump will actually display the values for any/all variables in the netcdf file which is also a way to check but not easy as there are a lot of numbers in the file Rosemarie Drummond这才是真正的治学精神，话说这个人真的是个大牛，都给我回复的这么详细基础！panoply很强悍，打开之后，可以直接plot出图：4、通过Fortran程序来读取，这个算是高