CN104166859B - 基于ssae和fsals‑svm极化sar图像分类 - Google Patents
基于ssae和fsals‑svm极化sar图像分类 Download PDFInfo
- Publication number
- CN104166859B CN104166859B CN201410395368.2A CN201410395368A CN104166859B CN 104166859 B CN104166859 B CN 104166859B CN 201410395368 A CN201410395368 A CN 201410395368A CN 104166859 B CN104166859 B CN 104166859B
- Authority
- CN
- China
- Prior art keywords
- mrow
- sae
- classification
- mtd
- input
- 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
Landscapes
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明的目的是提供一种基于SSAE和FSALS‑SVM的极化SAR图像分类方法。它利用栈式稀疏自动编码器(SSAE)的多隐层结构,获得具有对原始数据更本质的刻画能力并且更适合分类的深度特征,并用能够获得问题稀疏解的快速稀疏逼近最小二乘支持向量机(FSALS‑SVM)代替传统深度学习中常用的Softmax,与SSAE相结合,提高了对极化SAR图像的分类精度,并在一定程度上克服了传统基于像素的极化SAR图像分类方法受相干斑噪声影响大的缺点,保证了分类结果图中匀质区域的连贯性。
Description
技术领域
本发明涉及图像处理领域,针对极化SAR图像分类问题,提出了一种基于栈式稀疏自动编码器(SSAE)和快速稀疏逼近最小二乘支持向量机(FSALS-SVM)的极化SAR图像分类方法,可用于开展航空航天影像、天文学影像、军事等领域的数字图像预处理。
背景技术
合成孔径雷达(SAR)由于其全天时、全天候的工作能力,以及其分辨率高、能有效识别伪装和穿透掩盖物等特点而被广泛应用于遥感和地图测绘等领域。在近二十年来,极化SAR已被证明能够获得比传统的单极化SAR更加丰富的地物信息。目前,全球许多极化SAR系统,如TerraSAR-X、RADARSAT-2、ALOS-PALSAR等,已经提供了大量的极化SAR数据以供研究。然而,对这些大规模、复杂的数据进行人为手动研究是不现实的。因此,对极化SAR数据进行全自动或半自动研究是迫切需要的,其中极化SAR图像分类是极化SAR信息处理的一个重要分支。
极化SAR图像分类作为极化SAR图像解译的重要步骤,许多学者都对其进行了深入的研究,所提出的方法大致可以分为以下三类:1、基于极化散射机制的分类方法,它的特点是利用不同类的先验知识对图像进行分类;2、基于极化SAR图像统计特性的分类方法,以最大似然法和最大后验概率法为代表;3、基于二者结合的分类方法。然而,所有这些方法都可归结成基于像素的方法。基于像素的方法即:仅利用单个像素本身的特性对图像进行分类。这类方法在确定当前像素类别时,将每个像素都看成是独立的,不受周围像素的影响,能较好的保留图像细节。
发明内容
本发明的目的在于针对已有技术的不足,提出一种基于SSAE和FSALS-SVM的极化SAR图像分类方法,以提高分类效果,使分类结果更吻合真实地物。
为实现上述目的,本发明包括如下步骤:
步骤1,对待输入的所有极化SAR图像数据进行预处理,获得无标记训练样本和带标记训练样本;
步骤2,将每个稀疏自动编码器(SAE)均作为构成栈式稀疏自动编码器(SSAE)的一个基本单元,并将前一个SAE的输出作为后一个SAE的输入来逐个训练SAE,以获得每个SAE的模型参数,其中,步骤1中得到的所有无标记训练样本作为首个SAE的输入,优选地,本发明中的SAE个数为2;
步骤3,将步骤2中获得的多个SAE级联在一起构成SSAE,并将步骤1中得到的带标记训练样本输入该SSAE,进行前馈传导计算,获得这些带标记训练样本的深度特征;
步骤4,将深度特征进行归一化处理后,输入到快速稀疏逼近最小二乘支持向量机(FSALS-SVM),通过迭代训练获得训练后的FSALS-SVM分类器;
步骤5,将预处理后的待分类极化SAR图像数据,即步骤1中得到的无标记训练样本,输入到栈式稀疏自动编码器,获得数据的深度特征,并将这些特征输入到训练后的FSALS-SVM分类器,获得图像的最终分类结果。
本发明与现有的技术相比具有以下优点:
1.本发明在获得图像分类特征时,利用了栈式稀疏自动编码器,输入仅为经过简单预处理后的原始极化SAR协方差数据,避免了传统在对极化SAR图像进行分类前,复杂的基于极化目标的特征分解过程;
2.由于栈式稀疏自动编码器具有多隐层的人工神经网络结构,它具有优异的特征学习能力,学习得到的特征对数据有更本质的刻画,从而有利于可视化和分类;
3.本发明用FSALS-SVM代替了传统深度学习常用的Softmax分类器,它能够获得问题的稀疏解,克服了LS-SVM因计算复杂度高而不利于大规模数据分类的缺点,降低了分类所耗的时间复杂度,与栈式稀疏自动编码器结合,提高了分类精度,并在一定程度上加强了分类结果中匀质区域的连贯性。
附图说明
图1是本发明的主流程图;
图2是本发明仿真使用的极化SAR数据的PauliRGB合成图;
图3是用现有基于SVM的分类方法对图2的分类结果图;
图4是用现有基于Wishart分布的分类方法对图2的分类结果图;
图5是用本发明对图2的分类结果图。
具体实施方式
参照图1,本发明的具体实现步骤如下:
步骤1,对待输入的所有极化SAR图像数据进行预处理,获得无标记训练样本和带标记训练样本。
(1a)对待输入的所有极化SAR图像数据进行滑窗大小为7×7的Lee滤波;
(1b)Lee滤波后的极化SAR图像的任一像素点n均可表示为一个3×3协方差矩阵Mn:
其中,n∈[1,2,…,N],N为该极化SAR图像所包含的像素点个数,矩阵Mn中的大写字母A~I均为实数,可将这些字母组成列向量tn=[A B C D E F G H I]T,对每个像素点n所表示的列向量tn依次排放,构成整个待分类样本集;
(1c)对(1b)中得到的待分类样本集进行简单地处理使其利于分类,在本发明中,优选地,根据极化SAR图像数据本身的特点,直接对该待分类样本集统一乘以一个较大的整数Q就可获得比较好的分类结果,优选地,Q=500;
(1d)将(1c)中简单地处理后的待分类样本集中的所有样本均作为无标记训练样本,并根据真实地物类别,随机选取无标记样本中的一部分(优选地,2%)进行标记,作为带标记训练样本。
步骤2,将每个稀疏自动编码器(SAE)均作为构成栈式稀疏自动编码器(SSAE)的一个基本单元,并将前一个SAE的输出作为后一个SAE的输入来逐个训练SAE,以获得每个SAE的模型参数,其中,步骤1中得到的所有无标记训练样本作为首个SAE的输入,优选地,本发明中所需的SAE个数为两个。
具体实现步骤如下:
(2a)将步骤1中得到的所有无标记训练样本作为首个SAE的输入x,分别设置该SAE输入单元、隐藏单元和输出单元的个数。优选地,该SAE的输入单元和输出单元的个数均为9,隐藏单元的个数为50。
(2b)随机初始化模型参数W1 (1)、W1 (2)、其中,W1 (1)表示首个SAE中输入单元和隐藏单元间的连接权重,W1 (2)表示首个SAE中隐藏单元和输出单元间的连接权重,表示首个SAE中隐藏单元的偏置项,表示首个SAE中输出单元的偏置项。计算隐藏单元和输出单元的激活值,并利用梯度下降法,最小化SSAE的整体代价函数Jsparse,从而获得训练后的模型参数具体实现步骤如下:
①随机初始化连接权重W1 (1)、W1 (2),使它们均为范围在内的数,偏置项为0,其中nin表示当前SAE中输入单元的个数,nout表示当前SAE中输出单元的个数;
②利用如下前项传递函数F:
a(l+1)=F(W(l)a(l)+b(l))
进行前馈传导计算,分别得到隐藏单元和输出单元的激活值a(2)、a(3)。其中,对任一SAE均有,l∈{1,2}即l可以取值1或2,a(1)表示输入单元的激活值,a(2)表示隐藏单元的激活值,a(3)表示输出单元的激活值;
③定义SSAE的整体代价函数Jsparse:
其中,第一项为均方差项,第二项为权重衰减项,其目的在于减小权重幅度防止过度拟合,第三项 为稀疏性惩罚因子,用于将隐藏单元的平均活跃度保持在较小的范围内,W表示输入单元和隐藏单元之间的连接权重W(1)与隐藏单元和输出单元之间的连接权重W(2)的转置的集合,即 表示矩阵W(2)的转置,m表示输入SAE的样本个数,和y(i)分别表示输入第i个样本时得到的实际输出结果和理想输出结果,的上标2表示平方,表示实际输出结果与理想输出结果之间误差的平方,的上标2表示平方,下标F表示F范数,则表示W的F范数的平方,λ为权重衰减参数,β为控制稀疏性惩罚因子的权重,s为隐藏单元的个数,ρ为稀疏性参数,表示第j个隐藏单元的平均活跃度,它可由下式计算得到:
其中,表示在输入第i个样本x(i)的时隐藏单元j的激活度,上标2表示隐藏单元。利用反向传波算法,计算整体代价函数Jsparse的梯度,优选地,ρ=0.1、β=3、λ=0.003;
④通过梯度下降法求解代价函数Jsparse的最小化问题,从而获得训练后的模型参数
(2c)重新将步骤1中得到的所有无标记训练样本作为首个SAE的输入x,利用求得的模型参数可以计算出该SAE隐藏单元的激活值,并将该激活值作为第二个SAE的输入,设置该SAE输入单元、隐藏单元和输出单元的个数。优选地,该SAE的输入单元和输出单元的个数均为50,隐藏单元的个数为100;
(2d)按照步骤(2b)的方式对第二个SAE的输入进行训练,即可获得第二个SAE训练后的模型参数下标2代表第二个SAE;
(2e)令r表示第r个SAE,若r小于所需SAE的个数R,则对于第r个SAE的模型参数Wr (1)、Wr (2)、将第r-1个SAE中隐藏单元的激活值作为第r个SAE的输入,并设置第r个SAE输入单元、隐藏单元和输出单元的个数,按照步骤(2b)的方式对第r个SAE的输入进行训练,即可获得第r个SAE训练后的模型参数否则(即r等于所需SAE的个数R),停止训练。优选地,R=2。
步骤3,将步骤2中获得的R个稀疏自动编码器(SAE)级联在一起,构成栈式稀疏自动编码器(SSAE),并将步骤1中得到的带标记训练样本输入该SSAE,进行前馈传导计算,获得这些带标记训练样本的深度特征。
步骤4,将深度特征进行归一化处理后,输入到快速稀疏逼近最小二乘支持向量机(FSALS-SVM),通过迭代训练获得训练后的FSALS-SVM分类器。
(4a)设置FSALS-SVM分类器中所用的核函数,该发明中,优选地,核函数为径向基函数;
(4b)计算步骤3中获得的这些带标记训练样本的深度特征的均值mean和方差std,并将这些深度特征的取值范围规整化为零均值和单位方差,得到归一化的深度特征;
(4c)将第k类(k≥1)训练样本标记为1,其余类别均标记为-1,变为第k类的二分类问题;
(4d)输入归一化的深度特征,并利用(4c)中得到的第k类的二分类问题训练FSALS-SVM分类器,并在获得第k类参数后,令k的值加1;
(4e)若k大于训练样本类别数K,则停止训练,并将获得的第1类到第K类参数共K类参数依次排列后所构成的集合,作为训练后的FSALS-SVM分类器参数,获得训练后的FSALS-SVM分类器;否则(即k小于或等于训练样本类别数K),重复步骤(4c)~(4e),本发明的实施例中,根据真实地物的类别总数确定K=9。
步骤5,将预处理后的待分类极化SAR图像数据,即步骤1中得到的无标记训练样本,输入到SSAE,获得这些无标记训练样本的深度特征,并将这些深度特征输入到训练后的FSALS-SVM分类器,获得图像的最终分类结果。
(5a)将步骤1中得到的无标记训练样本输入到SSAE,通过前馈传导计算,获得这些无标记训练样本的深度特征;
(5b)将(5a)中获得的深度特征的取值范围规整化为零均值和单位方差,得到归一化的深度特征;
优选地,利用步骤(4b)中得到的均值mean和方差std,将(5a)中获得的深度特征的取值范围规整化为零均值和单位方差,得到归一化的深度特征;
(5c)将归一化深度特征输入到训练后的FSALS-SVM分类器,获得最终分类结果。
本发明的效果可以通过以下实验仿真进一步说明:
1、实验条件与方法
硬件平台为:Intel Core2 Duo CPU E6550@2.33GHZ、2GBRAM;
软件平台为:MATLAB R2010a;
实验方法:分别为本发明和现有的基于SVM以及基于Wishart分布的极化SAR图像分类方法,其中现有的这两种方法都是极化SAR图像分类中引用较多的经典方法。
2、仿真内容与结果
本发明将图2所示的荷兰Flevoland地区的地物分布情况作为测试图像,该图像来源于AIRSAR,大小为380×420。
仿真一,用现有的基于SVM的分类方法对图2进行分类,分类结果见图3。
由图3可见,地物类别虽然大致可以得到有效划分,但由于受相干斑噪声的影响,各个区域中均出现不同程度的混杂现象,特别是图3中圈出的区域错分现象严重。
仿真二,用现有的基于Wishart分布的分类方法对图2进行分类,分类结果见图4。
由图4可见,分类效果与图3相比,类别区分更合理,特别是图4中黑色圈包含的区域分类更加清晰,但白色圈包含的区域仍存在较严重的错分现象。
仿真三,用本发明对图2进行分类,分类结果见图5。
由图5可见,与图3和图4相比,本发明得到的分类结果不仅保持了图4中黑色圈包含的区域分类清晰的优势,同时还大大减少了其他区域中错分的情况,特别是图5中白色圈包含的区域相比图3和图4中的相应区域来说,错分情况明显减少。
综上所述,本发明的方法在极化SAR图像的有监督分类问题上能取得更好的分类结果,它在一定程度上保证了均质区域的连贯性,使分类结果更好地吻合真实地物。
Claims (8)
1.一种基于SSAE和FSALS-SVM极化SAR图像分类方法,包括如下步骤:
步骤1,对待输入的所有极化SAR图像数据进行预处理,获得无标记训练样本和带标记训练样本;
步骤2,将每个稀疏自动编码器SAE均作为构成栈式稀疏自动编码器SSAE的一个基本单元,并将前一个SAE的输出作为后一个SAE的输入来逐个训练SAE,以获得每个SAE的模型参数,其中,步骤1中得到的所有无标记训练样本作为首个SAE的输入;
步骤3,将步骤2中获得的多个SAE级联在一起,构成SSAE,并将步骤1中得到的带标记训练样本输入该SSAE,进行前馈传导计算,获得这些带标记训练样本的深度特征;
步骤4,将深度特征进行归一化处理后,输入到快速稀疏逼近最小二乘支持向量机FSALS-SVM,通过迭代训练获得训练后的FSALS-SVM分类器;
步骤5,输入预处理后的待分类极化SAR图像,经过栈式稀疏自动编码器,获得数据的深度特征,并将其输入到训练后的FSALS-SVM分类器,获得图像的最终分类结果。
2.根据权利要求1所述的极化SAR图像分类方法,其中,对待输入的所有极化SAR图像数据进行预处理,获得无标记训练样本和带标记训练样本包括:
(1a)对待输入的所有极化SAR图像数据进行滑窗大小为7×7的Lee滤波;
(1b)Lee滤波后的极化SAR图像的任一像素点n均可表示为一个3×3协方差矩阵Mn:
<mrow>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>A</mi>
</mtd>
<mtd>
<mrow>
<mi>D</mi>
<mo>+</mo>
<mi>i</mi>
<mi>E</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>F</mi>
<mo>+</mo>
<mi>i</mi>
<mi>G</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>D</mi>
<mo>-</mo>
<mi>i</mi>
<mi>E</mi>
</mrow>
</mtd>
<mtd>
<mi>B</mi>
</mtd>
<mtd>
<mrow>
<mi>H</mi>
<mo>+</mo>
<mi>i</mi>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>F</mi>
<mo>-</mo>
<mi>i</mi>
<mi>G</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>H</mi>
<mo>-</mo>
<mi>i</mi>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mi>C</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,n∈[1,2,…,N],N为该极化SAR图像所包含的像素点个数,矩阵Mn中的大写字母A~I均为实数,可将这些字母组成列向量tn=[A B C D E F G H I]T,对每个像素点n所表示的列向量tn依次排放,构成整个待分类样本集;
(1c)对(1b)中得到的待分类样本集进行简单地处理使其利于分类,根据极化SAR图像数据本身的特点,直接对该待分类样本集统一乘以一个大的整数Q就获得好的分类结果;
(1d)将(1c)中简单地处理后的待分类样本集中的所有样本均作为无标记训练样本,并根据真实地物类别,随机选取无标记样本中的一部分进行标记,作为带标记训练样本。
3.根据权利要求2所述的极化SAR图像分类方法,其中,Q=500。
4.根据权利要求1所述的极化SAR图像分类方法,其中,将每个稀疏自动编码器(SAE)均作为构成栈式稀疏自动编码器(SSAE)的一个基本单元进行训练的过程:
(2a)将步骤1中得到的所有无标记训练样本作为首个SAE的输入x,分别设置该SAE输入单元、隐藏单元和输出单元的个数,该SAE的输入单元和输出单元的个数均为9,隐藏单元的个数为50;
(2b)随机初始化模型参数W1 (1)、W1 (2)、其中,W1 (1)表示首个SAE中输入单元和隐藏单元间的连接权重,W1 (2)表示首个SAE中隐藏单元和输出单元间的连接权重,表示首个SAE中隐藏单元的偏置项,表示首个SAE中输出单元的偏置项,计算隐藏单元和输出单元的激活值,并利用梯度下降法,最小化SSAE的整体代价函数Jsparse,从而获得训练后的模型参数
(2c)重新将步骤1中得到的所有无标记训练样本作为首个SAE的输入x,利用求得的模型参数可以计算出该SAE隐藏单元的激活值,并将该激活值作为第二个SAE的输入,设置该SAE输入单元、隐藏单元和输出单元的个数,该SAE的输入单元和输出单元的个数均为50,隐藏单元的个数为100;
(2d)随机初始化第二个SAE的模型参数其中,下标2代表第二个SAE,重复步骤(2b)进行训练,即可获得第二个SAE训练后的模型参数
(2e)令r表示第r个SAE,若r小于所需SAE的个数R,则随机初始化第r个SAE的模型参数Wr (1)、Wr (2)、将第r-1个SAE中隐藏单元的激活值作为第r个SAE的输入,并设置第r个SAE输入单元、隐藏单元和输出单元的个数,重复步骤(2b)进行训练,即可获得第r个SAE训练后的模型参数否则,停止训练,R=2。
5.根据权利要求4所述的极化SAR图像分类方法,其中,随机初始化模型参数,计算隐藏单元和输出单元的激活值,并利用梯度下降法,最小化栈式稀疏自动编码器的整体代价函数,从而获得训练后的模型参数包括:
(3a)随机初始化连接权重W1 (1)、W1 (2),使它们均为范围在内的数,偏置项为0,其中nin表示当前SAE中输入单元的个数,nout表示当前SAE中输出单元的个数;
(3b)利用如下前项传递函数F:
a(l+1)=F(W(l)a(l)+b(l))
进行前馈传导计算,分别得到隐藏单元和输出单元的激活值a(2)、a(3);其中,对任一SAE均有,l∈{1,2}即l可以取值1或2,a(1)表示输入单元的激活值,a(2)表示隐藏单元的激活值,a(3)表示输出单元的激活值;
(3c)定义SSAE的整体代价函数Jsparse:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>J</mi>
<mrow>
<mi>s</mi>
<mi>p</mi>
<mi>a</mi>
<mi>r</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<msup>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfrac>
<mi>&lambda;</mi>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>W</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mi>F</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&beta;</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>s</mi>
</munderover>
<mo>&lsqb;</mo>
<mi>&rho;</mi>
<mi>log</mi>
<mfrac>
<mi>&rho;</mi>
<msub>
<mover>
<mi>&rho;</mi>
<mo>^</mo>
</mover>
<mi>j</mi>
</msub>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<mi>log</mi>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&rho;</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mover>
<mi>&rho;</mi>
<mo>^</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,第一项为均方差项,第二项为权重衰减项,其目的在于减小权重幅度防止过度拟合,第三项 为稀疏性惩罚因子,用于将隐藏单元的平均活跃度保持在小的范围内,W表示输入单元和隐藏单元之间的连接权重W(1)与隐藏单元和输出单元之间的连接权重W(2)的转置的集合,即 表示矩阵W(2)的转置,m表示输入SAE的样本个数,和y(i)分别表示输入第i个样本时得到的实际输出结果和理想输出结果,的上标2表示平方,表示实际输出结果与理想输出结果之间误差的平方,的上标2表示平方,下标F表示F范数,则表示W的F范数的平方,λ为权重衰减参数,β为控制稀疏性惩罚因子的权重,s为隐藏单元的个数,ρ为稀疏性参数,表示第j个隐藏单元的平均活跃度,它可由下式计算得到:
<mrow>
<msub>
<mover>
<mi>&rho;</mi>
<mo>^</mo>
</mover>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>m</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mo>&lsqb;</mo>
<msubsup>
<mi>a</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
其中,表示在输入第i个样本x(i)的时隐藏单元j的激活度,上标2表示隐藏单元,利用反向传波算法,计算整体代价函数Jsparse的梯度;
(3d)通过梯度下降法求解代价函数Jsparse的最小化问题,从而获得训练后的模型参数
6.根据权利要求5所述的极化SAR图像分类方法,其中,参数ρ=0.1、β=3、λ=0.003。
7.根据权利要求1所述的极化SAR图像分类方法,其中,将深度特征进行归一化处理后,输入到快速稀疏逼近最小二乘支持向量机(FSALS-SVM),通过迭代训练获得训练后的FSALS-SVM分类器包括:
(4a)设置FSALS-SVM分类器中所用的核函数;
(4b)计算步骤3中获得的这些带标记训练样本的深度特征的均值mean和方差std,并将这些深度特征的取值范围规整化为零均值和单位方差,得到归一化的深度特征;
(4c)将第k类(k≥1)训练样本标记为1,其余类别均标记为-1,变为第k类的二分类问题;
(4d)输入归一化的深度特征,并利用(4c)中得到的第k类的二分类问题训练FSALS-SVM分类器,并在获得第k类参数后,令k的值加1;
(4e)若k大于训练样本类别数K,则停止训练,并将获得的第1类到第K类参数共K类参数依次排列后所构成的集合,作为训练后的FSALS-SVM分类器参数,获得训练后的FSALS-SVM分类器;否则,重复步骤(4c)~(4e),根据真实地物的类别总数确定K=9。
8.根据权利要求1所述的极化SAR图像分类方法,其中,将步骤1中得到的无标记训练样本输入到SSAE,获得深度特征,并将这些特征输入到训练后的FSALS-SVM分类器,获得图像的最终分类结果包括:
(5a)将步骤1中得到的无标记训练样本输入到SSAE,通过前馈传导计算,获得这些无标记训练样本的深度特征;
(5b)将(5a)中获得的深度特征的取值范围规整化为零均值和单位方差,得到归一化的深度特征;
利用步骤(4b)中得到的均值mean和方差std,将(5a)中获得的深度特征的取值范围规整化为零均值和单位方差,得到归一化的深度特征;
(5c)将归一化深度特征输入到训练后的FSALS-SVM分类器,获得最终分类结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410395368.2A CN104166859B (zh) | 2014-08-13 | 2014-08-13 | 基于ssae和fsals‑svm极化sar图像分类 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410395368.2A CN104166859B (zh) | 2014-08-13 | 2014-08-13 | 基于ssae和fsals‑svm极化sar图像分类 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104166859A CN104166859A (zh) | 2014-11-26 |
CN104166859B true CN104166859B (zh) | 2017-09-22 |
Family
ID=51910661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410395368.2A Active CN104166859B (zh) | 2014-08-13 | 2014-08-13 | 基于ssae和fsals‑svm极化sar图像分类 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104166859B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104463213A (zh) * | 2014-12-09 | 2015-03-25 | 西安电子科技大学 | 基于小波核稀疏lssvm的极化sar图像分类方法 |
CN104751173B (zh) * | 2015-03-12 | 2018-05-04 | 西安电子科技大学 | 基于协同表示和深度学习的极化sar图像分类方法 |
CN105069479A (zh) * | 2015-08-19 | 2015-11-18 | 西安电子科技大学 | 基于在线序列极限学习机的极化sar图像分类方法 |
CN105389505B (zh) * | 2015-10-19 | 2018-06-12 | 西安电子科技大学 | 基于栈式稀疏自编码器的托攻击检测方法 |
CN105678340B (zh) * | 2016-01-20 | 2018-12-25 | 福州大学 | 一种基于增强型栈式自动编码器的自动图像标注方法 |
CN105825223A (zh) * | 2016-03-09 | 2016-08-03 | 西安电子科技大学 | 基于深度学习和度量学习的极化sar地物分类方法 |
CN106096651B (zh) * | 2016-06-12 | 2019-04-23 | 西安电子科技大学 | 基于超像素和度量学习的极化sar地物分类方法 |
CN106124212B (zh) * | 2016-06-16 | 2018-09-04 | 燕山大学 | 基于稀疏编码器和支持向量机的滚动轴承故障诊断方法 |
CN106251375B (zh) * | 2016-08-03 | 2020-04-07 | 广东技术师范学院 | 一种通用隐写分析的深度学习堆栈式自动编码方法 |
CN106323636A (zh) * | 2016-08-16 | 2017-01-11 | 重庆交通大学 | 栈式稀疏自动编码深度神经网络的机械故障程度特征自适应提取与诊断方法 |
CN106485353B (zh) * | 2016-09-30 | 2019-11-29 | 中国科学院遥感与数字地球研究所 | 空气污染物浓度预报方法及系统 |
CN107064894A (zh) * | 2017-01-11 | 2017-08-18 | 南京御达电信息技术有限公司 | 一种基于深度学习的杂波抑制方法 |
CN106897987A (zh) * | 2017-01-18 | 2017-06-27 | 江南大学 | 基于平移不变剪切波和栈式自编码的图像融合方法 |
CN106920243B (zh) * | 2017-03-09 | 2019-12-17 | 桂林电子科技大学 | 改进的全卷积神经网络的陶瓷材质件序列图像分割方法 |
CN107038421A (zh) * | 2017-04-17 | 2017-08-11 | 杭州电子科技大学 | 基于稀疏堆栈自编码的调制样式识别方法 |
CN107292885A (zh) * | 2017-08-08 | 2017-10-24 | 广东工业大学 | 一种基于自动编码器的产品缺陷分类识别方法及装置 |
CN109816002B (zh) * | 2019-01-11 | 2022-09-06 | 广东工业大学 | 基于特征自迁移的单一稀疏自编码器弱小目标检测方法 |
CN109949345A (zh) * | 2019-02-21 | 2019-06-28 | 西北大学 | 基于神经网络的获取非均衡sar图像变化检测图的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236063A (zh) * | 2013-05-03 | 2013-08-07 | 河海大学 | 基于多尺度谱聚类及决策级融合的sar图像溢油检测方法 |
CN103810699A (zh) * | 2013-12-24 | 2014-05-21 | 西安电子科技大学 | 基于无监督深度神经网络的sar图像变化检测方法 |
CN103955702A (zh) * | 2014-04-18 | 2014-07-30 | 西安电子科技大学 | 基于深度rbf网络的sar图像地物分类方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7171042B2 (en) * | 2000-12-04 | 2007-01-30 | Intel Corporation | System and method for classification of images and videos |
-
2014
- 2014-08-13 CN CN201410395368.2A patent/CN104166859B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103236063A (zh) * | 2013-05-03 | 2013-08-07 | 河海大学 | 基于多尺度谱聚类及决策级融合的sar图像溢油检测方法 |
CN103810699A (zh) * | 2013-12-24 | 2014-05-21 | 西安电子科技大学 | 基于无监督深度神经网络的sar图像变化检测方法 |
CN103955702A (zh) * | 2014-04-18 | 2014-07-30 | 西安电子科技大学 | 基于深度rbf网络的sar图像地物分类方法 |
Non-Patent Citations (1)
Title |
---|
基于多层编码器的SAR 目标及阴影联合特征提取算法;孙志军等;《雷达学报》;20130630;第2卷(第2期);第195-202页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104166859A (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104166859B (zh) | 基于ssae和fsals‑svm极化sar图像分类 | |
CN108038445B (zh) | 一种基于多视角深度学习框架的sar自动目标识别方法 | |
Zou et al. | Deep learning based feature selection for remote sensing scene classification | |
CN108416755A (zh) | 一种基于深度学习的图像去噪方法及系统 | |
CN110533606B (zh) | 基于生成式对抗网络的安检x光违禁品图像数据增强方法 | |
CN110210486A (zh) | 一种基于素描标注信息的生成对抗迁移学习方法 | |
CN103927551B (zh) | 基于超像素相关矩阵的极化sar图像半监督分类方法 | |
CN106250931A (zh) | 一种基于随机卷积神经网络的高分辨率图像场景分类方法 | |
CN104217438A (zh) | 基于半监督的图像显著性检测方法 | |
CN104156728B (zh) | 一种基于栈式编码和softmax的极化SAR图像分类方法 | |
CN106067042A (zh) | 基于半监督深度稀疏滤波网络的极化sar分类方法 | |
CN106651887A (zh) | 一种基于卷积神经网络的图像像素分类方法 | |
CN107292336A (zh) | 一种基于dcgan的极化sar图像分类方法 | |
CN105046268B (zh) | 基于Wishart深度网络的极化SAR图像分类方法 | |
CN107992891A (zh) | 基于光谱矢量分析多光谱遥感图像变化检测方法 | |
CN107273853A (zh) | 一种基于类心和协方差对齐的遥感图像迁移学习方法 | |
CN106683102A (zh) | 基于脊波滤波器和卷积结构模型的sar图像分割方法 | |
CN104573742A (zh) | 医学图像分类方法和系统 | |
CN106326938B (zh) | 基于弱监督学习的sar图像目标鉴别方法 | |
CN109948692A (zh) | 基于多颜色空间卷积神经网络及随机森林的计算机生成图片检测方法 | |
CN106611423A (zh) | 基于脊波滤波器和反卷积结构模型的sar图像分割方法 | |
US20210312243A1 (en) | Method for synthesizing image based on conditional generative adversarial network and related device | |
CN104992183A (zh) | 自然场景中的显著目标的自动检测方法 | |
CN106096660A (zh) | 基于独立成分分析算法的卷积神经网络 | |
CN109300128A (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 |