一种基于绿度水分综合指数的大豆制图方法

1.本发明涉及作物制图技术领域，具体为一种基于绿度水分综合指数的大豆制图方法。


背景技术：

2.基于遥感数据的农作物分类依赖于农作物独有的物候特征和光谱特征，由于种植系统的多样性以及不同国家地区气候环境等自然条件的巨大差异，许多农作物的物候和光谱特征存在交叠，导致单种作物分类精度不高，大豆作为重要的粮食和经济作物，通常与玉米水稻等粮食作物同时期间作种植且光谱特征较为相似，在进行大豆识别时极易与其它作物(尤其玉米)发生混淆，长时序遥感影像数据能够尽可能多地捕捉大豆的可区分特征，是当下进行大豆分类制图的主流数据源，机器学习等传统方法虽能获得较高大豆制图精度，但需要采集大量训练样本构建模型，耗费大量人力、物力和财力，方法难以在大范围应用推广；基于阈值的方法虽不需要采集样本，但需依赖长时序遥感影像开展作物制图，受天气等条件的干扰，高质量长时序遥感影像数据集往往难以获得，极大程度上限制了基于阈值的方法的应用，实际上，大豆制图的主要难点在于利用原始光谱信息识别大豆具有较大的局限性，如何最大化大豆的可区分特征，放大大豆与其它地物类型的信息差异(包括物候和光谱信息差异)是进行大豆高精度分类的主要难点所在。


技术实现要素：

