CN103064071A - 基于稀疏分解的雷达目标属性散射中心特征提取方法 - Google Patents
基于稀疏分解的雷达目标属性散射中心特征提取方法 Download PDFInfo
- Publication number
- CN103064071A CN103064071A CN201210414929XA CN201210414929A CN103064071A CN 103064071 A CN103064071 A CN 103064071A CN 201210414929X A CN201210414929X A CN 201210414929XA CN 201210414929 A CN201210414929 A CN 201210414929A CN 103064071 A CN103064071 A CN 103064071A
- Authority
- CN
- China
- Prior art keywords
- scattering center
- parameter
- dictionary
- scattering
- vector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种基于稀疏分解的雷达目标属性散射中心特征提取方法,主要解决现有基于对雷达图像进行图像分割提取属性散射中心方法的模型失配、特征易丢失以及参数估计精度较低的问题。其实现过程是:首先利用噪声样本建立散射中心强度门限,在雷达图像中进行强散射中心检测,并确定散射中心参数的取值集合;然后根据属性散射中心模型,利用坐标轮回下降技术构建超分辨字典,通过求解0范数优化问题得到目标属性散射中心参数超分辨估计集合,根据散射中心参数集合提取目标及其重要部件的几何尺寸特征。本发明能有效提取目标属性散射中心、超分辨估计散射中心参数、精确估计目标及其重要部件几何尺寸,可用于雷达目标分类识别。
Description
技术领域
本发明属于雷达技术领域,涉及一种雷达目标属性散射中心特征提取方法,可用于估计目标及其重要部件的几何尺寸,为目标分类识别提供重要的特征信息。
背景技术
雷达成像技术是20世纪50年代发展起来的,雷达图像是目标的二维散射图。传统雷达成像是以点散射模型为基础的,该模型只包含目标散射点位置信息,但仅利用目标散射点的位置信息构建的识别特征并不能完备表征雷达图像中目标的本质属性。在光学区,扩展目标的高频电磁散射响应可以用一组独立分布的散射体,或称散射中心的电磁散射响应之和近似表示。目标的散射中心主要产生于目标的边缘、拐点、棱角及尖端等不连续点部位,代表了目标的精细物理结构,所以散射中心模型能够更贴切地描述目标属性,也在雷达目标识别领域有着重要的应用。
基于几何绕射理论和物理光学理论,1999年Michael J.Gerry和Lee C.Potter提出了一个适用于合成孔径雷达(SAR)的参数化模型—属性散射中心模型,见[M.J.Gerry,L.C.Potter,I.J.Gupta,and A.van der Merwe,A parametric model for syntheticaperture radar measurements[J].IEEE Transactions on Antennas and Propagation,1999,Vol.47,NO.7,pp.1179-1188]。属性散射中心模型用一组参数描述每个散射中心的位置、形状、方向以及幅度等,这些属性都是关系目标的重要信息;同点散射模型相比,属性散射中心模型包含了更丰富的可用于目标分类识别的特征。
目标散射中心的特征提取本质上是一个从目标回波数据中估计各个散射中心参数的过程。由于属性散射中心模型结构复杂以及参数维数较高,增加了模型参数估计的复杂性。现有方法通过对以点散射模型为基础得到的雷达图像进行图像分割,得到阶数较低的目标散射区或者是孤立的散射中心,利用近似最大似然方法估计目标的属性散射中心参数。首先由于现有方法利用根据点散射模型得到的雷达图像提取属性散射中心,所以存在模型失配问题;其次由于该类方法是以图像分割为基础的,所以这类方法要求图像质量较高,此外当目标某些部件散射强度较弱时,通过图像分割方法是难以正确检测的,这就导致目标的一些重要特征容易丢失。除此以外由于该类方法的优化问题为非凸问题,并且存在很多局部最小解,所以存在参数初始化、模型阶数选择和散射中心结构类别判别的问题,使得最终参数估计精度较低。
发明内容
本发明的目的在于提出一种基于稀疏分解的雷达目标属性散射中心特征提取的方法,以解决现有方法中存在的模型失配、特征易丢失和参数估计精度较低问题。
本发明是这样实现的:
一.技术思路
雷达回波中,目标散射场绝大部分能量仅由少量强散射中心贡献,说明雷达回波在属性散射中心的参数空间具有很强的稀疏性。考虑属性散射中心参数空间维数较高,导致构造的冗余字典维数远远大于我们可处理的维数,本发明的方法基于稀疏分解理论,并利用坐标轮回下降技术降低问题维数并构建超分辨字典,通过求解0范数优化问题提取属性散射中心特征;建立散射中心强度门限并在雷达图像中进行强散射中心检测,结合先验信息确定散射中心参数的取值集合;根据属性散射中心模型,利用坐标轮回下降法技术构建超分辨字典,得到属性散射中心参数估计集合与目标及其重要部件几何尺寸特征。
二.技术方案
本发明的实现步骤包括如下:
1)根据雷达图像中噪声设定散射中心强度门限ξ,将雷达图像中强度大于ξ的散射中心确定为强散射中心,根据检测到的强散射中心坐标(x,y)确定坐标参数x,y的取值范围,由先验信息确定散射中心长度L取值范围,设定分布式散射中心方位角最终确定散射中心参数的集合
2b)根据属性散射中心模型,产生对应不同参数的原子di(f,φ)
i=1,…N0
其中i表示原子序号,N0表示初步估计(x,y,L)时字典的原子个数,exp(·)为自然指数函数,sinc(·)为辛克函数,为集合离散化后的第i组参数,f为雷达发射信号频率,φ为雷达波束方位角,c为光速;
将不同原子di(f,φ)列向量化,并对原子向量进行能量归一化,构建字典为:
其中vec(·)表示列向量化操作,||·||2为2范数算子;
其中,σ为待优化散射系数向量,σ1为优化的散射系数向量,s为频域观测信号矩阵的列向量化,||·||0为0范数算子,ε为能量误差约束因子,根据目标支撑区占整幅雷达图像的能量比确定;
4a)根据先验信息确定散射中心长度L取值范围,将相邻散射中心长度参数L的间隔设为ρa,得到字典的散射中心长度参数L的离散取值;根据雷达回波数据录取的方位角域确定分布式散射中心方位角取值范围,将相邻分布式散射中心方位角参数的间隔设为得到字典的分布式散射中心方位角参数的离散取值;散射中心坐标参数
其中σ2为优化的散射系数向量;
6a)在以中的坐标参数x为中心,长度为ρr的邻域内取值,将相邻参数x的间隔设为ρr/Ns,得到字典的坐标参数x的离散取值,Ns为超分辨倍数,一般Ns=2,4,8,...;在以中的y为中心,长度为ρa的邻域内取值,将相邻参数y的间隔设为ρa/Ns,得到字典的坐标参数y的离散取值;在以中的L为中心,长度为ρa的邻域内取值,将相邻参数L的间隔设为ρa/Ns,得到字典的散射中心长度参数L的离散取值;分布式散射中心方位角参数
其中σ3为优化的散射系数向量。
本发明与现有技术相比具有以下优点:
(1)现有提取属性散射中心特征的方法存在模型失配问题,并且由于图像分割会导致目标的一些重要散射中心易于丢失。本发明提出的基于稀疏分解的雷达目标属性散射中心特征提取方法从频域观测数据出发,通过求解0范数优化问题得到频域观测数据的稀疏表示,可以有效提取目标散射中心。
(2)现有方法基于对雷达图像分割的结果,进行参数初始化和模型阶数选择,通过近似最大似然方法估计雷达图像子区域中的属性散射中心的参数,这类方法估计的参数精度较低。本发明方法通过利用坐标轮回下降技术构建超分辨字典,对频域观测数据进行稀疏分解,得到雷达目标属性散射中心参数的超分辨估计。
(3)由于本发明方法能够有效提取目标散射中心,并得到散射中心参数的超分辨估计,所以本发明根据提取的散射中心参数集合,可以得到目标及其重要部件的精确几何尺寸特征。
附图说明
图1是本发明的流程图;
图2是属性散射中心模型中频率依赖因子α=1与频率依赖因子α=0时的距离频率域数据变化波形;
图3是属性散射中心模型中方位依赖因子γ=10-11与方位依赖因子γ=0时的方位频率域数据变化波形;
图4是雷达目标T72坦克的原始雷达图像;
图5是通过本发明方法对图4重构的雷达目标T72坦克的雷达图像;
图6是原始雷达目标D7铲车雷达图像;
图7是通过本发明方法对图6重构的雷达目标D7铲车雷达图像。
具体实施方式
一、技术原理
传统雷达图像是以点散射模型为基础得到的,该模型只包含目标散射点位置信息,但仅利用目标散射点的位置信息构建的识别特征并不能完备表征雷达图像中目标的本质属性。在光学区,扩展目标的高频电磁散射响应可以用一组独立分布的散射体,或称散射中心的电磁散射响应之和近似表示。属性散射中心模型用一组参数描述每个散射中心的位置、形状、方向以及幅度等,这些属性都是关系目标的重要信息;同点散射模型相比,属性散射中心模型包含了更丰富的可用于目标分类识别的特征。
根据属性散射中心模型可知,目标中第i个散射中心频率-方位二维回波信号为:
其中,i表示散射中心序号,f为雷达发射信号频率,φ为雷达方位角,exp(·)为自然指数函数,sin c(·)为辛克函数,c为光速,θi表示第i个散射中心的参数向量Ai为散射中心的散射强度,xi为距离维坐标,yi为方位维坐标,Li为分布式散射中心方位维的长度,为分布式散射中心的方位角,αi为频率依赖因子,一般αi∈{-1,-0.5,0,0.5,1},γi为局部式散射中心的方位依赖因子。
由各个散射中心的回波信号之和,即可构成目标频率-方位二维回波信号:
其中i表示散射中心序号,Ei(f,φ;θi)为第i个散射中心的回波信号,M为散射中心个数,θ表示M个散射中心参数矩阵,(·)T表示转置操作;
将式<9>的目标回波信号用矩阵形式表达,其表达式为:
s=D(θ)σ+n <10>
其中s为目标回波信号E(f,φ;θ)的列向量化,D(θ)为散射中心参数矩阵θ对应的字典,σ为散射系数向量,n表示高斯白噪声。
在雷达回波中,由于目标散射场的绝大部分能量仅由少量强散射中心贡献,故说明雷达回波在属性散射中心的参数空间具有很强的稀疏性。根据雷达回波在属性散射中心的参数空间的稀疏性,通过求解0范数优化问题即可得到观测数据s的稀疏表示,以及散射中心参数估计集合,通过散射中心参数集合可以估计目标及其重要部件的几何尺寸特征。
由于属性散射中心参数空间维数较高,导致构造的冗余字典维数远远大于需要处理的维数,因此可通过分析属性散射中心模型中的频率依赖因子α、方位依赖因子γ对回波信号的影响,对属性散射中心模型进行简化,降低参数空间维数。
当雷达的中心频率fc=9.6GHz,带宽B=591MHz,其它散射中心参数为0时,仅仅改变频率依赖因子α或者方位依赖因子γ,而频率依赖因子α与方位依赖因子γ对频域回波数据有如下影响:
a)由于雷达相对带宽较小,频率依赖因子α使散射中心的散射强度在距离频率维可发生微小变化,如图2所示,其中图2(a)中频率依赖因子α=1,图2(b)中频率依赖因子α=0。
b)由于雷达方位角φ变化范围较小,局部式散射中心方位依赖因子γ的变化能使局部式散射中心的散射强度在方位频率上稍有变化,如图3所示,其中图3(a)中方位依赖因子γ=10-11,图3(b)中方位依赖因子γ=0。
通过以上分析可知,频率依赖因子α与方位依赖因子γ对目标频域回波数据影响很小,因此估计参数时可以忽略频率依赖因子α与方位依赖因子γ,即认为频率依赖因子α=0,方位依赖因子γ=0,则原始属性散射中心模型式<8>可简化为:
根据简化后的属性散射中心模型,基于稀疏分解,利用坐标轮回下降技术,可以得到目标散射中心参数的超分辨估计,通过散射中心参数集合可以估计目标及其重要部件的几何尺寸特征。
二、实现步骤
参照图1,本发明的实现步骤如下:
步骤1,建立散射中心强度门限ξ进行强散射中心检测,确定散射中心参数取值范围。
1a)雷达成像:
由式<11>可知,雷达回波信号在频率-方位极坐标空间下均匀采样,雷达成像将回波信号重采样或者插值,得到在笛卡尔坐标系(fx,fy)下均匀采样的回波信号E(fx,fy),其中fx为距离频率,fy为方位频率;
fx=fcos(φ),fy=fsin(φ) <12>
其中,φ为雷达方位角;
对回波信号E(fx,fy)做二维逆傅里叶变换,得到目标的雷达图像;
1b)强散射中心检测,确定散射中心参数取值范围:
根据雷达图像中噪声设定散射中心强度门限ξ,将雷达图像中强度大于ξ的散射中心确定为强散射中心,根据检测到的强散射中心坐标(x,y)确定坐标参数x,y的取值范围,由先验信息确定散射中心长度L取值范围,设定分布式散射中心方位角最终确定散射中心参数的集合
例如,当雷达图像数据为公开MSTAR实测数据中方位角为80.774185°俯仰角为15°的T72坦克数据时,根据雷达图像中噪声设定散射中心强度门限ξ=0.28,确定的散射中心距离坐标参数x的取值范围为[11.1675,14.467]米,方位坐标参数y的取值范围为[13.40625,4.875]米,散射中心长度L取值范围为[0,4.0625]米,分布式散射中心方位角
例如,当雷达图像数据为公开MSTAR实测数据中方位角为80.774185°俯仰角为15°的T72坦克数据时,根据雷达图像中噪声设定散射中心强度门限ξ=0.28,确定的散射中心距离坐标参数x的取值范围为[11.1675,14.467]米,将相邻坐标参数x间隔为0.2538米,方位坐标参数y的取值范围为[13.40625,4.875]米,将相邻坐标参数y的间隔设为0.203125米,散射中心长度L取值范围为[0,4.0625]米,将相邻散射中心长度参数L的间隔设为0.203125米,分布式散射中心方位角
2b)根据属性散射中心模型,产生对应不同参数的原子di(f,φ):
i=1,…N0
其中,f为雷达发射信号频率,φ为雷达方位角,exp(·)为自然指数函数,sin c(·)为辛克函数,c为光速,i表示原子序号,N0表示初步估计(x,y,L)时字典的原子个数,为集合离散化后的第i组参数。
i=1,…N0
步骤3,利用正交匹配追踪OMP求解如式<15>的0范数优化问题:
其中σ为待优化散射系数向量,σ1为优化的散射系数向量,s为频域观测信号矩阵的列向量化,||·||0为0范数算子,ε为能量误差约束因子,其根据目标支撑区占整幅雷达图像的能量比确定。
本步骤的求解步骤如下:
3b)初始化,即将信号余量r初始化为频域观测信号矩阵的列向量s,将散射系数向量σ初始化为0,将原子索引向量a初始化为空,设初始迭代次数k=1,并开始迭代;
3c)计算字典与信号余量r的相关系数向量C
其中(·)H表示共轭转置,相关系数向量C中最大元素的序号为ak,可得原子索引向量a为:
a=[a1,…,ak] <17>
3d)利用最小二乘方法计算散射系数向量σ
更新余量r:
r=s-D0(:,a)·σ(a) <19>
其中σ(a)为散射系数向量σ中序号为ai的k个元素构成的列向量i=1,…,k;
3e)计算重构能量比ηk;
sk=D0(:,a)·σ(a)
其中sk为重构信号,表示2范数的平方;
3f)如果ηk-ηk-1≤δ停止迭代,σ1=σ,其中δ为相邻迭代能量比之差的门限,取δ=0.001;否则k=k+1并转至步骤3c)。
5a)根据先验信息确定散射中心长度L取值范围,将相邻散射中心长度参数L的间隔设为ρa,得到字典的散射中心长度参数L的离散取值;根据雷达回波数据录取的方位角域确定分布式散射中心方位角取值范围,将相邻分布式散射中心方位角参数的间隔设为得到字典的分布式散射中心方位角参数的离散取值;散射中心坐标参数
例如,当雷达图像数据为公开MSTAR实测数据中方位角为80.774185°俯仰角为15°的T72坦克数据时,散射中心长度L取值范围为[0,4.0625]米,将相邻散射中心长度参数L的间隔设为0.203125米,分布式散射中心方位角的取值范围为[-3°,3°],将相邻分布式散射中心方位角参数的间隔设为0.25°。
步骤6,利用正交匹配追踪OMP求解如式<22>的0范数优化问题:
其中σ2为优化的散射系数向量,其求解步骤如下:
6a)输入字典
6b)初始化,即将信号余量r初始化为频域观测信号矩阵的列向量s,将散射系数向量σ初始化为0,将原子索引向量a初始化为空,设初始迭代次数k=1,并开始迭代;
6c)计算字典与信号余量r的相关系数向量C
其中(·)H表示共轭转置,相关系数向量C中最大元素的序号为ak,可得原子索引向量a为:
a=[a1,…,ak] <24>
6d)利用最小二乘方法计算散射系数向量σ
更新余量r:
r=s-D1(:,a)·σ(a) <26>
其中σ(a)为散射系数向量σ中序号为ai的k个元素构成的列向量i=1,…,k;
6e)计算重构能量比ηk;
sk=D1(:,a)·σ(a)
6f)如果ηk-ηk-1≤δ停止迭代,σ2=σ,其中δ为相邻迭代能量比之差的门限,取δ=0.001;否则k=k+1并转至步骤6c)。
8a)确定字典的距离维坐标参数x离散取值:在以参数集合中的距离维坐标参数x为中心,长度为ρr的邻域内取值,将相邻参数x的间隔设为ρr/Ns,得到字典的距离维坐标参数x的离散取值,其中Ns为超分辨倍数,一般Ns=2,4,8,...;
8e)根据上述步骤确定的距离维坐标参数x,方位维坐标参数y,散射中心长度参数L,分布式散射中心方位角参数离散取值,由式<13>产生与离散化的字典参数组合对应的原子向量,并对原子进行能量归一化,构建超分辨字典为:
步骤9,利用正交匹配追踪OMP方法求解如式<29>的0范数优化问题:
其中σ3为优化的散射系数向量,其求解步骤如下:
9b)初始化,即将信号余量r初始化为频域观测信号矩阵的列向量s,将散射系数向量σ初始化为0,将原子索引向量a初始化为空,设初始迭代次数k=1,并开始迭代;
9c)计算字典与信号余量r的相关系数向量C
其中(·)H表示共轭转置,相关系数向量C中最大元素的序号为ak,得到原子索引向量a为:
a=[a1,…,ak] <31>
9d)利用最小二乘方法计算散射系数向量σ
更新余量r:
r=s-D2(:,a)·σ(a) <33>
其中σ(a)为散射系数向量σ中序号为ai的k个元素构成的列向量i=1,…,k;
9e)计算重构能量比ηk:
sk=D2(:,a)·σ(a)
9f)如果ηk-ηk-1≤δ停止迭代,σ3=σ其中δ为相邻迭代能量比之差的门限,取δ=0.001;否则k=k+1并转至步骤9c)。
例如,当雷达图像数据为公开MSTAR实测数据中方位角为80.774185°俯仰角为15°的T72坦克数据时,根据最终提取的散射中心及其参数集合估计T72坦克的几何尺寸为:T72坦克炮筒长度为6.4米,坦克体长度为6.8米,T72坦克的总长度为9.4米,T72坦克的宽度估值为3.8米。
本发明的效果通过以下对实测数据的实验进一步说明:
1)实验场景:
实验所用的实测数据为公开的MSTAR数据集中方位角为80.774185°俯仰角为15°的T72坦克的合成孔径雷达SAR数据以及方位角为0.307442°俯仰角为17°的D7铲车的SAR数据,雷达的中心频率fc=9.6GHz,带宽B=591MHz。T72坦克真实几何尺寸为:车体长6.41m,车宽3.52m,炮筒长6.155m,炮向前时炮筒伸出长度3.035m。D7铲车真实几何尺寸为:车宽2.56m,车长4.2m,车高3.294m,履带板宽度0.56m,履带中心距1.981m,接地履带长度2.72m。2)实验内容:
2a)对于方位角为80.774185°俯仰角为15°的T72坦克MSTAR数据,利用坐标轮回下降法构建超分辨字典,通过求解0范数优化问题得到散射中心参数估计集合。图4为原始T72坦克SAR图像,根据提取的散射中心及其参数重构的T72坦克图像如图5所示,由图4与图5可知本发明能够准确重构目标图像。
2b)对于方位角为0.307442°俯仰角为17°的D7铲车MSTAR数据,利用坐标轮回下降法构建超分辨字典,通过求解0范数优化问题得到散射中心参数估计集合。图6为原始D7铲车SAR图像,根据提取的散射中心及其参数重构的T72坦克图像如图7所示,由图6与图7可知本发明能够准确重构目标图像。
3)实验结果分析:
表1本发明方法处理T72坦克数据收敛过程
步骤3 | 步骤6 | 步骤9 | |
重构能量比 | 79.8% | 83.7% | 86.5% |
字典原子个数 | 9030 | 63000 | 47000 |
散射中心个数 | 126 | 94 | 92 |
表2本发明方法处理D7铲车数据收敛过程
步骤3 | 步骤6 | 步骤9 | |
重构能量比 | 84% | 86.7% | 88.7% |
字典原子个数 | 4800 | 31125 | 29925 |
散射中心个数 | 83 | 57 | 57 |
表1以及表2分别给出了本发明的方法处理T72坦克和D7铲车数据过程中重构能量比、字典原子个数以及散射中心个数的收敛情况。由于步骤3假定而仅仅对参数(x,y,L)进行初步估计,所以步骤3的重构能量较低;当步骤6在步骤3的基础上考虑分布式散射中心方位角时,重构能量比得到提高,提取的散射中心个数降低;步骤9中字典的参数分辨率增加,重构能量比得到提高,散射中心个数也进一步减少,字典原子个数得到降低,这是由于步骤3至步骤9之间利用了坐标轮回下降技术。表1、表2中数据收敛过程表明本发明方法能够有效提取散射中心。
根据最终提取的散射中心及其参数可以估计T72坦克与D7铲车的几何尺寸。
表3用于估计T72坦克炮筒长度的散射中心参数
表3表明,根据散射中心的坐标和散射中心长度估计坦克炮筒长度为6.4米。
表4用于估计T72坦克体长度的散射中心参数
表4表明,根据散射中心的坐标参数和散射中心长度估计坦克体长度为6.8米。
由目标支撑区可以得到T72坦克的总长度估值为9.4米;根据T72坦克上边缘的散射中心至炮筒散射中心的距离为T72坦克宽度的一半,得到T72坦克的宽度估值为3.8米。
表5用于估计D7铲车铲刀长度的散射中心参数
表5表明,根据散射中心的坐标参数和散射中心长度估计D7铲车铲刀长为3米。
表6用于估计D7铲车宽度的散射中心参数
表6表明,根据散射中心的坐标参数和散射中心长度估计D7铲车宽度为2.6米。
表7用于估计D7铲车履带宽度的散射中心参数
表7表明,根据散射中心的坐标参数和散射中心长度估计D7铲车履带宽度为0.6米。
由提取的履带散射中心可以估计D7铲车的接地履带长度为2.6米;由目标支撑区得到D7铲车车长估值为4.31米。
结合T72坦克与D7铲车的真实几何尺寸,可知本发明方法估计的目标及其重要部件的几何尺寸误差在6%以内,实验结果说明本发明方法能够有效地准确提取目标属性散射中心特征和精确估计目标及其重要部件的几何尺寸。
Claims (4)
1.一种基于稀疏分解的雷达目标属性散射中心特征提取方法,包括如下步骤:
1)根据雷达图像中噪声设定散射中心强度门限ξ,将雷达图像中强度大于ξ的散射中心确定为强散射中心,根据检测到的强散射中心坐标(x,y)确定坐标参数x,y的取值范围,由先验信息确定散射中心长度L取值范围,设定分布式散射中心方位角最终确定散射中心参数的集合
2b)根据属性散射中心模型,产生对应不同参数的原子di(f,φ)
i=1,...N0
其中i表示原子序号,N0表示初步估计(x,y,L)时字典的原子个数,exp(·)为自然指数函数,sinc(·)为辛克函数,为集合离散化后的第i组参数,f为雷达发射信号频率,φ为雷达波束方位角,c为光速;
其中vec(·)表示列向量化操作,||·||2为2范数算子;
其中,σ为待优化散射系数向量,σ1为优化的散射系数向量,s为频域观测信号矩阵的列向量化,||·||0为0范数算子,ε为能量误差约束因子,根据目标支撑区占整幅雷达图像的能量比确定;
4a)根据先验信息确定散射中心长度L取值范围,将相邻散射中心长度参数L的间隔设为ρa,得到字典的散射中心长度参数L的离散取值;根据雷达回波数据录取的方位角域确定分布式散射中心方位角取值范围,将相邻分布式散射中心方位角参数的间隔设为得到字典的分布式散射中心方位角参数的离散取值;散射中心坐标参数
其中σ2为优化的散射系数向量;
6a)在以中的坐标参数x为中心,长度为ρr的邻域内取值,将相邻参数x的间隔设为ρr/Ns,得到字典的坐标参数x的离散取值,Ns为超分辨倍数,一般Ns=2,4,8,...;在以中的y为中心,长度为ρa的邻域内取值,将相邻参数y的间隔设为ρa/Ns,得到字典的坐标参数y的离散取值;在以中的L为中心,长度为ρa的邻域内取值,将相邻参数L的间隔设为ρa/Ns,得到字典的散射中心长度参数L的离散取值;分布式散射中心方位角参数
其中,N2表示超分辨字典的原子个数;
其中σ3为优化的散射系数向量。
2.根据权利要求1所述的方法,其中步骤3所述的利用正交匹配追踪OMP方法求解如式<3>的0范数优化问题,按如下步骤进行:
3b)初始化,即将信号余量r初始化为频域观测信号矩阵的列向量s,将散射系数向量σ初始化为0,将原子索引向量a初始化为空,设初始迭代次数k=1,并开始迭代;
3c)计算字典与信号余量r的相关系数向量C
其中(·)H表示共轭转置,相关系数向量C中最大元素的序号为ak,可得原子索引向量a为:
a=[a1,...,ak] <9>
3d)利用最小二乘方法计算散射系数向量σ
其中表示伪逆,D0(:,a)为字典D0中列号为ai的k个列向量构成的矩阵i=1,...,k,k为当前迭代次数;
更新余量r:
r=s-D0(:,a)·σ(a) <11>
其中σ(a)为散射系数向量σ中序号为ai的k个元素构成的列向量i=1,...,k;
3e)计算重构能量比ηk;
sk=D0(:,a)·σ(a)
3f)如果ηk-ηk-1≤δ停止迭代,其中δ为相邻迭代能量比之差的门限,取δ=0.001;否则k=k+1并转至步骤3c)。
3.根据权利要求1所述的方法,其中步骤5所述的利用正交匹配追踪OMP方法求解如式<5>的0范数优化问题,按如下步骤进行:
5b)初始化,即将信号余量r初始化为频域观测信号矩阵的列向量s,将散射系数向量σ初始化为0,将原子索引向量a初始化为空,设初始迭代次数k=1,并开始迭代;
5c)计算字典与信号余量r的相关系数向量C
其中(·)H表示共轭转置,相关系数向量C中最大元素的序号为ak,可得原子索引向量a为:
a=[a1,...,ak] <14>
5d)利用最小二乘方法计算散射系数向量σ
更新余量r:
r=s-D1(:,a)·σ(a) <16>
其中σ(a)为散射系数向量σ中序号为ai的k个元素构成的列向量i=1,...,k;
5e)计算重构能量比ηk;
sk=D1(:,a)·σ(a)
5f)如果ηk-ηk-1≤δ停止迭代,其中δ为相邻迭代能量比之差的门限,取δ=0.001;否则k=k+1并转至步骤5c)。
4.根据权利要求1所述的方法,其中步骤7所述的利用正交匹配追踪OMP方法求解如式<7>的0范数优化问题,按如下步骤进行:
7b)初始化,即将信号余量r初始化为频域观测信号矩阵的列向量s,将散射系数向量σ初始化为0,将原子索引向量a初始化为空,设初始迭代次数k=1,并开始迭代;
7c)计算字典与信号余量r的相关系数向量C
其中(·)H表示共轭转置,相关系数向量C中最大元素的序号为ak,得到原子索引向量a为:
a=[a1,...,ak] <19>
7d)利用最小二乘方法计算散射系数向量σ
更新余量r:
r=s-D2(:,a)·σ(a) <21>
其中σ(a)为散射系数向量σ中序号为ai的k个元素构成的列向量i=1,...,k;
7e)计算重构能量比ηk:
sk=D2(:,a)·σ(a)
7f)如果ηk-ηk-1≤δ停止迭代,其中δ为相邻迭代能量比之差的门限,取δ=0.001;否则k=k+1并转至步骤7c)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210414929.XA CN103064071B (zh) | 2012-10-25 | 2012-10-25 | 基于稀疏分解的雷达目标属性散射中心特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210414929.XA CN103064071B (zh) | 2012-10-25 | 2012-10-25 | 基于稀疏分解的雷达目标属性散射中心特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103064071A true CN103064071A (zh) | 2013-04-24 |
CN103064071B CN103064071B (zh) | 2014-07-23 |
Family
ID=48106775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210414929.XA Active CN103064071B (zh) | 2012-10-25 | 2012-10-25 | 基于稀疏分解的雷达目标属性散射中心特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103064071B (zh) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558593A (zh) * | 2013-10-16 | 2014-02-05 | 北京航空航天大学 | 一种基于强散射点的舰船特定点的位置获取方法 |
CN103605121A (zh) * | 2013-11-18 | 2014-02-26 | 南京理工大学 | 基于快速稀疏贝叶斯学习算法的宽带雷达数据融合方法 |
CN103630886A (zh) * | 2013-12-05 | 2014-03-12 | 西安电子科技大学 | 基于属性散射中心模型的isar目标分解与重构方法 |
CN103760544A (zh) * | 2014-01-14 | 2014-04-30 | 北京环境特性研究所 | 雷达目标的散射中心提取方法及*** |
CN103969634A (zh) * | 2014-04-29 | 2014-08-06 | 西安电子科技大学 | 基于全极化属性散射中心模型的目标属性特征提取方法 |
CN104122554A (zh) * | 2014-07-31 | 2014-10-29 | 西安电子科技大学 | 高分辨率sar图像目标的属性散射中心特征提取方法 |
CN104182753A (zh) * | 2014-07-31 | 2014-12-03 | 西安电子科技大学 | 结合图像分割与子空间匹配追踪的目标散射中心提取方法 |
CN104732224A (zh) * | 2015-04-08 | 2015-06-24 | 重庆大学 | 基于二维泽尔尼克矩特征稀疏表示的sar目标识别方法 |
CN105044722A (zh) * | 2015-08-03 | 2015-11-11 | 西安电子科技大学 | 合成孔径雷达目标的全贝叶斯特征提取方法 |
CN105068062A (zh) * | 2015-08-19 | 2015-11-18 | 西安电子科技大学 | 基于稀疏散射中心提取的距离像数据外推方法 |
CN106154266A (zh) * | 2016-07-19 | 2016-11-23 | 西安电子科技大学 | 基于属性散射中心模型的多视角目标部件特征提取方法 |
CN106646395A (zh) * | 2016-09-30 | 2017-05-10 | 西安电子科技大学 | 一种飞行目标的雷达回波推演方法 |
CN107491831A (zh) * | 2017-07-09 | 2017-12-19 | 电子科技大学 | 一种多监视区域下自适应终止的mimo雷达优化布站方法 |
CN107784664A (zh) * | 2017-12-05 | 2018-03-09 | 韶关学院 | 一种基于k稀疏的快速鲁棒目标跟踪方法 |
CN108196256A (zh) * | 2013-10-25 | 2018-06-22 | 德克萨斯仪器股份有限公司 | 雷达中的角分辨率 |
CN108550145A (zh) * | 2018-04-11 | 2018-09-18 | 北京环境特性研究所 | 一种sar图像质量评估方法和装置 |
CN108562897A (zh) * | 2018-01-26 | 2018-09-21 | 桂林电子科技大学 | 一种mimo穿墙雷达的结构稀疏成像方法和装置 |
CN109061643A (zh) * | 2018-08-13 | 2018-12-21 | 南京理工大学 | 基于三维属性散射中心模型的多雷达融合高分辨成像方法 |
DE102018001499B3 (de) | 2018-02-26 | 2019-01-31 | Mbda Deutschland Gmbh | Radargerät mit Gruppenantenne und Verfahren zur Quelllokalisation durch eine zweidimensionale Radargruppenantenne |
CN109444844A (zh) * | 2019-01-04 | 2019-03-08 | 北京环境特性研究所 | 一种提取目标散射中心特征的方法和装置 |
CN109766899A (zh) * | 2018-12-28 | 2019-05-17 | 西安电子科技大学 | 物理特征提取和svm的sar图像车辆目标识别方法 |
CN111859266A (zh) * | 2020-07-30 | 2020-10-30 | 北京环境特性研究所 | 基于多源信息融合的空间目标结构反演方法和装置 |
CN113050050A (zh) * | 2021-03-03 | 2021-06-29 | 中国人民解放军63861部队 | 具有电磁散射特性的高分辨雷达目标复现方法及装置 |
CN114137521A (zh) * | 2021-11-26 | 2022-03-04 | 中山大学·深圳 | 基于属性散射中心模型的图像识别方法、装置及介质 |
CN115877340A (zh) * | 2022-12-22 | 2023-03-31 | 中国人民解放军军事科学院***工程研究院 | 一种用于参数估计的数据处理方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001194454A (ja) * | 2000-01-12 | 2001-07-19 | Mitsubishi Electric Corp | 合成開口レーダ装置および目標散乱点検出方法 |
JP4305148B2 (ja) * | 2003-12-02 | 2009-07-29 | 三菱電機株式会社 | Sar画像超解像装置 |
CN102222230A (zh) * | 2011-03-25 | 2011-10-19 | 中国人民解放军空军雷达学院 | 基于参数化词典的多输入多输出成像雷达目标识别方法 |
-
2012
- 2012-10-25 CN CN201210414929.XA patent/CN103064071B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001194454A (ja) * | 2000-01-12 | 2001-07-19 | Mitsubishi Electric Corp | 合成開口レーダ装置および目標散乱点検出方法 |
JP4305148B2 (ja) * | 2003-12-02 | 2009-07-29 | 三菱電機株式会社 | Sar画像超解像装置 |
CN102222230A (zh) * | 2011-03-25 | 2011-10-19 | 中国人民解放军空军雷达学院 | 基于参数化词典的多输入多输出成像雷达目标识别方法 |
Non-Patent Citations (3)
Title |
---|
KUSH R. VARSHNEY等: "Sparse Representation in Structured Dictionaries With Application to Synthetic Aperture Radar", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
汪雄良等: "基于紧致字典的基追踪方法在SAR图像超分辨中的应用", 《电子学报》 * |
王光新: "基于稀疏约束正则化模型的图像提高分辨率技术研究", 《万方学位论文数据库》 * |
Cited By (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558593B (zh) * | 2013-10-16 | 2015-08-26 | 北京航空航天大学 | 一种基于强散射点的舰船特定点的位置获取方法 |
CN103558593A (zh) * | 2013-10-16 | 2014-02-05 | 北京航空航天大学 | 一种基于强散射点的舰船特定点的位置获取方法 |
CN108196256B (zh) * | 2013-10-25 | 2021-10-29 | 德克萨斯仪器股份有限公司 | 雷达中的角分辨率 |
CN108196256A (zh) * | 2013-10-25 | 2018-06-22 | 德克萨斯仪器股份有限公司 | 雷达中的角分辨率 |
CN103605121A (zh) * | 2013-11-18 | 2014-02-26 | 南京理工大学 | 基于快速稀疏贝叶斯学习算法的宽带雷达数据融合方法 |
CN103630886A (zh) * | 2013-12-05 | 2014-03-12 | 西安电子科技大学 | 基于属性散射中心模型的isar目标分解与重构方法 |
CN103630886B (zh) * | 2013-12-05 | 2015-10-28 | 西安电子科技大学 | 基于属性散射中心模型的isar目标分解与重构方法 |
CN103760544B (zh) * | 2014-01-14 | 2016-05-11 | 北京环境特性研究所 | 雷达目标的散射中心提取方法及*** |
CN103760544A (zh) * | 2014-01-14 | 2014-04-30 | 北京环境特性研究所 | 雷达目标的散射中心提取方法及*** |
CN103969634A (zh) * | 2014-04-29 | 2014-08-06 | 西安电子科技大学 | 基于全极化属性散射中心模型的目标属性特征提取方法 |
CN103969634B (zh) * | 2014-04-29 | 2016-08-24 | 西安电子科技大学 | 基于全极化属性散射中心模型的目标属性特征提取方法 |
CN104182753A (zh) * | 2014-07-31 | 2014-12-03 | 西安电子科技大学 | 结合图像分割与子空间匹配追踪的目标散射中心提取方法 |
CN104122554A (zh) * | 2014-07-31 | 2014-10-29 | 西安电子科技大学 | 高分辨率sar图像目标的属性散射中心特征提取方法 |
CN104182753B (zh) * | 2014-07-31 | 2017-08-25 | 西安电子科技大学 | 结合图像分割与子空间匹配追踪的目标散射中心提取方法 |
CN104122554B (zh) * | 2014-07-31 | 2016-08-17 | 西安电子科技大学 | 高分辨率sar图像目标的属性散射中心特征提取方法 |
CN104732224A (zh) * | 2015-04-08 | 2015-06-24 | 重庆大学 | 基于二维泽尔尼克矩特征稀疏表示的sar目标识别方法 |
CN104732224B (zh) * | 2015-04-08 | 2017-11-03 | 重庆大学 | 基于二维泽尔尼克矩特征稀疏表示的sar目标识别方法 |
CN105044722A (zh) * | 2015-08-03 | 2015-11-11 | 西安电子科技大学 | 合成孔径雷达目标的全贝叶斯特征提取方法 |
CN105044722B (zh) * | 2015-08-03 | 2018-09-25 | 西安电子科技大学 | 合成孔径雷达目标的全贝叶斯特征提取方法 |
CN105068062A (zh) * | 2015-08-19 | 2015-11-18 | 西安电子科技大学 | 基于稀疏散射中心提取的距离像数据外推方法 |
CN106154266A (zh) * | 2016-07-19 | 2016-11-23 | 西安电子科技大学 | 基于属性散射中心模型的多视角目标部件特征提取方法 |
CN106154266B (zh) * | 2016-07-19 | 2019-02-19 | 西安电子科技大学 | 基于属性散射中心模型的多视角目标部件特征提取方法 |
CN106646395B (zh) * | 2016-09-30 | 2019-07-09 | 西安电子科技大学 | 一种飞行目标的雷达回波推演方法 |
CN106646395A (zh) * | 2016-09-30 | 2017-05-10 | 西安电子科技大学 | 一种飞行目标的雷达回波推演方法 |
CN107491831A (zh) * | 2017-07-09 | 2017-12-19 | 电子科技大学 | 一种多监视区域下自适应终止的mimo雷达优化布站方法 |
CN107784664A (zh) * | 2017-12-05 | 2018-03-09 | 韶关学院 | 一种基于k稀疏的快速鲁棒目标跟踪方法 |
CN107784664B (zh) * | 2017-12-05 | 2021-07-27 | 韶关学院 | 一种基于k稀疏的快速鲁棒目标跟踪方法 |
CN108562897A (zh) * | 2018-01-26 | 2018-09-21 | 桂林电子科技大学 | 一种mimo穿墙雷达的结构稀疏成像方法和装置 |
CN108562897B (zh) * | 2018-01-26 | 2022-01-11 | 桂林电子科技大学 | 一种mimo穿墙雷达的结构稀疏成像方法和装置 |
DE102018001499B3 (de) | 2018-02-26 | 2019-01-31 | Mbda Deutschland Gmbh | Radargerät mit Gruppenantenne und Verfahren zur Quelllokalisation durch eine zweidimensionale Radargruppenantenne |
CN108550145A (zh) * | 2018-04-11 | 2018-09-18 | 北京环境特性研究所 | 一种sar图像质量评估方法和装置 |
CN109061643B (zh) * | 2018-08-13 | 2022-05-27 | 南京理工大学 | 基于三维属性散射中心模型的多雷达融合高分辨成像方法 |
CN109061643A (zh) * | 2018-08-13 | 2018-12-21 | 南京理工大学 | 基于三维属性散射中心模型的多雷达融合高分辨成像方法 |
CN109766899A (zh) * | 2018-12-28 | 2019-05-17 | 西安电子科技大学 | 物理特征提取和svm的sar图像车辆目标识别方法 |
CN109766899B (zh) * | 2018-12-28 | 2021-05-04 | 西安电子科技大学 | 物理特征提取和svm的sar图像车辆目标识别方法 |
CN109444844A (zh) * | 2019-01-04 | 2019-03-08 | 北京环境特性研究所 | 一种提取目标散射中心特征的方法和装置 |
CN111859266A (zh) * | 2020-07-30 | 2020-10-30 | 北京环境特性研究所 | 基于多源信息融合的空间目标结构反演方法和装置 |
CN113050050A (zh) * | 2021-03-03 | 2021-06-29 | 中国人民解放军63861部队 | 具有电磁散射特性的高分辨雷达目标复现方法及装置 |
CN114137521A (zh) * | 2021-11-26 | 2022-03-04 | 中山大学·深圳 | 基于属性散射中心模型的图像识别方法、装置及介质 |
CN114137521B (zh) * | 2021-11-26 | 2024-07-19 | 中山大学·深圳 | 基于属性散射中心模型的图像识别方法、装置及介质 |
CN115877340A (zh) * | 2022-12-22 | 2023-03-31 | 中国人民解放军军事科学院***工程研究院 | 一种用于参数估计的数据处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103064071B (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103064071B (zh) | 基于稀疏分解的雷达目标属性散射中心特征提取方法 | |
CN107728142B (zh) | 基于二维卷积网络的雷达高分辨距离像目标识别方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN104122554B (zh) | 高分辨率sar图像目标的属性散射中心特征提取方法 | |
CN102914773B (zh) | 一种多航过圆周sar三维成像方法 | |
CN104950305A (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN103163523A (zh) | 基于压缩感知的低空风切变风速估计方法 | |
CN103886337A (zh) | 基于多稀疏描述的最近邻子空间sar目标识别方法 | |
CN107329137B (zh) | 基于box样条小波框架的探地雷达地下目标定位方法及*** | |
CN104182753B (zh) | 结合图像分割与子空间匹配追踪的目标散射中心提取方法 | |
CN105044722A (zh) | 合成孔径雷达目标的全贝叶斯特征提取方法 | |
CN104698431A (zh) | 模糊分量空间角的估计及多通道sar方位解模糊方法 | |
CN102628938A (zh) | 基于噪声先验的联合高斯模型雷达目标稳健识别方法 | |
CN103886606B (zh) | 基于联合广义伽玛分布参数的sar图像分割方法 | |
CN105068062A (zh) | 基于稀疏散射中心提取的距离像数据外推方法 | |
CN104200458A (zh) | 基于MeanShift的高分辨率遥感影像分割距离度量优化方法 | |
CN110596706B (zh) | 一种基于三维图像域投射变换的雷达散射截面积外推方法 | |
CN102254176B (zh) | 基于状态空间模型的雷达高分辨距离像目标识别方法 | |
CN103412306A (zh) | 基于多站多波段带宽融合的超宽带相参处理方法 | |
CN102928826A (zh) | 基于重构信号和时间采样的空中机动目标检测方法 | |
CN103631990A (zh) | Sar照射区域的仿真场景模型建立方法和*** | |
CN104021399A (zh) | 基于距离像时频图非负稀疏编码的sar目标识别方法 | |
CN109766899A (zh) | 物理特征提取和svm的sar图像车辆目标识别方法 | |
CN103777198B (zh) | 基于投影梯度的目标高度与反射面高度联合估计方法 | |
CN103605980A (zh) | 基于局部特征匹配的isar目标部件检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |