機載激光3D探測成像系統的關鍵技術

孟慶季，張續嚴，周 凌，王 超

(中國科學院長春光學精密機械與物理研究所，吉林長春130033)

1 引言

機載激光3D探測成像系統集光、機、電技術于一體，從組成上可分為激光測距系統、陀螺穩定平臺、GPS/INS(全球定位/慣性導航)組合導航系統、數據處理與3D顯示系統等。對比傳統探測載荷系統，機載激光3D探測成像系統具有以下優點：(1)主動照明，具備全天候能力;(2)電磁干擾能力強，對背景有極強抑制能力，不易受環境溫度及陽光的影響;(3)抗隱身能力強，能穿透一定的遮蔽物、偽裝和掩體;(4)具有高的距離、角度和速度分辨率，能同時獲得目標的多種圖像(如距離像、距離-角度像等)，圖像信息量豐富，易于目標識別等。

機載激光3D探測成像系統能夠快速精確地獲得地面的3D測量信息，與有人/無人機平臺相結合在地形測繪、森林調查、夜間軍事偵察、夜間搜尋與救援、軍事偽裝識別和水下目標探測諸多領域有著廣泛的應用。

2 國內外發展現狀

2.1 國外現狀

國外發達國家開展激光成像探測技術的研究較早，目前已經取得了相對豐碩的成果。

機載激光3D探測成像系統伴隨著激光的問世一同發展。美國麻省理工學院(MIT)的林肯實驗室是世界上激光探測成像領域比較前沿的研究中心，其主要從事激光3D成像系統的關鍵技術及系統技術研究，2004年已經成功研制出32×32焦平面成像陣列(FPA)激光3D成像系統，該系統采用固體光纖泵浦激光器(530和780 nm)，激光測距頻率為8～10 kHz，采取了掃描工作方式來擴大視場角［1］。

美國先進科技公司(Advanced Scientific Concepts Inc.)也公布了激光3D成像雷達的原型樣機(圖1)，該樣機采用128×128 InGaAs PIN陣列來實現激光飛行時間的測量［2］。

瑞典的CSEM開發了另外一種激光3D成像系統FPA傳感器，并致力于開發一種集成、低功耗、小型化的實時激光3D成像系統，其具有30 Hz的成像幀頻，124 pixel×160 pixel，每個像元均可以通過測量信號(激光調制產生)的相位差來得到激光的傳輸時間。目前該系統主要應用于無人車輛成像，測量距離較近［3］。

圖1 激光3D成像雷達的原型樣機Fig.1 Original prototype of laser 3D imaging system

另外，美國雷聲公司和洛克希德馬丁公司也都致力于下一代無人車激光3D成像雷達的研究，其采用的面陣雪崩光電二極管(APD)探測器等單元技術均極大地推進了激光3D成像雷達的發展［4］。

機載激光3D成像系統不但包括激光3D成像技術，還包括INS/GPS組合導航、陀螺穩定平臺、數據的處理與3D顯示等技術，其涉及學科更廣、技術難度更大。下面介紹幾個國外最新的機載激光3D成像系統：

(1)小型無人機載激光3D探測成像傳感(JIGSAW)系統

2007年資料報道的JIGSAW系統為美國高級防御研究計劃局(Defense Advanced Research Projects Agency，DARPA)和美國陸軍夜視和電子傳感器管理局(The U.S.Army Redcomcerdec Night Vision and Electronic Sensors Directorate，NVESD)研制的一種集成、輕型的，可以用來對植物遮蔽下的隱秘目標進行成像的激光3D成像系統，旨在得到更高分辨率的成像及對遮蔽物體的識別能力，計劃裝備于 DP-5X 型無人飛機［5，6］。系統基本結構如圖2所示。

JIGSAW系統可分為平臺上系統與平臺下系統。30 cm(12 in)的穩定平臺(內含光子計數直接探測型激光3D成像系統、傳感器控制電路及熱控等)采用4軸機構，結合INS/GPS組合導航系統實現激光3D成像系統視軸的高精度穩定、指向控制與感知。平臺下系統包括系統控制器、電源、高速數據處理器和高數數據傳輸系統(3D圖像向地面控制站的實時傳輸等)。

圖2 JIGSAW系統結構圖Fig.2 Structure diagram of JIGSAW system

(2)機載激光測深系統(LADS)

2008年的JANE防務周刊報道了機載激光測深系統(LADS)。LADS安裝在固定翼飛行器上用于探測深海，初期為澳大利亞政府的防務科學與技術部門設計，早于1998年已服役。系統參數如表1所示，其基本組成可分為［7］：

表1 系統參數Tab．1 System parameter

1)穩定激光平臺(Stabilized laser platform)：900 Hz的全固態紅外激光器安裝在穩定平臺上，以便隔離機體的俯仰、橫滾、航向的姿態運動。

2)機載捷聯慣導系統：為系統提供重要的位置、姿態基準信息。

3)機載激光控制與數據獲取軟件與設備：實時記錄激光測距信息、精確的位置信息與姿態信息(INS/GPS)、穩定平臺的角度信息等，為地面設備的后處理完成測量數據的收集。

