馬 娟,張鳴之,馮 振,黃 喆,薛躍明,石愛軍,邵 海
(中國地質環境監測院,北京 100081)
我國是地質環境脆弱、地質災害多發的國家之一,滑坡、崩塌、泥石流、地面塌陷等地質災害易發區面積約占國土面積的65%。國土資源部從1999年開始,在地質災害嚴重的縣(市),陸續部署開展了縣市地質災害調查與區劃工作[1],截至2017年,完成1 502個縣(市、區)的1∶5萬地質災害詳細調查和29 942處隱患點的勘查。我國共查明地質災害287 916處,其中,滑坡142 377處,崩塌68 419處,泥石流33 524處,地面塌陷、地裂縫、地面沉降等43 596處[2]。傳統的地質災害野外調查手段是使用羅盤、地質錘等工具,基于紙質地形圖進行野外數據采集,然后室內整理手工編錄等[3-5],這種工作方式已經不能滿足日益增長的快速獲取信息的需求,以及國家精細化地質調查的要求,地質災害調查的精度、廣度、深度有待于進一步提高。隨著信息技術的發展,掌上機、無人機、InSAR、高精度遙感等技術在地質災害領域已有了不同程度的應用[6-8],但也存在一些問題。以無人機遙感技術為例:一是地質災害調查人員通常為地質專業出身,學習無人機操控技術難度較大;二是市場無人機操作培訓價格高,難以在地質調查隊伍中普及;三是無人機設備搭載多鏡頭傾斜相機的價格較高,難以在野外地質調查小組中形成標配推廣。因此亟需設計研發一套適用于我國地質災害野外調查的無人機遙感調查技術方法。
筆者所在的團隊通過長期跟蹤地質調查隊伍,開發出一套地質災害調查微型無人機航線控制系統,用于解決無人機技術在地質災害調查中運用的普適性問題。系統基于市面現有成型微型無人機硬件技術,通過開發適用于地質災害調查無人機航線規劃及拍攝模式,使用無人機快速、及時獲取地質災害隱患點影像數據。通過數據處理軟件,生成地質災害隱患點及其周圍環境的正射影像數據、數字高程數據以及三維模型等,有利于地質災害調查人員獲取地質災害隱患點及其周邊精準數據。
根據無人機數字航空攝影測量與遙感外業技術規范要求,地質災害調查無人機數據采集流程[9-10]見圖1。

圖1 無人機數據采集流程Fig.1 Process of UAV data acquisition
資料收集:收集工作區地形圖、前期遙感影像及解譯圖和區域基本資料等,主要目的是為了解地質災害調查區域相對高差數據、周邊環境概況,以確保采集作業安全。如果調查區域缺乏以上資料,建議第二步現場踏勘時手動飛行至調查區域上空進行概況調查。
現場踏勘:依據收集的資料確定無人機起降場地及地面控制人員安全有效的控制位置,預規劃航線同時進行GPS信號檢測。
航線規劃:按照實際需要的地面分辨率進行設計,根據無人機的性能參數以及地質災害調查區域的地形特征進行航線規劃,航線能完整覆蓋整個航拍區域。
飛行作業:飛機組裝,檢查飛機狀態指標,按既定航線規劃進行飛行作業。
控制點測量:根據測區實際狀況和航線規劃情況進行控制點布設及測量。
影像下載:連接無人機的儲存設備到電腦,導出影像數據。
質量檢查:對獲取的測區影像進行質量檢查,確定影像質量是否有云、霧、雪,是否有變形、扭曲、發虛等現象,判斷是否需要補飛或重飛。
外業成果:經影像質量檢查無誤后形成的無人機數據采集結果。
面向地質災害調查業務需求以及地質災害無人機數據采集流程,設計微型無人機航線控制系統的總體架構(圖2),總體上包括支撐層、通信層、數據層、應用層以及UI層。系統基于無人機飛行器、遙控器以及移動設備等支撐層,通過通信手段,采集地質災害區域的影像數據;基于數據庫提供的同一數據模型和數據服務,構建微型無人機航線控制系統應用模塊;通過UI層實現各個應用模塊的集成應用,為地質災害調查用戶提供微型無人機航線控制服務。

