基于LIDAR數(shù)據(jù)的數(shù)字高程模型獲取

　　為了創(chuàng)建地面模型和建筑物模型而進(jìn)行的LIDAR數(shù)據(jù)的預(yù)處理包含了粗差剔除，數(shù)據(jù)的分類和分割。粗差實(shí)際上指的就是一些錯(cuò)誤的數(shù)據(jù)點(diǎn)，在數(shù)據(jù)獲取時(shí)由于不可預(yù)知的原因而造成，應(yīng)用之前應(yīng)當(dāng)先將這些壞點(diǎn)去掉。機(jī)載激光掃描測(cè)高激光腳點(diǎn)的分布并不規(guī)則，在三維空間的分布形態(tài)呈現(xiàn)隨機(jī)離散的數(shù)據(jù)點(diǎn)“點(diǎn)云”（Point cloud）。在這些點(diǎn)中，有些點(diǎn)位于真實(shí)地形表面有些點(diǎn)位于人工建筑物（房屋、煙囪、塔、輸電線等）或自然植被（樹、灌木、草）。從激光腳點(diǎn)數(shù)據(jù)點(diǎn)云中提取數(shù)字地形高程模型（DTM、DEM）需要將其中的地物數(shù)據(jù)點(diǎn)去掉，這就是所謂的激光測(cè)高數(shù)據(jù)的濾波。這里借用了數(shù)字信號(hào)處理中濾波的概念，即把地形表面當(dāng)作信號(hào)，而將地物（建筑物、樹木等）當(dāng)作噪聲，濾波算法就是從DSM得到DTM的過(guò)程。研究如何從數(shù)據(jù)點(diǎn)云中分離出地形表面激光腳點(diǎn)數(shù)據(jù)子集以及區(qū)分不同地物（包括房屋、道路、植被等）激光腳點(diǎn)數(shù)據(jù)子集，就是所謂的數(shù)據(jù)過(guò)濾和分類。

　　1、LIDAR數(shù)據(jù)濾波的原理

　　濾波的基本原理是基于鄰近激光腳點(diǎn)間的高程突變（局部不連續(xù)），一般不是由地形的陡然起伏所引起，更為可能的是較高點(diǎn)位于某些地物。即使高程突變是由地形變化所引起的，就一個(gè)區(qū)域來(lái)講，其表現(xiàn)形態(tài)也不會(huì)相同，陡坎只引起某個(gè)方向的高程突變，而房屋所引起的高程突變?cè)谒膫€(gè)方向都會(huì)形成階躍邊界。在同一區(qū)域，一定范圍大小內(nèi)地形表面激光腳點(diǎn)的高程和鄰近地物（房屋，樹木，電線桿等）激光腳點(diǎn)高程變化顯著，在房屋邊界處更為明顯。局部高程不連續(xù)的外圍輪廓就反映了房屋的形狀。當(dāng)激光掃描到枝葉繁茂的參天大樹時(shí)，激光腳點(diǎn)間的高程也會(huì)出現(xiàn)局部不連續(xù)的情況，但其表現(xiàn)形態(tài)卻與前者有顯著差異。

　　兩臨近點(diǎn)間的距離越近，兩點(diǎn)高差越大，較高點(diǎn)位于地形表面的可能性就越小，因此，判斷某點(diǎn)是否位于地形表面時(shí)，要顧及該點(diǎn)到參考地形表面點(diǎn)的距離，隨著兩點(diǎn)間距離的增加，判斷的閾值（threshold）也應(yīng)放寬，主要是為了同時(shí)考慮地形起伏產(chǎn)生的高程變化。兩地面點(diǎn)間的距離越遠(yuǎn)，自然高差（地形變化形成的高差）就會(huì)越大。

　　2、現(xiàn)有的濾波方法

