CN103323844A - 一种多通道干涉合成孔径雷达高程重建方法及装置 - Google Patents

Abstract

本发明公开了一种多通道干涉合成孔径雷达InSAR高程重建方法，该方法包括：根据干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为待测高程建模同余方程组；对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数；根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。本发明同时还公开了一种多通道InSAR高程重建装置。采用本发明的技术方案，能够利用较少的样本对地形进行高程重建。

Description

一种多通道干涉合成孔径雷达高程重建方法及装置

技术领域

本发明涉及多通道干涉合成孔径雷达的高程重建技术，尤其涉及一种多通道干涉合成孔径雷达（InSAR，Interferometric Synthetic Aperture Radar）高程重建方法及装置。

背景技术

InSAR技术是一种高精度对地三维观测技术。在InSAR技术中，获取地形表面的数字高程图（高程重建）是其中最重要的应用方向之一。而多通道InSAR高程重建则特别适用于解决较复杂地区的高程叠掩问题，而且有利于高度变化大于一个高度模糊数的情况下的相位解缠。

目前存在的多通道InSAR高程重建方法主要分为两大类：基于最大似然估计的高程重建和基于最大后验概率估计的高程重建。第一类高程重建算法将各通道干涉相位的概率密度函数的乘积作为似然函数，通过搜索使该似然函数的值最大的高度值而获得高程重建结果；第二类高程重建算法则利用了隐藏在干涉相位图中的相邻像素之间的关系，并以此作为先验信息，结合第一类算法中的似然函数，从而得到地表高度的最大后验概率估计。

以上两类高程重建算法均属于统计方法，需要较多的独立样本，而当独立样本的个数较少时，就不能通过上述方法进行高程重建了。

发明内容

有鉴于此，本发明的主要目的在于提供一种多通道干涉合成孔径雷达高程重建方法及装置，能够利用较少的样本对地形进行高程重建。

为达到上述目的，本发明的技术方案是这样实现的：

本发明提供了一种多通道干涉合成孔径雷达InSAR高程重建方法，通过设置于测量航天器上的至少两个测量雷达对待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系；该方法包括：

根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为所述待测高程建模同余方程组；

对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数；

根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。

上述方案中，所述根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为待测高程建模同余方程组，包括：

确定所述待测高程的缠绕干涉相位对应的斜距差；

根据所述斜距差、以及所述待测高程与所述至少两个测量雷达之间的位置关系，确定所述待测高程与干涉相位之间的关系；

根据设定的余数，对所述待测高程与干涉相位之间的关系进行等价变换，得到所述待测高程的同余方程组。

上述方案中，所述对所述同余方程组进行等价变换包括：

从同余方程组中选取同余方程作为参考方程；

将同余方程组中除参考方程以外的同余方程与所述参考方程进行相减，并将相减后的所有方程等价变化为一组同余方程组。

上述方案中，所述求解变换后的同余方程组包括：根据鲁棒性中国剩余定理，求解变换后的同余方程组。

上述方案中，根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程，包括：

将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中而确定出所述同余方程组的高程。

本发明还提供了一种多通道InSAR高程重建装置，该装置包括：第一确定单元、建模单元、获取单元以及第二确定单元；其中，

所述第一确定单元，用于通过设置于测量航天器上的至少两个测量雷达对待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系；

所述建模单元，用于根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为所述待测高程建模同余方程组；

所述获取单元，用于对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数；

所述第二确定单元，用于根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。

上述方案中，所述建模单元包括：第一确定子单元、第二确定子单元和第一变换子单元；其中，

所述第一确定子单元，用于确定所述待测高程的缠绕干涉相位对应的斜距差；

所述第二确定子单元，用于根据所述斜距差、以及所述待测高程与所述至少两个测量雷达之间的位置关系，确定所述待测高程与干涉相位之间的关系；

所述第一变换子单元，根据设定的余数，对所述待测高程与干涉相位之间的关系进行等价变换，得到所述待测高程的同余方程组。

上述方案中，所述获取单元包括：选取子单元和第二变换子单元；其中，

所述选取子单元，用于从同余方程组中选取同余方程作为参考方程；

所述第二等价变换子单元，用于将同余方程组中除参考方程以外的同余方程与所述参考方程进行相减，并将相减后的所有方程等价变换为一组同余方程组。

上述方案中，所述获取单元，还用于根据鲁棒性中国剩余定理，求解变换后的同余方程组。

