CN115856889B 一种误差自动校正的InSAR时序变形监测方法 （北京数慧时空信息技术有限公司）

一种误差自动校正的InSAR时序变形监测方法本发明提供一种误差自动校正的InSAR时序到解缠干涉图集合；S3识别目标区域内的稳定点并用于干涉对轨道误差与大气相位的联合校2S2设置时空基线阈值，基于时空基线阈值构建短基线集的差S3识别目标区域内的稳定点，基于稳定点对解缠S32计算解缠干涉图的相位标准差，并设置第一相位S34对质量筛选后的解缠干涉图集合进行堆叠的时序IS36根据稳定点的位置信息计算轨道误差和分层大气的联合相位，将其作为误差信号S42根据第二相干性阈值和第二相位标准差阈值对校正干涉图集合进行质量筛选，若S43采用三角相位闭环方法对质量筛选后的校正干涉图集合进行筛选，将质量筛选后i通过最小二乘法求解轨道误差与分层大气的联合相位3叠干涉图个数阈值Num_ifgthre，根据间歇性方式对质量筛选后的解缠干涉图集合进行堆S53根据滤波法对每个时刻的累积变形进行湍流大气校正和噪声去除处理，得到目标0MX1＝GMXNVNX1+AMX1Δh1X10(N-2)X1=αR(N-2)XNVNX1+0(N-2)X1Δh1X14对速率项和高程残差的进行求解，并对每个S531根据空间低滤波对每个时刻的累积变形进行处理，以去除累积变形中的噪声相S22将SAR影像的DEM数据转换为地形相位值，将干涉图集合和DEM生成地形相位做差5[0001]本发明涉及地形监测技术领域，尤其涉及一种误差自动校正的InSAR时序变形监[0002]干涉雷达指采用干涉测量技术的合成孔径雷达（InSAR是新近发展起来的空间[0008]S3识别目标区域内的稳定点，基于稳定点对解缠干涉图集合进行轨道误差和大[0010]S5基于形变算法对优质干涉图集合进行解算和处理，得到目标区域的地面形变[0014]S33根据第一相干性阈值和第一相位标准差阈值对解缠干涉图集合进行质量筛6[0015]S34对质量筛选后的解缠干涉图集合进行堆叠的时序InSAR处理，根据堆叠时序[0016]S35设定稳定点阈值区间，根据相位变化速率得到目标区域内[0020]p(x,y)=a,+azx+azy+a3xy+ax2+asy2+agh(x,y)+e(x,y)数阈值，根据间歇性方式对质量筛选后的解缠干涉图集合进行堆叠时序形[0029]S42根据第二相干性阈值和第二相位标准差阈值对校正干涉图集合进行质量筛7[0043]引入虚拟观测v,=0为起始数据，并附加加速度相等的平滑约束条件进行求解，V2-V1=V-VZ(K=1,2,3,..,N-2)[0048]其中a为平滑约束系数，GXW为速率项系数，AX1为高程残差项系数，[0051]S531根据空间低滤波对每个时刻的累积变形进行处理，以去除累积变形中的噪8[0056]S22将SAR影像的DEM数据转换为地形相位值，将干涉图集合和DEM生成地形相位[0075]S3识别目标区域内的稳定点，基于稳定点对解缠干涉图集合进行轨道误差和大9[0077]S5基于形变算法对优质干涉图集合进行解算和处理，得到目标区域的地面形变[0082]S22将SAR影像的DEM数据转换为地形相位值，将干涉图集合和DEM生成地形相位[0092]第三步：根据pm(l,p)和ps(l,p)来计算主影像和辅影像之间的偏移量口按照整像元偏移量在不同的行列之中来计算两[0105]H(u,v)=sfIF(u,v)le·F(u,v)(2)H(u,v)=s{IF(u,v)ly-F·F(u,v)(3)[0112]Y(i,i)=w(p(i,i))=p(i,)+2nn(i,i)(4)Y(i,i)为缠绕相位，p(i,)表示解缠相位，残差定义为：[0118]x=max(0,k,),x5=min(0,ka)(7)[0125]S33根据第一相干性阈值和第一相位标准差阈值对解缠干涉图集合进行质量筛[0127]S34对质量筛选后的解缠干涉图集合进行堆叠的时序InSAR处理，根据堆叠时序[0128]S35设定稳定点阈值区间，根据相位变化速率得到目标区域内用间歇性相干点进行堆叠时序形变解算。当堆叠干涉图个数Num_ifg>Num_ifg-。形变量的[50%-Ratioj-,50%+Ratiojp]区间时，该点被认为是稳定点。[0141]S42根据第二相干性阈值和第二相位标准差阈值对校正干涉图集合进行质量筛涉图的相位标准差Std2，并设置第二相位标准差阈值，若是或合条件的校正干涉图放入劣质校正干涉图集合2中进行筛除，符合条件的校正干涉图具有环，生成三个解缠干涉差分相位为unw12、unw23、unw13：设定闭合环相位阈值unw12、unw23、umnw13、unwy4、unw34，当解缠参考基准一致时，理论上存在于干涉图滤波及解缠误差的存在，等式右侧往往不为0。为了选取到解缠误差严重的干涉断标准。当满足RMS(unw1z+unwz有解缠相位；当时，则剔除闭合环相关的[0158]其中，为相干系数且ye[0,1l，M和S分别表示构成干涉对的两个SAR影像的局[0159]计算时间区间上每个像元的相干系数。假设有N+1幅SAR影像，可以形成N幅干涉图，每个像元就形成一个相干系数序列yz(k=1,2,..,N)。[0169]引入虚拟观测v,=0为起始数据，并附加加速度相等的平滑约束条件进行求解，[0170]VX2-V1=V1-V(K=1,2,3,...,N-2)(15)[0179]S531根据空间低滤波对每个时刻的累积变形进行处理，以去除累积变形中的噪共计23景哨兵SAR数据。卫星姿态为升轨方向，卫星航向角为一13.47°,雷达入射角为[0185]采用堆叠时序InSAR方法识别目标区域内的稳定点，并基于稳定点信息对所有解[0187]对优质的解缠