1、青海省气象灾害预报预警地理信息系统中的应用2010-12-01 作者：尹振良、蒋红、朱睿、王锡洁、王思维、陈芳、刘勇 来源：2009 年第八届 ESRI 中国用户大会论文集 1 前言最近几十年来，经济和社会的快速发展为气象信息服务提出了日益增长、不断提高的需求。气象信息服务通前言最近几十年来，经济和社会的快速发展为气象信息服务提出了日益增长、不断提高的需求。气象信息服务通常包括决策服务、公众服务、专业服务等几大类别，针对不同的受众，其信息服务的内容和要求呈现多样化、多层次的特点。气象灾害预报预警业已成为气象预报及其信息服务部门日益重要的职责组成部分，并在不断满足全社会各个领域日益增长的需求的基
2、础上向前发展。不同类型的气象灾害因其成因机制不同，影响的范围和频次都有很大的差异，针对不同的地理单元所造成的影响也有显著的不同。满足实际需求的气象信息服务普遍需要地理信息技术的支持。因此地理信息系统与气象信息服务业的有机结合成为目前发展的一个重要趋势。随着 Web Services 技术的兴起，地理信息服务技术应运而生，面向服务的系统架构被放在了一个十分突出的位置上，地理信息共享、分析与可视化在一个新的起点上继续向前发展。针对气象信息服务行业而言，一旦气象预报结果出来，气象信息服务的时效性要求该结果必须立刻提供给多个信息服务部门使用和进一步分析处理，以往基于单机的信息服务系统越来越无法满足实际
3、的需要，有必要发展基于 Web Services、具有网络信息交换与共享服务特点的气象灾害预报预警地理信息系统。本文拟结合青海省气象灾害预报预警地理信息系统探讨地理数据处理在气象灾害预报预警功能实现方面的应用。本系统采用 ArcGIS Server for Java 平台。2 地图服务及地理数据处理服务概述ArcGIS Server 是 ESRI 公司推出的一套后台基于 ArcObjects 搭建的强大的 B/S 开发工具，用于构建集中管理、支持多用户的企业级 GIS 应用平台。它不仅可以提供在线地图发布和在线地图浏览的功能,而且可以提供在线的地址编码服务(Geocode Service)、地
4、理数据处理服务(Geoprocessing Service)、地理数据服务(Geodata Service)和 3D 地图服务(Globe Service)等。本系统在建设过程中实际使用到地图服务和地理数据处理服务的网上发布和功能调用。地图服务是使用最多的一种 ArcGIS Server 服务。该服务可以支持发布二维地图(ArcMap 的 mxd 文档)。通过该服务，用户可以访问以 ArcMap 组织的地图数据和地图表现样式。地图服务中还支持建模操作，OGC WMS 和 KML 数据格式，以及在线编辑空间数据等功能。地图服务部署在 ArcGIS Server 上，运行中需要服务器对象管理器(S
5、OM)和服务器对象容器(SOC)的支持。SOM 主要在地图服务的管理、启动和关闭以及地图服务的运行服务器(SOC)的添加、删除及负载均衡方面起着专门的作用。SOC 主要负责运行地图服务并处理应用层提交的请求。地理数据处理服务就是将在 ArcGIS Server 服务器端建设好的地理数据处理模型或者包含一个地理数据处理模型工具层的地图文档发布为一个服务。地理数据处理模型是一个定义空间模型或地理数据处理工作流的工具，可以用可视化工具 Model Builder 创建地理数据处理模型，也可以用文本编程以脚本的方式创建模型。所有的地理数据处理工具可以用作创建模型中的处理，一个模型包含一个或多个处理，它
6、们可以链接在一起，也可以不链接在一起。一个处理包含一个工具、输入和输出。通过构建地理数据处理模型，可以自动完成地理数据处理工作流。地理数据处理服务是一个基于 Web 的地理数据处理工具，客户端提交处理请求，服务器执行空间分析和建模，然后把执行结果展现在客户端，它便于组织内部数据的集中管理和操作，实现了功能的共享。ArcGIS Server 的服务可以通过 Local 或者 Internet 两种方式进行连接。Local 方式直接连接到 SOM 上，通过 AO 进行交互，所以必须在本地有 AO 对象才可以进行连接。Internet 方式直接连接到 Web Service 的引用地址，它是通过本地
7、对象连接的。对于 Java ADF 而言，本地对象表示连接 ArcGIS Server 的类存在于本地 JRE 中。ArcGIS Server 在发布每一个服务时，也同时发布了一个 Web Service，因此我们可以通过 Web Services 的方式来直接访问 ArcGIS Server 上的服务，这样极大的提高了部署的灵活性。用户可以通过在 Java 中创建代理类的方式访问这个 Web 服务，如下图（图 1）所示。本系统采用Internet 方式在 Java 中创建代理类的方式来访问 ArcGIS Server 上的地图服务和地理数据处理服务。图 1.ArcGIS Server 服务的
8、两种连接方式3 系统体系结构本系统基于 Web Services 的理念，采用 ArcGIS Server 9.3 + SQL Server + J2EE 软件开发模式搭建了一个 B/S 架构的业务运行系统。整个系统自下而上可分为数据层、服务供给层和应用层三个部分（图 2）。数据层主要包括基础地理数据库和气象数据库，存储在数据库服务器中，采用 SQL Server 进行气象历史数据、预报数据的有效存储与管理。地理数据采用 ESRI Geodatabase 数据模型，通过空间数据库引擎 ArcSDE 保存在 SQL Server 中。服务供给层的作用是从数据库中提取所需要的气象数据、地理数据，然
9、后通过地理数据处理服务对数据进行处理分析，形成分别满足决策气象服务、公众气象服务、专业气象服务和气象信息发布需求的气象信息产品。这里需要用到地理信息系统空间数据及属性数据编辑功能、矢量数据叠置分析功能、地统计分析和栅格数据分析功能以及地图制图功能。这些功能通过 ArcGIS Server Web 应用程序定制或者地理数据处理服务来实现。应用层的主要内容就是定制一个易学易用的用户界面，用以实现有关气象信息与地理信息组合产品的网上发布，并进行用户访问控制，根据用户职能分别赋予不同的访问权限。图 2.系统体系结构示意图4 基础地理数据库本数据库包括基础地理数据、气象数据两个方面。其中气象预报结果需要
10、经过处理和变换生成的矢量图层数据库和气象灾情服务信息数据库。数据库的组织采用 ArcGIS 空间数据库 Geodatabase 体系结构。在 Geodatabase 中，所有基础地理矢量数据、数字高程模型、遥感影像和数据表单均保存在一个地理数据库中，并组成层次分支结构。在地理数据库中包含一个数据集、一个数字高程模型、若干不同分辨率遥感影像，以及若干数据表单。数据集中包括所有的基础地理矢量数据。本地理数据库内容包括基础地理数据、数字高程模型、卫星影像数据，以及有关表单。所有数据存储在一个 Geodatabase 中，位于分支树顶端，以下包含一个数据集 Albers，存储所有矢量数据。与该数据集并
11、列的是各种影像数据、数字高程模型和表单数据。这就是说，所有矢量、栅格数据和影像数据均进行了地图投影转换，采用亚伯斯（Albers）等积圆锥投影进行数据存储和管理。相应地，其中特征属性表中的面积项是有意义的。地理数据库包括以矢量方式分层存储和管理的境界线、水系、交通、居民地、地貌等要素。数据集的分层命名与国家地形数据库数据命名方式一致。相应建立分县、乡（镇）社会经济统计数据库并通过通用政区编码字段与地理数据库中的图层建立链接关系，并参与空间数据的分析。5 系统功能根据需求，本系统需要实现如下功能：(1) 基本信息查询检索功能：在 Web 页面上所有用户可以通过放大、缩小、漫游、查询、度量等工具进
12、行青海省基础地理图层的操作，可以将有关图层显示和关闭，可以查询特征图层的属性数据表。(2) 气象数据导入及其处理：系统在气象台局域网中与 MICAPS 系统文件夹下的数据集保持链接关系，有权限的用户可以通过菜单和弹出对话框方式选择要导入的气象数据类型（如气温）、设置插值方法、设置等值线间距等参数。在服务器端完成数据格式从文本文件到 ESRI Shapefile 数据格式的转换，并实现气温、降水量、气压等气象要素的站点记录在地图上的显示，以及风向风速的风向杆标注；系统还可以调用服务器上的地理数据处理服务对站点数据进行插值等操作，最终以等值线或者等值面的形式显示在地图上。(3) 灾情范围的确定与在
13、线勾绘：将数据编辑功能集成到网页中，授权用户根据天气形势分析结果，在网页上用鼠标手动编辑、交互式绘制灾害影响区域，用于发布灾害预警或者灾害损失评估分析。(4) 暴雨洪水预报：操作人员依据天气预报结果勾绘暴雨范围，系统将此暴雨区域与地理数据中的流域分区图层进行叠置分析，得到暴雨落地各个流域的范围大小，形成结果表单。进而如果知道暴雨量大小，就可以估计产生洪水径流量的多寡。(5) 灾害损失评估：操作人员依据天气预报或气象分析结果勾绘冰雹、干旱、暴雪等灾害性天气的影响范围，然后将该图层与政区图、草场资源分布图等进行叠置分析，结合政区图、草场资源分布图以及相关的统计数据分析灾害性天气可能造成的人口、财产
14、损失，结果以表单形式存储。(6) 制图打印功能：以图像形式保存结果分析地图页面，或直接送打印机输出。(7) 灾情信息管理与快报生成：对参与灾情分析的气象站点进行管理，添加新站点，删除不再发挥作用的站点；进行灾情信息资源管理，添加新的气象灾情报告记录到数据库中，并对已有报告记录进行数据库维护；添加或修改灾情防止对策数据库记录；在这些数据库建设的基础上根据灾情预警决策的需要自动完成气象灾情快报的制作。(8) 用户管理：将系统用户分成系统管理员、数据库管理员、灾情评估分析员和普通用户四种类型，不同类型的用户赋予不同的操作权限。本模块可以创建、删除用户，并赋予相应的权限。(9) 基于 Google E
15、arth API 的天气形势分析与预警模块：本模块基于 Google Earth API，可充分有效地利用 Google Earth 丰富的地理信息和遥感影像资源，授权用户在气象数据导入完毕、灾害影响区域手动勾绘完毕后可以通过点击菜单将气象数据和灾害预警区域导出到三维显示平台，进而可以在三维场景下浏览气象数据和灾害预警数据。本模块集成本系统所建立的青海省基础地理数据集，实时叠加和更新等温线、等压线、降雨色斑图、灾害预警区域等数据，以及实时显示 FY-2 卫星云图、环青海湖公路自行车赛天气预报等功能。6 地理数据处理建模及其服务发布6.1 数据格式转换及数据导入处理过程由于原始的气象观测数据是以
16、文本格式、单站点记录形式保存，无法实现在地图上标注及生成等值线等操作，进而也无法满足气象部门工作人员对气象数据的分析与判断，所以需要先将原始的气象数据转换成 ESRI Shapefile 数据格式，在此基础上再经过插值、平滑处理、定义投影坐标系、裁切、导入Geodatabase 等一些列地理数据处理过程实现气象数据以不同的方式在地图上标注与显示，数据格式转换处理流程如下图所示（图 3）。在本系统中，这一功能主要是通过建立地理数据处理模型发布地理数据处理服务来实现的。图 3. 数据格式转换处理流程在系统开发过程中，根据实际调研，需要将温度、降雨、气压、风向风速四类气象台站数据转换为地理数据在地图
17、上进行标注显示。这四类数据中温度、降雨、气压既需要在地图上标注出单个观测站点的数据也需要把经过地理数据处理生成的等温线、等压线和降雨量色斑图显示在地图上；风向风速数据只需要将单个站点的观测数据用风向杆根据风速的大小在地图上标注出来即可。所以，系统建立了站点观测数据导入模型、等值线数据导入模型、降雨量数据导入导入模型、风向风速数据导入模型，以满足四类气象数据的格式转换与在地图上标注显示的要求。站点观测数据导入模型（图 4）实现如下功能：把温度、降雨、气压三类气象数据的站点 Shapefile 数据作为模型的输入数据，经过模型的处理操作后将输出数据拷贝到显示图层，进而完成站点观测数据在地图上的标注
18、显示。系统在实现过程中，在地图文档中添加三个点状图层分别作为温度、降雨、气压三类气象数据的站点观测数据的显示图层。模型每次运行时，根据气象数据的类型决定把导入的站点观测数据拷贝到哪个显示图层。模型处理完成后，刷新地图页面即可实现新导入的数据在地图上显示。图 4. 站点观测数据导入模型等值线数据导入模型（图 5）实现如下功能：把温度、气压两类气象数据的站点 Shapefile 数据作为模型的输入数据，经过模型的处理操作后将输出数据拷贝到显示图层，进而完成温度、气压数据以等值线的形式在地图上显示。图 5. 等值线数据导入模型降雨量数据导入模型（图 6）实现如下功能：把 12 小时降雨、24 小时降
19、雨、过程降雨三类降雨数据的站点 Shapefile 数据作为模型的输入数据，经过模型的处理操作后将输出数据拷贝到显示图层，完成降雨数据以降雨量色斑图的形式在地图上显示。图 6. 降雨量数据导入模型风向风速数据导入模型（图 7）实现如下功能：把风向风速数据的站点 Shapiefile 数据作为模型的输入数据，经过模型的处理操作后将输出数据拷贝到显示图层，完成风向风速数据以风向杆形式在地图上显示。其中风向杆是通过风向杆字体 ESRI Weather 标注实现的，风向杆的旋转角度则是从原始数据中经过数据格式转换得到的旋转角度；风速的大小换算成风向杆字体大小的等级，再加上字体的其实 ASCII 值，便
20、得到风速所对应的风向杆字体的 ASCII 码值，最后将 ASCII 码值赋给 speed 字段，用以在地图上标注显示。图 7. 风向风速数据导入模型6.2 灾情分析处理过程本系统实现的灾情分析功能也是通过建立相应的地理处理模型实现的。灾情分析过程为：气象部门工作人员先手动在相应的灾害类型数据层上在线勾绘出灾害影响区域然后与 Geodatabase 中的行政区划数据层或者流域分区数据层进行叠置分析(Intersect)，再经过自动计算受灾面积、生成 DBF 数据表、删除临时数据等地理数据处理过程，最后通过程序访问 DBF 表中的相关数据列显示在客户端页面上，即完成了灾害影响区域分析。灾情分析处理
21、过程如下图所示（图8）。图 8. 灾情分析处理流程青海省境内经常发生的气象灾害及气象因素诱发的灾害主要有雪灾、暴雨洪涝及强降水、冰雹、干旱、雷电、低温冻害、大风、沙尘暴、山体滑坡和泥石流等九种灾害类型。系统在实现过程中根据这九类灾害在 Geodatabase 中建立九个多边形图层，每个多边形图层对应一种灾害类型。气象部门工作人员根据天气预报结果判断将要发生灾害时，即在相应的灾害类型图层上手动编辑出灾害影响区域，进而灾害区域在地图上根据不同的灾害类型用不同的颜色显示，达到灾害预警的目的。系统的灾情分析处理主要是通过灾情分析处理模型实现的，灾情分析处理模型如下图（图 9）所示。灾害影响区域作为模型
22、的一个输入数据，行政区划或者流域分区作为模型的另一个输入数据，将两个输入数据进行叠置分析。模型在运行时根据灾害类型通过程序控制参与叠置分析的灾害影响区域数据层，同理行政区划层和流域分区层也通过程序控制哪一层参与叠置分析。叠置分析后，经过自动计算受灾面积、生成 DBF 属性表、删除临时数据等操作完成灾情分析处理。最终灾情分析结果将以统计报表形式在客户端显示。图 9. 灾情分析处理模型6.3 气象数据及灾害预警区域导出处理过程系统本部分通过建立地理处理模型实现了将等温线、等压线、降雨量色斑图等气象数据和灾害预警区域数据导出成 KMZ 格式数据，以便气象数据和灾情预警区域在基于 Google Ear
23、th API 的三维显示平台上进行显示。数据导出过程如下图所示（图 10）。系统将经过数据导入处理的气象数据和灾害预警区域数据拷贝到显示图层，添加到地图文档中。在地图文档中设置好数据的显示风格，然后地图文档作为地理处理模型的输入数据，经过地理处理过程转成 KMZ 格式数据，即可在三维显示平台浏览导出的数据。数据导出处理模型如下图所示（图 11）。图 10. 数据导出处理过程 图 11. 数据导出处理模型6.4 地理数据处理服务发布过程在将地理数据处理模型发布为地理数据处理服务时，ArcGIS Server 对模型的输入、输出参数的数据类型有特殊的要求。以上所述的模型中，FeatureClass
24、 数据类型作为模型的输入参数，将模型发布为地理数据处理服务时系统就会报错，无法发布服务。本系统的解决方法是：先在 Model Builder中建立模型，将模型集成在工具箱中；第二步建立脚本模型，通过脚本调用工具箱中的模型，脚本模型的输入、输出参数全部设置为 String 类型；每一个模型对应建立一个脚本模型，将脚本模型同样集成在工具箱中，最后将脚本模型工具箱发布为地理数据处理服务。这样就满足了地理数据处理服务对输入、输出参数数据类型的特殊要求。以等值线数据导入脚本模型为例，脚本代码如下图所示（图 12）。图 12. 等值线数据导入脚本模型在 ArcCatalog 中，把建好的脚本模型工具箱，发
25、布为地理数据处理服务。系统在运行时，用户在客户端浏览器通过发送访问请求到 Web 服务器，由 Web服务器通过程序调用 GIS Server 上的地理数据处理服务，在 GIS Server 上完成数据的分析与处理，最后通过 Web 服务器把结果返回到客户端浏览器。这样就实现了地理数据处理服务的共享。图 13. 气象要素在地图上显示图 14. 灾情预警及评估分析图 15. 灾情预警区域三维显示7 结论本文基于 ArcGIS Server 平台开发了青海省气象灾害预报预警地理信息系统，对地理数据处理服务在气象灾害预报预警方面的应用进行了探索。探索表明，在把原始气象数据转换成地理信息数据并在地图上显
26、示的过程中，建立合适的地理数据处理模型并将模型发布成为地理数据处理服务是一种非常有效的方法。地理数据处理服务在 WebGIS 数据处理方面发挥了重要的作用。系统达到了预计的功基于 ArcGIS 的地质灾害预报预警及应急指挥系统1.前言地质灾害预报预警及应急指挥系统是面向国土资源部门应用，以 ArcGIS 为技术平台，集地质灾害预警管理、地质灾害应急管理、地质灾害危险评估、地质灾害点监测管理和地质灾害信息发布管理等于一体的地质灾害防治综合管理系统。2.总体设计1）系统目标根据国家和各地区地质灾害管理工作的有关规定，依照各地方地质灾害管理工作流程，以计算机技术为基础，以数据库为核心，以 ArcGI
27、S 为平台，以网络传输为通讯手段，建立地质灾害预报预警及应急指挥系统，实现地质灾害预报预警分析及成果签批发布的自动化、科学化；实现地质灾害调查评价、规划、监测、预报预警、应急处置、搬迁治理等地质灾害数据的存储和共享利用；对突发性地质灾害作出分析、评估，为各级政府快速组织相关资源，实施抢险救灾，开展地质灾害防灾减灾工作提供决策依据和技术服务。2）技术路线本系统对地质灾害的发育特点、地质灾害发生的内在（如地形地貌、地质构造、地层岩性等）和外在致灾因素（如降雨、人类工程活动、植被破坏等）进行研究分析，结合现有地质灾害调查基础资料及各类关联因素数据进行运算分析，建立地质灾害预报预警模型，并通过充分利用
28、 ArcGIS 平台提供的 GIS 技术手段，综合地质灾害预报预警模型的基础因子和诱发因子，生成地质灾害预报预警等级图。整个预报预警模型开发的技术路线如下图所示：图：总体技术路线图在本系统的研究开发过程中，对 GIS 技术的应用尤为重要，国土资源系统的核心是数据，海量的空间和属性数据势必要求通过具有海量数据管理能力的 GIS 软件来完成，由于 ArcSDE 本身所具有的海量数据存储、多用户并发访问、版本管理、长事务处理等强大优势，在本系统中引入 ArcSDE 作为空间数据存储和管理引擎，它利用 Oracle 关系数据库在数据存储、数据完整性等方面的先进技术手段，将海量空间数据（包括矢量数据和栅
29、格数据）有机地组织和管理起来，通过其内部异步缓冲、空间索引等先进的机制，提供对空间数据的多用户高效并发访问，为本系统的高性能品质提供了有力的保障。3.系统介绍1）系统总体架构系统采用“两网、一库、六系统”的应用架构，两网指利用现有的国土资源广域网和 Internet 外网网络资源；一库指建设地质环境与地质灾害信息库作为数据支撑；六系统指地质灾害数据库管理系统、地质灾害自动化预报预警系统、地质灾害综合管理信息系统、外网地质灾害信息发布系统、移动地质灾害信息查询系统及国土资源远程视频会商系统，系统间相互集成、共享数据、结合应用。系统总体架构如下图所示：图：系统总体架构图地质灾害数据库管理系统，提供
30、基础数据的加工处理以及与气象部门、水利部门的气象、雨量、水利数据共享等，同时为其它相关部门业务系统（如水利部门的山洪预报系统）的运行提供数据支撑；地质灾害自动化预报预警系统用于对地质灾害发生概率等级的分析计算，实现预报预警成果输出，为地质灾害应急指挥和受害群众的转移避让预案决策提供依据；地质灾害综合管理系统和外网地质灾害信息发布系统实现对地质灾害信息应用管理，以及提供灾害预警发布、应急指挥和应急处置等功能，两者面向的用户不同分别部署。地质灾害综合管理信息系统通过政务信息网联接各级国土资源部门及相关政府部门，外网地质灾害信息发布系统从地质灾害综合管理系统中功能提取，运行在外网 INTERNET
31、环境中，适用于乡（镇）国土所用户及社会公众。移动地质灾害信息查询系统通过无线网络满足地质灾害防治工作人员随时、随地了解最新地质灾害信息和及时上报地质灾害信息的需要。应用系统与远程视频会商系统相连，实现综合会商。2）系统功能 预报预警 综合有关地质和地理基础信息、地质灾害信息和降雨量、人类工程活动程度等因素，从不同空间和时间尺度上分析地质灾害的发生概率和分布范围，从而做出较为准确的区域地质灾害预报预警。 地质灾害信息管理 对各地质灾害隐患点信息进行全面管理；对地质灾害点灾情信息进行综合管理；提供地质灾害点分布图 GIS 功能。构建立体三维景观，逼真反映灾点所在位置的地形、地貌，并利用三维景观分析地质灾害点的空间信息。 地质灾害监测管理 对群测群防信息、群测群防监测数据进行管理。支持多种监测数据上报方式，支持以网络上报或短信上报监测数据。 预报预警信息发布 对预报预警信息进行发布。 应急处置管理 对应急预案、地质灾害防治信息、群众转移预案及地质灾害点临时避让的信息进行管理。 地灾评估治理 进行地灾评估管理和地灾治理。
