2. 中国海洋大学,青岛 266101;
3. 北京大学物理学院大气科学系, 北京 100081
2. Ocean University of China, Qingdao 266101;
3. Department of Atmospheric Sciences, Physics School, Peking University, Beijing 100871
气象数据分析和图形软件是现代气象业务、科研工作的必备工具。目前国内有较多气象数据和图形处理软件[1-5],其中气象信息综合分析处理系统(Meteoroloical Information Comprehensive Analysis and Processing System, MICAPS)是我国使用最为广泛的天气预报业务气象图形软件系统,目前最新版本是3.1[6]。MICAPS系统功能强大,数据兼容性好,使用方便简捷,是主要的天气预报业务工作平台。在气象科研工作上经常使用的气象绘图软件还有GRADS和Surfer。GRADS是一款免费跨平台气象绘图软件,获得了较为广泛的应用,目前它的最新版本为2.1,该软件主要使用命令行和脚本编程方式工作。美国Golden Software公司开发的Surfer软件是一款应用极为广泛的3D地质绘图软件,目前最新版本为9.0。Surfer操作界面友好,但不是专门的气象绘图软件,因此许多方面不能满足气象绘图的需要,例如Surfer本身不能生成地图、不能绘制流线图等。
北京大学大气科学系曾在Windows 32操作系统下开发了操作界面友好的客观分析绘图软件——客观分析诊断图形系统[1, 4],该软件主要面向天气分析科研工作,与天气分析业务软件系统MICAPS功能互补。该软件不仅能够对常规气象地面、探空数据进行客观分析诊断,且能够处理多种格式的气象数据和方便简捷地绘制各种常用气象图形,获得了较为广泛的应用。但随着大气科学探测技术的迅猛发展,新的气象探测资料在天气预报业务中广泛应用,特别是为我国新一代静止气象卫星和新一代天气雷达的业务运行提供了大量遥测气象数据,而旧版本的客观分析诊断图形系统已经不能适应目前气象业务和科研的需求,因此,在该背景下,我们对客观分析诊断图形系统软件做了进一步的改进与完善。
软件的改进主要包括增加新数据接口和新算法、改进算法、简化操作等。增加新数据接口和新算法是为了进一步满足气象业务和科研的需求;改进算法是为了提高算法的可靠性、效率与显示效果;简化操作是为了进一步提高用户工作效率。
1 软件新增数据接口新软件系统除了能够支持常规气象地面和探空观测原始数据、MICAPS格式数据、美国国家环境预报中心(National Centers for Enriromental Prediction,NCEP)再分析格点数据、VIS5D格式格点数据等数据外,也可以处理二进制格点加工数据(GRIdded Binary, GRIB)的GRIB2格式和GRADS格式的NCEP数据,FY2C、FY2D和MTSAT(Multi-functional Transport Satellite)地球静止卫星和雷达基数据,同时也能处理这些数据的压缩格式数据。
新软件能够读取的地球静止卫星数据包括国家卫星气象中心下发AWX(AWX是指由国家卫星气象中心所生成的卫星产品)格式的多种类型静止卫星云图资料与辐射亮温(blak body temperature, TBB)产品,国家卫星气象中心的HDF格式标称圆盘图像资料,星地通公司的GPF格式资料,北京大学大气科学系的GMS-5、GOES-9、MTSAT-1R格式资料,日本KoChi大学的GMS-5、GOES-9、MTSAT资料。读取的云图资料根据通道的不同既可以绘制云图(见图 1左下),也可以绘制TBB或者反照率等值线。
新软件可以读取敏视达公司SA/SB雷达基数据和3830 CINRAD C波段雷达基数据,并可以根据读取的基数据形成多个雷达产品和三维直角网格数据,也可以直接选择直线绘制垂直剖面图形(见图 1中下)。
新软件不仅可以处理以上列举的多种未压缩格式的数据,也可以读取以上数据的压缩格式数据。可处理的压缩格式主要有ZIP、RAR、7Z、BZ2、ARJ、ARC、CAB、GZIP、Z、TAR等格式。
不同数据接口绘制的具有相同地理坐标的图形不仅可以相互叠加显示,而且可便捷地编辑图形属性和调整显示次序。
2 软件新增算法为了适应天气分析科研工作的需要,新软件增加了多种算法,既包括客观分析诊断算法、雷达产品算法、强对流天气监测和临近预报算法,也包括新图形绘制和图形生成算法。
新软件中增加了等积Lambert地图投影算法、T-logP图中对流不稳定能量(CAPE和CIN)和各种指数(K指数、SI指数等)计算、等熵面物理量(300~350 K等熵面位势涡度、气压和风)客观分析、探空数据的时间剖面和垂直剖面制作等新算法。
新软件不仅能够读取雷达基数据,还能够计算组合反射率因子、垂直积分液态水含量(vertically integrated liquid, VIL)、回波顶高(echo tops, ET)等雷达数据产品,也能够利用雷达基数据进行交叉相关回波移动计算(CTREC)[7]矢量计算,雷暴识别、追踪、分析和临近预报(TITAN)[8]雷暴追踪和临近预报;并能够利用MICAPS格式数据和气象探测中心提供的闪电数据对冰雹、大风、闪电和短时强降水等强对流天气进行监测和统计。
CTREC算法[7]将反射率因子场分成若干个大小相当的“区域”。将这些在上一时刻的“区域”分别与下一时刻的各个“区域”做空间交叉相关,以找出此刻与上一个时刻的特定区域相关系数最大的“区域”,从而来确定整个区域中不同网格点处的回波移动矢量。
TITAN[8]使用单体质心算法来识别和追踪一个被看作三维实体的风暴。在风暴的检测上,TITAN使用类似三维聚类技术识别风暴,没有考虑风暴的虚假探测和漏测的问题。为了追踪风暴随时间的演变, TITAN将追踪问题转化为组合最优化问题来计算相邻两个雷达图像内风暴体的匹配关系,并且通过一种特殊的、基于风暴体质心位置的几何方法来识别风暴的分裂、合并。对于风暴体位置的预报,TITAN使用风暴体在不同时刻的质心位置进行外推预报。在一般情况下,TITAN可以取得较好的追踪结果。图 2为TITAN算法对2008年9月23日08:18 UTC热带气旋“黑格比”中对流风暴的识别与追踪结果。新软件也可以对TITAN算法的雷暴追踪结果进行统计分析。
韩雷,郑永光等[9-10]对TITAN算法进行了改进。改进的TITAN算法基于数学形态学的识别方法。该方法首先使用第一级阈值进行单阈值识别;其次,对识别得到风暴体执行基于动态卷积模板的腐蚀操作,以消除虚假合并;然后,使用高一级阈值进行识别,并对识别得到的风暴体进行膨胀操作,当风暴体的边界在膨胀的过程中相互之间接触,或接触到了原来较低阈值识别的风暴体的边界时,则停止膨胀过程;最后,逐次使用更高级别的阈值进行识别,并在每一个级阈值的识别过程中执行腐蚀和膨胀操作。试验的结果表明,改进的TITAN算法能够成功地识别出虚假合并风暴体,也能从风暴簇中分离出相距较近的风暴单体。
强对流天气的监测和统计功能不仅可以在天气业务中实现强对流天气的实时监测,也可以使用历史数据进行各类强对流天气时空分布统计[11]。目前该功能已经应用到国家气象中心强天气预报中心的实时业务中,取得了较好的业务效果。该功能输入和输出数据都为MICAPS格式数据文件,可直接与气象业务平台MICAPS衔接;该功能可以监测积云、雷暴、闪电、冰雹、龙卷、大风、短时强降水和深对流云的最近1、3、6、12和24小时的分布,该功能根据资料的观测时次进行时间滑动监测;可以进行月、旬、候等强对流天气分布监测;也可以设置为任意时次、任意天数的强对流实况监测。图 3为本软件监测的2009年6月3—4日大风、冰雹和短时强降水分布。
积云是大气不稳定状态的反映,其中以鬃积雨云发展最为强烈。但积云能否发展为强对流天气还与大气的不稳定程度、触发机制和大尺度环境条件密切相关。监测积云的分布可为强对流天气预报提供大气是否稳定的重要前期信息。积云的监测主要依靠常规地面观测资料完成,包括基本站、基准站和一般站的观测资料。
雷暴是最为普通的强对流天气。本监测系统监测的雷暴,既包括普通雷暴,也包括常规观测的闪电、冰雹、龙卷等强对流天气,是各种强对流天气的总和。雷暴的监测主要依靠常规地面观测资料和重要天气报告完成,包括基本站、基准站和一般站的观测资料。
闪电的监测区分为总地闪、负地闪和正地闪监测,并计算各类闪电的闪电密度。一般来说,较强的对流天气过程其闪电密度也较高。闪电资料的来源主要是中国气象局气象探测中心提供的全国闪电监测网闪电数据。
冰雹是我国暖季较常见的强对流天气,其监测资料来源是常规地面观测资料和重要天气报告,冰雹监测包括冰雹的发生时间和直径大小,监测结果如图 3。龙卷的监测资料来源也是常规地面观测资料和重要天气报告,但从2009年的监测结果来看,几乎没有监测到龙卷。大风的监测资料来源也是常规地面观测资料和重要天气报告,目前没有区分大风的性质,大风监测所取的风速阈值为17 m·s-1。
短时强降水的监测资料主要来源于自动站的1小时降水资料。短时强降水监测等级区分为三类,分别为:≥20 mm·h-1、≥30 mm·h-1、≥50 mm·h-1;系统也对每一个自动站进行最近3、6、12和24小时降水累积。
深对流云监测主要来自我国的风云静止卫星资料。2009年使用的是FY2C资料,2010年以来使用FY2E资料。深对流云的监测资料有两种,一种是使用国家卫星气象中心下发的云分类产品,提取其中的密卷云和积雨云生成深对流云产品;另一种使用红外1和水汽通道亮温资料,设定红外1通道的亮温阈值为-32 ℃,红外1通道和水汽通道的亮温差小于0。
新软件中增加了绘制风矢端图、格点图像、曲线、各种锋面和天气符号等新的图形类型,并在图形合并的基础上增加了图形组合功能。风矢端图可以更直观地了解垂直风切变大小和冷暖平流情况;格点图像直接使用格点数据文件绘制填充图像,不需要提取等值线,适用于格点数较多的数据绘图,比如卫星数据、雷达数据和高分辨率的数值模式数据等;格点图像可以根据需要选择图像透明度。图形组合功能可以把文字等图形对象与气象图形组合成一幅图形处理,方便编辑图形时的选定和移动操作。
新软件中用户可以根据需要选择地理区域裁剪图形(如图 1中右侧上部图形为对中国陆地区域的裁剪),也可以自定义裁剪范围;裁剪操作既可以保留多边形内图形(比如中国陆地区域裁剪),也可以保留多边形外图形(比如青藏高原大地形裁剪等)。
新软件在图形输出操作中增加了一次操作可以生成多幅EMF、GIF或者PNG等格式图像,或者生成矢量格式的可保存多个页面的Postscript图形,也可以利用时间序列数据生成的一系列图像制作动画GIF文件等。
3 软件改进的算法新软件中改进的算法主要有图形绘制、图形生成和离散点生成格点数据算法等。
新软件改进了图形的绘制方式。旧软件中只能使用Windows GDI进行绘图,新软件可使用GDI绘图也可使用GDI+绘图,默认使用GDI+绘图。GDI绘图的优势是绘图速度较快、需要内存少,但绘制的图形没有线条和字体“反走样”,在计算机屏幕上图形显得较为粗糙,不精美;GDI+绘图的优点是绘制的图形中线条和字体可以使用“反走样”功能,图形比较精美细腻,但绘图速度较慢、需要较多系统内存。由于目前计算机的速度和内存都有了很大提高,因此新软件推荐使用GDI+绘图。
新软件对不同类型图形的生成也做了较大改进。对于等值线图形,等值线生成时可以处理缺测数据,在缺测格点处等值线自动中断、填充颜色(默认使用白色显示)。软件读取格点数据时,某个格点数值不小于1×1035,则认为该格点数值缺测。对于生成的等值线图形,用户不仅可以进行等值线修改、等值线标值位置移动、格点数据显示和格点数据修改,也可以在编辑等值线图形时打开新的格点数据文件,提高软件对图形的编辑能力。
新软件对不规则分布离散点填图图形的生成也做了改进。离散点填图图形既可以标注数值,也可以标注文本,还可以绘制风羽或者箭头矢量,这些改进极大地提高了软件的功能和灵活性。新软件中用户可以根据需要对数据文件进行编辑,用户还可以根据数据文件自主选择数据所在列,也可以选择绘图的坐标范围(比如对于全球地面观测数据,可以只选择显示中国地理范围内的数据进行绘图)与填图数据的数值范围(比如可以选择只填图50 mm以上的降水数据等)。由于离散的数据点可能分布比较稠密,因此在显示时可能出现较多数据叠加在一起的情形,导致不能清楚显示离散点数据。新软件增加了自动设置离散点数据显示密度算法,用户放大图形时可显示较多的数据,缩小图形时可显示较少数据;用户可以根据需要自主选择是否使用此功能。
新软件除了在使用探空数据生成的T-logP图中增加了多个物理量值显示外,还增加了风矢端图、位温层结曲线、假相当位温层结曲线和饱和假相当位温层结曲线,并可以对这些曲线设置是否显示和线条类型。新软件可以任意选择和设置气块抬升曲线的起始气压、温度和露点温度。
新软件对地图图形进行了改进和完善,新软件中可以添加各国国界与国内的河流到地图中,也可以选择添加自定义的数据。此外,地图可以选择显示高分辨率或低分辨率地势,或者地势等值线(见图 1中上)。
本软件可以在图形页面中使用鼠标直接绘制多种图形形状(比如直线、折线、多边形、三角形、矩形等),新软件中增加了对这些鼠标绘制图形控制点的编辑修改功能,并增加可以直接与气象图形合并的功能,方便用户直接标注气象图形的重点特征。
新软件中增加了流线图和格点矢量图相互转换功能。流线图属性编辑功能中增加了可设置流线箭头疏密和角度大小的选项。格点矢量图形编辑功能中增加了可以根据图形缩放大小自动设置矢量显示密度的算法,用户放大图形时可显示较多的格点矢量,缩小图形时可显示较少格点矢量;用户也可以选择是否使用此功能。在编辑流线图和格点矢量图形属性时,用户可以根据需要选择是否绘制U、V分量和风速的等值线。
本软件具有对不规则分布离散点数据进行格点化、绘制等值线的功能,新软件改进了该算法。新软件中进行离散点数据格点化时可以选择是否为连续变量(比如温度、气压等)或者突变变量(比如降水、天气现象)。新软件中离散点格点化对话框中还提供了选择数据列的功能,用户可以对多列数据文件自主选择地理坐标和数值所在的数据列。
4 改进的软件操作改进软件操作是为了使软件操作更便捷,包括提供多种途径实现一种操作、批量生成图形图像等。
新软件主窗口工具条按钮中增加了多个线条属性、文字属性、窗口缩放、属性编辑、坐标轴编辑、图形大小编辑等按钮,并有多个按钮其备下拉菜单。在窗口查看方面,既可以选择页面显示,也可以选择不显示页面;既可以选择矩形区域放大到整个图形窗口,也可以拖动鼠标放大,还可以使用鼠标实时移动图形窗口显示内容。用户还可以选择最小化主窗口时隐藏主窗口,在双击Windows开始工具条右下方托盘中相应图标可以恢复主窗口。主窗口中还添加了资源管理器窗口和多剪切板窗口,用户可以在资源管理器窗口中直接启动客观分析进行数据处理或者NCEP数值模式资料处理等,多剪切板窗口可以记录用户的复制操作,类似微软Office中的多剪切板功能。
新软件中增加了用户是否选择把图形中的线条和文字随窗口进行缩放的功能,这是一项非常实用的功能。在用户只需要查看气象物理量场分布细节时,选择不缩放线条和文字功能,这样在放大图形后,线条和文字不跟随放大,可以更清楚地展示气象物理量场的分布特征,这跟MICAPS的缩放功能是一致的;当用户需要制作高分辨率的图像时,需要选择缩放线条和文字功能,这样在图形放大后,线条和文字一起跟随放大,可以获得更高分辨率的图像,这跟Surfer软件的缩放功能是一致的。
新软件中增加了特定等值线和格点图像绘制功能,用户可以根据需求预先配置等值线属性和地图属性,在绘图时直接进行等值线或者格点图像绘制,并自动叠加地图。
新软件中增加了等值线、格点图像、流线、格点矢量和离散点填图等图形的数据编辑功能,可以进行数据编辑或者更改为新数据文件,这样可以在不改变图形设置情况下绘制新图形,方便用户操作。新软件中对所绘制的图形可以保存设置到用户指定的文件中,用户可以在属性编辑对话框中直接打开用户所保存的图形设置文件。此外,在各种图形属性对话框中还增加了右键快捷菜单功能,方便用户选择绘图等值线间隔、填充颜色、线型、字体、颜色等。
在数据处理方面,新软件也做了较大改进。新软件增加了批处理NCEP数据和雷达基数据功能;在格点物理量运算方面增加了多种物理量计算功能,比如大于阈值和小于阈值频率统计、x与y方向平均、计算水汽通量、假相当位温等。
5 未来发展虽然新的客观分析诊断图形系统软件在天气分析方面功能已经较为成熟,但仍有许多不足之处。今后该软件还需要增加极轨气象卫星、MODIS数据和雷达产品数据处理能力;软件的绘制时间演变曲线功能尚不足,T-LogP图功能有待进一步增强;软件需要改进命令行操作,需要增加绘制及输出图形命令和增强图形批处理的能力;软件也需要增加类似MICAPS中的综合图功能以方便用户使用及其同MICAPS软件的兼容性;软件还可以进一步密切与天气业务数据环境的关系,以增强业务数据的支持能力等。
6 总结改进的客观分析诊断图形系统是一款具有较强数据处理能力、使用方便快捷的二维气象绘图和数据分析软件系统。目前该软件主要在国家气象中心强天气预报中心使用。
改进的新软件中增加了多种数据接口,能够处理常规气象观测数据、MICAPS格式数据、静止卫星、新一代天气雷达、VIS5D格式数据和NCEP数值模式数据,以及相应压缩格式的数据等。这进一步提高了该软件的数据处理能力,尤其是提高了强对流天气数据的分析处理能力。
改进的新软件功能得到了较大增强,增加了多种新算法,包括等熵面分析、CTREC矢量计算、TITAN对流风暴追踪、强对流天气的监测和统计、图形裁剪等。新软件中还增加了风矢端图、格点图像、曲线、各种锋面和天气符号等新的图形。
新软件也改进了原软件中的多种算法,包括等值线、流线、矢量、地图、不规则数据分布图等的生成、显示和编辑功能。新软件还改进了软件的数据处理功能和操作方式,使得软件的方便快捷性进一步增强。
[1] |
郑永光, 王洪庆, 陶祖钰, 等. Windows下二维气象绘图软件--客观分析诊断图形系统[J]. 气象, 2002, 28(3): 42-45. DOI:10.7519/j.issn.1000-0526.2002.03.009 |
[2] |
郑永光, 朱佩君, 白洁, 等. Windows下静止卫星云图处理软件[J]. 气象, 2003, 29(6): 16-21. DOI:10.7519/j.issn.1000-0526.2003.06.004 |
[3] |
刘淑媛, 孙建, 郭卫东, 等. 多普勒雷达数据处理显示系统[J]. 气象, 2004, 30(7): 44-46. DOI:10.7519/j.issn.1000-0526.2004.07.010 |
[4] |
郑永光, 陈炯, 王洪庆, 等. 一个气象数据分析绘图软件的设计与开发[J]. 应用气象学报, 2004, 15(4): 506-509. |
[5] |
郑永光, 陈炯, 朱佩君. 改进的静止卫星云图软件处理系统[J]. 气象, 2007, 33(12): 103-109. DOI:10.7519/j.issn.1000-0526.2007.12.016 |
[6] |
李月安, 曹莉, 高嵩, 等. MICAPS预报业务平台现状与发展[J]. 气象, 2010, 36(7): 50-55. DOI:10.7519/j.issn.1000-0526.2010.07.010 |
[7] |
Tuttle J D, Foote G B. Determination of the boundary layer airflow from a single doppler radar[J]. J Atmos Oceanic Technol, 1990, 7(2): 218-232. DOI:10.1175/1520-0426(1990)007<0218:DOTBLA>2.0.CO;2 |
[8] |
Dixon M, Wiener G. TITAN: Thunderstorm identification, tracking, analysis and nowcasting-A radar-based methodology[J]. J Atmos Oceanic Technol, 1993, 10(6): 785-797. DOI:10.1175/1520-0426(1993)010<0785:TTITAA>2.0.CO;2 |
[9] |
韩雷, 郑永光, 王洪庆, 等. 基于数学形态学的三维风暴体自动识别方法研究[J]. 气象学报, 2007, 65(5): 805-814. DOI:10.11676/qxxb2007.076 |
[10] |
郑永光, 张小玲, 周庆亮, 等. 强对流天气短时临近预报业务技术进展与挑战[J]. 气象, 2010, 36(7): 33-42. DOI:10.7519/j.issn.1000-0526.2010.07.008 |
[11] |
郑永光, 张小玲, 周庆亮. 强对流天气短时临近预报技术进展与监测技术[G]. //端义宏, 曲晓波. 2009年灾害性天气预报技术论文集. 北京: 气象出版社, 2010: 87-100.
|