蔣興明

摘要：根據(jù)森林調(diào)查標(biāo)準(zhǔn)地調(diào)查法所獲得的胸徑一樹高數(shù)據(jù)，用Excel擬合動(dòng)態(tài)樹高曲線，然后利用Excel的函數(shù)提取擬合樹高曲線方程中的參數(shù)來(lái)自動(dòng)計(jì)算林分平均高和徑階平均高，且計(jì)算結(jié)果會(huì)隨所選擇曲線方程類型的不同以及胸徑一樹高數(shù)據(jù)的增減變動(dòng)均能自動(dòng)更新計(jì)算結(jié)果，避免了因數(shù)據(jù)每次發(fā)生變化后都需要重復(fù)操作，省時(shí)高效。

關(guān)鍵詞：Excel;樹高曲線方程;林分平均高;徑階平均高

中圖分類號(hào)：TP311 文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1009-3044（2020）13-0247-03

1背景

林分中林木高度與胸徑之間存在著正相關(guān)關(guān)系，呈現(xiàn)出一定的變化規(guī)律，為了反映林分樹高與胸徑變化規(guī)律，以縱坐標(biāo)為樹高，橫坐標(biāo)為胸徑，即可建立胸徑一樹高曲線，即樹高曲線，二元材積表法計(jì)算林分蓄積量以及二元材積表導(dǎo)算一元材積表時(shí)都需要樹高曲線。圖解法繪制樹高曲線的傳統(tǒng)方法是用坐標(biāo)紙手工繪制樹高曲線，在樹高曲線上查取林分平均所對(duì)應(yīng)的樹高即是林分平均高，同樣在樹高曲線上可查得各徑階對(duì)應(yīng)的徑階平均樹高。求林分平均樹高還可使用數(shù)式法，即利用胸徑和樹高相關(guān)模型。此外可通過(guò)專業(yè)調(diào)查軟件完成，但其不夠普及和被廣泛應(yīng)用，基于Excel的大眾化和強(qiáng)大的數(shù)據(jù)處理功能，可很方便使用Excel擬合樹高曲線。本文以表1胸徑一樹高數(shù)據(jù)為例，介紹Excel（2016版）擬合樹高曲線方程參數(shù)的動(dòng)態(tài)提取及林分平均高及徑階平均高的自動(dòng)計(jì)算，解決了每次數(shù)據(jù)發(fā)生變動(dòng)都需要重復(fù)操作的難題。

2擬合樹高曲線

2.1錄入數(shù)據(jù)

新建“林分平均高自動(dòng)計(jì)算”的Excel工作薄，將sheetl工作表命名為“樹高曲線”，將胸徑一樹高調(diào)查整理表（見表1）中的數(shù)據(jù)錄入A、B、C三列（見圖1）。

2.2繪制散點(diǎn)圖

依據(jù)B、c列數(shù)據(jù)，利用菜單欄的“插人”圖表功能，建立胸徑一樹高的散點(diǎn)圖。

2.3添加趨勢(shì)線

為散點(diǎn)圖添加趨勢(shì)線，趨勢(shì)線包括指數(shù)、線性、對(duì)數(shù)、多項(xiàng)式、乘冪、移動(dòng)平均等6種，勾選“顯示公式”和“顯示R平方值”，則散點(diǎn)圖中并會(huì)添加相應(yīng)的趨勢(shì)線、回歸方程和擬合程度指標(biāo)R2。

2.4選擇樹高曲線方程

R2取值范圍在0～l之間，當(dāng)趨勢(shì)線的R²等于l或接近1時(shí)，其可靠性最高。繪制樹高曲線就是以平均胸徑為自變量、平均樹高為因變量所建立的回歸方程，所以選擇6種回歸方程中Rz最大者作為樹高曲線。并將“趨勢(shì)預(yù)測(cè)”向前、向后各推延1個(gè)周期。

2.5設(shè)置樹高曲線圖表格式

2.5.1設(shè)置水平軸的格式

對(duì)“坐標(biāo)軸選項(xiàng)”進(jìn)行設(shè)置，水平軸“邊界”的“最小值”與每木調(diào)查的起測(cè)徑階值相同，“最大值”不做設(shè)置，“單位”為徑階距l(xiāng)cm、2cm或4cm，本例設(shè)置為2，水平軸刻度值“數(shù)字”設(shè)為整數(shù)，即不保留小數(shù)，取消千分位分隔符。垂直（值）軸不做修改。

2.5.2設(shè)置曲線方程的格式

對(duì)圖表中的曲線方程進(jìn)行設(shè)置，數(shù)字類別選擇“數(shù)字”格式，小數(shù)位數(shù)為8位或與材積公式參數(shù)的小數(shù)位數(shù)一致，取消千分位分隔符。最后再適當(dāng)調(diào)整曲線方程的字體大小。

