CN108303043B - 多传感器信息融合的植物叶面积指数检测方法及系统 - Google Patents

Abstract

本发明公开了一种多传感器信息融合的植物叶面积指数检测方法及系统，所述方法包括：实现各传感器之间的数据同步；采用激光雷达传感器扫描植物冠层侧面，获得植物冠层切面轮廓的二维激光点云数据；采用惯性测量单元获取激光雷达传感器的姿态角，并根据姿态角对二维激光点云数据进行矫正；采用卫星导航定位系统获取激光雷达传感器的定位坐标，根据定位坐标计算两次激光扫描点之间的距离，融合定位坐标与矫正后的二维激光点云数据重构植物冠层三维激光点云图；获取叶面积指数反演计算所需的相关参数，构建叶面积指数估算模型，运算模型得到植物的叶面积指数。本发明与现有技术相比，改善了颠簸路面对叶面积指数测定影响，提高了叶面积指数测定效率。

Description

多传感器信息融合的植物叶面积指数检测方法及系统

技术领域

本发明涉及一种植物叶面积指数检测方法，尤其是一种多传感器信息融合的植物叶面积指数检测方法及系统，属于农业信息检测领域。

背景技术

传感器是获取信息的重要工具，其作用类似于人类的感知器官。多传感器信息融合(Multi-sensor Information Fusion，MSIF)技术是通过多类同构或异构传感器数据进行综合(集成或融合)获得比单一传感器更多的信息，形成比单一信源更可靠、更完全的融合信息，它突破单一传感器信息表达的局限性，避免单一传感器的信息盲区，提高了多源信息处理结果的质量，有利于对事物的判断和决策。

在使用车载激光雷达传感器检测植物叶面积指数时，由于路面不平整，很容易造成车上激光传感器的俯仰角与横滚角的变化，同时车体车在复杂地形条件下很难操控，很难保证车体的直线前进，又会造成激光传感器偏航角的变化，所以对由于地面不平整时激光传感器所测量的数据进行矫正具有重要的意义。

叶面积指数是指作物叶片的总面积与树冠投影占地面积的比值，是生态学研究的关键参数之一，具有重要的理论价值与实际价值。目前叶面积指数的测量方法大多效率低下，费时费力。例如专利号为201210004979.0的发明专利公开了一种叶面积指数快速自动测定装置及方法，其利用激光传感器测量树叶的层数，从而计算叶面积指数，这种方法在每次测量过程中需要按照每个树冠的形状调整测量时所走的路径，测量速度慢，效率低；专利号为201410746072.0的发明专利公开了一种基于激光雷达点云技术的植物叶面积指数测定方法，其利用激光雷达原始的数据点进行叶面积指数计算，没有考虑路面颠簸对检测结果的影响。

如何实现叶面积指数检测的高效率、高准确率已经成为急需解决的技术问题。

发明内容

本发明的目的是为了解决上述现有技术中叶面指数检测效率低、检测准确率低的问题，提供一种多传感器信息融合的植物叶面积指数检测方法。

本发明的另一目的在于提供一种多传感器信息融合的植物叶面积指数检测系统。

本发明的目的可以通过采取如下技术方案达到：

一种多传感器信息融合的植物叶面积指数检测方法，所述方法包括：

基于卡尔曼的滤波方法实现激光雷达传感器、惯性测量单元和卫星导航定位系统的数据同步；

采用激光雷达传感器扫描植物冠层侧面，获得植物冠层切面轮廓的二维激光点云数据；

采用惯性测量单元获取激光雷达传感器的姿态角，并根据姿态角对所述二维激光点云数据进行矫正；其中，所述姿态角包括横滚角、俯仰角和偏航角；

采用卫星导航定位系统获取激光雷达传感器的定位坐标，根据定位坐标计算两次激光扫描点之间的距离，融合定位坐标与矫正后的二维激光点云数据重构植物冠层三维激光点云图；