　　大部分激光掃描點(diǎn)濾波算法的思想基于如下假設(shè)：臨近激光腳點(diǎn)間的局部不連續(xù)（高程突變）一般不是由地形的起伏所造成，而是突變點(diǎn)位于非地面之上。高程突變的表現(xiàn)形態(tài)也不相同，例如陡坎只引起某個(gè)方向的高程突變，而房屋所引起的高程突變?cè)谒膫€(gè)方向都會(huì)形成階躍邊界。當(dāng)激光束投射到植被上時(shí)，激光點(diǎn)間的高程也會(huì)出現(xiàn)局部不連續(xù)的情況，但其表現(xiàn)形態(tài)與前者有顯著的差異。兩鄰近點(diǎn)間距離越近，兩點(diǎn)高差越大，較高點(diǎn)位于地形表面的可能性就越小。因此，判斷某點(diǎn)是否位于地形表面時(shí)，要顧及該點(diǎn)到參考地形表面點(diǎn)的距離，隨著兩點(diǎn)間距離的增加，判斷的閾值也應(yīng)放寬?？紤]到地形起伏產(chǎn)生的高程變化，兩地面點(diǎn)之間的距離越遠(yuǎn)，可能出現(xiàn)的高差會(huì)越大。目前用于機(jī)載LIDAR數(shù)據(jù)濾波的方法概括來(lái)講大致可以分為形態(tài)學(xué)濾波法，移動(dòng)窗口法，線性預(yù)測(cè)法，基于地形坡度濾波，移動(dòng)曲面擬合濾波方法，三角網(wǎng)迭代法，基于數(shù)據(jù)分割等幾種方法。

　　3、基于移動(dòng)曲面擬合濾波算法的濾波實(shí)驗(yàn)

　　本文采用移動(dòng)曲面擬合濾波算法進(jìn)行多次實(shí)驗(yàn)，包括平坦地形、水域、城區(qū)、山區(qū)等。實(shí)驗(yàn)區(qū)選擇我國(guó)湖北省武漢市地區(qū)，該區(qū)域地物比較豐富，有開(kāi)闊地、道路、居民住宅、大型建筑物、樹、森林和河流等。地形坡度在城市地區(qū)較平緩，部分地區(qū)為山區(qū)或小土墩。LIDAR數(shù)據(jù)由Leica公司的ALS50Ⅱ型航空激光雷達(dá)系統(tǒng)裝載在運(yùn)-5型飛機(jī)上獲得，相對(duì)航高1400m，激光間距大約為2m。

　　3.1 實(shí)驗(yàn)數(shù)據(jù)

　　為了便于直觀的顯示濾波效果，搜集了該地區(qū)的航空影像作為參考。實(shí)驗(yàn)區(qū)主要以工礦廠房為主，包含了4條互相交叉的高級(jí)公路和一小部分住宅。在實(shí)驗(yàn)區(qū)的西北方，有2個(gè)小型的人工湖泊，在西南方和東北方各有一座小山。原始點(diǎn)云數(shù)據(jù)。

　　3.2 濾波實(shí)驗(yàn)結(jié)果

　　首先將獲取的原始點(diǎn)云數(shù)據(jù)*.LAS格式使用LIDAR VIEWER軟件轉(zhuǎn)換為TXT文本格式，將轉(zhuǎn)出的TXT文件使用移動(dòng)曲面擬合算法濾波處理。獲得新的TXT文本。本文采取直接將TXT文本數(shù)據(jù)轉(zhuǎn)化為SHP，即可生成三維效果圖。

　　3.3 實(shí)驗(yàn)分析

　　移動(dòng)曲面擬合濾波算法能夠很好的過(guò)濾掉地面上的建筑物，原來(lái)在實(shí)驗(yàn)區(qū)內(nèi)部的建筑物點(diǎn)已經(jīng)全部濾掉，原建筑物部分的地面表現(xiàn)為空白，從點(diǎn)云分布圖可以看出是因?yàn)樵搮^(qū)域點(diǎn)已經(jīng)被濾掉的緣故。除建筑物外，移動(dòng)曲面濾波算法還過(guò)濾掉了地面的以及道路兩旁的樹木，公路現(xiàn)已經(jīng)表現(xiàn)為條帶的面狀物；西邊的2個(gè)小型人工湖泊和河流表現(xiàn)為洼地，西南部的和東北部的山地效果較為明顯，可以明顯的看到地勢(shì)的起伏。東部的大片建筑物和樹木已經(jīng)完全過(guò)濾，可以看到路基兩側(cè)較過(guò)濾之前平滑，和濾波閾值的選取有直接關(guān)系。閾值的大小，需要根據(jù)不同地區(qū)多次實(shí)驗(yàn)，方可找到合適的閾值。

　　3.4 綜合評(píng)價(jià)

　　移動(dòng)曲面擬合法濾波，算法簡(jiǎn)單明了。運(yùn)算速度快，自適應(yīng)性強(qiáng)，濾波性能基本上不受地形條件和地物數(shù)量的限制。但算法要求保證一定的數(shù)據(jù)密度，以保證趨勢(shì)面更新較快，更新的過(guò)程實(shí)際上就是濾波趨勢(shì)面自適應(yīng)地形起伏的過(guò)程。當(dāng)然算法還要保持局部第行數(shù)據(jù)的離散分布，避免數(shù)據(jù)點(diǎn)的畸形分布（比如數(shù)據(jù)幾乎共線）。對(duì)于當(dāng)前的機(jī)載激光掃描系統(tǒng)來(lái)說(shuō)，數(shù)據(jù)采集密度已經(jīng)不存在問(wèn)題，絕大多數(shù)系統(tǒng)提供的數(shù)據(jù)密度都能保證在1points/m2以上。

　　移動(dòng)曲面擬合算法實(shí)現(xiàn)過(guò)程中的另一關(guān)鍵之處是閾值的選取，這里的曲面擬合算法實(shí)質(zhì)是基于一種外推的算法，在數(shù)學(xué)上外推的精度要遠(yuǎn)遠(yuǎn)低于內(nèi)插的精度，所以閾值的設(shè)置不宜過(guò)大，這樣會(huì)積累外推的誤差，導(dǎo)致曲面無(wú)法移動(dòng)。解決的辦法是可以進(jìn)行多次濾波。

　　4、數(shù)字高程模型的創(chuàng)建

　　DEM是對(duì)一個(gè)地區(qū)的地表高程變化的表示，可以采用多種方式表達(dá)。地形表面的特征決定了地形表面表達(dá)的難度，因而在影響最終DEM表面的各種因素中扮演了重要的角色。隨著計(jì)算機(jī)應(yīng)用技術(shù)和GIS技術(shù)的發(fā)展，人們通常采用數(shù)字地面模型主要有三種基本形式：規(guī)則格網(wǎng)、不規(guī)則三角網(wǎng)和數(shù)字等高線。濾波后的點(diǎn)云數(shù)據(jù)可以認(rèn)定為一個(gè)地面特征點(diǎn)的集合，本文采取了不規(guī)則三角網(wǎng)的表示方法來(lái)表示DEM數(shù)據(jù)。將濾波后的點(diǎn)云先生成三角網(wǎng)進(jìn)而生成高程模型（DEM）。

　　雖然DEM的空間分辨率的最終影響因素是生成DEM的LIDAR數(shù)據(jù)中激光腳點(diǎn)的密度；但是若已經(jīng)存在的DEM的分辨率不能滿足應(yīng)用要求，則可通過(guò)重采樣來(lái)補(bǔ)充采樣點(diǎn)數(shù)量和密度上的不足，以達(dá)到提高空間分辨率的目的。

猜你喜歡：

機(jī)載激光雷達(dá)（LiDAR）點(diǎn)云數(shù)據(jù)結(jié)構(gòu)特點(diǎn)及處理流程

LIDAR技術(shù)的介紹及特點(diǎn)

版權(quán)聲明：文章來(lái)源于網(wǎng)絡(luò)，登載此文出于傳遞更多信息之目的，版權(quán)歸原作者及刊載媒體所有，如本文中圖片或文字侵犯您的權(quán)益，請(qǐng)聯(lián)系我們。