邢 浩， 李 昀， 潘遠(yuǎn)明， 陳名騰

(1.北京林業(yè)大學(xué) 信息學(xué)院,北京 100083;

2.中科方德軟件有限公司 基礎(chǔ)軟件國家工程研究中心, 北京 100190)

0 引 言

角規(guī)測樹是森林調(diào)查常用的一種地面調(diào)查方法，利用角規(guī)可以快速測得林分單位面積胸高斷面積，進(jìn)而結(jié)合林分高度就可以快速獲得林分單位面積蓄積量[1]。盡管近年來遙感、攝影測量等技術(shù)被廣泛應(yīng)用于大面積森林調(diào)查，但因其技術(shù)特點并不適用于林分級別的小面積數(shù)據(jù)獲取[2]，在這方面，地面調(diào)查和測量是不可替代的。

用于地面森林調(diào)查的便攜式測量設(shè)備主要有測高器、測徑儀、測距儀、角規(guī)等。根據(jù)文獻(xiàn)研究，可將其發(fā)展過程劃分為機(jī)械式測量設(shè)備、電子測量設(shè)備和智能測量設(shè)備 3個階段。機(jī)械式測量設(shè)備以1928年麻生測高器的發(fā)明為標(biāo)志，在此后的半個多世紀(jì)內(nèi)，出現(xiàn)了包括克里斯屯(Christen)測高器[3]、畢特利希桿式角規(guī)、光學(xué)測樹儀[4，5]等典型設(shè)備及其各類改型[6]。20世紀(jì)80年代開始，隨著遙測技術(shù)的發(fā)展以及傳感器、單片機(jī)等電子技術(shù)的成熟，對電子測量設(shè)備的研究占據(jù)主導(dǎo)，包括利用激光、超聲波對直徑[7]、高度[8]、距離[9]等參數(shù)進(jìn)行測量的研究，應(yīng)用數(shù)字?jǐn)z影技術(shù)的嘗試[10]和對電子角規(guī)[11]的研究，并出現(xiàn)了瑞典Vertex系列測高測距儀[9]和德國Leica、日本Nikon等公司生產(chǎn)的優(yōu)秀產(chǎn)品。21世紀(jì)以來，隨著信息化技術(shù)的加速發(fā)展，PDA、智能手機(jī)等智能便攜設(shè)備的出現(xiàn)，地面便攜式測量設(shè)備也呈現(xiàn)出與林業(yè)信息化相匹配的特點,如實現(xiàn)測量數(shù)據(jù)的數(shù)字化采集[12]和實時傳輸[13]；探索測算設(shè)備與PDA等智能便攜設(shè)備的協(xié)同工作模式[14]等。

目前，便攜式測量設(shè)備呈現(xiàn)出電子化、數(shù)字化、智能化的發(fā)展趨勢，這不但有利于提高測量效率[12]，更重要的是測量數(shù)據(jù)的數(shù)字化杜絕了人工輸入錯誤[16]，符合林業(yè)信息化發(fā)展需求。本文設(shè)計了集角規(guī)、測高器于一體的便攜電子測樹儀GForest，以適應(yīng)林業(yè)信息化發(fā)展的需要。

1 設(shè)計方案

1.1 計測思路

角規(guī)抽樣基于同心圓原理、三角函數(shù)原理和擴(kuò)大園原理[1]，在進(jìn)行角規(guī)抽樣調(diào)查時，根據(jù)林分平均直徑等特征選擇不同的斷面積系數(shù)Fg，通常為0.5，1，2，4。選定斷面積系數(shù)后，根據(jù)林分面積確定抽樣點的個數(shù)N，每個抽樣點使用角規(guī)產(chǎn)生的固定視角進(jìn)平繞測一周，樹干胸徑與視角兩條邊相割計數(shù)為1，相切計數(shù)為0.5，相離計數(shù)為0，相加得到該抽樣點的計數(shù)值Zi(i=1,2,3,…，N)。角規(guī)抽樣結(jié)束后，結(jié)合立木高度測量得到的標(biāo)準(zhǔn)木高度H，可以計算林分蓄積量M,如公式(1)所示

(1)

