李輝,林沂,孟祥爽,史振偉,蔡萬(wàn)園

基于地基激光雷達(dá)的欒樹(shù)分形特征分析

李輝,林沂*,孟祥爽,史振偉,蔡萬(wàn)園

北京大學(xué)地球與空間科學(xué)學(xué)院遙感與地理信息系統(tǒng)研究所空間信息集成與3S工程應(yīng)用北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100871

分枝結(jié)構(gòu)的幾何特征是樹(shù)木分枝按照一定規(guī)律組合而成的形態(tài)，體現(xiàn)了個(gè)體樹(shù)木在生理、生態(tài)角度對(duì)于環(huán)境的適應(yīng)能力，以及對(duì)光、水、養(yǎng)分、空間等生長(zhǎng)影響因子的競(jìng)爭(zhēng)能力，是生態(tài)領(lǐng)域研究的一個(gè)重要內(nèi)容。地基激光雷達(dá)（Terrestrial Laser Scanning, TLS）技術(shù)的發(fā)展為獲取單木結(jié)構(gòu)的參數(shù)提供了一種便捷、無(wú)損的綠色方式。本文使用基于TLS點(diǎn)云的方法來(lái)代替耗費(fèi)大量人力的方式來(lái)獲取樹(shù)木分枝的結(jié)構(gòu)參數(shù)，并進(jìn)行分形特征的分析，發(fā)現(xiàn)在欒樹(shù)的生長(zhǎng)發(fā)育階段，為了充分利用光照條件和生長(zhǎng)空間，隨著DBH的增加，分枝率和枝徑比同步上升，但分枝角度呈先上升后下降的趨勢(shì)。

地基激光雷達(dá); 欒樹(shù); 特征分析

樹(shù)木可以被看作各種組件組合而成的集合體[1]，不同分枝特征及樹(shù)葉、芽的動(dòng)態(tài)變化反映了樹(shù)木對(duì)光、水、養(yǎng)分與空間等資源的利用[2]，表達(dá)了樹(shù)木在不同生長(zhǎng)階段為了適應(yīng)環(huán)境的最終選擇的應(yīng)對(duì)策略[3]。對(duì)于樹(shù)木分形特征的研究，是為了了解樹(shù)木對(duì)于環(huán)境空間的利用能力，探討生長(zhǎng)環(huán)境與個(gè)體生長(zhǎng)之間的相互關(guān)系[4,5]。以往的研究表明，樹(shù)木的分形特征和冠層構(gòu)型反映了其對(duì)生存環(huán)境的妥協(xié)策略[6,7]。而經(jīng)過(guò)幾十年的發(fā)展，人們對(duì)不同樹(shù)木構(gòu)型和分枝特征的認(rèn)識(shí)和研究已逐步深入[8]。孫書(shū)存等通過(guò)實(shí)測(cè)數(shù)據(jù)分析了遼東栓皮櫟的分形特征[9]。張文輝等通過(guò)分析沙參的分枝特征，發(fā)現(xiàn)該理論可用于研究瀕危種群的分布[10]。馬克明等利用分形理論研究了東北草原格局，并提出了計(jì)算分支格局的模型[11]。趙相健等研究發(fā)現(xiàn)，競(jìng)爭(zhēng)地位不同的太白紅杉的分枝格局差異性與地形具備較強(qiáng)的關(guān)聯(lián)性[12]。林勇等通過(guò)分析不同離海距離下的木麻黃分支格局，得出持續(xù)性風(fēng)力對(duì)其影響較大，冠層形狀和分枝長(zhǎng)度都體現(xiàn)了對(duì)環(huán)境的響應(yīng)特征[13]。但由于測(cè)量方式的限制，該領(lǐng)域的研究也多基于草本植物，雖在之后擴(kuò)大到高大喬木，但依然存在樣本數(shù)據(jù)較少，難以大規(guī)模進(jìn)行研究的問(wèn)題。

