基于ArcObjectss中国大陆地质灾害气象预警系统GHR10.docVIP

基于ArcObjectss的中国大陆地质灾害气象预警系统GHR1.0 薛群威[1] 张徽[2] [1](中国地质环境监测院，北京，100081) [2](南京邮电学院，南京，210003) 一 中国大陆地质灾害气象预警概况 中国国土资源部与中国气象局联合开展的中国大陆地质灾害气象预报预警从2003年6月1日起启动。每年汛期（5-9月份），国土资源部与中国气象局通过中国中央电视台新闻联播之后的时段向社会发布中国大陆地质灾害气象预报预警信息。 中国大陆是地质灾害最为严重的国家之一，地质灾害种类多、分布广、危害大，严重制约着许多地质灾害多发地区的经济发展，威胁着人民生命财产安全。特别是在汛期，受气象因素的影响，崩塌、滑坡、泥石流等突发性地质灾害频繁发生，经常造成重大损失。据统计，近年来，除地震以外的各类地质灾害平均每年造成1000多人死亡、经济财产损失上百亿元。国土资源部与中国气象局将充分利用两部局在地质灾害防治和气象方面的技术力量及成果，提高中国大陆地质灾害监测预报工作的科技含量和预测预报精度，最大限度地避免和减少地质灾害造成的人员伤亡和财产损失。 地质灾害气象预报预警信息由各级国土资源主管部门和气象主管部门联合署名发布，通过当地电视台、广播电台、有线广播等新闻媒体在气象预报节目中播报。预报等级统一划分为5级：1级为可能性很小、2级为可能性较小、3级为可能性较大、4级为可能性大、5级为可能性很大，其中，3级在预报中为注意级，4级在预报中为预警级，5级在预报中为警报级。 二 系统框架 基于ArcObjects的中国大陆地质灾害气象预警系统GHR1.0完成对雨量监测数据、雨量预测数据、地质灾害区划数据、地质灾害隐患点数据、地形数据、地貌数据等的实时获取，并利用相应的算法对数据进行分析，得到地质灾害预警等级。系统的目标是实时生产出中国大陆地质灾害预警等值线。 系统由数据层、分析层、表达层构成，如图2-1。 图2-1 系统结构 基于ArcObjects的中国大陆地质灾害气象预警系统GHR1.0运行中涉及的数据源众多、数据流量大，运行时产生的数据流如图2-2。 图2-2 运行时数据流 三 数据分析 系统数据层数据全部由Microsoft SQL Server管理。雨量监测数据和雨量预测数据来源于气象局，气象局通过Internet把规定格式的雨量监测数据和雨量预测数据以文本形式发送到国土资源部门，然后数据被上载到Microsoft SQL Server进行管理。而地质灾害预警区划数据、地质灾害隐患点数据、地形数据、地貌数据、遥感影像数据、行政区划数据都通过ESRI ArcSDE存取和管理。预警成果数据由Microsoft SQL Server管理，以规定格式文本形式发送到气象局。 地质灾害预警等级等值线的生成依赖于面上的数据，而雨量监测数据、雨量预测数据来自于中国2435个雨量监测点，需要通过表面内插到面上。系统利用ESRI ArcObjects提供的IXYEventSource接口和IInterpolationOp接口完成了这一过程，图3-1显示2005年6月15日8：00和2005年6月16日8：00监测雨量通过IDW方法内插的结果。 图3-1 数据内插 在数据分析方面，系统更充分运用了ESRI ArcObjects提供的众多接口完成各类数据加工处理。 系统在利用地质灾害预警判据运算时利用Map Algebra比较圆满地完成了一系列运算。假设对栅格Δi进行计算，Qi={，，|∈{1,2,3,4,5,6,7},,∈N+,∈N+,j∈N+}构成地质灾害判据空间，其中表示累计雨量天数，，分别表示对应的判据，再令，并记累计雨量 = +++ 其中mk的当日预测雨量为，当日两次监测雨量分别记为和，前q天监测雨量记为，那么最终计算结果就是G(|∈{1,2,3,4,5,6,7})。 对Map Algebra结果的分析发现，由于计算结果分布的空间连续性还不够明显，不宜于进行Reclass计算，所以选择了INeighborhoodOp接口的FocalStatistics方法来进一步加强预警计算结果空间分布的连续性，达到了比较好的效果，图3-2以2005年8月13日预警结果为例，展现进行FocalStatistics前后的效果对比。 图3-2 邻域统计分析 为了对预警结果有一个趋势上的认识，利用上述计算结果通过IRasterConvertHelper接口的ToShapefile方法形成要素类，然后针对要素类通过IInterpolationOp接口的Trend方法进行了趋势分析，图3-3以2005年8月13日为例，展现趋势分析结果。 图3-3 趋势分析 四 预警成果 2005年6月17日GHR1.0预警成果图如图4-1，经过修正，实际预报为浙江大部，江西