首页 / 专利库 / 人工智能 / 机器学习 / 集成学习 / 随机森林 / 引入地形因子的黑土区大豆生物量多时相遥感反演方法

引入地形因子的黑土区大豆生物量多时相遥感反演方法

阅读:1008发布:2020-10-18

专利汇可以提供引入地形因子的黑土区大豆生物量多时相遥感反演方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了引入地形因子的黑土区大豆 生物 量多时相遥感反演方法,包括如下步骤:步骤一,地形特征信息获取;步骤二,生物量样本数据获取;步骤三,获取植被指数;步骤四,分段建模;步骤五,建立反演模型;步骤六,生物量预测;步骤七,生成图层;其中在上述的步骤一中,实地测量裸土期内 地 块 高 精度 DEM数据,使用 定位 精度为厘米级的海星达iRTK2对研究区地块进行实地测量;本发明,在传统植被指数 基础 上,创新性的加入地形相关因子,不仅提高了传统模型的精度和 稳定性 ,而且也能够对作物生长多时期反演估算;同时,解决了以往单纯用植被指数模型精度不高,稳定性不强的问题,也解决了生物量反演 时空 不连续的问题。,下面是引入地形因子的黑土区大豆生物量多时相遥感反演方法专利的具体信息内容。

1.引入地形因子的黑土区大豆生物量多时相遥感反演方法,包括如下步骤:步骤一,地形特征信息获取;步骤二,生物量样本数据获取;步骤三,获取植被指数;步骤四,分段建模;
步骤五,建立反演模型;步骤六,生物量预测;步骤七,生成图层;其特征在于:
其中在上述的步骤一中,实地测量裸土期内精度DEM数据,使用定位精度为厘米级的海星达iRTK2对研究区地块进行实地测量;将得到的780个精确坐标和高程点,在ArcGIS中生成TIN图层,并转成高精度DEM栅格数据,并提取地形特征信息;
其中在上述的步骤二中,分时段采集大豆整个生长期内,多期生物量样本数据,在大豆生长期批次均匀采样覆盖整个样区,每个样点面积为2.2m×2m;统计每个样点大豆株数,并在其中随机选取若干株长势均匀的植物样本,齐根剪取;在实验室中,将植物样品于80℃恒温烘干,获得植物地上部分干生物量;用以下公式:Wj=wjk×n/k,将其换算成采样点1m×
1m区域内大豆的地上干生物量;
其中在上述的步骤三中,订购研究区地块大豆生长各时期高精度遥感影像,进行正射校正和地形校正,提取植被指数;
其中在上述的步骤四中,根据作物生长不同时期状况,提取对应样点地形和对应时期的植被指数信息,进行分段建模;
其中在上述的步骤五中,在SPSS中,随机选取各时期样本数目的三分之二建模,三分之一验证;模型建立需要将植被指数、地形特征与作物生物量的相关性进行分析,选取相关系数最高的植被指数和海拔、坡度、坡向因子为备选估测因子,分别建立单一植被指数模型和引入地形因子的多元回归、神经网络、随机森林模型;模型构建前,为避免输入变量之间的相关性的影响,减少变维系数,先对上述变量进行主成分提取,建立基于主成分的大豆地上干生物量预测反演模型;
其中在上述的步骤六中,将其余样本带入各时期反演模型进行生物量预测,通过预测精度及均方根误差进行精度评价,可将样本数据分为为训练样本和验证样本,计算模型的决定系数R2和均方根误差RMSE,R2接近1,表示模型的相关性越好,拟合程度越高,方程越可靠;RMSE越小,表示方程的稳定性越高,结果越可靠;计算获得最佳模型为多元回归模型,试验区与验证区R2分别=0.91和0.89;RMSE=75.59和93.99,计算公式为:
其中,作物长势及产量是由作物种类、田间管理、气候土壤质地等多种因素综合作用的结果,考虑到田块尺度统一的耕种和施肥管理措施,气候在田块内部又相差甚微,地形导致的土壤分、温度、有机质等空间差异,就成为主导作物生长的重要因素;加入适当地形因子,有助于整体上提高模型的精度和稳定性;
其中在上述的步骤七中,选择模型精度和验证精度最高的模型作为各时期生物量最优反演模型,带入ArcGIS中生成图层。
2.根据权利要求1所述的引入地形因子的黑土区大豆生物量多时相遥感反演方法,其特征在于:所述步骤一中,高精度DEM栅格数据的分辨率为4米。
3.根据权利要求1所述的引入地形因子的黑土区大豆生物量多时相遥感反演方法,其特征在于:所述步骤二中,Wj(kg/m 2)表示第j个样单位面积内大豆地上干生物量;wjk表示第j个样单位面积内随机割取的大豆总干质量;n(株)为单位面积内大豆总株数;k(株)为样区内被割取称质量的大豆株数。
4.根据权利要求1所述的引入地形因子的黑土区大豆生物量多时相遥感反演方法,其特征在于:所述步骤三中,高精度遥感影像为SPOT-6影像,多光谱波段空间分辨率6m。
5.根据权利要求1所述的引入地形因子的黑土区大豆生物量多时相遥感反演方法,其特征在于:所述步骤六中,yi为生物量预测值,y为生物量样本值, 为实际样本值的平均值,n为样本个数。
6.根据权利要求1所述的引入地形因子的黑土区大豆生物量多时相遥感反演方法,其特征在于:所述步骤三和步骤四中,植被指数包括RVI比值植被指数、NDVI归一化差异植被指数、GNDVI绿波段归一化植被指数、DVI差值植被指数和EVI增强型植被指数。

