袁志超，張家盛，李士美，王永濤，崔孝平

(1.青島農(nóng)業(yè)大學 園林與林學院，山東 青島 266109；2.青島市林業(yè)局，山東 青島 266061；3. 嶗山區(qū)林業(yè)局，山東 青島 266061)

立木結(jié)構(gòu)特征決定著樹干材積的數(shù)量和質(zhì)量，是森林計測用表編制的基礎，對森林資源調(diào)查、評估及經(jīng)營管理都有著非常重要的意義。[1-3]森林計測常用的材積表主要為基于立木胸徑、樹高的一元材積表或二元材積表。然而，在林業(yè)資源管理部門處理盜伐案件以及處理林權(quán)糾紛中，經(jīng)常要對被伐木材積做出估計。對于被伐木，能測定的基本特征因子只有地徑，因而僅能通過對地徑的測定間接推算胸徑和樹高，然后用立木材積表計算被伐木的材積。[4,5]因而，研究地徑與胸徑、樹高的關(guān)系，構(gòu)建合適的數(shù)學模型，在林業(yè)資源管理上有重要的現(xiàn)實意義。

近年來，無人機逐漸成為森林資源調(diào)查與監(jiān)測的新手段，已應用于單木樹高、胸徑、冠幅等信息的獲取以及森林信息量化等方面的研究。[6,7]無人機影像在森林結(jié)構(gòu)參數(shù)提取方面主要包括單木樹高、單木冠幅、株數(shù)密度、郁閉度、林分平均高、生物量及蓄積量等因子，然而樹木胸徑的獲取主要是通過影像提取的樹木高度和冠幅來進行反演。[8,9]由于胸徑是立木結(jié)構(gòu)特征的基本參數(shù)，確立胸徑與樹高、冠幅等結(jié)構(gòu)參數(shù)的確切關(guān)系，是運用無人機對森林資源準確計測的前提和基礎。

黑松(Pinusthunbergii)系日本第一次侵占青島后引入的針葉樹種，目前已經(jīng)成為山東半島的重要基調(diào)樹種。然而，有關(guān)黑松立木結(jié)構(gòu)特征參數(shù)的相關(guān)研究鮮有報道。[10]本研究通過黑松立木結(jié)構(gòu)特征因子的實測，基于立木地徑、胸徑、冠幅、樹高等特征因子的關(guān)系分析，構(gòu)建黑松立木結(jié)構(gòu)特征參數(shù)模型，以期為區(qū)域森林資源調(diào)查與監(jiān)測、森林資源保護與管理等提供參考。

1 研究區(qū)概況

嶗山位于山東半島南部的黃海之濱，地處36°05′—36°19′ N，120°24′—120°42′ E，是中國海岸線第一高峰，主峰海拔1 132.7 m。嶗山屬暖溫帶季風氣候，年平均氣溫11.9 ℃，年日照平均時數(shù)2 622.3 h，年平均降雨量約800 mm。嶗山土壤的成土母巖主要是中生代花崗巖酸性巖類及噴發(fā)熔巖基性巖類，土壤多為棕色森林土，以沙壤為主。嶗山自然植被屬暖溫帶落葉闊葉林區(qū)，原以天然赤松林、櫟類次生林為主，后經(jīng)多年改造和更新造林，目前分布的主要森林類型為黑松林、赤松林、落葉松林、刺槐林、麻櫟林以及部分雜木林和灌叢等。黑松林多分布于山地陽坡、半陽坡、低山丘陵直至海岸線。

2 材料與方法

選擇無生長缺陷，樹干生長正常，胸徑大于5.0 cm的黑松作為樣木。樣木選定后，使用測徑圍尺測量地徑(D0)、距地表10 cm直徑(D0.1)和胸徑(D1.3)，運用布魯萊斯測高器測量樹高(H)，用塔尺測量枝下高(h)，用皮尺測量東西冠幅和南北冠幅(EW+NS)并計算其平均冠幅(CW)。共測定黑松樣木191株，地徑D0分布范圍7.1～44.1 cm，胸徑D1.3分布范圍5.1～38.5 cm，樹高H分布范圍3.1～20.2 m。

