[发明专利]一种合成孔径雷达旁瓣抑制方法

权利要求书

1.一种合成孔径雷达旁瓣抑制方法，包括以下步骤：

步骤一、获得SAR图像的大小和频谱，SAR图像距离向像素点数目为l_r，SAR图像方位向像素点数目为l_a，SAR图像距离向过采样率为f_sr，SAR图像方位向过采样率为f_sa；

步骤二、SAR图像距离向旁瓣抑制；步骤二包括以下步骤：

步骤2-1、按照公式(1)计算得到距离向加权窗函数门限：

式(1)中：K_r为对f_sr进行四舍五入运算得到的整数；

步骤2-2：将SAR图像沿着距离向看做f_sa个f_sr×1的一维数组，记作date_r(n)，其中，r(r＝1,2,3,L,f_sa)表示二维SAR图像距离向第r组像素点；

步骤2-3：从date_r(n)中第K_r+1个像素点开始迭代，每个像素点的实部和虚部均分别按照公式(2)计算距离向窗函数加权参数，

式(2)中：w_r(n)表示第n(n＝K_r+1,K_r+2,K_r+3,L)次迭代后当前窗函数加权参数，n表示迭代次数，w_r(n)的约束条件为：-w_rmax≤w_r(n)≤w_rmax；

步骤2-4：date_r(n)中每个像素点的实部和虚部均分别按照公式(3)取值范围判别条件得到经过旁瓣抑制后的SAR图像距离向像素点像素值：

公式(3)中，date'_r(n)表示第n(n＝K_r+1,K_r+2,K_r+3,L)次迭代后当前SAR图像距离向像素点像素值；

步骤2-5：所述步骤2-2中r分别取值1至f_sa，每个r值对应的date_r(n)循环重复步骤2-3以及步骤2-4中的迭代操作，直到迭代到SAR图像距离向第(f_sr-K_r)个像素点迭代结束，得到经过距离向旁瓣抑制后的SAR图像date'r(x,y)。

2.如权利要求1所述的合成孔径雷达旁瓣抑制方法，其特征在于，所述步骤一具体为：

步骤1-1、分别获取二维SAR图像date(x,y)距离向和方位向像素点数目，分别记作l_r、l_a；

步骤1-2、对SAR图像date(x,y)做二维傅里叶变换，得到SAR图像频谱Date(ω_x,ω_y)，即：

步骤1-3、分别根据频谱Date(ω_x,ω_y)距离向、方位向占空比，其倒数即可分别得到SAR图像在距离向和方位向过采样率，分别记作f_sr、f_sa。

3.如权利要求1或2所述的合成孔径雷达旁瓣抑制方法，其特征在于，还包括步骤三、SAR图像方位向旁瓣抑制；所述步骤三具体包括以下步骤：

步骤3-1：按照公式(4)计算得到方位向加权窗函数门限：

式(4)中：K_a为对f_sa进行四舍五入运算得到的整数；

步骤3-2：将SAR图像沿着方位向看做f_sr个f_sa×1的一维数组，记作date_a(m)，其中，a(a＝1,2,3,L,f_sr)表示二维SAR图像方位向第a组像素点；

步骤3-3：从date_a(m)中SAR图像方位向第K_a+1素点开始迭代，每个像素点的实部和虚部均分别按照公式(5)计算方位向窗函数加权参数，

式(5)中，w_a(m)表示第m(m＝K_a+1,K_a+2,K_a+3,L)次迭代后当前窗函数加权参数，m表示迭代次数，w_a(m)的约束条件为：-w_amax≤w_a(n)≤w_amax；

步骤3-4：按照公式(6)取值范围判别条件得到经过旁瓣抑制后的SAR图像方位向像素点像素值：

式(6)中，date'_a(m)表示第m(m＝K_a+1,K_a+2,K_a+3,L)次迭代后当前SAR图像方位向像素点像素值；

步骤3-5：步骤2-2中a分别取值1至f_sr，每个a值对应的date_a(m)循环重复步骤3-3以及步骤3-4中的迭代操作，直到迭代到SAR图像方位向第(f_sa-K_a)个像素点迭代结束，得到经过方位向旁瓣抑制后的SAR图像date'a(x,y)。