此外將圖表的標(biāo)題修改為“樹高曲線”，添加次要網(wǎng)格線，最后再根據(jù)自己的需要設(shè)置趨勢(shì)線的顏色、線寬線型，調(diào)整圖表的大小和圖表的位置。

2.6制作動(dòng)態(tài)樹高曲線

森林調(diào)查時(shí)，不同的林分，其徑階分布、株數(shù)多少是不固定的，為滿足不同林分都可使用，需要將樹高曲線制作成動(dòng)態(tài)的，讓樹高曲線可隨數(shù)據(jù)的增減而動(dòng)態(tài)變動(dòng)。

2.6.1定義名稱

2.6.2動(dòng)態(tài)樹高曲線

在散點(diǎn)圖表上單擊鼠標(biāo)右鍵，“選擇數(shù)據(jù)”快捷方式，進(jìn)入數(shù)據(jù)源對(duì)話框，點(diǎn)選“編輯”打開編輯數(shù)據(jù)系列對(duì)話框，重新編輯修改x軸系列值為“=林分平均高自動(dòng)計(jì)算.xlsxr！胸徑”;Y軸系列值為“=林分平均高自動(dòng)計(jì)算.xlsx'！樹高”。這樣，樹高曲線即可隨著B列和c列的數(shù)據(jù)增減而變動(dòng)，實(shí)現(xiàn)樹高曲線的動(dòng)態(tài)化。

3計(jì)算平均高

此處的平均高僅指林分條件平均高和林分徑階平均高。林分條件平均高是根據(jù)林分平均胸徑在樹高曲線上查得的樹高。林分徑階平均高是根據(jù)各個(gè)徑階中值在樹高曲線上所對(duì)應(yīng)的樹高，它是二元材積表法計(jì)算林分蓄積量的主要因子。傳統(tǒng)方法是在手繪樹高曲線上查得，費(fèi)時(shí)費(fèi)勁效率低，精度低還易出錯(cuò)。也可在Excel中利用擬合的曲線方程，代人徑階中值手工輸入計(jì)算公式計(jì)算徑階平均高，但是，胸徑一樹高數(shù)據(jù)發(fā)生變動(dòng)則擬合曲線方程的參數(shù)值也會(huì)發(fā)生改變，這時(shí)，就得重新輸入新的計(jì)算公式，不夠人性化。若能自動(dòng)提取曲線方程的參數(shù)，則可自動(dòng)利用曲線方程完成林分條件平均高和各徑階平均高的計(jì)算。

3.1提取曲線方程R2值和參數(shù)

3.1.1提取R2值

注：對(duì)提取的值用ROUNDO函數(shù)四舍五入保留8位小數(shù);公式中的符號(hào)均為英文輸入法狀態(tài)下的符號(hào)（下同），否則會(huì)出錯(cuò)。公式中的“樹高”是c列樹高的定義名稱，“胸徑”是B列胸徑的定義名稱，R²值提取結(jié)果見圖2-②。圖1中，E19-J19也是利用上述提取公式所得。

3.1.2提取參數(shù)

二次多項(xiàng)式（簡(jiǎn)稱二次式，下同）有A、B、C三個(gè)參數(shù)，三次多項(xiàng)式（簡(jiǎn)稱三次式，下同）有A、B、C、D四個(gè)參數(shù)，其他曲線方程只有A和B兩個(gè)參數(shù)。各參數(shù)提取結(jié)果見圖2一③，提取曲線方程參數(shù)的公式如下（保留8位小數(shù)）：

3.1.3列示曲線公式

用字符連接符（&）將所提取的參數(shù)按曲線方程表達(dá)式的形式連接成公式形式。公式連接后的結(jié)果見圖2一④，各個(gè)曲線

至此，完成了平均高的自動(dòng)計(jì)算。使用時(shí)，先錄入胸徑一樹高調(diào)查結(jié)果，之后在圖表中雙擊趨勢(shì)線打開“趨勢(shì)線選項(xiàng)”，選擇恰當(dāng)?shù)那€作為擬合的樹高曲線，然后，在D17單元格的下拉列表中選擇相應(yīng)的曲線類型即可計(jì)算林分條件平均高和徑階平均高。

4結(jié)論與討論

通過(guò)定義“樹高”和“胸徑”名稱，實(shí)現(xiàn)擬合的樹高曲線可隨數(shù)據(jù)的增減變化而實(shí)時(shí)變動(dòng)，利用Excel函數(shù)提取出擬合方程的參數(shù)，可得到了樹高曲線方程，這樣就可利用樹高曲線方程實(shí)現(xiàn)林分的平均高和徑階平均高計(jì)算的自動(dòng)化，解決了因不同林分調(diào)查的數(shù)據(jù)不同，重復(fù)擬合樹高曲線和計(jì)算的麻煩，提高了計(jì)算效率，在生產(chǎn)實(shí)踐中具有現(xiàn)實(shí)意義。