應(yīng)急響應(yīng)中PixelGrid無(wú)人機(jī)遙感數(shù)據(jù)處理的關(guān)鍵技術(shù)與應(yīng)用

賈 嬌，艾海濱，張 力，燕 琴

(1.蘭州交通大學(xué)，甘肅 蘭州730070;2.中國(guó)測(cè)繪科學(xué)研究院，北京100830)

一、引 言

現(xiàn)今，在應(yīng)急響應(yīng)和防災(zāi)減災(zāi)領(lǐng)域，通過(guò)低空無(wú)人機(jī)獲取災(zāi)區(qū)影像的方式發(fā)揮的作用越來(lái)越明顯，其最大的優(yōu)勢(shì)就是可以在惡劣的自然條件下，或者在人無(wú)法到達(dá)的地區(qū)，利用簡(jiǎn)單條件就可以低空飛行，獲取影像數(shù)據(jù)。同時(shí)，其也存在一個(gè)不容忽視的問(wèn)題，即與傳統(tǒng)的航空拍攝相比，應(yīng)用小型無(wú)人機(jī)獲取數(shù)據(jù)的方式受到氣候因素、重量因素、相機(jī)質(zhì)量因素等影響，導(dǎo)致飛機(jī)姿態(tài)不穩(wěn)定，從而使得通過(guò)無(wú)人機(jī)獲得的影像質(zhì)量會(huì)有一些不足，如影像航向重疊度和旁向重疊度都不夠規(guī)則、像幅較小、像片數(shù)量多、影像的傾角過(guò)大且傾斜方向沒(méi)有規(guī)律、航攝區(qū)域地形起伏大造成的影像旋偏角大、影像有明顯畸變等。這些情況下實(shí)現(xiàn)自動(dòng)空三是現(xiàn)有數(shù)字?jǐn)z影測(cè)量系統(tǒng)的主要挑戰(zhàn)，在大多數(shù)情況下都將導(dǎo)致錯(cuò)誤結(jié)果。

因此，完善的影像數(shù)據(jù)處理系統(tǒng)應(yīng)該既能快速處理數(shù)量龐大的無(wú)人機(jī)影像，又能保證處理結(jié)果準(zhǔn)確且滿足既定的要求。本文以影像數(shù)據(jù)測(cè)圖系統(tǒng)PixelGrid為平臺(tái)，介紹了無(wú)人機(jī)影像數(shù)據(jù)快速處理的關(guān)鍵算法與相關(guān)技術(shù)，以及PixelGrid在應(yīng)急響應(yīng)處理中的應(yīng)用。

二、攝影測(cè)量快速處理

針對(duì)無(wú)人機(jī)影像的特點(diǎn)對(duì)其進(jìn)行數(shù)據(jù)快速處理的流程如圖1所示。

圖1 無(wú)人機(jī)影像快速處理流程

獲得無(wú)人機(jī)數(shù)據(jù)以后，首先對(duì)影像作旋轉(zhuǎn)、主點(diǎn)修正、畸變改正或格式轉(zhuǎn)換等預(yù)處理;然后再進(jìn)行自動(dòng)相對(duì)定向、模型連接、航帶間轉(zhuǎn)點(diǎn)等完成空中三角測(cè)量;最后利用特征提取技術(shù)從影像中提取數(shù)字表面模型(DSM)，DSM經(jīng)濾波處理得到DEM。下面對(duì)原始影像作勻光勻色處理，再結(jié)合生成的DEM對(duì)其進(jìn)行數(shù)字微分糾正，然后利用正射影像自動(dòng)拼接技術(shù)得到整個(gè)測(cè)區(qū)的正射影像圖。

三、無(wú)人機(jī)影像處理的關(guān)鍵算法