地基激光雷達(dá)系統(tǒng)又稱(chēng)地面激光掃描（Terrestrial Laser Scanning, TLS）系統(tǒng)[14,15]，可通過(guò)發(fā)射激光脈沖獲取目標(biāo)表面的三維坐標(biāo)[16]。地基激光雷達(dá)掃描并結(jié)合相應(yīng)算法是獲取樹(shù)木的三維結(jié)構(gòu)的便捷方式，為實(shí)現(xiàn)可重復(fù)，無(wú)創(chuàng)且客觀的獲取樹(shù)木分枝結(jié)構(gòu)的提供了新的方式[15,17]。利用TLS掃描能夠便捷的估算樹(shù)木的結(jié)構(gòu)參數(shù)[18,19]，使用一些算法和定量結(jié)構(gòu)模型，樹(shù)木胸徑[20,21]、樹(shù)高[22-24]、樹(shù)冠直徑和面積[25,26]可以從單木點(diǎn)云中估算出來(lái)。Dassot M等獲取了42個(gè)不同樹(shù)種、不同大小的樣品的地基激光雷達(dá)數(shù)據(jù)，從點(diǎn)云數(shù)據(jù)中選取基徑大于7 cm的樹(shù)枝，通過(guò)柱面擬合方法進(jìn)行基徑預(yù)測(cè)建模[27]。R?der A等對(duì)13 m以下的31棵針葉樹(shù)每輪的枝高和枝數(shù)進(jìn)行了自動(dòng)識(shí)別和估計(jì)[28]。Raumonen P等提出了一種模型可以近似地得到整棵樹(shù)或部分樹(shù)的分枝結(jié)構(gòu)和大小特性[29]。Lau A等首次利用定量結(jié)構(gòu)模型提取樹(shù)木分枝參數(shù)[30]，Chen S等利用QSM（Quantitative Structure Model，QSM）模型進(jìn)行生物量估算[31]。倪文健等基于樹(shù)干點(diǎn)云垂直連續(xù)分布的特征提出了一種樹(shù)干識(shí)別方法[32]。王寧寧等由點(diǎn)云數(shù)據(jù)得到楊樹(shù)的樹(shù)冠特征參數(shù)并進(jìn)行分析[33]。劉金鵬等基于Skeltre算法提出一種利用柱形擬合法修正分支點(diǎn)處分支延伸方向的新方法[34]。截止目前，利用TLS點(diǎn)云對(duì)分枝分形特征的研究還比較少，本文旨在使用TLS點(diǎn)云的來(lái)代替耗費(fèi)人力的方式來(lái)獲取樹(shù)木分枝的結(jié)構(gòu)參數(shù)，并進(jìn)行分形特征的研究。

欒樹(shù)俗稱(chēng)“燈籠樹(shù)”，屬于落葉喬木，耐寒、耐旱[35]。在我國(guó)多作為園林綠化樹(shù)木，尤其在北方地區(qū)分布較廣[36]，同時(shí)其根系發(fā)達(dá)，抗污染性強(qiáng)，具備修復(fù)工業(yè)污染土地的能力[37]，葉子可作為染料的原材料[38]，欒樹(shù)花不僅可以藥用還能做油料[39]，具有較高的生態(tài)價(jià)值和經(jīng)濟(jì)價(jià)值。因此對(duì)于欒樹(shù)的相關(guān)研究不僅具備經(jīng)濟(jì)意義，在目前追求低碳環(huán)保、建立綠色中國(guó)的大背景下，還具備較強(qiáng)的生態(tài)學(xué)意義。本文以北京市昌平區(qū)的40棵10年生綠化欒樹(shù)為研究對(duì)象，同時(shí)獲取點(diǎn)云數(shù)據(jù)和人工測(cè)量數(shù)據(jù)，通過(guò)QSM定量結(jié)構(gòu)模型提取每棵樹(shù)木的分枝長(zhǎng)度、直徑等結(jié)構(gòu)參數(shù)，驗(yàn)證QSM算法可靠性的同時(shí)，針對(duì)欒樹(shù)研究其分枝格局和特征。

圖 1 技術(shù)路線圖

1 材料與方法

1.1 人工測(cè)量數(shù)據(jù)

研究區(qū)域位于北京市昌平區(qū)，將研究對(duì)象40棵欒樹(shù)按順序依次編號(hào)1至40號(hào)，測(cè)量的參數(shù)主要有枝下木的長(zhǎng)度、直徑，1級(jí)和2級(jí)分枝的長(zhǎng)度、直徑，對(duì)于每一級(jí)分枝按照逆時(shí)針進(jìn)行排序，子分枝與母分枝對(duì)應(yīng)記錄。為了便于統(tǒng)計(jì)，將枝下木定義為0級(jí)分枝，由0級(jí)分枝生長(zhǎng)出來(lái)的為1級(jí)分枝，源自1級(jí)分枝的定義為2級(jí)分枝，同樣，源自于2級(jí)分枝和3級(jí)分枝的定義為3級(jí)分枝和4級(jí)分枝。分枝長(zhǎng)度定義為基部到分枝末端之間的距離。直徑定義為基部直徑和末端直徑的平均值，通過(guò)分別測(cè)量分枝基部和末端的周長(zhǎng)后折算得出。由于隨著分枝層級(jí)的增加，分枝直徑逐漸遞減，人工測(cè)量誤差會(huì)相應(yīng)增大，因此只測(cè)量到2級(jí)分枝。即便如此，40棵樹(shù)木的人工測(cè)量在4人協(xié)作的情況下，共計(jì)耗費(fèi)2 d時(shí)間，且隨著測(cè)量時(shí)間增長(zhǎng)，測(cè)量誤差增大的概率也相應(yīng)增加。

表 1 實(shí)測(cè)數(shù)據(jù)統(tǒng)計(jì)

