地質災害三維遙感解譯成果的自動制圖方法

張明洋，馬維峰，唐湘丹，李仙偉

(1．中國地質大學(武漢)地球科學學院，武漢 430074;2．武漢地大信息工程有限公司，武漢 430074)

0 引言

地質災害遙感解譯應用于地質災害調查具有及時、快速、準確和經濟的特點，已廣泛應用于地質災害日常調查、監(jiān)測以及突發(fā)性地質災害的應急調查、動態(tài)監(jiān)測和損失評估中。在2008年“5·12”汶川地震災害中，地震誘發(fā)了大量的次生地質災害，中國地質調查局等單位在獲得的研究區(qū)震后多源遙感影像的基礎上進行了快速的遙感解譯，并繪制了研究區(qū)的地質災害分布圖，為進一步研究汶川地震次生地質災害的發(fā)育規(guī)律、發(fā)生機理、預測理論、災區(qū)恢復重建與選址等問題提供了重要的基礎信息［1］。在2010年“8·7”甘肅舟曲泥石流災害中，中國國土資源航空物探遙感中心利用災后航空攝影數據，結合國家基礎地理信息數據庫，快速進行了舟曲縣的地質災害遙感解譯，進而展開對比分析，獲得了倒塌城鎮(zhèn)房屋，受害公路、橋梁、農田、林區(qū)以及堰塞湖、泥石流、滑坡等各種災情的位置、類型、規(guī)模、分布特征等信息，為救援工作提供了有力依據，贏得了寶貴的搶險救災時間［2］。

地質災害遙感解譯成果的信息產品和表現形式主要是遙感影像專題圖，無論使用什么樣的制圖軟件，都需要將解譯結果進行繁瑣的編輯、渲染、符號化甚至數據格式轉換，需要熟悉各類軟件操作，費時費力［3］。相對于二維環(huán)境，三維環(huán)境下的地質災害遙感解譯能夠將高精度正射遙感影像與DEM數據疊加建立三維地形，從而提供真實、直觀的地理地貌環(huán)境，更有利于地質災害信息交互解譯和信息提取，已在具體工作中得到了廣泛應用［4－6］。而三維環(huán)境圖層多、信息量豐富，更難以通過第三方 GIS或CAD軟件進行制圖輸出。

本文以中國地質調查局遙感解譯專題圖相關制圖規(guī)范為依據［7－8］，研究在三維環(huán)境下地質災害遙感解譯成果的專題圖層設置、符號化、渲染和專題圖版面布局，實現遙感解譯成果的自動制圖和一鍵式輸出的方法。

1 技術框架

三維環(huán)境下地質災害遙感解譯成果的輸出圖件是以遙感圖像為底圖，附加上各種可選擇的GIS圖層、圖像圖層、標注圖層等柵格或矢量圖層，再根據一定的數學法則將解譯成果以符號化的點、線、多邊形等對象繪制于底圖上，從而形成能夠準確直觀表達地質災害信息及地理分布情況的專題地圖，其主要作用是為地質災害防治和應急業(yè)務提供基礎數據。本研究以中國地質調查局遙感解譯專題圖相關制圖規(guī)范為依據，根據地質災害三維遙感解譯的技術框架［9］，采用的地質災害遙感解譯成果自動制圖輸出的技術框架如圖1所示。

圖1 制圖技術框架Fig．1 Framework ofm apping technology

根據遙感解譯成果制圖規(guī)范，最終生成的成果專題圖的地理底圖所需的遙感圖像圖層以及附加的圖像圖層主要存儲于影像數據庫中;矢量圖層(即用戶根據解譯工作需求加載的專題數據)通過空間數據庫存儲在三維GIS圖層中;地質災害遙感解譯成果數據以三維GIS圖層的形式存儲于空間數據庫中。同時根據常規(guī)制圖規(guī)范，專題圖需按照成圖比例尺和圖件類型，分別設置圖名、圖例、比例尺、指北針、圖幅號、責任表和圖廓等地圖輔助要素，其樣式和版面布局如圖2所示。

圖2 版面布局Fig．2 Layout of printed sheet

確定了專題圖各要素的存儲方式、符號化方法以及地圖排版、樣式后，可以研制開發(fā)一鍵式制圖輸出功能，實現地質災害遙感解譯成果自動制圖，生成一張滿足地質災害調查、監(jiān)測和應急需求的遙感解譯專題圖。

2 關鍵技術

2．1 三維GIS平臺及系統(tǒng)框架

