一种SAR图像的旁瓣抑制方法、装置、电子设备及存储介质制造方法及图纸

本发明专利技术实施例涉及图像处理技术领域，特别涉及一种SAR图像的旁瓣抑制方法、装置、电子设备及存储介质。该方法包括：计算待处理的SAR图像的二维PSF估计矩阵；基于所述二维PSF矩阵和所述待处理的SAR图像的相关参数，求解解卷积函数，输出解卷积结果；其中，所述解卷积结果可还原为旁瓣抑制后的SAR图像，所述相关参数包括所述待处理的SAR图像的方位向采样点数、距离向采样点数、相位误差系数、模糊核大小、惩罚参数、平衡常数和最大迭代次数，所述方位向为雷达飞行轨道的方向，所述距离向为沿着雷达波发射的方向。本申请提供的SAR图像的旁瓣抑制方法旁瓣抑制性能好，并且能够有效滤除噪声，提高图像的质量。提高图像的质量。提高图像的质量。

【技术实现步骤摘要】
一种SAR图像的旁瓣抑制方法、装置、电子设备及存储介质


[0001]本专利技术涉及图像处理
，特别涉及一种SAR图像的旁瓣抑制方法、装置、电子设备及存储介质。

技术介绍

[0002]旁瓣抑制问题是传统合成孔径雷达图像处理中的一个重要问题。合成孔径雷达通过距离向脉冲压缩技术和方位向合成孔径技术获取二维高分辨率图像。但是，由于SAR系统接收的回波信号的频谱支撑域是有限分布的，SAR图像的点扩散函数(Point Spread Function，PSF)在距离向和方位向都是Sinc函数，存在不可忽略的旁瓣。
[0003]在SAR图像中，强目标，如建筑物和金属物体，不仅具有较强的主瓣能量，其旁瓣能量同样很高，极容易超过周围的弱目标的主瓣能量。弱目标被掩盖在强目标的旁瓣之下，从而无法识别。因此，需要对SAR图像进行旁瓣抑制，从而提升SAR图像的图像质量和视觉可读性，以便于后续的图像解译工作。
[0004]现有技术中，常用的SAR图像旁瓣抑制方法主要有频域加权方法、切趾滤波方法和空间变迹方法(Spatially Variant Apod本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种SAR图像的旁瓣抑制方法，其特征在于，包括：计算待处理的SAR图像的二维PSF估计矩阵；基于所述二维PSF估计矩阵和所述待处理的SAR图像的相关参数，求解解卷积函数，输出解卷积结果；其中，所述解卷积结果可还原为旁瓣抑制后的SAR图像，所述相关参数包括所述待处理的SAR图像的方位向采样点数、距离向采样点数、相位误差系数、模糊核大小、惩罚参数、平衡常数和最大迭代次数，所述方位向为雷达飞行轨道的方向，所述距离向为沿着雷达波发射的方向。2.根据权利要求1所述的方法，其特征在于，所述计算待处理的SAR图像的二维PSF估计矩阵，包括：获取所述待处理的SAR图像的散焦图像；基于所述待处理的SAR图像和所述待处理的SAR图像的散焦图像，计算所述待处理的SAR图像的第一因子；基于所述待处理的SAR图像和所述待处理的SAR图像的散焦图像，计算所述待处理的SAR图像的第二因子；基于所述第一因子和所述第二因子，计算所述待处理的SAR图像的辅助因子；对所述待处理的SAR图像的辅助因子做二维傅里叶逆变换，获得待处理的SAR图像的二维PSF估计矩阵。3.根据权利要求2所述的方法，其特征在于，所述获取所述待处理的SAR图像的散焦图像，包括：读入所述待处理的SAR图像和所述待处理的SAR图像的相关参数；其中，所述相关参数包括所述待处理的SAR图像的方位向采样点数N
a
、距离向采样点数N
r
、相位误差系数[Ep
a
,Ep
r
]、模糊核大小[h
r
,h
a
]、惩罚参数μ、平衡常数λ和最大迭代次数Iter
max
；对所述待处理的SAR图像进行二维傅里叶变换，得到第一频域图像；对所述第一频域图像添加相位误差，得到包含相位误差的第二频域图像；对所述第二频域图像进行二维傅里叶逆变换，得到所述待处理的SAR图像的散焦图像。4.根据权利要求3所述的方法，其特征在于，所述基于所述待处理的SAR图像和所述待处理的SAR图像的散焦图像，计算所述待处理的SAR图像的第一因子，包括：分别对所述待处理的SAR图像和所述待处理的SAR图像的散焦图像作距离向傅里叶变换，得到第一图像和第二图像；其中，所述第一图像中的每一行向量可以表示为所述第二图像中的每一行向量可以表示为式中，k＝1,2,
…
,N
a
、l＝1,2,
…
,N
r
；分别计算所述第一图像中每一行向量的辅助因子得到所述待处理的SAR图像的第一因子A(k,l)；其中，图像的第一因子A(k,l)；其中，5.根据权利要求4所述的方法，其特征在于，所述基于所述待处理的SAR图像和所述待
处理的SAR图像的散焦图像，计算所述待处理的SAR图像的第二因子，包括：分别对所述待处理的SAR图像和所述待处理的SAR图像的散焦图像作方位向傅里叶变换，得到第三图像和第四图像；其中，所述第三图像中的每一行向量可以表示为所述第四图像中的每一行向量可以表示为式中，k＝1,2,
…
,N
a
、l＝1,2,
…
,N
r
；分别计算所述第三图像中每一行向量的辅助因子得到所述待处理的SAR图像的第二因子B(k,l)；其中，图像的第二因子B(k,l)；其中，6.根据权利要求5所述的方法，其特征在于，所述基于所述第一因子和所述第二因子，计算所述待处理的SAR图像的辅助因子，包括：联立所述第一因子A(k,l)和所述第二因子B(k,l)，获得第一方程：A(k,l)a(k)
‑
B(k,...

【专利技术属性】
技术研发人员：杨威，刘艳芳，陈杰，李春升，曾虹程，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：