一种基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法与流程

本发明涉及一种基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法，属于农业信息技术领域。

背景技术：

作为提高农业用水效率的重要手段，调亏灌溉可以在不减产或少量减产的情况下节约大量农业用水，为了实现作物产量与灌溉用水之间的微妙平衡需要有效的作物水分胁迫状况监测手段。冠层温度作为作物冠层茎、叶表面温度的平均值，是衡量作物水分胁迫状况的重要指标，在国内外得到了广泛的应用。然而，传统的手持或地面红外测温传感器监测作物冠层温度时，具有范围较小、成本高等缺点，且没有考虑土壤、作物等的空间变异性，利用采样点的数据测量代表整个实验区域的水分胁迫状况，无法满足目前农场尺度的精细作业需求。

近年来，相比于卫星遥感和航空遥感技术，无人机凭借其平台凭借成本低、易于构建、运载便利、灵活性高、作业周期短、影像数据时空分辨率高等优势，可以低成本地多次重复采集农业遥感数据，已经成为农业遥感研究中快速、准确获取作物信息的主要工具，是当前研究的热点和趋势。由于受到微型无人机遥感平台载重能力的限制，目前搭载的热红外传感器多为重量较轻的非制冷热红外相机，存在冠层温度采集不稳定、拼接精度低、需要复杂预处理等问题，还需要对基于热红外遥感图像的作物冠层温度提取进行深入研究。

技术实现要素：

为了解决上述问题，本发明提供了一种基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法，即一种通过基于高分辨率可见光遥感数据获取大田玉米覆盖度，进行降维及二值化处理，然后采用矢量叠加的方式从热红外图像提取冠层温度的方法，可以低成本、快速获取大田作物水分胁迫状况分布图，为田间灌溉管理提供技术支持，促进农业信息化。

1、本发明提供一种基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法，包括以下步骤：

（1）无人机遥感数据采集航线设计

为了满足无人机可见光和热红外遥感数据的拼接要求，在进行航线设计时应满足航向及旁向重叠度大于85%，且可见光和热红外遥感图像的地面分辨率应满足以下要求：

(1)

式中gsdv为可见光遥感数据的地面分辨率（m），gsdt为热红外遥感数据的地面分辨率（m）。

无人机遥感图像的地面分辨率（gsd）是表征光谱数据的重要参数，由下式计算：

(2)

式中gsd为地面分辨率（m），a为多光谱传感器像元尺度（um），f为相机焦距（mm），h为无人机飞行高度（m）。

（2）地面像控点（groundcontrolpoint，gcp）布设：

为了实现无人机可见光和热红外遥感数据的几何校正和地理信息匹配处理，需要在采样区进行像控点布设，并利用全站仪进行三维坐标信息采集。像控点为正方形棋盘格形式，四个黑白格分别成对角放置，且尺寸大小与遥感图像地面分辨率满足以下关系：

(3)

式中l为正方形像控点边长，gsd为无人机遥感图像的地面分辨率。

（3）无人机遥感数据采集及预处理

选择晴朗天气，采用步骤（1）设计的航线以相机镜头垂直向下的方式分别进行无人机可见光和热红外遥感数据采集，基于同步采集、相同高度的黑（反射率3%）、白漫反射板（反射率58%）图像对可见光遥感图像进行辐射定标，获取地物波谱反射率。基于同步采集的包含黑、白漫反射板和水体的热红外图像提取上述三种地物的温度，并利用手持红外测温仪测量的上述三种地物温度建立热红外图像的温度校正模型并进行温度校正处理。

（4）无人机遥感影像的拼接、几何校正及地理信息匹配

基于步骤（3）采集的无人机可见光和热红外单张遥感数据及pos数据，利用遥感图像拼接软件得到采样区的可见光和热红外正射影像。在利用遥感图像拼接软件分别进行可见光和热红外遥感图像的拼接处理时，利用步骤（2）布设的地面控制点及测量的三维空间信息完成遥感图像的几何校正和地理信息匹配，使得获取的可见光和热红外遥感图像具有相同的地理信息。

（5）基于无人机可见光遥感的大田玉米覆盖度提取

在步骤（4）中获取的无人机可见光遥感数据上进行波段运算得到红绿比值指数rgri（red-greenratioindex）指数分布图，并应用otsu阈值算法获取大田玉米和地物的rgri分类阈值，小于阈值的为玉米，进而提取大田玉米覆盖度。所述红绿比值指数rgri由公式（4）计算得到。

