发明名称 测量和计算森林蓄积参数的方法
摘要 一种测量和计算森林蓄积参数的方法,属于用全景扫描的激光测距传感器对林地测量和计算森林蓄积参数的方法,按如下七个步骤进行:一是用激光测距传感器扫描林分,二是获取胸径,三是计算每公顷胸高断面积,四是计算平均胸径,五是树干密度,六是计算不同高度的植被覆盖率,七是求得林木分布空间的坐标格局;采用本非接触方式测量和计算森林蓄积参数,快捷、方便、准确、效率高、劳动强度低,同时能提供现有方法还不能提供的垂直方向不同高度的植被覆盖率,用于与航空遥感和航天遥感信息相结合估计森林蓄积和森林生物量等优点。
申请公布号 CN103616014B 申请公布日期 2016.03.30
申请号 CN201310633096.0 申请日期 2013.11.29
申请人 浙江农林大学 发明人 张茂震
分类号 G01C11/00(2006.01)I;G01B11/08(2006.01)I;G01B11/28(2006.01)I 主分类号 G01C11/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 周烽
主权项 一种测量和计算森林蓄积参数的方法,其特征是按如下步骤进行:(1)用激光测距传感器扫描林分:在森林中,将全景扫描仪安装于三脚架上1.3米高位置,调整至水平后,开始360度环绕扫描,得到周围所有树木与扫描仪之间的距离;(2)获取胸径:逐一扫描每一株树,定O为扫描仪位点,C为树干中心点,扫描仪匀速步进,其相邻两个扫描点之间的夹角ω相等,即∠EOG=∠GOC<sub>2</sub>=∠C<sub>2</sub>OH=∠HOF=ω,ω由步进电机转速和激光传感器扫描速率确定;以OE和OF之间的平分线OC为斜边构成的直角三角形EOC和COF相等,有∠EOF=α,<img file="FDA0000834389120000011.GIF" wi="605" he="135" />设R=OC,R'=OE=OF,则树木胸径<img file="FDA0000834389120000012.GIF" wi="374" he="133" />EC、FC相等,为<img file="FDA0000834389120000013.GIF" wi="85" he="148" />α=kω,k为在此树干外表实际扫描点数;对第i株树表面,从E点经G、C<sub>2</sub>、H到F逐点扫描,其返回距离从OE、OG、OC<sub>2</sub>到OH、OF,由大到小、再由小到大变化,当达到相等时表示对第i株树扫描结束;第i株树胸径d<sub>i</sub>、扫描点O距树木中心点C的距离R<sub>i</sub>以及扫描点O距树干上E点(或F点)的距离R'<sub>i</sub>的关系是:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><msub><mi>R</mi><mi>i</mi></msub><mo>&prime;</mo></msup><mo>=</mo><msqrt><mrow><msup><msub><mi>R</mi><mi>i</mi></msub><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>d</mi><mi>i</mi></msub><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000834389120000014.GIF" wi="453" he="94" /></maths>(3)、计算每公顷胸高断面积:根据角规原理,用d<sub>i</sub>和R<sub>i</sub>'计算每公顷胸高断面积;①计算判定标准:因行业内公用的角规杆长500mm,标准缺口宽10mm,即OC'=500mm,A'B'=10mm;根据相似三角形原理,<img file="FDA0000834389120000026.GIF" wi="356" he="126" />让视线从O点出发,通过缺口A'B'观察树干,能与缺口相切或相割的树干直径d应大于或等于一个标准直径d<sub>标</sub>,有:<img file="FDA0000834389120000021.GIF" wi="453" he="143" />因<img file="FDA0000834389120000022.GIF" wi="332" he="149" />则有:<img file="FDA0000834389120000023.GIF" wi="238" he="127" />这里的d<sub>标</sub>随距离OE变化;②角规计数:对于每一株被扫描样木,计算d<sub>标</sub>值,并用之与被扫描样木的实际胸径d<sub>i</sub>比较,当胸径d<sub>i</sub>大于d<sub>标</sub>时计1,当胸径d<sub>i</sub>等于d<sub>标</sub>时计0.5,否则计0,即:<img file="FDA0000834389120000024.GIF" wi="630" he="277" />角规计数汇总即为林分每公顷胸径断面积G;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>G</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>g</mi><mi>i</mi></msub></mrow>]]></math><img file="FDA0000834389120000025.GIF" wi="334" he="142" /></maths>单位:m<sup>2</sup>/hm<sup>2</sup>;(4)计算平均胸径:在步骤(2)的基础上,求算所测范围内林木的平均胸径d;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>d</mi><mo>&OverBar;</mo></mover><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>d</mi><mi>i</mi><mn>2</mn></msubsup></mrow></msqrt></mrow>]]></math><img file="FDA0000834389120000031.GIF" wi="378" he="159" /></maths>(5)计算树干密度:树干密度是林木胸高断面积与相应范围的土地面积之比;根据步骤(3)测得林分每公顷胸高断面积G,计算树干密度(以万分数表示):<img file="FDA0000834389120000032.GIF" wi="587" he="71" />(6)计算植被覆盖率:根据步骤(1)获得所有对象的距离,这些对象包括达到胸高位置的乔木、灌木植被,扫描返回植被数据的个数m<sub>1.3</sub>与全部扫描数据个数n<sub>1.3</sub>之比为胸高位置的植被覆盖率VC<sub>1.3</sub>,并以百分数表示;VC<sub>1.3</sub>=m<sub>1.3</sub>/n<sub>1.3</sub>*100%其中,下标1.3表示扫描高度为1.3米,即可在不同高度扫描获得不同高度的植被覆盖率;(7)求得林木分布空间坐标格局:根据步骤(2),在已有每株树木胸径和树木与扫描仪之间距离的基础上,按扫描仪的步进转速ω和激光传感器测量速率f,确定以扫描仪为圆心的每株树相对空间坐标(θ,R);其中,θ为树木中心点距扫描起始线的角度,R为树木中心点距扫描仪中心点的距离;由于步进电机每前进1步激光传感器扫描1次,水平转动角度为ω,如果扫描仪自扫描开始到完成第i株树干上的最后一个点扫描,共接收n次扫描数据,即步进电机运行了n步,则该树干表面最后一个扫描点相对于扫描起始点的角度为ω×n,该树干中心相对于起始点的角度为θ<sub>i</sub>=ω×n‑α;距离<img file="FDA0000834389120000041.GIF" wi="534" he="110" />第i株树的坐标为(θ<sub>i</sub>,R<sub>i</sub>)。
地址 311300 浙江省杭州市临安市环城北路88号