圖2 系統總體架構Fig.2 Diagram of the system framework
(1)支撐層
微型無人機航線控制系統支撐層主要包括無人機飛行器、遙控器、移動設備等硬件支撐設備。
(2)通信層
微型無人機航線控制系統通信層主要包括無人機飛行器、遙控器、移動設備、服務器之間的相互通信方式。其中,移動設備與遙控器通過USB通信協議連接,遙控器與飛行器通過無線電通信連接,移動設備與服務器通過網絡通信連接。通過多種通信方式,實現了地質災害數據在多平臺之間的傳輸。
(3)數據層
無人機地質災害調查數據庫面向地質災害調查業務應用的建設需求,構建統一的數據存儲、管理、應用和服務平臺,兼容多源、多尺度海量數據,實現各類地質災害調查業務數據的一體化存儲、管理和服務。
(4)應用層
業務應用層構建于無人機地質災害調查數據庫之上,提供面向地質災害調查業務應用、管理的一體化信息服務。業務應用主要功能包括:起飛條件檢查功能、飛行任務管理功能、航線規劃功能、任務模式功能、智能續飛功能、圖像快拼功能以及離線地圖功能等。
(5)UI層
UI層通過數據集成、業務集成和成果集成,實現各個應用模塊的集成應用,為用戶提供微型無人機航線控制服務。
微型無人機航線控制系統以手機或者平板電腦為終端,以Android Studio2.1.2和JDK1.8.0為開發平臺,使用SQLite和MongoDB數據庫,Java語言進行開發,運行環境Android4.4以上。

圖3 系統界面設計圖Fig.3 Design of the system interface
2.3.1起飛條件檢查
全面自動檢查無人機起飛條件,包括無人機與地面站(即微型無人機航線控制系統)連接是否正常、無人機電量是否充足(不低于30%)、GPS衛星數是否足夠(衛星數不小于6)、內存卡是否正常插入、航線是否繪制、設備是否已激活,只有當起飛條件都滿足時,才能一鍵起飛并執行后續的任務。

圖4 起飛檢查界面Fig.4 Interface of take-off check
2.3.2飛行任務管理
飛行任務管理分為已完成任務管理和未完成任務管理兩類。已完成任務指已經完成的飛行任務;未完成任務指已規劃完航線,但尚未進行飛行的任務。用戶可以進行飛行任務查看、命名、刪除以及數據下載和上傳等操作。數據下載指的是自動從無人機端下載無人機影像數據到手機端或平板電腦端,數據上傳指的是自動將已下載的無人機影像數據上傳到后方無人機數據處理平臺。
2.3.3航線規劃
航線規劃主要指設置航線參數,包括飛行高度、飛行速度(分為高、中、低三檔)、旁向(航向)重疊率(分為高、中、低三檔)、航線框類型等。在地理底圖上通過手勢拉框圈定范圍,通過飛行高度、數據重疊率、飛行區域、飛行速度等參數,自動計算生成無人機飛行航線及拍攝點(航線上灰色點即為拍攝點)。

圖5 航線規劃界面Fig.5 Route planning interface

參數可選范圍飛行速度低(1 m/s)、中(5 m/s)、高(10 m/s)旁向重疊率橫向:低(75%)、中(85%)、高(95%)航向重疊率縱向:低(60%)、中(70%)、高(80%)區域框矩形、圓形、自定義、橫線、折線
2.3.4任務模式
根據《滑坡崩塌泥石流災害調查規范(1∶50 000)》,滑坡、崩塌、泥石流等地質災害點需調查其分布范圍、規模、結構特征等。針對此要求,開發了任務模式。任務模式是本系統的核心功能,通過算法控制無人機自動旋轉機身以及自動調節鏡頭角度,實現單個鏡頭模擬多鏡頭傾斜相機拍攝效果,準確獲取地質災害及其發育環境要素信息,后期計算并生成地質災害三維模型,從而確定滑坡、崩塌、泥石流和不穩定斜坡的類型、規模及空間分布特征。任務模式共開發了4種適用于地質災害調查的數據采集方法:傾斜模式、正視模式、快速正射模式和全景模式。
傾斜模式適合坡度小于45°的斜坡,鏡頭垂直于水平面,通過旋轉機身及照相機鏡頭角度,每個拍攝點從1個垂直、2個傾斜3個角度拍攝3張照片:首先鏡頭垂直向下拍攝1張;接著機身順時針旋轉90°,鏡頭上仰45°拍攝1張;保持仰角不變,機身繼續順時針旋轉180°拍攝1張。
正視模式適合坡度大于45°的陡坡,鏡頭平行于水平面,通過旋轉機身及照相機鏡頭角度,每個拍攝點也從3個角度拍攝3張照片:首先鏡頭水平正對拍攝地物拍攝1張;然后機身逆時針旋轉45°拍攝1張;再順時針旋轉90°拍攝1張。
快速正射模式適用于地質災害應急調查,通過控制相機的曝光時間,鏡頭垂直于水平面定時拍照,實現無人機在執行航線飛行過程中不懸停拍照,極大縮短了調查時間。
全景模式適合地質災害點周圍概況調查,飛機飛到指定點上空,鏡頭分別呈水平方向、斜向下與水平面夾角35°、斜向下與水平面夾角70°,機身旋轉360°,轉3圈,每45°拍攝1張,最后鏡頭垂直向下拍攝1張,共拍照25張。
如遇陡緩結合型,比如后緣殘留危巖體與前緣崩滑(塌)堆積體構成的災害體,聯合正視模式、快速正射模式兩種模式進行拍攝,陡崖部分使用正視模式拍攝,斜坡崩塌堆積體使用傾斜模式拍攝,拍攝完成后進行聯合建模(圖6)。

