基于平面約束的機載LiDAR航帶平差方法研究

基于平面約束的機載LiDAR航帶平差方法研究

袁豹

(蘇州市測繪院有限責任公司,江蘇 蘇州 215006)

摘要：機載LiDAR技術作為一種新型的主動式空間觀測技術,在快速獲取空間三維地理信息方面具有不可估量的前景。然而，機載LiDAR系統獲得的相鄰航帶點云數據在重疊區存在“漂移”問題，需要采用航帶平差的方法實現不同航帶點云數據之間的“無縫”拼接。針對LiDAR點云數據難以獲取同名特征點來進行航帶平差迭代初始值計算的問題，文中在總結航帶平差迭代初始值求解方法的基礎上，提出一種基于點與平面約束關系的最小二乘求解方法，提高初始值求解的準確性，并通過具體實例數據進行驗證分析。

關鍵詞：機載LiDAR；航帶平差；點云數據；平面特征；迭代初始值

中圖分類號：P237

收稿日期：2014-07-27;修回日期：2015-01-10

作者簡介：袁豹(1988-)，男，碩士.

Research on the method of airborne LiDAR strip adjustment based on plane constraint

YUAN Bao

(Suzhou Surveying Institute Co.,Ltd.，Suzhou 215006, China)

Abstract：As a new type of active earth observation technology from space,the airborne LiDAR technology serves as an incalculable prospect in fast access to spatial three-dimensional geographic information.However,there has been a problem of data drifting in overlap region from two adjacent air-strips,and different strip adjustment methods are required to splice the data seamlessly.Aiming to obtain the initial value in strip adjustment iteration,on the basis of summarizing the related computing methods, it proposes a new algorithm using the least squares model to get a more accurate initial iteration value,which makes use of the constrains between the point and plane.Specific experiments are implemented to analyze and validate the effect of method presented.

Key words：airborne LiDAR；trip adjustment；point cloud data；plane feature；iterative initial value

機載LiDAR系統獲得的地面點云數據具有海量、散亂等特性，難以判別或提取重疊航帶點云之間的同名點，也很難在實際工程應用中很好處理航帶平差連接點。所以，在航帶平差過程中，迭代初始值的計算問題一直是學術界研究的熱點。在地面三維激光掃描點云粗配準中，有學者將圖像處理領域用于檢測圖像特征和邊緣、突變等性質的各種算子引入到粗配準的應用中，近年來，很多學者開始致力于將此方法應用于LiDAR數據處理中，該方法首先將點云轉換為深度圖像，利用各種作用不同的圖像檢測算子實現圖像特征點的提取，計算迭代初始值參數。然而，機載LiDAR點云數據平面精度較低，在轉為深度圖像時會由于插值處理導致誤差的累積和精度的損失。還有學者針對LiDAR點云數據本身進行處理以達到點云配準的目的，但目前不同方法的難點都集中在同名地物特征點的選取及配準方面，好的同名特征匹配方法是準確獲取航帶平差模型迭代初始值的關鍵。因此，針對迭代初始值的求解問題，本文提出一種基于點與平面約束關系的最小二乘模型來求解初始參數，并通過實例數據進行驗證分析。

1迭代初始值求解

1.1　算法流程設計

本文通過提取兩塊重疊點云同名幾何特征地物進行點云配準的方法，提出基于點與平面約束關系建立最小二乘模型的機載LiDAR點云初步配準方法。在機載LiDAR點云數據中，由于點云數據區域范圍大、點云數據散亂等特征的存在，線性特征地物的點云表示形式在機載LiDAR整個點云數據中并不明顯，也很難準確提取同名線性特征地物點云。而面狀地物點云在機載LiDAR整個點云數據中表現較為明顯，并且建筑平面等面狀地物點云高程屬性具有良好的判別性，易于從點云中提取，所以尋找兩重疊航帶點云的同名地物點云數據并進行匹配處理是一種較為可靠的理論方法。

基于面狀地物點云共面條件約束的配準方法的主要思想和操作流程如圖1所示。

1.2　基于點與平面約束關系的最小二乘模型

點云濾波可以方便快捷地獲得面狀地物的點云數據，假設兩面狀地物點云數據分別為S1和S2，根據S1和S2中的點云數據可以分別進行平面擬合獲得平面方程，定義平面方程為

Ax+By+Cz+D=0.

(1)

定義兩點云之間的轉換關系式為

(2)

(3)

其中，A=

(4)

將式(4)改寫為未知參數的誤差方程，則有：

(5)