获取叶面积指数反演计算所需的相关参数，构建叶面积指数估算模型，运算模型得到植物的叶面积指数。

进一步的，所述二维激光点云数据由一系列扫描角度间隔为α的数据点构成，每个数据点由三部分组成，包括扫描角度θ，该角度θ对应的径向距离radialDis，探测目标点的表面反射率RSSI，第n次激光扫描采样的点云数据以坐标(θ_k,radialDis_k,RSSI_k)_n表示，其中k在1到激光雷达传感器最大扫描线数之间取值，扫描角度间隔α等于(θ_k-θ_k-1)。

进一步的，所述根据姿态角对所述二维激光点云数据进行矫正，如下式：

其中，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样时激光雷达传感器的姿态角，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样的点云数据，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据，Fun1为矫正函数，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值。

进一步的，所述根据定位坐标计算两次激光扫描点之间的距离，如下式：

d＝Fun2(Coord_n,Coord_n-1)

其中，Coord_n为激光雷达传感器第n次激光扫描点的定位坐标，Fun2为两点定位坐标求解距离d的函数。

进一步的，所述获取叶面积指数反演计算所需的相关参数，构建叶面积指数估算模型，运算模型得到植物的叶面积指数，具体包括：

选择一株样本植物，对样本植物进行取样测量，获取样本植物准确的叶面积指数LAI_real；

获取样本植物的三维激光点云图数据(Coord,(θ’_k,radialDis’_k,RSSI’_k))_n；其中，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值，Coord_n为激光雷达传感器第n次激光扫描点的定位坐标，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据；

对样本植物构建叶面积指数估算模型，利用样本植物的叶面积指数LAI_real计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，如下式：

其中，LWA为植物的叶墙面积，P_hit为击中点的比例，P_through为击穿点的比例；

按照上述步骤，再选择一株样本植物，计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，取两次计算结果的平均值作为最终的计算值；

利用最终计算得到的叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，代入叶面积指数估算模型进行反演计算，得到单株植物的叶面积指数LAI，如下式：

进一步的，所述击中点是指小于植物行距的激光扫描测距点，所述击穿点是指大于植物行距的激光扫描测距点。

本发明的另一目的可以通过采取如下技术方案达到：

一种多传感器信息融合的植物叶面积指数检测系统，所述检测系统包括：

数据同步单元，用于基于卡尔曼的滤波方法实现激光雷达传感器、惯性测量单元和卫星导航定位系统之间的数据同步；

激光雷达传感器，用于扫描植物冠层侧面，获得植物冠层切面轮廓的二维激光点云数据；

惯性测量单元，用于获取激光雷达传感器的姿态角，并根据姿态角对所述二维激光点云数据进行矫正；其中，所述姿态角包括横滚角、俯仰角和偏航角；

卫星导航定位系统，用于获取激光雷达传感器的定位坐标，根据定位坐标计算两次激光扫描点之间的距离，融合定位坐标与矫正后的二维激光点云数据重构植物冠层三维激光点云图；

运算单元，用于获取叶面积指数反演计算所需的相关参数，构建叶面积指数估算模型，运算模型得到植物的叶面积指数。

进一步的，所述激光雷达传感器的角度分辨率不低于1°，测距距离不小于3m。

进一步的，所述惯性测量单元包括三轴加速度计和三轴角速度计，惯性测量单元的数据输出速率不低于10Hz，动态测角精度不低于0.1°。

进一步的，所述卫星导航定位系统的数据输出速率不小于10Hz，平面精度不低于1cm，高程精度不低于2cm。

本发明相对于现有技术具有如下的有益效果：

1、本发明融合了激光雷达传感器、卫星导航定位系统、惯性测量单元的检测数据，并利用卡尔曼滤波方法实现数据同步，获得比单一传感器更多的信息，形成比单一信源更可靠、更完全的融合信息，避免了单一传感器的信息盲区，有利于提高叶面积指数检测的准确率。

2、本发明利用惯性测量单元实时检测激光雷达传感器的姿态角，根据姿态角对激光扫描采样的点云数据进行数据矫正，这样即使路面不平整，在使用车载激光雷达传感器检测植物叶面积指数时，也能获得较高的检测准确率。

3、本发明利用抗干扰能力强的激光雷达传感器作为测距设备，探测实时性更强，探测准确率也更高。

4、本发明在构建叶面积指数估算模型后，计算两株样本植物的叶墙面积系数、击中点比例系数、击穿点比例因子、反射率系数和修正值，取两次计算结果的平均值作为最终的计算值，进一步提高叶面积指数检测的准确率。

附图说明

图1为本发明实施例1的植物叶面积指数检测系统结构框图。

图2为本发明实施例1的植物叶面积指数检测系统结构示意图。

图3为本发明实施例1的植物叶面积指数检测系统中卫星导航定位系统重构的三维激光点云图像示意图。

图4为本发明实施例2的植物叶面积指数检测方法流程图。

其中，1-数据同步单元，2-激光雷达传感器，3-惯性测量单元，4-卫星导航定位系统，5-运算单元，6-车体，7-植物，8-扫描面。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

实施例1：

如图1和图2所示，本实施例提供了一种多传感器信息融合的植物叶面积指数检测系统，该系统包括数据同步单元1、激光雷达传感器2、惯性测量单元3、卫星导航定位系统4和运算单元5，激光雷达传感器2、惯性测量单元3、卫星导航定位系统4装载在一车体6上，数据同步单元1、激光雷达传感器2、惯性测量单元3、卫星导航定位系统4和运算单元5的具体功能如下：

所述数据同步单元1，用于基于卡尔曼(Kalman)的滤波方法实现激光雷达传感器2、惯性测量单元3和卫星导航定位系统4之间的数据同步；其中，激光雷达传感器2获取的数据是二维激光点云数据，惯性测量单元3获取的数据是激光雷达传感器2的姿态角，卫星导航定位系统4获取的数据是激光雷达传感器2的定位坐标。

所述激光雷达传感器2采用二维(2-Dimension，2D)激光雷达传感器，用于扫描植物7冠层侧面，获得植物7冠层切面轮廓的二维激光点云数据，扫描面8可在图2中看到。

进一步地，激光雷达传感器2的角度分辨率不低于1°，测距距离不小于3m。

在本实施例中，二维激光点云数据用极坐标表示，其由一系列扫描角度间隔为α的数据点构成，每个数据点由三部分组成，包括扫描角度θ，该角度θ对应的径向距离radialDis，探测目标点的表面反射率RSSI，第n次激光扫描采样的点云数据以坐标(θ_k,radialDis_k,RSSI_k)_n表示，其中k在1到激光雷达传感器最大扫描线数之间取值，扫描角度间隔α等于(θ_k-θ_k-1)。

所述惯性测量单元3采用MEMS(Micro-Electro-Mechanical System)-IMU(Inertial Measurement Unit)微机电惯性测量单元，用于获取激光雷达传感器的姿态角，并根据姿态角对所述二维激光点云数据进行矫正；其中，所述姿态角包括横滚角、俯仰角和偏航角。

进一步地，惯性测量单元3包括三轴加速度计和三轴角速度计，优选地，惯性测量单元3还可以包括三轴磁力计。

进一步地，惯性测量单元3的数据输出速率不低于10Hz，可直接输出激光雷达传感器的横滚角

、俯仰角β和偏航角ψ的角度值，动态测角精度不低于0.1°。

在本实施例中，根据姿态角对所述二维激光点云数据进行矫正，如下式：

其中，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样时激光雷达传感器的姿态角，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样的点云数据，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据，Fun1为矫正函数，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值。

所述卫星导航定位系统4采用GNSS(Global Navigation Satellite System)卫星导航定位系统，用于获取激光雷达传感器的定位坐标，根据定位坐标计算两次激光扫描点之间的距离，融合定位坐标与矫正后的二维激光点云数据重构植物冠层三维激光点云图。

进一步地，卫星导航定位系统4的数据输出速率不小于10Hz，平面精度不低于1cm，高程精度不低于2cm。

在本实施例中，根据定位坐标计算两次激光扫描点之间的距离，如下式：

d＝Fun2(Coord_n,Coord_n-1)

其中，Coord_n为激光雷达传感器第n次激光扫描点的定位坐标，Fun2为两点定位坐标求解距离d的函数。

融合定位坐标Coord_n与矫正后的二维激光点云数据(θ’_k,radialDis’_k,RSSI’_k)_n重构植物冠层三维激光点云图，三维激光点云图数据用坐标(Coord,(θ’_k,radialDis’_k,RSSI’_k))_n表示，三维激光点云图像如图3所示。

所述运算单元5，用于获取叶面积指数反演计算所需的相关参数，构建叶面积指数估算模型，运算模型得到植物的叶面积指数。

优选地，叶面积指数反演计算所需的相关参数包括三维激光点云图数据(Coord,(θ’_k,radialDis’_k,RSSI’_k))_n。

上述数据同步单元1和运算单元5可以通过包含处理器的计算设备实现，如计算机、手机、平板电脑等。

实施例2：

如图2所示，本实施例提供了一种多传感器信息融合的植物叶面积指数检测方法，所述方法包括以下步骤：

S1、选择一株样本植物，对样本植物进行人工取样测量，获取样本植物准确的叶面积指数LAI_real(Leaf Area Index)。

S2、启动电源，初始化激光雷达传感器、卫星导航定位系统和惯性测量单元，保证各传感器工作正常。

S3、基于卡尔曼的滤波方法实现激光雷达传感器、惯性测量单元和卫星导航定位系统之间的数据同步。

S4、触发激光雷达传感器工作，扫描植物冠层侧面，获得植物冠层切面轮廓的二维激光点云数据。

本步骤中，二维激光点云数据由一系列扫描角度间隔为α的数据点构成，每个数据点由三部分组成，包括扫描角度θ，该角度θ对应的径向距离radialDis，探测目标点的表面反射率RSSI，第n次激光扫描采样的点云数据以坐标(θ_k,radialDis_k,RSSI_k)_n表示，其中k在1到激光雷达传感器最大扫描线数之间取值，扫描角度间隔α等于(θ_k-θ_k-1)。

S5、触发惯性测量单元，获取激光雷达传感器的姿态角，并根据姿态角对所述二维激光点云数据进行矫正；其中，所述姿态角包括横滚角、俯仰角和偏航角。

本步骤中，根据姿态角对所述二维激光点云数据进行矫正，如下式：

其中，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样时激光雷达传感器的姿态角，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样的点云数据，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据，Fun1为矫正函数，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值。

S6、触发卫星导航定位系统，获取激光雷达传感器的定位坐标，根据定位坐标计算两次激光扫描点之间的距离，融合定位坐标与矫正后的二维激光点云数据重构植物冠层三维激光点云图。

本步骤中，根据定位坐标计算两次激光扫描点之间的距离，如下式：

d＝Fun2(Coord_n,Coord_n-1)

其中，Coord_n为激光雷达传感器第n次激光扫描点的定位坐标，Fun2为两点定位坐标求解距离d的函数。

融合定位坐标Coord_n与矫正后的二维激光点云数据(θ’_k,radialDis’_k,RSSI’_k)_n重构植物冠层三维激光点云图，三维激光点云图数据用坐标(Coord,(θ’_k,radialDis’_k,RSSI’_k))_n表示。

S7、利用三维激光点云图数据(Coord,(θ’_k,radialDis’_k,RSSI’_k))_n，计算植物的叶墙面积LWA(Leaf Wall Area)，击中点比例P_hit，击穿点比例P_through。

本步骤中，以植物行距为分界点，小于植物行距的激光扫描测距点为击中点，是大于植物行距的激光扫描测距点为击穿点。

S8、对样本植物构建叶面积指数估算模型，利用样本植物的叶面积指数LAI_real计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，如下式：

其中，LWA(Leaf Wall Area)为植物的叶墙面积，P_hit为击中点的比例，P_through为击穿点的比例。

S9、按照上述步骤S1～S8，再选择一株样本植物，计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，取两次计算结果的平均值作为最终的计算值。

S10、利用最终计算得到的叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，代入叶面积指数估算模型进行反演计算，得到单株植物的叶面积指数LAI，如下式：

综上所述，本发明融合了激光雷达传感器、卫星导航定位系统、惯性测量单元的检测数据，并利用卡尔曼滤波方法实现数据同步，获得比单一传感器更多的信息，形成比单一信源更可靠、更完全的融合信息，避免了单一传感器的信息盲区，有利于提高叶面积指数检测的准确率。同时，本发明实时检测激光雷达传感器的姿态角，根据姿态角对激光扫描采样的点云数据进行数据矫正，这样即使路面不平整，在使用车载激光雷达传感器检测植物叶面积指数时，也能获得较高的检测准确率。

以上所述，仅为本发明专利较佳的实施例，但本发明专利的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明专利所公开的范围内，根据本发明专利的技术方案及其发明构思加以等同替换或改变，都属于本发明专利的保护范围。

Claims (7)

1.一种多传感器信息融合的植物叶面积指数检测方法，其特征在于：所述方法包括：

基于卡尔曼的滤波方法实现激光雷达传感器、惯性测量单元和卫星导航定位系统之间的数据同步；

采用激光雷达传感器扫描植物冠层侧面，获得植物冠层切面轮廓的二维激光点云数据；其中，所述二维激光点云数据由一系列扫描角度间隔为α的数据点构成，每个数据点由三部分组成，包括扫描角度θ，该角度θ对应的径向距离radialDis，探测目标点的表面反射率RSSI，第n次激光扫描采样的点云数据以坐标(θ_k,radialDis_k,RSSI_k)_n表示，其中k在1到激光雷达传感器最大扫描线数之间取值，扫描角度间隔α等于(θ_k-θ_k-1)；

采用惯性测量单元获取激光雷达传感器的姿态角，并根据姿态角对所述二维激光点云数据进行矫正，如下式：

其中，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样时激光雷达传感器的姿态角，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样的点云数据，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据，Fun1为矫正函数，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值；所述姿态角包括横滚角、俯仰角和偏航角；

采用卫星导航定位系统获取激光雷达传感器的定位坐标，根据定位坐标计算两次激光扫描点之间的距离，融合定位坐标与矫正后的二维激光点云数据重构植物冠层三维激光点云图；

获取叶面积指数反演计算所需的相关参数，构建叶面积指数估算模型，运算模型得到植物的叶面积指数，具体包括：

选择一株样本植物，对样本植物进行取样测量，获取样本植物准确的叶面积指数LAI_real；

获取样本植物的三维激光点云图数据(Coord,(θ’_k,radialDis’_k,RSSI’_k))_n；其中，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值，Coord_n为激光雷达传感器第n次激光扫描点的定位坐标，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据；

对样本植物构建叶面积指数估算模型，利用样本植物的叶面积指数LAI_real计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，如下式：

其中，LWA为植物的叶墙面积，P_hit为击中点的比例，P_through为击穿点的比例；

按照上述步骤，再选择一株样本植物，计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，取两次计算结果的平均值作为最终的计算值；

利用最终计算得到的叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，代入叶面积指数估算模型进行反演计算，得到单株植物的叶面积指数LAI，如下式：

2.根据权利要求1所述的一种多传感器信息融合的植物叶面积指数检测方法，其特征在于：所述根据定位坐标计算两次激光扫描点之间的距离，如下式：

