[发明专利]冰川探测雷达成像中杂波抑制的方法

摘要

本发明提供了一种冰川探测雷达成像中杂波抑制的方法。该方法采用最大似然估计算法来估计冰层表面杂波，从实际雷达回波数据中去除表面杂波影响，不仅能够有效抑制杂波，同时可以得到冰层表面回波信号的估计，是一种非常具有应用价值的杂波抑制算法。

主权项

1.一种冰川探测雷达成像中杂波抑制的方法，其特征在于，包括：步骤A，向待测冰川发射探测雷达波，由多通道正下视雷达跨航迹向分布的P个接收通道分别接收雷达回波数据C_i(t，η)，其中，i=1、2、……，P，t和η分别表示距离向快时间和方位向慢时间，其中P>2；步骤B，对于P个接收通道中每个接收通道接收的雷达回波数据C_i(t，η)，分别对其进行距离向和方位向处理，得到SAR图像D_i(r_l，a_m)，r_l表示第l个距离向单元，a_m表示第m个方位向单元；步骤D，设定冰表面是平表面，利用雷达高度H和雷达正下方冰层深度d，获取与雷达波穿透到雷达正下方冰层深度为d的地方经历时延相等的冰层左右表面杂波的波达角θ₁和θ₂；步骤E，将冰层左右表面杂波的波达角θ₁和θ₂代入下式，得到阵列控制矩阵A(Θ)：A(Θ)=ej2kb1sinθ1ej2kb1sinθ2ej2kb2sinθ1ej2kb2sinθ2······ej2kbPsinθ1ej2kbPsinθ2]]>其中，k为波数2π/λ，λ为雷达波的波长，b_i是接收通道i的基线长度；步骤F，提取每个SAR图像中同一方位向和距离向单元内的像素值形成雷达接收信号矢量x(l，m)：x(l，m)=[D₁(r_l，a_m)，D₂(r_l，a_m)，…，D_p(r_l，a_m)]步骤G，利用阵列控制矩阵A(Θ)和雷达接收信号矢量x(l，m)，得到冰层表面杂波s^(l,m)=(AH(Θ)A(Θ))-1AH(Θ)x(l,m)]]>其中，冰层表面杂波为一矢量值，包含天底点左边和右边的冰层表面杂波信号；以及步骤H，在P个接收通道中取某一个接收通道作为基准接收通道，从该基准接收通道对应的SAR图像中减去冰层表面杂波得到杂波抑制后的雷达回波数据。