本發(fā)明屬于機(jī)載激光雷達(dá)森林遙感應(yīng)用領(lǐng)域，具體涉及一種基于機(jī)載激光點(diǎn)云聚集關(guān)系的單木分割方法。

背景技術(shù)：

機(jī)載激光雷達(dá)(airbornelightdetectionandranging，a-lidar)是通過搭載在航空飛機(jī)上的傳感器所發(fā)出的激光來測定傳感器與目標(biāo)物之間距離的一種有效且主動(dòng)的遙感技術(shù)，不僅可以精確、快速的獲取地球表面幾何信息(如三維地理坐標(biāo))、回波地物的輻射信息，而且可以應(yīng)用于大尺度上，再加上其具有很強(qiáng)的穿透森林植被的能力,為森林結(jié)構(gòu)參數(shù)的提取以及森林生態(tài)系統(tǒng)功能參數(shù)的反演提供了全新的技術(shù)手段。

在林業(yè)應(yīng)用方面，利用機(jī)載激光雷達(dá)數(shù)據(jù)進(jìn)行研究大致包含森林參數(shù)反演(樹高、冠幅、葉面積指數(shù)等)、單木分割、林木判別以及森林3d模型重建。林木判別是從原始激光點(diǎn)云中識(shí)別出植被回波點(diǎn)的方法；而單木分割算法主要用來解決的問題是在林木判別后得到的植被回波點(diǎn)中分割出屬于特定單棵樹木的點(diǎn)云。在實(shí)際野外測量中，比如lai值、樹高、冠幅等實(shí)測參數(shù)往往都是針對單棵樹木進(jìn)行人工測量的，而在利用a-lidar或者其他遙感數(shù)據(jù)反演森林參數(shù)過程中，基本上都是以這些單木實(shí)測數(shù)據(jù)源作為驗(yàn)證數(shù)據(jù)，而且單木分割的位置精度、點(diǎn)云數(shù)據(jù)的密度等也影響著森林參數(shù)反演的精度,因而單木分割這個(gè)步驟往往是不能忽視的。

機(jī)載激光雷達(dá)數(shù)據(jù)不僅提供了高精度的地理位置數(shù)據(jù)，并且采用“自上而下”的方式掃描植被冠層，得到的回波數(shù)據(jù)中包含描述森林冠層的垂直分布特性的信息，使得在森林樣方下的單木分割成為可能。

目前國際上存在的單木分割方法中，最常用的方法是根據(jù)植被冠層的冠層高度模型(canopyheightmodel，chm)進(jìn)行單木分割，大致過程是它首先通過建立的數(shù)字表面模型(digitalsurfacemodel，dsm)與數(shù)字高程模型(digitalelevationmodel，dem)相減得到chm，然后從chm中定位出局部最大值，并自定義一個(gè)范圍來擴(kuò)展至被識(shí)別單棵樹的冠幅，以此冠幅為基礎(chǔ)提取出該棵樹。使用chm方法進(jìn)行單木分割時(shí)，由于建立的chm常常會(huì)使用插值點(diǎn)代替真實(shí)的三維空間點(diǎn)，會(huì)造成不可逆的信息丟失，再加上自定義的窗口大小對chm進(jìn)行處理，窗口大小與樹高和冠幅大小相聯(lián)系，因而在利用chm進(jìn)行分割時(shí)往往造成錯(cuò)誤分割以及不適當(dāng)?shù)男螤罘指睢?/p>

