孫哲趙曦賈寧（陜西省地震應急中心，西安 710069）

基于“天地圖”的地震專題圖自動出圖技術(shù)研究與實現(xiàn)——以陜西省為例1

孫哲趙曦賈寧
（陜西省地震應急中心，西安 710069）

孫哲，趙曦，賈寧，2016.基于“天地圖”的地震專題圖自動出圖技術(shù)研究與實現(xiàn)——以陜西省為例.震災防御技術(shù)，11（2）：396—402.doi：10.11899/zzfy20160223

結(jié)合地震應急中對專題圖產(chǎn)出的要求，研究“天地圖-陜西”的切片方式和服務方式，在ArcGIS中實現(xiàn)了天地圖瓦片數(shù)據(jù)與矢量數(shù)據(jù)相結(jié)合。通過研究數(shù)據(jù)庫監(jiān)聽技術(shù)、基于Python語言和ArcPy包的出圖技術(shù)，實現(xiàn)地震三要素的自動獲取和影響場產(chǎn)出，并最終實現(xiàn)了基于天地圖數(shù)據(jù)的地震專題圖自動產(chǎn)出，完全無需人工操作，提高了專題圖出圖效率。

ArcPy天地圖地震專題圖

引言

地震作為一種突發(fā)事件具有影響范圍廣、破壞性大、災害信息復雜、災害種類多樣及救援難度大的特點（魏艷旭等，2014），因此當?shù)卣鸢l(fā)生時，人們往往希望盡可能多地掌握震區(qū)人口、經(jīng)濟等情況，為隨后地震應急救援工作的開展提供準確和全面的參考。地震專題圖件因其直觀、信息量大的優(yōu)點，能及時有效地提供各類信息，在地震應急救援工作中發(fā)揮著重要作用（李金香等，2013）。尤其是根據(jù)地震參數(shù)采用適合震區(qū)的地震烈度衰減關(guān)系繪制成的地震影響范圍圖，為第一時間對災區(qū)的影響范圍和傷亡損失作出初步的評估提供了重要支撐（和朝霞等，2013），因此研究如何在地震發(fā)生后最短時間內(nèi)快速有效地產(chǎn)出地震專題圖具有重要意義。

目前普遍采用的專題圖出圖的方式有兩種，一種是通過ArcMap調(diào)用各種矢量數(shù)據(jù)和屬性數(shù)據(jù)，并在此基礎上實現(xiàn)快速出圖，但是這種出圖方式對制圖人員的專業(yè)要求較高，需要制圖人員具有一定的地理基礎知識并且熟悉相關(guān)出圖軟件，同時地震發(fā)生初期需要盡快了解震區(qū)相關(guān)情況，因此要求專題圖的產(chǎn)出盡可能快速，而人工出圖需要較長的時間制作專題圖，無法滿足快速出圖的要求。另一種是基于瀏覽器調(diào)用網(wǎng)絡地圖（例如Google等）并在此基礎上產(chǎn)出專題圖，這種方式提供了豐富的地理信息，更利于專題圖的產(chǎn)出，但是這種方式產(chǎn)出的專題圖在比例尺、圖例、符號化等圖面整飾方面有諸多的限制，同時還依賴于網(wǎng)絡的穩(wěn)定性。以上兩種方式各有優(yōu)缺點，使用時也各受其限。為此，本文開發(fā)出一種出圖方式以解決以上兩種方式的限制，并通過使用公共地圖擺脫保密數(shù)據(jù)限制的同時，再通過實現(xiàn)自動整飾、三要素自動獲取等技術(shù)減少人工干預，在滿足地震應急需求的基礎上實現(xiàn)地震專題圖的自動出圖，提高了專題圖的產(chǎn)出效率和出圖速度。

1　天地圖概述

2011年國家地理信息公共服務平臺——“天地圖”正式上線?！疤斓貓D”是基于互聯(lián)網(wǎng)的基礎地理信息服務網(wǎng)站，集成了海量基礎地理信息資源（陳鑫祥等，2013）。“天地圖”在服務器端采用金字塔式的瓦片存儲方式，將地圖按照不同級別切分為相同大小的切片，對相同級別的切片按照行列號規(guī)則統(tǒng)一存儲（圖1）（牛瑩等，2014）。

2013年陜西地震局與陜西地理信息中心簽署協(xié)議，在陜西地震局部署了“天地圖-陜西”地圖服務?！疤斓貓D-陜西”地圖服務分為切片服務和實體服務兩類，切片服務每一層級按照不同的比例尺存放，由上到下比例尺越來越大，最終構(gòu)成金字塔結(jié)構(gòu)進行切片組織和管理。每個切片都有對應行列值編號，系統(tǒng)通過調(diào)用相應編號來訪問對應的切片。實體服務分為行政區(qū)劃、交通、POI（興趣點）等服務。通過將這些服務與陜西地震局收集到的相關(guān)地理要素進行關(guān)聯(lián)和匹配，建立地震應急專題要素數(shù)據(jù)庫，為地震專題圖的產(chǎn)出提供非涉密的數(shù)據(jù)支撐。

