本发明涉及一种冬小麦叶绿素含量的预测系统及方法，其包括：预先设置的预测系统中的无人机按照设定的航线和高度飞行，由搭载在无人机上的多光谱相机对田间进行高空影像数据采集，同时采集冬小麦地里预先设定的小区内的冠层叶面积性状数据；通过影像携带的地理坐标信息对影像进行拼接，然后进行降噪及白板矫正，获取ROI的反射率；根据各个小区反射率，结合植被系数方程，计算出各个小区的不同植被指数的值；根据地表温度传感器获取到地表温度值，采集从种植当天到无人机的飞行数据的温度积累量；根据计算的不同植被指数，结合地表温度和测试的叶绿素表型数据，采用神经网络的算法构建叶绿素的估测模型，实现对冬小麦的叶绿素进行预测。
1.一种冬小麦叶绿素含量的预测方法，其特征在于，包括：
步骤S1、预先设置的预测系统中的无人机按照设定的航线和高度飞行，由搭载在无人机上的多光谱相机对田间进行高空影像数据采集，同时采集冬小麦地里预先设定的小区内的冠层叶面积性状数据；
步骤S2、通过影像携带的地理坐标信息对影像进行拼接，然后进行降噪及白板矫正，获取ROI的反射率；
步骤S3、根据各个小区反射率，结合植被系数方程，计算出各个小区的不同植被指数的值；
步骤S4、根据地表温度传感器获取到地表温度值，采集从种植当天到无人机的飞行数据的温度积累量；
步骤S5、根据计算的不同植被指数，结合地表温度和测试的叶绿素表型数据，采用神经网络的算法构建叶绿素的估测模型，实现对冬小麦的叶绿素进行预测。
2.如权利要求1所述预测方法，其特征在于，所述步骤S1中，多光谱相机包括5个通道。
3.如权利要求2所述预测方法，其特征在于，所述步骤S1中，打开多光谱相机的5个通道同时进行影像的拍摄，拍摄时镜头垂直向下；同时地面的冬小麦应按照种植图分别对每个小区进行地面叶绿素数据采集。
4.如权利要求1所述预测方法，其特征在于，所述步骤S2中，降噪采用椒盐噪声处理，获取平滑的影像。
5.如权利要求1所述预测方法，其特征在于，所述步骤S2中，利用ENVI掩膜处理获取各个小区的5个通道反射率；5个通道反射率分别为：红光反射率记做Rred，蓝光反射率记做Rblue，绿光反射率记做Rgreen，近红外光反射率记做Rnir，红边光谱反射率记做Rrededge。
6.如权利要求1所述预测方法，其特征在于，所述步骤S2中，采用白板进行校正从像素的DN值转换为反射率，反射率计算公式为：式中，R为该小区对应的光谱通道的反射率，DN多光谱为该小区所对应的多光谱单一通道的DN值，DN白板为该通道拍摄白板提取的DN值。
7.如权利要求1所述预测方法，其特征在于，所述步骤S3中，植被系数方程包括：式中，NDVI是归一化植被指数，GNDVI绿色归一化植被指数，EVI是增强型植被指数，NDRE是归一化差异红色边缘指数，NDRGI是归一化绿红差异指数。
8.如权利要求1所述预测方法，其特征在于，所述步骤S5中，叶绿素的估测模型构建方法包括以下步骤：步骤S51：根据计算的不同植被指数结合地面测试数据，采用随机森林的算法进行模型的构建并对冠层的叶面积指数进行估测；
步骤S52：计算出所有小区的植被指数，所有的植被指数和地面数据形成M*N的矩阵；M表示冬小麦地里所有的ROI小区，N表示包含地面数据和植被指数的列，其中小麦的地面叶面积属性表示分类属性；
步骤S53：将M*N的矩阵，先采用随机森林的算法选取贡献率最高的3个指数，然后随机选取70％的数据进行模型的建立，选取剩下的30％数据进行模型的验证。
9.一种冬小麦叶绿素含量的预测系统，其特征在于，该系统用于实现如权利要求1至8任一项所述预测方法，包括：地表温度传感器、处理器、无人机、多光谱传感器和白板；
所述地表温度传感器，设置在冬小麦种植区域内，用于将检测到的种植区域的地面温度传输至所述处理器；
所述无人机，用于搭载所述多光谱传感器，所述多光谱传感器通过所述无人机随时高空采集种植区域内的冬小麦影像数据；
所述白板，用于接收所述多光谱传感器传输至的影像数据，将图像的DN值转换为反射率，并传输至所述处理器；
所述处理器，根据接收到的数据得到叶绿素的估测模型，实现对冬小麦叶绿素含量的预测。
10.如权利要求9所述预测系统，其特征在于，所述地表温度传感器采用探针式的地表温度传感器，完全浸没在土壤中。
一种冬小麦叶绿素含量的预测系统及方法
技术领域
[0001] 本发明涉及一种智慧农业技术领域，特别是关于一种融合地表温度传感器和无人机多光谱遥感数据预测冬小麦叶绿素含量的系统及方法。
背景技术
[0002] 叶绿素是植物生长的必要色素，是农作物生长状况监测的重要指标。通过非破坏性的测量，研究叶绿素含量与冬小麦生长发育的关系，对于提高冬小麦质量，把握小麦在生长过程中的趋势具有重要参考价值。研究发现，作物叶片及冠层光谱指数与其叶绿素含量密切相关。
发明内容
[0003] 针对上述问题，本发明的目的是提供一种冬小麦叶绿素含量的预测系统及方法，其能得到不同生育期作物的叶绿素含量，获取小麦的表型，同时也实现了多传感融合的农业应用。
[0004] 为实现上述目的，本发明采取以下技术方案：一种冬小麦叶绿素含量的预测方法，其包括：步骤S1、预先设置的预测系统中的无人机按照设定的航线和高度飞行，由搭载在无人机上的多光谱相机对田间进行高空影像数据采集，同时采集冬小麦地里预先设定的小区内的冠层叶面积性状数据；步骤S2、通过影像携带的地理坐标信息对影像进行拼接，然后进行降噪及白板矫正，获取ROI的反射率；步骤S3、根据各个小区反射率，结合植被系数方程，计算出各个小区的不同植被指数的值；步骤S4、根据地表温度传感器获取到地表温度值，采集从种植当天到无人机的飞行数据的温度积累量；步骤S5、根据计算的不同植被指数，结合地表温度和测试的叶绿素表型数据，采用神经网络的算法构建叶绿素的估测模型，实现对冬小麦的叶绿素进行预测。
[0005] 进一步，所述步骤S1中，多光谱相机包括5个通道。
[0006] 进一步，所述步骤S1中，打开多光谱相机的5个通道同时进行影像的拍摄，拍摄时镜头垂直向下；同时地面的冬小麦应按照种植图分别对每个小区进行地面叶绿素数据采集。
[0007] 进一步，所述步骤S2中，降噪采用椒盐噪声处理，获取平滑的影像。
[0008] 进一步，所述步骤S2中，利用ENVI掩膜处理获取各个小区的5个通道反射率；5个通道反射率分别为：红光反射率记做Rred，蓝光反射率记做Rblue，绿光反射率记做Rgreen，近红外光反射率记做Rnir，红边光谱反射率记做Rrededge。
[0009] 进一步，所述步骤S2中，采用白板进行校正从像素的DN值转换为反射率，反射率计算公式为：
[0010]
[0011] 式中，R为该小区对应的光谱通道的反射率，DN多光谱为该小区所对应的多光谱单一通道的DN值，DN白板为该通道拍摄白板提取的DN值。
[0012] 进一步，所述步骤S3中，植被系数方程包括：
[0013]
[0014]
[0015]
[0016]
[0017]
[0018] 式中，NDVI是归一化植被指数，GNDVI绿色归一化植被指数，EVI是增强型植被指数，NDRE是归一化差异红色边缘指数，NDRGI是归一化绿红差异指数。
[0019] 进一步，所述步骤S5中，叶绿素的估测模型构建方法包括以下步骤：
[0020] 步骤S51：根据计算的不同植被指数结合地面测试数据，采用随机森林的算法进行模型的构建并对冠层的叶面积指数进行估测；
[0021] 步骤S52：计算出所有小区的植被指数，所有的植被指数和地面数据形成M*N的矩阵；M表示冬小麦地里所有的ROI小区，N表示包含地面数据和植被指数的列，其中小麦的地面叶面积属性表示分类属性；
[0022] 步骤S53：将M*N的矩阵，先采用随机森林的算法选取贡献率最高的3个指数，然后随机选取70％的数据进行模型的建立，选取剩下的30％数据进行模型的验证。
[0023] 一种冬小麦叶绿素含量的预测系统，该系统用于实现上述预测方法，包括：地表温度传感器、处理器、无人机、多光谱传感器和白板；
[0024] 所述地表温度传感器，设置在冬小麦种植区域内，用于将检测到的种植区域的地面温度传输至所述处理器；
[0025] 所述无人机，用于搭载所述多光谱传感器，所述多光谱传感器通过所述无人机随时高空采集种植区域内的冬小麦影像数据；
[0026] 所述白板，用于接收所述多光谱传感器传输至的影像数据，将图像的DN值转换为反射率，并传输至所述处理器；
[0027] 所述处理器，根据接收到的数据得到叶绿素的估测模型，实现对冬小麦叶绿素含量的预测。
[0028] 进一步，所述地表温度传感器采用探针式的地表温度传感器，完全浸没在土壤中。
[0029] 本发明由于采取以上技术方案，其具有以下优点：本发明采用的多光谱相机可以进行高精度，多时像的图像采集，并采用拼接算法对图像进行拼接，进行白板校正，把图像的DN值转换为反射率，通过反射率和光谱植被系数方程组运算得到多种植被指数。结合多种植被指数和地面的温度数据以及人工采集的冬小麦叶绿素表型数据，采用神经网络算法进行模型的生成，得到不同生育期作物的叶绿素含量模型，为小麦的表型获取提供了新的有效的技术手段，有效的提高了冬小麦质量。
附图说明
[0030] 图1是本发明实施例中的冬小麦叶绿素含量预测系统结构示意图；
[0031] 图2是本发明实施例中的冬小麦叶绿素含量预测方法流程示意图；
[0032] 图3是本发明实施例中的冬小麦多光谱影像拼接后的可见光ROI小区分布示意图。
具体实施方式
[0033] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。
[0034] 本发明提供一种融合地表温度传感器和无人机多光谱的冬小麦叶绿素含量的系统及方法，包含了地表的温度传感器以及四旋翼的无人机，其中该无人机一体式的多光谱成像系统，集成了1个可见光传感器及5个多光谱传感器(蓝光，绿光，红光，红边和近红外)。地表的温度传感器可以全天进行数据收集并保存。多光谱的相机可以进行高精度，多时像的图像采集，然后采用拼接算法对图像进行拼接，接着对拼接的图像数据进行白板校正，把图像的DN值转换为反射率，通过反射率和光谱植被系数方程组运算得到多种植被指数。结合多种植被指数和地面的温度数据以及人工采集的冬小麦叶绿素表型数据，采用神经网络算法进行模型的生成，找出不同生育期作物的叶绿素含量模型，为小麦的表型获取提供了新的有效的技术手段，同时也实现了多传感融合的农业应用。
[0035] 在本发明的第一实施方式中，如图1所示，提供一种冬小麦叶绿素含量的预测系统，其包括地表温度传感器、处理器、无人机、多光谱传感器和白板。
[0036] 地表温度传感器，设置在冬小麦种植区域内，用于将检测到的种植区域的地面温度传输至处理器；
[0037] 无人机，用于搭载多光谱传感器，多光谱传感器通过无人机随时高空采集种植区域内的冬小麦影像数据；在本实施例中优选四旋翼无人机；
[0038] 白板，用于接收多光谱传感器传输至的影像数据，将图像的DN值转换为反射率，并传输至处理器；
[0039] 处理器，根据接收到的数据得到叶绿素的估测模型，实现对冬小麦叶绿素含量的预测。
[0040] 上述实施例中，还包括可见光传感器。该可见光传感器搭载在无人机上，用于实时高空采集冬小麦的影像数据。
[0041] 上述实施例中，地表温度传感器优选在冬小麦种植以后按照种植图进行布局。在本实施例中，采用探针式的地表温度传感器，完全浸没在土壤中，该地表传感器可以实时采集地表温度并进行保存。
[0042] 上述实施例中，多光谱传感器为多光谱相机，该多光谱相机包括5个采集通道。在本实施例中，优选的，5个通道的中心波长分别为730nm，840nm，560nm，650nm，450nm。
[0043] 在本发明的第二实施方式中，提供一种冬小麦叶绿素含量的预测方法，该方法基于第一实施方式中的预测系统实现，在本实施例中，选择晴朗无云的天气，时间在上午的11点到下午的13点之间，进行冬小麦的监测。如图2所示，具体预测方法包括以下步骤：
[0044] 步骤S1、四旋翼无人机按照预先设置的航线和高度飞行，由多光谱相机对田间进行高空影像数据采集，同时采集冬小麦地里预先设定的小区内的叶绿素的性状数据；
[0045] 具体为：打开多光谱相机的5个通道同时进行影像的拍摄，拍摄时镜头垂直向下；同时地面的冬小麦应按照种植图分别对每个小区进行地面叶绿素数据采集；
[0046] 在本实施例中，预先设定的小区是根据地面冬小麦的种植图划分得到，采集的小区为ROI区域，如图3所示；
[0047] 步骤S2、通过影像携带的地理坐标信息对影像进行拼接，然后进行降噪及白板矫正，获取ROI的反射率；
[0048] 在本实施例中，降噪采用椒盐噪声处理，获取平滑的影像；利用ENVI掩膜处理获取各个小区的5个通道反射率；
[0049] 获取的5个通道反射率分别为：红光反射率记做Rred，蓝光反射率记做Rblue，绿光反射率记做Rgreen，近红外光反射率记做Rnir，红边光谱反射率记做Rrededge；
[0050] 步骤S3、根据各个小区反射率，结合植被系数方程，计算出各个小区的不同植被指数的值；
[0051] 步骤S4、根据地表温度传感器获取到地表温度值，采集从种植当天到无人机的飞行数据的温度积累量；
[0052] 步骤S5、根据计算的不同植被指数，结合地表温度和测试的叶绿素的性状数据，采用神经网络的算法构建叶绿素的估测模型，实现对冬小麦的叶绿素进行预测。
[0053] 上述各步骤中，本实施例选择上午的11点到下午的13点，这个时期是太阳直射，防止太阳倾斜造成误差，选择晴朗无云的天气是避免天气对航拍影响进行干扰。
[0054] 上述步骤S1中，无人机的飞行高度在10m～30m之间，飞行速度控制在1～4m/s,航向和旁向的重叠度均为70％。
[0055] 上述步骤S2中，采用白板进行校正从像素的DN值转换为反射率。具体的反射率计算公式如下：
[0056]
[0057] 式中，R为该小区对应的光谱通道的反射率，DN多光谱为该小区所对应的多光谱单一通道的DN值，DN白板为该通道拍摄白板提取的DN值。通过这样的方法依次转化5个通道对应小区的反射率。
[0058] 上述步骤S3中，植被系数方程包括以下等式：
[0059]
[0060]
[0061]
[0062]
[0063]
[0064] 式中，NDVI是归一化植被指数，GNDVI绿色归一化植被指数，EVI是增强型植被指数，NDRE是归一化差异红色边缘指数，NDRGI是归一化绿红差异指数。
[0065] 上述步骤S4中，温度积累量T的表达公式如下：
[0066]
[0067] 其中，种植下去的第一天i记做1，k是无人机飞行的时间距离种植第一天的天数，ti表示每天的平均温度。
[0068] 上述步骤S5中，叶绿素的估测模型构建方法包括以下步骤：
[0069] 步骤S51：根据计算的不同植被指数结合地面测试数据，采用随机森林的算法进行模型的构建并对冠层的叶面积指数进行估测；
[0070] 其中，植被指数由植被系数方程获得。
[0071] 步骤S52：计算出所有小区的植被指数，所有的植被指数和地面数据形成M*N的矩阵；M表示冬小麦地里所有的ROI小区，N表示包含地面数据和植被指数的列，其中小麦的地面叶面积属性表示分类属性；地面数据是指步骤S1中的每个小区所对应的叶绿素的性状数据；
[0072] 步骤S53：将M*N的矩阵，先采用随机森林的算法选取贡献率最高的3个指数，然后随机选取70％的数据进行模型的建立，选取剩下的30％数据进行模型的验证。
[0073] 本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD‑ROM、光学存储器等)上实施的计算机程序产品的形式。
[0074] 本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0075] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0076] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
