买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：电子科技大学

申请日：2024-10-30

公开（公告）日：2024-12-13

公开（公告）号：CN119126033A

专利技术分类：..抗干扰装置[2006.01]

专利摘要：本发明公开了一种基于可学习时频特征分离的SAR干扰抑制方法，首先建立基于可学习时频特征分离的SAR干扰抑制问题模型，在低秩稀疏分解模型框架中引入可学习可逆非线性变换，然后构建SAR干扰抑制网络并进行训练，将待处理的含干扰的SAR回波进行时频变换并裁剪后输入训练好的网络，输出SAR回波时频谱和干扰时频图，再将网络输出的SAR回波时频谱按裁剪方式重新拼接为完整时频数据，得到干扰抑制后的SAR回波时频谱，经逆短时傅里叶变换恢复成时域回波，最后将干扰抑制后的时域回波进行SAR成像处理，得到干扰抑制后的SAR图像。本发明的方法解决了在干扰环境下SAR图像质量受损或无法成像的问题，在时频域和图像域的干扰抑制和回波保真方面提供了可靠的性能。

专利权项：1.一种基于可学习时频特征分离的SAR干扰抑制方法，具体步骤如下：S1、建立基于可学习时频特征分离的SAR干扰抑制问题模型，在低秩稀疏分解模型框架中引入可学习可逆非线性变换，增强干扰和回波的特征表征能力；根据SAR回波和干扰的时频域分布特性，设置基于低秩和稀疏分解的干扰抑制问题基本模型，表达式如下： 其中，||·||*表示核范数，||·||1表示l1范数，λ表示加权因子，L表示低秩成分，S表示稀疏成分，X表示待处理的数据；然后在低秩稀疏分解模型基础上，引入两个可学习非线性变换和分别增强SAR回波和干扰的特征表征，并将SAR干扰抑制问题建模为基于可学习变换的低秩稀疏分解模型，表达式如下： 再构造增广拉格朗日函数求解式，表达式如下： 其中，ρ表示惩罚因子，Λ表示拉格朗日乘子，·表示内积，即矩阵对应元素乘积求和；再采用ADMM算法求解，将式3分解为三个子问题，每一步只更新一个变量而固定其余变量，交替重复更新；对于迭代次数k＝1,2,3,...，三个迭代子问题表达式如下： S2、基于步骤S1，构建基于可学习时频特征分离的SAR干扰抑制网络，包括：整体迭代算法的深度展开网络和可逆非线性变换学习网络；S21、整体迭代算法的深度展开网络；将求解式4的迭代过程展开为包含P个迭代层的深度学习的方案，3个子问题分别对应三个网络模块：LM、SM和ΛM，分别更新L，S和Λ，并将超参数ρ和λ设置为可学习的参数；其中，更新L，S和Λ对应的各网络模块具体如下：1第p层中低秩成分L的更新模块LMp表达式如下： 其中，表示增强干扰表征的可逆非线性变换，是其逆变换，p＝1,2,...,P表示展开网络的层数，对应于式4中的迭代次数k；表示核范数的近端算子，定义为W＝UΣVH即矩阵W的奇异值分解，U表示左奇异向量，V表示右奇异向量，H表示矩阵的共轭转置，Σ＝diagσ1,σ2,...,σr表示由奇异值组成的对角矩阵，r表示矩阵W的秩；用于对超过阈值的奇异值进行筛选；αp表示控制奇异值阈值的参数，被设计为可学习的网络参数，且随着迭代次数增加，αp逐层减小，则αp具体表达式如下： 其中，ωα表示可学习的系数因子，θα用于控制αp随层数p的变化程度，设置θα＞1来约束αp逐层减小；2第p层中稀疏成分S的更新模块SMp表达式如下： 其中，表示增强SAR回波时频表征的可逆非线性变换，是其逆变换；表示软阈值算子，控制软阈值的参数βp表达式如下： 其中，ωβ和θβ表示可学习的参数，θβ用于控制βp随层数p的变化程度，设置θβ＞1来约束βp随着p增加而递减；3第p层中拉格朗日乘子Λ的更新模块ΛMp表达式如下：Λp＝Λp-1+γpLp+Sp-X9其中，可学习的惩罚因子γp表达式如下： 其中，ωγ和θγ表示可学习的参数，指数p大于1，则当θγ＞1时，γp随着层数p的增长而增大；S22、可逆非线性变换学习网络；可逆非线性变换学习网络将输入分成两部分，通过两部分之间的交叉耦合来构造变换，且处理对象是复时频谱数据；首先使用重构算子将X的实部与虚部结合，将重组为具有两个通道的张量表达式如下： 其中，表示复数域，表示实数域，m表示离散的时间采样点，q表示离散的频率采样点，Nm表示时间采样点个数，Nq表示频率采样点个数；然后沿通道维将划分为两部分在正变换时执行前向过程，两路成分通过学习的函数实现交叉耦合；第i个耦合层的前向过程表达式如下： 其中，ψ·表示可逆下采样算子，用于将空间尺寸转化为通道维数，维持求和的两部分维度相等；耦合层的可逆性与和无关，采用卷积神经网络CNN模块，通过网络中的线性与非线性运算提取有效的深度特征；一个耦合层中对于X1,X2的操作并不对称，构建l个耦合层组合作为一个可逆变换网络模块，即可逆非线性变换学习网络；在耦合层之后，沿通道维合并两个输出成分，作为正变换的结果；其中，l3；在逆变换时执行反向过程，其中第i个耦合层的反向过程表达式如下： S3、将待处理的含干扰的SAR回波数据利用短时傅里叶变换逐脉冲转换至时频域，然后将时频域数据裁剪为步骤S2构建的SAR干扰抑制网络可处理尺寸；采用短时傅里叶变换STFT实现时频数据分析与处理，信号x[n]的局部时频特征Xtf[m,q]表达式如下： 其中，w[n]表示时域分析窗口，x[n]表示一个脉冲的待处理回波，n表示离散的时间变量，即分析窗口在时间轴的位置，R表示步长；在短时傅里叶变换后，将时频域数据裁剪为统一的网络可处理尺寸；S4、训练基于可学习时频特征分离的SAR干扰抑制网络；进行有监督训练，将不同场景SAR回波时频谱和不同类型、参数的干扰时频谱作为分离结果的真值来构建训练和验证数据集；对数据集中的干扰污染时频谱回波时频谱干扰时频谱进行幅度归一化处理；选择最大最小值归一化方法，且由于频率轴范围大于信号带宽，时频图中幅度最小值为0；为保证归一化后的三元组满足最大值采用的幅度最大值，则归一化操作表达式如下： 然后将归一化时频数据切割为统一尺寸的复数时频图，保存为三元组作为训练输入；网络输出为分离的SAR回波时频谱和干扰时频谱，且考虑网络输出的干扰和SAR回波时频谱的重构误差，定义损失函数表达式如下： 其中，L′和S′分别表示网络输出的干扰和SAR回波时频谱分解结果，σ表示调节权重的参数，||·||F表示矩阵的F范数；处理的数据均为复数据，F范数计算的是网络输出数据与标签间各像素复数值的差的模平方之和；S5、将待处理的含干扰的SAR回波按步骤S3进行时频变换并裁剪，然后将裁剪好的时频数据输入步骤S4训练好的SAR干扰抑制网络进行处理，得到网络输出的SAR回波时频谱和干扰时频图；S6、将步骤S5得到的网络输出的SAR回波时频谱按裁剪方式重新拼接为完整时频数据，得到干扰抑制后的SAR回波时频谱，经逆短时傅里叶变换恢复成时域回波，即干扰抑制后的时域回波；S7、将步骤S6干扰抑制后的时域回波进行SAR成像处理，得到干扰抑制后的SAR图像。

百度查询： 电子科技大学 一种基于可学习时频特征分离的SAR干扰抑制方法