劉海洋 馮仲科 呼 諾 劉金成 于新文

(1.北京林業(yè)大學(xué)精準(zhǔn)林業(yè)北京市重點(diǎn)實(shí)驗(yàn)室， 北京 100083； 2.中國林業(yè)科學(xué)研究院資源信息研究所， 北京 100091)

0 引言

樹高測量是森林資源調(diào)查過程中的重要環(huán)節(jié)，也是評價森林生長情況、材積量、碳儲量以及林分收獲出材量的重要數(shù)據(jù)依據(jù)，是我國森林資源連續(xù)清查中非常重要的調(diào)查因子[1-2]。由于林分情況復(fù)雜，現(xiàn)階段我國在森林資源連續(xù)清查中，樹高測量精度不能得到較好保障，引起了樹木樹高計算中增長速度過快或過慢，甚至出現(xiàn)樹高負(fù)增長情況，測量數(shù)據(jù)嚴(yán)重偏離了客觀規(guī)律，給森林調(diào)查、林業(yè)生產(chǎn)和林木生長研究工作帶來了諸多問題。由于長勢傾斜的立木樹高測量與實(shí)際樹長偏差較大，嚴(yán)重影響材積量和出材量的計算和評估，始終是林業(yè)測量中的難題。因此，研制一種操作簡單、測量精準(zhǔn)、適應(yīng)復(fù)雜林分測量和傾斜立木精準(zhǔn)測量的樹高測量裝置，解決林業(yè)樹高測量中的實(shí)際需求，提供可靠樹高數(shù)據(jù)，對發(fā)展精準(zhǔn)林業(yè)具有重要的意義[3-4]。

隨著科技不斷發(fā)展和林業(yè)生產(chǎn)的迫切需求，一大批國內(nèi)外專家學(xué)者提出了很多先進(jìn)的研究理論、技術(shù)和測量方法[5-7]，同時研制了一批新的測量設(shè)備，其新理論和測量方法主要體現(xiàn)在遙感反演[8-11]、激光測距、三維激光掃描[12-14]、攝影測量等領(lǐng)域。國外的樹高測量裝備價格昂貴，使用激光測距傳感器兩次或多次測距的方式進(jìn)行測量，樹頂部的激光點(diǎn)較難確定和測量，操作性差并且測量存在較大誤差，不適用于大面積的林業(yè)資源調(diào)查工作，如日本尼康公司生產(chǎn)的Forestry 550[15]型激光測距裝備等。國外的專家學(xué)者利用這類儀器進(jìn)行了大量的學(xué)術(shù)研究，如RUTTEN等[16]使用瑞典Hagl?f公司生產(chǎn)的Vertex IV Hypsometer有選擇地記錄研究山地森林的結(jié)構(gòu)和組成，測量裝備測量的局限性給森林研究工作帶來了限制。近年來國內(nèi)的專家為了發(fā)展精準(zhǔn)林業(yè)，適應(yīng)林業(yè)生產(chǎn)信息化和智能化的要求，對新型樹木測量設(shè)備與測樹方法進(jìn)行大量研究[17-21]。但目前國內(nèi)研制的儀器還存在模仿國外產(chǎn)品，在樹高測量方面仍使用傳統(tǒng)的測量原理和方法，無法解決操作性差、測量精度低等問題。與此同時眾多國內(nèi)外研究樹高測量儀器時，大多都未解決長勢傾斜立木樹高測量問題，對傾斜立木樹高測量誤差較大，不能測量樹干實(shí)際長度，無法精準(zhǔn)計算和評估傾斜立木的材積和出材量。

針對上述樹高測量存在的問題，本文以國家森林資源調(diào)查為實(shí)際需求，綜合考慮樹高測量過程中測量裝備的精度要求和作業(yè)環(huán)境復(fù)雜等特點(diǎn)，基于傳感器技術(shù)和嵌入式系統(tǒng)，研制一種便攜式精準(zhǔn)立木樹高測量裝置，實(shí)現(xiàn)復(fù)雜林分環(huán)境中立木樹高的精準(zhǔn)測量和傾斜立木樹長測量，為森林計測裝備提供新的測量方法和技術(shù)手段。