傳統(tǒng)攝影測(cè)量軟件在進(jìn)行相對(duì)定向之前要設(shè)定影像的航向重疊度，而低空無(wú)人機(jī)拍攝出的影像航向旁向重疊度都不規(guī)則，致使在金字塔影像最高層的搜索范圍不正確。無(wú)人機(jī)影像也存在旋偏角過(guò)大的問(wèn)題，傳統(tǒng)的灰度匹配算法會(huì)使影像不相關(guān)。以上兩個(gè)弊端導(dǎo)致了無(wú)人機(jī)影像的相對(duì)定向結(jié)果不理想。PixelGrid系統(tǒng)利用了SIFT算法進(jìn)行匹配，并基于尺度/旋轉(zhuǎn)不變性進(jìn)行全自動(dòng)相對(duì)定向［1］，圖2展示了在無(wú)人機(jī)影像旋偏角較大、紋理缺乏的情況下得到296個(gè)匹配點(diǎn)，σ值為1.12個(gè)像素，說(shuō)明此方案的匹配成功率和匹配點(diǎn)精度都較高，能滿足實(shí)際生產(chǎn)需求;得到相對(duì)定向元素以后，進(jìn)行模型連接，并構(gòu)成自由航帶網(wǎng)［2］。無(wú)人飛行器存在影像畸變大、飛行狀態(tài)不穩(wěn)定等不足，致使模型連接的中誤差超出限差過(guò)多。為了解決這個(gè)問(wèn)題，首先在影像重疊區(qū)域提取F?rstner特征，經(jīng)過(guò)灰度匹配獲得初始的模型連接點(diǎn)，再應(yīng)用直方圖統(tǒng)計(jì)法、模型連接比定義、幾何相似關(guān)系法逐步去除匹配錯(cuò)誤的點(diǎn)。此方法獲得的連接點(diǎn)分布均勻，轉(zhuǎn)點(diǎn)的成功率和可靠性高，為后續(xù)的航帶間轉(zhuǎn)點(diǎn)和區(qū)域網(wǎng)平差提供了保障。

完成自動(dòng)空三解算后，利用原始影像和內(nèi)外方位元素生成數(shù)字表面模型(DSM)或DEM。從遙感影像中自動(dòng)提取三維信息的方法很多，它們能滿足普遍的生產(chǎn)要求，但弊端是精度和可靠性不高［3］。針對(duì)受災(zāi)地區(qū)的復(fù)雜地況和無(wú)人機(jī)的非常規(guī)拍攝，可采用由粗到細(xì)(coarse-to-fine)的多級(jí)影像匹配策略［4］，融合特征點(diǎn)、特征線和格網(wǎng)點(diǎn)，集成多種成熟的、性能互補(bǔ)的影像匹配算法(基于面元的匹配、基于特征的匹配和關(guān)系匹配等)，同時(shí)自動(dòng)進(jìn)行匹配粗差定位和剔除，以保證匹配精度和可靠性。其核心算法為基于物方的多影像匹配方法［5］，可以同時(shí)匹配多于2景的影像，可有效地解決重復(fù)紋理或遮擋區(qū)域的匹配困難。同時(shí)采用兼顧地形特征線信息的全局概率松弛法影像匹配算法［5-6］，兼顧影像的全局和局部信息(如圖3所示)。

圖3 基于多基線、多匹配特征的地形信息自動(dòng)提取

由于拍攝相機(jī)和外部環(huán)境的影響，拼接后的影像容易出現(xiàn)色彩或灰度的差別。基于影像直方圖匹配的勻色算法［7］將縮小一定程度拼接后的大影像作為參考模板，使原大影像灰度直方圖變成規(guī)定形狀的直方圖而對(duì)圖像色彩作修正，這樣使得拼接的色調(diào)過(guò)渡平滑，也極大地提高了自動(dòng)化程度。另外，采用適當(dāng)調(diào)節(jié)影像的亮度和對(duì)比度屬性的方法［7］來(lái)去除影像上的大面積霧靄，增強(qiáng)影像對(duì)比度，彌補(bǔ)了以往的勻光算法(如自適應(yīng)模板法)［7］在此方面的缺陷。