針對森林調(diào)查中常用角規(guī)調(diào)查來測算林分蓄積量這一實際，本文考慮集成數(shù)字化測高和角規(guī)計測，設(shè)計一種便攜式產(chǎn)品提供角規(guī)繞測計數(shù)、林分蓄積量計算、計測數(shù)據(jù)存儲等功能。在實際設(shè)計中，通過硬軟件的協(xié)同來實現(xiàn)設(shè)計目標(biāo)。其中，測高器部分依據(jù)主流測高器所采用的三角幾何測高原理，采用相應(yīng)的傳感器采集數(shù)據(jù)，軟件上進(jìn)行處理計算得到立木高度值并自動存儲；角規(guī)硬件可以采用現(xiàn)在外掛常用的桿式角規(guī)的簡易產(chǎn)品——角規(guī)片，計數(shù)值和抽樣點個數(shù)通過軟件計算和控制；蓄積量計算可由軟件依據(jù)角規(guī)計數(shù)結(jié)果和樹高測定結(jié)果完成。GForest內(nèi)置3種測高方法，包括固定距離測高、輔助標(biāo)尺固定點測高、輔助標(biāo)尺不定點測高[8]，實際使用過程中用戶可以根據(jù)地形等實際情況選擇合適的方法測量立木高度，其中，固定距離測高法需提前量出測量者立足點與被測立木的距離，輔助標(biāo)尺固定點測高法和輔助標(biāo)尺不定點測高法則需要配合輔助標(biāo)尺實現(xiàn)。角規(guī)參數(shù)、測高方法等參數(shù)的設(shè)置可以通過軟件實現(xiàn)并保存設(shè)置方案，在設(shè)備開機(jī)時自動加載。

1.2 硬件選型與設(shè)計

為實現(xiàn)對立木高度的精確測量和對測量數(shù)據(jù)的計算處理，需要傾角傳感器和嵌入式處理器配合顯示屏、操作按鍵、蜂鳴器、激光發(fā)射器等硬件實現(xiàn);另外需存儲卡、USB接口等實現(xiàn)數(shù)據(jù)存儲和數(shù)據(jù)傳輸。為適應(yīng)林業(yè)野外調(diào)查的實際需要，硬件選型時，重點在保證測量精度的基礎(chǔ)上重點考慮低功耗、低成本、體積小等因素。經(jīng)過綜合比較，傾角傳感器選用芬蘭VTI公司生產(chǎn)的單軸傳感器SCA830—D07，工作電壓為3.0～3.6 V，測量范圍為±1gn，精度偏移為±30 mgn，具備高分辨率16位A/D轉(zhuǎn)換和SPI接口。MCU選型鎖定在32位STM32F103系列產(chǎn)品上，該系列產(chǎn)品基于ARM Cortex—M3內(nèi)核，具備高達(dá)72 MHz的主頻，32～512 k的Flash和最高達(dá)64k的SRAM，性價比很高。從性能、接口等方面的實際需求出發(fā)，本文最終選擇STM32F103RBT6作為微控單元。除核心硬件外，本文采用一塊以ST7565驅(qū)動芯片為核心的128×64 LCD顯示屏提供操作過程中所需的交互控制提示和結(jié)果顯示；采用Micro-SD卡作為數(shù)據(jù)存儲介質(zhì)；采用可充電鋰電池給設(shè)備供電并結(jié)合CN3052A芯片進(jìn)行充電管理；使用Mini-USB接口傳輸數(shù)據(jù)和給鋰電池充電。

設(shè)備硬件設(shè)計方案如圖1所示，整個硬件電路采用3.3 V直流供電。其中，SCA830—D07傾角傳感器將測量過程中的傾角值通過SPI接口傳輸給MCU處理；對Micro-SD卡進(jìn)行數(shù)據(jù)讀寫操作，以SPI模式進(jìn)行；因MCU只配置2個SPI接口，顯示屏同MCU之間采用GPIO模擬的SPI串行通信進(jìn)行控制顯示；GForest與PC之間的串口通信利用MAX3232芯片實現(xiàn)，結(jié)合JTAG接口進(jìn)行程序下載和設(shè)備調(diào)試。除上述主要模塊硬件和電路以外，整體硬件結(jié)構(gòu)還包括時鐘電路、復(fù)位電路等外圍電路。在實際應(yīng)用過程中，MCU通過處理SCA830—D07的傾角值計算得出樹木高度；進(jìn)行角規(guī)檢尺作業(yè)時，通過按鍵I/O信號進(jìn)行林分胸高斷面計數(shù)；根據(jù)標(biāo)準(zhǔn)木樹高和胸高斷面計數(shù)即可計算得到林分蓄積量。操作過程和測算結(jié)果通過顯示屏進(jìn)行顯示，相應(yīng)的測量數(shù)據(jù)會存儲在Micro-SD卡上。

圖1 GForest硬件結(jié)構(gòu)框圖

1.3 軟件設(shè)計