1 原理與方法

1.1 硬件設(shè)計

便攜式精準(zhǔn)立木樹高測量裝置由CPU、GPU、存儲設(shè)備RAM/ROM、激光測距儀、顯示屏、高清攝像頭、無線通信模塊、GPS模塊、高精度陀螺儀芯片、電源、天線等組成。其中中央處理器CPU采用Cortex-A9四核處理器，具有高速的運(yùn)算性能，并采用三星S5M8767電源管理模塊和10 000 mA·h鋰電池供電，可降低功耗、增加野外作業(yè)時長。

傳統(tǒng)的攝影測量裝備在手持測量樹高時誤差較大，為保證測量精度通常固定在三腳架上，在復(fù)雜林分的樹高測量中具有較大的局限性。本文便攜式精準(zhǔn)立木樹高測量裝置、激光測距儀和高清攝像頭均采用了前置設(shè)計，將激光測距儀和高清攝像頭并行排列置于裝置上方，測距方向垂直于顯示屏，這種設(shè)計可以在手持測量時，有效避免裝備位置偏移產(chǎn)生的測量誤差。硬件構(gòu)成和試驗(yàn)樣機(jī)如圖1所示。

1.2 軟件設(shè)計

樹高測量裝置軟件采用嵌入式系統(tǒng)和Android 7.1.1系統(tǒng)。樹高測量軟件APP采用Java語言匯編實(shí)現(xiàn)，在開發(fā)環(huán)境Android Studio 2.2下集成后載入設(shè)備，采用流程化和模塊化設(shè)計，包括樹高測量模塊、數(shù)據(jù)管理模塊和數(shù)據(jù)計算模塊，軟件流程圖如圖2所示。

圖1 硬件構(gòu)成和試驗(yàn)樣機(jī)圖Fig.1 Diagrams of hardware design and test prototype

圖2 樹高測量裝置APP軟件設(shè)計流程圖Fig.2 Flow chart of APP software design for tree height measuring device

樹高測量裝置APP主要實(shí)現(xiàn)樹高測量功能，此外還能對傾斜立木樹長進(jìn)行矯正測量。

在樹高測量基礎(chǔ)上，利用SQLite數(shù)據(jù)庫系統(tǒng)，實(shí)現(xiàn)測量數(shù)據(jù)管理和數(shù)據(jù)分析功能，其中數(shù)據(jù)管理包括：數(shù)據(jù)記錄功能記錄被測量樹GPS信息、胸徑、坡度、坡向等信息；數(shù)據(jù)查詢功能可根據(jù)數(shù)據(jù)屬性條件查詢，查詢同類別的數(shù)據(jù)(如地點(diǎn)、樹高取值范圍等)；數(shù)據(jù)編輯功能可將已測數(shù)據(jù)進(jìn)行編輯(包括增加、刪除、修改等)；數(shù)據(jù)導(dǎo)出功能可生成Excel數(shù)據(jù)文件并導(dǎo)出。數(shù)據(jù)分析包括：通過調(diào)用已測量數(shù)據(jù)，進(jìn)行樣地樹高統(tǒng)計分析(如最大值、最小值、平均值、標(biāo)準(zhǔn)方差等)；根據(jù)樹種、樹高、胸徑等信息實(shí)現(xiàn)單木材積計算分析；根據(jù)單木材積和樹種信息實(shí)現(xiàn)林木出材量計算分析；根據(jù)單木材積實(shí)現(xiàn)林分蓄積量計算分析。

1.3 樹高測量

1.3.1樹高測量原理

