于東海，馮仲科

?

基于無人機傾斜航空影像的樹冠體積測算方法

于東海1,2，馮仲科3※

（1. 甘肅省測繪工程院，蘭州 730000；2. 甘肅省應(yīng)急測繪工程研究中心，蘭州 730000； 3. 北京林業(yè)大學(xué)精準(zhǔn)林業(yè)北京市重點實驗室，北京 100083）

樹冠是結(jié)構(gòu)復(fù)雜的不規(guī)則體，對樹冠體積的精確測定一直是樹木測量研究中的難點問題。該文以消費級多旋翼無人機對目標(biāo)樹木進行傾斜攝影獲取的多角度航空影像為基礎(chǔ)，通過空三加密處理生成目標(biāo)樹木的三維點云模型；用等高線法分割樹冠點云，并確定樹冠最優(yōu)分割層數(shù)；用投影法對點云數(shù)據(jù)進行轉(zhuǎn)化，并選取測算點計算樹高和樹冠任意橫截面積；對分割后各規(guī)則體的體積進行累加獲得樹冠體積。結(jié)果表明：8棵目標(biāo)樹木的樹高測算值相對誤差為1.46%～4.10%，平均相對誤差為2.88%；樹冠體積測算值的相對誤差為6.95%～12.39%，平均相對誤差為9.42%；精度均可滿足林業(yè)調(diào)查中對于樹高和樹冠體積測量結(jié)果的要求。利用無人機傾斜航空影像建立單木的三維點云模型并進行樹冠體積測算的方法是可行且有效的，該方法可為研究單木樹冠幾何參數(shù)的提取提供參考。

無人機；圖像處理；林業(yè)；傾斜攝影；點云數(shù)據(jù)；單木參數(shù)提?。粯涔隗w積

0 引 言

樹冠指樹木在地面以上連同其生枝葉的部分，是樹木形態(tài)的重要組成部分。樹冠體積指樹冠輪廓以內(nèi)所占的體積，是表征樹冠三維空間大小的一項復(fù)合指標(biāo)[1]，也是監(jiān)測樹木生長狀況以及估算樹木生物量的重要依據(jù)。由于樹冠具有結(jié)構(gòu)復(fù)雜和形態(tài)不規(guī)則的特點，而樹冠體積又是一個無法通過直接測量獲取的復(fù)合因子，因此對于其值的精準(zhǔn)測定較為困難[2]。

為了能準(zhǔn)確高效地獲取測量數(shù)據(jù)，研究人員嘗試將測繪儀器和設(shè)備應(yīng)用于樹木測量中[3]，通過建立相應(yīng)的算法[4-6]，獲取樹木的多種幾何參數(shù)信息?，F(xiàn)有的樹冠體積測量方法，按作業(yè)方式不同可分為傳統(tǒng)手工量測法、全站儀測量法和三維激光掃描法3種。傳統(tǒng)手工量測法將冠幅和樹高等數(shù)據(jù)帶入已有經(jīng)驗?zāi)Ｐ椭泄浪鉡7]而獲得樹冠體積，雖然工作效率較高，但作業(yè)方式粗獷且結(jié)果無精度保證。三維激光掃法獲取的樹冠體積雖然精度較高，但設(shè)備價格昂貴體積龐大不便于野外作業(yè)[8-10]，且掃描獲取的大量點云數(shù)據(jù)[11-13]處理難度較大、耗時費力[14-15]。全站儀測量法獲取的樹冠體積精度遠高于傳統(tǒng)手工量測法，遜于三維激光掃描法，但對野外工作條件的要求和成本造價[16]都遠低于三維激光掃描法。

近年來，無人機(unmanned aerial vehicle, UAV)航空攝影測量[17]作為一項測繪新技術(shù)，在森林調(diào)查中的應(yīng)用也越來越廣泛。通過對無人機獲取的高分辨率數(shù)字正射影像（digital orthophoto map, DOM）和數(shù)字表面模型（digital surface model, DSM）進行處理[18-19]，提取冠幅、樹高和株數(shù)密度等參數(shù)[20-22]，再建立實地調(diào)查數(shù)據(jù)與影像提取數(shù)據(jù)之間的反演模型[23-25]，實現(xiàn)對森林蓄積量、生物量和碳儲量的估測[26-28]。可見，無人機航空攝影測量技術(shù)能夠提取不同尺度的森林參數(shù)信息，將該技術(shù)應(yīng)用于樹冠體積測量中，與現(xiàn)有3種方法相比具有以下優(yōu)勢：1）能夠快速獲取樹木的全方位、多角度影像，使樹木的細節(jié)紋理得以清晰呈現(xiàn)；2）影像經(jīng)過同名點匹配生成的點云數(shù)據(jù)量遠小于三維激光掃描法獲取的數(shù)據(jù)量，使得數(shù)據(jù)處理、分析更加易于操作；3）點云數(shù)據(jù)帶有空間坐標(biāo)信息，且能夠清晰的展現(xiàn)樹木的真實形態(tài)特征，克服了全站儀測量法不具備數(shù)據(jù)可視化的缺陷，使得各參數(shù)的提取計算更加直觀快捷高效。

本文以消費級多旋翼無人機對目標(biāo)樹木進行傾斜攝影獲取的多角度航空影像為基礎(chǔ)，通過空三加密處理生成目標(biāo)樹木的三維點云模型，并對樹冠點云進行投影和分割處理，提取樹高和樹冠任意橫截面積，再累加求和獲得樹冠體積。旨在為利用無人機傾斜航空影像研究單木樹冠幾何信息的提取提供參考。

1 材料與方法

1.1 試驗區(qū)概況

試驗區(qū)位于北京市海淀區(qū)西北部的鷲峰林場。地理坐標(biāo)為東經(jīng)116°28￠002，北緯39°54′00″，總面積832 hm2，海拔100～1 153 m。四季分明且溫差較大，全年溫度在-21.7～41.6 ℃范圍內(nèi)，年平均降水量628 mm。林場內(nèi)植被豐富茂盛、林木種類較多，主要樹種為油松（）、側(cè)柏（）、毛白楊（）、白蠟（）、榆樹（）等。

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

1.2.1 目標(biāo)樹木選取

通過實地踏勘，綜合考慮無人機飛行作業(yè)條件及樹木的大小規(guī)格，按照不同徑階分布[29]（中徑階13.0～24.9 cm、大徑階25.0～36.9 cm、特大徑階≥37.0 cm，由于小徑階樹木生長茂密，因而未能找到合適飛行作業(yè)的樹木），選取試驗區(qū)內(nèi)8棵不同樹種、樹冠形狀具有代表性的孤立單木作為目標(biāo)樹木（表1）。所選的目標(biāo)樹木須生長狀況良好且與周邊樹木樹冠無交叉和重疊，以確保無人機具有足夠充足的飛行作業(yè)空間。

表1 目標(biāo)樹木的基本信息

1.2.2 無人機數(shù)據(jù)獲取

數(shù)據(jù)獲取時間為2017年2月至4月間天氣晴好、風(fēng)力較小的時段。數(shù)據(jù)獲取平臺為大疆Inspire-1型消費級多旋翼無人機，該機型將Zenmuse-X3型普通數(shù)碼相機與云臺集成為一體，配有主、副2個遙控器，飛行時會自動記錄攝影瞬間像片的POS（position and orientation system）數(shù)據(jù)。無人機及相機的基本參數(shù)如表2所示。

使用該無人機采用圍繞目標(biāo)樹木盤旋上升的飛行方式采集多角度傾斜航空影像（圖1）。作業(yè)時，主遙控器控制無人機飛行狀態(tài)，實現(xiàn)爬升、下降、懸停旋轉(zhuǎn)及前后左右移動；副遙控器控制云臺狀態(tài)，實現(xiàn)相機在水平及垂直方位的多角度轉(zhuǎn)動，并設(shè)置自動曝光間隔為2 s；主、幅遙控器同時協(xié)同配合，實現(xiàn)對目標(biāo)樹木的傾斜航空影像采集。

1）在保證安全距離的情況下，以目標(biāo)樹木的樹干為中心進行近距離低速飛行攝影；

2）在盤旋上升的過程中實時調(diào)整相機姿態(tài)，以確保獲取目標(biāo)樹木不同位置和不同角度的影像；

3）確保在同一高度采集的相鄰影像重疊率超過90%，不同高度采集的影像重疊率超過60%。

依照此方案，手控操縱無人機分別對8棵目標(biāo)樹木進行傾斜攝影，共獲取2 083張有效航空影像及對應(yīng)的POS數(shù)據(jù)。

表2 無人機及相機的基本參數(shù)

圖1 無人機傾斜攝影示意圖

1.2.3 地面實測數(shù)據(jù)獲取

依據(jù)馮仲科等[30]提出的方法，利用南方NTS-372R型免棱鏡全站儀獲取目標(biāo)樹木的實測數(shù)據(jù)。實測數(shù)據(jù)為樹木不同位置處特征點的三維坐標(biāo)。如圖2所示，測量的關(guān)鍵步驟如下：

1）控制點布設(shè)。在目標(biāo)樹木周圍采用獨立坐標(biāo)系布設(shè)控制點：以目標(biāo)樹木為中心，在其周圍選擇相互通視且大致等距的某3點、和作為控制點；設(shè)定某點的三維坐標(biāo)為已知，在點架設(shè)全站儀并后視點，測量兩點間的距離和坐標(biāo)方位角以獲取點的三維坐標(biāo)；用支導(dǎo)線的方式，依次獲取點和點的三維坐標(biāo)。

2）樹冠表面特征點獲取。根據(jù)樹冠的大致形態(tài)，從東、南、西、北、東南、東北、西南、西北8個方向上選取樹冠特征點；特征點的位置需根據(jù)樹冠的伸展情況和彎曲的趨勢進行選取，一般取外凸和內(nèi)凹相交處，即形態(tài)變化的轉(zhuǎn)折點；將全站儀依次架設(shè)在3個控制點上，以碎部測量的方式分別測量獲取這些特征點的三維坐標(biāo)。

