【技术实现步骤摘要】
多轨道InSAR提取高精度三维形变方法
[0001]本专利技术涉及微波遥感测量领域，具体地涉及多轨道InSAR提取高精度三维形变方法。
技术介绍
[0002]干涉合成孔径雷达(InSAR)技术在过去二十年中一直是监测山体滑坡，地震，火山爆发，采矿沉降等地质灾害的有效方法。然而，InSAR技术仅提供雷达视线向(Line of sight，简称LOS)的一维变形。在实际应用中，多数情况下，LOS向变形不适用于评估灾害范围。
[0003]近年来，随着SAR卫星数量的增多，获取不同航向和入射角(升轨和降轨数据，或左视和右视数据)的多组LOS向形变数据成为可能。利用3组LOS向形变数据建立观测方程组，能够解得垂直、东西、南北向形变数据。但受SAR卫星近南北飞行轨道的限制，南北向解算结果精度远低于垂直、东西向。
技术实现思路
[0004]为解决现有技术中的技术问题，本专利技术提供多轨道InSAR提取高精度三维形变方法。包括以下内容:
[0005]本专利技术的第一方面，提供一种多轨道InSAR提取高精度三维
【技术保护点】
【技术特征摘要】
1.一种多轨道InSAR提取高精度三维形变方法，其特征在于，包括：从多轨道的SAR数据中获取每个轨道的LOS向形变；根据SAR传感器侧视成像几何关系，建立LOS向形变与三维方向形变关系方程组，建立LOS向形变与忽略南向北形变的二维方向形变关系方程组；根据LOS向形变与三维方向形变关系方程组和LOS向形变与二维方向形变关系方程组，利用最小二乘法，得到三维方向的形变结果以及垂直向和东西向形变；根据LOS向形变误差和三维方向的形变结果误差关系，建立三维形变误差传递模型，根据LOS向形变误差和垂直向和东西向形变误差关系，建立LOS向形变误差和二维形变误差传递模型；根据二维形变误差传递模型的垂直向及东西向系统误差和南北向形变分量的关系，建立额外系统误差传递模型；基于额外系统误差传递模型、三维解算南北向形变和忽略南北向解算的垂直向及东西向形变与LOS向形变的关系建立目标函数，并对目标函数添加正则化约束项；基于带有约束项的目标函数结合设定的误差传递模型，得到目标函数矩阵；基于目标函数矩阵得到满足目标函数极小值的LOS向形变误差，将LOS向形变误差代入三维形变误差传递模型，得到三维形变误差；根据三维形变误差得到三维形变结果。2.根据权利要求1所述的多轨道InSAR提取高精度三维形变方法，其特征在于，所述建立LOS向形变与三维方向形变关系方程组中，方程组为：式中，表示i轨道LOS向形变位移，i＝1,2,3；θ
i
，α
i
表示雷达卫星入射角和北向与雷达卫星飞行方向的顺时针方向夹角，i＝1,2,3，u
′
，e
′
，n
′
分别为垂直向、东西向和南北向三维形变。3.根据权利要求2所述的多轨道InSAR提取高精度三维形变方法，其特征在于，建立LOS向形变与忽略南向北形变的二维方向形变关系方程组中，方程组为：式中，u
″
，e
″
为忽略南北向时垂直向和东西向形变。4.根据权利要求3所述的多轨道InSAR提取高精度三维形变方法，其特征在于，所述根据LOS向形变与三维方向形变关系方程组和LOS向形变与二维方向形变关系方程组，利用最小二乘法，得到三维方向的形变结果以及垂直向和东西向形变中，公式为：D＝(A
T
WA)
‑1A
T
WL，式中A为：
L为雷达视线向形变，分别为：W为LOS向形变权重矩阵，采用等价权重矩阵，D为解算所得三维形变或忽略南北向得到的二维形变为：5.根据权利要求4所述的多轨道InSAR提取高精度三维形变方法，其特征在于，所述根据LOS向形变误差和三维方向的形变结果误差关系，建立三维形变误差传递模型为：式中，式中σ
li
为i轨道LOS向形变误差，a
i
，b
i
，c
i
分别为i轨道垂直向、东西向、南北向形变系数，σ
u
′
，σ
e
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。