4)數據處理軟件和設備。

(3)直升機3D-LZ Imaging LADAR

表2 系統參數Tab．2 System parameter

2010年的最新資料報道，直升機3D-LZImaging LADAR在2008被確定為解決直升機漂移控制與危險感知的核心技術［8］。該項目由美國空軍研究實驗室的相關部門負責研究，已經成功完成飛行實驗。由于系統工作距離短，直升機飛行姿態相對穩定，所以直升機3D-LZ Imaging LADAR沒有安裝在穩定平臺上，而是直接由INS慣導系統、激光3D成像系統、掃描系統組成，并結合了圖像處理系統與顯示系統。系統相應參數如表2所示，系統結構如圖3所示。

圖3 3D-LZ LADAR系統結構圖Fig.3 3D-LZ LADAR system block diagram

(4)機載對地高分辨詳查系統

記載高分辨詳查是高分辨縮放技術的基本前提。美國國家生態觀測網(National Ecological Observatory Network，NEON)為了掌握和預測大地氣候的變化、土地用途變化、植被生態情況等［9］，研制了適應于低空飛機的遙感詳查系統。系統載荷包括高分辨的機載成像光譜儀、機載激光成像系統和高分辨相機，輔助系統包括差分GPS和INS。其激光成像系統的主要指標如表3所示。

表3 系統參數Tab．3 System parameter

上面介紹的4個機載激光3D成像系統，由于工作任務的不同技術方案也不盡相同。JIGSAW項目的主攻目標是：(1)研制一種小型、集成、低功耗的3代無人機機載激光3D成像系統。(2)探討系統在遮蔽環境下進行目標的識別與探測的可行性。(3)推進數據的實時傳輸與處理技術在項目中的應用。

LADS系統采用特定波長的激光器實現對海洋深度的勘察，更具有大視場范圍，并且測量數據為后處理，不要求圖像信息的實時傳輸與處理。

直升機3D-LZ Imaging LADAR系統成像距離近，在系統穩定平臺結構選擇上大大簡化，直接由INS慣導系統、激光3D成像系統、掃描系統組成，并必須要求時時圖像處理與顯示系統。

機載對地高分辨詳查系統要求得到更加豐富的遙感信息，其激光成像系統要求更高的分辨率，載機工作高度較LADS高，對導航系統、掃描控制系統的要求也更高。

2.2 國內現狀

國內激光3D探測成像系統(LADAR)研究起步較晚，目前處于前期階段。主要有：電子部27所研制的直升機防撞激光3D成像系統，系統采用半導體泵浦的YAG激光器，利用兩個諧振鏡進行掃描;華中科技大學研制的海洋探測激光3D成像系統，系統采用YAG調Q倍頻激光器，利用卵形螺旋掃描方式;哈爾濱工業大學研制的障礙物回避用激光3D成像系統，已研制出實驗室樣機，采用1.06 μm半導體泵浦YAG激光器，利用兩個諧振鏡進行掃描成像，成像速率為7 frame/s，幀分辨率為32 pixel×32 pixel，作用距離為2 km，回波強度等級為16級。

整體而言，國內機載激光3D成像系統方面的研究已經取得一定的成果，但與國外相比存在較大差距，相關報道十分有限。

3 工作原理

以MIT林肯實驗室研制的基于蓋革模式雪崩光電二極管(APD)陣列的激光3D成像系統為例，原理如圖4所示［10，11］。高重頻脈沖激光發散光束，照明整個欲成像的場景;反射回的光照射到2D APD陣列上，APD陣列測得返回的光到達時間，而不是回波強度;每個像元給出距離值，從而對每一激光脈沖，激光3D成像系統獲得角-角-距3D像。

圖4 蓋革模式APD陣列的激光3D成像雷達原理圖Fig.4 Schematic diagram of Geiger-mode APD array laser 3D imaging radar

區別于傳統光學系統的探測光強，無論是基于焦平面成像(如APD陣列)還是單像元掃描，激光3D成像系統均是探測每個像元上激光發射到返回的傳播時間，結合激光光束的空間指向角信息，形成角-角-距測量數據，進而處理得到3D圖像。

機載激光3D成像系統在激光光軸的空間指向控制與獲取上與其它激光成像系統區別較大，也是機載成像系統的核心技術難點之一［12］。機載激光3D成像系統的工作原理如圖5所示。

圖5 機載激光3D成像系統原理圖Fig.5 Schematic diagram of airborne laser 3D imaging system

圖5中簡要介紹了兩個坐標系分別為地理坐標系(t系)、機體坐標系(b系)：

t系(O-xtytzt)：以地面某點 O 為原點，Oxt，Oyt，Ozt軸分別指向東、北、天空3個方向;

b 系(Ob-xb，yb，zb)，其中 Obyb為橫滾軸，Obxb為俯仰軸向，Obzb為航向軸。b系是個虛擬坐標系，為研究方便，可以認為捷聯慣導測量坐標系、穩定平臺測量坐標系與機體坐標系重合，統一為b 系 Ob-xbybzb。