3.本发明的目的在于提供一种基于绿度水分综合指数的大豆制图方法，以解决上述背景技术中提出的问题。
4.为实现上述目的，本发明提供如下技术方案：一种基于绿度水分综合指数的大豆制图方法，包括以下步骤：步骤一，计算时间窗口；步骤二，计算尺度因子；步骤三，构建绿度水分综合指数；步骤四，最优阈值确定；步骤五，大豆分类制图；
5.其中在上述步骤一中，下载4月份到11月份的多光谱遥感影像，计算每幅影像的ndvi，基于sg滤波法获得拟合后的ndvi时序剖面图，基于拟合后的ndvi时序剖面图，设定时间间距k，基于大豆整个生长季获取n个等间距时间节点及其对应的ndvi值，计算ndvi序列的变化率ρ，将ρ在生长季内首次小于0.1的时间节点t1作为大豆生长峰值时期的起点，首次小于-0.1的时间节点t2作为大豆生长峰值时期的终点，从而获得大豆生长峰值时期时间窗口(t1，t2)，即遥感影像数据获取的时间窗口；
6.其中在上述步骤二中，在设定的时间窗口(t1，t2)内任选一景无云的多光谱遥感影像m
t
，计算该影像的ndvi值并获取swir波段，由于此时间节点的ndvi值远大于swir波段的反射率，为保证ndvi值和swir处于同一可比较的水平，利用波段计算器计算整幅影像的尺度因子sf，选取部分代表性大豆样本计算相应尺度因子的平均值，将该值作为gwcci的尺度因子sf；
7.其中在上述步骤三中，由于大豆生长峰值时期大豆与其它地物类型的冠层含水量
差异最大，因此采用ndvi与swir波段反射率乘积的方式构建绿度水分综合指数(gwcci)：
[0008][0009]
其中在上述步骤四中，根据公式(1)计算整幅遥感影像的gwcci指数图，计算阈值区间，设定阈值步长，获得目标阈值集合，通过试错法，计算不同目标阈值下大豆的制图精度，将制图精度最高时对应的阈值确定为该指数的最优阈值；
[0010]
其中在上述步骤五中，认为当目标像素数值大于目标阈值时该像素被分类为大豆，反之为非大豆，基于最优阈值完成大豆分类制图。
[0011]
优选的，所述步骤一中，基于大豆整个生长季获取n个等间距时间节点(1，2，3...，i，...n)及其对应的ndvi值(ndvi1，ndvi2，...，ndvii，...，ndvin)。
[0012]
优选的，所述步骤一中，基于等间距时间节点集合，ndvi序列的变化率ρ的计算公式为
[0013]
优选的，所述步骤二中，整幅影像尺度因子sf的计算公式为sf＝2*(ndvi/swir))。
[0014]
优选的，所述步骤三中，t为获取的遥感影像的时间节点，即一年中的第t天(day of year，doy)，doy∈[0，365]，ndvi
t
为该时间点影像的ndvi值，为该时间点的swir波段的反射率，sf为上述计算得到的尺度因子。
[0015]
优选的，所述步骤四中，计算获得gwcci指数图后，在(1，2)阈值区间内设置步长(σ＝0.1)，获得目标阈值集合{1，1+σ，1+2σ，...，1+j*σ，...，2}，通过试错法，计算不同阈值下大豆的制图精度，将制图精度最高时对应的阈值确定为该指数的最优阈值。
[0016]
与现有技术相比，本发明的有益效果是：本发明相较于现有的大豆制图方法，创新性地结合ndvi和swir波段构建绿度水分综合指数(gwcci)，能够在大豆生长峰值时期显著增强大豆与其它地物类型的信息差异；本发明采用的基于gwcci的大豆制图方法，可实现基于单景遥感影像的、端到端的、全自动的大豆制图，不需要训练样本且不依赖长时间序列遥感数据集。
附图说明
[0017]
图1为本发明的步骤流程图；
[0018]
图2为本发明的gwcci方法流程图。
具体实施方式
[0019]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0020]
请参阅图1-2，本发明提供的一种实施例：一种基于绿度水分综合指数的大豆制图方法包括以下步骤：步骤一，计算时间窗口；步骤二，计算尺度因子；步骤三，构建绿度水分综合指数；步骤四，最优阈值确定；步骤五，大豆分类制图；
[0021]
其中在上述步骤一中，下载4月份到11月份的多光谱遥感影像，计算每幅影像的
ndvi，基于sg滤波法获得拟合后的ndvi时序剖面图，基于拟合后的ndvi时序剖面图，设定时间间距k，基于大豆整个生长季获取n个等间距时间节点(1，2，3...，i，...n)及其对应的ndvi值(ndvi1，ndvi2，...，ndvii，...，ndvin)，基于等间距时间节点集合，ndvi序列的变化率ρ的计算公式为将ρ在生长季内首次小于0.1的时间节点t1作为大豆生长峰值时期的起点，首次小于-0.1的时间节点t2作为大豆生长峰值时期的终点，从而获得大豆生长峰值时期时间窗口(t1，t2)，即遥感影像数据获取的时间窗口；
[0022]
其中在上述步骤二中，在设定的时间窗口(t1，t2)内任选一景无云的多光谱遥感影像m
t
，计算该影像的ndvi值并获取swir波段，由于此时间节点的ndvi值远大于swir波段的反射率，为保证ndvi值和swir处于同一可比较的水平，利用波段计算器计算整幅影像的尺度因子sf，选取部分代表性大豆样本计算相应尺度因子的平均值，将该值作为gwcci的尺度因子sf，其中，整幅影像尺度因子sf的计算公式为sf＝2*(ndvi/swir)；
[0023]
其中在上述步骤三中，由于大豆生长峰值时期大豆与其它地物类型的冠层含水量差异最大，因此采用ndvi与swir波段反射率乘积的方式构建绿度水分综合指数(gwcci)：
[0024][0025]
t为获取的遥感影像的时间节点，即一年中的第t天(day of year，doy)，doy∈[0，365]，ndvi
t
为该时间点影像的ndvi值，为该时间点的swir波段的反射率，sf为上述计算得到的尺度因子；
[0026]
其中在上述步骤四中，根据公式(1)计算整幅遥感影像的gwcci指数图，在(1，2)阈值区间内设置步长(σ＝0.1)，获得目标阈值集合{1，1+σ，1+2σ，...，1+j*σ，...，2}，通过试错法，计算不同阈值下大豆的制图精度，将制图精度最高时对应的阈值确定为该指数的最优阈值；
[0027]
其中在上述步骤五中，认为当目标像素数值大于目标阈值时该像素被分类为大豆，反之为非大豆，基于最优阈值完成大豆分类制图。
[0028]
基于上述，本发明的优点在于，该发明使用时，不需要依赖大量长时序遥感影像，规避了复杂且耗时的影像预处理过程，只需要获取时间窗口内任意一景多光谱遥感影像即可实现大豆的精确制图；第二，是一种端到端的、全自动的大豆制图方法，不需要训练样本构建模型，计算简单且方便快捷，特别适合于大尺度大豆制图应用需求。
[0029]
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。