本文地質災害遙感解譯成果自動制圖選用三維地理信息系統(tǒng)平臺iTelluro(網圖)作為三維環(huán)境進行研究。iTelluro是一款高效、穩(wěn)定的網絡三維GIS平臺，它采用面向Internet的分布式計算技術和三維可視化技術，支持跨區(qū)域、跨網絡的復雜大型三維應用系統(tǒng)集成［10］，為地質災害防治管理、決策和信息發(fā)布提供統(tǒng)一的解決方案和平臺技術。

iTelluro基于主流技術平臺．NET開發(fā)，產品開放性好、架構靈活、三維功能和GIS功能強大、支持TB級海量數據、多源數據(包括DEM，DOM，DLG、三維模型數據和其他專題數據)一體化管理和快速三維實時漫游功能，支持三維空間查詢、分析和運算，可與常規(guī)GIS軟件集成，提供全球范圍基礎影像資料，可方便快速地構建三維空間信息服務系統(tǒng)，亦可在二維GIS系統(tǒng)快速地完成向三維的擴展。

2．2 遙感影像底圖的獲取和繪制

系統(tǒng)解譯和制圖輸出所需的遙感圖像存儲于三維GIS平臺iTelluro的圖像數據庫中。iTelluro使用了基于哈希格網的改進橢球四叉樹方式來建立全球離散層次格網，即等邊長(度為單位)的切片剖分方式［10］。每一級切片都對應一個固定的圖像分辨率，iTelluro為應用模塊提供了WMS兼容的影像獲取接口，并可根據圖像分辨率獲取最佳分辨率的圖像數據。自動制圖模塊可根據制圖輸出的范圍、比例尺按需從遙感圖像數據庫獲取所需影像。制圖輸出成果的比例尺和圖像分辨率的關系為

式中:Res為圖像分辨率，單位是m/像素;scale為比例尺分母;dpi為每英寸的像素數，為了打印出圖的清晰，本文采用的是300 dpi;Rate為英寸和毫米的換算常數，為25．4。

根據式(1)，可由比例尺計算出所需遙感圖像的像素分辨率，從而得到對應的遙感圖像切片在iTelluro圖像數據庫中的級數，根據制圖范圍和級數可計算出對應的切片在iTelluro圖像數據庫中的索引，進而檢索出所需遙感圖像資料。如果圖像范圍過大，可從其左下角分塊從圖像數據庫返回多個圖像。獲取到遙感圖像之后，根據坐標對其進行拼接、渲染，生成遙感專題圖所需的地理底圖。這一過程是由制圖模塊自動完成的，無需人機交互。

2．3 專題圖層和解譯成果的符號化與繪制

地理和地質專題圖、地質災害遙感解譯成果，存儲在空間數據庫中，數據庫中除了遙感解譯信息表(屬性信息)外，還包括解譯目標的點、線和多邊形空間信息。從數據庫中讀取其空間對象，用于在地理底圖上繪制解譯目標;讀取其屬性信息，用于符號化和圖例、標注的繪制。基于空間分析，讀取制圖范圍內的空間對象，然后應用GDI+進行符號化、渲染和繪制。由于GDI+的繪制是以像素坐標(輸出坐標)為單位的，因此從數據庫中讀取的空間信息要先轉換成像素坐標，轉換方法根據其地理坐標和地理底圖的地理范圍的關系確定，并根據解譯成果制圖規(guī)范中的符號化規(guī)則，對解譯對象或專題對象的空間樣式進行符號化和繪制。

3 應用實例

基于iTelluro的地質災害三維遙感解譯系統(tǒng)解決了三維環(huán)境下的地質災害人機交互解譯，實現了基于空間數據庫的解譯目標空間、屬性數據一體化存儲。基于以上研究成果，本文在該系統(tǒng)平臺之上設計開發(fā)了遙感解譯成果一鍵式制圖輸出模塊。模塊主要界面如圖3所示。

圖3 模塊主要界面示意圖Fig．3 Schematic diagram of main interface of module

供選擇的出圖圖層有一般GIS圖層、本地影像圖層、本地矢量圖層和標注圖層4種。這些圖層的添加是為了豐富遙感底圖的信息，從而更好地描述專題內容。出圖的輔助要素主要有圖廓、圖名、圖件密級、標準圖幅號(非標準圖幅不需要)、指北針、比例尺、圖例和責任表等，地理底圖上需繪制經緯網。

制圖范圍確定有2種方式:自由選擇和標準圖幅。自由選擇即用戶可以直接輸入經度和緯度范圍，也可以直接在三維視圖中拉框選擇出圖范圍;標準圖幅則根據用戶輸入的坐標點和選擇的比例尺計算其所在的標準圖幅號，該標準圖幅即出圖范圍。添加的出圖圖層，可設置其透明度，以避免對基礎遙感影像的遮擋。經緯網的間隔可根據出圖范圍來自動確定，也可手工設置。

