一种基于微型光谱仪观测数据提取农作物叶面积指数和叶绿素含量的方法

1.本发明涉及反演技术领域，更具体地，涉及一种基于微型光谱仪连续观测提取农作物叶面积指数和叶绿素含量的方法。
2.

背景技术：

3.获取每日连续农作物参数，对于农田的智能化管理十分重要。农作物叶面积指数、叶绿素含量是其中最为重要的两个参数，他们与农田施肥密切相关。目前，可以使用手持仪器或者搭载在多种遥感平台上的传感器测量叶面积指数和叶绿素含量。但是，手持仪器需要耗费大量的人力，难以获取每日连续的数据。基于无人机、卫星等遥感平台的传感器，则会受到云、雨等天气状况的影响，难以获取连续数据。搭载在地面物联网iot上的传感器则可以突破天气状况的影响，可以获取每日连续的观测数据。目前常见的搭载在iot上的传感器主要是rgb相机，可以使用简单的孔隙率原理计算出叶面积指数，然而计算叶绿素含量需要用更多的波段，仅从rgb图像中很难准确计算出叶绿素含量，进而影响农作物氮含量的计算，影像施肥决策。因此，亟需发展可以从地面iot搭载的传感器上同步提取叶面积指数和叶绿素含量的算法。另外，当前的算法主要是从一天内某一时刻获取的数据提取参数，并未考虑一天之内冠层的辐射和反射率变化，会对叶面积指数和叶绿素含量的计算造成一定误差。
4.

技术实现要素：