(4)

式中和分别为红色波段和绿色波段的反射率。

然后采用最近邻算法对上述方法得到的大田玉米覆盖度分布图进行降维处理，得到热红外遥感影像分辨率尺度下的大田玉米覆盖度分布图并进行二值化。

（6）基于无人机热红外遥感影像的大田玉米冠层温度提取

利用步骤（5）获取的同分辨率尺度下的大田玉米覆盖度二值化图像，采用矢量叠加的方式从热红外遥感影像中提取大田玉米冠层温度。

本发明的优点和积极效果是：

本发明基于最新的无人机遥感技术，采用可见光和热红外遥感影像融合的方式实现了大田玉米的冠层温度的提取，可以有效的监测大田玉米的水分胁迫状况，为田间灌溉管理提供技术支持，促进农业信息化。

附图说明

图1基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法流程

图2实验地概况及样地划分

图3大田玉米冠层温度与实测温度的关系

图4大田玉米冠层温度与气孔导度的关系

图5大田玉米冠层温度分布图。

具体实施方式

以下结合附图对本发明做进一步详述。

如图1所示，这种基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法主要包括6个主要步骤：

（1）无人机遥感数据采集航线设计

（2）地面像控点（groundcontrolpoint，gcp）布设

（3）无人机遥感数据采集及预处理

（4）无人机遥感影像的拼接、几何校正及地理信息匹配

（5）基于无人机可见光遥感的大田玉米覆盖度提取

（6）基于无人机热红外遥感影像的大田玉米冠层温度提取

本发明基于最新的无人机遥感技术，采用可见光和热红外遥感影像融合的方式实现了大田玉米的冠层温度的提取，可以有效的监测大田玉米的水分胁迫状况，为田间灌溉管理提供技术支持，促进农业信息化。

以下结合实施例进行详细说明，本实施例的试验田域位于内蒙古鄂尔多斯市达拉特旗昭君镇（40°26'0.29"n，109°36'25.99"e，海拔高度1010m），种植作物为夏玉米，行距50cm，株距25cm，行向由东向西，种植面积20亩，采用中心轴式喷灌机进行灌溉，实验地被划分为5个扇形区域，每个区域布设1个土壤含水率采集节点（紫色三角形）和3个地面数据采样区（红色方框），如图2所示。以田间持水量的不同百分比在营养生长阶段（vegetativestage，v期）、生殖期（reproductivestage，r期）和成熟期（maturationstage，m期）进行不同梯度的水分处理，田间持水量的50%为充分灌溉，具体用水量设置梯度如表1所示。通过在喷灌机上安装流量计（mik2000h型）采集实际灌溉量，通过与实验地相邻的标准气象站采集实际降雨量。为了排除其它因素的干扰，根据种植经验进行统一施肥、施加除草剂。

表1用水量梯度设置

本发明步骤（1）具体操作如下：

根据试验田的区域大小（本实施例20亩）和续航时间大小进行了无人机可见光和热红外遥感数据的采集航线设计，无人机可见光遥感数据的采集高度为50m，航向及旁向重叠度为90%，地面分辨率为1.25cm；无人机热红外遥感数据的采集高度为60m，航向及旁向重叠度为85%，地面分辨率为7.8cm。满足航线设计的重叠要求，以及可见光遥感图像与热红外遥感图像的地面分辨率比例要求。

本发明步骤（2）具体操作如下：

由步骤（1）可知，无人机可见光遥感数据和热红外遥感数据的较低地面分辨率为7.8cm，为了方便计算取8cm，根据地面控制点的尺寸要求，正方形地面控制点的边长应大于等于1.6m，本实施例取1.6m。如图2所示，在试验田中布设了5个地面控制点并利用全站仪进行了三维坐标信息测量，满足地面控制点的数量布设要求。

本发明步骤（3）具体操作如下：

