[用NCL处理卫星数据一个简单的教程.docxVIP

用NCL处理卫星数据(一个简单的教程)2010年03月19日星期五 15:13这个教程主要介绍如何用NCL读取卫星数据，如何将数据存储为二进制格式，如何作图。本教程所用的数据：/file/f19ba404fc这个数据是AMSR-E的L2A产品，更多信息请查阅相关文档。在开始这个教程前，请确认已经安装好了NCL。1. 运行NCL和GrADS类似，NCL也可以运行在交互式环境中，即输入一条命令执行一条命令。交互式环境使用：在命令行下输入$ ncl退出请输入命令ncl 0 exit交互式环境不方便重复操作，一旦退出程序就需要全部重新来过，所以推荐编写NCL脚本。一个NCL脚本主要包括以下几个部分：load? 外部函数begin?;命令和操作endload用来载入外部函数。文件操作、数据操作、数据可视化等代码都放在load 语句之后。半角分号“;”后的内容是注释，NCL只有行注释功能。begin 和 end 可有可无，但如果出现了begin，就必须有end 与之对应。NCL脚本扩展名为.ncl。比如有一NCL脚本文件为 script.ncl，运行时可以用下面的方法：$ ncl script.ncl或 $ ncl script.ncl下面来写一个简单的NCL脚本。打开vi之类的编辑软件，将下面的内容复制进去，然后保存到卫星数据文件所在目录，文件名为test.ncl;********************************************************;? Load NCL scriptsload $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.nclload $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.nclload $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.nclload $NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.nclload $NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.nclload $NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl; load语句载入了NCL常用的几个外部程序，建议以后自己写脚本时也保留这6行的内容。;********************************************************? begin? fi = addfile(./AMSR_E_L2A.hdf,r)? print(fi)? end运行这个脚本$ ncl test.ncl log2. 熟悉数据打开log文件从第11行相关信息可以知道该文件为 HDFEOS_V2.5第663行开始可以得知该卫星数据的生成时间等信息第957行开始是卫星数据的维度信息(dimensions)第983行起是变量信息(variables)这些对文件内容进行描述的信息被称作元数据(metadata)，元数据对于数据的存储和发布很重要。它记录了数据集的主要特征，包括数据结构和内容描述，以及数据集开发的历史等。3. 读取数据下面以89V的数据为例89V的元数据为：short?89_0V_Res_1_TB__4?(?DataTrack_lo_Low_Res_Swath, DataXTrack_lo_Low_Res_Swath?)SCALE_FACTOR :??? 0.01OFFSET :??? 327.68UNIT :??? Kelvinhdf_name :??? 89.0V_Res.1_TBhdf_group :??? Low_Res_Swath/Data Fieldshdf_group_id :??? 4short：该数据是short格式89_0V_Res_1_TB__4：变量名DataTrack_lo_Low_Res_SwathDataXTrack_lo_Low_Res_Swath：数据大小，可以从维度信息(dimensions)了解数据大小。SCALE_FACTOROFFSET：为了减小存储数据所占的空间，一般都会对实型数据进行加减乘除处理，以整型来存储(pack)，用户在使用时，需要将整型转换为实型(unpack)，会用到scale_factor和offset信息。修改begin 和 end 之间的内容：? begin? fi = addfile(./AMSR_E_L2A.hdf,r)? data =new(dimsizes(fi-89_0V_Res_1_TB__4),float,-9999)? data = (fi-8