掃描過程中，光軸的實時指向PGb在Ob-xbybzb中可以得到(由陀螺穩定平臺的角度測量信息、擺掃鏡的角度測量信息通過坐標變化得到)。P點在地理系中的位置［XP，YP，XP］可由全球GPS或差分GPS(DGPS)精確獲得。

機載INS可以測量得到飛機實時的姿態，即b系到 t系的轉移矩陣 Ctb，由此可以得到PGt=CtbPGb即地理系下光軸的指向信息。P點到G點的距離可以由激光測距得到，所以由此可以解算得到G點的地理系坐標。

機載激光3D成像系統在此基礎上通過機載方式對目標進行持續的采樣，其實現方式主要分為掃描式激光3D成像技術與推掃式激光成像技術兩種。

4 單元關鍵技術

機載激光3D成像系統由于測距載體載機自身的運動(質點運動、姿態運動)，激光測距方向的掃描控制與感知，載機位置姿態的感知(INS/GPS)都是機載激光3D成像系統的關鍵技術。全系統關鍵技術一般可以分為總體關鍵技術和單元關鍵技術。由于機載激光3D成像系統分為穩定平臺、激光發射/接收測距系統、INS/DGPS組合導航系統等，所以系統必然設計分系統間誤差分配、多系統安裝誤差標定補償等系統集成的總體關鍵技術。

此外，研制過程中系統的單元技術更為核心，從工作目標上分為激光測距系統、激光光軸控制與指向測量系統、數據圖像處理與顯示系統。

4.1 激光測距系統

激光測距系統是激光3D成像系統的核心。可以概括為激光光源、探測器(G-APD陣列探測器)、激光發射/接收光學系統。

4.1.1 激光光源

激光發射系統設計的主要指標有：發射激光的峰值功率、脈沖寬度、重復頻率和發散角等。由激光3D成像系統距離方程可知，APD探測器上接收到的激光回波的峰值功率和發射激光的峰值功率成正比，即發射激光的峰值功率越高，APD探測器上接收到的激光回波的峰值功率也越高，從而作用距離也越遠。所以為了提高作用距離，應盡可能采用峰值功率高的激光器。

激光脈沖寬度減小有利于提高時刻鑒別精度，有利于減少半導體激光器發熱從而能提高重復頻率。但激光脈沖寬度越小意味著相應的電信號帶寬越寬，從而要求APD探測器、跨阻放大電路和壓控放大電路的帶寬越寬，半導體激光器的驅動電源輸出電流脈沖越窄，所以選擇激光脈沖時要綜合考慮。作為激光測距系統的關鍵部件，激光器的選擇十分關鍵。

(1)半導體激光器

半導體激光器的應用已經覆蓋了整個光電子學領域。美國是半導體激光器泵浦技術開發應用較早、進展最快的國家，日本、英國和德國的西門子公司早在1988～1989年就提供了高功率半導體激光器系列產品用作泵浦光源，目前國外用半導體激光器疊層陣列泵浦Nd∶YAG激光器，輸出功率可達到1 000 W，未來工業和國防應用的半導體激光泵浦技術將產生巨大變革［13］。

利用半導體激光器泵浦 Nd∶YAG(0.946 μm)，Nd ∶YAG(1.06 μm)，Nd ∶YVO4(1.3 μm)，Nd∶YAG(1.3 μm)，Nd∶BEL(LaBe2O3)(1.07 μm)，Ho ∶YAG(2.1 μm)，Er ∶YLF(2.3 μm)，Nd∶玻璃 (1.06 μm) 和 Er∶玻 璃(1.3 μm，1.5 μm)，近年來都有許多報道，有些已經商品化［14］。

二極管激光測距儀開發較早，小于1 km測量距離的商用測距儀已達到實用水平，用于測距報警系統、航海浮標測距、集裝箱檢查等。GaAs激光測距儀以每秒數千次的脈沖重復頻率工作，在距離幾公里內精度可達幾厘米。上世紀80年代，美國科頓公司的M931型半導體激光測距夜視儀將單目夜視裝置和GaAs半導體激光器集為一體，測程超過1 km，重量為1.3 kg。美國國際激光系統公司的GR500型激光測距機，采用GaAs激光器，重復頻率為2 kHz，脈寬為40 ns，發散角為5～50 mrad，測程為3 230 m，重量為10 kg。美國輕型反坦克武器激光測距機采用GaAs激光器，輸出功率為40 W，脈寬為70 nm，發射角為10 mrad，重復頻率為5.7 kHz，測距 ＞500 m。隨著半導體激光泵浦固體激光器的發展，美國麥道公司已將它引入軍用市場，90年代初開始在F/A-18戰斗機上進行試驗，并投入生產。用氣冷的激光二極管泵浦Nd∶YAG，輸出波長為1.064 μm或532 nm，脈沖能量達 200 mJ，工作溫度為-35℃～ +60℃，該裝置的重量為4.5～5.7 kg［15，16］。

