偏冠对树冠垂直投影面积计算的影响

孙 钊，潘 磊，乔晶晶，丁志丹，孙玉军，王轶夫

(北京林业大学森林资源和环境管理国家林业和草原局重点开放性实验室,北京 100083)

树冠信息是不可少的林分调查因子，反映了林木长期竞争水平[1-2]。无论是林业工作者还是研究人员越来越重视树冠信息的获取，随着航空航天遥感的发展，林业工作持续向精细化推进[3]。树冠冠幅是光学遥感测量的重点，准确测量地面树冠垂直投影面积对遥感冠幅提取的精度验证有较大帮助。

在一个开放、没有竞争、养分充足的理想环境中，树木通常会生长出非常对称的冠形，树冠的垂直投影区域会一直趋近于圆形[4-6]。但是实际环境中，为了获取更大的生长优势，无论是天然林还是人工林总会存在对生存环境的竞争现象，对于光线的竞争直接影响了树冠冠形，而此时的林木通常会形成不规则不对称的冠形，树冠偏向更为优势的生长空间，通常这种情况叫做树冠的偏冠现象[7]。由于偏冠现象的存在，林木不同方向上活枝长度会有较大差异，导致树冠垂直投影区域通常不会呈正圆形，传统的树冠投影面积计算方法是测量树冠活枝边缘投影到地面的位置到树干的水平距离[8-9]，计算树冠不同方向的平均冠幅，将树冠视作正圆和椭圆计算出树冠垂直投影面积，可能会与真实面积存在一定差异。地面测量不同方向冠幅后绘制的树冠垂直投影区域图像与现在光学遥感影像提取的树冠冠幅区域对比是当前判定树冠提取精度的重要手段。

怎样通过获取的森林遥感影像有效的进行森林参数提取是当前研究的重点。单木分割，树种分类等研究都需要对整幅影像进行分割处理，只有在地面测量精度较高的情况下，以此来进行林业遥感的地面验证才具有实际意义，所以怎样获取真实冠幅面积的重要性大大增加，无人机遥感进行树冠冠幅提取的精度验证通常是对比单木东西冠幅，南北冠幅[10-12]或利用东西冠幅，南北冠幅进行树冠垂直投影面积[13-15]计算，对树冠采用近似圆的方法进行树冠面积计算时，由于实际林木生长情况的不同，林内竞争压力的不同，从而产生树冠偏冠现象，可能导致最终计算的树冠投影面积与真实值相差较大，因此，直接对树冠投影形状进行面积测量更为准确，本研究提出一种直接对树冠投影区域面积计算的方法，并研究了偏冠对于树冠投影面积计算的影响，证明基于树冠投影区域形状的面积计算能够更加准确的计算树冠投影区域面积，为林业遥感的进一步发展提供一定的参考。

1 研究区概况

研究区为将乐国有林场（26°40′～26°50′ N，117° 25′～117°35′ E），位于福建省三明市将乐县，林场最高海拔1 203 m，最低海拔140 m，年均气温18.7℃，年均降 水1 669 mm，年均蒸发量1 204 mm。杉木（Cunninghamia lanceolata(Lamb.)Hook.）、马尾松（Pinus massonianaLamb.）是将乐国有林场主要针叶树种。

2 材料与方法

2.1 数据获取

选择生长状况良好的近熟和成熟杉木纯林样地（20 m×30 m）共5 块，样地平缓无坡度，测量样地内共453 株杉木不同方向冠幅数据。利用激光测距仪实测每木8 向（东、东南、南、西南、西、西北、北、东北）最大冠长到树干的距离，测量时，保持测量位置与活冠边界的充分一致，枯枝不记作树冠边界（图1）。

图1 树冠半径测量示意图Fig.1 Schematic diagram of crown radius measurement

在Visualstudio 中实现测量距离到点坐标的转换，便于在ArcGIS 中构建树冠垂直投影多边形。树冠测量半径坐标转换实现示意如图2。

图2 树冠半径转换为坐标示意图Fig.2 Diagram of transforming crown radius into coordinates

坐标转换公式：

以样地西南角为坐标原点，x为树干到样地西边界的距离，y为树干到样地南边界的距离，0′A为树冠东北向冠幅长度。

通过ArcGIS 点转面生成的树冠冠幅，点与点的连接线为直线，由于实际生长状况下，树冠边界基本不存在直线形式，为了使生成结果更加接近真实情况，需要对边界直线进行平滑，本研究选择Bezier 曲线进行树冠边界的平滑，Bezier 曲线是一种在图形平滑处理中广泛应用的曲线模型，在常用的图形处理软件中可以根据需求控制曲线形状[16]。Bezier 曲线通过调整图形边界上控制点的位置从而控制曲线弧度，将怎么绘制曲线的问题转换为控制点坐标的定位，首尾连接控制点构成光滑曲线，利用Bezier 公式趋近于原始多边形[17]（图3），从而得到所需的理想图形。随着控制点个数的增加，Bezier 曲线的复杂程度和灵活程度逐渐增加，通常使用的是包含3 个或4 个控制点的Bezier 曲线，N阶Bezier 曲线的拟合公式[18]为

图3 高阶Bezier 曲线示意图Fig.3 High order Bezier curve schematic

式中Pi为 Bezier 曲线的n+1个控制点，Bi,n(t)为Bernstein 基函数，有

2.2 树冠偏冠指数计算方法