在實際計算過程中，點云數據進行平面方程擬合后，并不完全符合該平面方程的函數關系，所以最小二乘平差模型在建立時以點與平面差距最小為準則進行約束。

2實例數據計算分析

2.1　數據來源

本文實驗數據來源于某庫區機載激光雷達航測及1∶2 000 比例尺DOM、DEM、DSM成果制作項目。測區位于黃河主干道，平均海拔350m，最高海拔800m，屬于山地地形。

數據采集使用的機載LiDAR系統為加拿大進口的國際上性能最為優越的機載LiDAR航測設備ALTMGemini。項目外業數據采集時間為2012年7月，飛行高度為1 000m(比例尺1∶2 000)，得到LAS1.0標準格式數據。鑒于數據保密性的要求，本文截取部分點云數據進行實驗分析，該部分點云包含有較多居民地、道路等地物信息，同時植被等地貌信息也較為豐富。相鄰航帶重疊率約為20%，平均掃描點密度為2.20點/m2，第一條航帶點云個數為1 830 091，數據量為48.8 MB，第二條航帶點云個數為1 598 130，數據量為42.6 MB。實驗區域點云高程圖如圖2所示。

2.2　實驗方案與結果

為了驗證本文提出的迭代初始值計算方法的有效性，比較采用點與平面約束關系計算得到的迭代初始值和根據人工選取的特征點計算迭代初始值兩種方法對航帶平差結果的影響，從航帶平差運行時間和模型精度等方面對結果進行對比分析。其中人工選取特征點的個數為四個，利用目視判斷選取的特征點作為控制，建立兩條航帶點云坐標轉換方程，求解出轉換參數作為迭代初始值。

根據濾波后獲得的非地面點云，在建筑物類別點云中選取同名平面屋頂表面的兩個航帶的點云，按照高程值不同對點云進行著色，顯示如圖3、圖4所示。

圖3　S 1″面狀地物點云

圖4　S 2″面狀地物點云

從上述兩幅平面屋頂點云高程屬性分布情況可以看出，點云高程具有一定的分散性。結合兩幅點云高程值屬性的具體分布表，可以發現：第一個屋頂點云數據的高程變化范圍為+124.57m~+126.21m，高程集中分布在+124.76m~+124.96m，最大最小高程差值為0.18m；而第二個屋頂點云數據的高程變化范圍在+123.31m~+126.14m之間，高程主要分布在+124.76m~+124.97m，最大高差為0.21m。所以獲取屋頂點云數據仍然存在粗差，這些粗差的存在對點云平面的擬合及初始參數的轉換具有很大的影響。因此，需要對點云數據進行粗差剔除處理。

(6)

首先根據面狀點云S1″進行平面擬合，以S2″中點云作為觀測值建立點與擬合平面方程之間的約束關系，建立最小二乘平差模型解算得到三個平移參數和三個旋轉角參數，比例因子初值設定為1.000，迭代初始參數值計算結果見表1，其中平移參數依次為Δx,Δy,Δz，旋轉角度參數表示為?x,?y,?z。

表1　初始迭代參數結果

根據相鄰兩航帶點云高程圖，根據點云顯示屬性及位置選取四對“概略同名控制點”，計算得出航帶平差迭代初始值如表2所示。

表2　人工選取控制點所得初始迭代參數

依據表1和表2中的迭代初始值，在相同的迭代限差條件下分別對相鄰航帶重疊區點云數據進行航帶平差處理，得到的平差結果如表3所示，其中方法 Ⅰ 表示人工選取特征點計算迭代初始值，而方法 Ⅱ 為本文提出的基于點與平面約束關系的計算方法。

表3　航帶平差結果對比

從表3的結果可以看出，兩種方法求得的參數結果具有一定的差異，采用基于點與平面約束關系建立模型求解得到的初始值與航帶平差后計算所得的轉換參數較為接近；而采用人工選取同名控制點計算迭代初始值進行航帶平差所得的轉換參數與給定的初值之間相差較大，說明僅憑人工選取同名控制點進行迭代初始值的計算具有不可靠性。兩種方法獲得航帶平差效果對比看出，本文采用的方法獲得的航帶平差模型中誤差較小，且迭代計算的次數和運算時間都遠遠小于人工選取同名控制點的平差方法，并且很好地保證航帶平差結果的可靠性和算法執行的效率。

3結論

針對航帶平差過程中的難點問題，提出一種依據點云平面特征進行航帶匹配，建立點與平面約束關系的最小二乘模型求解航帶平差迭代初始值。實例數據分析表明，本文中提出的方法可以有效地提高航帶平差效率和可靠性。如何從LiDAR點云數據中快速、有效地提取同名平面特征點云數據是需要進一步研究的方向。