上述方案中，所述第二确定单元包括：回调子单元和第三确定子单元；其中，

所述回调子单元，用于将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中；

所述第三确定子单元，用于在将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中之后，确定出所述同余方程组的高程。

本发明提供的多通道干涉合成孔径雷达高程重建方法及装置，通过设置于测量航天器上的至少两个测量雷达对所述待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系；根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为待测高程建模同余方程组；对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数；根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。如此，能够利用较少的样本对地形进行高程重建。

附图说明

图1为本发明多通道InSAR高程重建方法实现流程示意图一；

图2为本发明多通道InSAR高程重建方法实现流程示意图二；

图3为利用仿真数据模拟的高程示意图；

图4为利用干涉相位与高程之间的关系仿真的缠绕干涉相位示意图；

图5为利用本发明提供的方法得到的高程结果示意图；

图6a为本发明多通道InSAR高程重建装置的组成结构示意图；

图6b为本发明多通道InSAR高程重建装置中建模单元的组成结构示意图；

图6c为本发明多通道InSAR高程重建装置中获取单元的组成结构示意图；

图6d为本发明多通道InSAR高程重建装置中第二确定单元的组成结构示意图。

具体实施方式

为了能够更加详尽地了解本发明的特点与技术内容，下面结合附图对本发明的实现进行详细阐述，所附附图仅供参考说明之用，并非用来限定本发明。

本发明提供了一种多通道InSAR高程重建方法，如图1所示，该方法包括以下步骤：

步骤101：通过设置于测量航天器上的至少两个测量雷达对待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及待测高程与所述至少两个测量雷达之间的位置关系。

步骤102：根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为待测高程建模同余方程组。

本步骤具体包括：

确定所述待测高程的缠绕干涉相位对应的斜距差；

根据所述斜距差、以及所述待测高程与所述至少两个测量雷达之间的位置关系，确定所述待测高程与干涉相位之间的关系；

根据设定的余数，对所述待测高程与干涉相位之间的关系进行等价变换，得到所述待测高程的同余方程组。

步骤103：对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数。

这里，对所述同余方程组进行等价变换，具体为：

从同余方程组中选取同余方程作为参考方程；

将同余方程组中除参考方程外的同余方程与所述参考方程进行相减，并将相减后的所有方程等价变换为一组同余方程组。

本方案中，将同余方程组中除参考方程以外的同余方程与所述参考方程进行相减是为了消除未知参数—高程。

这里，所述求解变换后的同余方程组包括：根据鲁棒性中国剩余定理，求解变换后的同余方程组。

步骤104：根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。

本步骤具体为：将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中而确定出所述同余方程组的高程。

下面结合具体实施例对图1所示的方法做进一步详细描述，如图2所示，图2为本发明实施例提供的一种多通道InSAR高程重建实现方法，该方法包括以下步骤：

步骤201：通过设置于测量航天器上的两个测量雷达对待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及待测高程与所述至少两个测量雷达之间的位置关系。

步骤202：确定所述待测高程的缠绕干涉相位对应的斜距差。

这里，所述斜距是指测量雷达到待测高程位置之间的距离。

由于测量雷达接收到的电磁波的干涉相位与雷达到待测高程位置的斜距差有关，因此可以通过测量雷达接收到的干涉相位确定所述干涉相位对应的斜距差。

步骤203：根据所述斜距差、以及高程与测量雷达之间的位置关系，确定所述高程与干涉相位之间的关系。

本实施例中，获取的高程与干涉相位之间的关系是在去除平地效应后所获取的；平地效应是指水平地面上高度相同的两个物体由于卫星的距离不同所产生的相位差异。平地效应使得产生的干涉条纹过于密集，因此获取的高程与干涉相位之间的关系之前需要去除平地效应。

这里，干涉相位与高程之间的关系为：

$\frac{2\mathrm{\π}}{h_k}\×h={\mathrm{\φ}}_k+n_k\×2\mathrm{\π},k=\mathrm{1,2},...,K---\left(1\right)$

公式（1）中，h表示所求的地面上某个像素点的高程，k表示干涉通道的序号，φ_k表示第k组干涉通道上观测到的缠绕干涉相位值，n_k表示φ_k对应的模糊数，K表示干涉通道的总组数，h_k表示第k组干涉通道的高度模糊数。

优选地，高度模糊数h_k可以表示为：

