一种SAR雷达陆地坦克目标识别系统技术方案

本发明专利技术公开了一种SAR雷达陆地坦克目标识别系统，包括SAR雷达、数据库以及上位机；SAR雷达、数据库、上位机依次相连，所述SAR雷达对陆地进行实时监测，并将SAR雷达获得的图像数据存储到所述的数据库中，所述的上位机包括图像预处理模块、特征提取模块、特征选择模块、分类器训练模块以及结果显示模块。本发明专利技术提供一种实现在线识别、精度高的陆地坦克目标识别系统。

A SAR Radar Land Tank Target Recognition System

The invention discloses a SAR radar terrestrial tank target recognition system, which includes SAR radar, database and upper computer. SAR radar, database and upper computer are connected in turn. The SAR radar monitors land in real time and stores the image data obtained by SAR radar into the database. The upper computer includes: Image preprocessing module, feature extraction module, feature selection module, classifier training module and result display module. The invention provides a land tank target recognition system with online recognition and high accuracy.

【技术实现步骤摘要】
一种SAR雷达陆地坦克目标识别系统
本专利技术涉及图像处理领域，特别地，涉及一种SAR雷达陆地坦克目标识别系统。
技术介绍
数字图像与模拟图像相比具有存储，传输方便可压缩、传输过程中不易失真、处理方便等巨大优势，这些都为图像识别技术的发展提供了强大的动力。物体的识别主要指的是对三维世界的客体及环境的感知和认识，属于高级的计算机视觉范畴。它是以数字图像处理与识别为基础的结合人工智能、系统学等学科的研究方向，其研究成果被广泛应用在各种工业及探测机器人上。现代图像识别技术的一个不足就是自适应性能差，一旦目标图像被较强的噪声污染或是目标图像有较大残缺往往就得不出理想的结果。其中基于SAR图像的坦克目标检测已经有了广泛的研究，而坦克目标分类识别由于SAR图像分辨率的限制才刚刚起步，已有的一些研究成果也由于研究不够透彻，模型的效果并不是很好。因此，积极开展基于高分辨率SAR图像的坦克目标分类识别研究具有极其重大的意义。
技术实现思路
为了克服目前基于SAR图像的陆地坦克目标识别准确率不高的不足，本专利技术的目的在于提供一种实现实时分析的SAR雷达陆地坦克目标识别系统。本专利技术解决其技术问题所采用的技术方案是：一种SAR雷达陆地坦克目标识别系统，包括SAR雷达、数据库以及上位机，SAR雷达、数据库和上位机依次相连，所述SAR雷达对陆地进行实时监测，并将SAR雷达获得的图像数据存储到所述的数据库中，所述的上位机包括：图像预处理模块，用以进行SAR雷达图像数据预处理，采用如下过程完成：1)从数据库中传来的SAR图像灰度级为L,f(x0,y0)为像素点(x0,y0)处的灰度值，g(x0,y0)为像素点(x0,y0)的N×N邻域内像素的平均值，其中x0,y0分别表示像素点的横坐标和纵坐标；2)通过计算满足f＝m和g＝n的像素数目h(m,n)，得到二维联合概率密度pmn：pmn＝p(m,n)＝h(m,n)/M其中，M表示图像像素的总数目；3)计算二维直方图的均值向量μ：4)分别计算图像中目标和背景出现的概率P0,1和均值向量μ0,1：其中，t、s、下标0、下标1分别表示f分割阈值、g分割阈值、目标区域、背景区域；5)计算类间方差BCV：BCV＝P0(μ0-μ)(μ0-μ)′+P1(μ1-μ)(μ1-μ)′；其中，μ表示均值向量，上标’表示矩阵的转置。6)最佳阈值即为使得BCV为最大值时的二维阈值向量[s0,t0]：特征提取模块，用以进行坦克典型特征的提取，采用如下过程完成：1)从图像预处理模块传来的只包含一个坦克目标的SAR图像切片I(m,n)，其中只包含目标区域的二值图为B(m,n)，则只包含目标的图像T(m,n)：T(m,n)＝I(m,n)×B(m,n)其中，×表示对应像素相乘；2)在B(m,n)中根据坦克个体的主轴方向求得坦克主体区域的最小外接矩形，则该矩形的长边长度Length即为坦克个体的长度，矩形的短边长度Width即为坦克个体的宽度；3)计算得到几何结构特征，其中包括周长、面积、长宽比、形状复杂度、目标质心位置以及转动惯量：周长面积长宽比R＝Length/Width；形状复杂度C＝Length2/4πS；目标区域的质心位置转动惯量式中，r代表了目标像素点与质心之间的距离，4)计算得到灰度统计特征，其中包括质量、均值、方差系数、标准差、分形维数、加权填充比：质量均值方差系数标准差式中分别表示灰度对数和、灰度对数平方和；分形维数H＝(log10N1-log10N2)/(log10d1-log10d2)，该特征的计算方法是：用分割后的SAR图像切片构建一个保留了目标区域的K(这里取K＝50)个最亮像素点的二值图B2(m,n)，首先将一个大小为d1×d2的窗口在这个二值图中连续滑动，记下窗口中包含亮点的窗口总数记为N1，接着再用一个大小为d2×d2的窗口在这个二值图中连续滑动，记下在窗口中包含亮点的窗口总数记为N2；加权填充比特征选择模块，用以选择出最优特征子集，采用如下过程完成：1)计算类内距离类间距离以及类内类间距Ji:其中，i表示特征标号，ω表示坦克类别的标号，||Fi(ω)||2表示特征向量Fi(ω)的2范数，表示训练集样本的总体均值，Nω表示第ω类坦克的数量，N表示训练集中坦克总数，E表示期望，下标W、下标B分别表示类内、类间。2)计算得到归一化方差系数ρi(ω)：其中，i表示特征标号，ω表示坦克类别的标号，||Fi(ω)||2表示特征向量Fi(ω)的2范数，和E2[||Fi(ω)||2]分别表示特征的平方均值以及均值的平方。特征的方差系数ρi(ω)越小，表明该目标特征的稳定性越好；3)计算得到相关系数ri,j：其中，i,j表示特征标号，||Fi||2表示特征Fi的2范数，和分别表示Fi和Fj的均值，σi,i和σj,j分别表示Fi和Fj的标准差。由相关系数的性质可知，0≤ri,j≤1；如果两个特征完全不相关，ri,j＝0；如果两个特征完全相关，ri,j＝1；如果两个特征之间的相关性很低，即特征之间的信息冗余非常少，那么ri,j就会越接近0；反之，如果两个特征之间的相关性很高，即特征之间的信息冗余非常多，那么ri,j就会越接近1；4)通过上述得到的类内类间距、归一化方差系数、相关系数筛选出最优特征子集，构造最优输入特征向量；分类器训练模块，用以进行分类器训练，采用如下过程完成：1)从特征选择模块中采集N个SAR雷达图像xi作为训练样本，i＝1,2,…,N；2)对训练样本进行归一化处理，得到归一化样本3)将归一化后的训练样本重构，分别得到输入矩阵X和对应的输出矩阵Y:其中，D表示重构维数，D为自然数，且D<N,D的取值范围为50-70；4)将得到的X,Y代入如下分类器：Y＝wTΦ(x)+E其中，Φ(·)是核函数矩阵，w为权重系数矩阵，E为残差矩阵，核函数的似然函数为：其中，t为样本目标的分类值，σ2为坦克样本的方差，直接用最大似然估计的方法求解w,σ2，通常会导致严重的过拟合现象，使w中的元素大部分不是0，失去了稀疏性。为了避免这一现象，引入高斯先验来解决这一问题，即定义权重w的先验分布为依赖于超参数α的高斯分布：其中，α＝(α0,α1,…,αM)T为超参数，决定分类器的稀疏性。由贝叶斯准则，可得权重向量w的后验似然分布为：该权值的后验分布属于多变量高斯分布，其后验协方差和后验均值分别为：Σ＝(σ-2ΦTΦ+A)-1μ＝σ-2ΣΦTt其中，A＝diag(α0,α1,…,αM)为超参数的对角矩阵。通过最大化边缘似然函数化求得最大似然估计点αMP，如下：其中，I为单位矩阵。通过得到的α＝αMP代入式求得均值的点估计μMP，从而得到最后坦克分类的估计值t＝ΦμMP。采用自下向上系统来更新超参数α，即快速边缘最大化算法。在迭代过程中，大部分的αi会接近无穷大，即相应的wi为0，其基函数可以被删除，从而达到稀疏性。其他的αi会接近有限值，与之对应的非零μMP的数目也很少，这些元素称为相关向量。引入稀疏因子si和质量因子qi：其中，B＝σ2I+ΦA-1ΦT，表示该矩阵中去掉第i个基向量后得到的相应矩阵。si是衡量基向量是否在系统中存在的依据；qi是某一个基向量排除在系统外的误差调整尺度。可得：在训练的RVM坦克本文档来自技高网...

