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

本发明专利技术公开了一种基于绿度水分综合指数的大豆制图方法，包括步骤一，计算时间窗口；步骤二，计算尺度因子；步骤三，构建绿度水分综合指数；步骤四，最优阈值确定；步骤五，大豆分类制图；所述步骤一中，基于等间距时间节点集合，NDVI序列的变化率ρ的计算公式为所述步骤二中，整幅影像尺度因子sf的计算公式为sf＝2*(NDVI/SWTR)；本发明专利技术相较于现有的大豆制图方法，创新性地结合NDVI和SWIR波段构建绿度水分综合指数(GWCCI)，能够在大豆生长峰值时期显著增强大豆与其它地物类型的信息差异；本发明专利技术采用的基于GWCCI的大豆制图方法，可实现基于单景遥感影像的、端到端的、全自动的大豆制图，不需要训练样本且不依赖长时间序列遥感数据集NDVI。感数据集NDVI。感数据集NDVI。

【技术实现步骤摘要】
一种基于绿度水分综合指数的大豆制图方法


[0001]本专利技术涉及作物制图
，具体为一种基于绿度水分综合指数的大豆制图方法。

技术介绍

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

技术实现思路

[0003]本专利技术的目的在于提供一种基于绿度水分综合指数的大豆制图方法，以解决上述
技术介绍
中提出的问题。
[0004]为实现上述目的，本专利技术提供如下技术方案：一种基于绿度水分综合指数的大豆制图方法，包括以下步骤：步骤一，计算时间窗口；步骤二，计算尺度因子；步骤三，构建绿度水分综合指数；步骤四，最优阈值确定；步骤五，大豆分类制图；
[0005]其中在上述步骤一中，下载4月份到11月份的多光谱遥感影像，计算每幅影像的NDVI，基于SG滤波法获得拟合后的NDVI时序剖面图，基于拟合后的NDVI时序剖面图，设定时间间距k，基于大豆整个生长季获取n个等间距时间节点及其对应的NDVI值，计算NDVI序列的变化率ρ，将ρ在生长季内首次小于0.1的时间节点t1作为大豆生长峰值时期的起点，首次小于
‑
0.1的时间节点t2作为大豆生长峰值时期的终点，从而获得大豆生长峰值时期时间窗口(t1，t2)，即遥感影像数据获取的时间窗口；
[0006]其中在上述步骤二中，在设定的时间窗口(t1，t2)内任选一景无云的多光谱遥感影像M
t
，计算该影像的NDVI值并获取SWIR波段，由于此时间节点的NDVI值远大于SWIR波段的反射率，为保证NDVI值和SWIR处于同一可比较的水平，利用波段计算器计算整幅影像的尺度因子sf，选取部分代表性大豆样本计算相应尺度因子的平均值，将该值作为GWCCI的尺度因子sf；
[0007]其中在上述步骤三中，由于大豆生长峰值时期大豆与其它地物类型的冠层含水量
差异最大，因此采用NDVI与SWIR波段反射率乘积的方式构建绿度水分综合指数(GWCCI)：
[0008][0009]其中在上述步骤四中，根据公式(1)计算整幅遥感影像的GWCCI指数图，计算阈值区间，设定阈值步长，获得目标阈值集合，通过试错法，计算不同目标阈值下大豆的制图精度，将制图精度最高时对应的阈值确定为该指数的最优阈值；
[0010]其中在上述步骤五中，认为当目标像素数值大于目标阈值时该像素被分类为大豆，反之为非大豆，基于最优阈值完成大豆分类制图。
[0011]优选的，所述步骤一中，基于大豆整个生长季获取n个等间距时间节点(1，2，3...，i，...n)及其对应的NDVI值(NDVI1，NDVI2，...，NDVI
i
，...，NDVI
n
)。
[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，...，NDVI
i
，...，NDVI
n
)，基于等间距时间节点集合，NDVI序列的变化率ρ的计算公式为将ρ在生长季内首次小于0.1的时间节点t1作为大豆生长峰值时期的起点，首次小于
‑
0.1的时间节点t2作为大豆生长峰值时期的终点，从而获得大豆生长峰值时期时间窗口(t1，t2)，即遥感影像数据获取的时间窗口；
[0022]其中在上述步骤二中，在设定的时间窗口(t1，t2)内任选一景无云的多光谱遥感影像M
t
，计算该影像的NDVI值并获取SWIR波段，由于此时间节点的NDVI值远大于SWIR波段的反射率，为保证NDVI值和SWIR处于同一可比较的水平，利用波段计算器计算整本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于绿度水分综合指数的大豆制图方法，包括以下步骤：步骤一，计算时间窗口；步骤二，计算尺度因子；步骤三，构建绿度水分综合指数；步骤四，最优阈值确定；步骤五，大豆分类制图；其特征在于：其中在上述步骤一中，下载4月份到11月份的多光谱遥感影像，计算每幅影像的NDVI，基于SG滤波法获得拟合后的NDVI时序剖面图，基于拟合后的NDVI时序剖面图，设定时间间距k，基于大豆整个生长季获取n个等间距时间节点及其对应的NDVI值，计算NDVI序列的变化率ρ，将ρ在生长季内首次小于0.1的时间节点t1作为大豆生长峰值时期的起点，首次小于
‑
0.1的时间节点t2作为大豆生长峰值时期的终点，从而获得大豆生长峰值时期时间窗口(t1，t2)，即遥感影像数据获取的时间窗口；其中在上述步骤二中，在设定的时间窗口(t1，t2)内任选一景无云的多光谱遥感影像M
t
，计算该影像的NDVI值并获取SWIR波段，由于此时间节点的NDVI值远大于SWIR波段的反射率，为保证NDVI值和SWIR处于同一可比较的水平，利用波段计算器计算整幅影像的尺度因子sf，选取部分代表性大豆样本计算相应尺度因子的平均值，将该值作为GWCCI的尺度因子sf；其中在上述步骤三中，由于大豆生长峰值时期大豆与其它地物类型的冠层含水量差异最大，因此采用NDVI与SWIR波段反射率乘积的方式构建绿度水分综合指数(GWCCI)：其中在上述步骤四中，根据公式(1)计算整幅遥感影像的GWCCI指数图，计算阈值区间，设定阈值步长，获得目标阈值集合，通过试错法，计算不同目标阈值下大豆的制图精度，将制图精度...

【专利技术属性】
技术研发人员：陈慧，李华朋，张树清，
申请(专利权)人：中国科学院东北地理与农业生态研究所，
类型：发明
国别省市：