资源描述:
《grads绘图系统在站点等值线中的应用》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、GrADS绘图系统在站点等值线中的应用第2001姜期JOURNAlO广FG西UAN气GX象IMETEOROLOGY年6月V01.22No.2Jun.2001文章编号:1001—5191(2001)02-0053-04GRADS绘图系统在站点等值线中的应用况雪源(广西气候中心,广西南宁530021)摘要;姒一个实例详细介绍了利用GRADS绘图系统绘制具有不规则边界的站点资料等值线的方法,使站点等值线的绘制变得简单易行,图形输出更为美观.并对GRADS系统存在的一些问豚进行了讨论.荚镑调;GRADS系统;站点等值线绘图中囝分
2、类号:TP317,4文献标识码BApplicationofGRADSDrawingSysteminStationIsolineKUNXL!e—yuan(GuangxiClimateCenter,Nanning530022,China)Abstract:Bytakinganactualexample,thispaperintroducesindetailtheapplicationofGRADSdrawingsystemtomappinginformationisolineforthestationswithirregul
3、arboundary.Thismethodmakesthemappingofstationisolineeasyandsimple,andmorehandsomeimageoutput.ItalsodiscussestheexistingproblemsofGRADSdrawingsystem.Keywards:GRADSsystem;stationisoline;mapping引言1基本思路GrADS(GridAnalysisandDisplaySystem)绘图系统是一组多平台应用及气象数据图形分析和显示系统,该系统
4、具有气象数据分析功能强,地图投影丰富,高级编程语言容易,图形显示快速,并具有彩色动画功能等特点,并已成为国际气象数据显示的标准平台之一.GRADS系统提供了一个交互操作的气象格点数据和站点数据的分析和显示环境,但相对于格点资料,GRADS处理站点资料的能力较弱,且系统无法直接采用站点资料来绘制等值线.对于一般用户来说,如何利用GRADS绘制格点等值线功能及输出图形美观等优点,将一定区域的站点资料以等值线方式显示,并使等值线不越过不规则边界呢?刘宣飞Lj对此作了一些介绍,但较为简略.本文将针对这些问题作进一步探讨,并以绘制
5、广西88个站点2000年6月平均气温分布图为例作详细介绍.文中程序部分用Basic语言编写,对GRADS系统的基本使用可参看文献].GRADS系统无法采用站点资料来绘制等值线,而需将站点资料内插成格点资料后才能绘制等值线,GRADS提供的oacres内部函数可实现此内插功能,同时,内部函数maskout可将符合一定条件的格点值设为缺测,为此,我们采取以下方式进行站点等值线的绘制.首先,形成一份格点资料文件,此文件提供插值格点区域,并将边界外资料设为负值,边界内资料设为正值,然后,采用oacres函数将站点资料内插成格点资
6、料,并利用maskout函数将边界外的插值场设为缺测后绘制等值线,最后,利用广西区界,地区界经纬度资料在图上定位以绘制边界.2GRADS绘制站点等值线步骤2.1所需基本资料及格式为绘制广西88个站点2000年6月平均气温分布图形,需以下资料(表1);收稿日期:2000—07—31作者简升;况雪源(1971--)t女.贵州赫章人,硬士.广西气候中心预测工程师,现从事短期气候珂工作54广酉气象22卷表l绘制等值线圈资料2.2站点数据文件及其描述文件2.2.1站点文件格式站点资料采用二进制文件进行存储,按每个时次一个报告的顺序
7、直接记录各站资料,每个时次的站点报告组叉分成两部分:地面报和高空报.站点报数据的排放顺序为:①提供站点经纬度的头记录②地面报变量③高空报变量④一个时次完成后加上一个特殊头记录(没有数据组),表示本时次数据报结束.头记录共占28个字节,包括:STAT:站点标识,字符型(string或char),占8个字节;其赋值可以是任意数字,也可以是字母.RLON:经度,单精度型(single),占4个字节.RLAT:纬度,单精度型(single),占4个字节.TIM:本报的时次,单精度型(single),占4个字节,NLEV:头记录后
8、的数据组数,单精度型(single),占4个字节,包括一个地面报组加上各高度层数,取0时标志一个时次报文的结束.NFLAG:单精度型(single),占4个字节.取0表示头记录的报文中没有地面变量,取1时则有.在头记录之后,顺序记录本时次的各报文,首先写入地面变量,然后是高空变量.当本时次的所有报告写完后,写出一个特