基于Terrasolid軟件的激光雷達點云數(shù)據(jù)處理

李蓮芳，邵秋銘

(重慶市勘測院，重慶 400020)

基于Terrasolid軟件的激光雷達點云數(shù)據(jù)處理

李蓮芳?，邵秋銘

(重慶市勘測院，重慶 400020)

隨著空間數(shù)據(jù)應用領域的不斷擴展，空間數(shù)據(jù)的獲取和處理成為研究學者的熱點問題。20世紀70年代機載雷達測量技術以一種主動式、快速、精確的新型遙感測量方法應運而生，本文主要研究激光掃描數(shù)據(jù)后處理技術中的濾波分類方法，并通過Terrasolid軟件進行點云數(shù)據(jù)處理，結(jié)合具體實例進行粗差剔除，濾波分類等操作。在此過程中，通過設定不同的濾波分類參數(shù)進行精度比較，以實現(xiàn)最好的濾波分類效果。

機載LiDAR；點云；濾波算法；粗差剔除；地表面提取

1 引 言

機載LiDAR是遙感測量領域的一門新興技術，經(jīng)過幾十年的發(fā)展國內(nèi)外先后研制出了多種機載激光掃描測高系統(tǒng)，同時出現(xiàn)了一些成熟的機載LiDAR數(shù)據(jù)處理軟件，但其功能算法均由生產(chǎn)廠家單獨封裝保密提供，用戶看不到其具體的算法步驟，所以無法根據(jù)自身的應用需求對數(shù)據(jù)做出更精確的操作。本文針對點云數(shù)據(jù)后處理過程，嘗試尋找合理而高效的處理海量點云數(shù)據(jù)的算法來完成點云的自動濾波過程，本文的實驗操作是基于Terrasolid軟件下完成的，通過不同的算法來實現(xiàn)粗差剔除以及濾波分類參數(shù)設計實驗，給出了最佳的濾波效果和精度分析。

2 機載激光雷達對地定位原理

機載激光雷達系統(tǒng)所獲取的數(shù)據(jù)是地面點三維坐標[2]，它是通過一系列的坐標轉(zhuǎn)換求得的。首先利用測距儀測得的斜距P所在的坐標系作為初始坐標系，再結(jié)合慣性導航系統(tǒng)(INS)測定的飛行器的精確姿態(tài)信息，即側(cè)滾角，傾斜角，航向角(φ，ω，κ)進行坐標矩陣旋轉(zhuǎn)；最后聯(lián)合GPS提供的飛行器位置信息，就可以利已知GPS坐標加坐標旋轉(zhuǎn)矢量方法求出每個地面激光腳點的精確三維坐標(x，y，z)。

3 機載LiDAR數(shù)據(jù)后處理

機載LiDAR數(shù)據(jù)后處理主要進行濾波和分類兩個操作。從激光腳點數(shù)據(jù)點云中提取數(shù)字地面高程模型(DEM/DTM)需要將其中的地物數(shù)據(jù)腳點去掉，這就是所謂的激光雷達數(shù)據(jù)的濾波[4]。而若要進行地物提取以及三維重建等操作，還需對激光點云數(shù)據(jù)進行分類(分割)，區(qū)分植被數(shù)據(jù)點、人工地物點。

3.1 點云數(shù)據(jù)濾波方法及實驗

當前國內(nèi)外對機載LiDAR點云數(shù)據(jù)處理算法還在研究發(fā)展階段，現(xiàn)有的用于機載LiDAR測量數(shù)據(jù)的濾波算法都是基于三維激光腳點數(shù)據(jù)的高程突變等信息進行的。主要方法分為移動窗口法[3]，數(shù)學形態(tài)濾波法[5]，迭代線性最小二乘內(nèi)插法[6]，基于地形坡度濾波法[7]，移動曲面擬合法[3]等幾種。每種3方法都存在其自身的缺陷，仍有待改善。