$h_k=\frac{{\mathrm{\λ}}_k\×R_0\×\sin \mathrm{\θ}}{2\×B_{\⊥k}}---\left(2\right)$

公式（2）中，λ_k表示第k组干涉通道的系统波长，B_⊥k表示第k组干涉通道对应的垂直基线的长度，θ为系统的下视角，R₀为主图像的景中心斜距。这里，系统是指多通道InSAR系统。

步骤204：根据设定的余数，对所述高程与干涉相位之间的关系进行等价变换，得到所述高程的同余方程组。

这里，设定的余数可以表示为：

$r_k=\frac{{\mathrm{\φ}}_k}{2\mathrm{\π}}\×h_k,k=\mathrm{1,2},...,K---\left(3\right)$

对所述高程与干涉相位之间的关系进行等价变换，具体为：将公式（3）和公式（2）代入公式（1）中，并对公式（1）进行等价变换，得到下式：

h=n_k×h_k+r_k,k=1,2,...,K    （4）

公式（4）中，第k组干涉通道的高度模糊数h_k、以及余数r_k均为已知参数（可以通过多通道InSAR系统测出）；模糊数n_k未知参数；h为待求的同余方程组的解，为未知参数。这里，高度模糊数h_k和模糊数n_k是两个不同的参量。

这里，公式（4）为第k个同余方程的表达式，通过公式（4）可以得到K个同余方程，这K个同余方程进而组成同余方程组，即：

h=n₁×h₁+r₁

h=n₂×h₂+r₂

h=n₃×h₃+r₃   （5）

......

h=n_K×h_K+r_K

步骤205：从同余方程组中选取同余方程作为参考方程。

具体的，设定一组互质的数，这组互质的数中第k个数可以表示为：

Γ_k=h_k/M,1≤k≤K   （6）

公式（6）中，M是根据系统参数设定好的归一化因子，且M为实数。

根据公式（6），设定参数Γ和参数γ_k表示如下：

$\mathrm{\Γ}\stackrel{\mathrm{\Δ}}{=}{\mathrm{\Γ}}_1\×{\mathrm{\Γ}}_2\×...\×{\mathrm{\Γ}}_K$

（7）

${\mathrm{\γ}}_k\stackrel{\mathrm{\Δ}}{=}{\mathrm{\Γ}}_1\×...\×{\mathrm{\Γ}}_{k-1}\×{\mathrm{\Γ}}_{k+1}\×...\×{\mathrm{\Γ}}_K=\mathrm{\Γ}/{\mathrm{\Γ}}_k$

其中，

为定义符号，这里表示Γ的定义式为Γ₁Γ₂...Γ_K；同理，γ_k的定义式为Γ₁...Γ_k-1Γ_k+1...Γ_K。

通过公式（7）可知γ_k与Γ_k互质，因此γ_k关于Γ_k的乘法逆元

一定存在，也即满足以下同余方程的

一定存在：

$\stackrel{\‾}{{\mathrm{\γ}}_k}{\mathrm{\γ}}_k\≡1\mathrm{mod}{\mathrm{\Γ}}_K---\left(8\right)$

其中，

表示γ_k关于Γ_k的乘法逆元；≡为恒等符号。

将所有带噪声的余数

除以M取余，得到K个不同的新余数：

其中，

表示带噪声的余数；

表示由公式（9）所得的新余数；

表示下取整符号。

公式（9）中，

表示第k个新余数。

对公式（9）中的新余数进行平均，得到：

${\hat{r}}^c=\underset{x\∈R,0\≤x\≤M}{\arg \min }\underset{k=1}{\overset{K}{\mathrm{\Σ}}}{d}^2({\widehat{r_k}}^c,x)---\left(10\right)$

其中，∈表示属于符号；argmin表示取泛函的最小值，这里具体为取

的最小值。

公式（10）中，距离函数d(x,y)设定为：

$d(x,y)\stackrel{\mathrm{\Δ}}{=}x-y-j_{0}M---\left(11\right)$

公式（11）中，x，y为实数，且范围为大于等于0且小于M；j₀满足下式：

$j_0=\underset{j\∈\{-\mathrm{1,0,1}\}}{\arg \min }|x-y-j_0\×M|---\left(12\right)$

其中，||表示取绝对值符号，这里具体为取（x-y-j₀×M）的绝对值。

距离平均值

最近的所在的同余方程即为选取的参考方程。这里，距离平均值

最近的

也即

取最小值时的

步骤206：将同余方程组中的同余方程与所述参考方程进行等价变换。