首先，運用SPSS Statistics 19.0的相關(guān)分析功能，利用實測的191株樣木數(shù)據(jù)，確定林木結(jié)構(gòu)指標之間的相關(guān)關(guān)系。對191株樣木數(shù)據(jù)進行整理抽樣，按2 cm徑階進行整化，每個徑階按順序抽3個隔1個進行回歸模型的建立，遴選出建模樣本144株，檢驗樣本47株。然后，運用Origin Pro 8進行各指標之間的回歸分析，構(gòu)建各指標間的回歸模型。擬合過程中，根據(jù)繪制模型預估值的標準殘差圖，將超出±2倍標準差以外的數(shù)據(jù)作為極端觀測值予以剔除。原則上，選擇調(diào)整的決定系數(shù)(Adj.R2)最大、殘差平方和(RSS)和回歸均方(MSR)最小的模型，作為最優(yōu)樹高曲線模型。運用SPSS Statistics 19.0配對樣本t檢驗(Pair Sample t test)功能，對檢驗樣本的實測值與各模型推算的理論值進行差異顯著性檢驗，若檢驗結(jié)果t0.05，則說明模型的預測精度滿足要求。

3 結(jié)果與分析

3.1 各指標之間的相關(guān)性分析

相關(guān)分析結(jié)果表明(表1)，在P=0.01水平(雙側(cè))上，黑松立木結(jié)構(gòu)的7個指標之間均存在極顯著的正相關(guān)關(guān)系，其中以D0和D0.1的相關(guān)系數(shù)最大，其Pearson相關(guān)系數(shù)為0.986，其次為D0與D1.3，而h與WE的相關(guān)性相對較弱，其Pearson相關(guān)系數(shù)僅為0.531，但仍大于Pearson相關(guān)系數(shù)臨界值0.185(P=0.01，df=190)。

表1 黑松立木結(jié)構(gòu)指標之間的相關(guān)性

3.2 不同部位直徑之間的回歸模型

D1.3-D0建模樣木為144株，按照D1.3平均值±2倍標準差的標準未發(fā)現(xiàn)異常值。D1.3-D0之間的關(guān)系可以用線性式、指數(shù)式和邏輯斯蒂(Logistics)式等表示，其中邏輯斯蒂式的決定系數(shù)R2最小，為0.948 2，其他線性式和指數(shù)式的決定系數(shù)R2均為0.950 8(表2)。但線性式的RSS與MSR較小，且模型中所包含的參數(shù)較少，因而線性式是最優(yōu)的擬合方程(圖1)。模型預測精度檢驗結(jié)果表明，實測胸徑與預測胸徑的Pearson相關(guān)系數(shù)為0.982 5，差異顯著性檢驗t=-0.4750.05，說明線性式能準確地反應D1.3-D0之間的相關(guān)性。

表2 黑松D1.3-D0模型擬合結(jié)果比較

注：RSS-殘差平方和，MSR-回歸均方，下同

D1.3-D0.1建模樣木為144株，按照胸徑平均值±2倍標準差的標準未發(fā)現(xiàn)異常值。D1.3-D0.1之間的關(guān)系可用線性式、指數(shù)式和等表示。表3表明，在3種模型中，因多項式具有較小的RSS和MSR，且決定系數(shù)R2值較大，因而多項式的擬合優(yōu)度最佳(圖2)。對檢驗樣本的實測胸徑與預測胸徑進行配對樣本t檢驗，發(fā)現(xiàn)實測胸徑與預測胸徑的Pearson相關(guān)系數(shù)為0.969 6，差異顯著性檢驗t=0.112 70.05，說明多項式能反應D1.3-D0.1的相關(guān)性。

表3 黑松D1.3-D0.1模型擬合結(jié)果比較

圖2黑松D1.3-D0.1最優(yōu)模型擬合及殘差圖

D0-D0.1建模樣木為144株，按照胸徑平均值±2倍標準差的標準未發(fā)現(xiàn)異常值。D0-D0.1之間的關(guān)系可用線性式、指數(shù)式和多項式等表示。表4表明，在3種模型中，因多項式具有較小的RSS和MSE，且決定系數(shù)R2值較大，因而多項式為描述D0-D0.1關(guān)系的最佳模型(圖3)。對檢驗樣本的實測胸徑與多項式預測胸徑進行配對樣本t檢驗，發(fā)現(xiàn)實測胸徑與預測胸徑的Pearson相關(guān)系數(shù)為0.991，差異顯著性檢驗t=-0.358 20.05，說明多項式能反應D0-D0.1的相關(guān)性。但由于線性式對模型相對簡單，也可較準確地反應D0-D0.1之間的關(guān)系。

表4 黑松D0-D0.1模型擬合結(jié)果比較

圖3D0-D0.1最優(yōu)模型擬合及殘差圖

3.3 胸徑-樹高相關(guān)關(guān)系分析

