一种基于遥感的罂粟提取方法及系统[发明专利]

(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 (43)申请公布日 (21)申请号 201810049848.1
(22)申请日 2018.01.18
(71)申请人 中国科学院遥感与数字地球研究所
地址 100000 北京市朝阳区安定门外大屯
路甲20号北
(72)发明人 王思远　游永发　马元旭　申明　
刘卫华　
(74)专利代理机构 北京高沃律师事务所 11569
代理人 王戈
(51)Int.Cl.
G06T 7/00(2017.01)
G01C 11/00(2006.01)
(54)发明名称一种基于遥感的罂粟提取方法及系统(57)摘要本发明公开了一种基于遥感的罂粟提取方法及系统。

方法包括：首先，获取罂粟不同生长阶段的研究区域的多幅遥感影像，并计算多幅所述遥感影像的特征值，将所述特征值按照研究区域的罂粟生长阶段组成时间序列影像文件；然后，根据所述时间序列影像文件选取少量罂粟和研究区域的其它类型的背景植被的样本点，生成所述罂粟和其它各个类型的背景植被对应的时序变化曲线，并根据所述时序变化曲线建立罂粟提取模型；最后，设定罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果，实现了罂粟种植的有效监
测。

权利要求书3页 说明书8页 附图5页CN 108364284 A 2018.08.03
C N 108364284
A
1.一种基于遥感的罂粟提取方法，其特征在于，所述提取方法包括如下步骤：
获取罂粟不同生长阶段的研究区域的多幅遥感影像；
计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；
将所述特征值按照研究区域的罂粟生长阶段组成时间序列影像文件，所述时间序列影像文件包括归一化植被指数时间序列影像文件、增强型植被指数时间序列影像文件、近红外波段反射率时间序列影像文件和红波段反射率时间序列影像文件；
根据所述时间序列影像文件选取少量罂粟和研究区域的其它类型的背景植被的样本点，生成所述罂粟和其它各个类型的背景植被对应的时序变化曲线，所述时序变化曲线包括归一化植被指数时序变化曲线、增强型植被指数时序变化曲线、近红外波段反射率时序变化曲线和红波段反射率时序变化曲线；
根据所述时序变化曲线建立罂粟提取模型；
设定所述罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果。

2.根据权利要求1所述的一种基于遥感的罂粟提取方法，其特征在于，所述计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；具体包括：
对多幅所述遥感影像进行预处理，得到预处理后的遥感影像；
提取各所述预处理后的遥感影像的近红外波段反射率、红波段反射率和蓝波段反射率，获得所述遥感影像的近红外波段反射率和红波段反射率；
根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数。

3.根据权利要求2所述的一种基于遥感的罂粟提取方法，其特征在于，所述根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数，具体包括：
根据所述近红外波段反射率和红波段反射率，利用公式(1)计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
其中，ρNIR和ρRED分别表示近红外波段反射率和红波段反射率。

4.根据权利要求2所述的一种基于遥感的罂粟提取方法，其特征在于，所述根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数，具体包括：
根据所述近红外波段反射率、红波段反射率和蓝波段反射率，利用公式(2)计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数；
其中，ρNIR、ρRED和ρBLUE分别表示近红外波段反射率、红波段反射率和蓝波段反射率。

5.根据权利要求1所述的一种基于遥感的罂粟提取方法，其特征在于，所述设定所述罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果之后，还包括：
采用小斑消除算法对所述提取结果进行后处理，将零碎图斑去除，获得最终的提取结果。

6.一种基于遥感的罂粟提取系统，其特征在于，所述提取系统包括：
遥感影像获取模块，用于获取罂粟不同生长阶段的研究区域的多幅遥感影像；
特征值计算模块，用于计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；
时间序列影像建立模块，用于将所述特征值按照研究区域的罂粟生长阶段组成时间序列影像文件，所述时间序列影像文件包括归一化植被指数时间序列影像文件、增强型植被指数时间序列影像文件、近红外波段反射率时间序列影像文件和红波段反射率时间序列影像文件；
时序变化曲线建立模块，用于根据所述时间序列影像文件选取少量罂粟和研究区域的其它类型植被的样本点，生成所述罂粟和其它各个类型植被对应的时序变化曲线，所述时序变化曲线包括归一化植被指数时序变化曲线、增强型植被指数时序变化曲线、近红外波段反射率时序变化曲线和红波段反射率时序变化曲线；
提取模型建立模块，用于根据所述时序变化曲线建立罂粟提取模型；
提取结果获取模块，用于设定所述罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果。

7.根据权利要求6所述的一种基于遥感的罂粟提取系统，其特征在于，所述特征值计算模块，具体包括：
预处理子模块，用于对多幅所述遥感影像进行预处理，得到预处理后的遥感影像；
近红外波段反射率和红波段反射率获取子模块，用于提取各所述预处理后的遥感影像的近红外波段反射率、红波段反射率和蓝波段反射率，获得所述遥感影像的近红外波段反射率和红波段反射率；
归一化植被指数计算子模块，用于根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
增强型植被指数计算子模块，用于根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数。

8.根据权利要求7所述的一种基于遥感的罂粟提取系统，其特征在于，所述归一化植被指数计算子模块，具体包括：
归一化植被指数计算单元，用于根据所述近红外波段反射率和红波段反射率，利用公式(1)计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
其中，ρNIR和ρRED分别表示近红外波段的反射率和红波段的反射率。

9.根据权利要求7所述的一种基于遥感的罂粟提取系统，其特征在于，所述增强型植被指数计算子模块，具体包括：
增强型植被指数计算单元，用于根据所述近红外波段反射率、红波段反射率和蓝波段反射率，利用公式(2)计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数；
其中，ρNIR、ρRED和ρBLUE分别表示近红外波段反射率、红波段反射率和蓝波段反射率。

10.根据权利要求1所述的一种基于遥感的罂粟提取系统，其特征在于，所述提取系统还包括：
后处理模块，用于采用小斑消除算法对所述提取结果进行后处理，将零碎图斑去除，获得最终的提取结果。

一种基于遥感的罂粟提取方法及系统
技术领域
[0001]本发明涉及罂粟监测领域，特别涉及一种基于遥感的罂粟提取方法及系统。

背景技术
[0002]我国面临的禁毒形势严峻复杂，毒品犯罪多发高发，禁毒工作任务十分艰巨。

精确识别毒品原植物是禁毒工作的关键，而罂粟作为一种主要毒品原植物，对其进行识别至关重要。

非法罂粟种植的位置一般十分隐蔽，利用传统人工方式进行排查不仅费时费力，而且效率十分低下。

遥感技术的发展为大范围的罂粟种植监测提供了一条有效途径。

[0003]从遥感影像上提取罂粟的方法主要可分为监督分类与非监督分类两种。

基于监督分类方法提取罂粟的研究，如联合国毒品与犯罪问题办公室(The United Nations Office On Drugs and Crime)使用监督分类方法对高空间分辨率遥感影像进行分类，从而监测主要罂粟种植区域的罂粟种植情况。

该方法需要人工选取罂粟样本，其过程不仅费时费力，而且罂粟样本通常很难获得，因而限制了该方法的应用。

基于非监督分类方法提取罂粟的研究，如采用一种多基元混合像元分解的非监督分类方法从高光谱影像中识别罂粟，但该方法效率很低，针对该问题又提出了一种混合调谐匹配滤波的非监督分类方法从高光谱影像中检测罂粟。

现有的方法都是针对高光谱影像进行分析的，而高光谱影像虽然光谱分辨率很高，但其覆盖率很低同时幅宽小，因此采用高光谱影像实现大范围的罂粟检测也不太符合实际情况。

[0004]总体来看，以上这些方法都是基于单时相遥感影像进行分析的，未充分利用不同植被物候特征的差异进行分析，在有云的情况会导致影像成像质量较差，会导致分类结果出现很大误差，而且通常罂粟的生长期较短，而高光谱影像重返周期较长且地面覆盖率不高，会由于数据缺失导致该区域无法进行罂粟监测，因此现有的在遥感影像上提取罂粟的方法，无法实现罂粟种植的有效监测。

发明内容
[0005]本发明的目的是，为了实现罂粟种植的有效监测，提供一种基于遥感的罂粟提取方法及系统。

[0006]为实现上述目的，本发明提供了如下方案：
[0007]一种基于遥感的罂粟提取方法，所述提取方法包括如下步骤：
[0008]获取罂粟不同生长阶段的研究区域的多幅遥感影像；
[0009]计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；
[0010]将所述特征值按照研究区域的罂粟生长阶段组成时间序列影像文件，所述时间序列影像文件包括归一化植被指数时间序列影像文件、增强型植被指数时间序列影像文件、近红外波段反射率时间序列影像文件和红波段反射率时间序列影像文件；
[0011]根据所述时间序列影像文件选取少量罂粟和研究区域的其它类型植被的样本点，
生成所述罂粟和其它各个类型植被对应的时序变化曲线，所述时序变化曲线包括归一化植被指数时序变化曲线、增强型植被指数时序变化曲线、近红外波段反射率时序变化曲线和红波段反射率时序变化曲线；
[0012]根据所述时序变化曲线建立罂粟提取模型；
[0013]设定罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果。

[0014]可选的，所述计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；具体包括：
[0015]对多幅所述遥感影像进行预处理，得到预处理后的遥感影像；
[0016]提取各预处理后的遥感影像的近红外波段反射率、红波段反射率和蓝波段反射率，获得所述遥感影像的近红外波段反射率和红波段反射率；
[0017]根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
[0018]根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数。

[0019]可选的，所述根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数，具体包括：
[0020]根据所述近红外波段反射率和红波段反射率，利用公式(1)计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
[0021]
[0022]其中，ρNIR和ρRED分别表示近红外波段反射率和红波段反射率。

[0023]可选的，所述根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数，具体包括：[0024]根据所述近红外波段反射率、红波段反射率和蓝波段反射率，利用公式(2)计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数；
[0025]
[0026]其中，ρNIR、ρRED和ρBLUE分别表示近红外波段反射率、红波段反射率和蓝波段反射率。

[0027]可选的，所述设定罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果之后，还包括：
[0028]采用小斑消除算法对所述提取结果进行后处理，将零碎图斑去除，获得最终的提取结果。

[0029]一种基于遥感的罂粟提取系统，所述提取系统包括：
[0030]遥感影像获取模块，用于获取罂粟不同生长阶段的研究区域的多幅遥感影像；[0031]特征值计算模块，用于计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；
[0032]时间序列影像建立模块，用于将所述特征值按照研究区域的罂粟生长阶段组成时
间序列影像文件，所述时间序列影像文件包括归一化植被指数时间序列影像文件、增强型植被指数时间序列影像文件、近红外波段反射率时间序列影像文件和红波段反射率时间序列影像文件；
[0033]时序变化曲线建立模块，用于根据所述时间序列影像文件选取少量罂粟和研究区域的其它类型植被的样本点，生成所述罂粟和其它各个类型植被对应的时序变化曲线，所述时序变化曲线包括归一化植被指数时序变化曲线、增强型植被指数时序变化曲线、近红外波段反射率时序变化曲线和红波段反射率时序变化曲线；
[0034]提取模型建立模块，用于根据所述时序变化曲线建立罂粟提取模型；
[0035]提取结果获取模块，用于设定罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果。

[0036]可选的，所述特征值计算模块，具体包括：
[0037]预处理子模块，用于对多幅所述遥感影像进行预处理，得到预处理后的遥感影像；[0038]近红外波段反射率和红波段反射率获取子模块，用于提取各预处理后的遥感影像的近红外波段反射率、红波段反射率和蓝波段反射率，获得所述遥感影像的近红外波段反射率和红波段反射率；
[0039]归一化植被指数计算子模块，用于根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；[0040]增强型植被指数计算子模块，用于根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数。

[0041]可选的，所述归一化植被指数计算子模块，具体包括：
[0042]归一化植被指数计算单元，用于根据所述近红外波段反射率和红波段反射率，利用公式(1)计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
[0043]
[0044]其中，ρNIR和ρRED分别表示近红外波段反射率和红波段反射率。

[0045]可选的，所述增强型植被指数计算子模块，具体包括：
[0046]增强型植被指数计算单元，用于根据所述近红外波段反射率、红波段反射率和蓝波段反射率，利用公式(2)计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数；
[0047]
[0048]其中，ρNIR、ρRED和ρBLUE分别表示近红外波段反射率、红波段反射率和蓝波段反射率。

[0049]可选的，所述提取系统还包括：
[0050]后处理模块，用于采用小斑消除算法对所述提取结果进行后处理，将零碎图斑去除，获得最终的提取结果。

[0051]根据本发明提供的具体实施例，本发明公开了以下技术效果：
[0052]本发明公开了一种基于遥感的罂粟提取方法及系统，该方法，首先，生成时间序列影像文件，然后，通过选取少量的罂粟和研究区域的其它类型的背景植被的样本点，生成罂粟和每个类型的背景植被对应的时序变化曲线，并根据时序变化曲线建立罂粟提取模型，最后，根据所述罂粟提取模型，从所述时间序列影像文件中提取研究区域内种植的罂粟，获得提取结果。

本发明通过时间序列影像对罂粟与背景植被的物候及光谱变化特征进行分析，自动提取罂粟，实现了罂粟种植的有效监测。

附图说明
[0053]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

[0054]图1为本发明提供的一种基于遥感的罂粟提取方法的流程图；
[0055]图2为本发明提供的一种基于遥感的罂粟提取系统的结构框图；
[0056]图3为本发明提供的一个具体的实施例的基于遥感的罂粟提取方法的流程图；[0057]图4为本发明提供的一个具体的实施例的归一化植被指数时序变化曲线；[0058]图5为本发明提供的一个具体的实时例的增强型植被指数时序变化曲线；[0059]图6为本发明提供的一个具体的实时例的近红外波段反射率时序变化曲线；[0060]图7为本发明提供的一个具体的实施例的红波段反射率时序变化曲线。

具体实施方式
[0061]本发明的目的是，为了实现罂粟种植的有效监测，提供一种基于遥感的罂粟提取方法及系统。

[0062]为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对发明作进一步详细的说明。

[0063]如图1所示，本发明提供了一种基于遥感的罂粟提取方法，所述提取方法包括如下步骤：
[0064]步骤101，获取罂粟不同生长阶段的研究区域的多幅遥感影像。

[0065]步骤102，计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；具体包括：
[0066]对多幅所述遥感影像进行预处理，得到预处理后的遥感影像；提取各预处理后的遥感影像的近红外波段反射率、红波段反射率和蓝波段反射率，获得所述遥感影像的近红外波段反射率和红波段反射率；根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数。

[0067]所述根据所述近红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数，具体包括：
[0068]根据所述近红外波段反射率和红波段反射率，利用公式(1)计算各预处理后的遥
感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；
[0069]
[0070]其中，ρNIR和ρRED分别表示近红外波段反射率和红波段反射率。

[0071]所述根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数，具体包括：[0072]根据所述近红外波段反射率、红波段反射率和蓝波段反射率，利用公式(2)计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数；
[0073]
[0074]其中，ρNIR、ρRED和ρBLUE分别表示近红外波段反射率、红波段反射率和蓝波段反射率。

[0075]步骤103，将所述特征值按照研究区域的罂粟生长阶段组成时间序列影像文件，所述时间序列影像文件包括归一化植被指数时间序列影像文件、增强型植被指数时间序列影像文件、近红外波段反射率时间序列影像文件和红波段反射率时间序列影像文件。

[0076]步骤104，根据所述时间序列影像文件选取少量罂粟和研究区域的其它类型植被的样本点，生成所述罂粟和其它各个类型植被对应的时序变化曲线，所述时序变化曲线包括归一化植被指数时序变化曲线、增强型植被指数时序变化曲线、近红外波段反射率时序变化曲线和红波段反射率时序变化曲线。

[0077]步骤105，根据所述时序变化曲线建立罂粟提取模型。

[0078]步骤106，设定罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果。

[0079]步骤106所述设定罂粟提取模型的阈值，根据所述时间序列影像文件对研究区域内的罂粟进行提取，获得提取结果之后，还包括：采用小斑消除算法对所述提取结果进行后处理，将零碎图斑去除，获得最终的提取结果。

[0080]如图2所示，本发明还提供一种基于遥感的罂粟提取系统，所述提取系统包括：[0081]遥感影像获取模块201，用于获取罂粟不同生长阶段的研究区域的多幅遥感影像；[0082]特征值计算模块202，用于计算多幅所述遥感影像的特征值，所述特征值包括归一化植被指数、增强型植被指数、近红外波段反射率和红波段反射率；所述特征值计算模块202，具体包括：预处理子模块，用于对多幅所述遥感影像进行预处理，得到预处理后的遥感影像；近红外波段反射率和红波段反射率获取子模块，用于提取各预处理后的遥感影像的近红外波段反射率、红波段反射率和蓝波段反射率，获得所述遥感影像的近红外波段反射率和红波段反射率；归一化植被指数计算子模块，用于根据所述红外波段反射率和红波段反射率计算各预处理后的遥感影像的归一化植被指数，获得所述遥感影像的归一化植被指数；增强型植被指数计算子模块，用于根据所述近红外波段反射率、红波段反射率和蓝波段反射率计算各预处理后的遥感影像的增强型植被指数，获得所述遥感影像的增强型植被指数。

[0083]所述归一化植被指数计算子模块，具体包括：归一化植被指数计算单元，用于根据所述近红外波段反射率和红波段反射率，利用公式(1)计算各预处理后的遥感影像的归一。