依照此方法，分別獲取8棵目標(biāo)樹木的地面實測數(shù)據(jù)。

注: O為已知點；A、B、C為控制點。

1.3 研究方法

如圖3所示，本文對無人機獲取的數(shù)據(jù)進行空三加密處理，生成目標(biāo)樹木的三維點云模型；運用等高線法將樹冠點云進行分割，并確定樹冠最優(yōu)分割層數(shù)；采用投影法將點云數(shù)據(jù)進行轉(zhuǎn)化，并利用測算點分別獲得樹高和樹冠任意橫截面積的測算值；通過對分割后各規(guī)則體的體積進行累加，獲得整個樹冠體積的測算值；以全站儀實測值作為參考值，對樹高和樹冠體積的測算結(jié)果進行精度檢驗。

圖3 無人機數(shù)據(jù)測算樹冠體積技術(shù)路線

1.3.1 無人機數(shù)據(jù)處理

以Pix4D-mapper攝影測量軟件為工具，分別對8棵目標(biāo)樹木的傾斜航空影像進行處理。以其中的6號目標(biāo)樹木為例對主要處理步驟進行如下說明：

1）空三加密

導(dǎo)入影像數(shù)據(jù)和POS數(shù)據(jù)，將參考系設(shè)置為WGS-84/UTM zone 50N坐標(biāo)系，再導(dǎo)入相機檢校參數(shù)文件，完成影像屬性設(shè)置。通過空三加密處理后，軟件自動恢復(fù)攝影瞬間每張影像的方位姿態(tài)。

2）點云生成

為了能在準(zhǔn)確恢復(fù)目標(biāo)樹木三維形態(tài)的同時，又能避免因數(shù)據(jù)量過大而導(dǎo)致結(jié)果難以處理的問題，經(jīng)過反復(fù)試驗，最終確定將點云生成的“圖像比例”參數(shù)設(shè)置為“1/2”，“匹配最低數(shù)值”設(shè)置為“3”，“點密度”設(shè)置為“中等”。

3）點云編輯與輸出

利用人工手動框選的方式將生成的點云數(shù)據(jù)中地面、天空等多余噪聲點刪除，并輸出目標(biāo)樹木的三維點云模型，其多角度視圖如圖4所示。以LAS和XYZ兩種格式對包含坐標(biāo)和顏色信息的點云成果進行輸出，LAS格式的數(shù)據(jù)用于樹冠點云的分割處理，XYZ格式的數(shù)據(jù)用于樹冠點云的投影和篩選等處理。統(tǒng)計結(jié)果顯示，6號目標(biāo)樹木三維點云模型成果的點云平均密度為482.67個/m3，共生成355 350個三維點。

圖4 樹木三維點云模型

1.3.2 樹冠體積測算方法

利用三維點云模型對樹冠體積進行測算的基本思想是，將樹冠分割為多層規(guī)則體來模擬樹冠形態(tài)，然后求得樹干任意高度和樹冠任意橫截面積兩項幾何參數(shù)，再按照規(guī)則體的體積公式分別計算各層的體積，最后進行累加求和獲得樹冠體積的測算值。

1）樹干任意高度的計算

樹干任意高度指樹干上任意兩點間的長度，實際上是解算這兩點間的空間距離。樹冠的三維點云記錄了每個點的物方空間坐標(biāo)，通過獲取不同位置處的空間坐標(biāo)信息，就能實現(xiàn)樹干任意高度的計算。

式中表示高度，m；表示點號，取值為1, 2,…,。

2）樹冠任意橫截面積的計算

樹冠任意橫截面積的計算實際上是解算平面上各邊界點的連線所形成閉合圖形的面積。將樹冠任一特定高度值的三維點云展繪到軸和軸構(gòu)成的平面上，通過獲取各邊界點的平面坐標(biāo)信息，就能實現(xiàn)樹冠任意橫截面積的計算。

式中表示面積，m2；表示點號，取值為2,3,…,。

3）樹冠體積的計算

樹冠體積的計算實際上是解算樹冠分割后各層規(guī)則體的體積，進而累加求和得到樹冠的體積。首先利用式（1）計算分割后各層的高度；其次利用式（2）計算分割后各層的橫截面積；利用分層高度和橫截面積計算各層規(guī)則體的體積：最后將各層規(guī)則體的體積累加求和獲得樹冠體積。

將分割后的樹冠近似為由中間若干個臺體和頂部及底部兩個錐體構(gòu)成，臺體和錐體的計算公式分別為

式中表示截面號，取值為1, …,；其中，當(dāng)表示臺體截面時取值為2, …,1，當(dāng)表示錐體時取值為1或；S、S+1分別表示第層和第+1層的橫截面積，m2；H表示分割后相鄰兩層間的高度，m。

