关于nc文件的读取分析.docxVIP

关于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.nc?②q 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_250?Nc2TxT（版权归作者所有，本着学习的精神）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.shtml?ncview 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