在森林地面調(diào)查中，測高、角規(guī)抽樣和林分蓄積量計算存在操作流程上的先后順序和計測數(shù)據(jù)傳遞。為實現(xiàn)多功能集成，GForest軟件根據(jù)調(diào)查實際進(jìn)行抽象，將各項功能設(shè)計為獨立模塊，模塊間能夠進(jìn)行相應(yīng)的數(shù)據(jù)傳遞，以單任務(wù)調(diào)度模式實現(xiàn)對各功能模塊的調(diào)用。為減少高度測量過程中因手持設(shè)備抖動造成的測量誤差，提高測量精度，采樣控制算法設(shè)計為在每次測量傾角時在系統(tǒng)內(nèi)部進(jìn)行多次采樣，去除異常傾角數(shù)據(jù)后，在多個傾角樣本數(shù)據(jù)基礎(chǔ)上，計算傾角平均值，并以此來計算立木高度。

在選定硬件(HW)的基礎(chǔ)上，GForest操作程序可以分為硬件抽象層(HAL)和應(yīng)用程序?qū)?APP)兩層。硬件抽象層主要包括GPIO總線、RTC時鐘以及顯示屏(LCD)、接口(USB)、電源(PM)、傳感器(SEN)等主要設(shè)備的驅(qū)動程序，并提供面向應(yīng)用程序?qū)?APP)的接口。應(yīng)用程序?qū)又饕ㄕ{(diào)度菜單(MENU)、系統(tǒng)設(shè)置(SETTING)、高度測量(ALTIMETER)和角規(guī)抽樣(GAUGE)4個功能模塊和電源管理等獨立模塊。

GForest啟動后，MENU通過輪循方式監(jiān)聽按鍵中斷調(diào)用其他功能模塊，配合界面顯示實現(xiàn)用戶在設(shè)備使用過程中的各種操作。SETTING模塊包括樹種選擇和參數(shù)設(shè)定兩部分，在使用過程中，用戶需根據(jù)林分實際情況選擇對應(yīng)樹種和形數(shù)，用于林分蓄積量的計算；參數(shù)設(shè)定包括對測高方法、標(biāo)尺長度、角規(guī)常數(shù)、觀測者眼高等功能參數(shù)的選擇設(shè)定和對日期、自動關(guān)機(jī)時間、音效、語言等系統(tǒng)參數(shù)的選擇設(shè)定。ALTIMETER模塊根據(jù)用戶選定的測高方法和功能參數(shù)，通過傾角傳感器捕獲的傾角值計算立木高度。GAUGE模塊記錄每個抽樣點的數(shù)值，抽樣完畢后根據(jù)設(shè)定的樹種形數(shù)和測得的標(biāo)準(zhǔn)木高度計算林分蓄積量。各模塊在工作過程中的相互調(diào)用關(guān)系如圖2所示。

圖2 主要模塊調(diào)用關(guān)系

2 樣機(jī)實測

本文以赤峰市喀喇沁旗旺業(yè)甸林場主伐作業(yè)區(qū)內(nèi)的樟子松為測量對象對GForest樣機(jī)進(jìn)行實地測試，在伐區(qū)內(nèi)選取49株樟子松進(jìn)行實測分析。設(shè)備運行情況見圖3，測試過程中參數(shù)設(shè)置、高度測量、角規(guī)計數(shù)、蓄積量計算、數(shù)據(jù)查看等各項功能運行良好，相應(yīng)操作記錄均存儲在Micro-SD卡上。

圖3 GForest運行截圖

立木高度測量方面，采用瑞典Hagl?f生產(chǎn)的VertexIII測高儀(以下簡稱VertexIII)作為對比。立木伐倒前分別采用GForest和VertexIII測量立木高度得到2組測量值，伐倒后測量倒木得到高度真值。經(jīng)測量，49株樟子松實際立木高度分布在10.5～17.0 m之間，2種設(shè)備測量相對誤差分布如圖4所示。

圖4 GForest和VertexIII測量相對誤差散點分布圖

經(jīng)統(tǒng)計，GForest最大測量誤差為8 %，VertexIII最大測量誤差為10 %。在GForest與VertexIII各自測得的49組數(shù)據(jù)中，GForest有13組測量相對誤差超出5 %，超過73 %的數(shù)據(jù)落在[-5 %,5 %]區(qū)間內(nèi)；VertexIII有16組測量相對誤差超出5 %，67 %的數(shù)據(jù)落在[-5 %,5 %]區(qū)間內(nèi)。以上分布說明同等條件下，GForest獲得高精度數(shù)據(jù)的概率要比VertexIII稍高。此外，兩者的平均測量相對誤差均為4 %，能夠滿足林業(yè)外業(yè)測量需要。

通過對兩者測量誤差進(jìn)行頻率分布和和正態(tài)分布統(tǒng)計，得到兩者測量誤差的直方圖和正態(tài)分布圖(圖5)。從直方圖中可以看出：GForest產(chǎn)生的測量誤差多集中在[-7 %,7 %]，且在該區(qū)間內(nèi)基本呈現(xiàn)對稱分布。另一邊，VertexIII手持式測高器的測量誤差則呈比較明顯的偏態(tài)分布，誤差多集中在[0 %,10 %]的區(qū)域，但在高精度區(qū)間[-1 %,1 %]內(nèi)，VertexIII的數(shù)據(jù)量要明顯高于前者。通過正態(tài)分布擬合可以進(jìn)一步說明上述結(jié)論，盡管VertexIII的誤差分布的集中性(標(biāo)準(zhǔn)差3.7 %)較GForest(標(biāo)準(zhǔn)差4.1 %)要略好，但兩者的中位數(shù)(VertexIII為3 %，GForest為-1 %)說明GForest的誤差分布要更加符合隨機(jī)誤差分布規(guī)律。總體而言，兩者精度基本持平。

圖5 GForest和VertexIII測量相對誤差頻率分布直方圖和正態(tài)分布圖

3 結(jié) 論

本文根據(jù)嵌入式軟硬件開發(fā)方法，設(shè)計實現(xiàn)了一種林業(yè)手持式多功能測樹儀GForest，具備立木高度測量、角規(guī)抽樣計數(shù)以及數(shù)據(jù)處理能力，能夠計算林分蓄積量等林分指標(biāo)，具有數(shù)字化、多功能、成本低、適用性強(qiáng)等優(yōu)點。GForest結(jié)合配套的數(shù)據(jù)處理軟件，在野外使用過程中能夠脫離紙筆記錄的傳統(tǒng)方式，大幅提高測量效率。

參考文獻(xiàn):

[1]孟憲宇.測樹學(xué)[M].北京:中國林業(yè)出版社,2006:230-236.

[2]Petri R,Juha H,Hannu H,et al.Calibration of laser-derived tree height estimates by means of photogrammetric techniques[J].Journal of Forest Research,2004,19(6):524-528.

[3]李建華.基于三角原理的森林測高器研制與應(yīng)用[D].泰安：山東農(nóng)業(yè)大學(xué)，2011:19-25.

[4]唐義全. GCG—1型光學(xué)測高器的結(jié)構(gòu)原理與使用方法[J].林業(yè)資源管理，1987(2):86-91.

[5]張紹良,杜麗雁.DQL—1型測樹羅盤儀的設(shè)計與應(yīng)用[J].林業(yè)資源管理， 1982(4):48-50.

[6]陳愛國.多參數(shù)勃魯萊氏測高器刻度盤：中國，2347138[P].[1999—11—03].2013—01—20.

[7]Carr B.Using laser technology for forestry and engineering applications[J].Compiler,1992,10(4):5-16.

[8]鄢前飛.林業(yè)數(shù)字式測高測距儀的研制[J].中南林業(yè)科技大學(xué)學(xué)報，2007,27(5):67-70.

[9]Hagl?f Sweden AB .VertexIV[EB/OL].[2013—02—05].http:∥www.haglofcg.com/index.php.

[10] Clark N,Wynne R H,Schmoldt D L,et al.An assessment of the utility of a non-metric digital camera for measuring standing tree-s[J].Computers and Electronics in Agriculture,2000,28:151-169.

[11] 馮仲科,徐禎祥,杰林德·羅斯納爾.電子角規(guī)測樹儀及自動測樹方法：中國，1570557[P].2005—01—26.[2013—01—20].

[12] Inman-Narahari F,Giardina C,Ostertag R,et al.Digital data co-llection in forest dynamics plots[J].Methods in Ecology and Evolution,2010(1):274-279.

[13] J?nes J.Presenting an idea:Log-by-log real-time data transmission from harvester to server over GPRS[J].Forestry Studies,2006,45:177-181.

[14] Boston K.Precision log-making to maximise value recovery from plantation forests[C]∥The First International Precision Forestry Symposium,Seattle,2001:165-169.

[15] Guadagno L,Vandeweerd C,Stevens D,et al.Using PDAs for data collection[J].Applied Nursing Research,2004,17(4):283-291.