[发明专利]一种特殊地形条件下SAR干涉测量形变量解析方法

说明书

本发明涉及遥感影像处理领域，公开了一种特殊地形条件下SAR干涉测量形变量解析方法，该方法步骤如下：步骤1：对前后时相的SAR影像进行配准处理；步骤2：对配准后的影像对进行差分干涉处理；步骤3：对差分干涉处理后得到的相位差进行相位处理；将差分相位转换为视线方向的形变量；步骤4：对坡向视线角α进行计算。坡向视线角定义为点位梯度方向与雷达视线方向的夹角；步骤5：根据坡向视线角及形变量进行判断解析；步骤6：得到形变量解析结果。本发明提供的解析方法能实现特殊地形条件下，SAR干涉测量形变量的正确解析，对城市地面沉降监测有重要意义。

技术领域

本发明涉及遥感影像处理领域，具体说是一种用于特殊地形条件下SAR干涉测量形变量的解析方法。

背景技术

现有的SAR干涉测量形变量解析方法通常不考虑地面的地形条件。其通常采用前后时相的SAR影像进行配准、干涉处理、相位处理等步骤获取沿视线方向的点位形变量，然后根据其符号来判断点位的升降情况。但在特殊地形条件下，当地形梯度方向与雷达视线方向夹角为锐角时，传统的判断方法与真实情况呈现相悖的状态，使得传统方法应用局限极大。

发明内容

本发明所要解决的技术问题是：针对上述存在的问题，提供了一种用于特殊地形条件下的SAR干涉测量形变量解析方法。该方法在利用SAR干涉测量获得形变量的基础上，通过对梯度方向与雷达视线方向的位置关系的分析，实现点位沉降或抬升情况的判断。该发明克服了传统方法在特殊地形条件下不能正确判断点位升降情况的缺陷，具有更广泛的适用性。

本申请提供一种用于特殊地形条件下SAR干涉测量形变量的解析方法，采用的技术方案如下：

步骤1：对前后时相的SAR影像进行配准处理；

步骤2：对配准后得到的影像对进行差分干涉处理；差分干涉处理包括两个部分，首先对影像复数数据做干涉运算，然后将干涉相位与地形相位做差值运算；

步骤3：对差分干涉处理后得到的相位差进行相位处理，将差分相位转换为视线方向的形变量；

步骤4：对坡向视线角进行计算；其中，所述坡向视线角为点位梯度方向和雷达视线方向的夹角；

坡向视线角α的计算方式如下：

其中，α为坡向视线角，为点位的下视角，为点位切线的坡角。

步骤5：根据坡向视线角及形变量进行判断解析；

当α≥90°,d＜0或α＜90°,d＞0时，形变为沉降其中，d为形变量；当α≥90°,d＞0或α＜90°,d＜0时，形变为抬升；当d＝0时，不存在相对形变；

步骤6：得到形变量的分析结果。

进一步的，配准处理精度要达到亚像素级别，配准处理是进行干涉测量的基础，所以需确保其精度。

进一步的，所述步骤3中的相位处理包括：相位滤波、相位截仓和相高转换，其中，相位滤波采用自适应滤波方法，相位解缠使用最小费用流方法，通过相高转换后得到形变量。

与现有技术相比，采用上述技术方案的有益效果为：本发明能实现特殊地形条件下，SAR干涉测量形变量的正确解析，对城市地面沉降监测有重要意义。

附图说明

图1是本发明实施例提供的一种用于特殊地形条件下SAR干涉测量形变量的解析方法的流程示意图。

具体实施方式

下面结合附图对本发明做进一步描述。

实施例1

如图1所示，本发明实施例提供一种用于特殊地形条件下SAR干涉测量形变量的解析方法，该方法步骤包括以下步骤：