而在三維空間上直接對原始點(diǎn)云進(jìn)行分割發(fā)展的越來越好，比如h.lee等人(2010)，根據(jù)局部最大值的點(diǎn)與點(diǎn)之間的距離提出可適應(yīng)性聚類方法直接對植被點(diǎn)云數(shù)據(jù)進(jìn)行聚類，使其單木分割的總體精度達(dá)到了95.1％。lu等人(2014)基于強(qiáng)度和激光點(diǎn)云數(shù)據(jù)的在空間分布的3d結(jié)構(gòu)，提出了先識(shí)別樹干的自低向上的方式分割闊葉林，結(jié)果表明該方法能從原始點(diǎn)云中檢測出84％的樹木來，而在檢測出的樹木中進(jìn)行分割的正確率達(dá)到了97％。li等人(2012)使用小光斑離散回波激光雷達(dá)點(diǎn)云數(shù)據(jù)，通過考慮兩棵樹之間的相對間隔，以及形狀指數(shù)來判斷激光點(diǎn)，結(jié)果表明該方法能從原始點(diǎn)云中檢測出86％的樹木來，正確分割率達(dá)到了94％。

技術(shù)實(shí)現(xiàn)要素：

針對上述存在問題或不足，本發(fā)明提供了一種基于機(jī)載激光點(diǎn)云聚集關(guān)系的單木分割方法，基于原始點(diǎn)云數(shù)據(jù)，根據(jù)單棵樹內(nèi)的機(jī)載激光點(diǎn)云在三維空間上的聚集關(guān)系以及鄰近的樹與樹之間的點(diǎn)云在三維空間上的聚集關(guān)系進(jìn)行單木分割，提供了單木分割的新路徑。

本發(fā)明具體技術(shù)方案如下：

步驟1、林木判別：

由于單木分割算法只是針對植被點(diǎn)云數(shù)據(jù)的，把回波點(diǎn)分類為植被回波點(diǎn)、地面回波點(diǎn)和其余回波點(diǎn)。

步驟2、局部最大值法初始分割：

在未編碼(默認(rèn)編碼為0)的原始點(diǎn)云中尋找回波點(diǎn)類別屬于植被的且坐標(biāo)z值最大的點(diǎn)，并記錄該點(diǎn)的坐標(biāo)序號，同時(shí)把該點(diǎn)編碼設(shè)為k(k>＝1)，用半徑為r(實(shí)測目標(biāo)樹木數(shù)據(jù)優(yōu)勢植被的最小冠幅)的圓標(biāo)識(shí)該區(qū)域，并且把該區(qū)域內(nèi)的所有植被回波點(diǎn)都設(shè)為編碼k；然后在原始點(diǎn)云中去除已經(jīng)編碼了的點(diǎn)，樹編碼號變?yōu)閗+1，以此循環(huán)下去，直到原始點(diǎn)云中所有植被點(diǎn)的編碼都不為0。

步驟3、刪除錯(cuò)誤編碼樹：

當(dāng)同一編碼樹的點(diǎn)云數(shù)量小于激光點(diǎn)云密度或兩顆樹之間的間隔小于3r，修改該樹的編碼為0。

步驟4、聚類：

a、使用步驟3得到編碼數(shù)據(jù)，依據(jù)樹編碼號選取任意兩顆樹并確定兩樹之間的間隔i。

b、if判斷，如果i小于maxr，maxr為實(shí)測目標(biāo)樹木數(shù)據(jù)優(yōu)勢植被的最大冠幅半徑,則依次執(zhí)行①、②、③、④。

①、由于在樹冠邊緣的不規(guī)則，高樹相對矮樹樹高的邊緣植被點(diǎn)與其頂點(diǎn)形成的垂直夾角θ是不同的；首先比較兩樹的高矮，以低樹的中心點(diǎn)為原點(diǎn)，不考慮高程坐標(biāo)z，建立笛卡爾直角坐標(biāo)系，在高樹頂點(diǎn)與低樹頂點(diǎn)連線上的點(diǎn)云作為高樹的邊緣植被點(diǎn)(如圖4所示)。②、求最小高樹邊緣植被激光點(diǎn)與高樹頂點(diǎn)形成的閾值正切值mintanθ(原理如圖3),即以mintanθ為閾值tanθ。

③、記錄與mintanθ比較的激光點(diǎn)：

