[发明专利]干涉SAR叠掩区域数字高程模型重建的方法

权利要求书

1.一种干涉SAR叠掩区域数字高程模型重建的方法，包括：

检测由第一通道获取的第一干涉SAR复图像和由第二通道获取的第二干涉SAR复图像的叠掩区域；

将第一干涉SAR复图像和第二干涉SAR复图像进行共轭相乘，对共轭相乘的结果取相位，从而得到干涉相位图；

将叠掩区域作为掩膜，对干涉相位图上除叠掩区域外的非叠掩区域进行干涉相位解缠，得到非叠掩区域的解缠相位；

在每一方位向上，由所述非叠掩区域的解缠相位，分别得到叠掩区域左侧边界的非叠掩区域的第一解缠相位和右侧边界的非叠掩区域的第三解缠相位

在每一方位向，由所述非叠掩区域的解缠相位，计算叠掩区域的左侧沿距离向的前干涉相位频率f_before，和右侧沿距离向的后干涉相位频率f_after；

在每一方位向上，由第一解缠相位和第三解缠相位及前相位频率f_before和后相位频率f_after，计算出该方位向的叠掩区域左侧边界的解缠相位和右侧边界的解缠相位

在每一方位向上，由叠掩区域的左侧边界的解缠相位和右侧边界的解缠相位来拟合该方位向上叠掩区域内部的每一像素的解缠相位

在每一方位向的每一像素，由该像素的解缠相位和干涉SAR参数计算该像素的数字高程模型h。

2.根据权利要求1所述的干涉SAR叠掩区域数字高程模型重建方法，其中，所述由非叠掩区域的解缠相位，计算叠掩区域的左侧沿距离向的前干涉相位频率f_before，和右侧沿距离向的后干涉相位频率f_after的步骤包括：

在叠掩区域的左侧和右侧各选取(2M+1)×(2N+1)的估计窗口；

通过最大化对应估计窗口的相位模型的代价函数来估计前干涉相位频率f_before和后干涉相位频率f_after。

3.根据权利要求2所述的干涉SAR叠掩区域数字高程模型重建方法，其中，所述M＝16，N＝16。

4.根据权利要求1所述的干涉SAR叠掩区域数字高程模型重建方法，其中，根据以下公式，计算出该方位向的叠掩区域左侧边界的解缠相位和右侧边界的解缠相位

其中，N_pixel为叠掩区域沿距离向的像素数。

5.根据权利要求4所述的干涉SAR叠掩区域数字高程模型重建方法，其中，根据以下公式，由叠掩区域的左侧边界的解缠相位和右侧边界的解缠相位来拟合该方位向上叠掩区域内部的解缠相位

其中，n为待计算相位的像素相对于叠掩区域的左边界沿距离向偏移的像素数，n＝0，1，…N_pixel-1。

6.根据权利要求1所述的干涉SAR叠掩区域数字高程模型重建方法，其中，根据以下公式，由像素的解缠相位和干涉SAR参数计算该像素的数字高程模型h，则：

θ=α+asin(R12+B2-R222R1B)]]>

h＝H-R₁cosθ

其中，R₁为第一通道的斜距，λ为所用雷达的波长，R₂为第二通道的斜距，α为基线角，B为基线长度，θ为第一通道的视角，h为反演出目标的DEM，H为第一通道的高度。