发布信息

【联盟博览】无人机遥感在农田信息监测中的应用进展

作者:本站编辑      2023-06-08 13:13:45     28

摘要:快速实时地掌握农田信息是实施精准农作的基础。以无人机为平台的低空遥感探测技术具有空间分辨率高、时效性强和成本低等特点可填补地面监测和高空遥感间的测量尺度空缺因此在农田信息精准监测领域具有广泛的应用前景。近年来随着无人机飞行平台稳定性增强、操作难度降低机载遥感设备的轻量化和多样化以及遥感数据处理技术的进步无人机遥感在农田信息监测领域得到了快速发展。本文对国内外相关研究成果进行了总结对常用遥感技术类型和数据处理方法以及具体应用方向和实施效果进行了综述并提出了当前存在的突出问题和未来的发展方向以期为推动无人机遥感在农田信息监测和精准农业中更广泛的应用提供依据。

实施精准农作的前提是实时、准确获取农田信息和作物生长动态。无人机遥感(UAV-based remote sensing)是一种基于无人机载遥感传感器综合应用无人驾驶飞行、遥感、遥测遥控和空间定位等技术来获取空间信息的技术方法。在精准农业领域无人机遥感可用于监测作物生长参数:植被覆盖度、叶面积指数、株高以及这些参数与产量间的相关性在辅助农作物决策方面无人机遥感可用于诊断作物营养状况和田间墒情信息对于病虫害的监控也有较高的准确性此外无人机遥感还可快速获取农田空间位置信息划分土地利用类型、定位农田边界和基础设施、量算种植面积等的精度远高于传统测量方法。

与高空遥感相比无人机遥感具有空间分辨率高、应用成本低、时效性强、可重复使用等特点与传统的田间定位监测相比无人机遥感具有信息采集快捷、空间覆盖率高等特点。尤其是其监测范围覆盖了精准农业所关心的田块尺度且可方便地获取不同阶段的作物和土壤信息在农田信息监测领域有着传统点位监测、航空和卫星遥感无法比拟的优势得到越来越广泛的研究和应用。

1、无人机飞行平台与遥感测量技术

无人机遥感以无人驾驶飞机为探测平台利用搭载的各种任务负荷(通常是成像光谱仪等非接触式监测设备)获取有关农田和作物的遥感数字信息通过对数据的后期处理、挖掘和建模来获取农作物长势、农田环境等信息。

1.1无人机飞行平台

在农田信息监测领域人们关注的重点是监测范围、空间分辨率和测量精度等。无人机的机型、载重量、航行高度、续航时间、飞行稳定性、航线规划算法等均对探测效果有重要影响。目前常用的无人机可分为固定翼、单旋翼和多旋翼无人机等机型。固定翼无人机飞行速度快、续航时间长、载荷较大但飞行速度难以调节且需要较大的起飞着陆场在作物生长密集的农田常无法提供足够的起降场地单旋翼无人机稳定性较差会干扰传感器检测精度已逐步被多旋翼机型替代多旋翼无人机的航速姿态可调、飞行稳定、能够定点悬停适合获取多重复、定点、多尺度、高分辨率的植被信息因而在农田信息监测中应用最为广泛。为了保障探测精度和空间范围在无人机起飞前需要进行合理的航线规划。传统的无人机航线规划优化算法主要包括动态规划法、导数相关法、最优控制法、最速下降法、泰森多边形法(Voronoi); 现代智能算法主要包括遗传算法、人工神经网络、蚁群算法、粒子群算法等。蚁群算法在解决复杂航线规划方面效果良好但收敛效率低且容易陷入局部最优问题目前的研究多采用改进遗传算法。

1.2机载传感器与测量技术

无人机遥感传感器种类繁多感知原理和获取数据类型也各不相同。农田信息监测领域的传感技术大体可分为成像光谱和空间构型测量两大类。

成像光谱技术的原理是不同波段的光波作用于样本会产生不同的光谱特征由此可反映作物的生理生化指标多用于农田作物覆盖区的识别、叶片色素、养分含量等生物化学指标的预测建模。成像光谱的感光方式有棱镜/光栅色散型、干涉型、滤光片型和计算机层析型等感知光谱波段包括可见光、可见-近红外(380~2500 nm)和热红外波段(2.5~14 μm)。根据输出波段数量与连续性不同成像光谱又可分为高光谱和多光谱。目前的研究中普遍采用高清数码相机获取的农田正射影像、多光谱相机获取的多波段反射率、高光谱仪获取的连续光谱反射率、热红外传感器获取的田间温度信息。

