基于LiDAR的温室番茄冠层几何参数提取
Geometric Parameters Extraction of Tomato Canopy in Greenhouse Based on LiDAR
通讯作者: 王学农(1964-),男,陕西汉中人,研究员,硕士生导师,研究方向为农业机械化技术装备,(E-mail)xjwxn2010@sina.com
收稿日期: 2020-11-20
| 基金资助: |
|
Corresponding authors: WANG Xuenong (1964-), male, from Hanzhong, Shanxi Province, researcher, master supervisor, research direction : agricultural mechanization technology and equipment, (E-mail)xjwxn2010@sina.com
Received: 2020-11-20
| Fund supported: |
|
作者简介 About authors
杨征鹤(1994-),男,山东菏泽人,硕士研究生,研究方向为设施农业装备,(E-mail)247329956@qq.com
【目的】准确获取温室番茄作物行中单株冠层数据,为分析作物生长状态和为对靶喷药提供冠层数据支持。【方法】采用三维激光雷达(LiDAR)搭建番茄植株冠层检测平台,使用导轨以0.05 m/s的速度移动三维激光雷达,利用雷达上位机软件Ctrlview保存双侧扫描的A、B 2组共40株番茄植株点云。双侧点云使用ICP(迭代最近点)算法进行配准,利用基于特征值的平面拟合法去除地面,使用均值漂移算法(Meanshift)分割番茄行中的单株点云,获取冠层参数,与人工测量值比较验证精度,将单株点云在MATLAB中使用alpha shape算法进行重建并进行体积的获取,使用凸包算法作物参考值对比。【结果】该检测平台在激光雷达前进方向与垂直前进方向的测量误差分别为-2.65%、-3.95%;获取到的单株番茄植株高度与人工测量值相比,平均绝对误差分别为0.025和0.031 m;重建后求取的体积与凸包算法相比平均误差下降了约15.3%,与人工获取相比相差不大,各指标良好。【结论】番茄行点云分割结果与人工测量相比A、B 2组的均方根误差RMSE分别为0.039和0.043,冠层体积获取与参考值对比VRMSE为0.011 3,激光雷达在获取作物外形轮廓信息中具有一定的准确性和可靠性,该方法用于温室环境下单株作物冠层数据的获取。
关键词:
【Objective】 To achieve the accurate acquisition of the canopy data of single plant in the greenhouse tomato crop row, to analyze the crop growth status and to provide canopy data support for target spraying. 【Methods】 A 3 d LiDAR was used to build A detection platform for tomato plant canopy. A 3D LiDAR was used to move the 3D LiDAR at A speed of 0.05 m/s. A total of 40 tomato plant point clouds were saved by Ctrlview, the upper computer software of the radar. To bilateral point cloud using ICP (iterative closest point) algorithm for registration, the plane fitting method based on characteristic value is to remove the ground, using the mean shift algorithm (Meanshift) line segmentation of tomato yield point in the cloud, obtain canopy parameters, and comparing with manual measurement verification accuracy, finally will yield point cloud using alpha in MATLAB algorithm based on the shape reconstruction and volume, and compared using convex hull algorithm crops reference. 【Results】 The test results show that the measurement errors of the platform in the forward direction and the vertical direction of the lidar are -2.65% and -3.95% respectively. The average absolute error M 12 was 0.025 m and 0.031 m, respectively, compared with the measured height. The average error of the volume obtained by reconstruction using Alpha Shape algorithm is about 15.3% lower than that of the convex hull algorithm, which is not much different from that obtained manually, and the indexes are good. 【Conclusion】 The mean shift algorithm is adopted to tomato line point cloud segmentation result compared with the artificial measure A, B two groups of root mean square error of the RMSE are 0.039 and 0.043 respectively, using the alpha shape algorithm canopy volume access and reference comparison VRMSE 0.011,3, shows that the laser radar in crop silhouette information has A certain accuracy and reliability of this method can be used to plant crop canopy of greenhouse environment data acquisition.
Keywords:
本文引用格式
杨征鹤, 喻晨, 杨会民, 陈毅飞, 周欣, 马艳, 王学农.
Zhenghe YANG, Chen YU, Huimin YANG, Yifei CHEN, Xin ZHOU, Yan MA, Xuenong WANG.
0 引言
【研究意义】植物冠层数据是分析植物生长状态的重要参数,利用传感器技术检测植物位置信息及轮廓外形是对靶精准喷雾技术的基础[1],探测作物几何信息是对靶喷雾技术的重要研究内容[2],植物冠层体积特征也是调整农药用量的关键。温室中环境较为稳定,更加适合使用激光雷达采集作物信息。获取温室环境中番茄作物冠层数据,搭建基于激光雷达的温室番茄冠层检测平台,研究提出1种基于点云的单株番茄分割和体积获取方法,对实现从番茄作物行中获取单株番茄植株冠层数据与体积有实际意义。【前人研究进展】检测冠层表征采用的有超声波传感器技术、立体视觉技术、光传感器技术及激光雷达传感器技术[3]。超声波技术已经在计算作物冠层体积[4]和估算冠层密度[5]得到应用,Maghsoudi等[6]使用超声波传感器进行果树冠层量的有效估计,以减少果园中农药用量。翟长远等[7]基于超声波传感器设计了1种果树靶标体积的探测方法。但在实际使用中超声波传感器存在靶向性及采样频率较低等问题,存在一定的误差。Sylvain等[8]利用图像技术和SFM技术建立了作物行三维模型,预估植物高度和叶面积。立体视觉技术及光传感器受外界环境光影响较大[9],稳定性不足,在实际应用限制较多,相比之下利用激光雷达检测冠层表征是现阶段最高效的技术。激光雷达传感器具有速度快、精度高、靶向性强、非接触性、受外界环境光影响较小等优点[10],在预测冠层表征方面较为可靠,并已经广泛应用[11,12,13],且在不规则植株冠层特征信息的获取中有着难以比拟的优势[14]。Jordi等[15]验证了二维激光雷达测定植株体积和密度相关的冠层参数的精度,俞龙等[16]基于二维激光雷达设计了1种对丘陵山地果树冠层体积的激光测量方法,吴志鹏等[17]使用激光雷达对果树轮廓进行测量,试验结果证明,该方法具有一定的可行性。张美娜等[18]使用二维激光雷达设计了1种靶标探测系统试验平台,提出了1种基于点云信息的靶标叶面积密度的计算方法。三维雷达相较于二维雷达具有更高的信息获取能力,Ravi等[19]将LiDAR与无人机相结合,获取作物冠层覆盖情况并建立作物冠层地图,管贤平等[20]将三维激光雷达与导航系统、无人机相结合获取大田大豆单株几何冠层几何参数。【本研究切入点】激光雷达现阶段多用于针对森林或果园中林木及果树的冠层数据提取,很少在温室环境中对蔬菜类作物冠层进行分析。研究基于LiDAR的温室番茄冠层几何参数提取。【拟解决的关键问题】以番茄作物行为目标,设计基于LiDAR的温室番茄冠层检测平台,在番茄作物行两侧匀速运行激光雷达,获取番茄植株的三维点云信息,对点云进行去噪、配准处理,从行点云中分割获取单株点云,获取冠层信息,将点云重建并获取体积信息,验证获取的三维轮廓的准确性,比较与凸包算法和人工测量值,验证装置的精准性。
1 材料与方法
1.1 材料
1.1.1 番茄
选择新疆乌鲁木齐西山的一处温室番茄种植园内,选择的番茄作物生长期分别为40、80 d,高约0.6和1.5 m,株距同为0.3 m,分别定义为A组和B组。
1.1.2 仪器设备
试验平台主要由激光雷达部分、支架部分、导轨部分组成,其中导轨部分主要由滑块、线轨部分、步进电机、控制器及电源部分组成,激光雷达部分直接与PC机相连。图1
图1
图1
激光探测平台结构示意
注:1.激光雷达部分;2.支架部分;3.导轨部分
Fig.1
Structure diagram of detection platform
Note: 1. Lidar part; 2. Stents part ;3. Guide part
选用锐博公司生产的XG80导轨平台,该导轨平台具有承载能力强、精度较高、稳定性好的优点。滑轨总长度为1 500 mm,滑块通过步进电机驱动,精度良好,滑块移动速度为10~50 mm/s,滑块所连接的支架高度可调,以适应不同的作业环境。
选用北科天绘公司的R-Fans-32线程激光雷达。激光帧频在5~20 Hz可调,最大探测距离为200 m,距离分辨率为4 mm,激光波长为905 nm,测程远,回波强度较为准确,测量精度较高,角度覆盖和角分辨率,有效的抵抗背景光的干扰。
1.2 方法
番茄作物行点云通过获取、配准等一系列处理,最后获取单株点云并通过MATLAB进行重建与体积的计算。图2
图2
1.2.1 番茄行点云获取
图3
从XY平面方向对单株点云直接测量,计算点云中的最高点和最低点的差值作为激光雷达获取的番茄植株的冠层高度
点云中每个点的坐标计算:
式中:
图4
图5
1.2.2 点云去噪与匹配
点云去噪是通过为每个点搜索指定邻域点个数的相邻点,计算其与邻域点的平均距离,并计算这些平均距离的中值和标准差。如果该点的平均距离大于最大距离,则将其视为噪声点,将其去除。
利用ICP(Iterative Closest Points)即迭代最近点算法对双侧点云进行匹配。ICP配准算法是比较常用的点云配准方法之一,其基本原理为在P、Q两目标点云中,按照一定的约束找到最近的邻近点,计算出平移参数T和旋转参数R,对Q使用求得的平移和旋转参数,得到新的变换点集;如果新的变换点集与参考点集满足2点集的平均距离小于某一给定阈值,则停止迭代计算,否则新的变换点集作为新的Q继续迭代,直到达到目标函数的要求[22],目标函数为:
1.2.3 地面分割
式中:a、b、c为单位法向量,其中d为坐标原点到该拟合平面的距离,则对某一平面扫描的n个数据点的坐标为(
并且能够满足:
1.2.4 单株番茄点云分割算法
式中:以x为中心点,半径为h的高维度球区域,k为包含在
式中:
使用LiDAR360软件直接从XY平面方向对单株点云直接测量,计算点云中的最高点和最低点的差值作为激光雷达获取的番茄植株的冠层高度HL,人工测量值为从平行于导轨方向使用卷尺从每株番茄中获取冠层高度HM。图6
图6
1.2.5 单株番茄植株体积计算
选取的alp值为0.1,其中alp值越小,网格化体积越小,得到的点云体积更加精确。验证使用alpha shape算法求取体积的精度,使用凸包算法(Convex hull)求取体积,作为参考值对比。凸包算法能够对单株点云数据建立一个能够将所有点云都包含在内的最小整体凸包。
由于番茄植株具有并不规则的外形,枝叶相互遮挡,在计算体积时,需要对番茄植株点云数据进行三维重建。该次试验使用alpha shape算法对番茄植株的三维点云进行重建[28],alpha shape算法能够在离散的三维点云空间中提取边缘并建立对应最小包络。从A组和B组各选取10株番茄进行点云重建,选取的alp值为0.1,其中alp值越小,网格化体积越小,得到的点云体积更加精确,从而更精准的对番茄植株进行重建,最后使用Volume函数对所重建后包络进行体积的计算[29]。为验证使用alpha shape算法求取体积的精度,使用凸包算法(Convex hull)求取体积,作为参考值进行对比。凸包算法能够对单株点云数据建立一个能够将所有点云都包含在内的最小整体凸包,并进行体积的计算。
使用人工测量的方式对随机5株番茄植株进行体积获取。在单株体积计算时,获取的5个位置的冠层宽度,将5处位置的横截面近似为椭圆,求出椭圆面积与每个部分的冠层高度相乘,并通过累加得到番茄植株的冠层体积。
1.2.6 精度验证
使用人工测量的方式对随机5株番茄植株进行体积获取。激光雷达获取的冠层体积为测量值
式中:
2 结果与分析
2.1 规则框架人工测量值与系统检测值对比
图7
图7
规则框架示意
注:l1、l3方向为激光雷达前进方向,l2、l4方向为垂直于激光雷达前进方向
Fig.7
Schematic diagram of the rule framework
Note: l1 and l3 are the forward direction of the lidar; l2 and l4 are the direction of vertical lidar
表1 规则框架人工测量值与系统检测值对比
Table 1
| 参数 Parameter | 人工 Manual (cm) | 平台 Platform (cm) | 相对误差 Relative error (%) |
|---|---|---|---|
| 60.0 | 58.3 | -2.8 | |
| 80.0 | 76.4 | -4.5 | |
| 60.0 | 58.5 | -2.5 | |
| 80.0 | 77.3 | -3.4 |
2.2 冠层高度人工测量值与平台检测值对比
表2 冠层高度人工测量值与平台检测值对比
Table 2
| 编号 No. | 人工 Manual (m) | 平台 Platform (m) | 差值 Difference (m) |
|---|---|---|---|
| 1 | 0.70 | 0.73 | -0.03 |
| 2 | 0.73 | 0.75 | -0.02 |
| 3 | 0.59 | 0.64 | -0.05 |
| 4 | 0.50 | 0.56 | -0.06 |
| 5 | 0.66 | 0.73 | -0.07 |
| 6 | 0.75 | 0.76 | -0.01 |
| 7 | 0.55 | 0.58 | -0.03 |
| 8 | 0.45 | 0.49 | -0.04 |
| 9 | 0.57 | 0.54 | 0.03 |
| 10 | 0.81 | 0.84 | -0.03 |
| 11 | 0.69 | 0.74 | -0.05 |
| 12 | 0.57 | 0.58 | -0.01 |
| 13 | 0.70 | 0.70 | 0.00 |
| 14 | 0.60 | 0.63 | -0.03 |
| 15 | 0.61 | 0.58 | 0.03 |
| 16 | 0.49 | 0.45 | 0.04 |
| 17 | 0.53 | 0.59 | -0.06 |
| 18 | 0.62 | 0.65 | -0.03 |
| 19 | 0.60 | 0.65 | -0.05 |
| 20 | 0.62 | 0.65 | -0.03 |
| 21 | 1.67 | 1.68 | -0.01 |
| 22 | 1.59 | 1.63 | -0.04 |
| 23 | 1.6 | 1.71 | -0.11 |
| 24 | 1.74 | 1.76 | -0.02 |
| 25 | 1.60 | 1.58 | 0.02 |
| 26 | 1.65 | 1.60 | 0.05 |
| 27 | 1.45 | 1.45 | 0.00 |
| 28 | 1.60 | 1.64 | -0.04 |
| 29 | 1.72 | 1.75 | -0.03 |
| 30 | 1.63 | 1.64 | -0.01 |
| 31 | 1.50 | 1.56 | -0.06 |
| 32 | 1.54 | 1.56 | -0.02 |
| 33 | 1.60 | 1.65 | -0.05 |
| 34 | 1.54 | 1.57 | -0.03 |
| 35 | 1.67 | 1.69 | -0.02 |
| 36 | 1.69 | 1.65 | 0.04 |
| 37 | 1.68 | 1.67 | 0.01 |
| 38 | 1.70 | 1.74 | -0.04 |
| 39 | 1.68 | 1.69 | -0.01 |
| 40 | 1.74 | 1.73 | 0.01 |
图8
图9
图9
冠层高度检测值与测量值的线性模型
Fig.9
Linear model of detected and measured values of canopy height
2.3 番茄冠层重建与体积
表3 番茄冠层体积平台检测值与参考值
Table 3
| 编号 No. | 检测值 Detection value (m3) | 参考值 Reference value (m3) | 绝对误差 Absolute error (m3) | 相对误差 Relative error (%) |
|---|---|---|---|---|
| 1 | 0.019 0 | 0.022 3 | -0.003 3 | 14.8 |
| 2 | 0.017 2 | 0.019 2 | -0.002 0 | 10.4 |
| 3 | 0.022 1 | 0.025 2 | -0.003 1 | 12.3 |
| 4 | 0.027 6 | 0.031 7 | -0.004 1 | 12.9 |
| 5 | 0.016 1 | 0.019 0 | -0.002 9 | 15.3 |
| 6 | 0.020 3 | 0.024 2 | -0.003 9 | 16.1 |
| 7 | 0.016 6 | 0.018 9 | -0.002 3 | 12.2 |
| 8 | 0.019 7 | 0.023 4 | -0.003 7 | 15.8 |
| 9 | 0.023 4 | 0.026 3 | -0.002 9 | 11.0 |
| 10 | 0.026 9 | 0.030 6 | -0.003 7 | 12.1 |
| 11 | 0.074 6 | 0.090 2 | -0.015 6 | 17.3 |
| 12 | 0.073 9 | 0.087 1 | -0.013 2 | 15.2 |
| 13 | 0.069 2 | 0.083 1 | -0.013 9 | 16.7 |
| 14 | 0.076 7 | 0.094 1 | -0.017 4 | 18.5 |
| 15 | 0.081 0 | 0.094 6 | -0.013 6 | 14.4 |
| 16 | 0.071 4 | 0.087 4 | -0.016 0 | 18.3 |
| 17 | 0.065 7 | 0.079 8 | -0.014 1 | 17.7 |
| 18 | 0.069 4 | 0.087 6 | -0.018 2 | 20.8 |
| 19 | 0.083 2 | 0.097 6 | -0.014 4 | 14.8 |
| 20 | 0.076 4 | 0.095 1 | -0.018 7 | 19.7 |
| 平均值 Average value | -0.009 4 | 15.3 |
图10
获取的5个位置的冠层宽度,将5处位置的横截面近似为椭圆,求出椭圆面积与每个部分的冠层高度相乘,并通过累加得到番茄植株的冠层体积。人工测量值与平台检测值的
表4 人工测量与平台检测体积对比
Table 4
| 编号 No. | 人工 Manual (m3) | 平台 Platform (m3) | 相对误差 Relative error (%) |
|---|---|---|---|
| 1 | 0.022 6 | 0.019 0 | 15.9 |
| 2 | 0.025 9 | 0.022 1 | 14.6 |
| 3 | 0.095 2 | 0.076 7 | 19.4 |
| 4 | 0.098 6 | 0.081 0 | 17.8 |
| 5 | 0.097 8 | 0.076 4 | 21.9 |
3 讨论
利用激光雷达进行植物冠层数据的分析已经具备较大的发展潜力,激光雷达能够快速直接获取植物冠层形状,信息准确度较高,受外界环境影响较小,是其他技术无法比拟的优势。在实际的使用中,激光雷达在使用时需要保证一定的稳定性,减少因为振动而造成的噪点过多,从而影响数据的准确性。
采用平面拟合算法对地面进行分割时对地形的要求较高,若地形起伏较大,则产生的误差较大。采用均值漂移算法进行单株番茄植株的分割时,当植株枝叶过于茂密,相邻植株间交织严重,激光雷达获取的点云密集,从而无法通过点云特征处理,应根据不同作物的生长状况,采用不同的分类方法对单个植物进行分割[30]。通过alpha shape算法获得的植株体积小于使用凸包算法获取的体积,是因为alpha shape算法能够减少冠层间空隙带来的影响。
4 结论
4.1 温室环境中,作物状态较为稳定,受外界环境影响较小,相比于大田环境更加适合使用检测平台,检测平台在前进方向与垂直前进方向的测量误差分别为-2.65%、-3.95%。
4.2 采用基于平面拟合和均值漂移算法以实现从点云中分割地面和从番茄行中分割单株番茄点云,与人工测量冠层数据值相比,A、B 2组的均方根误差RMSE分别为0.039和0.043,平均绝对误差分别为0.025和0.031 m,均取得了较好的结果。
4.3 对单株点云使用alpha shape算法进行重建并进行体积的计算,平均绝对误差为0.009 4 m3,与凸包算法获得结果相比,平均相对误差下降了约15.3%,与人工获取相比,
参考文献
果树农药精确喷雾技术
[J].
Precision spray technology of fruit tree pesticide
[J].
变量喷雾技术研究进展分析
[J].
无人机遥感解析田间作物表型信息研究进展
[J].
Advances in the analysis of field crop phenotypic information by remote sensing using uavs
[J].
Investigation of laser and ultrasonic ranging sensors for measurements of citrus canopy volume
[J].
A Canopy Density Model for Planar Orchard Target Detection Based on Ultrasonic Sensors
[J].
Ultrasonic sensing of pistachio canopy for low-volume precision spraying
[J].
树型喷洒靶标外形轮廓探测方法
[J].
In-field crop row phenotyping from 3D modeling performed using Structure from Motion
[J].
西红柿采摘机器人目标识别、定位与控制技术研究
[D].
Research on target recognition, location and Control Technology of Tomato picking Robot
[D].
模拟复杂地形的喷雾靶标激光检测与三维重构
[J].
Laser detection and 3D reconstruction of spray targets for complex terrain
[J].
Estimation of the crop density of small grains using Li DAR sensors
[J].
Laser scanning based three dimensional measurement of vegetation canopy structure
[J].
Estimating wheat green area index from ground-based LiDAR measurement using a 3D canopy structure model
[J].
基于三维激光扫描仪的校园建筑物建模研究
[J].
Research on campus building modeling based on 3d laser scanner
[J].
Testing the Suitability of a Terrestrial 2D LiDAR Scanner for Canopy Characterization of Greenhouse Tomato Crops
[J].
丘陵山地果树冠层体积激光测量方法与试验
[J].
基于激光雷达的枣树轮廓测量平台的设计与试验
[J].
Design and test of a lidar based profile measurement platform for jujube
[J].
基于三维激光点云的靶标叶面积密度计算方法
[J].
Chinese journal of agricultural machinery
[J],
Implementation of UAV-Based lidar for high throughput phenotyping
[C].
基于机载三维激光扫描的大豆冠层几何参数提取
[J].
Geometric parameters extraction of soybean canopy based on airborne 3d laser scanning
[J].
3D激光扫描技术及其在森林测绘中的应用研究
[J].
Research on 3D laser scanning technology and its application in forest mapping
[J].
基于四参数-ICP的点云配准研究
[J].
Point Cloud Registration Based on four-parameter ICP
[J].
基于RANSAC和TLS的点云平面拟合
[J].
Point cloud plane fitting based on RANSAC and TLS
[J].
基于最小二乘拟合的墙面平整度检测方法
[J].
Wall Flatness Detection Method based on least-squares Fitting
[J].
三种点云数据平面拟合方法的精度比较与分析
[J].
Precision Comparison and Analysis of three plane fitting methods for point cloud data
[J].
基于水平集与多图谱的医学图像分割算法研究
[D].
Research on medical Image segmentation Algorithm based on Horizontal set and Multi-atlas
[D].
基于交通视频的目标检测与跟踪算法研究
[D].
Research on Target Detection and Tracking Algorithm based on traffic video
[D].
A Method to Obtain Orange Crop Geometry Information Using a Mobile Terrestrial Laser Scanner and 3D Modeling
[J].
The ball-pivoting algorithm for surface reconstruction
[J].
/
| 〈 |
|
〉 |