2 試驗與結(jié)果分析

2.1 樹干高度的提取

要從目標(biāo)樹木的三維點云模型上提取樹干任意位置的高度，就要運用投影法將三維空間的點云數(shù)據(jù)轉(zhuǎn)化至一維直線上。以XYZ格式的目標(biāo)樹木點云數(shù)據(jù)為基礎(chǔ)，將軸和軸所構(gòu)成的平面作為投影面，通過投影的方式轉(zhuǎn)化至二維平面上；然后以樹木根部點和頂點所在的直線為基準(zhǔn)，繪制出樹干的一維中心線；再在中心線上任意獲取兩個測算點的高程值，根據(jù)式（1）計算得到此段樹干的高度。

如圖5所示，以6號目標(biāo)樹木為例，在樹干中心線上，提取出樹木根部測算點所在直線的高程值為184.551 m，頂部測算點所在的直線的高程值為200.322 m，樹干第一枝處測算點所在直線的高程值為188.500 m，由此可得到樹高的值為15.77 m，樹干第一枝下高的值為3.95 m，樹冠（樹干第一枝至樹干頂部）高度的值為11.82 m。

圖5 樹干高度的提取

2.2 樹冠體積的計算

1）分割方法

運用等高線法對樹冠進行分割。等高線是對高程相等的相鄰各點所連成的閉合曲線，而目標(biāo)樹木的三維點云模型存儲了所有點的空間坐標(biāo)，因此對樹冠分割的過程實際上是篩選出點云中高程值相等的點的過程。

將LAS格式的目標(biāo)樹木點云數(shù)據(jù)導(dǎo)入到EPS軟件（北京山維科技股份有限公司）中，選取樹冠點云作為分割對象，將樹干第一枝處測算點的高程值作為等高線的起算基準(zhǔn)，通過“點云生成等高線”模塊設(shè)置相應(yīng)的等高距，生成樹冠點云的等高線。等高距即為樹冠的分層高度，等高線即為樹冠的分割線。

2）樹冠橫截面積的提取

樹冠點云被等高線分割后，落在同一等高線上的點都有相同的高程值，篩選出點云中有此高程值的所有點，并投影到以軸和軸所構(gòu)成的平面上。

如圖6所示，在投影后的圖形上，采用人機交互的方式，人為選取最外圍的邊界點作為測算點，選取原則為：①必須是投影后外圍的輪廓點；②要忽略樹冠離散點的影響；③依次連接后所形成的圖形必須包含投影后的大部分點（離散點除外）。

對測算點進行依次連接后所構(gòu)成的閉合平面就是此等高線所形成的橫截面，再根據(jù)式（3）計算得到此橫截面的面積。

圖6 樹冠橫截面積的提取

3）樹冠最優(yōu)分割層數(shù)的確定

在對樹冠分割時，分割層數(shù)的不同會影響樹冠體積的測算精度，因此需要尋找一個最優(yōu)的分割層數(shù)，使得在分割層數(shù)較少時體積測算值又相對準(zhǔn)確。

微積分求體積的思想[31]指出，當(dāng)隨著物體被分割的單元由大變小直至無窮小時，其求得的體積值會穩(wěn)定在一個常數(shù)附近，即有極限；此時的分割單元就是最優(yōu)分割尺寸，該常數(shù)就是物體的體積。依據(jù)此原理，對樹冠進行不同層數(shù)的多次分割，分別計算各次樹冠體積的測算值，并進行分析對比，確定最優(yōu)分割層數(shù)。

以6號目標(biāo)樹木為例，將樹干第一枝處測算點的高程值作為起點，將樹冠點云依次分割為2、3、4、5、6、7、8、9、10、11、12層，其對應(yīng)的等高距依次為6.00、4.00、3.00、2.50、2.00、1.70、1.50、1.40、1.30、1.10、1.00 m，分別計算不同分割層數(shù)時的樹冠體積測算值，結(jié)果如圖7所示。

圖7 不同分割層數(shù)下的樹冠體積測算結(jié)果

由圖7得出，隨著樹冠分割層數(shù)的增加，體積測算值逐漸增加，當(dāng)達到某一臨界值后基本趨于穩(wěn)定。當(dāng)分割層數(shù)＜8層時，隨著分割層數(shù)的增加，樹冠體積測算值也相應(yīng)增大，且波動較大，方差為3020.63 m6；當(dāng)分割層數(shù)≥8層時，隨著分割層數(shù)的增加，樹冠體積測算值基本趨于穩(wěn)定，方差為8.53 m6，且分層等高距由1.50 m遞減至1.00 m，變化范圍較小。綜上，將樹冠的分割層數(shù)確定為8層為最優(yōu)。如圖8所示，6號目標(biāo)樹木的樹冠點云被7條等高線分割為8層。

注: f1、f2、…、f7、f8表示樹冠分層；g1、g2、…、g6、g7表示樹冠橫截面。

2.3 測算結(jié)果精度檢驗

以全站儀實測獲取的樹木特征點的三維坐標(biāo)為基礎(chǔ)，根據(jù)文獻[2]的方法，分別求得8棵目標(biāo)樹木的樹高和樹冠體積的實測值，并以此為參考值，對本文方法測算樹高和樹冠體積的結(jié)果進行精度檢驗，結(jié)果如表3所示。

表3 目標(biāo)樹木測算值的精度檢驗結(jié)果

由表3知，8棵不同規(guī)格目標(biāo)樹木的樹冠最優(yōu)分割層數(shù)具有差異性。在確定樹冠最優(yōu)分割層數(shù)時，均以樹冠點云為對象，根據(jù)冠高進行多次分割，并計算不同分割層數(shù)下的樹冠體積值；樹冠體積值與分層后每層的橫截面積和分層高度有關(guān)，而橫截面積與樹冠形狀有關(guān)，分層高度又與冠高和分層數(shù)有關(guān)；且對于特定的單棵樹木而言，冠高一定，分層高度越小，分層數(shù)也就越多。由于目標(biāo)樹木的樹冠形狀和冠高都各不相同，因而在同樣的分析方法下確定的最優(yōu)分割層數(shù)也就不同。

樹高測算值與參考值之間相差較小，絕對誤差在0.18～0.52 m范圍內(nèi)，相對誤差在1.46%～4.10%范圍內(nèi)，8棵目標(biāo)樹木樹高的平均相對誤差為2.88%；這說明選取的樹高測算點與全站儀實測獲取的樹木頂、底兩端點位位置基本一致，能夠準(zhǔn)確計算出樹高。5號目標(biāo)樹木的樹冠體積測算值最小，為133.23 m3，7號目標(biāo)樹木的樹冠體積測算值最大，為447.65 m3；樹冠體積測算值的絕對誤差在14.14～43.77 m3范圍內(nèi)，相對誤差在6.95%～12.39%范圍內(nèi)，8棵目標(biāo)樹木樹冠體積的平均相對誤差為9.42%；測算樹冠體積的誤差相對較大，這是由于樹冠體積是一個復(fù)合因子，需要多種參數(shù)同時參與計算才能獲得，且各參數(shù)本身具有誤差，在運算時誤差進行傳播積累而造成的。結(jié)果表明，樹高和樹冠體積的測算精度均滿足林業(yè)調(diào)查中對于二者測量結(jié)果的要求。

3 討 論

針對樹冠是不規(guī)則體的特點，本文在利用生成的三維點云模型測算樹冠體積時運用等高線法將其分割為若干個規(guī)則體。分割樹冠點云所形成的平面就是規(guī)則體的上、下兩個橫截面，等高距就是規(guī)則體的高度；橫截面積和高度是通過在分割后的點云數(shù)據(jù)上人工選取測算點，根據(jù)其坐標(biāo)進行計算得到的。無人機傾斜航空影像和POS數(shù)據(jù)經(jīng)過空三加密處理生成的點云數(shù)據(jù)密度較高，且各點均帶有空間坐標(biāo)信息，在選取計算橫截面積和高度的測算點時又采用投影法進行轉(zhuǎn)化，能有效降低選點誤差對計算結(jié)果造成的影響，保證了人工選點的可靠性。

本文提出的利用無人機傾斜航空影像測算樹冠體積的方法可以直接通過獲取的影像重建出三維點云模型并進行參數(shù)提取，且消費級無人機相比三維激光掃描儀和地基激光雷達等設(shè)備成本更低、數(shù)據(jù)獲取速度更快、應(yīng)用前景更為廣闊，將其應(yīng)用到古樹名木調(diào)查保護工作中，最能體現(xiàn)實際使用價值。然而，由于受到無人機對飛行條件要求較高等因素的制約，目前只以孤立單木作為對象進行分析，如何實現(xiàn)結(jié)構(gòu)復(fù)雜林分內(nèi)多株樹木的三維點云模型重建并提高樹冠體積測算精度，有待進一步研究。

4 結(jié) 論

本文提出了一種基于無人機傾斜航空影像建立單木三維點云模型測算樹冠體積的方法。該方法以消費級多旋翼無人機搭載普通數(shù)碼相機獲取目標(biāo)樹木的多角度傾斜航空影像為基礎(chǔ)，根據(jù)傾斜攝影測量原理對影像進行空三加密處理，生成目標(biāo)樹木的三維點云模型，運用等高線法對樹冠點云進行分割，確定樹冠最優(yōu)分割層數(shù)，采用投影法提取樹高和樹冠任意橫截面積，通過計算獲得樹冠體積的測算值，并對結(jié)果進行精度檢驗。

8棵目標(biāo)樹木樹高測算值的平均相對誤差為2.88%，在最優(yōu)分割層數(shù)下計算的樹冠體積測算值的平均相對誤差為9.42%。利用生成的三維點云模型測算樹高和樹冠體積的精度均能滿足林業(yè)調(diào)查的要求。利用無人機傾斜航空影像進行樹冠體積測算是可行且有效的。