本文樹高測量主要利用攝影測量和三角函數(shù)原理進(jìn)行計算，通過改進(jìn)測量算法簡化作業(yè)操作，僅需2個操作步驟即可完成測量，測量方法如圖3所示。

圖3 樹高測量示意圖Fig.3 Diagram of tree height measurement

首先使用高清攝像頭拍攝樹根位置，獲取與激光測距傳感器、高精度陀螺儀芯片間的距離d和設(shè)備仰角α，然后利用高清攝像頭拍攝樹頂位置，獲取拍攝角β，計算得到立木樹高H，計算式為

(1)

1.3.2傾斜立木樹高矯正

使用設(shè)備貼緊樹干測量樹干傾斜角度，在樹木干型彎曲時利用圖像處理技術(shù)對傾斜立木樹高測量進(jìn)行矯正，利用邊緣檢測算法使用Canny算子進(jìn)行樹干邊緣檢測，首先進(jìn)行灰度化處理，計算公式為

f=0.299R+0.587G+0.114B

(2)

式中f——灰度

R、G、B——紅、綠、藍(lán)顏色分量

對圖像進(jìn)行高斯濾波，公式為

(3)

式中x、y——點(diǎn)坐標(biāo)值σ——標(biāo)準(zhǔn)差

用一階有限差分近似求取灰度梯度，計算圖像梯度幅值和方向，對梯度進(jìn)行非極大值抑制，在獲得梯度方向和大小后，對整幅圖像掃描，去除非邊界點(diǎn)，對每個像素進(jìn)行檢查，判斷該點(diǎn)的梯度是否為周圍具有梯度方向的點(diǎn)中最大。選取2個閾值，根據(jù)較高閾值得到一個邊緣圖像，這一圖像含有少量假邊緣，但由于閾值較高，產(chǎn)生的圖像邊緣可能不閉合， 為解決此問題采用另一個低閾值。在高閾值圖像中把邊緣連接成輪廓，當(dāng)達(dá)到輪廓的斷點(diǎn)時，在斷點(diǎn)的8鄰域點(diǎn)中尋找滿足低閾值的點(diǎn)，再根據(jù)此點(diǎn)收集新的邊緣，直到整個圖像邊緣閉合。樹干邊緣檢測如圖4所示。

圖4 邊緣檢測處理結(jié)果Fig.4 Edge detection processing result

利用輪廓近似法提取樹干輪廓邊緣點(diǎn)，獲得樹干邊緣離散點(diǎn)坐標(biāo)信息，使用最小二乘法對樹干邊緣信息進(jìn)行曲線擬合，處理后再進(jìn)行線性擬合，求得直線斜率后轉(zhuǎn)換為樹干傾斜角γ，對樹干高度進(jìn)行矯正計算。根據(jù)公式求得傾斜立木樹長H′和樹高H，傾斜立木樹高測量示意圖如圖5所示。

圖5 傾斜立木樹高測量示意圖Fig.5 Sketch of slant standing tree height measurement

在傾斜立木傾斜正面使用測距儀測量設(shè)備到樹根的距離d和設(shè)備仰角α，求出虛擬高度h

(4)

最后測量立木傾斜角γ，求出立木樹長

(5)

將式(4)代入式(5)，化簡得出

(6)

由樹長H′可求出傾斜立木樹高

H=H′sinγ

(7)

傾斜立木樹長H′可用于立木材積精準(zhǔn)計算，減少材積計算誤差，提高材積測量精度。

2 試驗(yàn)

2.1 樹高測量

試驗(yàn)在西山試驗(yàn)林場、松山國家自然保護(hù)區(qū)等地進(jìn)行，選取不同樹高的油松、雪松、側(cè)柏等針葉樹，蒙古櫟、白蠟、欒樹等闊葉樹進(jìn)行立木樹高測量。選用南方測繪NTS-282R6型全站儀進(jìn)行對比試驗(yàn)，并以全站儀測量數(shù)據(jù)為樹高測量真值進(jìn)行誤差計算，測量結(jié)果如圖6所示。