空间构型测量技术获取的是高精度的农田空间位置信息常用于识别农田表面的三维构型及其特征并以此来诊断作物株高、叶面积指数、地上部生物量等生物物理指标。根据测量方式和技术差异,又可分为基于激光直接测量的激光探测与测量(Light Detection and Ranging, LiDAR)和基于可见光摄影测量的基于运动的结构(Structure from Motion, SfM)方法。机载LiDAR主动发射脉冲遇到作物或地表后反射回波通过时间差计算距离结合惯性系统(Inertial Navigation System, INS)确定的飞机姿态和全球定位系统(Globe Positioning System, GPS)得到飞机位置解算农田各处的三维空间坐标。SfM方法以农田遥感影像为数据源基于像对间的特征匹配通过迭代光束平差过程(Iterative Bundle AdjustmentProcedure)求解相机方位及场景几何形态再通过引入地面控制点坐标(Ground Control Point, GCP)进行矩阵变换将农田空间点云转入现实世界坐标系。

2、无人机遥感数据处理与建模方法

成像光谱与空间构型数据在处理流程、建模方法和结果应用上均不相同前者以提取特征信息建立参数反演模型为目的后者则能够建立测区三维空间结构模型。成像光谱数据中的高光谱获取的是不限波段范围的窄间隔连续光谱数据因此数据量大波段提取和信息压缩是高光谱数据处理的首要任务。热红外获取的是2.5~14 μm波段的数据包括连续和非连续的数据。由于该波段与温度关系密切数据处理也以提取温度信息为主要目标。空间构型数据中, SfM方法基于可见光摄影测量解算后的数据是带有色彩信息的三维空间坐标基于激光直接测量的LiDAR技术获取的则仅有空间位置信息。

2.1成像光谱数据处理与建模

成像光谱技术中使用可见光成像和多光谱数据的通常直接选取特征波段或计算植被指数将主成分提取后的光谱特征信息与作物的各项生理生化指标建立反演模型。

高光谱因为波段密集数据量大需要先选择与农田信息密切相关的特征波段。信息量、投影、相似度计算及小波分解法均可作为特征波段选择的依据。基于信息量的波段选择可将高光谱影像灰度值(Digital Number, DN)的标准差可以看作最简单的信息量度量信息熵,波段熵值高则信息量丰富基于投影的选择包括主成分分析法、独立成分分析法等以压缩波段降低维数为目的相似度计算中波段的空间自相关程度高则包含的信息量大而相关程度高的波段间存在信息冗余小波分解法目前已有很多小波基可供选择(双正交样条、HararDaubechies、高斯、SymletMeyerColilet小波等), 选择某一小波基后即可对一维或二维信号进行离散或连续的小波变换。

提取出特征波段后的高光谱数据处理过程与可见光成像和多光谱类似可直接计算波段反射率或计算各项植被指数后采用一元回归和多元回归等方法建立其与作物生长参量间的线性、对数、二次、幂函数等模型;也可根据如植被光学辐射传输等模型通过全局敏感性分析法(Sobol法、二叉树遍历法、GLUE法、FAST)评价非线性模型中作物目标生长参数对于光谱信息的影响及相应的敏感波段通过查找表法、数值优化法、人工神经网络法等方法反演作物生长信息或者基于机器学习(基于神经网络、高斯过程回归、核岭回归、随机森林等)对农田信息反演建模。

热红外波段反演地温的方法主要有单窗算法、劈窗算法、多通道和多角度算法。单窗算法利用一个热红外波段反演地表温度劈窗算法利用相邻两个热红外波段是目前最为成熟的反演算法。研究者对热温数据分析处理通常使用相对热温指数包括目标器官或组织之间的热温差或热温比值目标区域监测员手指之间的指温差指数和指温比指数增加不同热像中热温数据的可比性。如林业中常用的边心温比就是伐桩横截面的边材与心材热温比值。由于红外辐射的衍射效应较可见光强所以图像对比度较弱边缘模糊目标区域分割边缘提取是一项难点。目前图像边缘检测一般利用图像边缘的一阶和二阶导数信息例如:梯度法模板算子法等被广泛使用的模板算子有Sobel算子、Prewitt算子和Laplacain算子等。