本文基于Terrasolid軟件進行濾波分類處理，該軟件是以Microstaion為平臺開發(fā)使用的。實驗所采用的數(shù)據(jù)是ISPRS第三工作組提供的一組森林地區(qū)的機載LiDAR數(shù)據(jù)。其濾波方法是Axelsson于2000年提出的一種基于不規(guī)則三角網(wǎng)[8](TIN)的漸進加密，該方法尋找到每個格網(wǎng)中的最低點，并將其作為初始種子點，去生成一個稀疏的TIN，再對其他各點進行判斷。這種算法的關鍵性問題在于濾波閾值的選取，主要的濾波參數(shù)有地形坡度角，迭代角，迭代距離，需根據(jù)經(jīng)驗選擇合適的參考值。本算法最大的一個優(yōu)勢就是可以有效地保留地面斷裂特征，處理曲面不連續(xù)的情況，適用于密集分布的城區(qū)。

3.2 利用Terrasolid進行噪聲剔除

常見的噪聲包括低位粗差和高位粗差。低位粗差是由于測量中多路徑誤差導致的極低點，而高位誤差則受到低空飛行物影響提前反射信號產(chǎn)生較大誤差。該算法是用一個點(即中心點)的高程值與給定距離范圍內(nèi)每一個點的高程值比較，如果中心點明顯低于其他點，這個點將被分離出來。如圖1所示，左圖設定參數(shù)剔除的是距離該點5 m之內(nèi)低于一般地面高度0.5 m的點，右圖設定的參數(shù)剔除的是半徑在5 m內(nèi)高位誤差大于5倍標準差的點。

圖1 分離低點算法以及剔除高位點的算法

一般情況粗差剔除前后對地形變化的影響是很大的。對于本實驗結(jié)果，低位點受多路徑影響較嚴重，剔除粗差后地形高程整體提高了。實際上，粗差對后期形成TIN的影響是非常大的，它是后期數(shù)據(jù)濾波的必要保障。

3.3 利用TerraScan進行濾波實驗

在Microstation環(huán)境下打開Terrascan。在實用工具菜單下選擇MDL應用，在打開的窗口中選中TSCAN，安裝并在TerraScan界面下點擊File菜單，讀入數(shù)據(jù)。在主菜單中選擇Classify菜單下的Routine選項→Ground項。

通過選取不同的參數(shù)值，可以得到如表1所示的實驗總結(jié)。利用不同的顏色可以區(qū)別哪些是地面點，哪些是地物點。

不同參數(shù)下的點云濾波 表1

3.4 植被的分層處理實驗

通常情況下，人們要將植被劃分為3種:即低矮植被，中等高度植被，高度植被。植被的分類算法使用的By height from算法，低矮植被距離地面0 m～2 m，中等植被一般距離地面2 m～5 m，高植被距離地面常為5 m～15 m。前面對原始激光腳點進行了濾波，基于它的基礎，就可利用此算法把符合要求的點分類出來，如圖2所示。

圖2 利用height from ground算法分類植被層

圖2 顯示的是從低植被層到中植被層的參數(shù)設置，由本算法得到了其最終的植被分層渲染效果如圖3所示。

圖3 植被分層效果圖

可以看出植被點高程是由矮到高逐步遞增的。藍色區(qū)域表達了地面點層，綠色代表了低矮植被層，黃色代表中度植被層，紅色代表高植被層。這樣的分層可以直觀感受到地形變化和特征，為更好地使用地形提供有利參靠。

4 濾波精度分析

通過設定不同的閾值參數(shù)，可以得到不同的濾波分類效果。此時需要對不同參數(shù)的濾波結(jié)果進行精度評定。評價精度的指標，是通過判斷兩類分類誤差的大小體現(xiàn)的。第Ⅰ類誤差是指拒絕了本屬于地形表面的激光腳點，第Ⅱ類誤差是指接受了不屬于地形表面的激光腳點。具體計算式如表2所示:

兩類誤差計算公式 表2

總誤差=(B+C)/(A+B+C+D)?100%

本文使用C#語言，編程比較實驗所得地面點點云與真實地面點點云數(shù)據(jù)間相應點的個數(shù)，看哪個參數(shù)閾值得到的真實點個數(shù)最多，與真實情況最符合，那么此后的設置就可以將所得閾值作為經(jīng)驗值使用。下面是部分比較程序:

class Point //點云數(shù)據(jù)基本結(jié)構(gòu)