為了減少聚類計(jì)算時(shí)間以及矮樹邊緣異值點(diǎn)對聚類結(jié)果的影響：記錄兩樹頂點(diǎn)連線為長邊長度及方向，矮樹樹冠幅為寬邊的平行四邊形所屬激光點(diǎn)point2peek，且該點(diǎn)的編碼不為高樹編碼的點(diǎn)(如圖4所示)。

④、根據(jù)多木判別原理進(jìn)行聚類：首先判斷正切值，如果步驟③中記錄的激光點(diǎn)對應(yīng)的正切值大于mintanθ，判斷該點(diǎn)編碼對應(yīng)的矮樹為高樹的一部分，并把該矮樹對應(yīng)的點(diǎn)編碼改為高樹的編碼，將平行四邊形內(nèi)所屬激光點(diǎn)遍歷(原理如圖3)。

⑤、根據(jù)單木歸屬原理進(jìn)行聚類：如果point2peek對應(yīng)的正切值小于mintanθ，那么判斷該點(diǎn)為高樹的冠幅點(diǎn)，修改該點(diǎn)的編碼為高樹的編碼，遍歷直到修改完該棵樹的所有回波點(diǎn)(原理如圖2)。

d、重復(fù)執(zhí)行a、b步驟，直到遍歷完所有組樹。

步驟5：最近鄰點(diǎn)法處理，把未被編碼的激光點(diǎn)分到距離最近的已編碼聚類的樹中。

步驟④和步驟⑤的聚類原理

單木歸屬原理：在同一棵樹中，假設(shè)在高度為h時(shí),邊緣植被點(diǎn)與頂點(diǎn)形成的垂直夾角為θ，邊緣植被點(diǎn)以內(nèi)的點(diǎn)與頂點(diǎn)形成的夾角為θ`；在高度低于h時(shí)，冠幅內(nèi)的點(diǎn)與頂點(diǎn)形成的夾角為a、b(如圖2所示)。在同一高度h，θ`大于θ角，在高度小于h時(shí)，a、b小于θ角。若已知高度大于h的點(diǎn)屬于該棵樹的冠層點(diǎn)，那么可以求出在高度為h時(shí)邊緣點(diǎn)與頂點(diǎn)形成的最小角度θ；根據(jù)高度小于h的其余點(diǎn)和頂點(diǎn)形成的夾角a、b與θ的關(guān)系，可以判斷其余點(diǎn)與該樹的所屬關(guān)系。

多木判別原理：在一組樹中，與矮樹頂點(diǎn)同一高度上，高樹的邊緣植被點(diǎn)與其頂點(diǎn)形成的垂直夾角為θ；高樹頂點(diǎn)與矮樹頂點(diǎn)形成的夾角為a；矮樹其余點(diǎn)與高樹頂點(diǎn)形成的夾角為b(如圖3所示)。如果在矮樹中的角a小于θ，且存在部分角b小于θ，那么矮樹不劃歸為高樹冠層部分；否則把矮樹劃歸為高樹。

本發(fā)明首先使用terrascan軟件進(jìn)行林木判別，把原始機(jī)載激光雷達(dá)數(shù)據(jù)分類為地面回波點(diǎn)、植被回波點(diǎn)和其余回波點(diǎn)；然后只針對植被回波點(diǎn)利用局部最值法進(jìn)行初始的單木判別，同時(shí)在此過程中刪除錯(cuò)誤判別的樹；再利用單棵樹內(nèi)的點(diǎn)云以及鄰近的樹與樹之間點(diǎn)云在三維空間上的聚集關(guān)系進(jìn)行聚類處理，從而得到較為精確的分割結(jié)果；最后跟據(jù)最近鄰點(diǎn)法把未分類的回波點(diǎn)劃分到離其最近的單木中，從而提取出整片區(qū)域的樹木來，其流程如圖1所示。