2.2空间构型数据处理与建模

SfM方法首先要对原始图像进行无畸变处理校正由相机镜头引起的失真获取无畸变图像序列相对定向过程选择匹配特征点多且居于飞行区域中央的两张影像作为初始像对通过特征匹配、前后方交会解算初始像对两幅影像的相机参数及同名坐标点初始值再利用光束法平差优化结果添加新影像利用已求得的三维点坐标,及影像间的匹配特征点空间后方交会估算新添加影像的相机参数、对应特征点坐标再次平差优化直到完成所有影像的添加最终完成由农田航拍影像到三维点云的构建。此过程生成的点云是离散的利用点云内插方法得到数字高程模型(Digital Elevation Model,DEM)后能覆盖整个测区。

LiDAR数据处理的首要工作是去除和平滑由于系统误差、目标物表面过亮或偶然因素导致的噪声点。随后修复点云漏洞及缺失部分在完整的点云基础上完成农田精细建模。

无人机遥感在农田信息监测中的应用

无人机可搭载的遥感传感器多种多样可以获取多维度、高精度的农田信息实现多类农田信息的动态监测。这些信息主要包括作物空间分布信息(农田定位、作物种类识别、面积估算及变化动态监测、田间基础设施提取)、作物生长信息(作物表型参数、营养指标、产量), 以及作物生长胁迫因子(田间墒情、病虫害)动态等。

3.1农田空间信息

农田空间位置信息包括田块的地理坐标及通过目视判别或机器识别得到的作物分类。通过地理坐标识别出田块边界还可以实现种植面积的估算。传统的方法通过以地形图作为底图进行数字化开展区域规划和面积测算时效性差边界位置与实际情况差异巨大且缺乏直观性不利于精准农业的实施。无人机遥感可以实时获取全面的农田空间位置信息具有传统方法无可比拟的优势。高清数码相机的航拍影像即可实现农田基本空间信息的识别和判定如白由路等[8]在上海将凤农场仅以数码相机即实现了作物种类的识别和地块面积的量算。空间构型技术的发展提高了农田位置信息研究的精度与深度在引入高程信息的同时提升了空间分辨率可实现更精细的农田空间信息监测。如张宏鸣等将无人机DEM数据用于农田灌溉渠系的提取渠系提取的完整度达85.61%

3.2作物生长信息

作物生长状况可以通过表型参数、营养指标以及产量等信息来表征。表型参数包括植被覆盖度、叶面积指数、生物量、株高等。这些参数相互关联、共同表征了作物的长势情况与最终产量直接相关。在农田信息监测研究中占有主导地位已经开展的研究相对较多。

3.2.1作物表型参数

叶面积指数(Leaf Area Index, LAI)是指单位地表面积上单面绿叶面积的总和可较好地表征作物对光能的吸收利用与作物的物质积累和最终产量关系密切。叶面积指数是目前无人机遥感监测的主要作物生长参数之一。以多光谱数据计算植被指数(比值植被指数、归一化植被指数、土壤调节植被指数、差值植被指数等)与地面实测数据建立回归模型是反演表型参数较为成熟的方法。高林等通过对多个生育期、多种植被指数和不同模型的比较选择鼓粒期(大豆主茎最上部4个具有充分生长叶片着生的节中任何一个节位上豆荚内绿色种子充满荚皮的时期)归一化差分植被指数(Normalized Difference Vegetation Index, NDVI)的线性回归模型反演大豆LAI, 决定系数R2=0.829, 均方根误差RMSE=0.301, 估测精度EA=85.4%。也有利用可见光图像估测LAI的研究如高林等[36]构建了基于可见光大气阻抗植被指数(Visible Atmospherically Resistant Index, VARI)原理的数字图像特征参数(UAV-based VARIRGB)的指数模型, R2也达到0.71。高光谱的高分辨率优势为研究者提供了更丰富和连续的数据。随着高光谱传感器的推广和高阶数据处理方法的发展应用高光谱估算LAI的研究逐渐增多。已有研究证实携式地物光谱仪(Analytica Spectra Devices, ASD)获取的地面高光谱比值植被指数(Ratio Vegetation Index, RVI)对数模型的LAI预测能力优于无人机多光谱的NDVI线性模型[37]; Cubert UHD 185-Firefly(UHD185)是新型的无人机载高光谱传感器研究者通过对冬小麦孕穗、开花、灌浆期的UHD185高光谱影像与冠层ASD反射率的比较发现其在第3波段~96波段(458~830nm)具有较好的光谱质量。采用偏最小二乘回归法(Partial Least Squares Regression, PLSR)与红边参数结合估算叶面积指数独立验证R2=0.757, RMSE=0.732;交叉验证R2=0.755, RMSE=0.762。田明璐等针对传统固定波段植被指数存在的波段范围问题通过动态搜索植被指数将波段范围内的反射率极值定义为极值植被指数提高了棉花LAI的估测精度(验证R2最大提高了0.11)