在玉米拔节期，即2018年7月12号（晴天）11:00-14:00，以flirvuepror640热红外相机镜头垂直向下，航向及旁向重叠度85%、飞行相对高度60m、速度5m/s的策略进行热红外遥感图像采集，图像地面分辨率7.8cm。飞行之前通过fliruavs软件将辐射率设置为0.97，图像格式设置为14位tiff，选择户外模式，根据具体天气状况设置天气条件、大气温度、空气温湿度等。在走航线采集热红外遥感数据之前，在60m高度采集包含两个黑（反射率3%）、白（反射率58%）漫反射板（尺寸3x3m，groupviii，usa）和水体的图像。同时，利用手持红外测温仪雷泰(raytek)st60+测量上述三种地物温度，建立温度校正模型并对热红外图像进行温度校正。以dji精灵4pro无人机镜头垂直向下，航向及旁向重叠度90%、飞行相对高度50m的策略进行可见光遥感图像采集，地面分辨率为1.25cm，根据具体天气条件设置相机iso值、快门时间（大于1250），白平衡设置为晴天。基于同步采集、相同高度的黑（反射率3%）、白漫反射板（反射率58%）图像对可见光遥感图像进行辐射定标，获取地物波谱反射率。

本实施例中热红外相机型号为：flirvuepror640，但不限于此；

本实施例中红外测温仪型号为：雷泰(raytek)st60+，但不限于此；

本实施例中无人机可见光遥感系统为：dji精灵4pro，但不限于此。

本发明步骤（4）具体操作如下：

遥感图像采集完成后，利用pix4dmapper软件进行图像拼接处理，并利用步骤（2）布设的地面控制点及其三维坐标信息对遥感图像进行了几何校正和地理信息匹配。

本发明步骤（5）具体操作如下：

在步骤（4）中获取的无人机可见光遥感数据上进行波段运算得到红绿比值指数rgri（red-greenratioindex）指数分布图，并应用otsu阈值算法获取大田玉米和其他地物的rgri分类阈值，小于阈值的为玉米，进而提取大田玉米覆盖度。所述红绿比值指数rgri由公式（4）计算得到。该大田玉米覆盖度分布图的地面分辨率为1.25cm，需要通过降维处理算法将地面分辨率降低为热红外遥感图像的地面分辨率7.8cm并进行二值化，本方法采用的降维算法为最近邻算法。

本发明步骤（6）具体操作如下：

利用步骤（5）获取的地面分辨率为7.8cm的大田玉米覆盖度二值化图像，采用矢量叠加的方式从热红外遥感影像中提取大田玉米冠层温度。

为了验证本发明基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法的提取效果，以及冠层温度对大田玉米水分胁迫状况的监测效果，本文利用手持红外测温仪雷泰(raytek)st60+采集的冠层温度数据以及利用植物叶片气孔计ap4测量的叶片气孔导度，分别对大田玉米冠层温度提取效果和水分胁迫状况监测效果进行验证。如图2所示，本实施例在试验田中五个灌溉梯度处理区分别均匀布设3个采样区，共15个采样区。每个样方布设3个数据采样点（红色圆圈），以三次测量的平均值作为采样点的数据测量，以三个采样点数据测量的平均值作为样方的数据测量。如图3所示，本发明方法提取的冠层温度与地面实测值具有高度的一致性（r²=0.94，rmse=0.70℃），可以有效的提取大田玉米冠层温度且具有较高精度。如图4所示，本发明方法提取的冠层温度与经典水分胁迫指数叶片气孔导度具有极显著（p<0.01）的相关性，决定系数r²为0.76，说明本方法所提冠层温度可以有效的监测大田玉米水分胁迫状况。

本实施例中基于无人机可见光和热红外遥感技术的提取的大田玉米水分胁迫状况分布图，如图5所示，可以通过冠层温度反映大田玉米水分胁迫状况的空间分布信息，为田间灌溉管理提供技术支持。

结合上述实施例可以说明本发明涉及的一种基于无人机可见光和热红外遥感的大田玉米冠层温度提取方法，可以有效提取大田玉米冠层温度且具有较高精度，同时通过大田玉米冠层温度分布图可以有效的监测大田玉米水分胁迫状况，解决了现有大田作物水分胁迫监测技术中费时、费力、成本较高，并且没有考虑土壤、作物等的空间变异性的问题，可以为田间灌溉管理提供技术支持，促进农业信息化。

上述实施例仅为本发明具体实施例，但并不局限于该实施例，凡在不脱离本发明构思的情况下，依据本申请所做的等效修饰和现有技术添加均视为本发明技术范畴。