1.本发明涉及一种基顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，属于微波遥感测量技术领域，主要是应用在多基线、高精度的相位信息解缠及数字高程模型重建领域。
背景技术：2.为了获取难以人工量测的陡峭、起伏地段的高程信息，合成孔径雷达干涉测量技术(insar)以信号处理及影像处理结合的形式实现了无需现场量测，就可获取地形测绘的数字高程模型(dem)，这是一项发展迅速的雷达微波遥感测量的现代新型测量技术之一。同时凭借其成像全天候，全天时，穿透力强的工作特点，为数字影像的处理提供了技术支持和技术手段。目前该测量手段为山区测绘、环境监测、地质勘探和气候分析等提供了数据支持及流程化管理，为其进一步硏究有着重要价值。
3.相位解缠技术一直是合成孔径雷达干涉测量技术中的关键。针对反演高程信息中的相位解缠技术，其解缠算法的性能好坏直接关系到高程信息的精确度和准确度。相位解缠在近几十年的研究进程中，已经发展为一门独特且多学科交融的技术，尤其随着深度学习算法的发展，对传统算法的改进和新算法的提出必然成为相位解缠技术新型趋势。
4.目前相位解缠的方法以基线的条数划分主要有单基线相位解缠算法和位移多基线相位解缠算法。单基线相位解缠算法由于受限于相位连续性的假设，例如枝切法，虽然设置了枝切线，有效的避开了含残差点的区域，但对相位欠采样区域的解缠效果不佳，多基线相位解缠算法通过引入多张相位干涉图求解缠绕相位微分的模糊数，并广泛应用于相位不连续的复杂地形，例如中国余数定理，借助求解同余方程组的形式精确同时求解多幅干涉图对应单一点位的模糊数，有较高的精确度，但其受噪声干扰程度大，其噪声鲁棒性差的特点使其不能在实测数据中大规模使用。
5.基于这两种方法的特点，可以借助枝切法设置有效的途径，将路径其扩展中国余数定理相位解缠中，解缠大范围区域的受噪声程度小的相位模糊数，再采用自适应矩形框的方式提取受噪声信息较大的点位，对其采用滤波的形式，再次利用中国余数定理解缠算法处理该区域的相位模糊数，最后拼接大范围区域及各个矩形框区域的解缠相位。
技术实现要素：6.本发明是为了克服现有技术中的相位解缠的方法，单基线枝切法算法对相位突变解缠存在孤岛现象以及多基线中国余数定理算法抗噪性弱的缺点，提供一种顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，采用单基线相位解缠算法结合多基线相位解缠方法的技术，既解决了孤岛区域难以解缠的难题，又增强了算法的噪声鲁棒性。
7.为实现上述目的，本发明采用以下技术方案：
8.本发明的一种顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，所述方法包括以下步骤：
9.针对同一感兴趣区域生成多基线insar相位干涉图；
10.采用枝切法对所述干涉图设置枝切线，并采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位；
11.针对干涉图中的枝切线分布区域，采用枝切法生成局部范围内的枝切线，并再次采用中国余数定理进行解缠，获得局部区域的解缠信息；
12.将大范围区域的相位解缠信息及局部区域的解缠信息结合，获得干涉图全部点位的相位解缠信息。
13.作为优选，所述针对同一感兴趣区域生成多基线insar相位干涉图，进一步包括：
14.步骤1，获取针对同一区域的多幅多基线sar影像，选择一幅主影像及多幅辅影像；
15.步骤2，所述主影像与辅影像逐一进行差分干涉处理，并按基线的长度排列，生成多幅待解缠的多基线insar相位干涉图。
16.作为优选，所述的采用枝切法对所述干涉图设置枝切线，并采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位，进一步包括：
17.步骤3，采用枝切法对所述的干涉图预处理，生成枝切线；
18.步骤4，根据所述的多幅干涉图建立多基线insar几何模型，根据多基线insar几何模型搭建中国余数定理数学模型；
19.步骤5，根据预处理的枝切线分布，避开枝切线采用合理的途径，使用中国余数定理对干涉图解缠，获得大范围区域的相位解缠信息。
20.作为优选，所述的针对干涉图中的枝切线分布区域，采用枝切法生成局部范围内的枝切线，并再次采用中国余数定理进行解缠，获得局部区域的解缠信息，进一步包括：
21.步骤6：针对干涉图中的枝切线分布区域采用自适应矩形框逐一提取；
22.步骤7：对框内的干涉相位采用枝切法重新选取展开点，生成局部范围内的枝切线；
23.步骤8：对局部范围内的枝切线上点的相位信息进行滤波处理；
24.步骤9：对矩形框内滤波处理后的干涉相位再次使用中国余数定理进行解缠，获得局部区域的解缠信息。
25.作为优选，所述步骤4，进一步包括：
26.对于多基线insar系统几何模型的参数表达式如下：
[0027][0028]
式中，h表示高差，λ表示波长，θ表示侧视角，α表示基线水平角，bi表示不同的基线长度，r表示相位中心相对于目标点的斜距，ki表示不同点位干涉相位的模糊数。
[0029]
若设置b0＝[b1,b2,
…
,bn]，其中，[]代表求解最小公倍数，同时设定模长mi＝b0/bi，干涉相位可以用高差表示：
[0030][0031]
若设置则不同基线下干涉图中对应于各个点位干涉相位的所要求解的模糊数满足以下等式：
[0032]
t＝wi+ki·
mi[0033]
对应于n组基线所形成的n幅干涉图，每个相对应的点位由n组等式组成，且每组等式中，mi是确定的，wi不同，t相同，根据数学理论，即构建了同余方程组，最终所要获取的值模糊数ki是求解中国余数定理数学模型的过程。
[0034]
作为优选，所述步骤5，进一步包括：根据预处理的枝切线分布，对于枝切线衍生较多的地形突变区域，采用中国余数定理对干涉图解缠，获得大范围区域的相位解缠信息。
[0035]
作为优选，所述步骤6中，干涉图的尺寸为m
×
n，矩形框的尺寸为u
×
v，所述自适应矩形框的尺寸满足公式：
[0036]
u＝floor[0.1
×
m+ω]
[0037]
v＝floor[0.1
×
n+ω]
[0038]
其中，ω为干涉图的信噪比。
[0039]
作为优选，所述步骤8中的滤波处理采用中值滤波处理。
[0040]
本发明提供的一种单基线相位解缠算法结合多基线相位解缠方法，考虑了枝切法孤岛现象的受限性及中国余数定理算法的抗噪能力弱，利用枝切线及自适应矩形框提取噪声强的区域，利用相关的滤波算法减弱噪声对干涉相位信息的干扰，利用两次中国余数定理算法分别解缠噪声信息弱的大面积区域及噪声信息强的局部区域，拼接其解缠的相位信息实现抗噪的多基线相位解缠算法。
附图说明
[0041]
图1是本技术实施例三待解缠的干涉图差分处理。
[0042]
图2是本技术实施例三多基线insar几何模型。
[0043]
图3是本技术实施例三中对应美国isolation peak国家公园真实dem数据三维及二维参考图。
[0044]
图4是本技术实施例三三条基线所对应的差分干涉图，其中按基线的长度排序。
[0045]
图5是本技术实施例三的流程图。
[0046]
图6是本技术实施例三对三幅干涉图采用枝切法处理全局范围所搜索的残差点。
[0047]
图7是本技术实施例三对三幅干涉图采用枝切法处理全局范围所形成的枝切线。
[0048]
图8是本技术实施例三采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位。
[0049]
图9是本技术实施例三采用枝切法对处理任意一处欠采样区域在矩形框内形成的局部枝切线。
[0050]
图10是本技术实施例三将所有矩形框内相位信息中值滤波之后再次采用中国余数定理所整合的解缠相位。
[0051]
图11是本技术实施例三将大范围解缠相位及局部区域所拼接的整幅影像的解缠相位。
[0052]
图12是本技术实施例三获取的解缠相位同参考相位作差的误差分布图。
[0053]
图13是本技术实施例三获取的解缠相位同参考相位作差的误差直方图。
具体实施方式
[0054]
下面结合附图和具体实施方式对本发明做进一步描述。
[0055]
实施例一
[0056]
如图1所示，本发明的一种顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，所述方法包括以下步骤：
[0057]
步骤s1，针对同一感兴趣区域生成多基线insar相位干涉图。
[0058]
步骤s2，采用枝切法对所述干涉图设置枝切线，并采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位。
[0059]
步骤s3，针对干涉图中的枝切线分布区域，采用枝切法生成局部范围内的枝切线，并再次采用中国余数定理进行解缠，获得局部区域的解缠信息。
[0060]
步骤s4，将大范围区域的相位解缠信息及局部区域的解缠信息结合，获得干涉图全部点位的相位解缠信息。
[0061]
本方案考虑了枝切法孤岛现象的受限性及中国余数定理算法的抗噪能力弱，利用枝切线及自适应矩形框提取噪声强的区域，利用相关的滤波算法减弱噪声对干涉相位信息的干扰，利用两次中国余数定理算法分别解缠噪声信息弱的大面积区域及噪声信息强的局部区域，拼接其解缠的相位信息实现抗噪的多基线相位解缠算法。
[0062]
实施例二
[0063]
本实施例是在实施例一基础上的优化方案，本实施例中分别对实施例一的步骤s1、s2、s3、s4进行了进一步扩展
[0064]
步骤s1具体进一步包括以下步骤：
[0065]
步骤1，获取针对同一区域的多幅多基线sar影像，选择一幅主影像及多幅辅影像。
[0066]
步骤2，所述主影像与辅影像逐一进行差分干涉处理，并按基线的长度排列，生成多幅待解缠的多基线insar相位干涉图。
[0067]
步骤s2具体进一步包括以下步骤：
[0068]
步骤3，采用枝切法对所述的干涉图预处理，生成枝切线。
[0069]
步骤4，根据所述的多幅干涉图建立多基线insar几何模型，根据多基线insar几何模型搭建中国余数定理数学模型；
[0070]
步骤5，根据预处理的枝切线分布，避开枝切线采用合理的途径，使用中国余数定理对干涉图解缠，获得大范围区域的相位解缠信息。
[0071]
具体的，所述步骤3采用枝切法生成枝切线涉及到二维相位解缠理论及残差点。
[0072]
首先，二维相位解缠是针对m
×
n大小的矩形空间而言的，其中，选取一个空间内的点作为起点向不同的方向展开相位，则二维的相位解缠路径分为水平向和垂直向，其中，设定起点位置的干涉相位φ与解缠相位存在以下关系：
[0073][0074]
其中，定义w[φ(i,j)]为相位缠绕函数，而上述提及的路径中，对于水平向及垂直向的相位梯度，可以定义为：
[0075]
△iφ(i,j)＝φ(i,j)-φ(i-1,j)
ꢀꢀ
(2)
[0076]
△jφ(i,j)＝φ(i,j)-φ(i,j-1)
ꢀꢀ
(3)
[0077]
根据路径的不同，利用相位缠绕函数，对相邻像元垂直向和水平向做缠绕相位梯
度：
[0078]
△iw[φ(i,j)]＝
△iφ(i,j)+2π
△ik(i,j)
ꢀꢀ
(5)
[0079]
△jw[φ(i,j)]＝
△jφ(i,j)+2π
△jk(i,j)
ꢀꢀ
(6)
[0080]
对上述不同方位的缠绕相位梯度再次利用相位缠绕函数，可以得到相邻像元垂直向和水平向做缠绕相位：
[0081]
w2[
△iw1[φ(i,j)]＝
△iφ(i,j)+2π[
△ik1(i,j)+
△ik2(i,j)]
ꢀꢀ
(7)
[0082]
w3[
△jw1[φ(i,j)]＝
△jφ(i,j)+2π[
△jk3(i,j)+
△jk4(i,j)]
ꢀꢀ
(8)
[0083]
要保障缠绕相位在-π到π之间，式中对应的各个参数要满足以下关系式：
[0084]
π≤
△iφ(i,j)《π
ꢀꢀ
(9)
[0085]
π≤
△jφ(i,j)《π
ꢀꢀ
(10)
[0086]
△ik1(i,j)+
△ik2(i,j)＝0
ꢀꢀ
(11)
[0087]
△jk3(i,j)+
△jk4(i,j)＝0
ꢀꢀ
(12)
[0088]
接下来，将会采取可靠的路径积分来实现解缠，先垂直向积分再水平向积分，先水平向积分再垂直向积分。
[0089]
先垂直向积分再水平向积分：
[0090]
首先进行垂直方向的首列进行像元解缠处理：
[0091][0092]
将首列的像元解缠之后，作为新的起始列，此时每一列依次对前一列进行水平向梯度积分处理，j列处理完毕则完成解缠过程
[0093][0094]
先水平向积分再垂直向积分：
[0095]
首先进行水平方向的首列进行像元解缠处理：
[0096][0097]
将首行的像元解缠之后，作为新的起始行，此时每一行依次对前一行进行水平向梯度积分处理，i行处理完毕则完成解缠过程
[0098][0099]
理想情况下，规定为相位连续的情形，根据上述两种路径积分的任意一种途径，就可从一点缠绕相位扩展到整幅影像的真实相位，但实际情况中，有噪声的存在，则路径的选取就会产生偏差，甚至会将点位误差积累，扩展为更大区间的误差。
[0100]
基于上述情况，需要引入残差点的概念，即通过求解相邻像元各个方位梯度的方式判定相位连续性：
[0101]
△1＝φ(i+1,j)-φ(i,j)
ꢀꢀ
(17)
[0102]
△2＝φ(i+1,j+1)-φ(i+1,j)
ꢀꢀ
(18)
[0103]
△3＝φ(i,j+1)-φ(i+1,j+1)
ꢀꢀ
(19)
[0104]
△4＝φ(i,j)-φ(i,j+1)
ꢀꢀ
(20)
[0105]
将求解的各个方位梯度求和即为残差：
[0106][0107]
利用求解的残差判定相位信息的连续性，即：t＝0，该点判定为相位连续点，t≠0，该点判定为残差点，通常定义为：t》0，该点判定为正残差点，t《0，该点判定为负残差点。
[0108]
残差点的引入对路径跟踪算法有着引导作用，尤其对于单基线解缠算法保证了大范围区域内相位连续性的假设。
[0109]
具体的，所述步骤4，进一步包括：
[0110]
对于多基线insar系统几何模型的参数表达式如下：
[0111][0112]
式中，h表示高差，λ表示波长，θ表示侧视角，α表示基线水平角，bi表示不同的基线长度，r表示相位中心相对于目标点的斜距，ki表示不同点位干涉相位的模糊数。
[0113]
若设置b0＝[b1,b2,
…
,bn]，其中，[]代表求解最小公倍数，同时设定模长mi＝b0/bi，干涉相位可以用高差表示：
[0114][0115]
若设置则不同基线下干涉图中对应于各个点位干涉相位的所要求解的模糊数满足以下等式：
[0116]
t＝wi+ki·
miꢀꢀ
(24)
[0117]
对应于n组基线所形成的n幅干涉图，每个相对应的点位由n组等式组成，且每组等式中，mi是确定的，wi不同，t相同，根据数学理论，即构建了同余方程组，最终所要获取的值模糊数ki是求解中国余数定理数学模型的过程。
[0118]
具体的，所述步骤5进一步包括，根据预处理的枝切线分布，对于枝切线衍生较多的地形突变区域，采用中国余数定理对干涉图解缠，获得大范围区域的相位解缠信息。
[0119]
步骤s3具体进一步包括以下步骤：
[0120]
步骤6：针对干涉图中的枝切线分布区域采用自适应矩形框逐一提取；
[0121]
步骤7：对框内的干涉相位采用枝切法重新选取展开点，生成局部范围内的枝切线；
[0122]
步骤8：对局部范围内的枝切线上点的相位信息进行滤波处理，所述步骤8中的滤波处理采用中值滤波处理；
[0123]
步骤9：对矩形框内滤波处理后的干涉相位再次使用中国余数定理进行解缠，获得局部区域的解缠信息。
[0124]
具体的，所述步骤6中，干涉图的尺寸为m
×
n，矩形框的尺寸为u
×
v，所述自适应矩形框的尺寸满足公式：
[0125]
u＝floor[0.1
×
m+ω]
ꢀꢀ
(25)
[0126]
v＝floor[0.1
×
n+ω]
ꢀꢀ
(26)
[0127]
其中，ω为干涉图的信噪比。
[0128]
实施例三
[0129]
本实施例是应用实施例二的方法，利用美国isolationpeak国家公园数据作为一
个具体操作实例，对本技术方案的详细说明。
[0130]
如图1、图2、图3、图4所示，本实施例采用的数据为一景主影像及三景辅影像，根据多基线几何模型及差分原理生成3组不同基线下所对应的干涉图作为待解缠的数据。如图5所示，以下本实施例的方法具体实施流程的进一步说明。
[0131]
步骤一：将四景影像区分出主影像与辅影像，进行差分干涉处理，并按基线的长度排列，形成多基线insar相位干涉图。
[0132]
步骤二：为了数据的准确性及实验对比分析，图3给出了美国isolation peak国家公园真实dem数据三维及二维参考图。
[0133]
步骤三：对三幅干涉图分别采用枝切法处理全局范围所搜索的正负残差点，采用公式计算，正残差点用1表示，负残差点用-1表示。
[0134]
步骤四：利用步骤三提取的残差点对全局相位设置枝切线，作路径跟踪指导作用。
[0135]
步骤五：基于干涉图中已经生成的枝切线，采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位。
[0136]
图6是对三幅干涉图采用枝切法处理全局范围所搜索的残差点，图7是对三幅干涉图采用枝切法处理全局范围所形成的枝切线。
[0137]
图8是采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位。
[0138]
步骤六：全局范围内的解缠相位生成，则需要考虑局部欠采样及相位突变处的区域，采用枝切法逐一对欠采样及相位突变处的区域处理，以形成矩形框内的局部枝切线，并标记该位置，其中，实验过程中有一处欠采样区域相位突变点较多，为了确保60％稳定的相位信息，自适应矩形框的尺寸扩展为98*52，展示为图9所示。
[0139]
步骤七：对枝切线上点的相位信息进行中值滤波处理，若矩形框内仅有1-2个相位突变点，则默认该矩形框内解缠效果良好，无需滤波处理再次解缠。
[0140]
步骤八：对矩形框内滤波处理后的干涉相位再次使用中国余数定理进行解缠，结果如图10所示。
[0141]
步骤九：将整体的相位解缠信息及小区域解缠信息结合，则得到了干涉图全部点位的相位解缠信息；如图11所示是大范围解缠相位及局部区域所拼接的整幅影像的解缠相位。
[0142]
步骤十：为验证该算法的抗噪性及相位突变处的解缠性能，作解缠相位同参考相位作差的误差分布图及误差直方图。
[0143]
由图11可见，该算法整合的多基线相位解缠图同原始的dem整体吻合程度较高，并没有出现孤岛现象。
[0144]
由图12可见，该算法作解缠相位同参考相位作差的误差分布图中，整体相位误差值较低，仅有几处亮点为忽略的矩形框内仅有1-2个没有经过滤波处理的相位突变点，其误差值较大。
[0145]
由图13可见，误差直方图分布也均在0rad左右，突变的点位占的频率很低，则该算法的鲁棒性有所提升，抗噪能力得到了验证。
技术特征：1.一种顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述方法包括以下步骤：针对同一感兴趣区域生成多基线insar相位干涉图；采用枝切法对所述干涉图设置枝切线，并采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位；针对干涉图中的枝切线分布区域，采用枝切法生成局部范围内的枝切线，并再次采用中国余数定理进行解缠，获得局部区域的解缠信息；将大范围区域的相位解缠信息及局部区域的解缠信息结合，获得干涉图全部点位的相位解缠信息。2.根据权利要求1所述的顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述针对同一感兴趣区域生成多基线insar相位干涉图，进一步包括：步骤1，获取针对同一区域的多幅多基线sar影像，选择一幅主影像及多幅辅影像；步骤2，所述主影像与辅影像逐一进行差分干涉处理，并按基线的长度排列，生成多幅待解缠的多基线insar相位干涉图。3.根据权利要求1所述的顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述的采用枝切法对所述干涉图设置枝切线，并采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位，进一步包括：步骤3，采用枝切法对所述的干涉图预处理，生成枝切线；步骤4，根据所述的多幅干涉图建立多基线insar几何模型，根据多基线insar几何模型搭建中国余数定理数学模型；步骤5，根据预处理的枝切线分布，避开枝切线采用合理的途径，使用中国余数定理对干涉图解缠，获得大范围区域的相位解缠信息。4.根据权利要求1所述的顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述的针对干涉图中的枝切线分布区域，采用枝切法生成局部范围内的枝切线，并再次采用中国余数定理进行解缠，获得局部区域的解缠信息，进一步包括：步骤6：针对干涉图中的枝切线分布区域采用自适应矩形框逐一提取；步骤7：对框内的干涉相位采用枝切法重新选取展开点，生成局部范围内的枝切线；步骤8：对局部范围内的枝切线上点的相位信息进行滤波处理；步骤9：对矩形框内滤波处理后的干涉相位再次使用中国余数定理进行解缠，获得局部区域的解缠信息。5.根据权利要求3所述的顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述步骤4，进一步包括：对于多基线insar系统几何模型的参数表达式如下：式中，h表示高差，λ表示波长，θ表示侧视角，α表示基线水平角，b
i
表示不同的基线长度，r表示相位中心相对于目标点的斜距，k
i
表示不同点位干涉相位的模糊数。若设置b0＝[b1,b2,
…
,b
n
]，其中，[]代表求解最小公倍数，同时设定模长m
i
＝b0/b
i
，干
涉相位可以用高差表示：若设置则不同基线下干涉图中对应于各个点位干涉相位的所要求解的模糊数满足以下等式：t＝w
i
+k
i
·
m
i
对应于n组基线所形成的n幅干涉图，每个相对应的点位由n组等式组成，且每组等式中，m
i
是确定的，w
i
不同，t相同，根据数学理论，即构建了同余方程组，最终所要获取的值模糊数k
i
是求解中国余数定理数学模型的过程。6.根据权利要求3所述的顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述步骤5，进一步包括：根据预处理的枝切线分布，对于枝切线衍生较多的地形突变区域，采用中国余数定理对干涉图解缠，获得大范围区域的相位解缠信息。7.根据权利要求4所述的顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述步骤6中，干涉图的尺寸为m
×
n，矩形框的尺寸为u
×
v，所述自适应矩形框的尺寸满足公式：u＝floor[0.1
×
m+ω]v＝floor[0.1
×
n+ω]其中，ω为干涉图的信噪比。8.根据权利要求4所述的顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，其特征是，所述步骤8中的滤波处理采用中值滤波处理。
技术总结本发明公开了一种顾及枝切线的鲁棒中国余数定理的多基线相位解缠方法，所述方法包括以下步骤：针对同一感兴趣区域生成多基线InSAR相位干涉图；采用枝切法对所述干涉图设置枝切线，并采用中国余数定理避开枝切线采取合适的途径获取的大范围解缠的相位；针对枝切线分布区域生成局部范围内的枝切线，并再次采用中国余数定理进行解缠；将大范围区域的相位解缠信息及局部区域结合获得最终相位解缠信息。本方案利用枝切线及自适应矩形框提取噪声强的区域，利用相关的滤波算法减弱噪声对干涉相位信息的干扰，利用两次中国余数定理算法分别解缠噪声信息弱的大面积区域及噪声信息强的局部区域，拼接其解缠的相位信息实现抗噪的多基线相位解缠算法。多基线相位解缠算法。多基线相位解缠算法。
技术研发人员：王启贵 周明 毛永生 王毅 陈兴红
受保护的技术使用者：浙江华东测绘与工程安全技术有限公司
技术研发日：2022.07.23
技术公布日：2022/11/1
