一种基于机器学习回归算法的土壤盐度预测方法技术

本发明专利技术公开了一种基于机器学习回归算法的土壤盐度预测方法，包括如下步骤：S1、以等边三角形为测点形状，采用EM38大地电导仪采集三角形的三个角点的实测盐度，测量时，读取垂直读数EMV和水平读数EMH，取平均值作为三角形区域的盐度，测量若干个三角形区域的EMH和EMV为基础真理训练集(TS)；S2、利用L波段雷达后向散射系数去除植被贡献来获得土壤组分；S3、通过直接栅格化将实测地块转换为栅格来创建训练集(TS)；S4、利用上述训练集(TS)，采用随机森林回归(RFR)算法建模应用于联合数据集进行盐度预测。

【技术实现步骤摘要】
一种基于机器学习回归算法的土壤盐度预测方法
本专利技术涉及土壤盐度检测领域，具体涉及一种基于机器学习回归算法的土壤盐度预测方法。
技术介绍
土壤盐渍化是农作物生长发育受抑制的重要土壤障碍因子之一。土壤盐渍化影响作物生产和粮食安全。绘制盐渍度的空间分布和土壤盐渍化严重程度对农业管理和发展至关重要。
技术实现思路
为解决上述问题，本专利技术提供了一种基于机器学习回归算法的土壤盐度预测方法。为实现上述目的，本专利技术采取的技术方案为：一种基于机器学习回归算法的土壤盐度预测方法，包括如下步骤：S1、以等边三角形为测点形状，采用EM38大地电导仪采集三角形的三个角点的实测盐度，测量时，读取垂直读数EMV和水平读数EMH，取平均值作为三角形区域的盐度，测量若干个三角形区域的EMH和EMV为基础真理训练集(TS)；S2、利用L波段雷达后向散射系数去除植被贡献来获得土壤组分；S3、通过直接栅格化将实测地块转换为栅格来创建训练集(TS)；S4、利用上述训练集(TS)，采用随机森林回归(RFR)算法建模应用于联合数据集进行盐度预测。进一步地，所述步骤S1中，所述等边三角形的边长15-20米，确保三角形可以近似表示一个TM(专题成像仪)像素；测量区块的规格1m×1m。进一步地，所述步骤S2具体包括如下步骤：S21、用无线电测量校准从欧洲航天局获取的LANDSAT5(陆地卫星5)TM图像，并采用散射法模型消除附加的大气影响，将所产生的反射率每一波段重新设定为0-1；S22、对从欧洲航天局获取的1.5级雷达结果进行几何校正，并将像素重新调整到12.5米的大小，对水平发射水平接收HH和水平发射垂直接收HV的数字量化值(DN)分别进行校正，按照公式转换为后向散射系数(σ0HH、σ0Hv)；然后应用增强Lee滤波器(EnhancedLeefilter3×3大小，Lee1980)去除斑点或噪声，导出σ0HH和σ0Hv，并重放大到30m像素用来匹配TM数据；S23、构建植被含水量(VWC)对雷达后向散射系数影响的水云模型：L2＝exp(-2BV2sec(θi))其中，σ0是来自植被冠层和土壤的总后向散射系数(σ0HH、σ0Hv)；σ0veg是植被的后向散射贡献，σ0soil是土壤的散射贡献；L2是双向植被衰减；θi雷达波束的入射角,A和B是植被参数；V1和V2是植被描述，V1适用叶面指数LAI(m2m-2)，V2适用土壤容积含水量VWC(kgm-2)；S24、取：LAI＝0.091exp(3.7579GDVI2)(R2＝0.932)VWC＝192.64NDVI5-417.46NDVI4+347.96NDVI3–138.93NDVI2+30.699NDVI–2.822(kg·m-2)(R2＝0.990)；A＝0.0045，B＝0.4179；以34.3°为θii雷达波束的入射角，计算植被移除后向散射系数(σ0HH、σ0Hv)。进一步地，所述步骤S3具体包括如下步骤：使用ArcGIS中的点到栅格转换工具，将平均野外测量图转换为光栅单元为30、60和90米大小，然后重放到30m像素，组成数据集。进一步地，所述步骤S4具体包括如下步骤：S41、采用随机森林回归算法RFR建模，将所有12个波段作为输入变量，利用观测值(训练集(TS))在EnMap-Box中的预设值100；以每个节点随机选择的特征或变量数的数目选取所有特征的平方根；结束标准(用于节点分裂)，其中预设值为节点中最小样本数，1，根据基尼指数计算的最小杂质，0；将光栅化EMV或EMH(EM38垂直读数或水平读数)为TS(训练集)参数化后，将生成的RFR模型应用于组合数据集，来预测土壤表观盐度(Eca)。进一步地，还包括将土壤表观盐度Eca转换为土壤电导率ece的步骤，其中，ECe-EM38读数与(Eca)之间的关系表示如下：ECe(dSm-1)＝0.0005EMV2-0.0779EMV+12.655(R2＝0.850)；ECe(dSm-1)＝0.0002EMH2+0.0956EMH+0.0688(R2＝0.791)。本专利技术具有以下有益效果：实现了土壤盐度的预测分析，对农业管理和发展有重要意义，随机森林回归算法RFR用于土壤盐度制图，有更高的精度和更低的归一化均方根误差。具体实施方式下面结合具体实施例对本专利技术进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本专利技术，但不以任何形式限制本专利技术。应当指出的是，对本领域的普通技术人员来说，在不脱离本专利技术构思的前提下，还可以做出若干变形和改进。这些都属于本专利技术的保护范围。本专利技术实施例的一种基于机器学习回归算法的土壤盐度预测方法，包括如下步骤：S1、以等边三角形为测点形状，采用EM38大地电导仪采集三角形的三个角点的实测盐度，所述等边三角形的边长15-20米，确保三角形可以近似表示一个TM(专题成像仪)像素；测量区块的规格1m×1m；测量时，读取垂直读数EMV和水平读数EMH，取平均值作为三角形区域的盐度，测量若干个三角形区域的EMH和EMV为基础真理训练集(TS)；S2、利用L波段雷达后向散射系数去除植被贡献来获得土壤组分；S21、用无线电测量校准从欧洲航天局LANDSAT5(陆地卫星5)TM图像，并采用散射法模型消除附加的大气影响，将所产生的反射率每一波段重新设定为0-1；S22、对从欧洲航天局1.5级雷达结果进行几何校正，并将像素重新调整到12.5米的大小，对水平发射水平接收HH和水平发射垂直接收HV的数字量化值(DN)分别进行校正，按照公式转换为后向散射系数(σ0HH、σ0Hv)；然后应用增强Lee滤波器(EnhancedLeefilter3×3大小，Lee1980)去除斑点或噪声，导出σ0HH和σ0Hv，并重放大到30m像素用来匹配TM数据；S23、构建植被含水量(VWC)对雷达后向散射系数影响的水云模型：L2＝exp(-2BV2sec(θi))其中，σ0是来自植被冠层和土壤的总后向散射系数(σ0HH、σ0Hv)；σ0veg是植被的后向散射贡献，σ0soil是土壤的散射贡献；L2是双向植被衰减；θi雷达波束的入射角,A和B是植被参数；V1和V2是植被描述，V1适用叶面指数LAI(m2m-2)，V2适用土壤容积含水量VWC(kgm-2)；S24、取：LAI＝0.091exp(3.7579GDVI2)(R2＝0.932)VWC＝192.64NDVI5-417.46NDVI4+347.96NDVI3–138.93NDVI2+30.699NDVI–2.822(kg·m-2)(R2＝0.990)；A＝0.0045，B＝0.4179；以34.3°为θi雷达波束的入射角，计算植被移除后向散射系数(σ0HH、σ0Hv)。S3、通过直接栅格化将实测地块转换为栅格来创建训练集(TS)；使用ArcGIS中的点到栅格转换工具，将平均野外测量图转换为光栅单元为30、60和90米大小，然后重放到30m像素，组成数据集；S4、利用上述训练集(TS)，采用随机森林回归(RFR)算法建模应用于联合数据集进行盐度预测。S41、采用随机森林回归算法RFR建模，将所有12个波段作为输入变量，利用24个观测值(训练集(TS))在EnMap-Box中的预设值100；以每个节点随机本文档来自技高网...