本發(fā)明是通過利用機(jī)載激光雷達(dá)獲取森林植被冠層三維激光點(diǎn)云數(shù)據(jù)，并根據(jù)這些三維點(diǎn)云數(shù)據(jù)在單棵樹內(nèi)的聚集關(guān)系以及鄰近的樹與樹之間的聚集關(guān)系，經(jīng)歷林木判別、局部最大值法初始分割、刪除錯(cuò)誤的分割單樹、循環(huán)聚類和最近鄰點(diǎn)法處理，建立基于原始三維激光點(diǎn)云數(shù)據(jù)的從森林中分割出單棵樹的單木分割方法。

綜上所述，本發(fā)明相比chm單木分割方法，能避免信息丟失、減少錯(cuò)誤分割和不適當(dāng)?shù)男螤罘指?，分割效果?yōu)異，并且計(jì)算時(shí)間也得到了大幅度的縮減，能快速、準(zhǔn)確的對森林激光點(diǎn)云數(shù)據(jù)分割出屬于特定單棵樹的回波激光點(diǎn)來，有助于單樹的參數(shù)反演(如葉面積指數(shù)、樹高、生物量等)、單樹的三維結(jié)構(gòu)特征重建以及對森林植被參數(shù)反演的驗(yàn)證。

附圖說明

圖1為本發(fā)明的流程示意圖；

圖2為單木歸屬原理圖；

圖3為多木判別原理圖；

圖4為聚類過程中的確定兩棵樹的相對方向以及取值范圍；

圖5為森林樣方林木判別結(jié)果圖；

圖6為本發(fā)明的操作界面圖；

圖7為初始識(shí)別和刪除錯(cuò)誤編碼樹后的結(jié)果圖；

圖8為單木分割最終結(jié)果圖；

圖9為使用實(shí)測樣方內(nèi)的林木棵數(shù)來驗(yàn)證分割精度的結(jié)果圖；

圖10為使用實(shí)測樹高來驗(yàn)證分割精度的結(jié)果圖。

具體實(shí)施方式

以下通過一個(gè)樣方的實(shí)施并結(jié)合附圖對本發(fā)明作進(jìn)一步詳細(xì)說明：

開發(fā)環(huán)境microsoftvisualstudio2010、開發(fā)平臺(tái)框架.netframework3.5、三維可視化引擎directx11、編程語言c#。

步驟1，以內(nèi)蒙古呼倫貝爾市西北部，額爾古納市東南部的上庫力農(nóng)場為研究區(qū)域，研究區(qū)域地理位置坐標(biāo)為東經(jīng)120°36′50.48″～120°52′56.53″，北緯50°21′11.08″～50°24′32.00″，并選取實(shí)地測量的30個(gè)樣方作為單木分割驗(yàn)證區(qū)，其優(yōu)勢植被為白樺樹。利用leicalas60掃描系統(tǒng)獲取該區(qū)域的機(jī)載激光雷達(dá)數(shù)據(jù)，航飛時(shí)間在2012年8月16日至2012年9月25日，系統(tǒng)具體參數(shù)參見表1。并基于terrascan軟件把原始激光回波點(diǎn)分類為植被、地面、其余回波點(diǎn)，分類結(jié)果如圖5。以下步驟為對樣方編號為yf_1，中心坐標(biāo)為(120.770998°，50.357450°)的單木分割詳細(xì)步驟。

表1機(jī)載激光雷達(dá)系統(tǒng)參數(shù)

步驟2，在原始激光點(diǎn)云中尋找局部最大高程值點(diǎn)，對該點(diǎn)進(jìn)行編碼并以優(yōu)勢植被最小冠幅的近似圓標(biāo)識(shí)該區(qū)域(樣方實(shí)測最小冠幅為1m)，以此循環(huán)下去，直到所有的原始點(diǎn)云都被編碼，最后得到初始識(shí)別的樹木棵樹為197(圖6所示)。采用上述局部最值法進(jìn)行單木分割，由于只是使用優(yōu)勢植被的最小冠幅r來標(biāo)識(shí)各個(gè)區(qū)域，而r并不能代表所有樹的冠幅，那么得到的分割結(jié)果將會(huì)存在多分現(xiàn)象、錯(cuò)誤識(shí)別現(xiàn)象以及形狀不規(guī)范性現(xiàn)象等。為了提高分割精度，需要其他的方法對初始分割結(jié)果進(jìn)行后處理。