人工實(shí)測(cè)樹(shù)木的分枝結(jié)構(gòu)參數(shù)統(tǒng)計(jì)見(jiàn)表1。從實(shí)測(cè)數(shù)據(jù)里看到，樹(shù)干（0級(jí)分枝）的直徑多在15 cm左右且分布相對(duì)均勻，1級(jí)分枝和2級(jí)分枝的長(zhǎng)度和直徑波動(dòng)均較大。從0級(jí)分枝到2級(jí)分枝，直徑和長(zhǎng)度的均值均呈現(xiàn)逐級(jí)遞減的特征。

1.2 TLS數(shù)據(jù)獲取

TLS數(shù)據(jù)使用徠卡BLK360地面激光掃描儀獲取。該儀器采用的激光波長(zhǎng)為830 nm，掃描速率360 000/s個(gè)點(diǎn)，掃描速度3 min內(nèi)可完成一站全景掃描，掃描距離0.6～60 m，點(diǎn)云精度距離10 m處的點(diǎn)間距為6 mm，測(cè)距精度10 m范圍內(nèi)誤差為4 mm。具體結(jié)構(gòu)參數(shù)見(jiàn)表2所示。為了獲取最佳效果的點(diǎn)云數(shù)據(jù)，本實(shí)驗(yàn)中每個(gè)測(cè)站距離樹(shù)干約2 m～3 m。采樣工作由1人獨(dú)立完成，共計(jì)耗費(fèi)6 h，相比人工測(cè)量大幅節(jié)約了人力和時(shí)間，且數(shù)據(jù)更加完善。

1.3 TLS數(shù)據(jù)預(yù)處理

點(diǎn)云數(shù)據(jù)預(yù)處理主要包括多站點(diǎn)數(shù)據(jù)的配準(zhǔn)和噪點(diǎn)的去除。掃描獲得的點(diǎn)云數(shù)據(jù)可通過(guò)對(duì)應(yīng)軟件讀取，使用CloudCompare點(diǎn)云處理軟件對(duì)樹(shù)木周邊地物進(jìn)行人工交互去除后，再用CSF濾波去除地面點(diǎn)、SOR濾波清理點(diǎn)云噪聲，并對(duì)每棵樹(shù)的點(diǎn)云進(jìn)行目視檢查，確保去除噪聲時(shí)不丟失樹(shù)木的任何部分?jǐn)?shù)據(jù)。最后進(jìn)行單木分割。

本文分離原始點(diǎn)云中地面點(diǎn)與非地面點(diǎn)時(shí)采用了布料模擬濾波算法（Cloth Simulation Filtering, CSF），算法的概述如圖2所示。該算法是由張吳明等于2016年提出，典型的優(yōu)點(diǎn)是適用性強(qiáng)和形象化，不需要根據(jù)地形來(lái)調(diào)整太多參數(shù)，簡(jiǎn)化地面濾波的處理流程[40]。Yilmaz CS等[41]、張昌賽等[42]通過(guò)將CSF算法[43]與Photoscan、Global Mapper和Lastools等多種濾波算法相比較，得出CSF算法在精度、穩(wěn)定性和適用性上均表現(xiàn)較強(qiáng)的結(jié)論。

1.4 樹(shù)木分枝結(jié)構(gòu)參數(shù)提取的QSM算法

通過(guò)定量結(jié)構(gòu)模型（Quantitative Structure Model, QSM）完成對(duì)欒樹(shù)分枝的重建，從而提取1、2級(jí)分枝的長(zhǎng)度、半徑等結(jié)構(gòu)參數(shù)，并與實(shí)測(cè)數(shù)值進(jìn)行匹配，驗(yàn)證其精度與可靠度的同時(shí)，為進(jìn)一步的分枝結(jié)構(gòu)分析提供數(shù)據(jù)基礎(chǔ)。

基于單木點(diǎn)云，通過(guò)QSM算法對(duì)分枝結(jié)構(gòu)進(jìn)行三維重建。QSM算法通過(guò)一系列圓柱來(lái)擬合分枝，將屬于同一分枝的圓柱長(zhǎng)度之和作為分枝長(zhǎng)度，將屬于同一分枝的第一個(gè)和最后一個(gè)圓柱的直徑平均值作為分枝直徑。完成對(duì)樹(shù)木三維結(jié)構(gòu)的重建之后，基于重建樹(shù)木的分枝長(zhǎng)度與實(shí)測(cè)樹(shù)木分枝長(zhǎng)度來(lái)進(jìn)行分枝匹配，對(duì)于每個(gè)分枝，將長(zhǎng)度與其最接近的重建分枝作為匹配對(duì)象，若重建分枝與實(shí)測(cè)分枝的長(zhǎng)度相對(duì)差異超過(guò)50%，則不對(duì)該分枝進(jìn)行配對(duì)。

研究采用Raumonen P等提出的用于重建樹(shù)木精細(xì)結(jié)構(gòu)的定量模型QSM[29]，是一種利用地面激光掃描得到的樹(shù)干和樹(shù)枝的點(diǎn)云，快速、自動(dòng)地建立高精度樹(shù)木模型的新方法。該方法通過(guò)輸入多站點(diǎn)掃描的一棵樹(shù)的點(diǎn)云，然后建立樹(shù)的柔性圓柱模型，對(duì)樹(shù)的可見(jiàn)部分進(jìn)行魯棒重構(gòu)，并能完整的獲取樹(shù)木拓?fù)浞种ЫY(jié)構(gòu)。主要步驟如圖3所示，其單樹(shù)點(diǎn)云及其重建結(jié)果如圖4所示，圖4（a）為其中1號(hào)欒樹(shù)預(yù)處理后點(diǎn)云，圖4（b）為初步重建結(jié)果。

