SAR图像的干涉相位 matlab_好文推荐：复杂环境下地基SAR粗差探测及应用

地基合成孔径雷达干涉技术(GB-InSAR)具有大范围、非接触、高精度与高时空分辨率、不受云雾天气影响等优点，广泛应用于矿区、滑坡、冰川和大坝等变形监测中。文献[1]将意大利Citrin山谷滑坡不同时期的地基SAR数据与同期GPS数据进行对比，得到了一致结果；文献[2]利用地基SAR系统成功对Alpine冰川位移进行实时监测，证明了该系统可用于冰川位移监测，且文献[3]利用IBIS-L系统对意大利Cancano大坝进行了位移测量，与正(倒)垂线方法测量结果的偏差均方值在0.2 mm内；文献[4]利用IBIS-L系统分别对紫坪铺大坝和滑坡体进行监测试验，获得了优于毫米级的监测结果；文献[5]利用IBIS-S系统对长江大桥进行桥梁静载试验，证明了该方法在实际桥梁结构安全监测中的可行性，且文献[6]提出了永久散射体网法，将复杂环境下的气象影响控制在亚毫米级；文献[7]采用IBIS-L系统对地面沉降进行了监测，表明该技术可用于毫米级精度的沉降监测；文献[8]利用地基SAR结合常规工程测量、岩土体探测技术及数值模拟技术，并联合监测巴西拉萨格勒拉火车站施工期间动态变化，得到了机械测量数据和数值力学模型变形解算结果一致的效果；文献[9]利用地基SAR对意大利阿尔卑斯山进行了不连续监测，证明地基SAR可用于平均每年几厘米位移的滑坡监测；文献[10]针对PS点选择的局限性，提出了多重阈值优化的PS点选择方法，其选取结果的稳定性和质量都优于单一模型。

以上地基SAR技术的应用研究均是在较为理想的状态下进行的，针对复杂环境下施工场地的地基SAR形变监测，雷达视线常被施工机械作业等间断性遮挡，导致部分影像出现相位奇异值，影响监测结果的精度与可靠性。本文提出改进的基于小波变换的信号奇异性检测，将遮挡影像识别转化为粗差探测问题，以解决外界干扰对监测结果带来的不利影响。

1 时序相位的奇异性及其检测 

1.1 地基SAR形变监测原理

雷达只能记录到目标回波相位中不满整周的相位信息，即所获取的相位是缠绕的，对同一目标不同时刻的两期雷达影像进行差分干涉可得其干涉相位为

(1)

式中，ϕ_ω为干涉相位；φ_disp为与实际位移有关的相位分量；φ_atm为大气相位分量；φ_noise为各种噪声相位分量，可选取高相关性的目标点及滤波处理进行削弱^[11]；k₁、k₂分别为目标点在两个时刻相位的整周模糊度。经消噪与大气校正，目标在雷达视线(LOS)向的真实形变值由解缠后的相位差来求解

(2)

式中，λ为波长。

1.2 时序相位的奇异性检测

雷达视线的遮挡会导致式(1)中得到错误的k_i值，而相邻影像目标点k_i值的大幅变化使其时序相位产生突变，直接由式(2)进行形变解算会得出错误的结果。时序相位在某时刻的突变所引起的信号非连续属于第一类间断点，称时序相位在此处具有奇异性^[12]。因此，为解决受遮挡影像对形变解算结果带来的不利影响，可通过检测PS点时序相位的奇异性，确定相位突变所发生的时空位置，剔除相应的受遮挡影像。基于奇异性检测的地基SAR受遮挡影像识别流程如图 1所示。

图 1 地基SAR遮挡影像识别流程