(2)微片激光器

美國的林肯實驗室已經開發出與3D雷達成像系統高度兼容的尖端激光器，其中用途最多的一種叫做微片激光器，如圖6所示。其主要采用ND-YAD泵浦和Cr4+泵浦，被動調Q方式，關鍵水平在于諧振腔的長度約為1 mm，僅支持小于增益帶寬的縱向模式激光器。小型激光器的另一個分支技術是它的輸出脈沖寬度＜1 ns，作為激光3D成像系統最佳的脈沖寬度［10］。

圖6 微片激光器Fig.6 Mini size laser

微型激光器發展的目標是增加單一脈沖的輸出能量。目前激光器可以達到脈沖重復頻率為1 kHz，輸出能量為250 J。而脈沖寬度低于ns量級(380 ps)的激光器正在研制中［10］。

如圖7所示，根據微片激光器重復頻率和單一脈沖輸出能量的經驗關系，最高的光子計數效率對應激光器有一個非常高的重復頻率。這樣一來，單一脈沖能量會非常小，并且主控振蕩器的放大率也應該適當對應。在這種情況下，整個系統的主要振蕩將來自于微片激光器［11］。

圖7 微片激光器的輸出性能Fig.7 Output performance of mini size laser

4.1.2 直接探測距離的陣列接收技術—G-APD陣列探測器

工作在蓋革模式下G-APD的反向偏置電壓高于雪崩電壓。如果接收到一個光子，會發生雪崩現象，使電流達到最大值，這個過程通常是瞬態的(不到1 ps)。G-APD陣列探測器兼具單光子探測靈敏度和皮秒級時間分辨率兩大特點，適用于對極微弱光目標的3D成像探測。同時，G-APD陣列探測器又是一種全固態的光電探測器件，不僅體積小、重量輕、可靠性高，而且還可利用現有的微光電子工藝實現規模化生產。

近來，林肯實驗室開發了一種32×32陣列的Si基G-APD探測陣列(如圖8所示)。鑒于此陣列探測器設計為用于微弱光目標角-角距3D成像，因此各像元電路的功能只為測量光子的到達時間［17］。G-APD的每個像元均采用平面結構。像元光敏面尺寸為 30～50 μm，像元間距為100 μm。每個像元均有相應的電路單元，電路單元是一種數字計時器，具有類似于秒表的功能;每個像元的計數器均按規定的程序對公共的時鐘脈沖計數，只要發現有單光子入射某像元，該像元探測器輸出的雪崩電流脈沖就將終止相應計數器單元的計數并保存結果，用于讀出。

圖8 32×32 APD陣列探測器Fig.8 32×32 APD array detector

4.1.3 激光發射/接收光學系統

光學系統由兩部分組成：發射光學單元和接收光學單元。其中，發射單元主要實現激光擴束、分光和光線的準直輸出;接收單元主要實現光線接收、準直和聚焦成像。激光發射系統對激光器的輸出光束進行準直，準直方法主要有：單透鏡法、組合透鏡法、漸變折射率透鏡法等。采用組合透鏡法可以保證一定精度的準直效果，同時避免過于復雜的結構。圖9給出了MIT林肯實驗室采用的光學系統原理圖［10］。

圖9 組合透鏡法的激光發射/接受系統光路圖Fig.9 Optical pathway of laser emission/receiving system

圖9中激光器波長為532 nm，32×32 APD探測器，發射系統焦距f=300 mm，口徑D=11 mm。該系統基于APD陣列的閃光雷達是一種非掃描激光3D成像系統，并設計了一套正交光柵衍射分光光學系統。利用伽利略望遠鏡對532 nm脈沖激光進行準直擴束;針對APD陣列的特點，使用2D正交光柵和聚焦透鏡對擴束光進行衍射分光，光斑經發射-接收分光鏡和發射透鏡后照明APD像元的瞬時視場;利用發射-接收分光鏡使發射光路與接收光路分開。

4.2 激光光軸控制與指向測量系統

激光光軸控制與指向測量系統主要由包括INS/GPS組合導航系統，陀螺穩定平臺和激光掃描控制系統等組成［18］。

4.2.1 INS/GPS組合導航系統

載機的精確位置和姿態是機載激光3D成像系統測量數據的重要組成部分，載機位置姿態的精確測量設備即為INS/GPS組合導航系統。

因INS/GPS組合導航系統的INS和GPS均為全球方位、全天候段的導航設備，能提供十分完整的導航數據，可優勢互補。INS主要具備實時性好、分辨率高、自主性強的優點，但也同時具有誤差隨時間發散、技術復雜、價格昂貴的缺點。GPS(或者格洛納斯、北斗)與INS相比具有價格低廉、誤差不發散、位置精度高的特點，但同時GPS的實時性、分辨率及自主性較差［19］。

組合導航系統可以綜合兩種導航系統的優點，結合先進的數據融合濾波技術可以大大提升單系統的能力。區別于傳統組合導航系統，成像補償用組合導航系統精度要求更高，我國在這方面技術差距很大。

表4給出了國內外應用于合成孔徑雷達(SAR)成像運動補償的組合導航系統參數。如加拿大APPLANIX公司的SAR運用撓性陀螺慣性測量單元(IMU)/GPS組合導航系統POS/AV510，POS/AV510應用于美國Sandia國家實驗室研制的Lynx SAR，使Lynx SAR實現了0.1 m分辨率［20］，具體參數見表4，其中 C/A與 DGPS是衛星導航中不同的定位解算方式。

表4 POS/AV510系統參數Tab．4 POS/AV510 system parameters

我國目前唯一用于成像運動補償類的組合導航系統是北航與航空618所聯合研制的基于SINS/DGPS的機載SAR運動補償系統SARMC，系統精度較國外差距較大，具體參數如表5所示。

表5SARMC系統參數Tab．5 SARMC system parameters

可見，INS/GPS組合導航系統是實現高分辨率成像系統的關鍵技術。

4.2.2 陀螺穩定平臺

陀螺穩定平臺用以隔離載機的俯仰、橫滾、航向姿態運動，為有效載荷(激光3D成像系統的測距系統)提供良好的工作環境，針對不同任務需求的激光3D成像系統，對陀螺穩定平臺的要求也不同。

2.1節提到的直升機激光3D著陸成像系統采用了單軸(俯仰軸)穩定，并且該俯仰軸集合了俯仰掃描功能，另有擺掃鏡進行水平擺掃。采用該結構可以有效降低系統的成本，光軸控制精度滿足直升機的工作要求。

澳大利亞政府支持的LADS系統，激光3D成像系統的測距系統安裝在3軸穩定平臺中(要求隔離俯仰、橫滾、航向運動)，擺鏡為橫向擺掃，隨著載機的前進推掃成像［21］。

美國的JIGSAW系統則采用高精度的4軸穩定平臺，可在GPS/INS組合導航系統的輔助下提供高精度的視軸穩定控制和高精度的激光指向信息測量［22，23］。

而用于航天器的激光3D成像系統，由于航天器的運動姿態穩定，一般激光3D成像系統的測距系統直接捷聯安裝在載體上，需要根據要求設計擺鏡的光軸擺掃方式，可以不用穩定平臺。

機載穩定平臺的結構形式有別于傳統航空吊艙穩定平臺的要求，相近于航空相機穩定平臺結構，旨在隔離載機的俯仰、橫滾運動。航空吊艙中常見的光電陀螺穩定平臺形式有：2軸2框架穩定平臺，2軸4框架穩定平臺，3軸3框架穩定平臺。文獻［24］給出了國內外機載光電平臺的主要性能指標對比列表。

比較適合機載激光3D成像系統的機構形式為3框架3軸穩定平臺系統(也可去掉方位軸)。該平臺是一個有3個自由度的系統：方位、俯仰、橫滾，可以隔離飛機方位、俯仰、橫滾方向的運動，即消除了機載環境對成像的干擾，提供了較好的成像環境。

為了進一步消除各種大氣現象和空氣動力現象顫振、擾流、抖振、結冰等引起的振動所產生的視軸瞄準線的晃動，提高圖像質量和指向精度，可采用一種內框架減振系統。該系統采用模塊式組合方式，具有靈活性和通用性，方便安裝調試［25］。另外為了解決2軸、3軸光電平臺中陀螺儀的解耦問題，提高光電平臺的穩定性，文獻［26］給出了2軸平臺正割解耦的算法，并給出了3軸平臺正弦解耦的算法以及基于ARM單片機的解耦軟硬件處理技術。

以文獻［27］給出的光電偵察吊艙的平臺結構為例(圖10)，光電偵察設備固定在方位框內，方位軸、俯仰軸與橫滾軸既起方位框、俯仰框與橫滾框及基座之間的連接作用，又同時實現了3個方向的運動，其中橫滾軸與俯仰軸消除機載環境的運動干擾，目標的穩定跟蹤由方位軸與俯仰軸來完成。光電探測器的角位置跟蹤誤差由位置外環來補償，速度內環用來消除環境干擾和角速度指令跟蹤。

圖10 3軸3框架示意圖Fig.10 Schematic diagram of three framework three axes

4.3 數據圖像處理與顯示系統

4.3.1 數據處理技術

機載激光3D成像系統測量的數據處理技術涵蓋多個方面，包括動態GPS數據后處理、INS和GPS組合姿態確定、不同傳感器觀測值的時間系統同步處理、激光腳點3D坐標計算、坐標系統的變換、系統誤差的改正、粗差的剔除、數據的濾波分類，還包括DEM/DTM的生成以及后續的地物提取、建筑物3D重建、3D模型建立等高級處理技術。其中激光腳點坐標的計算需要3組數據：系統檢校和安置參數、激光測距值及位置和姿態數據。為計算出數字高程模型及地物提取，先將激光腳點數據進行坐標轉換，然后進行濾波分類，即地面點與非地面點分開。通過數據濾波后，能內插生成數字高程模型(DEM)，以及后續的地物提取、建筑物3D重建、3D模型等高層次處理。

