CN102914501B - 一种利用激光点云计算三维森林冠层消光系数的方法 - Google Patents
一种利用激光点云计算三维森林冠层消光系数的方法 Download PDFInfo
- Publication number
- CN102914501B CN102914501B CN201210260635.6A CN201210260635A CN102914501B CN 102914501 B CN102914501 B CN 102914501B CN 201210260635 A CN201210260635 A CN 201210260635A CN 102914501 B CN102914501 B CN 102914501B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- line
- canopy
- projection
- blade
- 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.)
- Expired - Fee Related
Links
- 230000008033 biological extinction Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000009826 distribution Methods 0.000 claims abstract description 40
- 238000005070 sampling Methods 0.000 claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000008859 change Effects 0.000 claims abstract description 8
- 238000012935 Averaging Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 5
- 244000062793 Sorghum vulgare Species 0.000 claims description 3
- 235000019713 millet Nutrition 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000004088 simulation Methods 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 11
- 230000005855 radiation Effects 0.000 abstract description 10
- 230000006378 damage Effects 0.000 abstract description 2
- 238000007781 pre-processing Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 19
- 239000010410 layer Substances 0.000 description 17
- 238000000605 extraction Methods 0.000 description 7
- 241000196324 Embryophyta Species 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 241001269238 Data Species 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 241000208140 Acer Species 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 239000002356 single layer Substances 0.000 description 2
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 1
- 239000002028 Biomass Substances 0.000 description 1
- 235000014466 Douglas bleu Nutrition 0.000 description 1
- 240000001416 Pseudotsuga menziesii Species 0.000 description 1
- 235000005386 Pseudotsuga menziesii var menziesii Nutrition 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004173 biogeochemical cycle Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 108700002783 roundabout Proteins 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
Abstract
本发明提供了一种改进的利用计算几何算法来计算三维森林冠层任意给定入射光线下的消光系数的方法,属于森林冠层结构参数获取方法的研究领域。其步骤为:植被冠层的三维激光点云数据的获取及预处理;点云数据的三维网格化;基于体元数据结构的点云切片算法;点云切片的线采样分析;第k层切片的平均投影系数的计算;各层切片任意给定入射光线下的消光系数及整个冠层消光系数的计算。本发明与传统观测手段相比,工作量小,无需接触式观测,不破坏冠层结构和辐射特性,具有客观高效精确的特点;开发了从激光雷达数据中提取三维结构和生物物理多样性信息的方法,将叶片的水平和垂直分布变化规律特征化。
Description
一、技术领域
本发明涉及一种利用地面三维激光扫描仪获取的点云数据计算森林冠层消光系数的方法,具体地说,是指一种改进的利用计算几何算法来计算三维森林冠层任意给定入射光线下的消光系数的方法(流程如附图1所示)。
二、背景技术
森林冠层结构在土壤-冠层-大气的相互作用中占有重要的地位,并与生物地球化学循环中的物质能量交换密切相关。消光系数定义为单位叶面积在垂直于光线方向的平面上的平均投影面积,是定量的理解森林冠层结构和辐射传输的重要因子。他是基于叶倾角的方向性分布函数(包括倾角和方位角分布函数),对于任意给定角度的入射光线来计算的。
间接、准确地获取植被冠层的叶倾角分布密度函数一直以来是一个具有挑战性的工作,传统的方法主要是利用逐叶片直接观测的方法来获取给定单株树的叶倾角分布密度函数。比如:Norman等人利用量角器和指南针等设备制作了一种较为简单方便的测量叶片角度和空间分布的仪器(附图2左);Lang等人利用机械臂的灵活性和巧妙的几何传递关系制作了一种可以自由移动和伸缩的用来测量叶片空间和角度分布的仪器(附图2右)。
但由于其工作量大,主观性强,而且对所测量的目标树木的高度有限制,所以不能够广泛地应用于实际研究中。更为重要的一点是,直接观测往往会影响甚至破坏植被冠层的结构,进而破坏冠层内部和以下的辐射分布情况。而且就森林而言,实际中很难利用直接观测仪器来进行逐叶片测量。因此,需要一种新的快速、准确、间接的方法来获取冠层叶片的角度和空间分布情况。近年来,也有人开始尝试进行间接地获取叶片倾角分布。例如,2007年Shilbayama和Watanabe发表在期刊《Plant Production Science》第10卷的“Estimating the mean leafinclination angle of wheat canopies using reflected polarized light”一文中,提出利用反射极化光来进行平均叶倾角的估算方法。.
目前,还有一种较为流行的方法是利用数学几何模型来近似地逼近真实树木的叶片实际分布,其中被广泛接受的是Campbell在20世纪80年代提出的椭球体儿何模型。其核心思想是椭球体表面各切面倾角的不同分布来近似表达实际叶片的角度分布,并通过调整椭球体的长、短轴来调节不同倾角分布密度函数的变化。该方法可用以下公式表达:
其中
当χ=1时,该分布将成为球面分布,此时Λ=2,a、b分别是椭球体的半短、半长轴。但这种方法只是在理论上近似的逼近叶片的实际分布,不足以描述变化多端的真实情况。
地面激光雷达扫描仪的出现使得我们可以直接从三维角度获取高精度(毫米级)的森林冠层结构信息,进而计算冠层的消光系数。目前国内关于激光雷达在森林中的应用研究,特别是针对三维冠层结构的研究还处在起步和探索阶段。已有的研究大部分集中在利用地面激光雷达进行单株立木或森林样地水平的基本每木检尺的参数(包括胸径、树高和生物量等)提取方面。例如:2011年马利群等人在《世界林业研究》第1卷的“激光雷达在森林垂直结构参数估算中的应用”一文中介绍了激光雷达在提取森林垂直结构参数方面的应用;在激光雷达的森林应用方面,生成高精度的数字地面模型是其中关键的一步,周淑芳等人在2007年《遥感技术与应用》第22卷的“基于机载激光雷达数据的DEM获取及应用”一文中探讨了利用机载激光雷达***生成数字地面模型的方法和精确评价。除此之外,庞勇等人在2006-2009年间对利用激光雷达***提取森林单木和样方水平的树高估算进行了大量较深入的研究。但目前国内对利用激光雷达所获取的三维点云数据提取冠层消光系数的研究理论和技术还不够深入,需要进一步加强。
三、发明内容
本发明的目的是:
提供一套直接利用地面激光雷达***生成的三维点云数据,结合计算几何的方法,从三维角度来计算任意给定入射光的单株树木或森林样方消光系数的算法。并提供非接触式的快速准确观测森林冠层结构特征参数的方法。
本发明的原理如下:
利用较新的遥感技术手段(地面三维激光雷达扫描***),结合计算几何的技术手段,直接从三维的角度获取给定方向入射光下的单株树木和森林样方尺度叶片的消光系数,进而计算植被冠层的消光系数。首先建立三维空间的数学模型来精确定位单个叶片的位置,在此基础上采用计算几何方法,对三维点云数据进行网格化,利用点云切片算法对不同高度的点云进行分类。然后进行线采样分析,获取叶片截获光线的相关信息,进而计算每个切片层的投影系数,即消光系数,从而得到整个冠层的消光系数。完成给定入射光线下的消光系数及冠层整体的平均投影系数的计算。
本发明的技术方案主要包括以下步骤:
(1)首先利用地面激光雷达扫描***,获取植被冠层的三维点云数据。其中包含了扫描目标点的空间几何和激光束弹回的能量信息,三维坐标直接提供了任意一点的空间位置坐标信息,这也是传统光学理论中的数学模型的基础。获取的点云数据首先进行图像拼接,并手动去除地面点云。
(2)三维网格化。首先在点云区域定义一个X,Y,Z三个坐标轴的笛卡尔坐标系,将所有点云数据包含其中。并将点云数据分为有限的小区域,建立以体元(voxel)为基础的数据结构,这个过程叫做三维网格化。每个体元由长(l)宽(w)、高(h)三个参数决定其大小,所有的点云数据被分为m×n×p个体元,其中,m=(Xmax-Xmin)/w,n=(Ymax-Ymin)/l,P=(Zmax-Zmin)/h。如果,l=w=h,则体元为一个立方体。以这些点云数据的三维中心点,即[(Xmax-Xmin)/2,(Ymax-Ymin)/2,(Zmax-Zmin)/2]为新的原点,建立笛卡尔坐标系。坐标轴Z为树干的生长方向,垂直于水平面,X、Y轴位于与Z轴垂直的平面上,Z轴为纵向的轴,X、Y轴为横向的轴。如附图3(a,b,c)所示。
(3)点云切片算法。三维网格化过程结束后,点云数据被分为纵横方向的众多数据层,即切片,点云数据可视为纵横方向切片的叠加。如附图3(d)所示水平切片。所有体元设为(i,j,k)(i=1,2,...,m;j=1,2,...n;k=1,2,...,p)。当k=1,i和j为给定区域内的任意值,则可表示为第一层或第一横切面的所有体元。
为了近似模拟入射光线与森林冠层的相对位置关系,发展了另外一种旋转点云的切片方式。首先,将点云数据旋转至我们要模拟的相对位置,然后对其进行三维网格化。这样水平切片与XY轴所在的平面平行,垂直的切片与Z轴所在的平面平行。我们假设太阳直射光从天顶方向射入,与Z轴方向平行。当入射光方向改变时,可通过旋转点云数据来实现。
具体的说,三维点云数据的中心点为笛卡尔坐标系的原点,点云数据则根据点云与切片相对位置的需要围绕X、Y、Z轴旋转。例如,将点云数据围绕Y轴旋转30°,相当于点云被切片或太阳入射光以30°的倾斜角进行切片。因此,点云可在任意位置(水平0°-360°,垂直0°-90°)被切片,称作方向性切片。任何点云数据都可通过方向性切片生成薄板以提取三维冠层信息。
(4)线采样分析。叶片在切片中的空间分布通过线采样方法进行分析。点云数据每个位置的切片都是由许多样线组成的,这些样线垂直或平行于太阳入射光r(θ,β)。当研究沿r方向叶片的分布时,是考虑与r平行的样线。如果将每个体元的三维参数设置为激光雷达采样距离的1.5倍,我们可以保证每个体元最多只包含一个样点。那么,每条样线中非空的体元(N)个数代表样线与叶片相交的次数,且N服从随机分布。如果大于采样距离的1.5倍,一个体元内包含多个样点时,我们需要将样线上所有的点投影到与样线方向平行的面上,计算不同的N值,作为与叶片相交的次数。
定义Pn为样线穿过整个点云数据时,与叶片相交n次的可能性。P0表示样线穿越的都是空隙,叠加指数为0。通过求解所有样线穿过植被时N的平均值,可以求得整个点云数据中所有样线与叶片相交次数的平均值(m)和方差(σ2)。这个相交的次数并不是实际的叶片点云密度,而是投影到与样线方向垂直的平面上的投影。每条样线的方差(σ2)/均值(m)的相对变化可以反映叶片在真实空间的分布状况。单株树或林分的叶片分布可根据其相对变化(σ2/m=1,>1,和<1)相应的分为三种情况:规则分布、随机分布和聚集分布。
线采样分析是计算整个冠层和不同高度切片消光系数的关键,同时有利于点云密度、辐射通量和辐射分布状况等参数的获取。
(5)太阳入射光倾斜角(θ)固定时,计算第k层切片的平均投影系数。当只研究r(θ,β)一个方向的太阳入射光时,也就是说第k层切片的投影系数Gk(r)是指固定的r(θ,β)方向的投影系数。当不考虑入射光方位角,而倾斜角为θ时的k层平均投影系数Gk(θ),是通过方位角β在[0,2π]上的积分计算的:
(6)计算任意给定入射光线下的投影系数,即消光系数。首先考虑以固定方向r(θ,β)进行切片,在点云数据完成三维网格化之后,所有的三维点云都被包围在立体边界内,生成水平和垂直的切片及样线。如果所有的切片与方向r(θ,β)平行,那么所有的样线也平行于方向r(θ,β),且纵轴的方向即代表此方向。每个体元可以构建一个垂直于这个主轴的平面,所有的样点都可以投影到这个平面上。在三维网格框架中,这个投影面是两个垂直于样线方向的体元边界之。设,每个非空体元内包含的样点个数为n:
a.若n=1,每个二维的正方形以用户采样间距为边长,则样点代表这个正方形叶片区域的中心。若采样间距为s,则叶片区域面积为s×s。投影系数为以取样间距为边长的正方形面积(即s×s),与体元的长度(l)×宽度(w)的面积之比。
b.若n=2,可以构建一条样线来表示叶片区域,叶片投影面积为这条样线的投影长度乘以正方形的单位面积。投影系数为投影后和投影前的样线长度比。
c.若n≥3,首先检测是否所有的点都在一个面或是一条线上。如果在同一样线上,则将这些点分解到更小的尺度上(n=2)。如果在同一个面上,则构建一个三角形来表示叶片区域,叶片投影面积由三个顶点的投影坐标来计算。三个投影坐标设为(x1,y1),(x2,y2),(x3,y3),三角形的投影面积Ap即叶片投影面积,由以下矩阵计算:
d.若体元内的所有点都不共面或在一条线上,则构建一个三维的凸包,凸包表面积的一半即表示实际的叶片面积(A)。然后,将所有的点投影到与样线平行的平面上,用这些点建立二维的凸包,其面积即表示投影的叶片面积(Ap)
综上所述,投影系数Ap/A是由非空体元计算。第k层垂直于入射光的投影系数Gk(r)可由以下公式计算:
其中N为第k层切片非空体元的总个数,k为切片的层数。
求出Gk(r)后,代入公式(2),求出各切片层的平均投影系数,整个冠层的投影系数(即消光系数)为各层切片投影系数的总和。
与现有技术相比,本发明优势体现在:激光雷达技术为从三维角度获取森林冠层结构参数提供了有效的技术手段,结合计算几何方法进行间接观测,与传统观测手段相比,工作量小,无需接触式观测,不破坏冠层结构和辐射特性,具有客观高效精确的特点;开发了从激光雷达数据中提取三维结构和生物物理多样性信息的方法,将叶片的水平和垂直分布变化规律特征化。
具体有益效果如下:
本发明利用激光扫描仪获取的三维点云数据,以体元为基础计算森林冠层的消光系数,提供了一种间接的无需接触的观测森林冠层生物物理参数的方法。本发明通过旋转点云数据而保持太阳光从树冠上方沿纵轴平行方向入射,来模拟现实情况中太阳光在不同角度照射树冠的情景,既简便又高效,且不会对森林结构和辐射特性造成任何不良影响。本发明不仅可以计算冠层整体的消光系数,还可以计算树冠不同高度切片的消光系数,这有利于进一步研究光照辐射在冠层内部的截留和分布机制,以及其他生物物理参数的研究。本发明从三维空间对森林冠层机构进行观测,其关键是利用几何投影技术将三维点云数据转换为二维栅格图像,从不同方向计算冠层不同高度平面的消光系数,不仅极大提高了估算精度,而且为冠层内部的辐射传输机制的研究提供了大量信息。
实际应用表明,本发明提供了间接观测森林冠层消光系数的有效方法,无论是单株树木还是森林样方冠层的消光系数,均可直接使用本方法求得。该方法克服了传统方法需花费大量时间和大量人力物力,影响森林冠层结构特性,且误差大的缺陷。提高了植被生物物理参数估算的效率,增强了三维激光扫描技术应用的普及性和有效性。可以更好地服务于森林资源调查、植被生态遥感等资源环境研究项目。
四、附图说明
图1为本发明的流程示意图;
图2为传统的叶片角度和空间分布观测仪器;
图3为体元数据结构示意图和水平切片示意图
a.两个单独的体元,以及决定其大小的长、宽、高三个参数;
b.一个2×2×2的三维栅格空间;
c.体元结构的三维空间网格化示意图;
d.基于体元结构的点云水平切片算法示意图;
图4为一株北美花旗松的三维点云数据;
图5为三维网格化和点云切片示意图
a.针对单株北美花旗松的基于体元结构的点云水平切片算法示意图;
b.单层切片的体元三维空间网格化(5×5×5个体元)示意图;
图6为单株立木的三维点云数据在三维网络化的基础上的线采样示意图
a.“线采样”(图中红色长方体)的提取和分析;
b.线采样提取样本的横断面;
c.线采样提取样本的直方图分析。
图7为单个体元内部叶片消光系数计算方法示意图
a.单个体元内部叶片三维点云示意图;
b.叶片三维凸包结果及面表达;
c.叶片三维凸包结果及点、线表达;
d.叶片在水平面上投影点集的二维凸包计算结果及面表达。
图8为14m-16m高度的点云切片层消光系数分布特征示意图
a.本切片层点云分布示意图;
b.本切片层三维消光系数空间分布特征示意图;
c.本切片层二维点云密度空间分布特征示意图;
d.本切片层二维消光系数空间分布特征示意图。
五、具体实施方式
以下通过实例对本发明作进一步解释:
以一株北美花旗松(Douglas-fir)为研究对象(树高约25m),使用地面三维激光扫描仪LeicaScanStation 2(其参数如表1所示)和高精度GPS,在树的一个侧面进行三维点云数据的采集,采样间距为2cm。手动去除地面点云及其他噪声点云,得到单株北美花旗松的三维点云数据,如附图4所示。
表1 三维激光扫描仪Leica ScanStation 2参数
在取得花旗松的点云数据并预处理后,如附图1所示,采用计算几何技术对三维点云数据进行处理,第一步为三维网格化,第二步为基于体元的点云切片,第三步为线采样分析。
根据技术方案步骤(2)所述,对点云数据进行网格化,建立体元数据结构,每个体元大小为0.5m×0.5m×1m。如附图5(b)中所示,为单层切片的体元三维空间网格化(5×5×5个体元)示意图,图中的点云颜色代表不同高度,所示点云为附图5(a)中高度为17至18米间的点云。
对三维网格化后的点云数据进行切片,如附图5(b)为针对一株北美花旗松的基于体元结构的水平点云切片算法示意图,灰色为水平的7个切片平面,点云的不同颜色代表不同的高度。为了模拟太阳光从不同方向照射到树冠的情形,我们将太阳入射光固定在与Z轴平行的正上方,通过旋转花旗松的点云数据,来模拟太阳入射光与树冠相对角度的变化,以计算任意角度入射光下的投影系数。
根据技术方案步骤(4)所述,对三维网格化的点云数据进行线采样提取和分析。我们对花旗松点云数据进行垂直方向的线采样提取分析,如附图6(a)所示,红色长方体为“线采样”,附图6(b)为线采样提取样本的横断面,附图6(c)为线采样提取样本的直方图分析,可见14m-18m间的点云数据最多。
为更好地描述单个体元内部叶片的消光系数计算方法,下面以大叶枫树的一个叶片为例进行说明,如附图7所示,7(a)为枫树叶片的三维点云。根据技术方案步骤(6)中d所述,对叶片点云构建一个三维凸包,如附图7(b)所示,为凸包的计算结果与面表达。附图7(c)为凸包的计算结果及点、线表达。计算三维凸包的表面积,其表面积的一半即表示实际的叶片面积(A)。然后将三维点云投影到与太阳入射光垂直的平面上(即与Z轴垂直、与XY轴平行的平面),用这些点建立二维的凸包,其面积即表示投影的叶片面积(Ap),如附图7(d)所示。Ap/A即该叶片的消光系数。
根据本发明提出的算法,对花旗松点云进行分析,依照技术方案步骤(5)和(6)所述,获取不同高度切片中每个体元的点云密度空间分布和消光系数。以高度为14m-16m处的切片为例,如附图8所示:图8(a)为本切片层点云分布示意图,图8(b)为本切片层三维消光系数空间分布特征示意图,图8(c)和8(d)分别为本切片层二维点云密度和消光系数空间分布特征示意图。
从以上结果可以看出,本发明还可以提供树冠叶片密度垂直变化及其特定高度水平分布的三维信息,是其他方法难以比拟的。
Claims (2)
1.一种利用激光点云计算三维森林冠层消光系数的方法,其主要包括以下步骤:
(1)植被冠层的三维激光点云数据的获取及预处理;
(2)点云数据的三维网格化:定义一个以X,Y,Z为轴的笛卡尔坐标系,将点云数据分为有限的小区域,建立以体元为基础的数据结构,每个体元由长1、宽w、高h三个参数决定其大小;
(3)基于体元结构的点云切片算法:三维网格化过程结束后,点云数据被分为纵横方向的多个切片,点云数据为纵横方向切片的叠加;所有体元设为(i,j,k),其中i=1,2,...,m;j=1,2,...n;k=1,2,...,p;当k=1,i和j为给定区域内的任意值,则表示为第一横切面的所有体元;当太阳直射光从天顶方向射入,与Z轴方向平行;通过旋转点云,近似模拟入射光线与森林冠层的相对位置的变化,点云在水平0°-360°、垂直0°-90°被切片,称作方向性切片;
(4)线采样分析:叶片在切片中的空间分布通过线采样方法进行分析,每条样线中非空的体元N个数代表样线与叶片相交的次数,且N符合随机分布;定义Pn为与太阳入射光方向平行的样线穿过整个点云数据时,与叶片相交n次的可能性;P0表示样线穿越的都是空隙,叠加指数为0;通过求解所有样线穿过植被时N的平均值,求得整个点云数据中所有样线与叶片相交次数的平均值m和方差σ2;每条样线的方差σ2/均值m的相对变化反映叶片在真实空间的分布状况;单株树或林分的叶片分布根据其相对变化σ2/m=1、>1和<1相应的分为三种情况:规则分布、随机分布和聚集分布;
(5)太阳入射光倾斜角θ固定时,计算第k层切片的平均投影系数:第k层切片的投影系数Gk(r)是指固定的r(θ,β)方向的投影系数;当不考虑入射光方位角,而倾斜角为θ时的第k层平均投影系数Gk(θ),是通过方位角β在[0,2π]上的积分计算的:
(6)计算任意给定入射光线下的投影系数,即消光系数:将方向r(θ,β)的太阳入射光固定为与Z轴平行,对依照入射光与冠层的相对位置进行旋转的点云数据进行切片,生成水平和垂直的切片及样线;每个非空体元内包含的样点个数为n
a.若n=1,每个二维的正方形以用户采样间距为边长,则样点代表这个正方形叶片区域的中心;若采样间距为s,则叶片区域面积为s×s;投影系数为:以取样间距为边长的正方形面积即s×s与体元的长度1×宽度w的面积之比;
b.若n=2,构建一条样线来表示叶片区域,叶片投影面积为这条样线的投影长度乘以正方形的单位面积;投影系数为投影后和投影前的样线长度比;
c.若n≥3,首先检测是否所有的点都在一个面或一条线上;如果在同一样线上,则将这些点分解到更小的尺度上,即n=2;如果在同一个面上,则构建一个三角形来表示叶片区域,叶片投影面积由三个顶点的投影坐标来计算;三个投影坐标设为(x1,y1),(x2,y2),(x3,y3),三角形的投影面积Ap即叶片投影面积,由以下矩阵计算:
d.若体元内的所有点都不共面或在一条线上,则构建一个三维的凸包,凸包表面积的一半即表示实际的叶片面积A;然后,将所有的点投影到与样线平行的平面上,用这些点建立二维的凸包,其面积即表示投影的叶片面积Ap;
投影系数Ap/A是由非空体元计算;
第k层垂直于入射光的投影系数Gk(r)由以下公式计算:
其中N为第k层切片非空体元的总个数,k为切片的层数;
求出Gk(r)后,代入公式(1),求出各切片层的平均投影系数,整个冠层的投影系数,即消光系数,为各层切片投影系数的总和。
2.根据权利要求1所述的一种利用激光点云计算三维森林冠层消光系数的方法,其特征在步骤(1)中,所述三维激光点云数据是由地面三维激光扫描仪获取的森林植被冠层点云,其中包含了扫描目标点的空间几何和激光束弹回的能量信息,以及每个点的空间位置坐标信息,对获取的点云进行图像拼接,并手动去除地面点云,作为提取冠层结构信息的数据源。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210260635.6A CN102914501B (zh) | 2012-07-26 | 2012-07-26 | 一种利用激光点云计算三维森林冠层消光系数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210260635.6A CN102914501B (zh) | 2012-07-26 | 2012-07-26 | 一种利用激光点云计算三维森林冠层消光系数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102914501A CN102914501A (zh) | 2013-02-06 |
CN102914501B true CN102914501B (zh) | 2015-01-14 |
Family
ID=47612963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210260635.6A Expired - Fee Related CN102914501B (zh) | 2012-07-26 | 2012-07-26 | 一种利用激光点云计算三维森林冠层消光系数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102914501B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201303540D0 (en) * | 2013-02-27 | 2013-04-10 | Arithmetica Ltd | Image processing |
CN105371789B (zh) * | 2015-10-09 | 2019-02-19 | 南京大学 | 一种利用航空激光点云计算有效叶面积指数的方法 |
CN106383998B (zh) * | 2016-09-09 | 2019-03-12 | 厦门大学 | 一种基于地面激光雷达扫描的树木胸径自动计算方法 |
CN107643048B (zh) * | 2017-07-28 | 2020-09-29 | 北京农学院 | 基于点云数据的测树因子自动提取方法 |
CN108171761B (zh) * | 2017-12-13 | 2020-10-16 | 北京大学 | 一种基于傅里叶图变换的点云帧内编码方法及装置 |
CN109472850B (zh) * | 2018-11-26 | 2020-12-04 | 广东精鹰传媒集团股份有限公司 | 一种三维立体电流光效果的实现方法 |
CN110702166B (zh) * | 2019-09-29 | 2021-09-28 | 北京农业信息技术研究中心 | 一种用于植物冠层参数测量的装置 |
CN111006645A (zh) * | 2019-12-23 | 2020-04-14 | 青岛黄海学院 | 基于运动与结构重建的无人机测绘方法 |
CN111580131B (zh) * | 2020-04-08 | 2023-07-07 | 西安邮电大学 | 三维激光雷达智能车在高速公路上识别车辆的方法 |
CN111551530B (zh) * | 2020-04-23 | 2022-07-22 | 江苏大学 | 一种作物群体冠层叶绿素荧光三维分布信息获取方法与装置 |
CN111581710B (zh) * | 2020-05-19 | 2021-04-13 | 北京数字绿土科技有限公司 | 架空输电线路杆塔挠度自动获取方法及装置 |
CN112254820B (zh) * | 2020-10-14 | 2021-04-27 | 中国科学院空天信息创新研究院 | 一种离散森林场景热红外辐射传输模拟方法 |
CN112287859A (zh) * | 2020-11-03 | 2021-01-29 | 北京京东乾石科技有限公司 | 物体识别方法、装置和***,计算机可读存储介质 |
CN112964172B (zh) * | 2020-12-08 | 2022-08-26 | 聚时科技(上海)有限公司 | 基于结构光相机的航空叶片表面量测方法及量测设备 |
CN112946669B (zh) * | 2021-02-05 | 2023-04-25 | 国际竹藤中心 | 基于地基激光雷达的毛竹林数量识别方法及装置 |
CN114663786B (zh) * | 2022-03-21 | 2024-03-22 | 南京林业大学 | 基于点云数据和计算机图形学的林分辐射通量计算方法 |
CN117434517B (zh) * | 2023-12-21 | 2024-02-27 | 西南林业大学 | 一种基于消光系数优化的森林冠层高度估测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520307A (zh) * | 2008-02-26 | 2009-09-02 | 中国计量科学研究院 | 应用三维激光影像扫描***测定树冠体积分形维数方法 |
CN102305622A (zh) * | 2011-06-14 | 2012-01-04 | 北京林业大学 | 一种基于三维激光扫描仪的乔木三维绿量测定方法 |
CN102608620B (zh) * | 2012-03-12 | 2013-09-18 | 北京北科安地科技发展有限公司 | 一种依据反射强度与地形形态的激光扫描点云植被过滤方法 |
-
2012
- 2012-07-26 CN CN201210260635.6A patent/CN102914501B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102914501A (zh) | 2013-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102914501B (zh) | 一种利用激光点云计算三维森林冠层消光系数的方法 | |
Gastellu-Etchegorry et al. | Discrete anisotropic radiative transfer (DART 5) for modeling airborne and satellite spectroradiometer and LIDAR acquisitions of natural and urban landscapes | |
CN103413133B (zh) | 无序激光点云数据中自动提取电力线方法 | |
CN102314546B (zh) | 基于虚拟植物的植物生长生物量变化估算方法 | |
Zheng et al. | Computational-geometry-based retrieval of effective leaf area index using terrestrial laser scanning | |
Boudreault et al. | A LiDAR method of canopy structure retrieval for wind modeling of heterogeneous forests | |
Ishimoto et al. | Irregularly shaped ice aggregates in optical modeling of convectively generated ice clouds | |
CN104699952B (zh) | 一种湿地水生植被冠层brdf蒙特卡洛模型 | |
CN103983230B (zh) | 一种地面叶面积指数间接测量验证方法 | |
CN101673413B (zh) | 基于表面数据的植物枝体三维形态建模***及方法 | |
CN105371789A (zh) | 一种利用航空激光点云计算有效叶面积指数的方法 | |
CN102401898A (zh) | 一种合成孔径雷达森林遥感数据的定量化模拟方法 | |
CN105389538A (zh) | 一种基于点云半球切片估算森林叶面积指数的方法 | |
Nevalainen et al. | Estimating the rut depth by UAV photogrammetry | |
CN106909750B (zh) | 一种阔叶植被冠层反射率的计算方法 | |
CN107479065A (zh) | 一种基于激光雷达的林窗立体结构量测方法 | |
CN102034266A (zh) | 激发荧光断层成像的快速稀疏重建方法和设备 | |
CN105844057A (zh) | 基于光束和三角面片求交的激光扫描成像快速仿真方法 | |
CN114663786B (zh) | 基于点云数据和计算机图形学的林分辐射通量计算方法 | |
Liang et al. | Application and algorithm of ground-penetrating radar for plant root detection: a review | |
Kallel | FluLCVRT: Reflectance and fluorescence of leaf and canopy modeling based on Monte Carlo vector radiative transfer simulation | |
Sun et al. | Feasibility study on the estimation of the living vegetation volume of individual street trees using terrestrial laser scanning | |
Lou et al. | Growth parameter acquisition and geometric point cloud completion of lettuce | |
CN102865814A (zh) | 植物群体三维重建误差测量方法 | |
Jiang et al. | Simulation on different patterns of mobile laser scanning with extended application on solar beam illumination for forest plot |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150114 Termination date: 20210726 |
|
CF01 | Termination of patent right due to non-payment of annual fee |