圖3 定量模型QSM主要流程圖

圖4 QSM重建結(jié)果

1.5 基于TLS的欒樹(shù)分形特征分析

為進(jìn)行欒樹(shù)的分枝特征分析，本文使用改進(jìn)后QSM算法提取前四個(gè)級(jí)別的分枝，針對(duì)這些分枝分別提取長(zhǎng)度、直徑、枝徑比（Ratio of Branch Diameter，RBD）、分枝率（Bifurcation Ratio，BR）和分枝角度等參數(shù)，進(jìn)行分形特征分析。枝徑比的計(jì)算公式如下：RBD=D+1/D(1)

其中RBD為第級(jí)分枝的枝徑比，D+1和D分別為第+1級(jí)和第級(jí)分枝的直徑。

分枝率包括總體分枝率（Overall Bifurcation Ratio，OBR）和逐步分枝率（Stepwise Bifurcation Ratio，SBR），總體分枝率計(jì)算公式如下：=(N-N)/(N-1) (2)

其中N為所有級(jí)別分枝總數(shù)，1為1級(jí)分枝個(gè)數(shù)，N為最后1級(jí)分枝個(gè)數(shù)。

逐步分枝率（SBR）則可根據(jù)下式計(jì)算：SBR=N/N+1(3)

其中SBR為第級(jí)分枝的分枝率，N為第級(jí)分枝的個(gè)數(shù)，N+1為第+1級(jí)分枝的個(gè)數(shù)。

2 結(jié)果與分析

2.1 分枝重建精度評(píng)價(jià)

表3 分枝重建精度評(píng)價(jià)

表3顯示運(yùn)用QSM算法對(duì)不同直徑級(jí)別的分枝重建精度，可以發(fā)現(xiàn)對(duì)于分枝直徑小于5 cm、大于5 cm且小于10 cm和大于10 cm三個(gè)級(jí)別，分枝長(zhǎng)度提取的誤差均較小，且提取結(jié)果均小于實(shí)際分枝長(zhǎng)度，相對(duì)誤差分別為-2.88%、-4.05%和-7.15%。QSM算法對(duì)分枝直徑提取的誤差較大，對(duì)于三個(gè)直徑級(jí)別而言，提取結(jié)果均小于實(shí)際的分枝直徑，且相對(duì)誤差均在50%左右。該結(jié)論也與前人研究相一致，即QSM對(duì)于直徑小于20 cm的分枝直徑估算誤差相對(duì)較大[44]。研究通過(guò)QSM算法完成了欒樹(shù)分枝結(jié)構(gòu)重建，驗(yàn)證了QSM算法在直徑偏?。ㄖ睆叫∮?0 cm）的分枝重建時(shí)，對(duì)直徑產(chǎn)生明顯的低估，且重建直徑均小于實(shí)測(cè)直徑。

2.2 不同分枝級(jí)別的分枝個(gè)數(shù)和長(zhǎng)度統(tǒng)計(jì)分析

表4 不同DBH組別欒樹(shù)各級(jí)分枝個(gè)數(shù)與長(zhǎng)度

圖5 不同DBH級(jí)別欒樹(shù)各級(jí)分枝個(gè)數(shù)

表4展示了不同DBH組別欒樹(shù)各級(jí)分枝的個(gè)數(shù)與長(zhǎng)度，可以發(fā)現(xiàn)對(duì)于不同組別，隨著DBH的增加，1級(jí)分枝的長(zhǎng)度有所增加，而其他級(jí)別的分枝沒(méi)有隨著DBH的變化產(chǎn)生明顯變化。隨著分枝級(jí)別的加深，不同DBH組別的分枝長(zhǎng)度均呈現(xiàn)下降趨勢(shì)。就分枝個(gè)數(shù)而言，隨著DBH的增加，2級(jí)及以后級(jí)別的分枝個(gè)數(shù)均呈現(xiàn)上升趨勢(shì)；而隨著分枝級(jí)別的增加，分枝個(gè)數(shù)呈現(xiàn)先上升后下降的趨勢(shì)，如圖5所示。由此可以得出，研究區(qū)欒樹(shù)的DBH越大，1級(jí)分枝的長(zhǎng)度越長(zhǎng)，同時(shí)也有助于各級(jí)別分枝個(gè)數(shù)的增長(zhǎng)。研究區(qū)欒樹(shù)正處于生長(zhǎng)階段，因此隨著DBH的增長(zhǎng)，1級(jí)分枝的長(zhǎng)度增加，同時(shí)各級(jí)分枝的個(gè)數(shù)也有所上升，以此促進(jìn)上部枝條對(duì)光照條件的競(jìng)爭(zhēng)，從而促進(jìn)欒樹(shù)的生長(zhǎng)發(fā)育。