圖1　天地圖切片組織結(jié)構(gòu)Fig.1　The organization structure of TianDiTu（Map World）tiles

2　系統(tǒng)設計

本文涉及的專題圖自動出圖系統(tǒng)基于MyEclipse平臺，系統(tǒng)分為3個組成部分：數(shù)據(jù)庫采用Oracle11g，數(shù)據(jù)庫監(jiān)聽采用Java語言實現(xiàn)，出圖腳本采用Python腳本語言和ArcPy包構(gòu)建。

2.1ArcPy語言簡介

ArcPy是ESRI公司從ArcGIS10開始發(fā)布的基于Python腳本語言的站點包，目的是為以實用高效的方式通過Python執(zhí)行地理數(shù)據(jù)分析、數(shù)據(jù)轉(zhuǎn)換、數(shù)據(jù)管理和地圖自動化產(chǎn)出等功能。ArcPy集成了Python的特點：易于跨平臺操作、占有資源較少、速度快，在自動產(chǎn)出、快速生成方面具有較大優(yōu)勢（周津津等，2014），因此本文用Python語言構(gòu)建設計專題圖自動出圖系統(tǒng)，可縮短出圖時間以滿足地震應急的需求。

2.2技術(shù)路線

地震發(fā)生后，陜西省地震局臺網(wǎng)中心在進行地震速報的同時，將速報結(jié)果寫入到地震速報信息中轉(zhuǎn)表中。系統(tǒng)通過數(shù)據(jù)庫監(jiān)聽技術(shù)定時監(jiān)控數(shù)據(jù)庫中的記錄，當發(fā)現(xiàn)一條新紀錄后觸發(fā)出圖腳本。腳本被觸發(fā)后從數(shù)據(jù)庫中讀取地震三要素，首先根據(jù)地震位置進行是否需要生成地震影響場的判斷，并確定相關(guān)影響場參數(shù)后生成地震影響場和震中。其次再按照地震位置調(diào)用省內(nèi)、全國和全球不同專題圖模板，對圖名、比例尺、出圖范圍等要素進行修改整合后，導出專題圖完成出圖，并上傳產(chǎn)出圖件至共享FTP用于分發(fā)使用（圖2）。

圖2　系統(tǒng)技術(shù)路線圖Fig.2　Technical route in the system

2.3系統(tǒng)平臺構(gòu)建

系統(tǒng)編程基于MyEclipse平臺，數(shù)據(jù)庫使用Oracle11g，GIS平臺使用ArcGIS軟件。數(shù)據(jù)庫監(jiān)聽基于Java語言通過TimerTask函數(shù)調(diào)用實現(xiàn)，自動出圖使用Python腳本語言編寫，影響場的生成、專題圖整飾、專題圖產(chǎn)出等GIS功能通過調(diào)用ArcPy包進行二次開發(fā)完成，通過cx_Oracle插件的調(diào)用實現(xiàn)對數(shù)據(jù)庫的訪問和操作，產(chǎn)出的專題圖通過Commnos包上傳至共享FTP。

3　系統(tǒng)實現(xiàn)

3.1模板分類

因為地震震級不同、專題圖是否產(chǎn)出影響場以及產(chǎn)出影響場的范圍不同，地震專題圖所涉及的專題圖符號、地圖比例尺、專題內(nèi)容的詳細程度等也有不同，因此需要按照地震震級和專題圖的用途對專題圖模板進行分類，提前制作好專題圖模板。本系統(tǒng)根據(jù)中國地震局《破壞性地震應急專題地圖產(chǎn)出流程與制圖規(guī)范（試行）》對專題圖的要求，結(jié)合陜西省天地圖自身地理信息情況完成專題圖模板的分類和制作，具體分類和模板樣例見圖3和圖4。

圖3　模板分類Fig.3　The classification of thematic map

圖4　專題圖模板Fig.4　The template of thematic map

3.2數(shù)據(jù)庫的建立

建立地震速報數(shù)據(jù)庫表，用于獲取和存儲地震三要素，數(shù)據(jù)庫表結(jié)構(gòu)如表1所示。

表1　地震速報數(shù)據(jù)庫表Table 1　Data table of the earthquake quick report

3.3地震三要素自動獲取

系統(tǒng)使用TimeTask實現(xiàn)對地震速報庫循環(huán)掃描，并判斷是否有新的記錄。初步定義系統(tǒng)每隔三分鐘執(zhí)行一次數(shù)據(jù)庫掃描。當?shù)卣鸢l(fā)生后由陜西省地震局臺網(wǎng)中心確定地震三要素并寫入地震速報數(shù)據(jù)庫中，系統(tǒng)監(jiān)測到數(shù)據(jù)庫有新記錄后調(diào)用Python腳本進行判斷是否需要產(chǎn)出專題圖，對于正式速報地震進行專題圖產(chǎn)出工作，對于余震或自動速報則不觸發(fā)響應。