圖6 不同地形任務模式示例Fig.6 Examples of different terrain mission modes
2.3.5智能續飛
針對微型無人機供電電池容量小、飛行時間短的問題,微型無人機航線控制系統提供智能續飛功能:后一架次的起點自動設為航線未完成部分的第一個點,依此類推,從而滿足較大區域連續拍攝需求。
2.3.6圖像快拼
當任務模式為快速正射模式時,數據采集任務完成后,APP端通過遙控器與無人機連接全自動下載該任務所有拍攝照片,并對下載照片進行自動拼接,快速獲得現場正射影像圖。

圖7 圖像快拼示意圖Fig.7 Diagram of image quick spelling
2.3.7離線地圖
離線地圖下載可以框選所需調查區域進行地圖下載,分為正在下載和已完成兩類,可實時顯示下載進度。用戶可對地圖下載任務進行暫停、繼續、刪除、查看等操作。
試驗對象選擇伊寧縣喀拉亞尕奇鄉潘津村皮里青河3號滑坡。實地調查發現該點位于皮里青河主溝右岸斜坡,屬山前黃土丘陵,坡體裸露,植被覆蓋率低其剖面形態呈凹形,平面呈不規則,坡度約40°。斜坡地表出露第四系風積黃土(Qeol),土體內部垂直節理裂隙較發育,下伏石炭系伊什基里克組砂礫巖(C2y)。斜坡坡腳為河流階地,主河道距離坡址約55 m。滑坡滑動地層為黃土,滑體厚約12 m,滑動后在坡腳堆積。
本試驗使用大疆公司的精3專業版(Phantom-3-Pro)無人機,其具有重量輕、成本低、靈活、便攜、快速等優勢。實地踏勘該滑坡后,認為該滑坡地形使用傾斜攝影任務模式較為適合,規劃航線后,使用續飛功能飛行2個架次,完成該滑坡無人機調查,共取得照片84張,設置飛行高度260 m,飛行速度高,旁向重疊率中,航向重疊率中。
飛行獲得的數據通過上傳自主研發的iFlyCloud無人機三維數據處理云平臺(http://fly.cigem.cn),生成正射影像圖、三維模型以及數字表面模型等一系列數據產品。從生成的數據產品質量可以看出,微型無人機航線控制系統所采用的數據采集方法完全可行。
圖8為制作生成的皮里青河3號滑坡三維模型,從圖中可以明顯劃分滑坡的滑源區、滑動區、灑落區以及堆積區,同時可通過高程變化判別出滑坡的滑動方向。通過自主開發iFlyCloudViewer三維模型瀏覽軟件可快速測量出整個滑坡的寬度67 m,長度340 m,面積2.2×104m2,體積117.8×104m3,坡度為40°。滑坡堆積體延伸至主河道中間,前緣被河流沖刷侵蝕形成深約5.2 m的陡坎。滑坡后緣可見明顯的陡坎和拉張裂縫,陡坎高2~3 m,大型拉張裂縫約有4條,最長的達到65 m。滑坡堆積體松散,后緣陡坎形成新的臨空面,在降雨、冰雪融水,河流沖刷等作用下,易發生局部滑動,滑坡范圍將向山頂進一步擴展。

圖8 三維立體解譯(攝于2017年7月19日)Fig.8 3-D interpretation(taken on July 19, 2017)
目前,微型無人機的硬件技術已經成熟,但缺乏自動控制微型無人機對地質災害調查進行數據采集的航線控制系統軟件。本項目基于Android系統設計了適用于地質災害野外調查微型無人機航線控制系統,以微型無人機為調查工具,實現了地質災害野外數據采集流程全自動化,極大降低了地質災害調查人員的工作強度,提高了地質災害調查數據采集的效率。
此系統能夠廣泛運用于地質災害日常調查以及突發性地質災害的應急調查中,特別是4種無人機遙感數據采集任務模式能夠快速科學規劃航線以適應斜坡、陡坡、陡緩結合型等不同地形需求,采集獲取地質災害及其發育環境要素遙感影像數據。通過三維建模軟件平臺計算,生成二、三維數據產品,在產品的基礎上,通過遙感信息提取技術與地質災害分析相結合,提供研究地質災害調查點的基礎數據,對提高地質災害調查的精度有著重要意義。
致謝:武漢地大信息工程有限公司提供了相關材料和技術實現,在此表示感謝!