发明名称 强杂波背景下目标信号的提取和超分辨率增强处理方法
摘要 本发明公开了一种强杂波背景下目标信号的提取和超分辨率增强处理方法,其主要思路是:得到地面场景中点目标的成像结果后,选择一个包括该点目标和地杂波的目标矩形区域切片数据,并提取该目标矩形区域切片数据中与傅里叶基函数相关度最大值对应的散射点二维位置对应的强散射点,采用梯度下降法对该强散射点邻近的弱散射点进行提取,得到满足梯度下降法条件的散射点二维位置对应的强散射点邻近的弱散射点区域,将目标矩形区域切片数据减去满足梯度下降法条件的散射点二维位置对应的强散射点邻近的弱散射点区域,进而得到一幅完整SAR图像,并对该幅完整SAR图像利用正则化方法进行超分辨率增强处理,得到一幅最终的超分辨率SAR图像。
申请公布号 CN104991241A 申请公布日期 2015.10.21
申请号 CN201510371744.9 申请日期 2015.06.30
申请人 西安电子科技大学 发明人 孙光才;景国彬;盛佳恋;邢孟道;保铮
分类号 G01S7/48(2006.01)I 主分类号 G01S7/48(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种强杂波背景下目标信号的提取和超分辨率增强处理方法,其特征在于,包括以下步骤:步骤1,对地面场景的SAR雷达回波信号依次进行运动补偿和常规线频调变标成像处理,得到第p个点目标的成像结果<img file="FDA0000749035940000011.GIF" wi="174" he="78" />进而得到一幅包含P个点目标的成像结果Img;其中,p∈{1,2,…,P},P表示地面场景中的点目标总个数;步骤2,在一幅包含P个点目标的成像结果Img中,提取包含P个点目标和地杂波的目标矩形区域切片数据D<sub>a</sub>,分别得到该目标矩形区域切片数据D<sub>a</sub>的二维频域离散表达式S(n,m)和该目标矩形区域切片数据D<sub>a</sub>的时域表达式G(k,h);其中,n∈[1,N],m∈[1,M],k∈[1,N],h∈[1,M],N表示目标矩形区域切片数据D<sub>a</sub>的距离采样单元个数,M表示目标矩形区域切片数据D<sub>a</sub>的方位采样单元个数,D<sub>a</sub>表示包括P个点目标和地杂波的目标矩形区域切片数据,P表示地面场景中的点目标总个数;步骤3,将目标矩形区域切片数据D<sub>a</sub>的四个大小一致的角落杂波矩形区域与傅里叶基函数W<sub>kh</sub>(n,m)分别进行相关性处理,计算得到自适应杂波门限ε;其中,n∈[1,N],m∈[1,M],N表示目标矩形区域切片数据D<sub>a</sub>的距离采样单元个数,M表示目标矩形区域切片数据D<sub>a</sub>的方位采样单元个数,D<sub>a</sub>表示包括P个点目标和地杂波的目标矩形区域切片数据,P表示地面场景中的点目标总个数;步骤4,将目标矩形区域切片数据D<sub>a</sub>与傅里叶基函数W<sub>kh</sub>(n,m)进行相关性处理,提取得到目标矩形区域切片数据D<sub>a</sub>中与傅里叶基函数W<sub>kh</sub>(n,m)相关度最大值对应的散射点二维位置(k,h)<sup>q</sup>对应的强散射点G(k,h)<sup>q</sup>;其中,k∈[1,N],h∈[1,M],N表示目标矩形区域切片数据D<sub>a</sub>的距离采样单元个数,M表示目标矩形区域切片数据D<sub>a</sub>的方位采样单元个数,D<sub>a</sub>表示包括P个点目标和地杂波的目标矩形区域切片数据,P表示地面场景中的点目标总个数,q表示迭代次数;步骤5,采用梯度下降法对目标矩形区域切片数据D<sub>a</sub>与傅里叶基函数W<sub>kh</sub>(n,m)相关度最大值对应的散射点二维位置(k,h)<sup>q</sup>对应的强散射点G(k,h)<sup>q</sup>邻近的弱散射点区域进行提取,得到满足梯度下降法条件的强散射点G(k,h)<sup>q</sup>邻近的弱散射点区域<img file="FDA0000749035940000012.GIF" wi="199" he="90" />其中,n∈[1,N],m∈[1,M],k∈[1,N],h∈[1,M],<img file="FDA0000749035940000013.GIF" wi="412" he="78" />N表示目标区域切片数据D<sub>a</sub>的距离采样单元个数,M表示目标区域切片数据D<sub>a</sub>的方位采样单元个数,D<sub>a</sub>表示包括P个点目标和地杂波的目标区域切片数据,P表示地面场景中的点目标总个数,q表示迭代次数;步骤6,将目标矩形区域切片数据D<sub>a</sub>减去满足梯度下降法条件的强散射点<img file="FDA0000749035940000021.GIF" wi="170" he="72" />邻近的弱散射点区域<img file="FDA0000749035940000022.GIF" wi="199" he="85" />得到剩余目标矩形区域切片数据<img file="FDA0000749035940000023.GIF" wi="83" he="70" />并对该剩余目标矩形区域切片数据<img file="FDA0000749035940000024.GIF" wi="60" he="69" />利用步骤4和步骤5进行迭代操作,直到剩余目标矩形区域切片数据<img file="FDA0000749035940000025.GIF" wi="51" he="67" />中所有散射点与傅里叶基函数W<sub>kh</sub>(n,m)相关度最大值低于自适应杂波门限ε,迭代停止,得到P个强散射点和该P个强散射点各自邻近的弱散射区域,进而组合成一幅包含P个点目标的完整SAR图像G,且<img file="FDA0000749035940000026.GIF" wi="377" he="101" />其中,<img file="FDA0000749035940000027.GIF" wi="177" he="90" />表示满足梯度下降法条件的散射点二维位置<img file="FDA0000749035940000028.GIF" wi="105" he="92" />对应的所有散射点,<img file="FDA0000749035940000029.GIF" wi="415" he="78" />N表示目标矩形区域切片数据D<sub>a</sub>的距离采样单元个数,M表示目标矩形区域切片数据D<sub>a</sub>的方位采样单元个数,D<sub>a</sub>表示包括P个点目标和地杂波的目标矩形区域切片数据,P表示地面场景中的点目标总个数,P也表示目标矩形区域切片数据D<sub>a</sub>中的强散射点总个数,q表示迭代次数;步骤7,利用正则化方法对一幅包含P个点目标的完整SAR图像G进行超分辨率增强处理,得到一幅最终的超分辨率SAR图像<img file="FDA00007490359400000210.GIF" wi="66" he="66" />其中,P表示地面场景中的点目标总个数。
地址 710071 陕西省西安市太白南路2号