3.4影響場產(chǎn)出

本系統(tǒng)使用的地震影響場模型采用橢圓衰減模型（見公式1），并收集不同地區(qū)相關(guān)的衰減參數(shù)。系統(tǒng)獲取地震三要素后首先判斷地震落入到哪個分區(qū)，提取相應的衰減參數(shù)后按照公式計算生成影響場，并保存在fileGDB中，出圖時將其導入到MXD模板中，最后成圖。

3.5專題圖整飾和產(chǎn)出

系統(tǒng)首先對地震三要素進行空間分析判斷，以此根據(jù)提前規(guī)劃好模板的分類方式確定產(chǎn)出專題圖的模板。然后對于每一個模板通過ListLayoutElements函數(shù)獲取到模板圖名、圖例、產(chǎn)出時間等圖面要素并對其進行修改。通過extent函數(shù)確定專題圖的出圖范圍，具有影響場的專題圖出圖范圍為影響場最外圈外切矩形的范圍，對于無影響場產(chǎn)出的地震，其出圖范圍設定為模板原始范圍。最后通過ExportToJPEG函數(shù)設定出圖路徑、文件名稱、分辨率等要素并完成出圖。

4　系統(tǒng)應用與測試

最后通過測試地震數(shù)據(jù)，對自動出圖系統(tǒng)進行測試。測試地震數(shù)據(jù)為：2015年6月11 日9時00分00秒陜西省西安市周至縣6.5級地震（108.02°E，34.13°N）。圖5為9點00分00秒模擬臺網(wǎng)記錄在數(shù)據(jù)庫中寫入地震信息。

圖5　在數(shù)據(jù)庫中模擬一個地震信息Fig.5　Simulating a seismic information in the database

運行專題圖自動出圖軟件，監(jiān)視系統(tǒng)運行結(jié)果如圖6，09:00:20系統(tǒng)觸發(fā)成功，自動開始制作專題圖，09:05:20專題圖制作完成（圖7），產(chǎn)出A1、A3、A4共3個幅面40張專題圖。

圖6　系統(tǒng)運行信息Fig.6　System running information

圖7　產(chǎn)出的地震專題圖Fig.7　The outputting seismic thematic map

5　結(jié)語

本文通過對天地圖-陜西地圖服務、數(shù)據(jù)庫監(jiān)聽和Python自動出圖相關(guān)技術(shù)的研究，構(gòu)建基于天地圖的地理要素數(shù)據(jù)庫，并在此基礎上編程實現(xiàn)了基于天地圖的地震專題圖快速出圖。系統(tǒng)避免了人工干預，在應急人員到崗之前就可以完成專題圖的制作，提高了專題圖產(chǎn)出效率。同時使用天地圖數(shù)據(jù)避免了使用保密數(shù)據(jù)帶來的專題圖在共享和分發(fā)方面的限制，擴充了專題圖的使用范圍，方便了專題圖的共享。

陳鑫祥，吳錦超，2013.利用ArcGIS Engine實現(xiàn)天地圖瓦片數(shù)據(jù)的切割.北京測繪，（4）：17—20.

和朝霞，賈寧，2013.基于ArcGIS的地震專題圖快速生成技術(shù)研究.四川地震，（1）：45—47.

李金香，溫少妍，張琳琳等，2013.新疆地震應急救援快速圖研究.高原地震，25（2）：56—59.

?，?，劉亮，2014.基于C/S的網(wǎng)絡地圖的研究與實現(xiàn)——以天地圖為例.甘肅科技，30（18）：38—41.

魏艷旭，王寶坤，賈軍鵬等，2014.基于AE的地震應急專題圖快速生成系統(tǒng)設計與實現(xiàn).價值工程，（8）：229—230.

周津津，陳少鋒，劉曉娟，2014.淺談基于ArcGIS的影像處理工具集開發(fā).37（11）：155—157.

Research and Implementation ofAutomatic Mapping for Seismic Thematic Map Based on TianDiTu（Map World）of Shaanxi Province

Sun Zhe，Zhao Xi and Jia Ning
（Earthquake Administration of Shaanxi Province，Xi’an710069，China）

For the requirements of earthquake emergency response with seismic thematic mapping，it’s necessary to combine the tile map data and vector data in ArcGIS for studying the slice and the service of the“Map World-Shaanxi”.Here，we can automatically achieve the three elements of earthquake and influence field output through the technique of database monitoring technology and thematic mapping technology based on the python language and arcpy package，and finally achieve the seismic thematic mapping based on（TianDiTu）.This study tends to the improvement of the production efficiency without manual operation.

Arcpy；TianDiTu（Map World）；Seismic thematic map

地震應急青年重點任務（CEA_EDEM_201508）、2014年陜西省地震局啟航與創(chuàng)新基金（201306）資助

2015-07-06

孫哲，男，生于1984年。碩士，工程師。主要從事地震應急研究。E-mail：sz@eqsn.gov.cn