2.3 不同分枝級(jí)別的枝徑比

表5 不同DBH級(jí)別欒樹(shù)各級(jí)分枝枝徑比

圖6 不同DBH級(jí)別欒樹(shù)各級(jí)分枝枝徑比

表5和圖6展示了不同DBH級(jí)別欒樹(shù)各級(jí)分枝枝徑比的變化，可以發(fā)現(xiàn)對(duì)于不同的DBH組別，在前四級(jí)分枝，每一級(jí)分枝的枝徑比隨著DBH的增加均呈現(xiàn)下降趨勢(shì)，由此可以得出，隨著欒樹(shù)發(fā)育階段的深入，后一級(jí)分枝與前一級(jí)分枝的直徑差異在增大。在同一DBH組別內(nèi)，前三級(jí)分枝直徑比逐漸增大，在DBH小于15 cm的范圍內(nèi)，前四級(jí)分枝直徑比均隨著級(jí)別增加而增大，但DBH在15 cm以上的的組別里，4級(jí)分枝的枝徑比均開(kāi)始下降，即3級(jí)分枝的枝徑比達(dá)峰值。通過(guò)對(duì)比不同分枝級(jí)別的枝徑比可以發(fā)現(xiàn)，1級(jí)分枝的枝徑比均較低，后三級(jí)分枝的枝徑比相差不大，隨著分枝級(jí)別的加深，枝徑比呈現(xiàn)先上升后平緩的變化趨勢(shì)。由此可以得出，1級(jí)分枝通常具有較大的直徑，2級(jí)分枝的直徑與1級(jí)分枝直徑間的差異較大，隨著分枝級(jí)別的上升，后一級(jí)別與前一級(jí)別分枝的直徑差距在逐漸減小。

2.4 不同分枝級(jí)別的分枝率

表 6 不同DBH級(jí)別欒樹(shù)各級(jí)分枝率

圖 7 不同DBH級(jí)別欒樹(shù)各級(jí)分枝率

表6和圖7展示了不同DBH級(jí)別欒樹(shù)各級(jí)分枝率變化，可以發(fā)現(xiàn)對(duì)于不同級(jí)別，1級(jí)分枝與其他幾個(gè)級(jí)別的分枝相比具有較高的分枝率，2、3、4級(jí)分枝的分枝率相互差別比較小，但隨著分枝級(jí)別的上升，分枝率總體呈現(xiàn)下降的趨勢(shì)。對(duì)于不同的DBH組別，隨著發(fā)育階段的加深，1、2、3、4級(jí)分枝的分枝率隨著DBH的增加整體呈現(xiàn)上升趨勢(shì)，尤其是1級(jí)分枝的分枝率增長(zhǎng)幅度較大。由此可見(jiàn)，在欒樹(shù)的生長(zhǎng)發(fā)育階段，隨著DBH的增長(zhǎng)，有助于上級(jí)分枝分離出更多的下級(jí)分枝，從而增強(qiáng)對(duì)光照的競(jìng)爭(zhēng)，促進(jìn)樹(shù)木整體的發(fā)育。總體上，隨著DBH的增加，研究區(qū)欒樹(shù)的分枝率呈下降趨勢(shì)。

2.5 不同分枝級(jí)別的分枝角度

表7 不同DBH級(jí)別欒樹(shù)各級(jí)分枝角度

圖8 不同DBH級(jí)別欒樹(shù)各級(jí)分枝角度

分枝角度對(duì)于樹(shù)冠的形態(tài)具有重要作用，樹(shù)木分枝在樹(shù)冠中的分布方向影響著樹(shù)冠對(duì)光照和CO2等條件的利用。圖8和表7展示了不同DBH級(jí)別欒樹(shù)各級(jí)分枝角度的變化，可以發(fā)現(xiàn)對(duì)于不同的分枝級(jí)別，隨著級(jí)別的增加，分枝角度呈現(xiàn)上升趨勢(shì)，到高級(jí)別分枝的變化則較平緩，以便更好地利用生長(zhǎng)空間并吸收光照。對(duì)于不同的DBH組別，能夠發(fā)現(xiàn)隨著DBH的增加，分枝角度呈現(xiàn)先上升后下降的趨勢(shì)。對(duì)于不同的分枝界別而言，DBH在15～17 cm的組別具有最大的分枝角度。由此可以得出，在樹(shù)木生長(zhǎng)的前期分枝角度逐漸增大，隨著樹(shù)木生長(zhǎng)過(guò)程的深入，分枝角度又逐漸變小，這有助于分枝更好地利用空間并吸收光照。

3 討論