说明书全文

引入地形因子的黑土区大豆生物量多时相遥感反演方法

技术领域

[0001] 本发明涉及大豆生物量多时期精准反演与监测技术领域,具体为引入地 形因子的黑土区大豆生物量多时相遥感反演方法。

背景技术

[0002] 生物量是作物长势监测的重要指标,单位面积生物量反映了作物群体长 势,是作物产量估算的基础。生物量估算是服务现代农业的一项重要内容, 及时准确的生物量模拟为田尺度作物精准管理提供技术支持,对国家农业 决策、农田生产管理、粮食仓储安全等都有重要意义。
[0003] 传统地面调查监测的统计模型与物理模型难以实用化,无论是从时间还 是从空间度来获取生物量,都有一定难度和局限;随着遥感技术的发展, 较高空间分辨率遥感数据可以准确捕获农田土壤与作物信息,能更为快速、 准确、高效、无破坏地对生物量进行估算。目前基于传统遥感植被指数的估 算模型方法因素单一,精度稳定性低,且无法对作物生长过程连续性表征, 难以满足田块尺度精准管理的要求,因此,设计引入地形因子的黑土区大豆生 物量多时相遥感反演方法是很有必要的。

发明内容

[0004] 本发明的目的在于提供引入地形因子的黑土区大豆生物量多时相遥感反 演方法,以解决上述背景技术中提出的问题。
[0005] 为实现上述目的,本发明提供如下技术方案:
[0006] 引入地形因子的黑土区大豆生物量多时相遥感反演方法,包括如下步骤: 步骤一,地形特征信息获取;步骤二,生物量样本数据获取;步骤三,获取 植被指数;步骤四,分段建模;步骤五,建立反演模型;步骤六,生物量预 测;步骤七,生成图层;
[0007] 其中在上述的步骤一中,实地测量裸土期内地块高精度DEM数据,使用 定位精度为厘米级的海星达iRTK2对研究区地块进行实地测量;将得到的780 个精确坐标和高程点,在ArcGIS中生成TIN图层,并转成高精度DEM栅格数 据,并提取地形特征信息;
[0008] 其中在上述的步骤二中,分时段采集大豆整个生长期内,多期生物量样 本数据,在大豆生长期批次均匀采样覆盖整个样区,每个样点面积为 2.2m×2m;统计每个样点大豆株数,并在其中随机选取若干株长势均匀的植 物样本,齐根剪取;在实验室中,将植物样品于80℃恒温烘干,获得植物 地上部分干生物量;用以下公式:Wj=wjk×n/k,将其换算成采样点1m×1 m区域内大豆的地上干生物量;
[0009] 其中在上述的步骤三中,订购研究区地块大豆生长各时期高精度遥感影 像,进行正射校正和地形校正,提取植被指数;
[0010] 其中在上述的步骤四中,根据作物生长不同时期状况,提取对应样点地 形和对应时期的植被指数信息,进行分段建模;
[0011] 其中在上述的步骤五中,在SPSS中,随机选取各时期样本数目的三分之 二建模,三分之一验证;模型建立需要将植被指数、地形特征与作物生物量 的相关性进行分析,选取相关系数最高的植被指数和海拔、坡度、坡向因子 为备选估测因子,分别建立单一植被指数模型和引入地形因子的多元回归、 神经网络、随机森林模型;模型构建前,为避免输入变量之间的相关性的影 响,减少变维系数,先对上述变量进行主成分提取,建立基于主成分的大豆 地上干生物量预测反演模型;
[0012] 其中在上述的步骤六中,将其余样本带入各时期反演模型进行生物量预 测,通过预测精度及均方根误差进行精度评价,可将样本数据分为为训练样 本和验证样本,计算模2 2
型的决定系数R和均方根误差RMSE,R接近1,表 示模型的相关性越好,拟合程度越高,方程越可靠;RMSE越小,表示方程的 稳定性越高,结果越可靠;计算获得最佳模型为多元回归模型,试验区与验 证区R2分别=0.91和0.89;RMSE=75.59和93.99,计算公式为:
[0013]其中,作物长势及产量是 由作物种类、田间管理、气候、土壤质地等多种因素综合作用的结果,考虑 到田块尺度统一的耕种和施肥管理措施,气候在田块内部又相差甚微,地形 导致的土壤分、温度、有机质等空间差异,就成为主导作物生长的重要因 素;加入适当地形因子,有助于整体上提高模型的精度和稳定性;
[0014] 其中在上述的步骤七中,选择模型精度和验证精度最高的模型作为各时 期生物量最优反演模型,带入ArcGIS中生成图层。
[0015] 根据上述技术方案,所述步骤一中,高精度DEM栅格数据的分辨率为4 米。
[0016] 根据上述技术方案,所述步骤二中,Wj(kg/m 2)表示第j个样单位面 积内大豆地上干生物量;wjk表示第j个样单位面积内随机割取的大豆总干 质量;n(株)为单位面积内大豆总株数;k(株)为样区内被割取称质量的大豆 株数。
[0017] 根据上述技术方案,所述步骤三中,高精度遥感影像为SPOT-6影像,多 光谱波段空间分辨率6m。
[0018] 根据上述技术方案,所述步骤六中,yi为生物量预测值,y为生物量样 本值,为实际样本值的平均值,n为样本个数。
[0019] 根据上述技术方案,所述步骤三和步骤四中,植被指数包括RVI比值植 被指数、NDVI归一化差异植被指数、GNDVI绿波段归一化植被指数、DVI差 值植被指数和EVI增强型植被指数。
[0020] 与现有技术相比,本发明的有益效果是:本发明,在传统植被指数基础 上,创新性的加入地形相关因子,提高了传统模型的精度和稳定性,并能够 对作物生长多时期反演估算;解决了以往单纯用植被指数模型精度不高,稳 定性不强的问题,同时也解决了生物量反演时空不连续的问题;本发明对田 块尺度精准管理及生物量遥感估算及监测具有实际意义。附图说明
[0021] 图1是本发明的方法流程图
[0022] 图2是生物量6月份遥感反演图;
[0023] 图3是生物量7月份遥感反演图;
[0024] 图4是生物量8月份遥感反演图;