5.针对背景技术中的问题，本发明提出一种基于微型光谱仪连续观测提取农作物叶面积指数和叶绿素含量的方法，包括：s1，基于多年历史总par观测数据和iot上的par传感器的实测入射par数据，估算iot实时par的散射比例；s2，针对s1得到的iot实时par的散射比例，建立每个波段的散射比例和par的散射比例的多元回归模型，使用所述多元回归模型估算微型光谱仪每个波段的散射比例；s3，耦合冠层辐射传输模型和线性半经验核驱动模型，估算每天的绿色叶片叶面积指数、叶片叶绿素含量和冠层叶绿素含量。
6.本发明还提出一种基于微型光谱仪观测数据提取农作物叶面积指数和叶绿素含量的系统，其特征在于，包括处理器，所述处理器能够实现如前所述的方法。
7.在本发明中，使用微型光谱仪，在不使用辐射定标板的情况下，提取相对辐射数据进行参数的反演，突破了以往必须使用辐射定标板计算绝对反射率进行反演的局限。在参数反演过程中，考虑了一天内不同天气情况下所有的太阳入射辐射，即直射辐射和散射辐射，而不是假设散射辐射为0。本发明利用气象站提供的历史多年的太阳总辐射、散射辐射
数据，创新性地使用反向传播神经网络训练学习方法，输入iot实测的太阳总辐射，反演计算得到实时的散射辐射，进而计算出散射比例。
8.本发明可以同时获取每日叶面积指数和叶绿素含量，不仅时间连续，而且精度高。更具体的创新点还包括：1) 使用微型光谱仪，在不使用辐射定标板的情况下，提取相对辐射数据进行参数的反演，突破了以往必须使用辐射定标板计算绝对反射率进行反演的局限。
9.2) 参数反演过程中，考虑了一天内不同天气情况下所有的太阳入射辐射，即直射辐射和散射辐射，而不是假设散射辐射为0。
10.3) 利用气象站提供的历史多年的太阳总辐射、散射辐射数据，创新性地使用反向传播神经网络训练学习方法，输入iot实测的太阳总辐射，反演计算得到实时的散射辐射，进而计算出散射比例。
11.4) 改进已有的冠层辐射传输模型，建立考虑散射比例的辐射传输模型，并计算相对反射率。
12.5) 耦合冠层辐射传输模型和线性半经验核驱动模型，将模型训练的相对反射率和地面观测的相对辐射数据，分别由线性半经验核驱动模型表达，通过反演权重系数的方式，考虑一天内的反射率和辐射变化。
13.6)基于耦合后的模型考虑每日的太阳辐射变化，建立反向传播神经网络，实现每日连续叶面积指数和叶绿素含量的同步估算，为农田智慧管理提供必要的输入数据。
14.附图说明
15.为了更容易理解本发明，将通过参照附图中示出的具体实施方式更详细地描述本发明。这些附图只描绘了本发明的典型实施方式，不应认为对本发明保护范围的限制。
16.图1为本发明的方法的一个实施例的流程图。
17.图2为本发明的方法的另一个实施例的流程图。
18.图3为本发明的方法的实验效果图。
19.具体实施方式
20.下面参照附图描述本发明的实施方式，以便于本领域的技术人员可以更好的理解本发明并能予以实施，但所列举的实施例不作为本发明的限定，在不冲突的情况下，下述的实施例及实施例中的技术特征可以相互组合，其中相同的部件用相同的附图标记表示。
21.在本发明中，采用了地面物联网iot设备。iot上搭载有：微型光谱仪和测量par（光合有效辐射）的par传感器。
22.如图1-2所示，本发明的方法包括s1-s5。
23.s1，基于多年历史总par观测数据和iot上par传感器的实测入射par数据，估算iot实时par的散射比例。更具体地，步骤s1包括s11-s12。
24.s11，太阳入射辐射分为直射par和由云等散射引起的散射par两部分。通过地面气象观测站获取多年总辐射par和散射par，作为输入数据来训练第一反向传播神经网络模
型。优选地，还将如下数据输入到所述第一反向传播神经网络模型来对该模型进行训练：太阳天顶角和日地距离，多年历史par的散射比例。
25.反向传播神经网络模型是遥感领域常用的神经网络模型之一，适用于非线性问题的处理，具有功能稳定、结构原理易于理解、精度较高和容错率强的优势，因而本发明选择其进行参数反演。
26.s12，以地面iot（物联网）上par传感器实测的入射par数据作为输入，通过已训练的第一反向传播神经网络模型反演得到散射par数据，通过散射par数据和所述入射par数据计算iot实时par的散射比例。即，散射比例=散射par/入射par。
27.优选地，还将太阳天顶角和日地距离输入到以训练的第一反向传播神经网络模型进行反演。
28.s2，针对s1得到的实时par的散射比例，建立每个波段的散射比例和par的散射比例的多元回归模型，使用所述多元回归模型估算iot的微型光谱仪每个波段的散射比例。更具体地，步骤s2包括s21-s23。
29.s21，使用大气辐射传输模型（优选使用6s大气辐射传输模型），模拟计算iot上搭载的微型光谱仪的每个波段的散射比例，以及对应的par的散射比例。
30.s22，建立每个波段的散射比例和par的散射比例的一元线性或非线性回归模型。
31.具体来说，对每个波段的散射比例，以par的散射比例为自变量x，以该波段的散射比例为因变量y，使用最小二乘法建立二者间的正规方程，求解方程组系数得到回归方程的表达式。
32.通常是分别建立一元线性回归方程（y=ax b）、一元非线性回归方程（y=anxn an−1xn−1 ...... a0x0，其中n代表阶数，大于等于2），进行拟合，将拟合后的估计的y值与实际的因变量y进行对比，计算二者间的相关系数，取其中相关系数最大的那个回归方程作为最终的方程。
33.如下为得到的一组一元回归模型关系实例，针对某研究区，经过测试，确定了与之间的一元五次关系式（如下表所示）。
34.s23，将此多元回归关系运用到步骤s12得到的实测par的散射比例上，进而计算得到实时的每个波段的散射比例。
35.s3，耦合冠层辐射传输模型和线性半经验核驱动模型，估算每天的绿色叶片叶面积指数gai、叶片叶绿素含量lcc和冠层叶绿素含量ccc。更具体地，步骤s3包括s31-s35。可用的冠层辐射传输模型包括了一维冠层辐射传输模型prosail、acrm和slc等。也可以使用三维冠层辐射传输模型，例如dart、raytran等。
36.s31，按照模型的输入参数（例如冠层结构、生理参数与太阳、观测角度），自行设定每个参数的输入范围及步长，如叶面积指数设为0到7，步长是1，叶倾角设为30
°‑
60
°
，步长为10
°
等，输入参数间自由组合，生成一个训练数据集，将每一组训练数据集输入到冠层辐射传输模型，模拟出iot微型光谱仪对应波段的反射率。进一步，以多波段的反射率均值为分母，分别计算每个波段的相对反射率。
37.其中，代表波段上的冠层辐射传输模型模拟反射率，或者iot微型光谱仪实测的辐射，为太阳天顶角，为观测天顶角，为相对方位角，即太阳方位角减去观测方位角。代表相对值，即冠层辐射传输模型模拟的相对反射率，或者iot微型光谱仪实测的相对辐射值。
38.s32，基于s31中得到的每日反射率，建立第一核驱动模型，并得到第一核驱动权重系数。
39.1）同时模拟一天内的多组数据，对每组数据给定一个par波段的散射比例，建立第一核驱动模型：其中，代表由线性半经验核驱动模型模拟得到的数据，和分别为各向同性散射、体散射与几何光学散射对应的权重系数。分别为体散射核与几何光学核。
40.2）对每一天的模拟数据，对第一核驱动模型使用优化算法，计算核驱动模型中的权重系数：其中，j代表优化算法的非线性函数关系表达式。
41.s33，以s32模拟得到的每个波段的第一核驱动模型权重系数为输入，以对应的gai、lcc和ccc为输出，分别训练第二反向传播神经网络模型。
42.s34，对iot地面实测的数据，采用s31-s32同样的方法，计算每日核驱动模型权重系数。具体包括：1）基于iot微型光谱仪每个波段散射比例，建立第二核驱动模型。第二核驱动模型与第一核驱动模型类型是相同的，区别是可以根据需要k值不同。
43.2）基于iot微型光谱仪实测的辐射数据，对第二核驱动模型使用优化算法，计算每日核驱动模型权重系数。
44.s35，以s34得到的iot实测的每日核驱动权重系数为输入，使用s33训练后的第二反向传播神经网络模型进行反演，分别得到每日gai、lcc和ccc。
45.其中，第二反向传播神经网络模型可以与第一反向传播神经网络模型相同或不同。
46.本发明耦合了冠层辐射传输模型和线性半经验核驱动模型，将辐射传输模型训练的相对反射率和地面观测的相对辐射数据，分别由线性半经验核驱动模型表达，通过反演权重系数的方式，考虑一天内的反射率和辐射变化。基于耦合后的模型考虑每日的太阳辐射变化，建立反向传播神经网络，实现每日连续叶面积指数和叶绿素含量的同步估算，为农田智慧管理提供必要的输入数据。
47.图3中横坐标为地面验证数据，纵坐标为iot反演得到的结果，三个子图分别代表了gai、lcc和ccc的iot反演结果与地面验证数据的对比。可以看出：经过地面验证数据集的验证，本发明可以得到高精度、连续的gai、lcc和ccc。
48.以上所述的实施例，只是本发明较优选的具体实施方式，本说明书使用词组“在一种实施例中”、“在另一个实施例中”、“在又一实施例中”或“在其他实施例中”，其均可指代根据本公开的相同或不同实施例中的一个或多个。本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。