步驟3，在初始識(shí)別得到的結(jié)果中，如果所識(shí)別的樹包含的激光回波點(diǎn)數(shù)小于5個(gè)點(diǎn)(這個(gè)點(diǎn)數(shù)閾值根據(jù)點(diǎn)云密度和最小冠幅面積而設(shè)定)，就把屬于該棵樹的所有編碼設(shè)為未編碼狀態(tài)，編碼值設(shè)為0。當(dāng)兩棵樹之間的距離小于1.5米時(shí)，取消矮樹的所有編碼，最終得到修改后的樹木棵樹為45棵，如圖6、圖7所示。

步驟4，在實(shí)測中該樣方單棵樹的最大冠幅為4.2m，這里把最大半徑閾值設(shè)為2.1m，根據(jù)激光點(diǎn)云在單樹或者相鄰樹的聚集情況進(jìn)行聚類處理，最后得到的分割棵樹為38棵，而實(shí)測棵樹為33棵。聚類結(jié)果如圖7所示。

步驟5，根據(jù)距離把未被編碼的回波點(diǎn)分到最近的已編碼已聚類的樹中，如圖7所示。

步驟6，精度檢驗(yàn)。分別從單木分割后得到的樣方內(nèi)樹木棵樹和分割后得到的點(diǎn)云數(shù)據(jù)進(jìn)行單木樹高的反演兩個(gè)方面對分割結(jié)果的精度進(jìn)行了檢驗(yàn)。如圖8、圖9所示。在圖8中，對實(shí)地測量的30個(gè)樣方進(jìn)行分割，分割后得到的樹木棵樹的精度分別是，絕對系數(shù)r²為0.7398，均方根誤差rmse為6.845棵。在圖9中，對分割后得到數(shù)據(jù)進(jìn)行了單木提取以及單木樹高反演，樹高反演精度分別是，絕對系數(shù)0.7912，均方根誤差1.894米。

該實(shí)例根據(jù)本發(fā)明提出的方法，對白樺樹樣方的機(jī)載激光雷達(dá)點(diǎn)云數(shù)據(jù)進(jìn)行分析，依照技術(shù)方案所述，實(shí)現(xiàn)了對30個(gè)樣方的單木分割，并對分割后的結(jié)果進(jìn)行了兩種方法的精度檢驗(yàn)，精度檢驗(yàn)結(jié)果分別為(r²＝0.7398，rmse＝6.845棵)、(r²＝0.7912，rmse＝1.894米)，都取得了較好的精度，充分說明該方法的可適用性。由檢驗(yàn)結(jié)果可知，在一般情況下，分割后的單木棵樹比實(shí)測的稍多，是因?yàn)樵诜指顣r(shí)，如果單棵樹存在比較大的分支，那么單棵樹內(nèi)的回波點(diǎn)云聚集情況將會(huì)改變?yōu)橄噜彉涞狞c(diǎn)云聚集情況，造成多分，尤其是多枝丫的樹出現(xiàn)多分的頻率會(huì)更高。樹高反演的均方根誤差稍微顯大，是因?yàn)闃涓叻囱莸臄?shù)據(jù)來源是經(jīng)過單木分割后得到的點(diǎn)云數(shù)據(jù)，單木分割中本身就存在一定的誤差(如錯(cuò)誤分割、少分、多分等)，這樣就導(dǎo)致單木分割中的誤差傳遞到單木樹高反演中，進(jìn)一步帶來偏差。該分割算法更適用于不存在大的分支的樹種(針葉林)或者樹木冠層形狀比較規(guī)范的闊葉林等。