CN103175491B - 一种树木测量方法及树木测量系统 - Google Patents
一种树木测量方法及树木测量系统 Download PDFInfo
- Publication number
- CN103175491B CN103175491B CN201210249851.0A CN201210249851A CN103175491B CN 103175491 B CN103175491 B CN 103175491B CN 201210249851 A CN201210249851 A CN 201210249851A CN 103175491 B CN103175491 B CN 103175491B
- Authority
- CN
- China
- Prior art keywords
- point
- axis coordinate
- segmental arc
- tree
- coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
本发明公开了一种树木测量方法,包括:处理器获取无棱镜全站仪实时采集的采样点信息,所述采样点信息包括最高点高程Ht、最低点高程Hb、采样点坐标;所述处理器根据所述最高点高程Ht及最低点高程Hb计算树高H,所述树高H=Ht-Hb;所述处理器根据所述采样点坐标计算树冠投影面积。本发明还公开了一种树木测量系统。采用本发明,可在地形复杂的条件下,利用无棱镜全站仪高精度、高效率地采集目标树木的采样点信息,并通过处理器计算树高、树冠投影面积,同时,计算树冠投影面积时,根据树冠总体生长态势,结合凹凸多边形,把非线性的树冠投影线性化,将树冠投影转换为具有凹弧、凸弧的弧段多边形进行计算,精确度高。
Description
技术领域
本发明涉及森林资源监测领域,尤其涉及一种树木测量方法及树木测量系统。
背景技术
在传统的森林资源监测中,精确测量树高和树冠投影面积一直是一个难以解决的问题。
通常,树高主要使用测高器、罗盘仪等进行测量,但使用该方法测量时,数据无法实时存储与解算,仅能依靠手工记录数据,精度低、工作量大。
而树冠投影面积则多将其视为圆形,通过拉尺测量冠径来求算面积,但该测量方法采用目估法,手工操作,精度较低。
另外,北京林业大学的冯仲科等人曾提出利用有棱镜全站仪,采用三角高程法、三维前方交会法、全站仪解析法等方法进行树高和树冠的测量。但该测量方法中,有棱镜全站仪必须在树底下设置棱镜,通过反射获得数据。而在野外测量中,树下环境不一定适宜架设棱镜,所以这种方法耗时,效率低,并且受地形影响大。在计算树冠投影面积时,利用近圆或者近椭圆模拟不规则的树冠投影面,通过有棱镜全站仪测量树冠的冠幅以得到冠径,再求出树冠面积。因此,利用此类测量方法得出来的树冠投影面积与实际面积误差较大,精度低。
所以,在森林资源监测过程中,传统的树高、树冠投影面积测量方法易受森林中复杂地形因素的限制,且工作量大、自动化程度低、精度不高。因此,实现树高和树冠投影面积的高精度、高效率的测量成为本领域科技人员急待解决的问题。
发明内容
本发明所要解决的技术问题在于,提供一种树木测量方法及树木测量系统,可在地形复杂的条件下,利用无棱镜全站仪高精度、高效率地采集目标树木的采样点信息,并通过处理器计算树高、树冠投影面积,同时,计算树冠投影面积时,根据树冠总体生长态势,结合凹凸多边形,把非线性的树冠投影线性化,将树冠投影转换为具有凹弧、凸弧的弧段多边形进行计算,精确度高。
为了解决上述技术问题,本发明提供了一种树木测量方法,包括:处理器获取无棱镜全站仪实时采集的采样点信息,所述采样点信息包括最高点高程Ht、最低点高程Hb、采样点坐标,所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标;所述处理器根据所述最高点高程Ht及最低点高程Hb计算树高H,所述树高H=Ht-Hb;所述处理器根据所述采样点坐标计算树冠投影面积。
作为上述方案的改进,所述处理器根据采样点坐标计算树冠投影面积的步骤包括:提取所述采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形;根据所述基于线段的多边形生成基于弧段的多边形;计算所述基于弧段的多边形的面积。
作为上述方案的改进,所述提取采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形的步骤包括:提取所述采样点坐标中的X轴坐标及Y轴坐标,根据所述X轴坐标及Y轴坐标将所述采样点绘于同一图层内;相邻的采样点两两相连,组合成基于线段的多边形。
作为上述方案的改进,所述根据基于线段的多边形生成基于弧段的多边形的步骤包括:获取预设的样品参数N;计算所述相邻的采样点间的弦长L;根据所述样品参数N及所述相邻的采样点间的弦长L,计算所述相邻的采样点间对应的弧段的弦高D,所述弦高D=N×L;计算由相邻的三个采样点依次相连而成的内角的角度,所述三个采样点依次为起点、中点、终点;判断所述内角的角度是否大于180°,判断为是时,则所述中点与终点间对应的弧段为凹弧,判断为否时,则所述中点与终点间对应的弧段为凸弧;根据所述弦高D、弦长L及弧段的凹凸性生成所述相邻的采样点间的弧段。
作为上述方案的改进,所述处理器获取无棱镜全站仪实时采集的采样点信息的步骤之后还包括存储所述采样点信息。
相应地,本发明还提供了一种树木测量系统,包括:无棱镜全站仪,用于实时采集目标树木的采样点信息,所述采样点信息包括最高点高程Ht、最低点高程Hb、采样点坐标,所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标;处理器,用于根据所述无棱镜全站仪实时采集的采样点信息计算树高及树冠投影面积。
作为上述方案的改进,所述处理器包括:获取单元,用户获取所述无棱镜全站仪实时采集的采样点信息;与所述获取单元相连的树高计算单元,用于根据所述无棱镜全站仪实时采集的采样点信息计算树高;与所述获取单元相连的树冠投影面积计算单元,用于根据所述无棱镜全站仪实时采集的采样点信息计算树冠投影面积。
作为上述方案的改进,所述树冠投影面积计算单元包括:线段多边形单元,提取所述采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形;与所述线段多边形单元相连的弧段多边形单元,根据所述基于线段的多边形生成基于弧段的多边形;与所述弧段多边形单元相连的计算单元,用于计算所述基于弧段的多边形的面积。
作为上述方案的改进,所述弧段多边形单元包括:参数获取子单元,用于获取预设的样品参数N;弦长计算子单元,用于计算所述相邻的采样点间的弦长L;与所述参数获取子单元及弦长计算子单元分别相连的弦高计算子单元,用于根据所述样品参数N及所述相邻的采样点间的弦长L,计算所述相邻的采样点间对应的弧段的弦高D,所述弦高D=N×L;内角计算子单元,用于计算由相邻的三个采样点依次相连而成的内角的角度;与所述内角计算子单元相连的判断子单元,用于判断所述内角的角度是否大于180°并做相应处理;与所述弦高计算子单元及判断子单元分别相连的弧段生成子单元,用于根据所述弦高D、弦长L及判断子单元的判断结果生成所述相邻的采样点间的弧段。
作为上述方案的改进,所述处理器还包括与所述获取单元相连,用于存储所述采样点信息的存储单元。
实施本发明的有益效果在于:可在地形复杂的条件下,利用无棱镜全站仪高精度、高效率地采集目标树木的采样点信息,通过处理器计算树高、树冠投影面积。其中,无棱镜全站仪具有测量精度高、测量速度快、测量距离远等优点,特别适合于人员无法到达的地区,受地形影响小,实现了“所瞄即所测”,效率高,适用范围广,安全性好,操作简单。另外,处理器自动计算,免去了手工计算出图,减少了工作量,提高了工作效率。
计算树高时,仅需选取目标树木顶端枝条或者与该枝条处于同一水平面上的树叶为参考点,测量最高点高程Ht,同时,选取与目标树木根部处于同一水平面的有效测区为参考点,测量最低点高程Hb,处理器获取最高点高程Ht及最低点高程Hb以计算树高H,所述树高H=Ht-Hb。因此,不用测量无棱镜全站仪至树干中心的平距,克服了在测区内有陡坡或沟壑无法测量水平距离的缺点。
同时,计算树冠投影面积时,根据树冠总体生长态势,选取树冠在各个方向的最大伸展处作为采样点,利用无棱镜全站仪测量采样点的三维坐标。处理器把非线性的树冠投影线性化,根据采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形,并结合凹凸多边形的特点,将树冠投影转换为具有凹弧、凸弧的弧段多边形进行计算,精确度高。
附图说明
图1是本发明一种树木测量方法的流程示意图;
图2是图1中处理器根据采样点坐标计算树冠投影面积的流程示意图;
图3是图1中处理器根据采样点坐标计算树冠投影面积的另一流程示意图;
图4是本发明一种树木测量方法的第一实施例中基于线段的多边形;
图5是本发明一种树木测量方法的第一实施例中基于圆弧的多边形;
图6是本发明一种树木测量方法的第二实施例中基于线段的多边形;
图7是本发明一种树木测量方法的第二实施例中基于圆弧的多边形;
图8是本发明一种树木测量方法的第三实施例中基于线段的多边形;
图9是本发明一种树木测量方法的第三实施例中基于圆弧的多边形;
图10是本发明一种树木测量系统的结构示意图;
图11是本发明一种树木测量系统中处理器的结构示意图;
图12是本发明一种树木测量系统中处理器的另一结构示意图;
图13是本发明一种树木测量系统的处理器中树冠投影面积计算单元的结构示意图;
图14是图13中弧段多边形单元的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
图1是本发明一种树木测量方法的流程示意图,包括:
S100,处理器获取无棱镜全站仪实时采集的采样点信息。
所述采样点信息包括最高点高程Ht、最低点高程Hb、采样点坐标。
所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标。
S101,所述处理器根据所述最高点高程Ht及最低点高程Hb计算树高H。
需要说明的是,处理器利用公式H=Ht-Hb,计算树高。
计算树高时,仅需选取目标树木顶端枝条或者与该枝条处于同一水平面上的树叶为参考点,并利用无棱镜全站仪测量最高点高程Ht。同时,选取与目标树木根部处于同一水平面的有效测区为参考点,并利用无棱镜全站仪测量最低点高程Hb。无棱镜全站仪将测量所得的最高点高程Ht及最低点高程Hb上传至处理器。处理器获取最高点高程Ht及最低点高程Hb,并依据公式H=Ht-Hb计算树高H。
需要说明的是,利用无棱镜全站仪测量时,需要考虑风速、植物通视程度、枝条粗细等因素。当枝条受风速影响小而且通视情况良好时,可直接测量植物顶端枝条以确定最高点高程Ht;否则,应测量与该枝条处于同一水平面上的叶子,且应选取枝条摆动弧度上限位置为测量点以确定最高点高程Ht。
S102,所述处理器根据所述采样点坐标计算树冠投影面积。
测量树冠投影面积时,根据树冠总体生长态势,初步选取树冠在各个方向的最大伸展处作为采样点,再根据树冠投影形状和凹凸点分布情况最终确定采样点,利用无棱镜全站仪顺时针或逆时针测量采样点的三维坐标。
处理器根据采样点坐标中的X轴坐标及Y轴坐标,描绘采样点并根据采样点模拟出与树冠投影相似的不规则多边形,此时,计算出的多边形面积可近似为树冠投影的面积,大大提高了测量的精度。
需要说明的是,由于地形、通视情况等因素影响,往往需要布设多个控制点,并在近地面风速小于0.5m/s时,实现激光对准目标树木的情况下,从不同角度对树冠进行测量,部分叶片厚重的树木,可以在风速小于1m/s测量。
需要说明的是,步骤S101与步骤S102之间没有必然的先后顺序,在进行步骤S101的同时也可以进行步骤S102。
图2是图1中处理器根据采样点坐标计算树冠投影面积的流程示意图,包括:
S200,提取所述采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形。
测量树冠投影面积时,根据树冠总体生长态势,初步选取树冠在各个方向的最大伸展处作为采样点,再根据树冠投影形状和凹凸点分布情况最终确定采样点,利用无棱镜全站仪顺时针或逆时针测量采样点的三维坐标。
所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标。
处理器获取无棱镜全站仪实时采集的采样点坐标,依次提取采样点坐标中的X轴坐标及Y轴坐标,根据X轴坐标及Y轴坐标,描绘采样点并连接成线段多边形。
S201,根据所述基于线段的多边形生成基于弧段的多边形。
在基于线段的多边形的基础上,将基于线段的多边形转换为具有凹弧、凸弧的基于弧段的多边形,使模拟的多边形更接近树冠投影的形状,大大提高了测量的精度。
S202,计算所述基于弧段的多边形的面积。
计算基于弧段的多边形的面积时,可利用积分原理对弧段多边形进行面积计算,即通过计算上下曲线/折线的方程,再求上下曲线积分得到面积。
也可以使用多边形像素映射方法计算基于弧段的多边形面积,即先计算多边形内的像素数量,再按比例换算成面积。
还可以采用基于扫描线方法的面积算法计算基于弧段的多边形,即通过从每一个采样点引出的横向扫描线,把整个多边形区域划分为若干个三角形、梯形(只要存在两条对边平行的四边形就称为梯形)、或者是弧形(在存在圆弧的情况下,先将圆弧的弦连起来作为一个弧形单独考虑),从而分别求算各个分块的面积,各个分块累加即可得到整个基于弧段的多边形的面积。
例如(参见图5),多边形为含有弧形边的不规则图形,采用基于扫描线方法计算面积时,先对所有的弦和直线边组成的多边形进行划分,计算获得对应的面积块累加公式。其中,对多边形排序后的边集合为{(A,B)(B,G)(G,C)(C,D)(D,E)(E,F)(F,A)},另外,弧(A,B)(B,C)(D,E)(E,A)围成的面积为凸出于多边形的弧形,以及弧(C,D)围成的面积凹入多边形的弧形。然后,根据弧所对应的凹凸关系,对应的总面积减去或加上圆弧围成的面积。
图3是图1中处理器根据采样点坐标计算树冠投影面积的另一流程示意图,包括:
S300,提取所述采样点坐标中的X轴坐标及Y轴坐标,根据所述X轴坐标及Y轴坐标将所述采样点绘于同一图层内。
所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标。
确定采样点后,利用无棱镜全站仪顺时针或逆时针测量采样点的三维坐标。处理器获取无棱镜全站仪实时采集的采样点坐标,并按照无棱镜全站仪的测量顺序依次提取采样点坐标中的X轴坐标及Y轴坐标,根据X轴坐标及Y轴坐标将所述采样点绘于同一图层内。
S301,相邻的采样点两两相连,组合成基于线段的多边形。
相邻的采样点,依照测量顺序通过线段形式两两相连,形成基于线段的多边形。
S302,获取预设的样品参数N。
需要说明的是,样品参数N为相邻样品点的弦高与弦长之比的平均值。
在太阳直射点上,工作人员进行实地测量,选定样品树木,并在地面上标出样品树木的树冠投影,分别选取树冠投影上的两个相邻的样品点,沿树冠投影勾画出两个相邻的样品点间的弧段。用皮尺测量出两样品点的距离(即弦长)以及两样品点间连线的垂直平分线的长度(即弦高),此时,将两样品点间连线的垂直平分线的长度与两样品点的距离之比称为固定比。经上千次的反复测量,分别选取不同种类的样品树木及不同的样品点,计算可得出上千个固定比。此时,计算各固定比的平均值即称为样品参数N。
优选地,所述样品参数N的值为1/4.15。
S303,计算所述相邻的采样点间的弦长L。
根据步骤S301组合成的基于线段的多边形,依照测量顺序计算基于线段的多边形内相邻的采样点间的弦长L。
S304,根据所述样品参数N及所述相邻的采样点间的弦长L,计算所述相邻的采样点间对应的弧段的弦高D。
所述弦高D=N×L,即相邻的采样点间对应的弧段的弦高D为样品参数N与相邻的采样点间的弦长L的乘积。
S305,计算由相邻的三个采样点依次相连而成的内角的角度。
所述三个采样点依次为起点、中点、终点。
根据步骤S301组合成的基于线段的多边形,依照测量顺序计算基于线段的多边形内相邻的三个采样点依次相连而成的内角的角度。
例如(参见图5),相邻的三个采样点依次相连而成的内角共有5个。分别为:采样点A、B、C组成的∠ABC,此时A为起点,B为中点,C为终点;采样点B、C、D组成的∠BCD,此时B为起点,C为中点,D为终点;采样点C、D、E组成的∠CDE,此时C为起点,D为中点,E为终点;采样点D、E、A组成的∠DEA,此时D为起点,E为中点,A为终点;采样点E、A、B组成的∠EAB,此时E为起点,A为中点,B为终点。处理器根据各采样点的X轴坐标及Y轴坐标计算上述各角的内角角度。
S306,判断所述内角的角度是否大于180°。
判断为是时,则所述中点与终点间对应的弧段为凹弧。
判断为否时,则所述中点与终点间对应的弧段为凸弧。
例如(参见图5),采样点A、B、C组成的∠ABC小于180°,此时A为起点,B为中点,C为终点,则点B与点C间对应的弧段为凸弧,同理可得,点D与点E间对应的弧段为凸弧,点E与点A间对应的弧段为凸弧,点A与点B间对应的弧段为凸弧。另外,采样点B、C、D组成的∠BCD大于180°,此时B为起点,C为中点,D为终点,则点C与点D间对应的弧段为凹弧。
S307,根据所述弦高D、弦长L及弧段的凹凸性生成所述相邻的采样点间的弧段。
更佳地,处理器获取无棱镜全站仪实时采集的采样点信息后,处理器存储所述采样点信息,实现了采样点信息的实时存储。
下面结合具体的实施例,进一步说明本方法。
实施例1
根据树冠总体生长态势,选取目标树木的树冠在各个方向的采样点,分别为A、B、C、D、E,共5点,利用无棱镜全站仪逆时针测量采样点的三维坐标,采样点的三维坐标数据如表1所示。
表1
名称 | X轴坐标 | Y轴坐标 | Z轴坐标 |
A | 1 | 4 | 10.4 |
B | 3.732050808 | 3.267949192 | 10.4 |
C | 2.732050808 | 5 | 10.4 |
D | 4.464101615 | 6 | 10.4 |
E | 1 | 6 | 10.3 |
处理器获取无棱镜全站仪实时采集的采样点坐标,并按照无棱镜全站仪的测量顺序依次提取采样点A、B、C、D、E的坐标中的X轴坐标及Y轴坐标,根据X轴坐标及Y轴坐标将所述采样点绘于同一图层内。
如图4所示,以线段形式依次连接采样点A与B,B与C,C与D,D与E,E与A。
获取预设的样品参数N=1/4.15。根据X轴坐标及Y轴坐标依次计算采样点A与B,B与C,C与D,D与E,E与A间的弦长,同时,根据公式D=N×L,代入数据,可知各采样点间弦长、弦高如表2所示。
表2
名称 | 弦长L | 弦高D=N×L |
弦AB | 2.828427125 | 0.681548705 |
弦BC | 2 | 0.481927711 |
弦CD | 2 | 0.481927711 |
弦DE | 3.464101615 | 0.834723281 |
弦EA | 2 | 0.481927711 |
然后,根据各采样点的X轴坐标及Y轴坐标,计算相邻的三个采样点依次相连而成的内角的角度,并判断出构成内角的中点与终点所对应的弧段凹凸性,如表3所示。
表3
名称 | 角度 | 弧段凹凸性 |
∠ABC | 45° | 点B与点C间对应的弧段为凸弧 |
∠BCD | 270° | 点C与点D间对应的弧段为凹弧 |
∠CDE | 30° | 点D与点E间对应的弧段为凸弧 |
∠DEA | 90° | 点E与点A间对应的弧段为凸弧 |
∠EAB | 105° | 点A与点B间对应的弧段为凸弧 |
此后,根据弦高D、弦长L及弧段的凹凸性生成所述相邻的采样点间的弧段,形成基于弧段的多边形(参见图5)。
计算基于弧段的多边形面积,采用基于扫描线方法计算面积,先将所有的弦和直线边组成的多边形划分为梯形EFCD、梯形FCGA、三角形AGB、弧(A,B)(B,C)(D,E)(E,A)围成的面积为凸出于多边形的弧形,以及弧(C,D)围成的面积凹入多边形的弧形。此时,结合采样点的X轴坐标及Y轴坐标计算出基于线段的多边形总面积(即梯形EFCD、梯形FCGA、三角形AGB的总面积)为5.464101615,而弧(A,B)所对弧形面积为2.557,弧(B,C)(C,D)(E,A)所对弧形面积均为1.278,而弧(E,D)所对弧形面积为3.835,
根据弧所对应的凹凸关系,对应的总面积减去或加上圆弧围成的面积可等基于弧段的多边形面积为13.134。
实施例2
根据树冠总体生长态势,选取目标树木的树冠在各个方向的采样点,分别为g030111、g030112、g030113、g030114、g030115、g030116、g030117,共7点,利用无棱镜全站仪顺时针测量采样点的三维坐标,采样点的三维坐标数据如表4所示。
表4
名称 | X轴坐标 | Y轴坐标 | Z轴坐标 |
g030111 | 49083.932 | 19649.779 | 10.59 |
g030112 | 49081.282 | 19648.144 | 10.532 |
g030113 | 49078.449 | 19649.429 | 10.513 |
g030114 | 49079.183 | 19645.928 | 10.537 |
g030115 | 49082.082 | 19645.26 | 10.427 |
g030116 | 49083.306 | 19646.534 | 10.46 |
g030117 | 49085.179 | 19647.824 | 10.48 |
处理器获取无棱镜全站仪实时采集的采样点坐标,并按照无棱镜全站仪的测量顺序依次提取采样点g030111、g030112、g030113、g030114、g030115、g030116、g030117的坐标中的X轴坐标及Y轴坐标,根据X轴坐标及Y轴坐标将所述采样点绘于同一图层内。
如图6所示,以线段形式依次连接采样点g030111与g030112,g030112与g030113,g030113与g030114,g030114与g030115,g030115与g030116,g030116与g030117,g030117与g030111。
获取预设的样品参数N=1/4.15。根据X轴坐标及Y轴坐标依次计算各采样点间的弦长,同时,根据公式D=N×L,计算弦高。
然后,代入各采样点的X轴坐标及Y轴坐标,计算相邻的三个采样点依次相连而成的内角的角度,并判断出构成内角的中点与终点所对应的弧段凹凸性。
此后,根据弦高D、弦长L及弧段的凹凸性生成所述相邻的采样点间的弧段,形成基于弧段的多边形(参见图7)。
最后,计算得出基于弧段的多边形的面积为21.4220914335。
实施例3
根据树冠总体生长态势,选取目标树木的树冠在各个方向的采样点,分别为g030211、g030212、g030213、g030214、g030215、g030216,共6点,利用无棱镜全站仪逆时针测量采样点的三维坐标,采样点的三维坐标数据如表5所示。
表5
名称 | X轴坐标 | Y轴坐标 | Z轴坐标 |
g030211 | 49079.167 | 19657.373 | 10.408 |
g030212 | 49081.621 | 19657.199 | 10.42 |
g030213 | 49082.56 | 19657.954 | 10.473 |
g030214 | 49081.524 | 19659.795 | 10.44 |
g030215 | 49080.704 | 19657.884 | 10.314 |
g030216 | 49079.041 | 19658.288 | 10.357 |
处理器获取无棱镜全站仪实时采集的采样点坐标,并按照无棱镜全站仪的测量顺序依次提取采样点g030211、g030212、g030213、g030214、g030215、g030216的坐标中的X轴坐标及Y轴坐标,根据X轴坐标及Y轴坐标将所述采样点绘于同一图层内。
如图8所示,以线段形式依次连接采样点g030211与g030212,g030212与g030213,g030213与g030214,g030214与g030215,g030215与g030216,g030216与g030211。
获取预设的样品参数N=1/4.15。根据X轴坐标及Y轴坐标依次计算各采样点间的弦长,同时,根据公式D=N×L,计算弦高。
然后,代入各采样点的X轴坐标及Y轴坐标,计算相邻的三个采样点依次相连而成的内角的角度,并判断出构成内角的中点与终点所对应的弧段凹凸性。
此后,根据弦高D、弦长L及弧段的凹凸性生成所述相邻的采样点间的弧段,形成基于弧段的多边形(参见图9)。
最后,计算得出基于弧段的多边形的面积为6.36798183354。
图10是本发明一种树木测量系统3的结构示意图,包括:
无棱镜全站仪1,用于实时采集目标树木的采样点信息。
所述采样点信息包括最高点高程Ht、最低点高程Hb、采样点坐标,所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标。
处理器2,用于根据所述无棱镜全站仪1实时采集的采样点信息计算树高及树冠投影面积。
需要说明的是,无棱镜全站仪1的数量可以为一个或多个,处理器2能海量地处理无棱镜全站仪1采集的采样点信息,实现树高及树冠投影面积的实时计算。
图11是本发明一种树木测量系统3中处理器2的结构示意图,包括:
获取单元21,用户获取所述无棱镜全站仪1实时采集的采样点信息。
与所述获取单元21相连的树高计算单元22,用于根据所述无棱镜全站仪1实时采集的采样点信息计算树高。
需要说明的是,树高计算单元22利用公式H=Ht-Hb,计算树高。利用无棱镜全站仪1测量时,需要考虑风速、植物通视程度、枝条粗细等因素。当枝条受风速影响小而且通视情况良好时,可直接测量植物顶端枝条以确定最高点高程Ht;否则,应测量与该枝条处于同一水平面上的叶子,且应选取枝条摆动弧度上限位置为测量点以确定最高点高程Ht。
计算树高时,仅需选取目标树木顶端枝条或者与该枝条处于同一水平面上的树叶为参考点,并利用无棱镜全站仪1测量最高点高程Ht。同时,选取与目标树木根部处于同一水平面的有效测区为参考点,并利用无棱镜全站仪1测量最低点高程Hb。无棱镜全站仪1将测量所得的最高点高程Ht及最低点高程Hb上传至处理器2。获取单元21获取最高点高程Ht及最低点高程Hb,并由树高计算单元22依据公式H=Ht-Hb计算树高H。
与所述获取单元21相连的树冠投影面积计算单元23,用于根据所述无棱镜全站仪1实时采集的采样点信息计算树冠投影面积。
测量树冠投影面积时,根据树冠总体生长态势,初步选取树冠在各个方向的最大伸展处作为采样点,再根据树冠投影形状和凹凸点分布情况最终确定采样点,利用无棱镜全站仪1顺时针或逆时针测量采样点的三维坐标。获取单元21获取无棱镜全站仪1上传的采样点三维坐标,并由树冠投影面积计算单元23根据采样点坐标中的X轴坐标及Y轴坐标,描绘采样点并根据采样点模拟出与树冠投影相似的不规则多边形,此时,计算出的多边形面积可近似为树冠投影的面积,大大提高了测量的精度。
图12是本发明一种树木测量系统3中处理器2的另一结构示意图,与图7不同的是,处理器2还包括:
与所述获取单元21相连,用于存储所述采样点信息的存储单元24。
存储单元24实现了采样点信息的实时存储,计算时,树高计算单元22及树冠投影面积计算单元23可随时提取存储单元24内的采样点信息,对树高、树冠投影面积进行计算。
图13是本发明一种树木测量系统3的处理器2中树冠投影面积计算单元23的结构示意图,包括:
线段多边形单元231,提取所述采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形。
获取单元21获取无棱镜全站仪1实时采集的采样点坐标后,线段多边形单元231依次提取采样点坐标中的X轴坐标及Y轴坐标,并根据X轴坐标及Y轴坐标,将采样点描绘于同一图层,并连接成线段多边形。
与所述线段多边形单元231相连的弧段多边形单元232,根据所述基于线段的多边形生成基于弧段的多边形。
弧段多边形单元232在基于线段的多边形的基础上,将基于线段的多边形转换为具有凹弧、凸弧的基于弧段的多边形,使模拟的多边形更接近树冠投影的形状,大大提高了测量的精度。
与所述弧段多边形单元232相连的计算单元233,用于计算所述基于弧段的多边形的面积。
计算基于弧段的多边形的面积时,可利用积分原理对弧段多边形进行面积计算,即通过计算上下曲线/折线的方程,再求上下曲线积分得到面积。
也可以使用多边形像素映射方法计算基于弧段的多边形面积,即先计算多边形内的像素数量,再按比例换算成面积。
还可以采用基于扫描线方法的面积算法计算基于弧段的多边形,即通过从每一个采样点引出的横向扫描线,把整个多边形区域划分为若干个三角形、梯形(只要存在两条对边平行的四边形就称为梯形)、或者是弧形(在存在圆弧的情况下,先将圆弧的弦连起来作为一个弧形单独考虑),从而分别求算各个分块的面积,各个分块累加即可得到整个基于弧段的多边形的面积。
图14是图13中弧段多边形单元232的结构示意图,包括:
参数获取子单元2321,用于获取用户预设的样品参数N。优选地,所述样品参数N的值为1/4.15。
需要说明的是,样品参数N为相邻样品点的弦高与弦长之比的平均值。在太阳直射点上,工作人员进行实地测量,选定样品树木,并在地面上标出样品树木的树冠投影,分别选取树冠投影上的两个相邻的样品点,沿树冠投影勾画出两个相邻的样品点间的弧段。用皮尺测量出两样品点的距离(即弦长)以及两样品点间连线的垂直平分线的长度(即弦高),此时,将两样品点间连线的垂直平分线的长度与两样品点的距离之比称为固定比。经上千次的反复测量,分别选取不同种类的样品树木及不同的样品点,计算可得出上千个固定比。此时,计算各固定比的平均值即称为样品参数N。
弦长计算子单元2322,用于计算所述相邻的采样点间的弦长L。
弦长计算子单元2322根据采样点的X轴坐标及Y轴坐标计算相邻的采样点间的弦长L。
与所述参数获取子单元2321及弦长计算子单元2322分别相连的弦高计算子单元2323,用于根据所述样品参数N及所述相邻的采样点间的弦长L,计算所述相邻的采样点间对应的弧段的弦高D。
所述弦高D=N×L。
例如,相连的采用点A、B,其中采样点A的坐标为(1,3),采样点B的坐标为(4,7),弦长计算子单元2322根据采样点A及采样点B的X轴坐标及Y轴坐标可计算相邻的采样点间的弦长L为5。相应地,弦高计算子单元2323利用公式D=N×L,可计算出采样点A及采样点B间对应的弧段的弦高D为5/4.15。
内角计算子单元2324,用于计算由相邻的三个采样点依次相连而成的内角的角度。
所述三个采样点依次为起点、中点、终点。如图5所示,相邻的三个采样点依次相连而成的内角共有5个。分别为:采样点A、B、C组成的∠ABC,此时A为起点,B为中点,C为终点;采样点B、C、D组成的∠BCD,此时B为起点,C为中点,D为终点;采样点C、D、E组成的∠CDE,此时C为起点,D为中点,E为终点;采样点D、E、A组成的∠DEA,此时D为起点,E为中点,A为终点;采样点E、A、B组成的∠EAB,此时E为起点,A为中点,B为终点。
内角计算子单元2324根据采样点的X轴坐标及Y轴坐标计算由相邻的三个采样点依次相连而成的内角的角度。
例如,相连的采用点A、B、C组成的∠ABC,其中采样点A的坐标为(1,2),采样点B的坐标为(1,1),采样点C的坐标为(1+,1),此时,内角计算子单元2324可直接计算出∠ABC的内角为90°。
与所述内角计算子单元2324相连的判断子单元2325,用于判断所述内角的角度是否大于180°并做相应处理。
判断子单元2325根据内角计算子单元2324的计算结果判断所述内角的角度是否大于180°,当判断出所述内角的角度大于180°时,则所述中点与终点间对应的弧段为凹弧;当判断出所述内角的角度小于180°时,则所述中点与终点间对应的弧段为凸弧。
例如(参见图5),判断子单元2325判断出采样点A、B、C组成的∠ABC小于180°,此时A为起点,B为中点,C为终点,则点B与点C间对应的弧段为凸弧,同理可得,点D与点E间对应的弧段为凸弧,点E与点A间对应的弧段为凸弧,点A与点B间对应的弧段为凸弧。另外,判断子单元2325判断出采样点B、C、D组成的∠BCD大于180°,此时B为起点,C为中点,D为终点,则点C与点D间对应的弧段为凹弧。
与所述弦高计算子单元2323及判断子单元2325分别相连的弧段生成子单元2326,用于根据所述弦高D、弦长L及判断子单元2325的判断结果生成所述相邻的采样点间的弧段。
由上可知,无棱镜全站仪1根据操作者的需求实时测量所需的采样点信息,进行树高计算时,无棱镜全站仪1在操作者的控制下采集高点高程Ht、最低点高程Hb并上传至处理器2,处理器2的获取单元21获取相应的高点高程Ht及最低点高程Hb,利用树高计算单元22,结合公式H=Ht-Hb,计算树高。
另外,进行树冠投影面积计算时,无棱镜全站仪1在操作者的控制下按照逆时针或顺时针的采集采样点坐标,并上传至处理器2。处理器2的获取单元21按顺序依次获取相应的采样点坐标,线段多边形单元231依次提取所述采样点坐标中的X轴坐标及Y轴坐标,并根据X轴坐标及Y轴坐标,描绘采样点并连接成线段多边形。相应地,参数获取子单元2321取用户预设的样品参数N,弦长计算子单元2322根据采样点的X轴坐标及Y轴坐标计算相邻的采样点间的弦长L,弦高计算子单元2323根据公式D=N×L,计算所述相邻的采样点间对应的弧段的弦高D,内角计算子单元2324计算由相邻的三个采样点依次相连而成的内角的角度,并由判断子单元2325判断内角的角度是否大于180°,当判断出内角的角度大于180°时,则三个采样点中的中点与终点间对应的弧段为凹弧;当判断出所述内角的角度小于180°时,则三个采样点中的中点与终点间对应的弧段为凸弧。然后,弧段生成子单元2326根据弦高D、弦长L及判断子单元2325的判断结果生成基于弧段的多边形。最后,计算单元233计算所述基于弧段的多边形的面积。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (4)
1.一种树木测量方法,其特征在于,包括:
处理器获取无棱镜全站仪实时采集的采样点信息,所述采样点信息包括最高点高程Ht、最低点高程Hb、采样点坐标,所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标;
所述处理器根据所述最高点高程Ht及最低点高程Hb计算树高H,所述树高H=Ht-Hb;
所述处理器根据所述采样点坐标计算树冠投影面积;
所述处理器根据采样点坐标计算树冠投影面积的步骤包括:提取所述采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形;根据所述基于线段的多边形生成基于弧段的多边形;计算所述基于弧段的多边形的面积;
所述提取采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形的步骤包括:提取所述采样点坐标中的X轴坐标及Y轴坐标,根据所述X轴坐标及Y轴坐标将所述采样点绘于同一图层内;相邻的采样点两两相连,组合成基于线段的多边形;
所述根据基于线段的多边形生成基于弧段的多边形的步骤包括:获取预设的样品参数N;计算所述相邻的采样点间的弦长L;根据所述样品参数N及所述相邻的采样点间的弦长L,计算所述相邻的采样点间对应的弧段的弦高D,所述弦高D=N×L;计算由相邻的三个采样点依次相连而成的内角的角度,所述三个采样点依次为起点、中点、终点;判断所述内角的角度是否大于180°,判断为是时,则所述中点与终点间对应的弧段为凹弧,判断为否时,则所述中点与终点间对应的弧段为凸弧;根据所述弦高D、弦长L及弧段的凹凸性生成所述相邻的采样点间的弧段。
2.如权利要求1所述的测量方法,其特征在于,所述处理器获取无棱镜全站仪实时采集的采样点信息的步骤之后还包括存储所述采样点信息。
3.一种树木测量系统,其特征在于,包括:
处理器,用于根据无棱镜全站仪实时采集的采样点信息计算树高及树冠投影面积;
所述无棱镜全站仪可实时采集目标树木的采样点信息,所述采样点信息包括最高点高程Ht、最低点高程Hb、采样点坐标,所述采样点坐标包括X轴坐标、Y轴坐标、Z轴坐标;
所述处理器包括:获取单元,用户获取所述无棱镜全站仪实时采集的采样点信息;与所述获取单元相连的树高计算单元,用于根据所述无棱镜全站仪实时采集的采样点信息计算树高;与所述获取单元相连的树冠投影面积计算单元,用于根据所述无棱镜全站仪实时采集的采样点信息计算树冠投影面积;
所述树冠投影面积计算单元包括:线段多边形单元,提取所述采样点坐标中的X轴坐标及Y轴坐标,生成基于线段的多边形;与所述线段多边形单元相连的弧段多边形单元,根据所述基于线段的多边形生成基于弧段的多边形;与所述弧段多边形单元相连的计算单元,用于计算所述基于弧段的多边形的面积;
所述弧段多边形单元包括:参数获取子单元,用于获取预设的样品参数N;弦长计算子单元,用于计算相邻的采样点间的弦长L;与所述参数获取子单元及弦长计算子单元分别相连的弦高计算子单元,用于根据所述样品参数N及所述相邻的采样点间的弦长L,计算所述相邻的采样点间对应的弧段的弦高D,所述弦高D=N×L;内角计算子单元,用于计算由相邻的三个采样点依次相连而成的内角的角度;与所述内角计算子单元相连的判断子单元,用于判断所述内角的角度是否大于180°并做相应处理;与所述弦高计算子单元及判断子单元分别相连的弧段生成子单元,用于根据所述弦高D、弦长L及判断子单元的判断结果生成所述相邻的采样点间的弧段。
4.如权利要求3所述的树木测量系统,其特征在于,所述处理器还包括与所述获取单元相连,用于存储所述采样点信息的存储单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210249851.0A CN103175491B (zh) | 2012-07-19 | 2012-07-19 | 一种树木测量方法及树木测量系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210249851.0A CN103175491B (zh) | 2012-07-19 | 2012-07-19 | 一种树木测量方法及树木测量系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103175491A CN103175491A (zh) | 2013-06-26 |
CN103175491B true CN103175491B (zh) | 2016-01-20 |
Family
ID=48635524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210249851.0A Active CN103175491B (zh) | 2012-07-19 | 2012-07-19 | 一种树木测量方法及树木测量系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103175491B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105571561A (zh) * | 2015-12-15 | 2016-05-11 | 浙江农林大学 | 立木信息采集方法和装置 |
CN105513078A (zh) * | 2015-12-15 | 2016-04-20 | 浙江农林大学 | 基于图像的立木信息采集方法和装置 |
CN109099856A (zh) * | 2018-07-12 | 2018-12-28 | 河北农业大学 | 一种基于方位角与距离的树冠投影测量方法及系统 |
CN112461185A (zh) * | 2020-12-03 | 2021-03-09 | 深圳市测友科技有限公司 | 一种不规则物体的测量装置和测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2192877Y (zh) * | 1993-09-17 | 1995-03-22 | 赵世凯 | 不定距离测高器 |
JP2001133255A (ja) * | 1999-11-04 | 2001-05-18 | Hitachi Zosen Corp | 建物高さ測量方法およびその装置 |
CN101021416A (zh) * | 2007-01-18 | 2007-08-22 | 北京林业大学 | 全站仪量测树高方法 |
-
2012
- 2012-07-19 CN CN201210249851.0A patent/CN103175491B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2192877Y (zh) * | 1993-09-17 | 1995-03-22 | 赵世凯 | 不定距离测高器 |
JP2001133255A (ja) * | 1999-11-04 | 2001-05-18 | Hitachi Zosen Corp | 建物高さ測量方法およびその装置 |
CN101021416A (zh) * | 2007-01-18 | 2007-08-22 | 北京林业大学 | 全站仪量测树高方法 |
Non-Patent Citations (1)
Title |
---|
无棱镜全站仪测量树高及树冠的方法研究;谢鸿宇等;《中南林业科技大学学报》;20111130;第31卷(第11期);2-3节,附图4-7 * |
Also Published As
Publication number | Publication date |
---|---|
CN103175491A (zh) | 2013-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107063197B (zh) | 一种基于空间信息技术的水库特征曲线提取方法 | |
CN113177744B (zh) | 一种城市绿地系统碳汇量估算方法及系统 | |
CN103389136B (zh) | 基于三维激光扫描技术的外浮顶立式金属罐容积测量方法 | |
CN103175491B (zh) | 一种树木测量方法及树木测量系统 | |
CN102252612B (zh) | 一种基于高程等值线法量测树冠体积的方法 | |
CN103616014B (zh) | 测量和计算森林蓄积参数的方法 | |
CN104331927B (zh) | 一种储层数字露头表层模拟的方法和系统 | |
CN105066877A (zh) | 基于智能终端镜头的树木测量方法 | |
CN101021416A (zh) | 全站仪量测树高方法 | |
WO2022062305A1 (zh) | 一种基于激光点云的线路塔基地形自动提取与成图方法 | |
CN110440762B (zh) | 一种多旋翼无人机山区航测成图的网格化像控点布设方法 | |
CN107833281A (zh) | 一种快速建立地质临空面gocad三维地质模型的方法 | |
CN101158576A (zh) | 内置程序式电子经纬仪角规自动立体测树技术 | |
CN105571561A (zh) | 立木信息采集方法和装置 | |
CN110207676A (zh) | 一种田沟塘参数的获取方法及装置 | |
CN105806266A (zh) | 基于激光扫描数据的树木冠层叶面积计算方法 | |
CN105806269B (zh) | 一种便携式非接触的树木参数测量装置与测量方法 | |
CN105466391B (zh) | 利用数字高程模型与外业数据修正生成塔基断面的方法 | |
Fan et al. | New method for forest resource data collection based on smartphone fusion with multiple sensors | |
CN107255446B (zh) | 一种矮化密植果树冠层三维地图构建系统和方法 | |
CN103185557A (zh) | 一种树木三维绿量测量方法及系统 | |
CN102207363A (zh) | 一种树高测定三角尺及其使用方法 | |
CN104898106A (zh) | 一种面向复杂地形机载激光雷达数据的地面点提取方法 | |
CN115984721A (zh) | 基于倾斜摄影和图像识别技术实现乡村风貌管理的方法 | |
CN102692190A (zh) | 一种基于数字高程模型的树木体积计测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |