本发明涉及一种反演方法，尤其涉及一种利用布谷鸟搜索算法优化bp神经网络微波遥感反演土壤水分的方法，属于遥感反演技术领域。
背景技术：
土壤水分在多种水文模型、气候模型和生态模型中都是及其重要的输入参数。在干旱、洪涝、泥石流滑坡等自然灾害中，土壤水分的分布和时空变化可以建立相应模型，提供更好的帮助和支持。微波遥感能够全天候、全天时进行对地观测，并且对植被有一定的穿透能力，弥补了传统土壤水分监测方法的不足，在干旱区、植被覆盖区土壤水分监测、地表参数反演等研究中具有极大的优势。自1970年以来，雷达土壤水分反演算法不断发展，国内外学者相继提出了多种算法进行土壤水分反演。反演算法物理过程清晰、明确但很复杂，还需付出大量迭代运算代价或输入大量的辅助数据，并且大部分未有效利用主动微波空间分辨率高的特性。尤其对于地面站野外实测数据的需求，在大区域应用时很难满足，土壤参数不易获取，影响了现有模型的实现以及土壤水分反演方法的适用性。
技术实现要素：
为了解决上述现有技术中的不足之处，本发明提供了一种基于优化的bp神经网络微波遥感反演土壤水分的方法。
为了解决上述技术问题，本发明采用的技术方案是：
一种基于优化的bp神经网络微波遥感反演土壤水分的方法，包括以下步骤：
步骤1：获取与雷达数据相同或相近时间内的gf-1光学数据，同时对gf-1数据进行辐射定标，大气校正等预处理，后计算归一化植被指数；
步骤2：对alos-2双极化雷达影像进行滤波，多视，地理编码等预处理，通过预处理分别得到hh和hv两个极化的总后向散射系数；
步骤3：利用水云模型，去除植被和土壤共同产生的后向散射对土壤后向散射系数的影响，分别得到hh和hv两个极化的土壤后向散射系数；
步骤4：在alos-2卫星过境前后，进行实地土壤样本采集，样本采集点均匀分布并且间隔大于3m，土壤样本带回实验室进行处理，采用实验室烘干法获得实测土壤质量含水量和体积含水量；
步骤5：根据bp神经网络建立三层拓扑结构，同时建立对应数据集，包括训练集和验证集，利用布谷鸟搜索算法优化bp神经网络，从而反演得到土壤水分，并与实测土壤水分进行精度验证；
所述的土壤水分反演方法，所述步骤5中，包括以下步骤：
步骤501、确定bp神经网络拓扑结构，根据输入、输出数据的个数，hh和hv作为输入值，对应经纬度实测土壤水分数据作为输出值，建立对应的神经网络模型，确定初始bp神经网络权值、阈值和长度，随后cs对初始值进行编码；
步骤502、根据cs算法，初始化鸟巢数n、维度d、发现概率pa、鸟窝界值以及最大迭代次数maxn，后计算初始化鸟窝适应度值，更新最优鸟窝位置和最优解；
步骤503、不断进行优化，更新适应度较优的鸟窝和最优解，如果满足最大迭代次数，输出全局最优鸟窝和最优值，即反演后的土壤水分值，如果没有满足最大迭代值，即重复步骤502继续优化，直到得到最优解；
步骤504、根据步骤503，获取到最优解后，即反演后的土壤水分值，同实测的土壤水分计算rmse、mae以及mape，进行精度评价，rmse、mae和mape值越小，说明结果越好，得到反演后的土壤水分越准确。
所述的土壤水分反演方法，所述步骤1中，包括以下步骤：
步骤101、从陆地观测卫星数据服务中心订购gf-1影像，gf-1数据时间与alos-2数据相同或相近，并且筛选出云量小于20％的gf-1数据；
步骤102、对gf-1数据的进行预处理，gf-1wfv数据，空间分辨率16m,在使用之前需要进行辐射校正，几何校正，利用envi5.3进行操作；
步骤103、利用gf-1数运算，计算归一化植被指数；
所述的土壤水分反演方法，所述步骤2中，包括以下步骤：
步骤201、根据雷达数据处理软件envi5.3.1 sarscape5.2.1，读取alos-2原始数据，首先设置系统参数，导入数据，生成的数据包括slc数据、地理范围矢量数据、以及相应的kml文件；
步骤202、根据步骤201生成的slc数据，进行多视处理，通过多视处理的slc数据，空间分辨率降低，提升了数据的辐射分辨率；
步骤203、根据步骤202多视处理得到的sar数据继续进行滤波处理，选择frost滤波方法，抑制斑点噪声；
步骤204、根据步骤203滤波后得到的sar数据继续进行地理编码和辐射定标处理，下载对应区域的dem文件，定标结果是无量纲的(单位是linear)，由于后向散射的量级都比较小，实际工作中便于分析往往会转成db的单位，如果要输出db为单位的定标结果(10*log10ofthelinearvalue)，可以选择相应单位的定标结果输出。
所述的土壤水分反演方法，所述步骤3中，包括以下步骤：
步骤301、在水云模型中，由于植被体散射贡献和下垫面散射贡献共同组成了冠层的总后向散射系数，且植被层在一定程度上削减了下垫面散射的贡献，如果在不考虑雷达阴影与地形起伏的影响下，水云模型表达为如下公式：
其中，σo为植被覆盖地表雷达获得的总体后向散射系数(m2/m2)；
τ2＝exp(-2bmv/cosθ)
其中，θ为入射角；mv为冠层中的水分含量(kg/m3)；a、b为经验常数；
步骤302、在alos-2说明文件中，找到水云模型所需要的参数，中心入射角；
步骤303、根据步骤1中gf-1数据计算得到的归一化植被指数作为植被含水量输入；
步骤304、通过最小二乘法计算得到水云模型中的经验系数a、b；
步骤305、确定所有模型需要的参数，编程计算得到去除植被影响后的土壤后向散射系数。
所述的土壤水分反演方法，所述步骤4中，包括以下步骤：
步骤401、根据地形地况在卫星过境前后进行土壤样本采集，采样点的分布考虑植被覆盖、地形、交通通达情况等因素，尽量在研究区域的alos-2整景影像上均布布设；
步骤402、利用环刀法采取土壤样本，取土深度为0-10cm，记录采样点的经纬度、植被覆盖、植被类型等相关信息；
步骤403、将土壤样本带回实验室，采用实验室烘干法获得土壤质量含水量以及体积含水量。
本发明具有以下有益效果：提供了一种利用布谷鸟搜索算法优化bp神经网络反演微波遥感土壤水分反演方法。该方法基于传统的水云模型，去除植被后向散射影响，但由于植被覆盖区土壤粗糙度等因素的限制，加入了布谷鸟搜索算法优化bp神经网络模型，无需考虑土壤粗糙度等因素影响的土壤水分反演方法。该方法放宽了现有方法需要的土壤粗糙度等因素影响，不在受土壤因素限制。为在缺少土壤粗糙度、极化方式少等情况下，依靠布谷鸟搜索算法优化bp神经网络，根据步骤5，通过bp神经网络自身训练及学习规则，以网络误差平方为目标函数、采用梯度下降法来计算目标函数的最小值，计算得到土壤水分，bp神经网络在计算过程中易陷入局部最小值情况，加入布谷鸟算法优化bp神经网络可以增加精确度和预测速度，在一定程度上防止多峰问题陷入局部最优的不足，提高了搜索能力，能获得更优的结果，通过布谷鸟搜索算法优化神经网络可以有效的进行土壤水分反演，为微波遥感土壤水分反演提供了新的解决思路。
附图说明
图1为神经网络结构；
图2为bp训练结果误差；
图3为神经网络训练状况；
图4为神经网络预测误差曲线；
图5为cs网络预测误差；
图6为实测土壤水分和反演土壤水分精度验证曲线；
具体实施方式
以下结合具体实施例，对本发明进行详细说明。
步骤1：获取2020年9月28日与雷达数据相同或相近时间内的gf-1光学数据，利用软件envi5.3，对gf-1数据进行预处理，包括辐射定标，大气校正，后计算归一化植被指数；
步骤101、从陆地观测卫星数据服务中心订购gf-1影像，gf-1数据时间与alos-2数据相同或相近，alos-2过境时间为2020年9月28日，同时筛选出云量小于20％的gf-1数据；
步骤102、对gf-1数据的进行预处理，gf-1wfv数据，空间分辨率16m,在使用之前需要进行辐射校正，几何校正，利用envi5.3进行操作；
步骤103、根据近红外4，红外3波段，计算归一化植被指数ndvi，并裁剪出同alos-2同区域图像；
步骤2：对alos-2双极化雷达影像进行滤波，多视，辐射定标和地理编码等预处理，通过预处理分别得到hh和hv两个极化的总后向散射系数；
步骤201、根据雷达数据处理软件envi5.3.1 sarscape5.2.1，读取alos-2原始数据，首先设置系统参数，导入数据，生成的数据包括slc数据、地理范围矢量数据、以及相应的kml文件；
步骤202、根据步骤201生成的slc数据，进行多视处理，通过多视处理的slc数据，空间分辨率降低，提升了数据的辐射分辨率；
步骤203、根据步骤202多视处理得到的sar数据继续进行滤波处理，选择frost滤波方法，抑制斑点噪声；
步骤204、根据步骤203滤波后得到的sar数据继续进行地理编码和辐射定标处理，下载对应区域的dem文件，定标结果是无量纲的(单位是linear)，由于后向散射的量级都比较小，实际工作中便于分析往往会转成db的单位，为了后续计算，得到以db为单位的后向散射系数。
步骤3：利用水云模型，去除植被、植被和土壤共同产生的后向散射对土壤后向散射系数的影响，分别得到hh和hv两个极化的土壤后向散射系数；
步骤301、在水云模型中，由于植被体散射贡献和下垫面散射贡献共同组成了冠层的总后向散射系数，且植被层在一定程度上削减了下垫面散射的贡献，如果在不考虑雷达阴影与地形起伏的影响下，水云模型表达为如下公式：
其中，σo为植被覆盖地表雷达获得的总体后向散射系数(m2/m2)；
τ2＝exp(-2bmv/cosθ)
其中，θ为入射角；mv为冠层中的水分含量(kg/m3)；a、b为经验常数；
步骤302、在导入sar数据后生成的头文件slc.sml文件中，找到水云模型所需要的参数，中心入射角；
步骤303、根据步骤1中gf-1数据计算得到的归一化植被指数作为植被含水量输入；
步骤304、通过最小二乘法计算得到水云模型中的经验系数a、b；
步骤4：在alos-2卫星过境前后，进行实地土壤样本采集，样本采集点均匀分布并且间隔大于3m，土壤样本带回实验室进行处理，采用实验室烘干法获得实测土壤质量含水量和体积含水量；
步骤401、根据地形地况在卫星过境前后进行土壤样本采集，采样点的分布考虑植被覆盖、地形、交通通达情况等因素，尽量在研究区域的alos-2整景影像上均布布设；
步骤402、利用环刀法采取土壤样本，取土深度为0-10cm，记录采样点的经纬度、植被覆盖、植被类型等相关信息，采样点共120个；
步骤403、将土壤样本带回实验室，采用实验室烘干法获得土壤质量含水量以及体积含水量。
步骤5：根据bp神经网络建立三层拓扑结构，同时建立对应数据集，包括训练集和验证集，利用布谷鸟搜索算法优化bp神经网络，从而反演得到土壤水分，并与实测土壤水分进行精度验证。
步骤501、确定bp神经网络拓扑结构，hh和hv作为输入值，对应经纬度实测土壤水分数据作为输出值，故输入层为2，输出层为1，隐含层为5，其中设置输入值的三分之二为训练集，三分之一为验证集，建立对应的神经网络模型，确定初始bp神经网络权值、阈值和长度，随后cs对初始值进行编码；
步骤502、根据cs算法，初始化鸟巢数n、维度d、发现概率pa、鸟窝界值以及最大迭代次数maxn，后计算初始化鸟窝适应度值，更新最优鸟窝位置和最优解；
步骤503、不断进行优化，更新适应度较优的鸟窝和最优解，如果满足最大迭代次数，输出全局最优鸟窝和最优值，即反演后的土壤水分值，如果没有满足最大迭代值，即重复步骤502继续优化，直到得到最优解；
步骤504、根据步骤503，获取到最优解后，即反演后的土壤水分值，同实测的土壤水分计算rmse、mae以及mape，进行精度评价，rmse、mae和mape值越小，说明结果越好，得到反演后的土壤水分越准确。
根据2020年9月28日alos-2双极化l波段雷达影像，通过预处理提取影像总后向散射，结合gf-1数据，通过水云模型，得到土壤后向散射系数，hh和hv两极化数据。实地采样数据点120个，经过实验室烘干法处理得到土壤质量含水量和体积含水量，通过查阅相关资料，确定当地土壤土质多为黏土、壤土，容重范围多为1.0～1.2g/cm3，故设本实验数据土壤容重为1.0g/cm3。根据实测采样经纬度提取出相同经纬度的hh和hv两极化的土壤后向散射系数，经过筛选，剔除极大、极小的土壤水分数据以及不在范围内的后向散射数据，确定最终样本数量为93。后建立算法输入数据集，训练集个数为62个，验证集个数为31个，神经网络参数设置：2个输入1个输出，隐层节点数为5，最大训练次数100，训练精度0.00001，学习率0.1；鸟巢个数100，发现概率pa为0.7。后进行训练，训练次数100 ，最终确定最优结果。将归一化的数据进行算法训练，得出神经网络结构如图1，图2为bp训练结果误差，图3为神经网络训练状况，图4为神经网络预测误差曲线，图5为cs网络预测误差，图6为实测土壤水分和反演土壤水分精度验证曲线。
实验结果表明，通过光学数据和雷达数据结合，根据cs算法优化神经网络，反演得到土壤水分，均方根误差为2.881，平均绝对误差为2.487，平均绝对百分比误差为0.062，实测土壤水分和反演土壤水分验证精度都较高，说明此方法能有效的提高土壤水分的反演精度。传统方法的建立都是基于土壤粗糙度、相关长度等实测土壤相关参数完整的情况下进行计算，但由于实际情况的影响，在实测采样不能有效的进行情况下，缺失实测土壤参数，不能用传统方法进行土壤水分反演，在仅有的相关数据上，提出的利用cs算法优化神经网络来进行主动微波土壤水分反演，减少了由于数据不足影响反演方法的进行。本发明利用的cs算法优化神经网络来进行主动微波土壤水分反演的方法。本发明放宽了现有方法的前提假设条件，可以弥补数据缺失的不足，更准确的进行主动微波土壤水分反演。
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
