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

摘要

本发明公开了一种用于InSAR地形测量的方位向预滤波处理方法，本发明的实现步骤是：（1）输入SAR图像数据；（2）生成复干涉图并对干涉图进行分块；（3）估计相位梯度；（4）计算各分块方位向预滤波处理参数；（5）多普勒域滤波；（6）输出方位预滤波处理后的主、辅SAR图像数据。本发明具有在长基线情况下，能够对配准后的主、辅SAR图像完成高质量的方位向预滤波处理，满足高质量干涉合成孔径雷达InSAR处理的实际工程性能要求，获得高质量的测绘产品。

主权项

用于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)/2f′_d2＝(f_d1‑Δf_d+f_d2)/2B′_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图像数据。