利用已有的地質災害遙感解譯數據進行制圖輸出，系統(tǒng)出圖過程和所得結果如圖4所示。

圖4 遙感解譯專題圖Fig．4 Thematic map of remote sensing interpretation

從生成的地質災害遙感解譯成果圖來看，本系統(tǒng)實現了三維環(huán)境下的地質災害遙感解譯成果一鍵化自動制圖輸出，省時省力，無需借助于第三方軟件，對地質災害遙感調查、監(jiān)測和應急業(yè)務應用具有重要的現實意義。

4 結論

本研究解決了三維環(huán)境下地質災害遙感解譯成果制圖中需借助第三方軟件進行編輯、格式轉換等繁瑣操作才能輸出成果產品的難題，實現了解譯成果根據需求，自動完成符號化、渲染和專題圖輸出，提高了地質災害遙感應用的工作效率，對促進和深化遙感技術在地質災害調查、監(jiān)測和應急救援中的應用具有一定現實意義。

［1］ 許 沖，戴福初，陳 劍，等．汶川Ms8．0地震重災區(qū)次生地質災害遙感精細解譯［J］．遙感學報，2009，13(4):754－762．Xu C，Dai FC，Chen J，et al．Identification and analysis of secondary geological hazards triggered by a Magnitude 8．0 Wenchuan Earthquake［J］．Journal of Remote Sensing，2009，13(4):754－762．

［2］ 童立強，張曉坤，程 洋，等．“8·7”甘肅舟曲縣特大泥石流災害遙感解譯與評價研究［J］．遙感信息，2011(5):109－113．Tong L Q，Zhang X K，Cheng Y，et al．Remote sensing interpretation and evaluation of super－large geological disaster caused by“8·7”debris flow of Zhouqu［J］．Remote Sensing Information，2011(5):109－113．

［3］ 屈新原，李家存，趙曉麗，等．面向災害遙感快速制圖的關鍵技術［J］．地理空間信息，2011(5):112－114．Qu X Y，Li JC，Zhao X L，et al．Critical techniques of disaster oriented remote sensing quick mapping［J］．Geospatial Information，2011(5):112－114．

［4］ 唐湘丹，馬維峰，楊強根，等．基于MAS的達州地質災害遠程監(jiān)控系統(tǒng)［J］．計算機技術與應用，2012(6):148－154．Tang X D，MaW F，Yang Q G，et al．The remote sensing system of Dazhou geological disaster based on MAS［J］．Computer Technology and Its Application，2012(6):148－154．

［5］ 烏云其其格，馬維峰，張時忠，等．基于三維的地質災害遙感解譯標志管理系統(tǒng)設計與實現［J］．國土資源遙感，2012，24(2):148－151．Wu Y Q Q G，MaW F，Zhang SZ，et al．The design and development of the interpretation marks management system for geological disasters based on 3D［J］．Remote Sensing for Land and Resources，2012，24(2):148－151．

［6］ 尹江濤，楊 斌，何政偉，等．GIS三維可視化技術在天山公路地質災害評價與決策信息系統(tǒng)中的應用研究［J］．遙感信息，2008(6):81－84．Yin JT，Yang B，He ZW，et al．Application and research of 3D GIS visualization in Tianshan highway geo－h(huán)azard appraisal and decision information system［J］．Remote Sensing Information，2008(6):81－84．

［7］ 中國地質調查局．DD2011－02—2011遙感解譯地質圖制圖規(guī)范(1∶250000)［S］．China Geological Survey．DD2011－02—2011 mapping specification of geological maps of remote sensing interpretation(1∶250000)［S］．

［8］ 中國地質調查局．DD2011－01—2011遙感影像地圖制圖規(guī)范(1 ∶50000、1 ∶250000)［S］．China Geological Survey．DD2011－01—2011 mapping specification of maps of remote sensing interpretation(1 ∶50000、1 ∶250000)［S］．

［9］ 唐湘丹，馬維峰．基于Virtual Globes的地質災害遙感解譯應用框架［C］//Proceedings of 2010 International Conference on Remote Sensing(ICRS 2010)Volume 4．2010．Tang X D，MaW F．Application framework for remote sensing interpretation of geological disaster based on Virtual Globes［C］//Proceedings of 2010 International Conference on Remote Sensing(ICRS 2010)Volume 4．2010．

［10］ 武漢地大信息工程有限公司．iTelluro Globe Engine SDK使用說明書 2．5［Z］．武漢:［s．n．］，2011．Wuhan Infoearth Information Co．，Ltd．The manual of iTelluro Globe Engine SDK［Z］．Wuhan:［s．n．］，2011．