[发明专利]一种基于图像后处理的太赫兹SAR图像副瓣抑制方法

权利要求书

1.一种基于图像后处理的太赫兹SAR图像副瓣抑制方法，其特征在于，包括：设SAR图像为I，距离向尺寸为N_r，即图像的行数，方位向尺寸为N_a，即图像的列数；若强点的高副瓣贯穿SAR图像的整个方位向，则选取一定方位宽度的清晰区域，计算噪底数值和每列像素的数值抑制函数，生成数值抑制的模板图像，抑制方位高副瓣；若强点的高副瓣延伸到有限距离范围，则计算图像每个强点的数值抑制函数，生成数值抑制的模板图像，抑制强点高副瓣。

2.根据权利要求1所述的基于图像后处理的太赫兹SAR图像副瓣抑制方法，其特征在于，所述抑制方位高副瓣，包括：在SAR图像清晰区域选择一块区域，在方位向求均值，得到距离向的一维向量，将一维向量的中位数作为噪底数值；将噪底数值的1.2倍作为阈值，将清晰区域中大于该阈值的距离门位置作为需要高副瓣抑制的位置，选取每列像素前后一定宽度的区域，在方位向求均值，作为该列像素的投影向量，删除其需要高副瓣抑制的位置对应的像素，重新插值，得到与投影向量长度相同的向量，将两个向量的对应元素相除，得到抑制倍数向量；生成与SAR图像尺寸相同的抑制模板图像，根据抑制倍数向量，采用钟形函数定义，得到每个需要高副瓣抑制的位置对应的钟形函数，取两者对应元素的最小值，向抑制模板图像赋值；将抑制模板图像与SAR图像的对应元素相乘，得到方位高副瓣抑制后的SAR图像。

3.根据权利要求2所述的基于图像后处理的太赫兹SAR图像副瓣抑制方法，其特征在于，所述计算噪底数值，包括：在图像I的一侧提取行数为N_r、列数为N_{a_part}的区域I_part，在方位向求均值，得到长度为N_r的一维向量V_part，作为清晰区域的投影向量，用中位数函数mid(·)计算中位数T_noise＝mid(V_part)，作为SAR图像的噪底数值。

4.根据权利要求3所述的基于图像后处理的太赫兹SAR图像副瓣抑制方法，其特征在于，所述计算抑制倍数向量，包括：设置高副瓣区域的抑制阈值为T_lobe＝1.2T_noise，V_part中数值大于T_lobe对应的距离门位置定义为高副瓣区域，得到高副瓣的距离位置信息向量S_lobe，不同方位用不同倍数抑制SAR图像中距离门在S_lobe处的数值。

5.根据权利要求4所述的基于图像后处理的太赫兹SAR图像副瓣抑制方法，其特征在于，所述计算抑制倍数向量，还包括：设图像I的第i列像素为I(:,i)，i的取值范围为[N_p,N_a-N_p]，在位置i的两边各扩展N_p个像素，得到距离向尺寸为N_r、方位向尺寸为i+2N_p的图像块I(:,i-N_p:i+N_p)，在其方位向求均值，得到长度为N_r的距离向向量V_r，作为列像素I(:,i)的投影向量，删除V_r中的高副瓣位置S_lobe，采用线性插值方法填充原位置，保持元素数量与V_r相同，得到修正后的一维向量V′_r，用对应元素相除./将V_r和V′_r相除N_depress＝V_r./V′_r，得到I(:,i)处的抑制倍数向量N_depress，计算图像I的杂波区域每列像素对应的抑制倍数向量。