圖6 測量結(jié)果Fig.6 Measurement results

試驗(yàn)分別測量56棵樹，樹高在5.6～15.4 m之間，從測量數(shù)據(jù)分析得出其測量相對誤差范圍為-0.4～1.1 m，平均絕對誤差為0.19 m。誤差分布如圖7所示。

圖7 誤差分布Fig.7 Error distribution map

根據(jù)測量誤差分別計算乖離率為0.05(Bias)、0.58(Bias%)、均方根誤差為0.40、4.30 m[22]。

2.2 不同樹種分析

將測量數(shù)據(jù)根據(jù)不同樹種劃分，并分別計算測量誤差，結(jié)果如表1所示。

表1 不同樹種樹高測量結(jié)果Tab.1 Measurement results of different tree species

由表1可以看出，本測量裝置所測量不同樹種樹高平均絕對誤差范圍在0.12～0.28 m之間，測量不同樹種樹高相對誤差范圍為0.97%～2.79%，不同樹種平均測量精度達(dá)98.04%，不同樹種的測量誤差均能滿足林業(yè)中樹高測量要求，試驗(yàn)裝置適用于多樹種立木樹高測量。

2.3 傾斜立木樹長測量試驗(yàn)

在松山國家自然保護(hù)區(qū)對不同傾斜程度的立木分別進(jìn)行測量，得出立木樹長，并使用全站儀分別測量樹頂和樹根空間位置坐標(biāo)，計算得出樹高真值，試驗(yàn)結(jié)果如表2所示。

測量設(shè)備在不同樹種間測量平均絕對誤差在0.1～0.5 m之間，各種不同傾斜立木間測量結(jié)果的相對誤差為1.09%～4.94%,平均測量誤差為3.11%，測量矯正效果較好，能有效測量樹木樹長，提高傾斜立木材積計算精度。

3 結(jié)論

(1)利用電子信息技術(shù)、傳感器技術(shù)和圖像處理技術(shù)，在樹高測量研究上提出新的測量方法并改進(jìn)了樹高測量技術(shù)，研制出一種便攜式精準(zhǔn)立木樹高測量儀，實(shí)現(xiàn)了復(fù)雜林分環(huán)境下的精準(zhǔn)立木樹高測量和傾斜立木樹長測量，同時此設(shè)備還實(shí)現(xiàn)了數(shù)據(jù)記錄和存儲功能，在記錄樣地信息的同時可根據(jù)輸入樹種、胸徑等其他信息，計算該樣地中的材積、出材量、蓄積量等信息?？筛鶕?jù)需求，與手機(jī)數(shù)據(jù)接收端協(xié)同使用，實(shí)現(xiàn)對測量數(shù)據(jù)無線實(shí)時接收和讀取，并實(shí)現(xiàn)測量數(shù)據(jù)的調(diào)用、編輯、導(dǎo)出Excel等功能。

表2 傾斜立木樹長測量結(jié)果Tab.2 Measurement results of tilted tree

(2)試驗(yàn)表明本文設(shè)計的精準(zhǔn)立木樹高測量設(shè)備在樹高測量中測量精度達(dá)到98.04%，傾斜立木樹長測量平均誤差為3.11%，符合國家森林資源連續(xù)清查等森林調(diào)查的樹高測量精度要求。經(jīng)測量作業(yè)實(shí)踐證明，利用攝影方式解決了傳統(tǒng)激光測高儀器在復(fù)雜林分環(huán)境下樹冠激光點(diǎn)難確定的問題，優(yōu)化了森林計測儀器的測量方法，利用圖像處理技術(shù)計算傾斜立木樹長，精準(zhǔn)估算傾斜立木材積和出材量，可為林業(yè)調(diào)查工作人員提供可靠數(shù)據(jù)支持，同時為樹高測量設(shè)備的研制提供新的研究方向。