在树冠偏冠的研究中通常采用数学“圆度”来衡量冠幅是否对称，冠幅投影形状越接近于圆则认为该树冠越对称[19]。以树干位置为中心，由测量的一组树冠冠幅长度构成的不规则图形被视为树冠投影区域，该图形与圆的偏离程度视作树冠偏冠指数CAI（Crown Asymmetry Index）[20]。在数学图像学领域，大量学者作了任意图形“圆度”测量的研究，本研究结合已有的树冠冠幅数据，参考孔繁琳[16]的研究，选择了其结合Herrera-Navarro[21]等人在2013年提出的数学图形平均圆度测量方法改编的树冠偏冠指数计算方法。CAI的计算方法如下：

Ri为 第i个树冠半径的长度，R为N个树冠半径的 平均长度，N为所测量树冠半径的总数。

3 结果与分析

3.1 树冠垂直投影面积计算

Bezier 曲线现在广泛运用于各类图形处理软件中，在ArcGIS 中的Smooth polygon 命令中选择Bezier 算法实现树冠投影面积平滑，结果如图4。

由图4 可以看出，相较于直线连接的树冠边界，经过Bezier 曲线平滑之后的树冠垂直投影面积更加真实树冠投影形状。所以选择Bezier 曲线平滑后的图形面积作为真实树冠垂直投影面积。

图4 树冠垂直投影面积计算结果（a.冠幅直线连接,b.Bezier 曲线平滑连接，c.两种冠幅叠加比较）Fig.4 Calculation results of vertical projection area of tree crown (a.Straight line connection of crown width,b.smooth connection of Bezier curve,c.superposition comparison of two kinds of crown width)

3.2 树冠垂直投影面积和传统计算方法的比较

在ArcGIS 中生成的树冠面图形，可以通过属性中的面积字段，直接计算当前面的面积。传统的树冠垂直投影面积计算，是将树冠投影图像作为圆或者椭圆，以平均冠幅作为圆的半径，以东西或南北冠幅作为长轴或短轴。本研究分别选择4 向和8 向冠幅的正圆和椭圆面积计算方法，结果见图5。

由图5 可以看出，8 向冠幅正圆面积计算方法和椭圆面积计算方法和树冠垂直投影面积所得RMSE分别为0.423 0 m2和0.743 3 m2，4 向冠幅正圆面积计算方法和椭圆面积计算方法和树冠垂直投影面积所得RMSE分别为1.368 5 m2和1.390 5 m2。可以看出，测量冠幅半径数量越多，计算结果越接近真实冠幅面积，同时，相比于椭圆的计算方法，把树冠视为正圆计算的结果更加接近真实树冠垂直投影面积。

图5 不同树冠垂直投影面积计算方法和树冠垂直投影面积差异Fig.5 Different calculation methods of canopy vertical projection area and difference of canopy vertical projection area

3.3 树冠偏冠指数和冠幅面积计算差值的相关性

根据式5 计算出每个树冠的CAI值，由CAI计算方法可以看出，树冠不同方向的冠幅长度直接影响了树冠偏冠指数的结果，将不同树冠垂直投影面积计算方法所得结果和树冠垂直投影面积差值与CAI进行相关性分析，结果见图6，表1。

图6 每木CAI 计算结果Fig.6 CAI results of each tree

由表1 可以看出，不同树冠垂直投影面积计算方法所得结果和树冠垂直投影面积差值与CAI呈显著正相关。在计算树冠垂直投影面积时，正是由于树冠偏冠的存在，才导致了不同面积计算方式与真实面积的差异。

表1 CAI 和树冠垂直投影面积计算差值相关性Table 1 Correlation between CAI and vertical projection area of tree crown

4 结论

树冠垂直投影面积的精准测量对无人机调查精度意义重大[22]。本研究以研究区453 棵杉木为研究对象，通过实地测量单木8 向最大冠长到树干的距离，以此作为单木8 向冠幅，以Visualstudio 和ArcGIS 混合编程，实现树冠测量半径的坐标转换，为了更加贴近真实树冠形状，采用Bezier 曲线进行边界平滑，得到树冠垂直投影面积，分别与将树冠视作正圆和椭圆的计算方法进行对比分析。结果表明，本研究提出的树冠垂直投影面积计算方法能直接提取出树冠垂直投影面积；树冠垂直投影面积计算方法中，测量冠幅半径数量越多，计算结果越接近真实冠幅面积；相比于椭圆的计算方法，把树冠视为正圆所得面积更接近真实面积。

树冠不对称性是近年来森林生态研究的热点，本研究参考前人研究，以“圆度”来衡量树冠偏冠，并计算了CAI与树冠垂直投影面积计算差值的相关性。结果表明，不同树冠垂直投影面积计算方法所得结果和树冠垂直投影面积差值与CAI呈显著正相关。因此认为，在计算树冠垂直投影面积时，由于树冠偏冠现象的存在，导致了传统的不同面积计算方式与树冠真实面积的差异。

本研究中，虽然比较了4 向和8 向冠幅结果，但为了更进一步证实冠幅测量数量对树冠垂直投影面积计算的影响，应进一步测量16 向或32 向冠幅结果，以确定冠幅测量数量对面积计算的影响。本研究应用Bezier 曲线进行树冠边界平滑，但以此代表实际结构更加复杂，形状不规则的树冠远远不够，如何更加真实还原真实树冠，更贴近森林真实样貌有待进一步研究。树冠不对称性是树冠形态的一个特征属性，树冠垂直投影面积的精准测量，将会极大的促进林业遥感和人工智能在林业领域的发展。