[1] 閆飛. 森林資源調(diào)查技術(shù)與方法研究[D]. 北京：北京林業(yè)大學(xué)，2014. Yan Fei. Research of Technology and Methord of Forest Resource Inventory[D]. Beijing: Beijing Forestry University, 2014. (in Chinese with English abstract)

[2] 何誠，馮仲科，袁進軍，等. 基于數(shù)字高程模型的樹木三維體積測量[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(8)：195－199. He Cheng, Feng Zhongke, Yuan Jinjun, et al. Three- dimensional volume measurement of trees based on digital elevation model[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(8): 195－199. (in Chinese with English abstract)

[3] 馮仲科，黃曉東，劉芳. 森林調(diào)查裝備與信息化技術(shù)發(fā)展分析[J]. 農(nóng)業(yè)機械學(xué)報，2015，46(9)：257－265. Feng Zhongke, Huang Xiaodong, Liu Fang. Forest survey equipment and development of information technology[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(9): 257－265. (in Chinese with English abstract)

[4] 趙芳，馮仲科，高祥，等. 樹冠遮擋條件下全站儀測量樹高及材積方法[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(2)：182－190. Zhao Fang, Feng Zhongke, Gao Xiang, et al. Measure method of tree height and volume using total station under canopy cover condition[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactionsof the CSAE), 2014, 30(2): 182－190. (in Chinese with English abstract)

[5] 焦有權(quán)，馮仲科，趙禮曦，等. PSO嵌入SVM算法的活立木材積預(yù)報研究[J]. 光譜學(xué)與光譜分析，2014，34(1)：175－179. Jiao Youquan, Feng Zhongke, Zhao Lixi, et al. Research on living tree volume forecast based on PSO embedding SVM[J]. Spectroscopy and Spectral Analysis, 2014, 34(1): 175－179. (in Chinese with English abstract)

[6] 于東海，馮仲科，曹忠，等. 全站儀測量立木胸徑樹高及材積的誤差分析[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(17)：160－167. Yu Donghai, Feng Zhongke, Cao Zhong, et al. Error analysis of measuring diameter at breast height and tree height and volume of standing tree by total station[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(17): 160－167. (in Chinese with English abstract)

[7] 吳明欽，孫玉軍，郭孝玉，等. 長白落葉松樹冠體積和表面積模型[J]. 東北林業(yè)大學(xué)學(xué)報，2014，42(5)：1－5. Wu Mingqin, Sun Yujun, Guo Xiaoyu, et al. Predictive models of crown volume and crown surface area for Korean Larch[J]. Journal of Northeast Forestry University, 2014, 42(5): 1－5. (in Chinese with English abstract)

[8] Hildebrandt R, Iost A. From points to numbers: A database-driven approach to convert terrestrial LiDAR point clouds to tree volumes[J]. European Journal of Forest Research, 2012, 131(6): 1857－1867.

[9] 劉魯霞，龐勇，李增元，等. 用地基激光雷達提取單木結(jié)構(gòu)參數(shù)：以白皮松為例[J]. 遙感學(xué)報，2014，18(2)：365－377. Liu Luxia, Pang Yong, Li Zengyuan, et al. Retrieving structural parameters of individual tree through terrestrial laser scanning data[J]. Journal of Remote Sensing, 2014, 18(2): 365－377. (in Chinese with English abstract)

[10] Asner G P, Knapp D E, Boardman J, et al. Carnegie Airborne Observatory-2: Increasing science data dimensionality via high-fidelity multi-sensor fusion[J]. Remote Sensing of Environment, 2012, 124: 454－465. doi:10.1016/j.rse.2012. 06.012

[11] Zhen Z, Quackenbush L, Zhang L. Trends in automatic individual tree crown detection and delineation-evolution of LiDAR data[J]. Remote Sensing, 2016, 8(4): 333. doi:10. 3390/rs8040333

[12] 徐偉恒，馮仲科，蘇志芳，等. 一種基于三維激光點云數(shù)據(jù)的單木樹冠投影面積和樹冠體積自動提取算法[J]. 光譜學(xué)與光譜分析，2014，34(2)：465－471. Xu Weiheng, Feng Zhongke, Su Zhifang, et al. An automactic extraction algorithm for indvidual tree crown projection area and volume based on 3D point cloud data[J]. Spectroscopy and Spectral Analysis, 2014, 34(2): 465－471. (in Chinese with English abstract)

[13] 李增元，劉清旺，龐勇. 激光雷達森林參數(shù)反演研究進展[J]. 遙感學(xué)報，2016，20(5)：1138－1150. Li Zengyuan, Liu Qingwang, Pang Yong. Review on forest parameters inversion using LiDAR[J]. Journal of Remote Sensing, 2016, 20(5): 1138－1150. (in Chinese with English abstract)

[14] 郭彩玲，宗澤，張雪，等. 基于三維點云數(shù)據(jù)的蘋果樹冠層幾何參數(shù)獲取[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(3)：175－181. Guo Cailing, Zong Ze, Zhang Xue, et al. Apple tree canopy geometric parameters acquirement based on 3D point clouds[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(3): 175－181. (in Chinese with English abstract)

[15] 王佳，張芳菲，高赫，等. 地基激光雷達提取單木冠層結(jié)構(gòu)因子研究[J]. 農(nóng)業(yè)機械學(xué)報，2018，49(2)：199－206. Wang Jia, Zhang Fangfei, Gao He, et al. Extracting crown structure parameters of individual tree by using Ground- based Laser Scanner[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(2): 199－206. (in Chinese with English abstract)

[16] 謝鴻宇，趙耀龍，楊木壯，等. 基于無棱鏡全站儀的樹冠體積算法[J].中南林業(yè)科技大學(xué)學(xué)報，2014，34(1)：12－17. Xie Hongyu, Zhao Yaolong, Yang Muzhuang, et al. Tree crown volume algorithm based on non-prism total station[J]. Journal of Central South University of Forestry & Technology, 2014, 34(1): 12－17. (in Chinese with English abstract)

[17] 李德仁，李明. 無人機遙感系統(tǒng)的研究進展與應(yīng)用前景[J]. 武漢大學(xué)學(xué)報：信息科學(xué)版，2014，39(5)：505－513，540. Li Deren, Li Ming. Research advance and application prospect of unmanned aerial vehicle remote sensing system[J]. Geomatics and Information Science of Wuhan University, 2014, 39(5): 505－513, 540. (in Chinese with English abstract)

[18] Zarco-Tejada P J, Diaz-Varela R, Angileri V, et al. Tree height quantification using very high resolution imagery acquired from an unmanned aerial vehicle (UAV) and automatic 3D photo-reconstruction methods[J]. European Journal of Agronomy, 2014, 55: 89－99. doi:10.1016/ j.eja.2014.01.004

[19] 劉清旺，李世明，李增元，等. 無人機激光雷達與攝影測量林業(yè)應(yīng)用研究進展[J]. 林業(yè)科學(xué)，2017，53(7)：134－148. Liu Qingwang, Li Shiming, Li Zengyuan, et al. Review on the applications of UAV-based LiDAR and Photogrammetry in forestry[J]. Scientia Silvae Sinicae, 2017, 53(7): 134－148. (in Chinese with English abstract)

[20] Gatziolis D, Lienard J, Vogs A, et al. 3D tree dimensionality assessment using photogrammetry and small unmanned aerial vehicles[J]. PLOS ONE, 2015, 10(9). doi:10.1371/journal. pone.0137765

[21] Gaetano R, Masi G, Poggi G, et al. Marker-controlled watershed-based segmentation of multiresolution remote sensing images[J]. IEEE Transactions on Geoscience＆Remote Sensing, 2015, 53(6): 2987－3004.

[22] 陳崇成，李旭，黃洪宇. 基于無人機影像匹配點云的苗圃單木冠層三維分割[J]. 農(nóng)業(yè)機械學(xué)報，2018，49(2)：149－155，206. Chen Chongcheng, Li Xu, Huang Hongyu. 3D segmentation of individual tree canopy in forest nursery based on drone image-matching point cloud[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(2): 149－155, 206. (in Chinesewith English abstract)

[23] Zahawi R A, Dandois J P, Holl K D, et al. Using lightweight unmanned aerial vehicles to monitor tropical forest recovery[J]. Biological Conservation, 2015, 186: 287－295. doi:10.1016/j.biocon.2015.03.031

[24] 史潔青，馮仲科，劉金成. 基于無人機遙感影像的高精度森林資源調(diào)查系統(tǒng)設(shè)計與試驗[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(11)：82－90. Shi Jieqing, Feng Zhongke, Liu Jincheng. Design and experiment of high precision forest resource investigation system based on UAV remote sensing images[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(11): 82－90. (in Chinese with English abstract)

[25] 劉文萍，仲亭玉，宋以寧. 基于無人機圖像分析的樹木胸徑預(yù)測[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(21)：99－104. Liu Wenping, Zhong Tingyu, Song Yining. Prediction of trees diameter at breast height based on unmanned aerial vehicle image analysis[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(21): 99－104.(in Chinese with English abstract)

[26] Puliti S, Orka H O, Gobakken T, et al. Inventory of small forest areas using an unmanned aerial system[J]. Remote Sensing, 2015, 7(8): 9632－9654.

[27] Kachamba D, Orka H O, Gobakken T, et al. Biomass estimation using 3D data from unmanned aerial vehicle imagery in a tropical woodland[J]. Remote Sensing, 2016, 8(11): 968. doi:10.3390/rs8110968

[28] 何游云，張玉波，李俊清，等. 利用無人機遙感測定岷江冷杉單木樹干生物量[J]. 北京林業(yè)大學(xué)學(xué)報， 2016，38(5)：42－49. He Youyun, Zhang Yubo, Li Junqing, et al. Estimation of stem biomass of individual Abies faxoniana through unmanned aerial vehicle remote sensing[J]. Journal of Beijing Forestry University, 2016, 38(5): 42－49. (in Chinese with English abstract)

[29] 國家林業(yè)局. 國家森林資源連續(xù)清查技術(shù)規(guī)定[S/OL]. 2014: 40－41. https://wenku.baidu.com/view/e0b332884bfe 04a1b0717fd5360cba1aa9118c35.html

[30] 馮仲科，何誠，姚山，等. 一種基于高程等值線法量測樹冠體積的方法：201110164615.4[P]. 2011-11-23.

[31] 同濟大學(xué)數(shù)學(xué)系. 高等數(shù)學(xué)（下冊） [M]. 第六版. 北京：高等教育出版社，2007：132－135.

Tree crown volume measurement method based on oblique aerial images of UAV

Yu Donghai1,2, Feng Zhongke3※

(1.730000,; 2.730000,; 3.100083,)

Tree crown volume is an important basis for monitoring tree growth and estimating tree biomass. Accurate measurement of tree crown volume has always been a difficult issue in forest measurement research. Traditional method of estimating tree crown volume is to bring crown breadth and tree height measurements to inherent empirical models, but it’s faced with problems such as rough operation mode and no precision guarantee of measurement results. In recent years, the emergence of the modern equipments and technologies lay the foundation to achieve high precision tree crown volume measurements. Unmanned aerial vehicle (UAV) oblique aerial photography technology with high-resolution images changes traditional measurement ways, which can use oblique aerial images to generate point cloud data and extract different types of tree geometry parameters by point cloud information. In this paper, a consumer-level multi-rotor UAV named DJI Inspire-1 was used as data acquisition platform, which was equipped with RGB band of ordinary digital camera named Zenmuse-X3. In the Beijing Jiufeng Forest Farm, comprehensively considering the flight operating conditions of the UAV and tree size specifications, we selected eight target trees with different types and sizes. Using the UAV in the manner of spiral flying, we obtained multi-angle oblique aerial images of these target trees. During operation, the main remote controllor controlled the flight status of the UAV and the auxiliary remote controllor controlled the status of the camera haeundae. These two were operated at the same time to collect the oblique aerial images. The following points were the schemes for collecting UAV data: a) In the case of ensuring a safe distance, taking the trunk of the target tree as the center for low-speed flight photography. b) Adjusting the camera pose in real-time during the hovering process so that making sure obtain images of the target tree at different positions and angles. c) Ensuring the overlapping rate of adjacent images collected at the same height exceeded 90%, and the overlapping rate of images collected at different height exceeded 60%. The acquired images were processed through the principle of aerial triangulation for generating three-dimensional point cloud models of target trees. Based on three-dimensional point cloud models, the research segmented the tree crown point cloud by contour lines method and determined the optimal segment number of tree crown point cloud. To extract tree measurement factors, projection method was used to reduce the dimension of the point cloud data. And the measured values of tree height and the arbitrary cross-sectional area of tree crown were calculated by using the key points. According to the established algorithm, the measured volume of the entire tree crown was calculated by accumulating the volume of each rule body after segmentation. Taking actual values by total station as reference, the accuracy of the tree height and tree crown volume measurement results was examined. The results showed that it was a feasible and effective method that the oblique aerial images of UAV were used to establish the three-dimensional point cloud models of single trees and to calculate the tree crown volume. In addition, the average relative error of tree height and tree crown volume of eight target trees was 2.88% and 9.42%, respectively. The accuracy met the standard for tree height and tree crown volume measurement resulted in forestry surveys. In conclusion, three-dimensional point cloud models generated by oblique aerial images of multi-rotor UAV can realize the extraction of measurement factors of single trees, which could be applied to the investigation and protection of ancient and famous trees. This method provides a new approach for the extraction of single trees geometry parameters.

unmanned aerial vehicle; image processing; forestry; oblique photography; point cloud data; parameters extraction of single trees; tree crown volume

2018-07-17

2018-11-19

國家自然科學(xué)基金（U1710123）

于東海，甘肅平?jīng)鋈?，主要研究方向測繪與3S技術(shù)應(yīng)用。 Email：yudh1006@163.com

馮仲科，教授，博士生導(dǎo)師，主要研究方向為精準(zhǔn)林業(yè)、測繪與3S集成研究。Email：fengzhongke@126.com

10.11975/j.issn.1002-6819.2019.01.011

TP79；S771

A

1002-6819(2019)-01-0090-08

于東海，馮仲科. 基于無人機傾斜航空影像的樹冠體積測算方法[J]. 農(nóng)業(yè)工程學(xué)報，2019，35(1)：90－97. doi：10.11975/j.issn.1002-6819.2019.01.011 http://www.tcsae.org

Yu Donghai, Feng Zhongke. Tree crown volume measurement method based on oblique aerial images of UAV[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(1): 90－97. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.01.011 http://www.tcsae.org