本文基于欒樹(shù)的TLS點(diǎn)云使用QSM算法提取前四級(jí)分枝的結(jié)構(gòu)參數(shù)，進(jìn)行了欒樹(shù)分枝的分形特征分析。結(jié)果表明，研究區(qū)的欒樹(shù)正處于生長(zhǎng)階段，隨著發(fā)育階段的深入，欒樹(shù)的DBH增長(zhǎng)，1級(jí)分枝的長(zhǎng)度也隨之增加，同時(shí)各級(jí)分枝的個(gè)數(shù)也有所上升。后一級(jí)分枝與前一級(jí)分枝的直徑差異在增大，枝徑比上升。欒樹(shù)前四級(jí)分枝的分枝率隨著DBH的增加表顯示出上升趨勢(shì)，其中1級(jí)分枝的分枝率具有較大的增長(zhǎng)幅度，分枝角度則隨著DBH的增加表現(xiàn)出先上升后下降的趨勢(shì)。以上分枝格局有助于欒樹(shù)分離出更多的下級(jí)分枝，增強(qiáng)對(duì)光照的競(jìng)爭(zhēng)和對(duì)發(fā)育空間的利用，從而促進(jìn)對(duì)樹(shù)木整體的發(fā)育。探討欒樹(shù)的分枝格局將有助于增進(jìn)其對(duì)光照和生長(zhǎng)空間的利用情況，并揭示其應(yīng)對(duì)環(huán)境變化的響應(yīng)，從而有助于進(jìn)行對(duì)欒樹(shù)的精確管理和科學(xué)經(jīng)營(yíng)。應(yīng)當(dāng)注意的是，本文選擇的欒樹(shù)為同齡10年樹(shù)，均處于生長(zhǎng)發(fā)育階段，不同樹(shù)木之間結(jié)構(gòu)的差異性較小，因此對(duì)比度更高，但是在更大范圍進(jìn)行數(shù)據(jù)采集時(shí)，可能就會(huì)面臨對(duì)比難度增大的問(wèn)題。同時(shí)，盡管QSM算法在分枝參數(shù)提取上有很強(qiáng)的能力，但是實(shí)驗(yàn)證明其在分枝直徑方面的誤差較大，因此后續(xù)研究需要對(duì)如何提升直徑的提取能力以及對(duì)長(zhǎng)度等參數(shù)的精度提高上作出努力。再者，如果可以針對(duì)同采樣區(qū)的欒樹(shù)進(jìn)行年度的跟蹤研究，則可以更加系統(tǒng)科學(xué)的了解不同年齡的欒樹(shù)分形特征以及隨樹(shù)齡增長(zhǎng)所產(chǎn)生的動(dòng)態(tài)變化，后從而在更廣泛的角度討論TLS在進(jìn)行樹(shù)木分枝結(jié)構(gòu)研究方面的應(yīng)用潛力。

4 結(jié)論

本實(shí)驗(yàn)在開(kāi)展過(guò)程中，發(fā)現(xiàn)人工測(cè)量需要人員多、耗時(shí)長(zhǎng)、測(cè)量級(jí)別少的缺點(diǎn)，且如果對(duì)較大樹(shù)木進(jìn)行測(cè)量，非破壞難以完成，而地基激光雷達(dá)相對(duì)節(jié)約大量人力和時(shí)間。本文基于QSM算法的重建結(jié)果，進(jìn)行了欒樹(shù)冠層分形特征的分析，通過(guò)對(duì)比不同組別欒樹(shù)分枝特征的差異，發(fā)現(xiàn)在欒樹(shù)的生長(zhǎng)發(fā)育階段，隨著DBH的增加，欒樹(shù)1級(jí)分枝的長(zhǎng)度上升，同時(shí)具有更多的次級(jí)分枝數(shù)目，分枝率和枝徑比也同步上升。而隨著DBH的增加，各級(jí)分枝角度均呈現(xiàn)先上升后下降的趨勢(shì)，以更充分地利用光照條件和生長(zhǎng)空間，促進(jìn)欒樹(shù)的發(fā)育。因此，本研究證明基于TLS點(diǎn)云能夠?qū)崿F(xiàn)對(duì)欒樹(shù)分枝結(jié)構(gòu)參數(shù)的精確獲取，從而開(kāi)展多級(jí)別分枝的分形特征研究，大幅節(jié)省人工測(cè)量所需要的大量人力物力，同時(shí)增進(jìn)對(duì)欒樹(shù)分枝結(jié)構(gòu)與生態(tài)屬性聯(lián)系的理解，為樹(shù)木的精細(xì)化管理與科學(xué)經(jīng)營(yíng)提供關(guān)鍵的決策依據(jù)。

[1] Ohsawa M, Nitta I. Patterning of subtropical warm-temperate evergreen broad-leaved forests in east Asian mountains with special reference to shoot phenology [J]. Tropics, 1997,6(4):317-334

[2] Wilkerson GG. Plant growth modeling and applications [J]. Agricultural Systems, 2004,82(1):93-94

[3] 孫棟元,趙成義,王麗娟,等.荒漠植物構(gòu)型研究進(jìn)展[J].水土保持研究,2011(5):281-287