{ public double X；public double Y；public double Z；

｝

static void Compute()

{

//讀取真實地面點數(shù)據(jù)文件，并轉(zhuǎn)化成Point類型列表

var RGP=File.ReadAllLines("RealGroundPoint.txt").Select(p=＞

{

var Pstr=p.Split('').Cast ＜double＞().ToList ()；

return new Point(){

X=Pstr[0]，

Y=Pstr[1]，

Z=Pstr[2]｝；

｝).ToList()；

……

……同理讀取真實地物點數(shù)據(jù)RFP，實驗地面點數(shù)據(jù)TGP，實驗地物點數(shù)據(jù)TFP

//計算真實地面點同實驗地物點交集(共同數(shù)據(jù))，得到第一類錯誤

var FirstTypeError=TFP.Intersect(RGP).Count()/ (double)RGP.Count()；

//計算真實地物點同實驗地面點交集(共同數(shù)據(jù))，得到第二類錯誤

var SecondTypeError=TGP.Intersect(RFP).Count ()/(double)RFP.Count()；

實驗得到的兩類錯誤統(tǒng)計如表3所示:

兩類誤差結(jié)果統(tǒng)計表 表3

從誤差結(jié)果表中可看出，最大建筑物尺寸、最大地形角、迭代角的選擇對地表面構(gòu)建有著很重要的參考價值。實驗1中的參數(shù)，常作為人們長期實踐得到的適應地形條件的標準參數(shù)而使用，選擇稍大一點的迭代角和迭代距離可以有效地保存真實地形特征。其次再進行程序的比較和篩選，就可將真實地面展現(xiàn)出來，減少了兩類錯誤的發(fā)生。

5 結(jié) 語

本文的三個實驗是基于深度理解濾波分類原理后，利用軟件進行算法操作的過程，主要進行了以下幾項操作:粗差剔除，地面點濾波，植被點分類處理。從圖像上來看，每個算法都能夠基本符合現(xiàn)實的情況，把大部分的點分類到相應的層中。但由于Terrasolid軟件源代碼不公開，所以對于數(shù)據(jù)的準確性和真實度還有待商榷。本文基于該軟件開發(fā)了一個精度估算的功能，將實驗結(jié)果與參考數(shù)據(jù)比較，為實驗參數(shù)的設置提供了一個很大的技術參考。機載LiDAR是攝影測量技術的一種有力補充，未來若能融合其他數(shù)據(jù)源進行地物識別，提出更加完美的濾波算法，必將是機載Li-DAR技術的一個革命性轉(zhuǎn)折。

[1] 劉經(jīng)南，張小紅.激光掃描測高技術的發(fā)展與現(xiàn)狀[J].武漢大學學報·信息科學版，2003(2).

[2] 劉經(jīng)南，張小紅.機載激光雷達測高數(shù)據(jù)濾波[J].測繪科學，2004，29(6):158～160.

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

[4] 賴旭東.機載激光雷達數(shù)據(jù)處理中若干關鍵技術的研究[D].武漢:武漢大學，2006.

[5] 崔屹.圖像處理與分析:數(shù)學形態(tài)學方法及應用[M].北京:科學出版社，2002.

[6] Kraus K，Pfeifer N.Determination of Terrain Models in Wooded Areaswith Airborne Laser Scanner Data[J].ISPRS Journal of Photogrammetry and Remote Sensing，1998，53 (4):193～203.

[7] Vosselman G.Slope based filtering of laser altimetry data [A].In:International Archives of Photogrammetry and Remote sensing，vol.XXXIII，Part B3，Amsterdam，2000:935～942.

[8] 史建青，史照良.TerraScan在LIDAR點云濾波中的應用研究[J].淮海工學院學報，2011(3).

Processing of the Laser Radar Point Cloud Data Based on the Terrasolid Software

Li Lianfang，Shao Qiuming
(Chongqing Survey Institute，400020，Chongqing，China)

With the developing of the spatial data application，spatial data acquisition and processing become a hot issue for researchers.In 1970s the airborne LiDAR technology is a new remote sensingmeasurementmethod with functions of active，fast and accurate.This paper expounds the gross error elimination and filter classification method towards the point cloud data through the Terrasolid software.In this process，it chieved the best filer effect by setting different filter-classification parameters for precision comparation.

airborne LiDAR；point cloud；filtering algorithm；ground extraction

1672-8262(2013)05-67-04

P237

B

2013—02—18

李蓮芳(1982—)，女，工程師，主要從事工程測量工作。

十二五國家科技支撐計劃基金(2011BAH12B07-03)