為得到精確的DEM，要對(duì)色彩調(diào)整后的原始影像進(jìn)行數(shù)字微分糾正以生成正射影像圖。正射影像拼接需要做到無(wú)縫鑲嵌、無(wú)影像錯(cuò)位，清晰準(zhǔn)確地呈現(xiàn)測(cè)區(qū)的實(shí)際情況。可運(yùn)用動(dòng)態(tài)規(guī)劃算法［8］，根據(jù)多影像的整理分布情況自動(dòng)生成區(qū)域范圍內(nèi)的接縫線網(wǎng)絡(luò)，進(jìn)行整體接縫線自動(dòng)編輯，同時(shí)對(duì)接縫線進(jìn)行實(shí)時(shí)的羽化處理，以保證鑲嵌處理的效率。生成的拼接線能夠繞開(kāi)房屋、樹(shù)木、橋梁等，而且一般不需要人工手動(dòng)修改。

四、高效能處理技術(shù)

數(shù)字?jǐn)z影測(cè)量軟件可處理的數(shù)據(jù)量越來(lái)越大，在應(yīng)急響應(yīng)中越快獲知災(zāi)區(qū)情況越好，為了提高效率，并行處理方式在數(shù)據(jù)處理中廣泛應(yīng)用起來(lái)。并行處理是計(jì)算機(jī)系統(tǒng)中同時(shí)執(zhí)行多個(gè)處理機(jī)的一種計(jì)算方法，目的就是快速解決大型且復(fù)雜的計(jì)算問(wèn)題。單機(jī)多核是最簡(jiǎn)單的一種并行方式，即一臺(tái)計(jì)算機(jī)擁有多個(gè)內(nèi)核，允許多個(gè)進(jìn)程在一臺(tái)機(jī)器上同時(shí)進(jìn)行處理。經(jīng)測(cè)試驗(yàn)證，一臺(tái)8核PC機(jī)只需75 min即可完成856個(gè)像對(duì)UCXP航空影像的自動(dòng)相對(duì)定向，200個(gè)像對(duì)用單機(jī)多核的方式只需3 h即可完成DSM的自動(dòng)提取。

為了獲得更優(yōu)質(zhì)的效率，還可以采取多機(jī)集群分布式的數(shù)據(jù)處理方法，即將局域網(wǎng)中互聯(lián)的所有工作站(包括PC機(jī)和高性能集群計(jì)算機(jī))通過(guò)軟件的方式進(jìn)行通信和協(xié)作，從而以一定的任務(wù)調(diào)度策略共同完成影像數(shù)據(jù)的分布式處理［9］。圖4展示了PixelGrid系統(tǒng)的分布式處理架構(gòu)，基于高速局域網(wǎng)絡(luò)，基本上實(shí)現(xiàn)了基于松散耦合并行服務(wù)中間件的影像數(shù)據(jù)集群分布式并行計(jì)算［10］。當(dāng)接到任務(wù)時(shí)，任務(wù)調(diào)度模塊首先根據(jù)性能檢測(cè)模塊的報(bào)告，按照負(fù)載均衡的方式將待處理的任務(wù)發(fā)送到相應(yīng)的處理節(jié)點(diǎn)上;然后操作員通過(guò)軟件界面實(shí)時(shí)了解任務(wù)進(jìn)展，接收遠(yuǎn)端處理完的成果數(shù)據(jù)，并在本機(jī)上儲(chǔ)存該數(shù)據(jù)。影像正射糾正中涉及的重采樣操作往往耗時(shí)巨大，在高數(shù)據(jù)處理量的系統(tǒng)中，采用分布式計(jì)算，將系統(tǒng)計(jì)算功能分塊并行計(jì)算，可大大提高處理效率。

