[发明专利]用于InSAR地形测量的方位向预滤波处理方法

权利要求书

1.用于InSAR地形测量的方位向预滤波处理方法，包括如下步骤：

（1）输入SAR图像数据；

（1a）输入干涉合成孔径雷达InSAR主天线获取的主图像数据；

（1b）输入干涉合成孔径雷达InSAR辅天线获取的已完成与主图像完全配准的辅图像数据；

（1c）输入干涉合成孔径雷达InSAR成像处理的处理参数和系统参数；

（2）复干涉图分块：

（2a）对输入的主、辅SAR图像进行干涉处理，得到复干涉图；

（2b）对复干涉图沿方位向进行分块处理，得到复干涉图分块；

（3）估计相位梯度：

采用相位梯度估计方法，分别对每个复干涉图分块的数据进行方位向干涉相位梯度估计，得到每个复干涉图分块各个距离单元的相位梯度估计值；

（4）计算各分块方位向预滤波处理参数：

（4a）利用步骤（3）得到的每个复干涉图分块各个距离单元的相位梯度估计值，计算每个复干涉图分块对应的主辅SAR图像区域的多普勒中心偏移量；

（4b）利用下式，估计预滤波后每个复干涉图分块对应的主、辅SAR图像区域数据的多普勒中心和多普勒带宽：

f′_d1＝(f_d1+Δf_d+f_d2)/2

f′_d2＝(f_d1-Δf_d+f_d2)/2

B′_a＝B_a-|f_d1-f_d2-Δf_d|

其中，f′_d1和f′_d2分别表示预滤波后复干涉图分块对应的主、辅SAR图像区域数据的多普勒中心，d表示多普勒中心，f_d1和f_d2分别表示预滤波前复干涉图分块对应的主、辅SAR图像区域数据的多普勒中心，Δf_d表示复干涉图分块对应的主辅SAR图像区域数据多普勒中心的偏移量，B′_a表示预滤波后复干涉图分块对应的主、辅SAR图像的多普勒带宽，a表示多普勒域，B_a表示预滤波前复干涉图分块对应的主、辅SAR图像区域数据的多普勒带宽，|·|表示取绝对值操作；

（5）多普勒域滤波：

（5a）采用傅里叶变换，将每一个复干涉图分块对应的主、辅SAR图像区域的数据，变换至距离多普勒域，获得距离多普勒的SAR数据；

（5b）采用去除窗函数的方法，获得除去窗函数后的距离多普勒的SAR数据；

（5c）利用步骤（4b）获取的预滤波参数，对每一复干涉图分块对应的主、辅SAR图像区域的数据沿距离向分别构造窗函数，进行滤波处理；

（5d）采用傅里叶逆变换，将滤波后的数据变换回二维时域；

（6）输出方位预滤波处理后的主、辅SAR图像数据。

2.根据权利要求1所述的用于InSAR地形测量的方位向预滤波处理方法，其特征在于：步骤（2a）所述的干涉处理是利用下式实现：

I＝R·*conj(S)

其中，I表示主、辅SAR图像的干涉图，R表示主图像数据，S表示辅图像数据，conj(·)表示取共轭操作。

3.根据权利要求1所述的用于InSAR地形测量的方位向预滤波处理方法，其特征在于：步骤（2b）所述的分块处理是指，将复干涉图从起始方位时刻开始，以2的整数次幂点的分块大小对复干涉图进行分块。