【技术保护点】
1.一种基于机器学习回归算法的土壤盐度预测方法，其特征在于：包括如下步骤：S1、以等边三角形为测点形状，采用EM38大地电导仪采集三角形的三个角点的实测盐度，测量时，读取垂直读数EMV和水平读数EMH，取平均值作为三角形区域的盐度，测量若干个三角形区域的EMH和EMV为基础真理训练集(TS)；S2、利用L波段雷达后向散射系数去除植被贡献来获得土壤组分；S3、通过直接栅格化将实测地块转换为栅格来创建训练集(TS)；S4、利用上述训练集(TS)，采用随机森林回归(RFR)算法建模应用于联合数据集进行盐度预测。

【技术特征摘要】
1.一种基于机器学习回归算法的土壤盐度预测方法，其特征在于：包括如下步骤：S1、以等边三角形为测点形状，采用EM38大地电导仪采集三角形的三个角点的实测盐度，测量时，读取垂直读数EMV和水平读数EMH，取平均值作为三角形区域的盐度，测量若干个三角形区域的EMH和EMV为基础真理训练集(TS)；S2、利用L波段雷达后向散射系数去除植被贡献来获得土壤组分；S3、通过直接栅格化将实测地块转换为栅格来创建训练集(TS)；S4、利用上述训练集(TS)，采用随机森林回归(RFR)算法建模应用于联合数据集进行盐度预测。2.如权利要求1所述的一种基于机器学习回归算法的土壤盐度预测方法，其特征在于：所述步骤S1中，所述等边三角形的边长15-20米，确保三角形可以近似表示一个TM像素；测量区块的规格1m×1m。3.如权利要求1所述的一种基于机器学习回归算法的土壤盐度预测方法，其特征在于：所述步骤S2具体包括如下步骤：S21、用无线电测量校准从欧洲航天局获取的LANDSAT5TM图像，并采用散射法模型消除附加的大气影响，将所产生的反射率每一波段重新设定为0-1；S22、对从欧洲航天局获取的1.5级雷达结果进行几何校正，并将像素重新调整到12.5米的大小，对水平发射水平接收HH和水平发射垂直接收HV的数字量化值分别进行校正，按照公式转换为后向散射系数(σ0HH、σ0Hv)；然后应用增强Lee滤波器(EnhancedLeefilter3×3大小，Lee1980)去除斑点或噪声，导出σ0HH和σ0Hv，并重放大到30m像素用来匹配TM数据；S23、构建植被含水量(VWC)对雷达后向散射系数影响的水云模型：L2＝exp(-2BV2sec(θi))其中，σ0是来自植被冠层和土壤的总后向散射系数(σ0HH、σ0Hv)；σ0veg是植被的后向散射贡献，σ0soil是土壤的散射贡献；L2是双向植被衰减；θi雷达波束的入射角...

【专利技术属性】
技术研发人员：吴伟成，周晓亭，祝民强，刘光萍，
申请(专利权)人：东华理工大学，
类型：发明
国别省市：江西,36