d＝Fun2(Coord_n,Coord_n-1)

其中，Coord_n为激光雷达传感器第n次激光扫描点的定位坐标，Fun2为两点定位坐标求解距离d的函数。

3.根据权利要求1所述的一种多传感器信息融合的植物叶面积指数检测方法，其特征在于：所述击中点是指小于植物行距的激光扫描测距点，所述击穿点是指大于植物行距的激光扫描测距点。

4.一种多传感器信息融合的植物叶面积指数检测系统，其特征在于：所述检测系统包括：

数据同步单元，用于基于卡尔曼的滤波方法实现激光雷达传感器、惯性测量单元和卫星导航定位系统之间的数据同步；

激光雷达传感器，用于扫描植物冠层侧面，获得植物冠层切面轮廓的二维激光点云数据；其中，所述二维激光点云数据由一系列扫描角度间隔为α的数据点构成，每个数据点由三部分组成，包括扫描角度θ，该角度θ对应的径向距离radialDis，探测目标点的表面反射率RSSI，第n次激光扫描采样的点云数据以坐标(θ_k,radialDis_k,RSSI_k)_n表示，其中k在1到激光雷达传感器最大扫描线数之间取值，扫描角度间隔α等于(θ_k-θ_k-1)；

惯性测量单元，用于获取激光雷达传感器的姿态角，并根据姿态角对所述二维激光点云数据进行矫正，如下式：

其中，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样时激光雷达传感器的姿态角，(θ_k,radialDis_k,RSSI_k)_n为第n次激光扫描采样的点云数据，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据，Fun1为矫正函数，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值；所述姿态角包括横滚角、俯仰角和偏航角；

卫星导航定位系统，用于获取激光雷达传感器的定位坐标，根据定位坐标计算两次激光扫描点之间的距离，融合定位坐标与矫正后的二维激光点云数据重构植物冠层三维激光点云图；

运算单元，用于获取叶面积指数反演计算所需的相关参数，构建叶面积指数估算模型，运算模型得到植物的叶面积指数，具体包括：

选择一株样本植物，对样本植物进行取样测量，获取样本植物准确的叶面积指数LAI_real；

获取样本植物的三维激光点云图数据(Coord,(θ’_k,radialDis’_k,RSSI’_k))_n；其中，θ为扫描角度，radialDis为角度θ对应的径向距离，RSSI为探测目标点的表面反射率，k在1到激光雷达传感器最大扫描线数之间取值，Coord_n为激光雷达传感器第n次激光扫描点的定位坐标，(θ’_k,radialDis’_k,RSSI’_k)_n为根据姿态角采样矫正后的第n次激光扫描采样的点云数据；

对样本植物构建叶面积指数估算模型，利用样本植物的叶面积指数LAI_real计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，如下式：

其中，LWA为植物的叶墙面积，P_hit为击中点的比例，P_through为击穿点的比例；

按照上述步骤，再选择一株样本植物，计算出叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，取两次计算结果的平均值作为最终的计算值；

利用最终计算得到的叶墙面积系数λ₁、击中点比例系数λ₂、击穿点比例因子λ₃、反射率系数λ₄和修正值R，代入叶面积指数估算模型进行反演计算，得到单株植物的叶面积指数LAI，如下式：

5.根据权利要求4所述的一种多传感器信息融合的植物叶面积指数检测系统，其特征在于：所述激光雷达传感器的角度分辨率不低于1°，测距距离不小于3m。

6.根据权利要求4所述的一种多传感器信息融合的植物叶面积指数检测系统，其特征在于：所述惯性测量单元包括三轴加速度计和三轴角速度计，惯性测量单元的数据输出速率不低于10Hz，动态测角精度不低于0.1°。

7.根据权利要求4所述的一种多传感器信息融合的植物叶面积指数检测系统，其特征在于：所述卫星导航定位系统的数据输出速率不小于10Hz，平面精度不低于1cm，高程精度不低于2cm。

Images