作物生长后期地上部生物量与产量和品质的关系均很密切。目前农业上用无人机遥感进行生物量估测仍多使用多光谱数据提取光谱参数、计算植被指数进行建模空间构型技术在生物量的估算方面有一定优势如王东亮以呼伦贝尔草地为研究对象提出基于无人机LiDAR技术的草层高和盖度提取方法并用这两项参数反演了地上生物量(R2=0.784, RMSE=108.9 g·m-2)。该研究还探讨了无人机飞行高度对草层高度和盖度提取结果的影响并应用镶嵌算法提升了图像拼接的效率和效果对于农田作物生物量的估算具有参考意义。Bendig等利用SfM算法获取作物表面模型(Crop Surface Models, CSM)提取作物冠层高度结合3种可见光区植被指数来估算大麦生物量发现该方法在抽穗前期可靠但生长后期预测效果不佳。可见生育期对于建模参数的选择有很大影响。陆国政等在对大豆生物量的反演过程中采取了分段建模的方式。在开花结荚期以NDVI705和优化土壤调节植被指数(Optimization of Soil Adjusted Vegetation Index, OSAVI), 红边位置辅以株高为自变量通过最小二乘法建立多元线性回归模型独立验证R2=0.727, RMSE0.145;交叉验证R20.714,RMSE0.393;在生长后期(即鼓粒成熟期),由于株高稳定、对生物量影响小不再作为建模参数4种高光谱植被指数建立的生物量回归模型独立验证R20.698, RMSE0.238;交叉验证R2.697,RMSE0.386

3.2.2作物营养指标

传统的作物营养状态监测需要通过田间取样、室内化学分析以诊断营养物质或指标(叶绿素、氮素等)的含量而无人机遥感则依据不同物质具有特异的光谱反射吸收特征进行诊断。叶绿素的监测依据是其在可见光波段有两个强吸收区640~663 nm的红光部分和430~460nm的蓝紫光部分而在550 nm处吸收很弱。作物缺素时叶片颜色、纹理特征均会变化发掘不同缺素情况对应的颜色和纹理的统计特征及相关特性是营养监测的关键。与生长参数监测类似特征波段、植被指数和预测模型的选择依旧是研究的主要内容。肖宇钊对多光谱植被指数、纹理特征建立不同叶绿素的相对含量值(Soil and Plant Analyzer Development, SPAD)预测模型比较得出纹理特征易受成像质量影响稳定性差于优选植被指数该研究还发现延后采集时间、增加采集高度、降低飞行速度均能提高模型预测精度。在高值经济作物研究中, Arai等利用无人机近红外影像监测茶树叶片氮含量优化茶叶采摘时间在保持茶叶口感的同时提高收获量显著提高了经济效益。植被辐射传输机理模型可描述光在作物叶片和冠层吸收、反射的物理过程模型以作物生理信息为输入参数输出模拟的冠层光谱信息。通过查找表法、数值优化法、人工神经网络等方法可以反演作物的生长信息。高精度植被辐射传输机理模型被越来越多的研究者所使用。

3.2.3作物产量