在國外，機載激光3D成像系統測量技術已經相當成熟，硬件技術與系統問題研究基本解決，但相應的數據處理的算法還處于前期研究階段，其中如何消除系統誤差對3D激光腳點坐標的影響、激光3D成像系統數據的濾波分類以及地物識別和提取是關鍵研究內容［28］。

目前中國科學院遙感應用研究所采用一種組合CCD圖像和稀疏激光測距數據的建筑物3D信息提取的方法。該方法采用高分辨率的航空CCD圖像和機載激光掃描測距數據準確而自動地提取了大部分建筑物的3D信息，并充分考慮了CCD圖像的高分辨率特性和激光掃描測距數據的3D特性，做到互相補充。通過CCD圖像提取建筑物的輪廓信息，并根據激光測距數據提取建筑物的高度信息，二者結合起來重建了建筑物的 3D 信息［29］。

4.3.2 測量誤差處理

機載激光3D成像系統的測量誤差主要包括：激光掃描測距誤差、GPS動態定位誤差、INS姿態測量誤差、動態時延誤差、二類高程誤差、掃描角誤差、激光束發散角所產生的系統誤差和系統集成綜合誤差。為了提高機載激光3D成像系統的測量精度，最大可能地降低各種系統誤差的影響，通常采用建立誤差改正模型、儀器檢校和條帶平差改正的方法［30，31］。

4.3.3 數據濾波和分類

從激光腳點數據點云中提取數字地面高程模型(DTM/DEM)，需要將其中的地物數據腳點去掉，即進行激光3D成像系統的濾波。在進行地物提取和物體的3D重建過程中，植被數據點和人工地物點的區分，以及城市建筑物數據點云的提取，即為激光3D成像系統的數據分類。研究如何從激光數據點云中分離出地形表面激光腳點數據子集，以及區分不同地物激光腳點數據子集，就是數據的濾波和分類，有時濾波和分類是同時進行的。

目前，用于機載激光3D成像系統測量數據濾波的方法多是基于3D激光腳點數據的高程突變信息進行的，大致可以分為形態學濾波法、移動窗口法、高程紋理分析法、迭代線性最小二乘內插法、基于地形坡度濾波等幾種方法。但是幾乎每種方法都有自身的缺陷，有待改進，有些方法需要分步反復迭代進行。通常所有的方法都強調高程低的腳點屬于地面腳點，高程較高的腳點不在地面上，這樣會帶來一定的系統誤差。所以濾波算法一般采用半自動或人工操作方式，能否設計出一種簡單、實用、穩健的濾波算法，是機載激光3D成像系統數據處理中必須解決的關鍵問題之一。圖11給出的是城區內機載激光3D成像系統數據濾波前后的對比圖，結果表明，濾波算法能夠優先過濾出大型建筑物［32］。

圖11 數據濾波前后對比圖Fig.11 Data comparison before and after filtering

如果要進行地物提取，就必須在濾波的基礎上對地物腳點進行分類，以實現機載激光3D成像系統自動提取地面目標。目前絕大多數算法都是先將原始數據直接內插成規則格網的距離圖像，在此基礎上進行基于高程紋理特征的分類。即利用2D GIS數據進行分類;基于小波變換和尺度空間理論進行分類以及利用局部直方圖分析技術并融合激光強度信息進行分類。圖12是機載激光3D成像系統數據經過濾波分類后的效果圖，圖中所示的車輛和植被被過濾了，保留了地形的特征線［32］。

圖12 數據濾波分類效果圖Fig.12 Classification effect of data filtering

4.3.4 激光回波信號強度與激光腳點高程的融合

機載激光3D成像系統的測量系統不僅能提供數據點的高程信息，且越來越多的系統能夠提供激光回波信號的強度信息。激光脈沖打到相同的物質表面時，其回波信號的強度較為接近，而每種物質對激光信號的反射特性是不一樣的。利用介質材料的反射率能夠確定激光回波的能量值，從而對不同介質材料進行標定，根據這一特點，能夠方便地區分不同的地物，如道路、沙石、草地和灌木叢等。

圖13右側為掃描區域真實的相片，左側為激光3D成像系統數據融合激光強度信息和高程信息進行分類后的效果圖。將系統的激光回波信號的強度信息與激光腳點高程信息融合，進而實現數據的分類，能夠提高系統對地形的測量精度，目前此問題在國際上還處于研究階段［33］。

圖13 融合強度信息濾波分類結果Fig.13 Classification results of aggregation strength information

4.3.5 建立3D模型

機載激光3D成像系統技術的發展已有20多年的歷史，但利用機載激光3D成像系統數據進行3D重建的研究并不多。目前主要有以下幾種方法：1995年Weidner等人從高空分辨率的數字表面模型中自動提取規則棱柱形房屋模型;1999年Haala等人聯合多光譜影像和激光3D成像系統數據進行建筑物和樹的模型提取;1999年Mass等人基于分析激光腳點點云的矩不變提出了從機載激光3D成像系統數據中自動提取建筑物模型的方法。