【技术保护点】
1.一种SAR雷达陆地坦克目标识别系统，其特征在于：包括SAR雷达、数据库以及上位机，SAR雷达、数据库和上位机依次相连，所述SAR雷达对陆地进行实时监测，并将SAR雷达获得的图像数据存储到所述的数据库中，所述的上位机包括图像预处理模块、特征提取模块、分类器训练模块、分类器训练模块和结果显示模块，所述图像预处理模块、特征提取模块、分类器训练模块、分类器训练模块和结果显示模块依次相连。

【技术特征摘要】
1.一种SAR雷达陆地坦克目标识别系统，其特征在于：包括SAR雷达、数据库以及上位机，SAR雷达、数据库和上位机依次相连，所述SAR雷达对陆地进行实时监测，并将SAR雷达获得的图像数据存储到所述的数据库中，所述的上位机包括图像预处理模块、特征提取模块、分类器训练模块、分类器训练模块和结果显示模块，所述图像预处理模块、特征提取模块、分类器训练模块、分类器训练模块和结果显示模块依次相连。2.根据权利要求1所述SAR雷达陆地坦克目标识别系统，其特征在于：所述图像预处理模块用以进行SAR雷达图像数据预处理，采用如下过程完成：1)从数据库中传来的SAR图像灰度级为L,f(x0,y0)为像素点(x0,y0)处的灰度值，g(x0,y0)为像素点(x0,y0)的N×N邻域内像素的平均值，其中x0,y0分别表示像素点的横坐标和纵坐标；2)通过计算满足f＝m和g＝n的像素数目h(m,n)，得到二维联合概率密度pmn：pmn＝p(m,n)＝h(m,n)/M其中，M表示图像像素的总数目；3)计算二维直方图的均值向量μ：4)分别计算图像中目标和背景出现的概率P0,1和均值向量μ0,1：其中，t、s、下标0、下标1分别表示f分割阈值、g分割阈值、目标区域、背景区域；5)计算类间方差BCV：BCV＝P0(μ0-μ)(μ0-μ)′+P1(μ1-μ)(μ1-μ)′；其中，μ表示均值向量，上标’表示矩阵的转置。6)最佳阈值即为使得BCV为最大值时的二维阈值向量[s0,t0]：3.根据权利要求1所述SAR雷达陆地坦克目标识别系统，其特征在于：所述特征提取模块用以进行坦克典型特征的提取，采用如下过程完成：1)从图像预处理模块传来的只包含一个坦克目标的SAR图像切片I(m,n)，其中只包含目标区域的二值图为B(m,n)，则只包含目标的图像T(m,n)：T(m,n)＝I(m,n)×B(m,n)其中，×表示对应像素相乘；2)在B(m,n)中根据坦克个体的主轴方向求得坦克主体区域的最小外接矩形，则该矩形的长边长度Length即为坦克个体的长度，矩形的短边长度Width即为坦克个体的宽度；3)计算得到几何结构特征，其中包括周长、面积、长宽比、形状复杂度、目标质心位置以及转动惯量：周长面积长宽比R＝Length/Width；形状复杂度C＝Length2/4πS；目标区域的质心位置转动惯量式中，r代表了目标像素点与质心之间的距离，4)计算得到灰度统计特征，其中包括质量、均值、方差系数、标准差、分形维数、加权填充比：质量均值方差系数标准差式中分别表示灰度对数和、灰度对数平方和；分形维数H＝(log10N1-log10N2)/(log10d1-log10d2)，该特征的计算方法是：用分割后的SAR图像切片构建一个保留了目标区域的K(这里取K＝50)个最亮像素点的二值图B2(m,n)，首先将一个大小为d1×d2的窗口在这个二值图中连续滑动，记下窗口中包含亮点的窗口总数记为N1，接着再用一个大小为d2×d2的窗口在这个二值图中连续滑动，记下在窗口中包含亮点的窗口总数记为N2；加权填充比4.根据权利要求1所述SAR雷达陆地坦克目标识别系统，其特征在于：所述特征选择模块用以选择出最优特征子集，采用如下过程完成：1)计算类内距离类间距离以及类内类间距Ji:其中，i表示特征标号，ω表示坦克类别的标号，||Fi(ω)||2表示特征向量Fi(ω)的2范数，表示训练集样本的总体均值，Nω表示第ω类坦克的数量，N表示训练集中坦克总数，E表示期望，下标W、下标B分别表示类内、类间。2)计算得到归一化方差系数ρi(ω)：ρi(ω)＝E[||Fi(ω)|...

【专利技术属性】
技术研发人员：刘兴高，吴俊，孙元萌，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33