提高农作物产量是农业活动的主要目标准确估测产量对农业生产和管理决策部门均有重要意义。众多研究者通过多因素分析试图建立更高预测精度的估产模型。如许童羽等从监测时间频率入手利用二元定距变量方法对单天、各旬、各月冠层NDVI与产量进行相关性分析相关性最优的五个时间范围内的数据用线性回归、二次、三次曲线分别与产量建模比较独立建模与组合建模的结果得到了预测最优水稻产量的模型是6月中旬和8月上旬的组合模型R2=0.771, 相对误差为4.06%, RMSE0.474 t·hm-2, 对应的北方粳稻生育期为分蘖盛期和抽穗期。赵晓庆等[14]的研究则着眼于监测的采样范围基于不同空间尺度提取的高光谱数据采用偏最小二乘回归构建植被指数与产量模型通过模型方程估算精度的曲线变化趋势分析最优空间尺度面积及采样小区的长宽比最终确定采样空间长宽与种植小区长宽比例介于4.25:54.5:5估算模型最为准确(r=0.8117)

3.3 作物生长胁迫因子

3.3.1农田墒情

农田墒情常通过热红外法进行监测。植被高度覆盖区因叶片气孔的关闭可降低因蒸腾造成的水分损失使地表潜热通量降低地表感热通量增加进而引起冠层温度升高地表温度可认为是植物冠层温度。由于反映作物能量平衡的水分胁迫指数可量化作物含水量与冠层温度的关系因此通过热红外传感器获取的冠层温度即可反映农田水分状况裸土或植被覆盖度小的地区可以用下垫面温度间接反演土壤水分其原理是:水的比热大受热温度变化慢因此白天下垫面温度的空间分布可间接反映土壤水分的分布。在对冠层温度的监测中裸露的土壤是重要的干扰因素。有研究者对裸土温度与作物地表覆盖度的关系进行了研究明确了由裸土导致的冠层温度测量值与真值间的差距将修正结果用于农田水分的监测提高了监测结果的精度。在实际农田生产管理中田间水分渗漏情况也是被关注的重点已有研究利用红外成像仪监测灌溉渠道水分渗漏精度可达93%

3.3.2病虫害

利用近红外光谱反射率监测植物病虫害其依据是:叶片在近红外区的反射受海绵组织和栅栏组织控制健康植物这两种组织间隙充满水分而膨胀是各种辐射良好的反射体当植物受害,叶片遭到破坏组织萎蔫水分随之减少红外反射减少直到丧失。

热红外监测的温度也是反映作物病虫害的重要指标。植株在健康条件下主要通过控制叶片气孔开闭进行蒸腾作用的调节维持自身温度的稳定在遭受病害后会发生病理变化病原物-寄主互作中病原物对植物的影响特别是对蒸腾有关方面的影响会决定侵染部分温度的升降。一般而言植物受感会导致气孔的开张调节失调因而病变区域的蒸腾作用高于健康区域旺盛的蒸腾作用会导致感染区域温度的下降叶面温差较正常叶片高直到叶片表面出现坏死斑点。坏死部位的细胞完全死亡该部分的蒸腾作用完全丧失温度开始升高但是由于叶片其他部分开始感染所以叶面温差始终高于健康植株。

美国农业部使用机载MS4100多光谱相机获取棉花影像采用多元线性回归和统计学方法分析、处理生成处方图暗紫色和蓝色区域棉花长势良好粉红色和红色区域棉花生病或者死亡监测结果直观可以快速投入应用辅助农艺管理如根据棉田病害处方图进行喷药管理。

3.4 其他信息

在农田信息监测领域无人机遥感数据还有更广泛的应用。如:利用多项纹理特征提取玉米倒伏面积[56]; 利用NDVI指数反映棉花成熟期叶片老熟程度生成的脱落酸施用处方图可有效指导棉花脱落酸的喷施避免农药过度施用等。根据农田监测与管理的需求不断挖掘无人机遥感数据信息、拓展其应用领域是未来信息化、数字化农业发展的必然趋势。

存在问题与发展方向

目前利用无人机监测农田信息还处于起步阶段不仅在飞行平台的研发、应用和管理方面存在不少问题而且在遥测数据的挖掘、处理、和综合应用方面也具有极大的提升空间。首先无人机的一次性投入成本较高飞行稳定性受气象条件、操控人员能力、甚至飞行政策等因素的制约真正投入生产实践的案例还比较少针对大田作物生长的阶段性特点现有的研究对于周期性的农田信息获取仍急需完善多载荷无人机虽然已经投入使用但获取的数据尚未充分实现多源数据的综合应用由于获取遥感数据缺少快速的处理方法和管理技术指导平台高精度、高阶数据的有效信息未能完全发掘。