[4] 徐程揚(yáng).不同光環(huán)境下紫椴幼樹(shù)樹(shù)冠結(jié)構(gòu)的可塑性響應(yīng)[J].應(yīng)用生態(tài)學(xué)報(bào),2001,12(3):339-343

[5] 李俊清,臧潤(rùn)國(guó),蔣有緒.歐洲水青岡(L.)構(gòu)筑型與形態(tài)多樣性研究[J].生態(tài)學(xué)報(bào),2001(1):151-155

[6] 陳波,宋永昌,達(dá)良俊.木本植物的構(gòu)型及其在植物生態(tài)學(xué)研究的進(jìn)展[J].生態(tài)學(xué)雜志,2002(3):53-57,29

[7] 鐘章成,曾波.植物種群生態(tài)研究進(jìn)展[J].西南師范大學(xué)學(xué)報(bào):自然科學(xué)版,2001,26(2):230-236

[8] 黃文娟,李志軍,楊趙平,等,胡楊異形葉結(jié)構(gòu)型性狀及其與胸徑關(guān)系[J].生態(tài)學(xué)雜志,2010,29(12):2347-2352

[9] 孫書(shū)存,陳靈芝.遼東櫟植冠的構(gòu)型分析[J].植物生態(tài)學(xué)報(bào),1999,23(5):433-440

[10] 張文輝,祖元?jiǎng)?馬克明.裂葉沙參與泡沙參種群分布格局分形特征的分析[J].植物生態(tài)學(xué)報(bào),1999,23(1):31-39

[11] 馬克明,祖元?jiǎng)?興安落葉松分枝格局的分形特征[J].植物研究,2000,20(2):235-241

[12] 趙相健,王孝安.太白紅杉分枝格局的可塑性研究[J].西北植物學(xué)報(bào),2005,25(1):113-117

[13] 林勇明,俞偉,劉奕,等.不同距海處木麻黃分枝格局特征及冠形分析[J].應(yīng)用與環(huán)境生物學(xué)報(bào),2013,19(4):587-592

[14] Grau E. Estimation of 3D vegetation density with Terrestrial Laser Scanning data using voxels. A sensitivity analysis of influencing parameters [J]. Remote Sensing of Environment, 2017,191:373-388

[15] Malhi , Jackson, Bentley,. New perspectives on the ecology of tree structure and tree communities through terrestrial laser scanning [J]. Interface Focus, 2018,8(2):20170052

[16] Bazezew MN, Hussin YA, Kloosterman EH. Integrating Airborne LiDAR and Terrestrial Laser Scanner forest parameters for accurate above-ground biomass/carbon estimation in Ayer Hitam tropical forest, Malaysia [J]. International Journal of Applied Earth Observation & Geoinformation, 2018,73:638-652

[17] Wilkes P, Lau A, Disney M,. Data acquisition considerations for terrestrial laser scanning of forest plots [J]. Remote Sensing of Environment, 2017,196:140-153

[18] Takoudjou SM, Ploton P, Sonk? B,. Using terrestrial laser scanning data to estimate large tropical trees biomass and calibrate allometric models: A comparison with traditional destructive approach [J]. Methods in Ecology & Evolution, 2017,9(4):905-916

[19] Stovall, Vorster, Anderson,. Non-destructive aboveground biomass estimation of coniferous trees using terrestrial LiDAR [J]. Remote Sensing of Environment, 2017,200:31-42

[20] Tansey K, Semles N, Anstee A,. Estimating tree and stand variables in a Corsican Pine woodland from terrestrial. laser scanner data [J]. International Journal of Remote Sensing, 2009,30(19):5195-5209

[21] Huang H, Li Z, Gong P,. Automated methods for measuring DBH and tree heights with a commercial scanning。Lidar [J]. Photogrammetric Engineering & Remote Sensing, 2011,77(3):219-227

[22] Burt A, Disney MI, Armston P,. Rapid characterisation of forest structure from TLS and 3D modelling. in. Geoscience and Remote Sensing Symposium (IGARSS) [J]. IEEE International. 2013, PP: 3387-3390

[23] Krooks A , Kaasalainen S , Kankare V, e. Predicting tree structure from tree height using terrestrial laser scanning. and quantitative structure models [J]. Silva Fenn, 2014, 48: 1125

[24] Brede B, Lau A, Bartholomeus HM,. Comparing RIEGL RiCOPTER UAV LiDAR derived canopy height and DBH with terrestrial LiDAR [J]. Sensors, 2017,17:2371

[25] Zhao F, Guo Q, Kelly M. Allometric equation choice impacts lidar-based forest biomass estimates: A case study from the Sierra National Forest, CA [J]. Agricultural & Forest Meteorology, 2012,165:64-72

[26] Srinivasan S, Popescu SC, Eriksson M, et al. Terrestrial laser scanning as an effective tool to retrieve tree level height, crown width, and stem diameter [J].Remote Sensing, 2015,7(2): 1877-1896

[27] Dassot M, Colin A, Santenoise P,. Terrestrial laser scanning for measuring the solid wood volume, including branches, of adult standing trees in the forest environment [J]. Computers & Electronics in Agriculture, 2012,89:86-93