激光腳點數據經過濾波和分類后，形成一系列彼此分開的點群信息標定其地物屬性，這樣就確定了屋頂激光腳點點集，并采用掃描帶法分離出每個建筑物的激光腳點集。通過分析這些激光腳點間的幾何拓撲關系，利用3D Hough變換就能重建出屋頂的表面形態。將分類得到的房屋激光腳點垂直投影到xy平面，這些投影點的幾何構形就反映了建筑的邊界、外形、大小、方位等幾何參數［34］。

5 結束語

本文對機載激光3D探測成像系統的國內外發展現狀進行了總結，并介紹了激光3D成像系統的工作原理，重點對系統的單元關鍵技術進行了詳細的剖析。機載激光3D成像技術能夠提供極高的角分辨率和距離分辨率，還可以合成目標的高分辨率3D圖像，而且3D圖像的生產周期短，因而在軍事(例如導彈的精確制導等)和民用領域(例如公路鐵路勘探等)具有極高的應用價值。除了文中探討的單元關鍵技術，還需考慮系統的重量、體積以及成本等多方面因素。機載激光3D成像作為一種新興技術，具有傳統技術無法比擬的優點，必將在未來的國防和國民經濟中占有越來越重要的地位。

［1］ JUBERTS M，BARBERA A.Status report on next generation LADAR for driving unmanned ground vehicles［J］.SPIE，2004，5609：1-12.

［2］ MAARTEN U H，VENABLE D，SMEARCHECK M.Use of 3D laser radar for navigation of unmanned aerial and ground vehicles in urban and indoor environments［J］.SPIE，2007，6550：65500C.

［3］ STONE W C，JUBERTS M，DAGALAKIS N，et al..Performance analysis of next-generation LADAR for manufacturing，construction，and mobility［G］.Gaithersburg：National Institute of Standards and Technology(NIST)，2004.

［4］ JUBERTS M，BARBERA A.Status report on next generation LADAR for driving unmanned ground vehicles［J］.SPIE，2004，5609：1-12.

［5］ PEREIRA J，MOEBIUS B，PFENNIGBAUER M，et al..Imaging LIDARs for space applications［J］.SPIE，2008，7061：706100J.

［6］ VAIDYANATHAN M，BLASK S，HIGGINS T，et al..Jigsaw phase Ⅲ：a miniaturized airborne 3-D imaging laser radar with photon-counting sensitivity for foliage penetration1［J］.SPIE，2007，6550：65500N.

［7］ ANTHES J P.Non-scanned LADAR imaging and applications［J］.SPIE，1993，1936：11-22.

［8］ SAVAGE J，HARRINGTON W，MCKINLEY A R，et al..3D-lz helicopter ladar imaging system［J］.SPIE，2010，7684：768407.

［9］ BRIAN R J，THOMAS U K，MICHELE K.Development of airborne remote sensing instrumengtations for NEON［J］.SPIE，2010，7809：78090I.

［10］ RICHARD M M，TIMOTHY S，ROBERT E H，et al..A compact 3D imaging laser radar system using Geiger-mode APD arrays：system and measurements［J］.SPIE，2003，5086：1-12.

［11］ RICHARD M M，BRIAN F A，RICHARD M M，et al..Three-dimensional laser radar with APD Array［J］.SPIE，2001，4377：217-228.

［12］ LI S K，XUE Y Q，LIU ZH.A new airborne remote sensing platform for generating geocoding image without ground control point［J］.IEEE，1997：1987-1989.

［13］ 李學千，張興德.高功率半導體激光器的進展及其在軍事上的應用［J］.光學精密機械，1993(2)：1-5.LI X Q，ZHANG X D.Development and application in military of high power semiconductor laser［J］.Opt．Precision Eng.，1993(2)：1-5.(in Chinese)

［14］ 王德，李學千.半導體激光器的最新進展及應用現狀［J］.光學 精密工程，2001，9(3)：279-283.WANG D，LI X Q.New progress in semiconductor lasers and their applications［J］.Opt．Precision Eng.，2001，9(3)：279-283.(in Chinese)

［15］ 付有余.激光測距與跟蹤系統低噪聲電子設計方法［J］.光學 精密工程，2000，8(3)：250-253.FU Y Y.Low noise electronic design method for laser range finding and tracking system［J］.Opt．Precision Eng.，2000，8(3)：250-253(in Chinese).

［16］ 何興仁.大功率半導體激光器及其在軍事上的應用［J］.半導體情報，1995，32(3)：26-33.HE X R.High power semiconductor laser and its application in military［J］.Semiconductor Information，1995，32(3)：26-33.(in Chinese)

［17］ MUGUIRA M R.Range determination for scannerless uImaging：US，6088086［P］.2000-11-11.

［18］ SLOBODAN N V.Suppression of torsional oscillations in a high-performance speed servo drive［J］.IEEE，1998，45(1)：1-8.

