一种基于互质多基线的高程重建方法及装置

摘要

本发明公开了一种基于互质多基线的高程重建方法，包括：基于干涉合成孔径雷达InSAR系统给定参数，设计各天线间形成的基线长度为互质值；根据观测干涉相位的概率密度函数建立最大似然函数，并利用最大似然估计算法确定地形高程。本发明还同时公开了一种基于互质多基线的高程重建装置，采用本发明的技术方案，能够使高度不连续地区的高程重建具有唯一性，并提高了测量精度。

主权项

一种基于互质多基线的高程重建方法，其特征在于，所述方法包括：基于干涉合成孔径雷达InSAR系统给定参数，使各天线间形成的基线长度为互质值，从而确定观测干涉相位；根据所述观测干涉相位的概率密度函数建立最大似然函数，并利用最大似然估计算法确定地形高程；其中，所述基于InSAR系统给定参数，使各天线间形成的基线长度为互质值，从而确定观测干涉相位，包括：依据互质性要求将所述各天线进行放置，并使各天线间形成的基线长度的互质精度在设定阈值范围内；将所述各天线所在通道所采集的回波数据及直通信号进行同步，并形成合成孔径雷达SAR图像；获得SAR图像后，以其中一幅SAR图像作为基准进行配准，并形成干涉图，通过去平地效应处理、干涉图滤波处理，得到观测干涉相位；其中，所述根据观测干涉相位的概率密度函数建立最大似然函数，利用最大似然估计算法确定地形高程，包括：分别确定每幅干涉图的高度模糊数h_ak，其中，k＝1,...,K，K表示干涉图的总幅数；确定h_ak两两互质时的最大高度模糊数H_m，其中H_m为：确定观测干涉相位的概率密度函数f(φ_k|h)，其中，所述观测干涉相位的概率密度函数f(φ_k|h)为：$f\left(\phi \right(i,j\left)\right|h\left(i,j\right))=\frac{1}{2\pi }\frac{1-|\gamma {|}^2}{1-|\gamma {|}^2{\cos }^2\left(\phi \right(i,j)-\frac{2\pi }{h_a}h(i,j\left)\right)}·\{1+\frac{\left|\gamma \right|cos\left(\phi \right(i,j)-\frac{2\pi }{h_a}h(i,j\left)\right){\cos }^{-1}[-|\gamma \left|cos\left(\phi \right(i,j)-\frac{2\pi }{h_a}h(i,j\left)\right)\right]}{{[1-|\gamma {|}^2{\cos }^2\left(\phi \right(i,j)-\frac{2\pi }{h_a}h(i,j\left)\right)]}^{1/2}}\},$φ(i,j)表示干涉缠绕相位信号，φ(i,j)∈(‑π,π]，h(i,j)表示地表高程，γ为相干系数，h_a表示高度模糊数；将使似然函数为最大值时的h确定为地形高程。