高光谱监测方法、装置、设备及存储介质与流程

1.本技术涉及遥感影像处理技术领域，尤其是涉及一种高光谱监测方法、装置、设备及存储介质。


背景技术：

2.随着人类活动轨迹的发展和扩大，一些外来植物侵占了当地原有植物的栖息地，改变了当地的生态系统，对生态环境造成了一定程度的破坏。
3.相关技术中，采取中低分辨率时序数据分析法对入侵植物的分布趋势进行预测，然而，中低分辨率时序数据分析法的分辨率低，导致入侵植物的分布趋势预测的准确度不够，单个像元的尺寸太大，在一些小尺度区域内不能使用。


技术实现要素：

4.本技术旨在至少解决现有技术中存在的技术问题之一。为此，本技术提出一高光谱监测方法、装置、设备及存储介质，能够提高入侵植物区域识别的准确度。
5.为实现上述目的，本公开实施例的第一方面提出了一种高光谱监测方法，包括：
6.获取待检测区域内多个预设的时间周期的高光谱影像数据集；其中，所述高光谱影像数据集包括多个高光谱影像数据图，所述时间周期内的所述入侵植物具有突出的光谱特征；
7.对所述高光谱影像数据集中的每一高光谱影像数据图进行预处理，得到每一高光谱影像数据图的标准影像图；
8.获取预设的入侵植物的光谱曲线图；其中，所述光谱曲线图与所述时间周期内的入侵植物的所述光谱特征相匹配；
9.根据每一高光谱影像数据图的标准影像图和所述光谱曲线图得到分布图；其中，所述分布图用于表征所述待检测区域内每一时间周期内所述入侵植物的分布区域；
10.根据多个分布图评估所述待检测区域内的所述入侵植物的分布趋势。
11.根据本技术实施例的高光谱监测方法，至少具有如下有益效果：通过使用高光谱影像，提高了入侵植物区域识别的准确度，通过多张入侵植物的分布图，能够直观地呈现外来入侵植物的蔓延趋势，反映出入侵植物的蔓延方向，并对入侵植物的未来蔓延趋势做出预测。
12.根据本技术的一些实施例，
13.所述预处理包括：几何校正处理和大气校正处理；
14.所述对所述高光谱影像数据集中的每一高光谱影像数据图进行预处理，得到每一高光谱影像数据图的标准影像图，包括：
15.根据预设的第一校正算法对所述高光谱影像数据集中的每一高光谱影像数据图几何校正处理，得到几何标准图；
16.根据预设的第二校正算法对所述几何标准图进行大气校正处理，得到每一高光谱
影像数据图的标准影像图。
17.根据本技术的一些实施例，所述根据每一高光谱影像数据图的标准影像图和所述光谱曲线图得到分布图，具体包括：
18.对标准影像图中每一像元进行转换处理，得到每一像元的像元特征值；
19.对所述光谱曲线图进行转换处理，得到光谱特征值；
20.根据多个所述像元特征值和所述光谱特征值得到分布图。
21.根据本技术的一些实施例，所述根据多个所述像元特征值和光谱特征值得到分布图，具体包括：获取预设的像元比较阈值；
22.计算每一所述像元特征值和所述光谱特征值的差值，得到特征差值；
23.根据所述像元比较阈值和所述特征差值得到每一像元的比较结果；
24.根据多个所述比较结果得到分布图。
25.根据本技术的一些实施例，所述根据所述像元比较阈值和所述特征差值得到每一像元的比较结果，具体包括：
26.若所述特征差值小于所述像元比较阈值，则所述比较结果为所述像元相似；
27.若所述特征差值大于或等于所述像元比较阈值，则所述比较结果为所述像元不相似。
28.根据本技术的一些实施例，所述方法还包括：
29.记录获取待检测区域内高光谱影像数据图的时间顺序。
30.根据本技术的一些实施例，所述根据多个分布图评估所述待检测区域内的所述入侵植物的分布趋势，具体包括：
31.根据所述时间顺序对多个所述分布图进行排序处理，得到变化趋势图集；
32.对所述变化趋势图集进行动态化处理，得到动态变化图；
33.根据动态变化图评估所述待检测区域内的所述入侵植物的分布趋势。
34.为实现上述目的，本公开实施例的第二方面提出了一种高光谱监测装置，包括：
35.获取模块，用于获取待检测区域内多个预设的时间周期的高光谱影像数据集；其中，所述高光谱影像数据集包括多个高光谱影像数据图，所述时间周期内所述入侵植物具有突出的光谱特征；
36.预处理模块，用于对所述高光谱影像数据集中的每一高光谱影像数据图进行预处理，得到每一高光谱影像数据图的标准影像图；
37.入侵植物光谱获取模块，用于获取预设的入侵植物的光谱曲线图，其中，所述光谱曲线图与所述时间周期内的入侵植物的所述光谱特征相匹配；
38.分布图生成模块，用于根据每一高光谱影像数据图的标准影像图和所述光谱曲线图得到分布图；其中，所述分布图用于表征所述待检测区域内每一时间周期内所述入侵植物的分布区域；
39.评估模块，用于根据多个分布图评估所述待检测区域内的所述入侵植物的分布趋势。
40.为实现上述目的，本公开实施例的第三方面提出了一种计算机设备，包括：
41.至少一个存储器；
42.至少一个处理器；
43.至少一个程序；
44.所述程序被存储在所述存储器中，处理器执行所述至少一个程序以实现本公开实施例第一方面所述的方法。
45.为实现上述目的，本公开实施例的第四方面提出了一种存储介质，该存储介质是计算机可读存储介质，所述计算机可读存储介质存储有计算机可执行指令，所述计算机可执行指令用于使计算机执行：
46.本公开实施例第一方面所述的方法。
47.本技术的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本技术的实践了解到。
附图说明
48.下面结合附图和实施例对本技术做进一步的说明，其中：
49.图1是本技术实施例提供的高光谱监测方法的流程图；
50.图2是图1中所示步骤s200的流程图；
51.图3是图1中所示步骤s400的流程图；
52.图4是图3中所示步骤s430的流程图；
53.图5是图1中所示步骤s500的流程图；
54.图6是本技术实施例提供的高光谱监测装置的模块框图；
55.图7是本技术实施例提供的电子设备的硬件结构示意图。
具体实施方式
56.下面详细描述本技术的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
57.在本技术的描述中，需要理解的是，涉及到方位描述，例如上、下、前、后、左、右等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本技术和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本技术的限制。
58.在本技术的描述中，若干的含义是一个或者多个，多个的含义是两个以上，大于、小于、超过等理解为不包括本数，以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。
59.本技术的描述中，除非另有明确的限定，设置、安装、连接等词语应做广义理解，所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本技术中的具体含义。
60.随着人类活动轨迹的发展和扩大，一些外来植物侵占了当地原有植物的栖息地，改变了当地的生态系统，对生态环境造成了一定程度的破坏。
61.相关技术中，采取中低分辨率时序数据分析法对入侵植物的分布趋势进行预测，然而，中低分辨率时序数据分析法的准确度不够，在一些小尺度区域内不能使用。
62.在外来入侵植物的监测方面，对于入侵植物的扩展活动规律呈现仍然比较抽象，不能直观的显示外来入侵植物的蔓延趋势；对于外来入侵植物的蔓延趋势以理论预测为主，只是对外来入侵植物的入侵概率进行估测，对物种或区域环境依赖程度高，普适性差。
63.基于此，本技术提供了一种高光谱监测方法，通过利用高光谱影像数据，提高了对外来入侵植物辨别的准确性。
64.下面结合附图，对本技术实施例的高光谱监测方法作进一步阐述。
65.请参照图1，本技术的一些实施例提供了一种高光谱监测方法，该方法包括但不限于步骤s100、步骤s200、步骤s300、步骤s400、步骤s500，下面对这五个步骤进行详细介绍。
66.步骤s100，获取待检测区域内多个预设的时间周期的高光谱影像数据集；其中，高光谱影像数据集包括多个高光谱影像数据图，预设的时间周期内的入侵植物具有突出的光谱特征，便于后续对高光谱影像数据图的分析和处理。
67.在步骤s100中，可以通过无人机获取高光谱影像数据；也可以通过卫星，还可以直接互联网上的相关资源，对此，本技术不作具体限制。
68.步骤s200，对高光谱影像数据集中的每一高光谱影像数据图进行预处理，得到每一高光谱影像数据图的标准影像图。
69.步骤s300，获取预设的入侵植物的光谱曲线图，其中，光谱曲线图与预设时间周期内的入侵植物的光谱特征相匹配。
70.需要说明的是，光谱特征指入侵植物的光谱曲线与地物背景光谱曲线差值最大的波段，在差值最大的波段上，可以实现入侵植物和背景的区分。
71.步骤s400，根据每一高光谱影像数据图的标准影像图和光谱曲线图得到分布图；其中，分布图用于表征所述待检测区域内每一时间周期内入侵植物的分布区域。
72.步骤s500，根据多个分布图评估所述待检测区域内的所述入侵植物的分布趋势。
73.本技术实施例的高光谱监测方法，通过在入侵植物具有突出的光谱特征的时间周期内，获取待检测区域的高光谱影像数据，得到高光谱影像数据集，然后对高光谱影像数据集进行预处理，得到每一高光谱影像数据图的标准影像图，再通过入侵植物的光谱特征生成分布图，根据多张分布图量化计算得出入侵植物的分布趋势，从而能够准确地评估入侵植物的蔓延过程以及未来的蔓延趋势。
74.参照图2，在一些实施例中，步骤s200包括但不限于步骤s210和步骤s220。
75.步骤s210：根据预设的第一校正算法对高光谱影像数据集中的每一高光谱影像数据图几何校正处理，得到几何标准图。
76.需要说明的是，在步骤s210中，第一校正算法主要用到二元n次多项式变换坐标的方法，具体变换公式如下所示：
[0077][0078]
其中，在公式(1)中，x，y为变换前图像坐标，u，v为变换后图像坐标，a
ij
，b
ij
为多项式系数，n与几何畸变的复杂程度密切相关，畸变程度越高，n取值越高，一般情况下n取值小于等于3。经过几何校正得到几何标准图，能够消除图像中的几何畸变。
[0079]
步骤s220，根据预设的第二校正算法对所述几何标准图进行大气校正处理，得到每一高光谱影像数据图的标准影像图。
[0080]
需要说明的是，在步骤s220中，第二校正算法主要用到黑暗像元法，黑暗像元法是一种大气校正处理用到的算法，其具体方法如下：
[0081]
使用envi(the environment for visualizing images，可视化图像处理软件)找到几何标准图中的黑暗像元并确定其各个波段反射率，将其他像元各个波段的反射率减去黑暗像元对应的各个波段的反射率，得到标准的像元的反射率，根据标准的像元所成图像即为标准影像图。
[0082]
需要说明的是，黑暗像元指反射率或辐射亮度很小的像元，其判定阈值可在envi中自由设定。
[0083]
通过预处理获得标准影像图，消除了图像中的几何畸变和大气畸变，提高了处理结果的精确度。
[0084]
参照图3，图3是步骤s400的细化步骤流程图，在一些实施例中，步骤s400包括但不限于以下步骤：步骤s410、步骤s420、步骤s430。
[0085]
步骤s410，对标准影像图中每一像元进行转换处理，得到每一像元的像元特征值。
[0086]
获取每一像元在特定波段上的反射率，该反射率即为像元的像元特征值。
[0087]
具体地，在该特定波段上，入侵植物具有突出的光谱特征。
[0088]
步骤s420，对所述光谱曲线图进行转换处理，得到光谱特征值。
[0089]
具体地，获取入侵植物在特定波段上的反射率，该反射率即为入侵植物的光谱特征值。
[0090]
步骤s430，根据多个所述像元特征值和所述光谱特征值得到分布图。
[0091]
将入侵植物在特定波段的光谱特征作为识别手段，并量化其数值，实现了对入侵植物的精确判断。
[0092]
参照图4，图4是步骤s430的细化步骤流程图，在一些实施例中，步骤s430包括但不限于以下步骤：步骤s431、步骤s432、步骤s433、步骤s434。
[0093]
步骤s431，获取预设的像元比较阈值。
[0094]
需要说明的是，像元比较阈值是在入侵植物与目标区域像元的光谱曲线的差值最大波段上，可接受的能将像元识别为入侵植物的最大差值，需要人为设定，以用于区别入侵植物与其他区域。
[0095]
步骤s432，计算每一所述像元特征值和所述光谱特征值的差值，得到特征差值。
[0096]
需要说明的是，像元特征值和光谱特征值均为光谱反射率的值，光谱反射率用于表示某波段该像元的光波吸收情况，可以反映出地物的颜色信息，特征差值为像元特征值和光谱特征值的差值的绝对值。
[0097]
步骤s433，根据所述像元比较阈值和所述特征差值得到每一像元的比较结果。
[0098]
具体地，比较过程包括：若所述特征差值小于所述像元比较阈值，则所述比较结果为像元相似，认定为入侵区域；若所述特征差值大于或等于所述像元比较阈值，则所述比较结果为像元不相似，认定为非入侵区域。
[0099]
步骤s434，根据多个所述比较结果得到分布图。
[0100]
具体地，在同一个标准影像图中采取两种不同的颜色分别表示入侵区域和非入侵
区域，得到分布图。入侵区域选用较为醒目的颜色如红色、紫色，非入侵区域选用较为平淡的颜色如绿色、蓝色。
[0101]
在本技术的一些实施例中，高光谱监测方法还包括：
[0102]
记录获取待检测区域内高光谱影像数据图的日期和时间顺序。
[0103]
具体地，在本实施例中，将高光谱影像数据图命名为获取日期加时间先后序号的格式，如某高光谱影像数据图是2020年1月1日获取的1号图片，则命名为1-20200101，方便后续的使用和整理。
[0104]
参照图5，图5是步骤s500的细化步骤流程图，在一些实施例中，步骤s500包括但不限于以下步骤：步骤s510、步骤s520、步骤s530。
[0105]
步骤s510：根据时间顺序对多个分布图进行排序处理，得到变化趋势图集；
[0106]
步骤s520：对变化趋势图集进行动态化处理，得到动态变化图；
[0107]
步骤s530：根据动态变化图评估待检测区域内的入侵植物的分布趋势。
[0108]
具体地，变化趋势图集是按时间先后顺序排列好的分布图，使用常用的图像处理软件，如pr(adobe premiere pro，视频编辑软件)，经过补帧得到动态变化图，动态变化图是短视频或者gif(graphics interchange format，图形交换格式)图的形式，动态的展现出待检测区域在预设时间周期内的入侵植物蔓延情况。
[0109]
具体地，根据动态变化图评估待检测区域内的入侵植物的分布趋势包括：
[0110]
使用谷歌地图应用接口提供的方法计算按时间顺序最早的分布图内入侵区域的中心点，计算每个固定角度入侵区域边缘点到中心点的距离，用该距离的变化值和高光谱影像数据图成像间的时间间隔，得到固定角度的入侵植物蔓延速度。根据蔓延速度即可定量地预测入侵植物未来的蔓延趋势。
[0111]
本公开实施例还提供一种高光谱监测评估装置，参照图6，该装置包括：第一获取模块100、预处理模块200、第二获取模块300、分布图生成模块400、评估模块500。
[0112]
第一获取模块100，用于获取待检测区域内多个预设的时间周期的高光谱影像数据集；
[0113]
预处理模块200，用于对所述高光谱影像数据集中的每一高光谱影像数据图进行预处理；
[0114]
第二获取模块300，用于获取预设的入侵植物的光谱曲线图；
[0115]
分布图生成模块400，用于根据每一高光谱影像数据图的标准影像图和所述光谱曲线图得到分布图；
[0116]
评估模块500，用于根据多个分布图评估所述待检测区域内的所述入侵植物的分布趋势。
[0117]
本公开实施例的高光谱监测装置，通过在入侵植物具有突出的光谱特征的时间周期内，获取待检测区域的高光谱影像数据，得到高光谱影像数据集，然后对高光谱影像数据集进行预处理，得到每一高光谱影像数据图的标准影像图，再通过入侵植物的光谱特征生成分布图，根据多张分布图量化计算得出入侵植物的分布趋势，从而能够准确地评估入侵植物的蔓延过程以及未来的蔓延趋势。
[0118]
需要说明的是，本实施例的高光谱监测装置的具体实施方式与上述高光谱监测方法的具体实施方式基本一致，在此不再赘述。
[0119]
本公开实施例还提供了一种计算机设备，包括：
[0120]
至少一个存储器；
[0121]
至少一个处理器；
[0122]
至少一个程序；
[0123]
所述程序被存储在所述存储器中，处理器执行所述至少一个程序以实现：本公开实施例上述的高光谱监测方法。该计算机设备可以为包括手机、平板电脑、个人数字助理、车载电脑等任意智能终端。
[0124]
请参阅图7，图7示意了另一实施例的电子设备的硬件结构，电子设备包括：
[0125]
处理器701，可以采用通用的cpu(central processing unit，中央处理器)、微处理器、应用专用集成电路(application specific integrated circuit，asic)、或者一个或多个集成电路等方式实现，用于执行相关程序，以实现本公开实施例所提供的技术方案；
[0126]
存储器702，可以采用rom(read only memory，只读存储器)、静态存储设备、动态存储设备或者ram(random access memory，随机存取存储器)等形式实现。存储器702可以存储操作系统和其他应用程序，在通过软件或者固件来实现本说明书实施例所提供的技术方案时，相关的程序代码保存在存储器702中，并由处理器701来调用执行本公开实施例的评估方法；
[0127]
输入/输出接口703，用于实现信息输入及输出；
[0128]
通信接口704，用于实现本设备与其他设备的通信交互，可以通过有线方式(例如usb、网线等)实现通信，也可以通过无线方式(例如移动网络、wifi、蓝牙等)实现通信；和
[0129]
总线705，在设备的各个组件(例如处理器701、存储器702、输入/输出接口703和通信接口704)之间传输信息；
[0130]
其中处理器701、存储器702、输入/输出接口703和通信接口704通过总线705实现彼此之间在设备内部的通信连接。
[0131]
本公开实施例还提供了一种存储介质，该存储介质是计算机可读存储介质，该计算机可读存储介质存储有计算机可执行指令，该计算机可执行指令用于使计算机执行上述高光谱监测方法。
[0132]
上面结合附图对本技术实施例作了详细说明，但是本技术不限于上述实施例，在所述技术领域普通技术人员所具备的知识范围内，还可以在不脱离本技术宗旨的前提下作出各种变化。
[0133]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0134]
尽管已经示出和描述了本技术的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本技术的范围由权利要求及其等同物限定。