[28] R?der A,Klemmt HJ., Seifert S,Integration of the forest growth simulator SILVA in practical forest management. EnviroInfo (1). 2007:239-246

[29] Raumonen P, Kaasalainen M, ?kerblom M,. Fast automatic precision tree models from terrestrial laser scanner data [J]. Remote Sensing, 2013,5(2):491-520

[30] Lau A, Martius C., Bartholomeus H, et. Estimating architecture-based metabolic scaling exponents of tropical trees using terrestrial LiDAR and 3D modelling [J]. Forest Ecology and Management, 2019,439:132-145

[31] Chen S, Feng Z., Chen P,. Nondestructive estimation of the above-ground biomass of multiple tree species in boreal forests of China using Terrestrial Laser Scanning [J]. Forests, 2019,10(11):936-961

[32] 倪文儉,過(guò)志峰,孫國(guó)清,等.基于地基激光雷達(dá)數(shù)據(jù)的單木結(jié)構(gòu)參數(shù)提取研究[J].高技術(shù)通訊,2010,20(2):191-198

[33] 王寧寧,尹文廣,黃秦軍,等.三維掃描技術(shù)在獲取楊樹(shù)樹(shù)冠結(jié)構(gòu)特征參數(shù)上的應(yīng)用[J].林業(yè)科學(xué),2015,51(5):108-116

[34] 劉金鵬,張懷清,劉閩,等.基于點(diǎn)云骨架的單木分枝參數(shù)提取研究[J].北京林業(yè)大學(xué)學(xué)報(bào),2016,38(10):15-20

[35] 楊吉華,張永濤,王貴霞,等.欒樹(shù)、黃連木、黃櫨水分生理生態(tài)特性的研究[J].水土保持學(xué)報(bào),2002(4):152-154

[36] 張勝,湯少展,彭兆豐.欒樹(shù)中有效成分的提取、應(yīng)用及其生物修復(fù)綜述[J].江蘇林業(yè)科技,2020,47(2):48-51

[37] 馮獻(xiàn)賓,董倩,李旭新,等.黃連木和黃山欒樹(shù)的抗寒性[J].應(yīng)用生態(tài)學(xué)報(bào),2011,22(5):1141-1146

[38] 馬廣恩,申雅維,魯學(xué)照,等.欒樹(shù)抗菌有效成分的研究[J].中草藥,1998(2):84-85

[39] 翟梅枝,郭景麗,王磊,等.欒樹(shù)花黃酮類(lèi)化合物的提取工藝研究[J].西北林學(xué)院學(xué)報(bào),2010,25(2):136-139

[40] Zhao JQ, Ledoux H, Stoter J,. HSW: Heuristic Shrink-wrapping for automatically repairing solid-based City GML LOD2 building models [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2018,146:289-304

[41] Yilmaz CS, Yilmaz VO, Güng R. Investigating the performances of commercial and non-commercial software for ground filtering of UAV-based point clouds [J]. International Journal of Remote Sensing, 2018:5016-5042

[42] 張昌賽,劉正軍,楊樹(shù)文,等.基于LiDAR數(shù)據(jù)的布料模擬濾波算法的適用性分析[J].激光技術(shù),2018,42(3):410-416

[43] Zhang WM, Qi JB, Wan P,. An easy-to-use airborne LiDAR data filtering method based on cloth simulation [J]. Remote Sensing, 2016,8(6):501

[44] Lau A, Bentley LP, Martius C,. Quantifying branch architecture of tropical trees using terrestrial LiDAR and 3D modelling [J]. Trees, 2018,32(5):1219-1231

Fractal Characteristics ofon Ground-based LiDAR

LI Hui, LIN Yi*, MENG Xiang-shuang, SHI Zhen-wei, CAI Wan-yuan

100871,

The geometric characteristics of tree meristems describes the structures that branches form according to certain laws. It reflects the adaptability of individual tree to the environment, and the competitive power to factors that affect growth such as light, water and growth space.The advancement of Terrestrial Laser Scanning (TLS) technology has made it possible to obtain parameters of single wood constructions in a convenient, non-destructive, and environmentally friendly manner. In this study, we used TLS point cloud based method instead of original labor-intensive methods to estimate structural parameters of tree branches, and perform further analysis of fractal characteristics. The results showed that to make better use of light and growth space during the growth stage, branching ratio and branch diameter ratio rised synchronously as DBH icreasd，while branching angles increased first, and then decreased。

LiDAR; Koelreuteria paniculata; characteristics analysis

P237

A

1000-2324(2022)03-0475-09

10.3969/j.issn.1000-2324.2022.03.021

2021-12-14

2022-02-07

國(guó)家自然科學(xué)基金面上項(xiàng)目(32171782,31870531)

李輝(1983-),男,碩士研究生,講師,主要從事激光雷達(dá)遙感研究. E-mail:Hui_li@pku.edu.cn

Author for correspondence. E-mail:yi.lin@pku.edu.cn