另外本文提出的基于點與平面約束關系建立最小二乘模型進行航帶平差迭代初始值的求解方法，對航帶重疊部分點云數據的地物屬性有一定的要求，需要存在一些較為平整的地物，能夠進行平面擬合處理，從而建立良好的數學模型。對于平面特征缺失的航帶點云平差問題，研究其初始值求解方法還需進一步深入。

參考文獻：

[1]張小紅.機載激光雷達測量技術理論與方法[M].武漢：武漢大學出版社,2007.

[2]孫步陽.機載激光雷達航帶拼接技術研究[D].武漢:中國地質大學,2009.

[3]WANGLY,LIUZhJ,SONGWD,etal.AirborneLiDARstripadjustmentbasedonLSM[C].2010 3rdInternationalCongressonImageandSignalProcessing.China：ShangHai,October16-18,2010.

[4]王麗英.面向航帶平差的機載LiDAR系統誤差處理方法研究[D].遼寧:遼寧工程技術大學,2011.

[5]劉春,陳華云,吳杭彬，等.激光三維遙感的數據處理與特征提取[M].北京：科學出版社,2010.

[6]梁靜,范洪冬,鄧喀中.基于多分辨率分析的LiDAR數據處理與可視化[J].測繪工程,2011,20(4)：50-54.

[7]尚大帥,馬東洋,高振峰等.機載LiDAR點云數據與影像數據融合[J].測繪工程,2012,21(1)：18-20.

[8]CROMBAGHSM,MINED,BRUEGELMANNR.Ontheadjustmentofoverlappingstripsoflaseraltimeterheightdata[J].InternationalArchivesofPhotogrammetryandRemoteSensing,2000(33)：230-237.

[9]尹丹.古建筑三維重建中的深度圖像配準技術研究與應用[D].西安：西北大學,2008.

[10]田慧,周紹光,李浩.深度圖像自動配準點云的方法研究[J].測繪科學,2012,37(3)：41-43.

[11]馮聰慧.機載激光雷達系統數據處理方法的研究[D].鄭州:解放軍信息工程大學,2007.

[12]張熠斌.機載LiDAR點云數據處理理論及技術研究[D].西安：長安大學,2010.

[13]張靖,江萬壽.基于虛擬連接點模型的機載LiDAR系統安置角誤差自檢校[J].測繪學報,2011,40(6)：762-769.小紅.機載激光雷達測量技術理論與方法[M].武漢：武漢大學出版社,2007.

[2]孫步陽.機載激光雷達航帶拼接技術研究[D].武漢:中國地質大學,2009.

[3]WANGLY,LIUZhJ,SONGWD,etal.AirborneLiDARstripadjustmentbasedonLSM[C].2010 3rdInternationalCongressonImageandSignalProcessing.China：ShangHai,October16-18,2010.

[4]王麗英.面向航帶平差的機載LiDAR系統誤差處理方法研究[D].遼寧:遼寧工程技術大學,2011.

[5]劉春,陳華云,吳杭彬，等.激光三維遙感的數據處理與特征提取[M].北京：科學出版社,2010.

[6]梁靜,范洪冬,鄧喀中.基于多分辨率分析的LiDAR數據處理與可視化[J].測繪工程,2011,20(4)：50-54.

[7]尚大帥,馬東洋,高振峰等.機載LiDAR點云數據與影像數據融合[J].測繪工程,2012,21(1)：18-20.

[8]CROMBAGHSM,MINED,BRUEGELMANNR.Ontheadjustmentofoverlappingstripsoflaseraltimeterheightdata[J].InternationalArchivesofPhotogrammetryandRemoteSensing,2000(33)：230-237.

[9]尹丹.古建筑三維重建中的深度圖像配準技術研究與應用[D].西安：西北大學,2008.

[10]田慧,周紹光,李浩.深度圖像自動配準點云的方法研究[J].測繪科學,2012,37(3)：41-43.

[11]馮聰慧.機載激光雷達系統數據處理方法的研究[D].鄭州:解放軍信息工程大學,2007.

[12]張熠斌.機載LiDAR點云數據處理理論及技術研究[D].西安：長安大學,2010.

[13]張靖,江萬壽.基于虛擬連接點模型的機載LiDAR系統安置角誤差自檢校[J].測繪學報,2011,40(6)：762-769.

[責任編輯：李銘娜]