CN104331706B - 基于rbm和svm的极化sar图像分类 - Google Patents
基于rbm和svm的极化sar图像分类 Download PDFInfo
- Publication number
- CN104331706B CN104331706B CN201410597098.3A CN201410597098A CN104331706B CN 104331706 B CN104331706 B CN 104331706B CN 201410597098 A CN201410597098 A CN 201410597098A CN 104331706 B CN104331706 B CN 104331706B
- Authority
- CN
- China
- Prior art keywords
- polarization
- classification
- sar image
- polarimetric sar
- pixel
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V30/00—Character recognition; Recognising digital ink; Document-oriented image-based pattern recognition
- G06V30/10—Character recognition
- G06V30/19—Recognition using electronic means
- G06V30/192—Recognition using electronic means using simultaneous comparisons or correlations of the image signals with a plurality of references
- G06V30/194—References adjustable by an adaptive method, e.g. learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于RBM和SVM的极化SAR图像分类方法,主要解决现有极化SAR图像分类方法分类精度不高的问题。其实现步骤为:(1)输入待分类的极化SAR图像,并进行精致极化Lee滤波;(2)基于极化相干矩阵、极化协方差矩阵和Cloude分解提取每个像素点的原始特征并归一化;(3)初始化并训练RBM,得到相关参数;(4)根据RBM重新学习到的特征利用SVM进行分类,得到分类结果。本发明与现有的一些方法相比,不仅充分考虑到图像的空间相关性,而且能够提取出更有利于极化SAR图像分类的特征,明显地提高了极化SAR图像分类精度,可用于对极化SAR图像的地物分类和目标识别。
Description
技术领域
本发明属于图像处理技术领域,涉及图像分类领域的应用,具体的说是一种基于RBM和SVM的极化SAR图像分类方法,可用于对极化SAR图像的地物分类和目标识别。
背景技术
极化SAR是一种高分辨主动式相干多通道微波遥感成像雷达,它是SAR的一个重要分支,具有全天候、全天时、分辨率高、可侧视成像等优点,可广泛应用于军事、农业、导航、地理监视等诸多领域。极化SAR能够获得更加丰富的目标信息,在国际遥感领域受到高度重视,因此极化SAR图像分类已成为极化SAR信息处理的一个重要研究方向。
现有的极化SAR图像分类方法可以分为有监督分类和无监督分类。
有监督分类方法包括:Lee等提出的基于复Wishart分布的极化协方差矩阵监督分类方法,由于极化协方差矩阵C可以通过线性变换得到极化相干矩阵T,这样就得到Wishart分类器,这种方法要求C或T矩阵的概率密度分布函数服从复Wishart分布,对数据分布要求严格;Heermann等提出的基于后向传播神经网络的极化SAR图像分类方法,这种方法收敛速度慢,且容易陷入局部最优。
无监督分类方法包括:Cloude等提出的H/α非监督分类,它是通过Cloude目标分解得到散射熵H和平均散射角α特征参数后,根据这两个参数的范围对目标进行八分类,这种方法分类边界固定导致区域的划分过于武断,且只利用H和α这两个参数,极化信息没有得到充分利用,导致分类准确度低;Lee等提出基于Cloude目标分解和Wishart分类器的H/α-Wishart非监督分类方法,它是在原始H/α分类基础上增加了Wishart迭代,弥补了H/α分类固定边界的缺陷,但这种方法不能很好的保持各类的极化散射特性;Lee等基于Freeman分解提出了一种极化SAR图像分类方法,它主要是根据Freeman分解获得平面散射功率、二面角散射功率和体散射功率的大小对极化数据进行划分,并对初始划分进行类别合并,最后再利用Wishart分类器重新划分,这种方法保持了各类的散射特性,但存在多类别的划分及合并,计算复杂度较高。
发明内容
本发明的目的在于针对已有技术的不足,提出一种基于RBM和SVM的极化SAR图像分类方法,充分利用图像的空间相关性和极化信息,对原始特征空间重新学习,提取更有效的特征进行分类,提高了分类精度。
为实现上述目的,本发明包括如下步骤:
步骤1,输入待分类的极化SAR图像,进行精致极化Lee滤波,去除斑点噪声,得到滤波后的极化SAR图像;
步骤2,基于滤波后的极化SAR图像每个像素点的极化相干矩阵T和极化协方差矩阵C,得到图像的原始特征空间;
步骤3,根据极化SAR图像的地物分布参考图中已标注的信息,随机选择每类的训练样本,得到训练样本集合:
本发明中采样得到的每类样本数目为600个;
步骤4,初始化限制玻耳兹曼机(Restricted Boltzmann Machine,RBM)的相关参数;
步骤5,将选取的训练样本分批次并归一化到[0.1,0.9]后训练RBM网络;
步骤6,重复步骤5,直至满足终止条件,本方法中选取的终止条件为达到最大迭代次数200次,然后得到RBM的参数θ=(W,a,b)以及隐藏层的输出特征p(h=1|v(0)):
其中,W是可见层与隐藏层之间的连接权重,a是可见层偏置,b是隐藏层偏置,p(h=1|v(0))表示在已知可见层状态v(0)时采样出隐藏层各个单元状态h=1的概率,并将此作为RBM的输出特征;
步骤7,将步骤6得到的隐藏层特征(即极化SAR图像训练样本的原始数据经过RBM重新学习得到的隐藏层特征)以及训练样本标签去训练SVM分类器;
步骤8,利用训练好的SVM分类器预测分类:
将待分类的极化SAR图像原始的测试数据归一化到[0.1,0.9]后,输入训练好的RBM,得到RBM隐藏层的特征,然后利用这些特征输入到训练好的SVM对待分类的极化SAR图像进行分类,得到每个像素点的类别;
步骤9,输出图像并计算分类精度。
本发明首先对极化SAR图像滤波,减少噪声影响;其次,充分提取图像的原始特征并进行归一化;然后,初始化并训练RBM网络,对图像的原始特征重新学习;最后,利用RBM网络隐藏层输出的特征训练SVM分类器,并预测分类以及计算分类精度。本发明与现有的技术相比具有以下优点:
1.本发明采用精致极化Lee滤波对极化相干矩阵和极化协方差矩阵中所有元素进行相干斑抑制,减少噪声影响,同时避免了极化通道之间的串扰,又保持了极化通道之间的极化信息和统计相关性;
2.本发明考虑到图像空间相关性,将邻域特征也作为中心像素的特征,不仅充分利用了极化信息,而且能够自适应的将相似度高的分为一类,从而提高分类精度;
3.本发明采用新的特征学习方法—RBM,将极化SAR图像的原始特征重新进行学习,分级抽象化,提取出更有利于极化SAR图像分类的特征;
4.本发明实现了在极化SAR图像原始特征的基础上,重新学习并提取新的特征,再引入SVM分类器进行分类,实验结果表明这种策略得到了更好的分类结果,明显地提高了分类精度。
附图说明
图1是本发明的总流程图;
图2是本发明仿真使用的一幅极化SAR图像;
图3是本发明仿真使用的极化SAR图像的地物分布参考图;
图4是对图2精致极化Lee滤波后,用现有的Wishart分类器的方法,利用训练样本得到各类的聚类中心,再利用Wishart距离重新聚类得到的分类仿真结果图;
图5是对图2精致极化Lee滤波后,用现有的基于目标分解—Cloude分解和Freeman分解提取特征,再利用SVM分类的方法得到的分类仿真结果图;
图6是本发明对图2的分类仿真结果图。
具体实施方式
参照图1,本发明的具体实现步骤如下:
步骤1,输入待分类的极化SAR图像,进行精致极化Lee滤波,去除斑点噪声,得到滤波后的极化SAR图像。
优选地,采用如下步骤:
1a)边缘检测及方向窗口选取:在待分类的极化SAR分类图像的span总功率图像(图像每一像素点的值等于该像素点的极化协方差矩阵对角线元素之和)上,设置滤波窗口大小为7*7,按照像素空间位置,从左到右、从上到下将滤波窗口依次分解为9个子窗口,其中子窗口大小为3*3,子窗口之间有重叠,计算每个子窗口的均值,获得3*3大小的均值窗口;在均值窗口内,用4方向(水平,垂直,45度,135度)边缘模板进行检测以确定边缘方向,即将均值窗口分别与4个模板进行加权,对得到的加权结果求绝对值,选出所有绝对值中的最大值,将该最大值对应的方向作为边缘方向;在3*3的均值窗口内,比较边缘方向两侧的像素元素与中心像素(位于3*3均值窗口中心的那个像素)的邻近性,即均值窗口在边缘方向两侧的值减去中心像素的值,选取结果绝对值最小的一侧作为方向窗口(即其滤波同质区);
1b)权值参数计算及滤波处理:在方向窗口内,采用下式:
计算精致极化Lee滤波的滤波系数b,(1)式中的y表示方向窗口内极化SAR总功率图像的像素值,var(y)表示方向窗口内像素值的总功率方差,表示y的数学期望,σv表示相干斑噪声的标准差。由滤波系数b,根据下式
计算滤波后中心像素的极化相干矩阵或极化协方差矩阵(2)式中的表示方向窗口内像素的极化相干矩阵或极化协方差矩阵的数学期望,Z表示中心像素的极化相干矩阵或极化协方差矩阵。
步骤2,基于滤波后的极化SAR图像每个像素点的极化相干矩阵T和极化协方差矩阵C,得到图像的原始特征空间。
2a)由于极化相干矩阵T和极化协方差矩阵C具有丰富的关于雷达目标的相位和幅度信息,而且都是哈密顿半正定矩阵,所以可提取维数为3*3的极化相干矩阵T的上三角位置的6个元素的模值,维数为3*3的极化协方差矩阵C的上三角位置的6个元素的实部和虚部(对角线元素虚部为0,不包括在内),以及对角线元素之和得到功率span;
2b)对极化相干矩阵T进行Cloude分解:
其中,TD=diag{λ1,λ2,λ3},λi(i=1,2,3)是极化相干矩阵的特征值,且λ1≥λ2≥λ3,U=[v1,v2,v3],vi是与λi相对应的单位特征矢量,()H表示这个矩阵的共轭转置;
接着,得到相应的特征:散射熵其中
平均散射角其中αi为散射体的内部自由度,
反熵A=(λ2-λ3)/(λ2+λ3),极化相干矩阵T的特征值λ1,λ2,λ3,继续作为每个像素点的原始特征。基于上述特征提取方法,每个像素点共22维原始特征[T11,T22,T33,abs(T12),abs(T13),abs(T23),C11,C22,C33,real(C12),imag(C12),real(C13),imag(C13),real(C23),imag(C23),span,H,α,A,λ1,λ2,λ3];
2c)考虑图像具有空间相关性,所以本方法将每个像素点的4邻域的每个像素点的22维特征也都作为其特征,图像边界采用镜像对称进行填充,加上步骤2b)中得到的22维原始特征,每个像素点共提取了110维的原始特征。
步骤3,根据极化SAR图像的地物分布参考图中已标注的信息,随机选择每类的训练样本,得到训练样本集合。
本发明中采样得到的每类样本数目为600个。
步骤4,初始化限制玻耳兹曼机(Restricted Boltzmann Machine,RBM)的相关参数。
优选地,采用如下步骤:
4a)为了并行运算和提高效率,一般将选取的训练样本分成包含10到100个样本的小批量数据,然后逐个运算,本方法中每批次数据包含每类的5个训练样本数据;设置可见层单元数目为极化SAR图像每个像素点原始特征的维数110,隐藏层单元数目设置为30,即nv=110,nh=30;
4b)初始化相关参数:本方法中将可见层单元偏置和隐藏层单元偏置、可见层与隐藏层之间的连接权重的学习率均设置为0.01,ηvb=ηhb=ηW=0.01;权重损失系数λ设置为0.0002;初始动量学习率和最终动量学习率分别设置为0.5和0.9,即ρini=0.5,ρfinal=0.9;k步对比散度(Contrastive Divergence,CD)算法中参数k=1;初始化可见层与隐藏层的连接权重W是服从均值为0,标准差为0.1的正态分布的随机数,隐藏层偏置b初始化为0,可见层偏置a初始化为[0,1]内的随机数,相应参数的增量均设置为0,即ΔW(0)=Δa(0)=Δb(0)=0。
步骤5,将选取的训练样本分批次并归一化到[0.1,0.9]后训练RBM网络。
5a)对于每次迭代,在每小批量数据上,调用CD-k算法,进行k次交替Gibbs采样,本方法只需1次交替Gibbs采样即可保证采集到的样本符合目标分布,具体为:
其中(3)式表示在已知可见层状态v(0)时采样出隐藏层的第j个单元的状态等于1的概率,v(0)表示可见层的输入状态,即训练样本原始特征向量归一化后的数据,表示可见层第i个单元的状态,表示隐藏层第j个单元的状态,bj表示隐藏层第j个单元的偏置,wj,i表示可见层第i个单元与隐藏层第j个单元的连接权重;
接着,产生[0,1]内的随机数rj,对进行二值化,即:若则否则,
其中(4)式表示在已知隐藏层状态时采样出可见层各个单元状态以此来代替目标分布,其中ai表示可见层第i个单元的偏置;
利用(5)、(6)、(7)式分别计算每次迭代时目标函数在每小批量数据上关于连接权重W、可见层偏置a、隐藏层偏置b的偏导数;
Δwi,j (iter)=ρ*Δwi,j (iter-1)+ηW*(Δwi,j (iter)/nblock-λ*wi,j (iter-1)) (8)
Δai (iter)=ρ*Δai (iter-1)+ηvb*(Δai (iter)/nblock) (9)
Δbi (iter)=ρ*Δbi (iter-1)+ηhb*(Δbi (iter)/nblock) (10)
本方法采用(8)、(9)、(10)式更新当前迭代(iter≥1)每小批量数据上的平均偏导数,其中nblock为小批量数据的数目,ρ为动量学习率,本方法中,当每次迭代次数大于30次时,ρ=ρfinal=0.9,否则,ρ=ρini=0.5;
5b)利用梯度上升法更新RBM参数θ=(W,a,b),具体为:
W(iter)=W(iter-1)+ΔW(iter)
a(iter)=a(iter-1)+Δa(iter) (11)
b(iter)=b(iter-1)+Δb(iter)
按上式更新当前迭代每小批量数据上的连接权重、可见层偏置和隐藏层偏置。
步骤6,重复步骤5,直至满足终止条件,本方法中选取的终止条件为达到最大迭代次数200次,然后得到RBM的参数θ=(W,a,b)以及隐藏层的输出特征p(h=1|v(0)):
其中,W是可见层与隐藏层之间的连接权重,a是可见层偏置,b是隐藏层偏置,p(h=1|v(0))表示在已知可见层状态v(0)时采样出隐藏层各个单元状态h=1的概率,并将此作为RBM的输出特征。
步骤7,将步骤6得到的隐藏层特征(即极化SAR图像训练样本的原始数据经过RBM重新学习得到的隐藏层特征)以及训练样本标签去训练SVM分类器。
支持向量机(Support Vector Machine,SVM)是一种有效的监督分类器,根据间隔最大化准则构造最优分类超平面进行分类,具有良好的泛化能力;
针对输入空间线性可分问题,SVM通过求解其对偶优化问题:
s.t.αi≥0
其中,αi是Lagrange乘子,xi,xj是训练样本数据,yi,yj是训练样本标签,m是训练样本数目,<·,·>表示欧式内积;
求解出Lagrange乘子α,然后根据求解出超平面的法向量w,最后在α和w以求解出的情况下,可通过求解出超平面的截距b,其中T表示转置;
这样,最优分类超平面可表示为
该方法引入核函数可将输入空间中线性不可分的问题转化到高维空间,使其在新空间中线性可分,典型的核函数包括多项式核函数、径向基核函数和S型核函数等,本方法采用的核函数为径向基核函数;对于多类的分类问题,可以在2类问题的基础上进行推广,其中具有代表性的2种算法为一对多算法和一对一算法,本方法采用的分类算法为一对一算法。
步骤8,利用训练好的SVM分类器预测分类。
将待分类的极化SAR图像原始的测试数据归一化到[0.1,0.9]后,输入训练好的RBM,得到RBM隐藏层的特征,然后利用这些特征输入到训练好的SVM对待分类的极化SAR图像进行分类,得到每个像素点的类别。
步骤9,输出图像并计算分类精度。
9a)利用SVM预测分类的像素类别,将R(红色)、G(绿色)、B(蓝色)作为三基色,按照三基色上色法进行上色,得到上色后的极化SAR图像,然后将其输出;
9b)将极化SAR图像得到的像素类别与真实地物类别进行比较,将类别一致的像素个数与全部像素个数的比值作为极化SAR图像的分类精度。
本发明的效果可以通过下面的实验仿真进一步说明:
1、实验条件与方法
硬件平台为:Intel Core2 Duo CPU [email protected]、1.95GB RAM;
软件平台为:MATLAB R2009a;
实验方法:分别为本发明和现有的Wishart分类器的方法以及基于Cloude分解和Freeman分解提取特征,再利用SVM分类的方法,其中现有的这两种方法都是极化SAR图像分类中经典的方法。
2、仿真内容与结果
本发明和现有的Wishart分类器的方法以及基于Cloude分解和Freeman分解提取特征,再利用SVM分类的方法分别对图2所示的Flevoland农田极化SAR图像进行分类仿真,分类结果见图4到图6,其中,图4为Wishart分类器方法分类结果图,图5为基于Cloude分解和Freeman分解提取特征,再利用SVM分类的结果图,图6为本发明的分类结果图。
仿真实验中,根据图3所示的地物分布参考图,随机选取每类600个训练样本,在剩下的样本上计算精度,应用各类分类精度和总分类精度作为评价指标。
评价结果如表1所示,其中,Alg1是Wishart分类器的方法,Alg2是基于Cloude分解和Freeman分解提取特征,再利用SVM分类的方法,Alg3是本发明的方法。
表1.本发明和两种对比方法在仿真实验中得到的各类分类精度和总分类精度
3.实验结果分析
对比Wishart分类器方法得到的分类结果图4、基于Cloude分解和Freeman分解提取特征,再利用SVM分类的方法得到的分类结果图5和本发明得到的分类结果图6以及表1可以看出,图6所示的本发明相对于其他两种方法得到了更好的实验结果,图6不仅杂散点少,而且保持了区域一致性与连贯性,每类分类精度都高于两种对比方法,总分类精度明显地提高;图4所示的Wishart分类器方法得到的分类结果虽然区域边缘划分相对平滑,但杂散点较多,而且有严重的错分现象;图5所示的基于Cloude分解和Freeman分解提取特征,再利用SVM分类的方法能够有效的正确分类,但是杂散点较多,丢失了图像的细节信息。
综上所述,本发明提出的基于RBM和SVM的极化SAR图像分类方法能够明显地提高极化SAR图像分类精度,不仅有效的去除了噪声影响,保持了区域一致性与连贯性,而且可用于对各种极化SAR图像分类。
Claims (3)
1.一种基于RBM和SVM的极化SAR图像分类方法,包括如下步骤:
步骤1,输入待分类的极化SAR图像,进行精致极化Lee滤波,去除斑点噪声,得到滤波后的极化SAR图像;
其中步骤1所述的精致极化Lee滤波,按照如下步骤进行:
1a)边缘检测及方向窗口选取:在待分类的极化SAR分类图像的span总功率图像上,设置滤波窗口大小为7*7,按照像素空间位置,从左到右、从上到下将滤波窗口依次分解为9个子窗口,其中子窗口大小为3*3,子窗口之间有重叠,计算每个子窗口的均值,获得3*3大小的均值窗口;在均值窗口内,用水平、垂直、45度、135度这4方向边缘模板进行检测以确定边缘方向,即将均值窗口分别与4个模板进行加权,对得到的加权结果求绝对值,选出所有绝对值中的最大值,将该最大值对应的方向作为边缘方向;在3*3的均值窗口内,比较边缘方向两侧的像素元素与中心像素的邻近性,即均值窗口在边缘方向两侧的值减去中心像素的值,选取结果绝对值最小的一侧作为方向窗口;
1b)权值参数计算及滤波处理:在方向窗口内,采用下式:
计算精致极化Lee滤波的滤波系数b,(1)式中的y表示方向窗口内极化SAR总功率图像的像素值,var(y)表示方向窗口内像素值的总功率方差,表示y的数学期望,σv表示相干斑噪声的标准差,由滤波系数b,根据下式
计算滤波后中心像素的极化相干矩阵或极化协方差矩阵(2)式中的表示方向窗口内像素的极化相干矩阵或极化协方差矩阵的数学期望,Z表示中心像素的极化相干矩阵或极化协方差矩阵;
步骤2,基于滤波后的极化SAR图像每个像素点的极化相干矩阵T和极化协方差矩阵C,提取图像的原始特征空间;
步骤3,根据极化SAR图像的地物分布参考图中已标注的信息,随机选择每类的训练样本,得到训练样本集合:
步骤4,初始化限制玻耳兹曼机(Restricted Boltzmann Machine,RBM)的相关参数,包括:
4a)将选取的训练样本分成包含10到100个样本的小批量数据,然后逐个运算,本方法中每批次数据包含每类的5个训练样本数据;设置可见层单元数目为极化SAR图像每个像素点原始特征的维数110,隐藏层单元数目设置为30,即nv=110,nh=30;
4b)初始化相关参数:本方法中将可见层单元偏置和隐藏层单元偏置、可见层与隐藏层之间的连接权重的学习率均设置为0.01,ηvb=ηhb=ηW=0.01;权重损失系数λ设置为0.0002;初始动量学习率和最终动量学习率分别设置为0.5和0.9,即ρini=0.5,ρfinal=0.9;k步对比散度(Contrastive Divergence,CD)算法中参数k=1;初始化可见层与隐藏层的连接权重W是服从均值为0,标准差为0.1的正态分布的随机数,隐藏层偏置b初始化为0,可见层偏置a初始化为[0,1]内的随机数,相应参数的增量均设置为0,即;
步骤5,将选取的训练样本分批次并归一化到[0.1,0.9]后训练RBM网络,包括:
5a)对于每次迭代,在每小批量数据上,调用CD-k算法,进行k次交替Gibbs采样,本方法只需1次交替Gibbs采样即可保证采集到的样本符合目标分布,具体为:
其中(3)式表示在已知可见层状态v(0)时采样出隐藏层的第j个单元的状态等于1的概率,v(0)表示可见层的输入状态,即训练样本原始特征向量归一化后的数据,表示可见层第i个单元的状态,表示隐藏层第j个单元的状态,bj表示隐藏层第j个单元的偏置,wj,i表示可见层第i个单元与隐藏层第j个单元的连接权重;
接着,产生[0,1]内的随机数rj,对进行二值化,即:若则否则,
其中(4)式表示在已知隐藏层状态时采样出可见层各个单元状态以此来代替目标分布,其中ai表示可见层第i个单元的偏置;
利用(5)、(6)、(7)式分别计算每次迭代时目标函数在每小批量数据上关于连接权重W、可见层偏置a、隐藏层偏置b的偏导数;
△wi,j (iter)=ρ*△wi,j (iter-1)+ηW*(△wi,j (iter)/nblock-λ*wi,j (iter-1)) (8)
△ai (iter)=ρ*△ai (iter-1)+ηvb*(△ai (iter)/nblock) (9)
△bi (iter)=ρ*△bi (iter-1)+ηhb*(△bi (iter)/nblock) (10)
本方法采用(8)、(9)、(10)式更新当前迭代(iter≥1)每小批量数据上的平均偏导数,其中nblock为小批量数据的数目,ρ为动量学习率,本方法中,当每次迭代次数大于30次时,ρ=ρfinal=0.9,否则,ρ=ρini=0.5;
5b)利用梯度上升法更新RBM参数θ=(W,a,b),具体为:
按上式更新当前迭代每小批量数据上的连接权重、可见层偏置和隐藏层偏置;
步骤6,重复步骤5,直至满足终止条件,然后得到RBM的参数θ=(W,a,b)以及隐藏层的输出特征p(h=1|v(0)):
其中,W是可见层与隐藏层之间的连接权重,a是可见层偏置,b是隐藏层偏置,p(h=1|v(0))表示在已知可见层状态v(0)时采样出隐藏层各个单元状态h=1的概率,并将此作为RBM的输出特征;
步骤7,将步骤6得到的极化SAR图像训练样本的原始数据经过RBM重新学习得到的隐藏层特征以及训练样本标签去训练SVM分类器,包括:
7a)支持向量机(Support Vector Machine,SVM)是一种有效的监督分类器,根据间隔最大化准则构造最优分类超平面进行分类,具有良好的泛化能力;
针对输入空间线性可分问题,SVM通过求解其对偶优化问题:
s.t.αi≥0
其中,αi是Lagrange乘子,xi,xj是训练样本数据,yi,yj是训练样本标签,m是训练样本数目,<·,·>表示欧式内积;
求解出Lagrange乘子α,然后根据求解出超平面的法向量w,最后在α和w以求解出的情况下,可通过求解出超平面的截距b,其中T表示转置;
这样,最优分类超平面可表示为
7b)该方法引入核函数可将输入空间中线性不可分的问题转化到高维空间,使其在新空间中线性可分;
步骤8,利用训练好的SVM分类器预测分类:
将待分类的极化SAR图像原始的测试数据归一化到[0.1,0.9]后,输入训练好的RBM,得到RBM隐藏层的特征,然后利用这些特征输入到训练好的SVM对待分类的极化SAR图像进行分类,得到每个像素点的类别;
步骤9,输出图像并计算分类精度。
2.根据权利要求1所述的基于RBM和SVM的极化SAR图像分类方法,其中步骤2所述的提取原始特征,按照如下步骤进行:
2a)由于极化相干矩阵T和极化协方差矩阵C具有丰富的关于雷达目标的相位和幅度信息,而且都是哈密顿半正定矩阵,所以可提取维数为3*3的极化相干矩阵T的上三角位置的6个元素的模值,维数为3*3的极化协方差矩阵C的上三角位置的6个元素的实部和虚部,以及对角线元素之和得到功率span;
2b)对极化相干矩阵T进行Cloude分解:
其中,TD=diag{λ1,λ2,λ3},λi(i=1,2,3)是极化相干矩阵的特征值,且λ1≥λ2≥λ3,U=[v1,v2,v3],vi是与λi相对应的单位特征矢量,()H表示这个矩阵的共轭转置;
接着,得到相应的特征:散射熵其中
平均散射角其中αi为散射体的内部自由度,
反熵A=(λ2-λ3)/(λ2+λ3),极化相干矩阵T的特征值λ1,λ2 ,λ3,继续作为每个像素点的原始特征,基于上述特征提取方法,每个像素点共22维原始特征[T11,T22,T33,abs(T12),abs(T13),abs(T23),C11,C22,C33,real(C12),imag(C12),real(C13),imag(C13),real(C23),imag(C23),span,H,α,A,λ1,λ2,λ3];
2c)考虑图像具有空间相关性,所以本方法将每个像素点的4邻域的每个像素点的22维特征也都作为其特征,图像边界采用镜像对称进行填充,加上步骤2b)中得到的22维原始特征,每个像素点共提取了110维的原始特征。
3.根据权利要求1所述的基于RBM和SVM的极化SAR图像分类方法,其中步骤9所述的输出图像及计算精度,按照如下方式得到:
9a)利用SVM预测分类的像素类别,将红色、绿色、蓝色作为三基色,按照三基色上色法进行上色,得到上色后的极化SAR图像,然后将其输出;
9b)将极化SAR图像得到的像素类别与真实地物类别进行比较,将类别一致的像素个数与全部像素个数的比值作为极化SAR图像的分类精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410597098.3A CN104331706B (zh) | 2014-10-29 | 2014-10-29 | 基于rbm和svm的极化sar图像分类 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410597098.3A CN104331706B (zh) | 2014-10-29 | 2014-10-29 | 基于rbm和svm的极化sar图像分类 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104331706A CN104331706A (zh) | 2015-02-04 |
CN104331706B true CN104331706B (zh) | 2018-03-16 |
Family
ID=52406427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410597098.3A Active CN104331706B (zh) | 2014-10-29 | 2014-10-29 | 基于rbm和svm的极化sar图像分类 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104331706B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104751183B (zh) * | 2015-04-08 | 2018-03-06 | 西安电子科技大学 | 基于张量mpca的极化sar图像分类方法 |
CN104965822B (zh) * | 2015-07-29 | 2017-08-25 | 中南大学 | 一种基于计算机信息处理技术的中文文本情感分析方法 |
CN105184309B (zh) * | 2015-08-12 | 2018-11-16 | 西安电子科技大学 | 基于cnn和svm的极化sar图像分类 |
CN105069459B (zh) * | 2015-08-18 | 2018-06-19 | 电子科技大学 | 一种针对高分辨率sar图像地物类型提取方法 |
CN105512661B (zh) * | 2015-11-25 | 2019-02-26 | 中国人民解放军信息工程大学 | 一种基于多模态特征融合的遥感影像分类方法 |
CN105574616A (zh) * | 2015-12-16 | 2016-05-11 | 国网山东省电力公司电力科学研究院 | 一种架空输电线路动态容量预测方法 |
CN105894013B (zh) * | 2016-03-29 | 2019-04-23 | 西安电子科技大学 | 基于cnn和smm的极化sar图像分类方法 |
CN106485673B (zh) * | 2016-09-19 | 2018-12-18 | 电子科技大学 | 一种对海面sar图像的滤波方法 |
CN106909939B (zh) * | 2017-02-20 | 2020-09-11 | 中国人民解放军国防科学技术大学 | 一种结合旋转域极化零角特征的极化sar地物分类方法 |
CN106897741B (zh) * | 2017-02-20 | 2019-06-14 | 中国人民解放军国防科学技术大学 | 结合旋转域极化相干特征的极化sar地物分类方法 |
CN107528824B (zh) * | 2017-07-03 | 2020-08-04 | 中山大学 | 一种基于二维度稀疏化的深度信念网络入侵检测方法 |
CN107563447B (zh) * | 2017-09-08 | 2020-11-10 | 哈尔滨工业大学 | 一种遥感图像中目标到目标部位的分级识别方法 |
CN108509989B (zh) * | 2018-03-26 | 2020-04-21 | 西安电子科技大学 | 基于高斯选控玻尔兹曼机的hrrp识别方法 |
CN108898101B (zh) * | 2018-06-29 | 2021-09-28 | 西安电子科技大学 | 基于素描图和先验约束的高分辨sar图像路网检测方法 |
CN110472627B (zh) * | 2019-07-02 | 2022-11-08 | 五邑大学 | 一种端到端的sar图像识别方法、装置及存储介质 |
CN110703221A (zh) * | 2019-10-16 | 2020-01-17 | 艾索信息股份有限公司 | 基于极化特征的城市低空小目标分类识别*** |
CN111915559B (zh) * | 2020-06-30 | 2022-09-20 | 电子科技大学 | 基于svm分类可信度的机载sar图像质量评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722883A (zh) * | 2012-04-16 | 2012-10-10 | 上海交通大学 | 一种具有空间自适应性的极化sar图像分割方法 |
CN102402685B (zh) * | 2011-11-10 | 2013-06-26 | 西安电子科技大学 | 基于Gabor特征的三马尔可夫场SAR图像分割方法 |
CN103366184A (zh) * | 2013-07-23 | 2013-10-23 | 武汉大学 | 基于混合分类器的极化sar数据分类方法及*** |
CN104112144A (zh) * | 2013-12-17 | 2014-10-22 | 深圳市华尊科技有限公司 | 人车识别方法及装置 |
CN104112138A (zh) * | 2013-12-17 | 2014-10-22 | 深圳市华尊科技有限公司 | 物体颜色分类方法及装置 |
-
2014
- 2014-10-29 CN CN201410597098.3A patent/CN104331706B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102402685B (zh) * | 2011-11-10 | 2013-06-26 | 西安电子科技大学 | 基于Gabor特征的三马尔可夫场SAR图像分割方法 |
CN102722883A (zh) * | 2012-04-16 | 2012-10-10 | 上海交通大学 | 一种具有空间自适应性的极化sar图像分割方法 |
CN103366184A (zh) * | 2013-07-23 | 2013-10-23 | 武汉大学 | 基于混合分类器的极化sar数据分类方法及*** |
CN104112144A (zh) * | 2013-12-17 | 2014-10-22 | 深圳市华尊科技有限公司 | 人车识别方法及装置 |
CN104112138A (zh) * | 2013-12-17 | 2014-10-22 | 深圳市华尊科技有限公司 | 物体颜色分类方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104331706A (zh) | 2015-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104331706B (zh) | 基于rbm和svm的极化sar图像分类 | |
CN105184309B (zh) | 基于cnn和svm的极化sar图像分类 | |
CN105469034B (zh) | 基于加权式鉴别性稀疏约束非负矩阵分解的人脸识别方法 | |
Zhao et al. | Spectral-spatial classification of hyperspectral imagery based on stacked sparse autoencoder and random forest | |
CN107085716A (zh) | 基于多任务生成对抗网络的跨视角步态识别方法 | |
Wang et al. | Hyperspectral anomaly detection using ensemble and robust collaborative representation | |
Montazer et al. | An improved radial basis function neural network for object image retrieval | |
Chen et al. | Spectral–spatial preprocessing using multihypothesis prediction for noise-robust hyperspectral image classification | |
Liu et al. | Enhancing spectral unmixing by local neighborhood weights | |
CN103971123B (zh) | 基于线性回归与Fisher判别字典学习的高光谱图像分类方法 | |
CN102982338B (zh) | 基于谱聚类的极化sar图像分类方法 | |
CN104268593B (zh) | 一种小样本情况下多稀疏表示的人脸识别方法 | |
CN108537102B (zh) | 基于稀疏特征与条件随机场的高分辨sar图像分类方法 | |
CN105894018B (zh) | 基于深度多示例学习的极化sar图像分类方法 | |
CN105894013B (zh) | 基于cnn和smm的极化sar图像分类方法 | |
CN105138970A (zh) | 基于空间信息的极化sar图像分类方法 | |
Ghosh et al. | A novel approach to neuro-fuzzy classification | |
Manivannan et al. | Particulate matter characterization by gray level co-occurrence matrix based support vector machines | |
CN110458192B (zh) | 基于视觉显著性的高光谱遥感图像分类方法及*** | |
CN103793908A (zh) | 一种基于脑核磁共振图像多维度纹理建立预测模型的方法 | |
CN106096652A (zh) | 基于稀疏编码和小波自编码器的极化sar图像分类方法 | |
CN102999762B (zh) | 基于Freeman分解和谱聚类的极化SAR图像分类方法 | |
CN104282008B (zh) | 对图像进行纹理分割的方法和装置 | |
CN103632160A (zh) | 一种融合多尺度形态学特征的组合核函数rvm高光谱分类方法 | |
CN102999761A (zh) | 基于Cloude分解和K-wishart分布的极化SAR图像分类方法 |
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 |