［19］ 房建成.先進導航技術講座［Z］.北京航空航天大學，2006：23-32.FANG J CH.Advanced navigation technology lecture［Z］.Beijing：Beijing Aeronautics and Astronautics University，2006：23-32.(in Chinese)

［20］ MARCELO C，ALGRAIN D，EHLERS E.Suppression of gyroscope noise effect in pointing and tracking system［J］.SPIE，1994，2221：1-8.

［21］ BO L.Nonlinear induced disturbance rejection in inertial stabilization systems［J］.SPIE，1996，2739：1-8.

［22］ 石紅生，盧廣山.一種新型狀態觀測器在陀螺穩定平臺中的應用［J］.電光與控制，1999(1)：24-28.SHI H SH，LU G SH.Application of a new type state observer in gyro stabilized platform［J］.Electronics Opt．Control，1999(1)：24-28.(in Chinese))

［23］ 王合龍，朱培申，姜世發.陀螺穩定平臺框架伺服系統變結構控制器的設計和仿真［J］.電光與控制，1998(2)：24-29.WANG H L，ZHU P SH，JIANG SH F.Design and simulation of framework servo system variable structure controller of gyro stabilized platform［J］.Electronics Opt．Control，1998(2)：24-29.

［24］ 賈平，張葆.航空光電偵察平臺關鍵技術及其發展［J］.光學 精密工程，2003，11(1)：82-88.JIA P，ZHANG B.Critical technologies and their development for airborne opto-electronic reconnaissance platforms［J］.Opt．Precision Eng.，2003，11(1)：82-88.(in Chinese)

［25］ 甘至宏.光電吊艙內框架減振系統設計［J］.光學 精密工程，2010，18(9)：2036-2043.GAN ZH H.Design of inner frame vibration absorbing system for optoelectronic pod［J］.Opt．Precision Eng.，2010，18(9)：2036-2043.(in Chinese)

［26］ 陳占軍，葛文奇，李英，等.光電平臺中陀螺解耦問題的研究［J］.光學 精密工程，2008，16(9)：1712-1715.CHEN ZH J，GE W Q，LI Y，et al..Research of gyro decoupling problem on photo-electric platform［J］.Opt．Precision Eng.，2008，16(9)：1712-1715.(in Chinese)

［27］ 鄒志峰，閆峰.機載光電偵察穩定平臺系統分析［J］.機械工程師，2007，16(4)：39-40.ZOU ZH F，YAN F.Analysis for the system of air borne optoelectronic reconnaissance platforms［J］.Mechanical Eng．，2007，16(4)：39-40.(in Chinese)

［28］ 李清泉，李必軍，陳靜.激光雷達測量技術及其應用研究［J］.武漢測繪科技大學學報，2000，25(5)：387-392.LI Q Q，LI B J，CHEN J.Research on laser range scanning and its application［J］.J．Wuhan Technical University Surveying and Mapping，2000，25(5)：387-392.(in Chinese)

［29］ 尤紅建，張世強.組合CCD圖像和稀疏激光測距數據的建筑物三維信息提取［J］.光學 精密工程，2006，14(2)：297-302.YOU H J，ZHANG SH Q.Reconstructing 3D buildings based on airborne CCD image and laser scanning rangef inder data［J］.Opt．Precision Eng．，2006，14(2)：297-302.(in Chinese)

［30］ 章大勇，吳文起，吳美平.機載激光雷達系統標定方法［J］.光學 精密工程，2009，17(11)：2806-2813.ZHANG D Y，WU W Q，WU M P.Calibration technology of airborne lidar［J］.Opt．Precision Eng.，2009，17(11)：2806-2813.

［31］ 劉經南，張小紅，李征航.影響機載激光掃描測高精度的系統誤差分析［J］.武漢大學學報·信息科學版，2002，27(2)：111-117.LIU J N，ZHANG X H，LI ZH H.Analysis of systematic error influencing the accuracy of airborne laser scanning altimetry［J］.Geomatics and Information Sci．Wuhan University，2002，27(2)：111-117.

［32］ 劉經南，張小紅.利用激光強度信息分類激光掃描測高數據［J］.武漢大學學報·信息科學版，2005，30(3)：189-192.LIU J N，ZHANG X H.Classification of laser scanning altimetry data using laser intensity［J］.Geomatics and Information Sci．Wuhan University，2005，30(3)：189-192.

［33］ 張小紅，劉經南，李征航.3-D機載激光掃描測高數據濾波和分類研究［J］.武漢大學學報·信息科學版，2002，27(增刊)：181-183.ZHANG X H，LIU J N，LI ZH H.Research on the filter and classification of 3-D airborne laser scanning altimetry Data［J］.Geomatics and Information Sci．Wuhan University，2002，27(s)：181-183.

［34］ 張玉方，程新文，歐陽平，等.機載Lidar數據處理及其應用綜述［J］.工程地球物理學報，2008(1)：119-124.ZHANG Y F，CHENG X W，OUYANG P，et al..The data processing technology and application of airborne LIDAR［J］.zChinese J．Eng．Geophysics，2008(1)：119-124.