本发明涉及雷达遥感以及图像处理技术领域，尤其涉及一种联合ps目标与ds目标的双极化相位优化地表形变监测方法。
背景技术：
合成孔径雷达干涉测量(insar)技术是20世纪60年代产生的一项空间对地观测技术，它能够不受光照和天气条件的限制实现全天候、全天时对地观测，还可以穿透地表和植被获取地表下信息。1998年，产生了极化合成孔径雷达干涉测量技术，该技术研究同一目标的不同极化方式，通过统计分析同一地区多幅sar影像时间序列上幅度与相位信息稳定性，探测保持高相干性的点为ps点。通过这些稳定ps点目标，可以估计并消除大气效应相位的贡献，获得更加精确的三维信息和形变量。
然而，传统的ps点选取方法往往因为自身的局限性，导致选点数量与选点精度的不足。如何解决上述方法的不足，是本发明主要解决的问题。
技术实现要素：
为了提高sentinel-1a数据中相干点的数量和优化相干点目标(ps与ds目标)相位，本发明提出了一种极化干涉雷达领域中的ds目标识别以及极化相位优化方法，该方法利用vv、vh极化幅度信息选取ps与ds目标，并基于ps目标与ds目标的时空间相干矩阵，利用特征值分解极化优化技术，抑制相位中的噪声分量，提高相位质量。
bgsm算法是双极化中ps点目标与ds目标优化算法中的方法之一。通过在投影矢量双变量范围内以一定的步长取值，寻找最优变量优化目标函数。
为了实现上述目的，本发明提供了一种联合ps目标与ds目标的双极化相位优化地表形变监测方法，包括如下步骤：
s1、获取n景原始的vv极化、vh极化slc影像，并对所述slc影像进行裁剪、配准、多视预处理，得到多视后的slc影像；
s2、基于所述多视后的slc影像分别计算vv极化数据和vh极化数据的平均幅度值，并设置阈值选取ps目标，将vv极化数据和vh极化数据选取的ps目标的总和作为最终的ps目标；采用anderson-darling检验和平均幅度阈值法识别得到ds目标；
s3、利用投影矢量和pauli矢量将vv极化、vh极化的复数像元组合成新的slc影像，并对所述新的slc影像中的ds目标进行自适应均值滤波处理，得到滤波后的ds目标；
s4、分别计算得到所述新的slc影像中的ps目标集与ds目标集的时空间相干矩阵；
s5、通过特征值分解方法分解所述时空间相干矩阵，得到ps目标集的特征值和ds目标集的特征值；
s6、通过bgsm算法搜索最优投影矢量使得最大特征值占所有特征值总和的比重最大，得到优化后的干涉相位；
s7、根据所述优化后的干涉相位联合ps目标与ds目标反演地表形变，得到地表形变速率图。
本发明提供的技术方案带来的有益效果是：本发明利用平均幅度阈值法分别选取vv、vh双极化数据的ps目标，将两种极化数据的ps目标合并，利用anderson-darling检验和平均幅度阈值法识别ds目标，并对ds目标进行自适应均值滤波处理；分别针对ps目标与ds目标生成时空间相干矩阵，采用相干矩阵特征值分解方法分离出最大特征值，并采用bgsm算法寻找最优投影矢量使得第一主成分贡献率达到最大；最终相较于单极化数据能有效增加相干目标点数量，抑制干涉相位中的噪声成分，提高相干点目标相位质量。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明实施例的可见光遥感影像；
图2是本发明实施例联合ps目标与ds目标的双极化相位优化地表形变监测方法的流程图；
图3是本发明实施例ps目标与ds目标相位优化与其他方法的结果对比；
其中，图3(a)为ps目标与ds目标时空间相干矩阵特征值分解极化相位优化结果；
图3(b)为esm极化干涉相位；
图3(c)为vv极化干涉相位。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
参考图1，图1是本发明实施例的可见光遥感影像，本发明实施例采用覆盖青藏铁路部分路段的sentinel-1avv、vh双极化sar数据来完成提出的相干点目标识别与相位优化算法，该数据集包含17景sar影像，时间范围为2017年03月21日至2018年04月21日，重访周期为12天，距离向与方位向视数比为4：1，研究区大小为10001000像元。
参考图2，图2是本发明实施例联合ps目标与ds目标的双极化相位优化地表形变监测方法的流程图；
本发明实施例提供了一种联合ps目标与ds目标的双极化相位优化地表形变监测方法，包括如下步骤：
s1、获取n景原始的vv极化、vh极化slc影像，并对所述slc影像进行裁剪、配准、多视预处理，得到多视后的slc影像。
本实施中采用的slc数据为原始burst条带，但是本发明不限于此，预处理的方法如下：
s11：在sar影像范围内建立感兴趣区，将所有slc原始影像裁剪至与感兴趣区同一范围，此时burst条带的数据分辨率没有损失；
s12：选取vv、vh极化数据中时间较为居中的一幅影像为主影像，将其余从影像配准到主影像空间，配准精度为亚像元级，保证不同影像上同一地物对应的像元一一对应；
s13：根据配准结果进行8：1多视处理，使之大小与实际地物形状一致。
s2、基于所述多视后的slc影像分别计算vv极化数据和vh极化数据的平均幅度值，并设置阈值选取ps目标，将vv极化数据和vh极化数据选取的ps目标的总和作为最终的ps目标；采用anderson-darling检验和平均幅度阈值法识别ds目标。
本实施例中，s2的具体步骤如下：
s21、对vv极化数据、vh极化数据分别计算平均幅度值，并根据阈值选取ps目标，最终的ps目标为两种极化数据得到ps目标的总和，计算公式如下：
其中，n为影像数量；
s22、采用anderson-darling检验和平均幅度阈值法识别ds目标。
s3、利用投影矢量和pauli矢量将vv极化、vh极化的复数像元组合成新的slc影像，并对所述新的slc影像中的ds目标进行自适应均值滤波处理，得到滤波后的ds目标。
本实施例中，s3的具体步骤如下：
s31、利用投影矢量和pauli矢量将vv、vh两种的复数像元组合成新的slc影像，组合方式如下：
k＝[svv2svh]t
s＝ω *k
其中，svv、shv分别表示vv、vh极化对应的复数像元，k为pauli基矢量化结果，ω为由α、φ两个待优化变量表达的极化散射机制，α与目标的散射类型有关，φ与目标的散射方向有关，s为两种极化状态组合而成的新的slc影像的像元值， 表示共轭转置；
s32、对新的slc影像中的ds目标进行自适应均值滤波处理，计算公式如下：
其中，ψn为中心像元相位；nn为估计窗口内像元数总和；a(i)为第i个像元的振幅；ψ(i)为第i个像元的相位。
s4、分别计算所述新的slc影像中的ps目标集与ds目标集的时空间相干矩阵。
本实施例中，s4的具体步骤如下：
s41：对新的slc影像中的ps目标集和ds目标集计算时空间协方差矩阵，计算公式如下：
其中，x＝[x1,x2…,xn]t表示分布式目标的同质点在n景slc影像上的复观测量；nshps为相干目标点总数量；ω表示相干目标点集合；h表示复共轭转置；
s42：将时空间协方差矩阵进行归一化处理，得到时空间相干矩阵，计算公式如下：
其中，
s5、通过特征值分解方法分解ps目标集和ds目标集的时空间相干矩阵，得到ps目标集的特征值和ds目标集的特征值。
根据特征值分解方法，时空间相干矩阵分解为：
其中，
s6、通过bgsm算法搜索最优投影矢量使得最大特征值(第一主成分)占所有特征值总和的比重最大，进行干涉相位优化。
本实施例中，s6的具体步骤如下：
s61、取待优化变量α、φ初始值分别为0、-π；
s62、通过bgsm算法搜索全局最优投影矢量组(α,φ)使得最大特征值贡献率最大，计算公式如下：
其中，
s7、根据优化后的干涉相位反演形变，得到地表形变速率图。
本实施例中，s7的具体步骤如下：
s71、建立稀疏格网相位解缠，利用delaunay不规则三角网建立ps候选点之间的连接关系，然后基于加权的最小二乘法进行相位解缠；
s72、aps(大气相位)的估计和去除，通过时间维和空间维的滤波分离非线性形变和大气相位；
s73、相干点目标(ps目标与ds目标)形变估计，在去除大气相位后，重新相位解缠，得到更加精确的形变量估计。
请参考图3，图3是本发明实施例ps目标与ds目标相位优化与其他方法的结果对比，图3(a)为ps目标与ds目标时空间相干矩阵特征值分解极化相位优化结果，图3(b)为esm极化干涉相位，图3(c)为vv极化干涉相位；图3(a)中相干目标数量为276362，图3(b)中ps目标数量为168164，图3(c)中ps目标数量为30502，可见本发明相较于单极化数据选取的相干点目标数量明显增加。此外，就图3中的相位质量而言，本发明相较于单极化方法的干涉相位更加紧密和连续，有效抑制了噪声信号。
本发明提供了一种联合ps目标与ds目标的双极化相位优化地表形变监测方法，利用平均幅度阈值法分别选取vv、vh双极化数据的ps目标，将两种极化数据的ps目标合并，利用anderson-darling检验和平均幅度阈值法识别ds目标，并对ds目标进行自适应均值滤波处理；分别针对ps目标与ds目标生成时空间相干矩阵，采用相干矩阵特征值分解方法分离出最大特征值，并采用bgsm算法寻找最优投影矢量使得第一主成分贡献率达到最大。最终相较于单极化数据能有效增加相干点目标(ps目标与ds目标)数量，还能提高主导信号比重，从而抑制噪声成分，提高相干点目标相位质量。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。
1.一种联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，包括如下步骤：
s1、获取n景原始的vv极化、vh极化slc影像，并对所述slc影像进行裁剪、配准、多视预处理，得到多视后的slc影像；
s2、基于所述多视后的slc影像分别计算vv极化数据和vh极化数据的平均幅度值，并设置阈值选取ps目标，将vv极化数据和vh极化数据选取的ps目标的总和作为最终的ps目标；采用anderson-darling检验和平均幅度阈值法识别得到ds目标；
s3、利用投影矢量和pauli矢量将vv极化、vh极化的复数像元组合成新的slc影像，并对所述新的slc影像中的ds目标进行自适应均值滤波处理，得到滤波后的ds目标；
s4、分别计算得到所述新的slc影像中的ps目标集与ds目标集的时空间相干矩阵；
s5、通过特征值分解方法分解所述时空间相干矩阵，得到ps目标集的特征值和ds目标集的特征值；
s6、通过bgsm算法搜索最优投影矢量使得最大特征值占所有特征值总和的比重最大，得到优化后的干涉相位；
s7、根据所述优化后的干涉相位联合ps目标与ds目标反演地表形变，得到地表形变速率图。
2.根据权利要求1所述的联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，所述s1具体包括以下步骤：
s11、获取n景原始的vv极化、vh极化slc影像；
s12、在所述slc影像范围内建立感兴趣区，将所有slc原始影像裁剪至与感兴趣区同一范围；
s13、选取vv极化、vh极化的某一景slc影像为主影像，将其余从影像配准到主影像空间；
s14、根据配准结果进行距离向和方位向视数比为8：1的多视处理，得到多视后的slc影像。
3.根据权利要求1所述的联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，所述s2具体包括以下步骤：
s21、对vv极化数据、vh极化数据分别计算平均幅度值，并根据阈值选取ps目标，最终的ps目标为两种极化数据得到ps目标的总和，计算公式如下：
其中，n为影像数量；
s22、采用anderson-darling检验和平均幅度阈值法识别ds目标。
4.根据权利要求1所述的一种联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，所述s3具体包括以下步骤：
s31、利用投影矢量和pauli矢量将vv、vh两种的复数像元组合成新的slc影像，组合方式如下：
k＝[svv2svh]t
s＝ω *k
其中，svv、shv分别表示vv、vh极化对应的复数像元，k为pauli基矢量化结果，ω为由α、φ两个待优化变量表达的极化散射机制，α与目标的散射类型有关，φ与目标的散射方向有关，s为两种极化状态组合而成的新的slc影像的像元值， 表示共轭转置；
s32、对新的slc影像中的ds目标进行自适应均值滤波处理，计算公式如下：
其中，ψc为中心像元相位；nn为估计窗口内像元数总和；a(i)为第i个像元的振幅；ψ(i)为第i个像元的相位。
5.根据权利要求1所述的联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，所述s4具体包括如下步骤：
s41：对新的slc影像中的ps目标集和ds目标集计算时空间协方差矩阵，计算公式如下：
其中，x＝[x1，x2...，xn]t表示分布式目标的同质点在n景slc影像上的复观测量；nshps为相干目标点总数量；ω表示相干目标点集合；h表示复共轭转置；
s42：将时空间协方差矩阵进行归一化处理，得到时空间相干矩阵，计算公式如下：
其中，
6.根据权利要求1所述的联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，所述s5中，根据特征值分解方法，时空间相干矩阵分解为：
其中，
7.根据权利要求1所述的联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，所述s6具体步骤如下：
s61、取待优化变量α、φ初始值分别为0、-π；
s62、通过bgsm算法搜索全局最优投影矢量组(α，φ)使得最大特征值贡献率最大，以得到优化后的相位，计算公式如下：
其中，
8.根据权利要求1所述的联合ps目标与ds目标的双极化相位优化地表形变监测方法，其特征在于，所述s7具体步骤如下：
s71、基于ps目标与ds目标优化后的干涉相位建立稀疏格网相位解缠，得到地表；
s72、在相位中估计和去除aps；
s72、根据去除aps后的相位对ps目标和ds目标形变估计。
技术总结