圖4 PixelGrid系統(tǒng)分布式處理架構(gòu)

在實(shí)際生產(chǎn)中，可以根據(jù)需要選擇使用單機(jī)多核或者多機(jī)分布式處理。表1以匹配DSM為例比較了兩種處理方式的效率，具體處理內(nèi)容是從36張8 cm分辨率的UCD航空影像生成格網(wǎng)間隔為40 cm的DSM。由表1可知多機(jī)分布式并行處理比單機(jī)處理更能提高效率，減輕了人員的工作量，實(shí)現(xiàn)了作業(yè)步驟的高度自動(dòng)化，為災(zāi)害的應(yīng)急處理爭(zhēng)取了時(shí)間。

表1 PixelGrid兩種處理方式在匹配DSM中的效率對(duì)比

GPU技術(shù)以其卓越的圖形處理功能，在數(shù)字?jǐn)z影測(cè)量領(lǐng)域的應(yīng)用越來(lái)越重要，根據(jù)其并行結(jié)構(gòu)和硬件特點(diǎn)，使利用GPU實(shí)現(xiàn)通用計(jì)算和圖像處理的高性能并行計(jì)算成為可能，并且發(fā)展成為趨勢(shì)。對(duì)于一幅分辨率為0.2 m、大小為756 MB的正射影像，采用CPU糾正大約需要8 min，而基于FX1800顯卡的GPU精糾正處理只需30 s，加速比為16倍。由于正射糾正是針對(duì)單個(gè)像元的操作，隨著影像像元數(shù)的增多，計(jì)算量也急劇增加。PixelGrid系統(tǒng)利用OpenGL作為手段，修改三維繪制管線的可編程Vertex和 Fragment Shader來(lái)實(shí)現(xiàn)快速的糾正［11］。以21張UCD航空遙感影像(測(cè)區(qū)面積約73 km2，數(shù)據(jù)量9.2 GB)為例，利用單機(jī)單線程GPU處理模式生成分辨率為0.2 m的正射影像，在NVIDIA Quadro 4000顯卡和1000轉(zhuǎn)速SAS盤(pán)的機(jī)器上運(yùn)行，只需要約8 min。通過(guò)GPU并行正射糾正的實(shí)例說(shuō)明，

綜上所述，PixelGrid在玉樹(shù)地震后僅用了7個(gè)小時(shí)即完成了從自動(dòng)空三到生成正射影像的全部流程，制作出玉樹(shù)結(jié)古鎮(zhèn)震后0.4 m分辨率的正射影像圖，圖5的正射影像顯示出全鎮(zhèn)的受災(zāi)情況，圖6是結(jié)古鎮(zhèn)震后2 m分辨率DSM彩色暈渲圖。另外，對(duì)舟曲泥石流后近200幅無(wú)人機(jī)航空影像數(shù)據(jù)進(jìn)行處理，從自動(dòng)空三到正射影像制作完成共用近6 h，獲得了0.2 m分辨率的舟曲縣泥石流災(zāi)害影像圖，如圖7所示，能顯示泥石流經(jīng)過(guò)地區(qū)的受災(zāi)狀況。這些成果可以幫助救災(zāi)人員判斷受災(zāi)的具體區(qū)域和嚴(yán)重程度，為調(diào)派人力物力前往救助提供保障。另外，如圖8所示，正射影像圖也能為監(jiān)測(cè)災(zāi)后地質(zhì)形態(tài)(滑坡群、堰塞湖等)變化、關(guān)注受災(zāi)后城鎮(zhèn)的恢復(fù)與重建提供保障。利用GPU對(duì)攝影測(cè)量中相關(guān)圖像處理算法的并行化，可以極大地提高攝影測(cè)量處理的效率。

圖5 結(jié)古鎮(zhèn)震后正射影像圖

圖6 結(jié)古鎮(zhèn)震后DSM彩色暈渲圖