具体的，将同余方程组中除参考方程以外的同余方程与参考方程相减，并将相减后的所有方程等价变换为一组同余方程组：

${\hat{n}}_1={\hat{n}}_2\×\frac{{\mathrm{\Γ}}_2}{{\mathrm{\Γ}}_1}+\frac{{\hat{q}}_{\mathrm{2,1}}}{{\mathrm{\Γ}}_1}$

${\hat{n}}_1={\hat{n}}_3\×\frac{{\mathrm{\Γ}}_3}{{\mathrm{\Γ}}_1}+\frac{{\hat{q}}_{\mathrm{3,1}}}{{\mathrm{\Γ}}_1}---\left(13\right)$

...

${\hat{n}}_1={\hat{n}}_K\×\frac{{\mathrm{\Γ}}_K}{{\mathrm{\Γ}}_1}+\frac{{\hat{q}}_{K,1}}{{\mathrm{\Γ}}_1}$

公式（13）中，

这里，[]表示四舍五入取整符号。

步骤207：根据鲁棒性中国剩余定理，求解变换后的同余方程组。

具体的，求解公式（13）中的同余方程组，得到的解以及与该同余方程组对应的模糊数为：

${\hat{n}}_1=\underset{k=2}{\overset{K}{\mathrm{\Σ}}}{\widehat{\mathrm{\ϵ}}}_{k,1}\×b_{k,1}\×\frac{{\mathrm{\γ}}_1}{{\mathrm{\Γ}}_k}\mathrm{mod}{\mathrm{\γ}}_1$

（14）

${\hat{n}}_k=\frac{{\hat{n}}_1\×{\mathrm{\Γ}}_1-{\hat{q}}_{k,1}}{{\mathrm{\Γ}}_k},k=\mathrm{1,2},...K$

公式（14）中，

b_k,1为γ₁/Γ_k关于模Γ_k的乘法逆元。

步骤208：将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中进而确定出所述同余方程组的高程。

具体的，将步骤205中求得的解以及与该同余方程组对应的模糊数（如公式（14）所示）代入到原来的同余方程组（如公式（5）所示），求得各个同余方程的解。优选地，并取各个同余方程中高程的平均值作为多通道InSAR高程的估计值。该平均高程表示为：

$\hat{h}=\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}({\hat{n}}_k\×h_k+{\hat{r}}_k)---\left(15\right)$

下面结合图3至图5对本发明的多通道InSAR高程重建方法做进一步描述。

图3为利用仿真数据模拟的高程示意图，如图3所示，右图为三维立体图，不同的形状具有不同的高度；左图为右图的俯视图，通过左图灰度的明暗程度可以看出实际情况下不同形状的相对于基准面的高度。

图4为利用干涉相位与高程之间的关系仿真的加噪声后的缠绕干涉相位图，该图为三个干涉通道所得到的干涉相位图。从图中可以明显看出图像受到噪声的干扰，显示不清晰；要想得到实际情况下不同形状的相对于基准面的高度，需要综合考虑三幅干涉相位图。

图5为利用本发明提供的方法得到的高程结果示意图，如图5所示，通过该图灰度的明暗程度可以看出实际情况下不同形状的相对于基准面的高度，与图3所示的高程示意图极为接近，因此本发明提供的多通道InSAR高程重建方法具有较高的精确度。

本发明实施例还提供了一种多通道InSAR高程重建装置，如图6a所示，该多通道InSAR高程重建装置包括：第一确定单元61、建模单元62、获取单元63以及第二确定单元64；其中，

所述第一确定单元61，用于通过设置于测量航天器上的至少两个测量雷达对所述待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系；

所述建模单元62，用于根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为待测高程建模同余方程组；

所述获取单元63，用于对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数；

所述第二确定单元64，用于根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。

优选地，所述获取单元63，还用于根据鲁棒性中国剩余定理，求解变换后的同余方程组。

如图6b所示，所述建模单元62包括：第一确定子单元621、第二确定子单元622和第一变换子单元623；其中，

所述第一确定子单元621，用于确定所述待测高程的缠绕干涉相位对应的斜距差；

所述第二确定子单元622，用于根据所述斜距差、以及所述待测高程与所述至少两个测量雷达之间的位置关系，确定所述待测高程与干涉相位之间的关系；

所述第一变换子单元623，根据设定的余数，对所述待测高程与干涉相位之间的关系进行等价变换，得到所述待测高程的同余方程组。

