基于雷达后向散射实测数据对单层地表介电参数与粗糙度参数快速反演的联合优化算法
【专利摘要】本发明公开了一种遗传与多输出支持向量机的联合优化算法，结合雷达后向散射实测数据，快速反演单层地表(均方根高度kσ＜1.5，均方根斜率s＜0.3)介电参数与粗糙度参数，其步骤包括：利用单层地表HH和VV极化雷达后向散射系数实测数据，获得同极化比值；根据粗糙地表电磁散射的小斜率近似方法计算后向同极化比；采用遗传算法结合小斜率近似同极化比以及实测数据反演地表介电常数；将其反演结果带入电磁散射积分方程模型，生成后向散射系数随粗糙度变化的数据文件；结合两种极化的雷达后向散射实测数据，构成目标函数，利用M-SVR优化算法反演地表粗糙度参数，并评估反演误差和效率。本发明在保证反演精度的同时，可实现对地表参数的实时预测。
【专利说明】基于雷达后向散射实测数据对单层地表介电参数与粗糙度参数快速反演的联合优化算法
【技术领域】
[0001]本发明涉及微波遥感【技术领域】，尤其是涉及一种基于雷达实测数据对单层地表介电参数和粗糙度参数快速反演的联合优化算法。
【背景技术】
[0002]电磁散射和电磁逆散射研究在微波遥感中起着极为重要的作用，对地微波遥感的雷达技术近些年来取得了长足的进步，其能够全天候、全天时观测，具有广泛的应用前景。在陆地和海洋等典型地表的微波遥感中，其后向散射回波与地表的物理参数以及土壤水分等密切相关，其包含了众多的地表信息，例如利用地表微波遥感技术可以获得土壤的含水量、农作物成熟情况以及海水的海态、风速、介电常数等信息，因此逆散射相关领域的科学研究将对国民技术的发展产生重大影响，对地表参数的反演研究工作势在必行。
[0003]反演地表的介电参数和粗糙度参数是微波遥感的重要应用之一，近二十多年来，地表参数的反演已发展为多波段、多极化、多角度地表数据的反演，地表介电常数和粗糙度的反演取得了显著的发展。目前，反演方面的研究主要有四条路径:一是经验公式法，如Oh模型，Dubois模型和Shi模型等，但其各自适用范围有限。二是优化算法结合地表电磁散射模型和雷达实测数据，反演粗糙面的介电参数和粗糙度参数，当电磁散射模型较为复杂时，计算速度相当缓慢。一般常采用遗传算法或粒子群算法分别结合积分方程(IEM)模型、微扰法模型以及双尺度模型联合反演电磁参数。王悦泉和金亚秋等基于裸土后向散射系数实测数据结合遗传算法和双尺度模型反演了地表粗糙度和湿度参数，反演结果和实测值吻合良好，但该反演只局限于满足两尺度模型的地表。电子科技大学贾明全利用L / S / C /X波段裸土单频全极化和双频同极化后向散射实测数据，结合神经网络方法和AIEM法反演了裸土的介电常数和粗糙度参数，其反演精度较高，但耗时较长，无法实现地表参数的实时预测。三是优化算法结合地表电磁散射同极化比值和雷达实测数据，反演介电参数。不同地表电磁模型的同极化比适用范围不同。该反演方法在实验数据较为理想的情况下，运算速度快，反演精度高。Ceraldi采用地表电磁散射模型的微扰法、基尔霍夫近似法的镜像同极化比消去粗糙度的影响，结合遗传算法反演了粗糙面的介电常数，但其只适用于微起伏、曲率不大的粗糙面。四是基于机器学习算法的地表参数反演。传统的支持向量机(SVM)只能进行地表单参数反演，因此无法保证多参数的总反演精度，对于地表多参数同时未知的情况，SVM变得无能为力。M-SVR相比于SVM具有多参数反演能力，且具有精度高、速度快等优点，根据目前检索国内外资料表明，M-SVR应用于地表电磁参数反演还未见报道。
【发明内容】
[0004]本发明所要解决的技术问题是:对于雷达波照射下的地表回波电磁散射模型，如常规的微扰法、基尔霍夫近似法、双尺度法、小斜率近似以及积分方程法等，其地表适用范围不同。已有的地表后向散射模型同极化比仅针对微扰法和基尔霍夫近似法，但这两者的适用范围较窄，使得在未知地表粗糙度的情况下采用同极化比方法反演该粗糙面介电常数，在模型选用上具有一定的盲目性；而传统的IEM地表散射模型在较宽的粗糙度适用范围内和实测数据吻合良好，与GA方法联合反演电磁参数的精度较高，但其公式繁琐，反演时间漫长，难以实现地表参数的实时预测。
[0005]为解决上述技术问题，本发明的技术方案是:在地表粗糙度满足k O <2，均方根斜率s < 0.3的情况下，推导SSA近似同极化比公式，拓宽所反演地表粗糙面的适用范围，在其均方根斜率较小的情况下均具有精确解；采用IEM法生成训练数据，由M-SVR反演参数，解决反演过程中的小样本、非线性以及局部极小等问题，在保证反演精度的同时，减少了反演时间。基于雷达后向散射实测数据对地表介电参数与粗糙度参数反演的联合优化算法，包括如下步骤:
[0006]I)根据单层地表雷达后向散射系数实测数据，获得同极化比值；
[0007]2)计算粗糙地表电磁散射的小斜率近似方法(SSA)后向同极化比；
[0008]3)采用GA算法结合SSA近似同极化比与实测数据同极化比反演地表介电常数；
[0009]4)将地表介电常数反演结果带入IEM电磁散射模型，降低未知数维度，生成后向散射系数随均方根高度和相关长度变化的文件(Al)；
[0010]5)将数据文件(Al)作为训练模型，采用M-SVR建立目标函数E 2并优化为模型Model-1，将地表HH和VV极化雷达后向散射系数实测数据作为测试样本，优化反演地表粗糙度信息(均方根高度和相关长度)；
[0011]6)评估反演误差RMS和反演所需时间。
[0012]进一步的，其中，步骤2)所述的SSA近似方法的同极化比公式为:
【权利要求】
1.基于雷达后向散射实测数据对单层地表介电参数与粗糙度参数快速反演的联合优化算法，其特征在于，包括如下步骤: 1)获取地表HH和VV极化雷达后向散射系数实测数据，处理实测数据获得同极化比值； 2)根据小斜率近似方法(SSA)，理论计算粗糙地表雷达后向散射系数的同极化比； 3)根据I)和2)，建立同极化比值的目标函数EI，采用遗传算法优化反演地表介电常数； 4)将反演的地表介电常数，代入粗糙面电磁散射积分方程(IEM)，降低未知数维度，生成HH和VV极化雷达后向散射系数随均方根高度和相关长度变化的数据文件(Al)； 5)将数据文件(Al)作为训练模型，采用M-SVR方法建立目标函数E2，其最优模型为Model-1，将地表HH和VV极化雷达后向散射系数实测数据作为测试样本，利用模型Model-1优化反演地表粗糙度信息(均方根高度和相关长度)； 6)评估地表介电参数与粗糙度散射的反演误差和反演所需时间。
2.根据权利要求1所述的基于雷达后向散射实测数据对单层地表介电参数与粗糙度参数快速反演的联合优化算法，其特征在于，其中，步骤2)所述的SSA近似方法的同极化比公式为:
3.根据权利要求1所述的基于雷达后向散射实测数据对单层地表介电参数与粗糙度参数快速反演的联合优化算法，其特征在于，其中，步骤3)所述的采用遗传算法结合SSA近似同极化比理论值和实测数据同极化比值，反演地表介电常数的目标函数E I取为:
4.根据权利要求1所述的基于雷达后向散射实测数据对单层地表介电参数与粗糙度参数快速反演的联合优化算法，其特征在于，其中，步骤5)所述的采用M — SVR方法结合单层地表电磁散射IEM模型反演粗糙度参数(均方根高度和相关长度)，按如下步骤进行: 4a)选择入射角和IEM模型HH与VV极化后向散射系数理论值作为训练输入向量X(mX2 维)； 4b)将训练输入向量对应的均方根高度和相关长度作为训练输出向量Y(mX2维)；4c)采用M-SVR方法结合训练输入向量X和训练输出向量Y，建立目标函数E 2并优化为模型Model-1 ； 对于非线性问题，支持向量机利用非线性映射函数将样本映射到高维线性空间，并建立一个线性模型来估计回归函数 Y=f (X) =W ? (X) +b(3) 其中W=[wi，w2，...，wk]T，b=[bl，b2，...，bk]T为高维特征空间的K维线性回归量，所有的训练数据在精度e下无误差地用线性函数拟合。其中，w是可调的权值向量，b是偏置.采用e-支持向量机，可转化为如下约束优化问题
5.评估其反演误差RMS和反演所需时间。
【文档编号】G06F19/00GK103617344SQ201310544995
【公开日】2014年3月5日 申请日期:2013年11月2日 优先权日:2013年11月2日
【发明者】吴振森, 张元元, 苏翔, 李海英, 令狐龙翔 申请人:西安电子科技大学