由影像自動(dòng)匹配生成的DSM包含地面建筑物、樹(shù)木和橋梁等高度，需要轉(zhuǎn)換為DEM以供后續(xù)流程應(yīng)用。一般是靠單人手工編輯，耗費(fèi)人力和時(shí)間巨大。為了提高效率，運(yùn)用多人網(wǎng)絡(luò)分布式DSM在線編輯系統(tǒng)，即由服務(wù)器端將一片測(cè)區(qū)分成若干塊區(qū)，若干客戶端在網(wǎng)絡(luò)內(nèi)同時(shí)進(jìn)行DSM編輯，處理時(shí)間可成倍減少，為后續(xù)流程的順利進(jìn)行提供了方便。

五、應(yīng) 用

圖7 舟曲縣泥石流正射影像圖

圖8 堰塞湖正射影像圖

六、結(jié)束語(yǔ)

本文基于PixelGrid系統(tǒng)的數(shù)據(jù)處理流程和方法，針對(duì)低空無(wú)人機(jī)拍攝方式的現(xiàn)狀，介紹了遙感影像快速處理技術(shù)的關(guān)鍵算法，詳細(xì)闡述了并行計(jì)算技術(shù)與數(shù)字?jǐn)z影測(cè)量相關(guān)算法相結(jié)合的模式，并通過(guò)PixelGrid在玉樹(shù)地震和舟曲泥石流的災(zāi)后影像處理中的應(yīng)用，驗(yàn)證了這兩種技術(shù)相結(jié)合可以充分發(fā)揮應(yīng)急響應(yīng)中數(shù)據(jù)處理的速度和效率。

［1］王琳，張力，艾海濱.無(wú)人飛行器遙感影像的自動(dòng)相對(duì)定向研究［J］.測(cè)繪通報(bào)，2011(11)：11-14.

［2］王琳.高精度、高可靠的無(wú)人機(jī)影像全自動(dòng)相對(duì)定向及模型連接研究［D］.北京：中國(guó)測(cè)繪科學(xué)研究院，2011.

［3］張力，張繼賢.基于多基線影像匹配的高分辨率遙感影像DEM 的自動(dòng)生成［J］.測(cè)繪科學(xué)，2008，33(S2)：35-39.

［4］袁修孝，明洋.一種綜合利用像方和物方信息的多影像匹配方法［J］.測(cè)繪學(xué)報(bào)，2009，38(3)：216-222.

［5］張祖勛，張劍清.數(shù)字?jǐn)z影測(cè)量學(xué)［M］.武漢：武漢大學(xué)出版社，2010：148-213.

［6］ZHANG Li，GRUEN A.Multi-image Matching for DSM Generation from IKONOSImagery［J］.ISPRSJournal of Photogrammetry and Remote Sensing，2006，60(3)：195-211.

［7］王邦松，艾海濱，安宏，等.航空影像色彩一致性處理算法研究［J］.遙感信息，2011(1)：45-49.

［8］CHON J，KIM H.Determination of the Optimal Seamlines in Image Mosaicking with the Dynamic Programming(DP)on the Converted Cost Space［J］.Lecture Notes in Computer Science，2006(4029)：750-757.

［9］王彥敏，盧剛.基于PixelGrid實(shí)現(xiàn)DOM的快速更新［C］∥地理信息與物聯(lián)網(wǎng)論壇暨江蘇省測(cè)繪學(xué)會(huì)2010年學(xué)術(shù)年會(huì).無(wú)錫：江蘇省測(cè)繪學(xué)會(huì)，2010.

［10］艾海濱，張劍清.基于中間件的高分辨率衛(wèi)星影像正射糾正的分布式處理［J］.測(cè)繪科學(xué)，2009，34(4)：158-160.

［11］韓元利.基于GPU編程的虛擬自然環(huán)境技術(shù)研究［D］.武漢：武漢大學(xué)，2007.