如图6c所示，所述获取单元63包括：选取子单元631和第二变换子单元632；其中，

所述选取子单元631，用于从同余方程组中选取同余方程作为参考方程；

所述第二等价变换子单元332，用于将同余方程组中除参考方程以外的同余方程与所述参考方程进行相减，并将相减后的所有方程等价变化为一组同余方程组。

如图6d所示，所述第二确定单元64包括：回调子单元641和第三确定子单元642；其中，

所述回调子单元641，用于将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中；

所述第三确定子单元642，用于在将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中之后，确定出所述同余方程组的高程。

以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。

Claims (10)

1.一种多通道干涉合成孔径雷达InSAR高程重建方法，其特征在于，通过设置于测量航天器上的至少两个测量雷达对待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系；该方法包括：

根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为所述待测高程建模同余方程组；

对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数；

根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。

2.根据权利要求1所述的多通道InSAR高程重建方法，其特征在于，所述根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为待测高程建模同余方程组，包括：

确定所述待测高程的缠绕干涉相位对应的斜距差；

根据所述斜距差、以及所述待测高程与所述至少两个测量雷达之间的位置关系，确定所述待测高程与干涉相位之间的关系；

根据设定的余数，对所述待测高程与干涉相位之间的关系进行等价变换，得到所述待测高程的同余方程组。

3.根据权利要求1所述的多通道InSAR高程重建方法，其特征在于，所述对所述同余方程组进行等价变换包括：

从同余方程组中选取同余方程作为参考方程；

将同余方程组中除参考方程以外的同余方程与所述参考方程进行相减，并将相减后的所有方程等价变化为一组同余方程组。

4.根据权利要求1所述的多通道InSAR高程重建方法，其特征在于，所述求解变换后的同余方程组包括：根据鲁棒性中国剩余定理，求解变换后的同余方程组。

5.根据权利要求1至4任一项所述的多通道InSAR高程重建方法，其特征在于，根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程，包括：

将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中而确定出所述同余方程组的高程。

6.一种多通道InSAR高程重建装置，其特征在于，该装置包括：第一确定单元、建模单元、获取单元以及第二确定单元；其中，

所述第一确定单元，用于通过设置于测量航天器上的至少两个测量雷达对待测高程所在的位置分别发射电磁波，接收各发射电磁波的反射电磁波，并确定反射电磁波之间的干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系；

所述建模单元，用于根据所述干涉相位、以及所述待测高程与所述至少两个测量雷达之间的位置关系为所述待测高程建模同余方程组；

所述获取单元，用于对所述同余方程组进行等价变换，并求解变换后的同余方程组，获取所述变换后的同余方程组的解及所述变换后的同余方程组对应的模糊数；

所述第二确定单元，用于根据所述变换后的同余方程组的解及所述模糊数，确定出所述同余方程组的高程。

7.根据权利要求6所述的多通道InSAR高程重建装置，其特征在于，所述建模单元包括：第一确定子单元、第二确定子单元和第一变换子单元；其中，

所述第一确定子单元，用于确定所述待测高程的缠绕干涉相位对应的斜距差；

所述第二确定子单元，用于根据所述斜距差、以及所述待测高程与所述至少两个测量雷达之间的位置关系，确定所述待测高程与干涉相位之间的关系；

所述第一变换子单元，根据设定的余数，对所述待测高程与干涉相位之间的关系进行等价变换，得到所述待测高程的同余方程组。

8.根据权利要求6所述的多通道InSAR高程重建装置，其特征在于，所述获取单元包括：选取子单元和第二变换子单元；其中，

所述选取子单元，用于从同余方程组中选取同余方程作为参考方程；

所述第二等价变换子单元，用于将同余方程组中除参考方程以外的同余方程与所述参考方程进行相减，并将相减后的所有方程等价变换为一组同余方程组。

9.根据权利要求6所述的多通道InSAR高程重建装置，其特征在于，所述获取单元，还用于根据鲁棒性中国剩余定理，求解变换后的同余方程组。

10.根据权利要求6至9任一项所述的多通道InSAR高程重建装置，其特征在于，所述第二确定单元包括：回调子单元和第三确定子单元；其中，

所述回调子单元，用于将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中；

所述第三确定子单元，用于在将所述变换后的同余方程组的解以及所述变换后的同余方程组对应的模糊数代入到所述同余方程组中之后，确定出所述同余方程组的高程。

Images