4.1简单易用、低成本的无人机遥感平台

稳定的无人机平台与协调的传感器配置对于提高农田信息监测水平意义重大但是目前无人机存在成本过高、稳定性较差、受大风、阴雨等恶劣天气的限制等问题遥感平台的操作复杂且过度依赖于人工设置制约了无人机遥感在农田监测管理中的广泛应用。目前的研究表明无人机稳定性受到阵风的影响最大在飞行控制中也均考虑了风场扰动。滑模变结构控制和基于自抗干扰的飞行控制均取得了一定进展。商用无人机也通过改进自稳系统和航线规划算法简化了无人机的操作[58-59]。目前我国对于7 kg以下、飞行半径500 m范围内、高度120 m以下的无人机未作明确规定但飞行器的飞行空域依然要接受相关部门的监管。随着无人机的普及相关政策也有待明确。

4.2作物生长情况的周期性动态监测

对于农田信息监测缺少周期性和连续性是当前研究存在的普遍问题。由于作物生长情况处于动态变化中每一生长阶段的生长参数与遥感信息间的关系均不相同单一或少数几个生育期的作物生长模型不具有普适性和代表性。例如由于不同生育期的适用模型不同有大豆的鲜生物量的研究只能采取分段建模的方式。而且对于同一生长参数不同时期监测的准确性也存在较大差异因此对于作物生长状态的连续监测探寻同一生长参数不同生育期的变化趋势并构建适用于多生育期的普适模型十分必要。

4.3多源数据的综合应用

单一来源的遥感数据难以全面反应农田信息。随着传感器的轻型化和无人机载重及续航时间的增加已经实现多源数据同步监测农田信息。如何将不同遥感信息综合使用提高监测精度拓展监测范围是需要进一步思考的内容。目前已有以激光雷达为核心集成高分辨率相机、热成像仪、高光谱成像仪等传感器的Crop3D平台实现株高、株幅、叶长、叶宽、叶倾角及叶面积等作物各生长时期多源表型数据的提取实现高通量作物表型的测量为基因组学和生物学分析提供数据支持。也有将高清数码影像与高光谱数据结合提取株高和光谱参数综合应用于生物量估算的案例。但多源数据综合应用的方法和范围还很局限如何实现空间构型数据与成像光谱数据的融合无人机遥感数据与地面、卫星遥感数据的融合仍是需要深入研究的问题。

4.4高精度高阶数据实时处理与管理指导技术

相比于整个田块农田中植被覆盖区的遥感信息与作物状态的相关性更高。基于植被指数的遥感处理方法并未充分利用无人机遥感超高的空间分辨率。因此利用高空间分辨率的优点在田块内精细区分作物、垄间裸土获取单纯植物覆盖区的遥感信息对于提升作物模型的预测水平有实际意义。目前虽然基于可见光或光谱数据即可实现物体识别与分类但由于其受阴影等影响实际应用中存在不少问题。特别是作物生长前期由于地表覆盖度低遥感数据易受土壤介质的影响生长后期易受阴影等因素干扰将物体的空间特征也同时应用于农田物体识别分类将会有助于提高无人机遥感对作物的分辨能力实现对农田区域更精准的分割。

如何处理无人机遥感获取的海量数据也是亟待解决的问题。无法实现即时的数据解译指导农艺操作无人机遥感将丧失实时性的优势。因此无人机数据处理平台的开发尤为必要生成相应的农业管理措施更好地指导农业生产。

结语

无人机遥感作为一种新型的高空间、时间分辨率的信息监测技术,在农田信息监测中的研究和应用方兴未艾目前尚存在不少的问题。但与传统地面定点监测和高空遥感监测相比无人机遥感监测的空间尺度和精度优势显著尤其适合中尺度农田的数字信息快速获取。在可以预见的未来通过无人机飞行平台、机载多源信息采集技术、数据挖掘和建模技术、决策支持技术平台等方面的发展完善无人机遥感技术有望在农业和相关领域中得到更广泛、更深入的应用。

本文来源:土壤学报

相关内容 查看全部