按照樹高平均值±2倍標準差的準則剔除異常數(shù)據(jù)2株，最終保留142株黑松胸徑樹高數(shù)據(jù)用于建模。表5說明，D1.3-H之間的相關(guān)關(guān)系可用指數(shù)式、邏輯斯蒂(Logsitic)式和岡珀茨(Gompertz)式表示，其中指數(shù)式的擬合優(yōu)度最佳，其決定系數(shù)R2為0.768 9，樹高殘差均在2倍標準差-5.859 8～5.859 8范圍內(nèi)(圖4)。對檢驗樣本的實測樹高與指數(shù)式預測樹高進行配對樣本t檢驗，發(fā)現(xiàn)實測樹高與預測樹高的Pearson相關(guān)系數(shù)為0.872 9，差異顯著性檢驗t=0.802 40.05，說明指數(shù)式能反應胸徑樹高的相關(guān)性。

表5 黑松D1.3-H曲線擬合模型結(jié)果比較

圖4黑松D1.3-H最優(yōu)模型擬合及殘差圖

3.4 胸徑-冠幅相關(guān)關(guān)系分析

對建模樣本的胸徑-平均冠幅數(shù)據(jù)回歸分析，發(fā)現(xiàn)線性式、指數(shù)式和多項式均能較好地反應二者之間的線性關(guān)系(表6)。由于線性式模型相對簡單，且決定系數(shù)R2值較大，而RSS值與其他兩個模型基本一致，因此線性式是反應D1.3-CW之間關(guān)系的最佳模型，殘差圖結(jié)果顯示D1.3殘差在2倍標準差即-15.335～15.335范圍內(nèi)(圖5)。對47株檢驗樣本的實測胸徑與線性式預測胸徑進行配對樣本t檢驗，發(fā)現(xiàn)實測胸徑與預測胸徑相關(guān)系數(shù)為0.898 4，差異顯著性檢驗t≈00.05，說明指數(shù)式能反應胸徑樹高的相關(guān)性。

表6 黑松D1.3-CW曲線擬合模型結(jié)果比較

圖5黑松D1.3-CW最優(yōu)模型擬合及殘差圖

4 結(jié)論與討論

立木結(jié)構(gòu)特征是森林資源調(diào)查、評估及經(jīng)營管理的基礎，涵蓋立木不同部位的直徑、樹高、冠幅、干形等指標和因子。本研究通過對191株黑松立木結(jié)構(gòu)特征因子的實測，基于立木地徑、胸徑、冠幅、樹高等特征因子的關(guān)系分析，構(gòu)建了黑松立木結(jié)構(gòu)特征參數(shù)模型。研究發(fā)現(xiàn)，黑松立木結(jié)構(gòu)地徑(D0)、距地面10 cm直徑(D0.1)、胸徑(D1.3)、樹高(H)、枝下高(h)、冠幅(CW)等指標之間均存在極顯著的正相關(guān)關(guān)系，其中以不同部位直徑之間的相關(guān)性最強。

本研究擬合并比較了黑松立木結(jié)構(gòu)指標之間的回歸模型。首先，利用建模樣本數(shù)據(jù)，分別擬合了3種常用的模型，通過比較決定系數(shù)R2、殘差平方RSS和和回歸均方MSR，篩選出擬合效果最佳的模型。然后，利用檢驗樣本數(shù)據(jù)，對模型的擬合效果進行配對樣本t檢驗，驗證模型的預測效果。本研究發(fā)現(xiàn)，D1.3-D0、D1.3-D0.1、D0-D0.1之間的相關(guān)模型可分別用線性式、多項式和多項式表示，模型的決定系數(shù)R2分別為0.950 8、0.939 3和0.969 5。以往對針葉樹種不同部位直徑關(guān)系的研究結(jié)果也表明，胸徑D1.3、地徑D0之間存在顯著的線性關(guān)系，如杉木地徑D0與胸徑D1.3間存在顯著的線性關(guān)系[11,12]、沙地樟子松地徑D0.1與胸徑D1.3之間存在顯著的多項式線性關(guān)系[13]。

由于樹木生長受多種因素的制約，如樹種生物學特性、立地條件、林分密度、管理措施等因素的影響，不同樹種的胸徑樹高曲線模型存在較大差別，如七姊妹山黃杉胸徑樹高曲線為冪函數(shù)式[14]、帽兒山林場20個闊葉樹種的胸徑樹高曲線模型分別為冪函數(shù)式、Schumacher式、單分子式、Richards式、Gompertz式和Logistic式等[15]。本研究黑松胸徑樹高(D1.3-H)之間的關(guān)系用指數(shù)式擬合優(yōu)度最佳，決定系數(shù)R2為0.768 9，這與赤松的胸徑樹高模型比較一致[16]，可能是由于兩個樹種在生物學特性比較相似。