具体实施方式

[0025] 下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行 清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而 不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做 出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
[0026] 请参阅图1-4,本发明提供一种技术方案:
[0027] 引入地形因子的黑土区大豆生物量多时相遥感反演方法,包括如下步骤: 步骤一,地形特征信息获取;步骤二,生物量样本数据获取;步骤三,获取 植被指数;步骤四,分段建模;步骤五,建立反演模型;步骤六,生物量预 测;步骤七,生成图层;
[0028] 其中在上述的步骤一中,实地测量裸土期内地块高精度DEM数据,使用 定位精度为厘米级的海星达iRTK2对研究区地块进行实地测量;将得到的780 个精确坐标和高程点,在ArcGIS中生成TIN图层,并转成分辨率为4米的高 精度DEM栅格数据,在ArcGIS中提取海拔、坡向、坡度、剖面曲率,利用 ERDAS 9.2建模,计算坡长、地形湿度等地形因子;
[0029] 其中在上述的步骤二中,在大豆生长期内,从6月15日开始到8月13 日大豆生长期内,分3批次均匀采样60个,覆盖整个样区,每个样点面积为 2.2m×2m;统计每个样点大豆株数,并在其中随机选取若干株长势均匀的植 物样本,齐根剪取;在实验室中,将植物样品于80℃恒温烘干,获得植物 地上部分干生物量;用以下公式:Wj=wjk×n/k,将其换算成采样点1m ×1m区域内大豆的地上干生物量;
[0030] 其中在上述的步骤三中,编程订购2016-6-9、2016-7-18、2016-8-8、3 期SPOT-6影像,多光谱波段空间分辨率6m;利用ENVI5.1对影像进行辐射 定标、大气校正、正射校正的处理,在ArcGIS 10.1中对影像进行裁剪;然 后对裁剪后的影像进行C模型地形校正;
[0031] 其中在上述的步骤四中,基于国内外植被指数研究状况和本研究区实地 状况,选择已构建的被认为有潜在预测作物生物量能的植被指数,根据作 物生长不同时期状况,提取对应样点地形和对应时期的植被指数信息,进行 分段建模;
[0032] 其中在上述的步骤五中,在SPSS中,根据作物生长状态不同将60个样 本分为两部分,7、8月份一起建模,6月份单独建模;随机选取各部分三分 之二样点用来建模,三分之一验证;以地形参数和植被指数参数作为输入量, 建立不同时期生物量反演模型;
[0033] 其中在上述的步骤六中,将其余样本带入各时期反演模型进行生物量预 测,通过预测精度及均方根误差进行精度评价,可将样本数据分为为训练样 本和验证样本,计算模型的决定系数R2和均方根误差RMSE,R2接近1,表 示模型的相关性越好,拟合程度越高,方程越可靠;RMSE越小,表示方程的 稳定性越高,结果越可靠;计算获得最佳模型为多元回归模型,试验区与验 证区R2分别=0.91和0.89;RMSE=75.59和93.99,计算公式为:
[0034]
[0035] 其中在上述的步骤七中,选择模型精度和验证精度最高的模型作为各时 期生物量最优反演模型,选择多元回归模型,在ArcGIS中反演出田块内大豆 生物量年内多时期空间分布情况,该生物量反演图能够准确反演大豆生长期 内时空格局连续性变化,为田间精准管理提供指导。
[0036] 根据上述技术方案,步骤一中,高精度DEM栅格数据的分辨率为4米。
[0037] 根据上述技术方案,步骤二中,Wj(kg/m 2)表示第j个样单位面积内 大豆地上干生物量;wjk表示第j个样单位面积内随机割取的大豆总干质量; n(株)为单位面积内大豆总株数;k(株)为样区内被割取称质量的大豆株数。
[0038] 根据上述技术方案,步骤三中,高精度遥感影像为SPOT-6影像,多光谱 波段空间分辨率6m。
[0039] 根据上述技术方案,步骤六中,yi为生物量预测值,y为生物量样本值,  为实际样本值的平均值,n为样本个数。
[0040] 根据上述技术方案,步骤三和步骤四中,植被指数包括RVI比值植被指 数、NDVI归一化差异植被指数、GNDVI绿波段归一化植被指数、DVI差值植 被指数和EVI增强型植被指数。
[0041] 表1为大豆生物量模型预测与验证表:
[0042]
[0043] 表2为大豆生物量反演模型结果:
[0044]
[0045] 基于上述,本发明的优点在于,本发明,实地测量裸土期内地块高精度 DEM数据,使用定位精度为厘米级的海星达iRTK2对研究区地块进行实地测 量;将得到的780个精确坐标和高程点,在ArcGIS中生成TIN图层,并转成 分辨率为4米的高精度DEM栅格数据,在ArcGIS中提取海拔、坡向、坡度、 剖面曲率,利用ERDAS 9.2建模,计算坡长、地形湿度等地形因子;在大豆 生长期内,从6月15日开始到8月13日大豆生长期内,分3批次均匀采样 60个,覆盖整个样区,每个样点面积为2.2m×2m;统计每个样点大豆株数, 并在其中随机选取若干株长势均匀的植物样本,齐根剪取;在实验室中,将 植物样品于80℃恒温烘干,获得植物地上部分干生物量;用以下公式: Wj=wjk×n/k,Wj(kg/m 2)表示第j个样单位面积内大豆地上干生物量; wjk表示第j个样单位面积内随机割取的大豆总干质量;n(株)为单位面积内 大豆总株数;k(株)为样区内被割取称质量的大豆株数,将其换算成采样点1 m×1m区域内大豆的地上干生物量;编程订购2016-6-9、2016-7-18、 2016-8-8,3期SPOT-6影像,多光谱波段空间分辨率6m;利用ENVI5.1对影 像进行辐射定标、大气校正、正射校正的处理,在ArcGIS 10.1中对影像进 行裁剪;然后对裁剪后的影像进行C模型地形校正;基于国内外植被指数研 究状况和本研究区实地状况,选择已构建的被认为有潜在预测作物生物量能 力的植被指数,植被指数包括RVI比值植被指数、NDVI归一化差异植被指数、 GNDVI绿波段归一化植被指数、DVI差值植被指数和EVI增强型植被指数,根 据作物生长不同时期状况,提取对应样点地形和对应时期的植被指数信息, 进行分段建模;在SPSS中,根据作物生长状态不同将60个样本分为两部分, 7、8月份一起建模,6月份单独建模;随机选取各部分三分之二样点用来建 模,三分之一验证;以地形参数和植被指数参数作为输入量,建立不同时期 生物量反演模型;将其余样本带入各时期反演模型进行生物量预测,通过预 测精度及均方根误差进行精度评价,可将样本数据分为为训练样本和验证样 本,计算模型的决定系2 2
数R 和均方根误差RMSE,R 接近1,表示模型的相 关性越好,拟合程度越高,方程越可靠;
RMSE越小,表示方程的稳定性越高, 结果越可靠;计算获得最佳模型为多元回归模型,试验区与验证区R2分别 =0.91和0.89;RMSE=75.59和93.99,计算公式为:
[0046]其中,yi为生物量预测值, y为生物量样本值, 为实际样本值的平均值,n为样本个数,作物长势及 产量是由作物种类、田间管理、气候、土壤质地等多种因素综合作用的结果, 考虑到田块尺度统一的耕种和施肥管理措施,气候在田块内部又相差甚微, 地形导致的土壤水分、温度、有机质等空间差异,就成为主导作物生长的重 要因素;加入适当地形因子,有助于整体上提高模型的精度和稳定性;选择 模型精度和验证精度最高的模型作为各时期生物量最优反演模型,带入 ArcGIS中生成图层,选择多元回归模型,在ArcGIS中反演出田块内大豆生 物量年内多时期空间分布情况,该生物量反演图能够准确反演大豆生长期内 时空格局连续性变化,为田间精准管理提供指导。
[0047] 需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来 将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示 这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、 “包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系 列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明 确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有 的要素。
[0048] 尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而 言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行 多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限 定。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