CN102737253B - 一种sar图像目标识别方法 - Google Patents
一种sar图像目标识别方法 Download PDFInfo
- Publication number
- CN102737253B CN102737253B CN201210201460.1A CN201210201460A CN102737253B CN 102737253 B CN102737253 B CN 102737253B CN 201210201460 A CN201210201460 A CN 201210201460A CN 102737253 B CN102737253 B CN 102737253B
- Authority
- CN
- China
- Prior art keywords
- coefficient
- vector
- target data
- alpha
- training sample
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开了一种SAR图像目标识别方法,本发明的方法利用稀疏表示理论将目标数据表示为训练样本的线性组合,通过求解最优化问题得到了具有可区分能力的近似非负稀疏系数,然后基于各类别系数和的大小确定样本的类别。本发明的方法训练样本系数值的大小反应了目标数据与该训练样本的相似性程度,系数值越大相似性越高,反之则相似性越低,因此实际上目标数据对应系数的大小在某种程度上可以体现目标数据的真实类别;而且通过非负约束保证测试图像均为各训练样本的非负加权和,使其具有更加可以解释的意义,同时也更有利于识别。
Description
技术领域
本发明属于合成孔径雷达(SAR,Synthetic Aperture Radar)自动目标识别技术领域,具体涉及一种SAR图像目标识别方法。
背景技术
SAR相对于普通的光学成像等***,对地面目标特别是地面静止目标具有强大的成像优势,作为一种重要的目标探测手段受到广泛的关注和研究。随着SAR成像体制方式的不断创新和发展,各种不同应用的SAR***开始丰富起来,且成像功能及效果也不断增强和完善,利用这些强大的成像***我们获得了大量的探测数据。这些数据在为我们提供对目标更加准确全面的探测结果的同时,也使人们面临着如何从大量探测数据中及时而准确地提取有用信息的困境,单纯依靠人工辨识已经无法跟上数据源膨胀的节奏。当前SAR图像解译***的发展已经远滞后于信息源的获取能力,SAR图像自动目标识别作为SAR图像解译和分析的重要组成部分,具有重要的民用和军事价值,日益成为国内外图像处理和模式识别领域的研究热点。
SAR图像与平台的俯仰角、成像的方位角、目标的电磁散射特性、回波的噪声分布等诸多因素有关,所以SAR图像本身具有复杂的图像特点,将稀疏表示应用于SAR图像目标识别时需要设计针对SAR图像特点的识别算法;另外目前采用的稀疏表示重建算法大多没有对求解的结果增加非负的约束,导致系数向量中存在负系数,例如基追踪(Basis Pursuit,BP)、匹配追踪(Matching Pursuit,MP)、正交匹配追踪(Orthogonal Matching Pursuit,OMP)、梯度投影(Gradient Projection for Sparse Reconstruction,GPSR)等算法得到的系数向量中都可能存在负系数,虽然负系数也体现了该训练样本与目标数据之间的某种对应关系,但是这些负系数在样本之间的稀疏表示过程中将可能会产生复杂的减性关系,而不具备目标识别中由局部累加为一个整体的直观性,且应用中难以对这些负系数赋予明确的意义使之便于理解,这些不利因素综合起来将最终降低SAR图像目标识别性能。
发明内容
本发明的目的是为了解决现有的SAR图像目标识别存在的上述问题,提出一种SAR图像目标识别方法。
本发明的技术方案是:一种SAR图像目标识别方法,具体包括如下步骤:
S1.将k类训练样本的特征向量组成矩阵Θ=[Θ1,Θ2,...,Θi,...,Θk],其中,第i类目标的ni个训练样本的特征向量组成矩阵的列向量集令n=k×ni,可得Θ∈Rm×n,矩阵的列向量为SAR图像训练样本的m维特征向量,任一目标数据的特征向量y∈Rm,其中表示m×ni维实数数组,Rm×n表示m×n维实数矩阵,Rm表示m维实数向量;
S2.1.将目标数据表示成训练样本的线性组合,得到目标数据和训练样本之间的稀疏表示模型,即:
y=Θα
其中,α∈Rn且α为线性表示的系数向量;
S2.2.将式y=Θα的求解转化为最小l1范数凸优化问题,并对系数向量α施加非负约束,则求解y=Θα的非负稀疏重建模型为:
s.t.||Θα-y||2≤ε
α≥0
由于限制α≥0,α直接表示为:
α=u,u≥0
所以 可表示为:
s.t.u≥0
有:
将 重写为标准的BCQP形式有:
s.t.z≥0
其中,
z=u,b=ΘTy,c=τ1n-b,B=ΘTΘ
目标函数F(z)的梯度方程为:
S4.输出目标数据类别identity(y):将测试样本归类到样本系数之和最大的类别:
本发明的有益效果:本发明的方法利用稀疏表示理论将目标数据表示为训练样本的线性组合,通过求解最优化问题得到了具有可区分能力的近似非负稀疏系数,然后基于各类别系数和的大小确定样本的类别,从以上介绍可以看出该方法有四个特点:其一,通过非负约束保证测试图像均为各训练样本的非负加权和,使其具有更加可以解释的意义,同时也更有利于识别;其二,训练样本系数值的大小反应了目标数据与该训练样本的相似性程度,系数值越大相似性越高,反之则相似性越低,因此实际上目标数据对应系数的大小在某种程度上可以体现目标数据的真实类别;其三,目标数据的类别不再依赖于单个或小部分训练样本,而是由所有训练样本系数的相对大小共同决定,避免了只由单个或小部分训练样本最终决定目标数据类别的偶然性;其四,当训练样本的种类或数量越多时,往往能够真正决定目标数据类别的训练样本数量相对却减少,也就是说与目标数据相关的训练样本在所有训练样本中所占的比例相对减少,因此稀疏表示越容易取得更稀疏的样本表示结果,也越有利于实现对样本的正确识别。这些特点为本发明设计的识别方法的可行性提供了保证,有利于本发明取得更好的识别率和鲁棒性。
附图说明
图1本发明的实施流程图。
图2总体系数分布情况。
图3单个目标数据的系数分布。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的说明。
利用稀疏表示的重建算法可以求解稀疏表示的系数向量α0,理想情况下,系数向量α0中只有与目标数据同类的第i类训练样本系数可能非0,其它类别训练样本对应的系数都为0,也就是说,目标数据基于全体训练样本的线性组合实际上只依赖于部分同类型的训练样本。
但是实际应用中受成像角度、噪声干扰、模型误差等因素影响,使得许多相对较小的非0系数分散到其它多个类上。但是总体上目标数据y所属类中仍然集中了多数相对较大的非0系数,也就是说对于可正确分类的情况,目标数据对应类别中的非0系数在所有非0系数中仍然占有主要的部分。由此可见稀疏表示的系数向量本身对目标数据类别具有可区分能力,这种系数向量的分布特点为识别方法的设计提供了指导。因此本发明利用该分布特点设计的分类算法为:首先计算每一类训练样本的系数之和,然后将目标数据归类到样本系数之和最大的类别。
本发明仅涉及SAR图像目标识别阶段的技术内容,但识别前进行简单的SAR预处理和特征提取,将更有利于提高目标识别的性能。因此在本发明中当提到训练样本和目标数据时,均指代其对应的特征向量。另外,训练样本中必须包含目标数据的相关类别。
下面结合附图,给出本发明的具体实施过程:
S1.根据目标数据选取训练样本,训练样本包含目标数据类型。将k类训练样本的特征向量组成矩阵Θ=[Θ1,Θ2,...,Θi,...,Θk],其中第i类目标的ni个训练样本的特征向量组成矩阵的列向量集令n=k×ni,可得Θ∈Rm×n,矩阵的列向量为SAR图像训练样本的m维特征向量,任一目标数据的特征向量y∈Rm,其中表示m×ni维实数数组,Rm×n表示m×n维实数矩阵,Rm表示m维实数向量;
S2.1.目标数据y在全体训练样本下的线性表示为:
即:y=Θα,其中,是理想情况下稀疏表示的系数向量,α中只有与目标数据同类的第i类训练样本系数可能非0,其它类系数都为0,也就是说,目标数据基于全体训练样本的线性组合实际上只依赖于部分同类的训练样本。实际应用中受成像角度、噪声干扰、模型误差等因素影响,使得许多相对较小的非0系数分散到其它多个类上。但是总体上目标数据y所属类中仍然集中了多数相对较大的非0系数,也就是说对于可正确分类的情况,目标数据对应类别的训练样本中的非0系数在所有非0系数中仍然占有主要的部分。
s.t.||Θα-y||2≤ε (1)
α≥0
本发明在GPSR算法的基础上推导出非负梯度投影稀疏重建算法(Non-negative GradientProjection for Sparse Reconstruction,GPSR)对该问题进行求解。通过NGPSR算法就可以获得非负系数向量
该方法对公式(1)的求解形式可以转化为:
式中α∈Rn,||z||2=(∑(zi)2)1/2表示z的欧几里德范数,||z||1=∑|zi|表示z的l1范数。是为了优化稀疏表示的误差,||α||1是为了优化系数向量的稀疏性,τ为预先设置的调节系数,具体为非负常数,其作用是为了调节两者在整体优化结果中的比重。l1范数形式将引导向量α中值较小的系数逐渐趋向于零,最终促进解向量α的稀疏性,这种性质使得公式在统计学和信号处理中稀疏性问题方面占有重要地位。
由于限制α≥0,α可以直接表示为:
α=u,u≥0 (3)
因此其中,为n个1组成的向量,由于α=u,为了实现正确的向量运算,这里的n与向量α∈Rn中的n取同样大小。所以式(2)又可以描述为如下边界受限的二次方程(Bound-Constrained Quadratic Program,BCQP):
s.t.u≥0
经推导,有:
则式(4)等价于:
s.t u≥0
式(6)也可重写为如下更标准的BCQP形式:
s.t.z≥0
其中,
z=u,b=ΘTy,c=τ1n-b,B=ΘTΘ
式(7)中目标函数的梯度方程为:
下面将利用梯度投影算法(Gradient Projection,GP)求解式(7),进而求得稀疏表示的非负系数向量的
GP算法按照下面描述的公式从z(t)迭代到z(t+1)。首先选择参数ρ(t)>0,并设:
其中,(·)+为取正算子,定义为(α)+=max{0,α}。然后选择第二个标量参数λ(t)∈[0,1],并设:
z(t+1)=z(t)+λ(t)(w(t)-z(t)) (10)
该算法的关键在于参数ρ(t)和λ(t)的计算,下面将详细描述。
首先定义向量g(t)为:
然后估算ρ(t)的初始值:
或者通过下列公式准确计算:
为了防止ρ0过大或过小,将ρ0限定在区间[ρmin,ρmax](0<ρmin<ρmax)内,并定义算子mid(a,b,c)为这三个标量值的中值。
综上所述,利用GP算法的求解步骤描述如下:
(2)利用式(13)计算ρ0,并用mid(ρmin,ρ0,ρmax)替代ρ0。
(3)在数字序列ρ0,βρ0,β2ρ0,…选择第一个满足下列不等式的常数等于ρ(t)
并令
(4)利用终止条件判断是否收敛:如果是,则结束循环,并得到稀疏解z(t+1);否则令t←t+1并返回步骤(2)。
S4.输出目标数据类别:将测试样本归类到样本系数之和最大的类别。
本发明的实验数据是MSTAR计划录取的实测SAR地面静止军用目标数据。包括3大类:BMP2(装甲车)、BMP70(装甲车)、T72(主战坦克),成像分辨率是0.3m×0.3m,方位角覆盖范围是0~360°,图像大小为128×128。实验使用该计划推荐的训练样本是目标在俯仰角为17°时的成像数据,目标数据是目标在俯仰角为15°时的成像数据。其中,目标数据集中的3大类包含6个型号,相同类型不同信号的目标在配置上有些差别,如表1所示为本实验中训练样本和目标数据的类别及其对应的图像数量。
表1
图2表明当目标数据的类别分别为BMP2、BTR70和T72时的总体系数分布情况,这里对目标数据的系数向量的绝对值之和做了归一化处理。由图可以看出,目标数据在相同类别训练样本上的系数分布数值更大,密度更集中,从总体上直观地反应出本发明可以实现对目标的有效识别。
图3是在三类目标中任意抽取的单个目标数据的系数分布。图3(a)是预处理后分割出的图像切片,图3(b)是图像切片基于训练样本的系数分布。可以看出单个目标数据时,系数向量的分布同样集中于与之类别相同的训练样本上,因此将目标数据归类为样本系数之和最大的类别将有利于实现对目标的正确识别。另外注意到系数分布图内的图像切片对应于系数值最大的训练样本,可以看出该训练样本与目标数据在目标形状和散射分布上都很相似,这说明系数向量中的系数值大小反应了目标数据与训练样本之间的相似程度。
表2给出了本发明基于表1实验数据的识别率情况,可以看出本发明对三类典型目标均取得了很高的识别率,同时与其它几种识别方法作对比可看出其平均识别率明显高于其他几种。
表2
本发明利用稀疏表示理论将目标数据表示为训练样本的线性组合,通过求解最优化问题得到了具有可区分能力的非负稀疏系数,分析了该非负稀疏系数在训练样本上的分布特性,从理论和实验两方面研究了该识别方法的原理和可行性。实验结果表明,本发明的可以实现对目标的有效识别,能明显提高对目标的正确识别率。
综上可以看出本发明在GPSR算法的基础上对目标数据的稀疏系数施加非负的约束,推导了一种非负梯度投影稀疏重建算法以保证测试图像均为各训练样本的非负加权和,使其具有更加可以解释的意义,同时也更有利于识别,最后利用重建的非负系数设计了识别算法实现对SAR图像的有效识别。
本领域的普通技术人员将会意识到,这里所述的基于MSTRA的实验是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实验。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (1)
1.一种SAR图像目标识别方法,具体包括如下步骤:
S1.将k类训练样本的特征向量组成矩阵Θ=[Θ1,Θ2,...,Θi,...,Θk],其中,第i类目标的ni个训练样本的特征向量组成矩阵的列向量集令n=k×ni,可得Θ∈Rm×n,矩阵的列向量为SAR图像训练样本的m维特征向量,任一目标数据的特征向量y∈Rm,其中表示m×ni维实数数组,Rm×n表示m×n维实数矩阵,Rm表示m维实数向量;
S2.1.将目标数据表示成训练样本的线性组合,得到目标数据和训练样本之间的稀疏表示模型,即:
y=Θα
其中,α∈Rn且α为线性表示的系数向量;
S2.2.将式y=Θα的求解转化为最小l1范数凸优化问题,并对系数向量α施加非负约束,则求解y=Θα的非负稀疏重建模型为:
s.t.||Θα-y||2≤ε
α≥0
由于限制α≥0,α直接表示为:
α=u,u≥0
所以 可表示为:
s.t.u≥0
有:
s.t.z≥0
其中,
z=u,b=ΘTy,c=τ1n-b,B=ΘTΘ
目标函数F(z)的梯度方程为:
S4.输出目标数据类别identity(y):将测试样本归类到样本系数之和最大的类别:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210201460.1A CN102737253B (zh) | 2012-06-19 | 2012-06-19 | 一种sar图像目标识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210201460.1A CN102737253B (zh) | 2012-06-19 | 2012-06-19 | 一种sar图像目标识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102737253A CN102737253A (zh) | 2012-10-17 |
CN102737253B true CN102737253B (zh) | 2014-03-05 |
Family
ID=46992717
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210201460.1A Expired - Fee Related CN102737253B (zh) | 2012-06-19 | 2012-06-19 | 一种sar图像目标识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102737253B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103426001B (zh) * | 2013-08-30 | 2016-08-10 | 西安电子科技大学 | Sar图像目标识别方法 |
CN104091335B (zh) * | 2014-07-04 | 2017-01-18 | 西安电子科技大学 | 极化sar图像舰船目标检测方法 |
CN104732224B (zh) * | 2015-04-08 | 2017-11-03 | 重庆大学 | 基于二维泽尔尼克矩特征稀疏表示的sar目标识别方法 |
CN105224621B (zh) * | 2015-09-22 | 2019-03-01 | 中国科学技术大学 | 一种压缩感知域的图像筛选、检索和识别方法 |
CN110472627B (zh) * | 2019-07-02 | 2022-11-08 | 五邑大学 | 一种端到端的sar图像识别方法、装置及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6272246B1 (en) * | 1997-02-05 | 2001-08-07 | Sharp Kabushiki Kaisha | Positioning device |
CN101551856A (zh) * | 2009-05-22 | 2009-10-07 | 西安电子科技大学 | 基于稀疏最小二乘支撑向量机的sar目标识别方法 |
CN102073875A (zh) * | 2011-01-06 | 2011-05-25 | 西安电子科技大学 | 基于稀疏表示的背景杂波量化方法 |
CN102122355A (zh) * | 2011-03-15 | 2011-07-13 | 西安电子科技大学 | 基于核稀疏表示的sar目标识别方法 |
-
2012
- 2012-06-19 CN CN201210201460.1A patent/CN102737253B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6272246B1 (en) * | 1997-02-05 | 2001-08-07 | Sharp Kabushiki Kaisha | Positioning device |
CN101551856A (zh) * | 2009-05-22 | 2009-10-07 | 西安电子科技大学 | 基于稀疏最小二乘支撑向量机的sar目标识别方法 |
CN102073875A (zh) * | 2011-01-06 | 2011-05-25 | 西安电子科技大学 | 基于稀疏表示的背景杂波量化方法 |
CN102122355A (zh) * | 2011-03-15 | 2011-07-13 | 西安电子科技大学 | 基于核稀疏表示的sar目标识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102737253A (zh) | 2012-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kim et al. | Color–texture segmentation using unsupervised graph cuts | |
CN102737253B (zh) | 一种sar图像目标识别方法 | |
CN109682382A (zh) | 基于自适应蒙特卡洛和特征匹配的全局融合定位方法 | |
Zhao et al. | Automatic recognition of loess landforms using Random Forest method | |
CN106096506A (zh) | 基于子类类间判别双字典的sar目标识别方法 | |
CN102129573A (zh) | 基于字典学习和稀疏表示的sar图像分割方法 | |
CN102592134A (zh) | 一种高光谱与红外数据多级决策融合分类方法 | |
CN109740588A (zh) | 基于弱监督和深度响应重分配的x光图片违禁品定位方法 | |
CN110533069B (zh) | 一种基于支持向量机算法的二维箔条分布特性识别方法 | |
CN103345757A (zh) | 多层次多特征约束下的光学和sar影像自动配准方法 | |
CN105354841A (zh) | 一种快速遥感影像匹配方法及*** | |
CN103440513A (zh) | 基于稀疏非负张量分解的大脑特定视觉认知状态判定方法 | |
CN106203483A (zh) | 一种基于语义相关多模态映射方法的零样本图像分类方法 | |
CN105447490A (zh) | 基于梯度回归树的车辆关键点检测方法及装置 | |
Pepik et al. | 3d object class detection in the wild | |
CN104021559A (zh) | 基于互信息和Harris角点检测的图像配准方法 | |
CN104657717A (zh) | 一种基于分层核稀疏表示的行人检测方法 | |
CN104732224B (zh) | 基于二维泽尔尼克矩特征稀疏表示的sar目标识别方法 | |
CN101196564A (zh) | 拉普拉斯正则化最小二乘合成孔径雷达自动目标识别方法 | |
CN104616022B (zh) | 一种近红外光谱的分类方法 | |
CN108629297A (zh) | 一种基于空域自然场景统计的遥感图像云检测方法 | |
CN103295031A (zh) | 一种基于正则风险最小化的图像目标计数方法 | |
CN103854290A (zh) | 一种结合骨架特征点和分布场描述子的扩展目标跟踪方法 | |
CN102314610A (zh) | 一种基于概率潜语义分析模型的面向对象影像聚类方法 | |
CN113723572B (zh) | 船只目标识别方法、计算机***及程序产品、存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140305 Termination date: 20150619 |
|
EXPY | Termination of patent right or utility model |