利用無人機的航空攝影測量影像可以有效分析常綠森林的樹冠結(jié)構(gòu)，張凝等[17]利用的無人機的攝影測量影像分割方法，發(fā)現(xiàn)在快速獲得單木冠幅信息的前提下冠幅提取精度達90.05%。然而，樹木胸徑數(shù)據(jù)無法直接測量，只能通過反演方法進行估測。本研究發(fā)現(xiàn)，黑松的胸徑D1.3與冠幅CW之間存在顯著的線性相關(guān)關(guān)系，其決定系數(shù)Adj.R2為0.805 8。檢驗樣本的驗證結(jié)果也表明，該模型能準確地反應胸徑D1.3與冠幅CW之間的關(guān)系。

本研究基于191株樣木的實測數(shù)據(jù)，通過擬合與比較分析，篩選出黑松立木結(jié)構(gòu)特征參數(shù)模型。但是，由于未對樣木所處的立地條件區(qū)別分類，如坡度、坡位、坡向等因子，對模型的預測精度有一定影響[18,19]。下一步的研究工作中，將不斷豐富完善數(shù)據(jù)，把立地條件如坡度、坡向、海拔等因子與樹木直徑、樹高和冠幅等結(jié)構(gòu)指標聯(lián)合建立估測模型來提高預測精度。

：

[1] Maltamo M,Mustonen K,Hyypp? J,et al. The accuracy of estimating individual tree variables with airborne laser scanning in a boreal nature reserve[J]. Canadian Journal of Forest Research,2004,34(9):1791-1801

[2] Weaver S A,Ucar Z,Bettinger P,et al. Assessing the Accuracy of Tree Diameter Measurements Collected at a Distance[J]. Croatian Journal of Forest Engineering,2015,36(1):73-83

[3] 姚愛靜,朱清科,張宇清,等.林分結(jié)構(gòu)研究現(xiàn)狀與展望[J].林業(yè)調(diào)查規(guī)劃,2005,30(2):70-76

[4] 閉海秀,黃清經(jīng).桂東區(qū)杉木地徑與胸徑、樹高、材積相關(guān)分析[J].廣西林業(yè)科學,2014(4):444-449

[5] 曾慎松.巨尾桉人工林地徑與胸徑、樹高相關(guān)模型的研究[J].福建林業(yè)科技,2007,34(1):40-42

[6] 王偉.無人機影像森林信息提取與模型研建[D].北京：北京林業(yè)大學,2015

[7] 李祥,鄭淯文,戴楚彥,等.基于無人機影像的森林信息獲取研究進展[J].世界林業(yè)研究,2017,30(4):41-46

[8] 李亞東,馮仲科,明海軍,等.無人機航測技術(shù)在森林蓄積量估測中的應用[J].測繪通報,2017(4):63-66

[9] 劉清旺,李世明,李增元,等.無人機激光雷達與攝影測量林業(yè)應用研究進展[J].林業(yè)科學,2017,53(7):134-148

[10] 許景偉,李傳榮,李琪,等.沿海黑松防護林立木材積表的編制[J].山東林業(yè)科技,2003(4):22-25

[11] 王金榮.杉木天然林林木地徑與胸徑的關(guān)系研究——以福建五一國有林場為例[J].林業(yè)資源管理,2015(3):121-124

[12] 曾偉生.利用誤差變量聯(lián)立方程組建立南方杉木一元立木材積模型和胸徑地徑回歸模型[J].中南林業(yè)調(diào)查規(guī)劃,2012,31(4):1-4

[13] 張日升,劉廣,于洪軍,等.遼西北沙地樟子松地徑與胸徑相關(guān)關(guān)系及其應用研究[J].防護林科技,2006(3):19-21

[14] 熊斌梅,汪正祥,李中強,等.七姊妹山自然保護區(qū)黃杉年齡胸徑樹高的相關(guān)性研究[J].林業(yè)資源管理,2016(4):41-46

[15] 王小明,李鳳日,賈煒瑋,等.帽兒山林場天然次生林闊葉樹種樹高-胸徑模型[J].東北林業(yè)大學學報,2013(12):116-120

[16] 徐悅,陳昌華,蔣之富,等.天然赤松胸徑與樹高相關(guān)模型的研究[J].林業(yè)調(diào)查規(guī)劃,2008,33(3):56-58

[17] 張凝，馮躍文，張曉麗，等.結(jié)合航空影像紋理和光譜特征的單木冠幅提取[J].北京林業(yè)大學學報，2015，37(3):13-19

[18] 袁曉紅,李際平.杉木人工林南北坡向樹高-胸徑生長曲線研究[J].西北林學院學報,2012,27(2):180-183

[19] 杜紀山,李希菲.杉木樹高曲線參數(shù)穩(wěn)定性的研究[J].林業(yè)科學研究,1998,11(2):169-174