申请/专利权人：大连理工大学;国网吉林省电力有限公司检修公司;国网吉林省电力有限公司;国家电网公司

申请日：2018-08-27

公开（公告）日：2021-05-07

公开（公告）号：CN109145993B

主分类号：G06K9/62(20060101)

分类号：G06K9/62(20060101)

优先权：

专利状态码：有效-授权

法律状态：2021.05.07#授权;2019.01.29#实质审查的生效;2019.01.04#公开

摘要：一种基于多特征与非负自动编码器的SAR图像分类方法，属于图像处理技术领域。基于灰度梯度共生矩阵提取SAR图像的图像块空间域特征；基于二维Gabor变换提取SAR图像的图像块变换域特征；将图像块的空间域特征与变换域特征组合；选取SAR图像块的训练样本集和测试样本集；利用训练样本集对多层非负自动编码器与softmax分类器进行训练；采用训练好的非负自动编码器网络进行分类；得到分类结果图。本发明结合了SAR图像的空间信息与变换域信息，获得了SAR图像的多维特征，并利用非负自动编码器对特征进行优化，提升了特征的区分性，进而有效提高了分类的准确度，可用于高分辨率SAR图像地物分类与目标识别等。

主权项：1.一种基于多特征与非负自动编码器的SAR图像分类方法，步骤如下：1基于灰度梯度共生矩阵的SAR图像空间域特征提取：1a输入一幅a×b的SAR图像，按照w×w的窗大小进行分块，得到图像块；1b基于灰度梯度共生矩阵，提取每个图像块的空间域特征，计算公式为： 其中，Hij表示SAR图像块的灰度值为i且对应梯度图的梯度值为j的像素个数，表示归一化的Hij，Nh和Nt分别表示灰度级与梯度级；1c将每个图像块的空间域特征向量化，具体如下：T＝[T1,T2,T3,T4,…,T15]；2基于二维Gabor变换的SAR图像变换域特征提取：2a对分块后得到的SAR图像的图像块，提取基于二维Gabor变换的变换域特征，计算公式为： 其中，Ix,y表示SAR图像，x和y表示SAR图像像素点的位置，x0＝xcosθ+ysinθ，y0＝-xsinθ+ycosθ，w0表示Gabor滤波器的中心频率，取c个不同尺度，θ表示Gabor滤波器的方向，取d个方向，分别为σ表示高斯函数的标准差，表示卷积运算；提取c×d维的Gabor特征，按照w×w的窗大小进行分块，并计算块的均值特征来代表图像块的c×d维空间域特征G；2b对每个图像块的变换域特征向量化，具体如下：G′＝[G1,G2,G3,G4,…,Gc×d]；3将空间域特征与变换域特征进行组合：将每个图像块的空间域特征向量和变换域特征向量组合成一个特征向量，具体如下：f＝[T,G′]；4选取SAR图像块的训练样本集和测试样本集：从SAR图像所有带标签的图像块中，随机选取p％的图像块形成训练样本集，其余的图像块作为测试样本集；5利用训练样本集对多层非负自动编码器与softmax分类器进行训练：5a先初始化三层非负自动编码器的权重和偏置参数，再利用训练样本集逐层预训练各层网络；预训练第k层非负自动编码器的权重和偏置参数的目标方程为： 其中，上式第一项为解码输出与编码前输入的平均重建误差，第二项为负权重的惩罚项，表示第k层网络的解码输出，表示第k层网络的输入同时也表示第k-1层网络的隐含输出，N表示训练样本数量，λ表示负权重惩罚因子；利用反向传播算法求解上式；5b预训练完各层非负自动编码器的权重和偏置参数后，第三层非负自动编码器的隐含输出作为优化后的样本特征；5c将优化后的训练样本特征与标签输入到softmax分类器，训练得到softmax分类器的权重和偏置参数；5d从softmax分类器、第三层非负自动编码器、第二层非负自动编码器、第一层非负自动编码器的自上而下的顺序，对整个网络的参数进行微调；反向微调整个网络权重和偏置参数的目标方程为： 其中，上式第一项为所有编码层总的重建误差，第二项为各个编码层负权重的惩罚项，K表示网络层数，表示第k层网络的隐含输出；利用反向传播算法求解上式，得到微调的各层非负自动编码器与softmax分类器的权重和偏置参数；6对测试样本集进行分类：6a将测试样本xi输入到三层非负自动编码器中，第三层非负自动编码器的隐含输出作为测试样本优化后的样本特征；6b将优化后的测试样本特征输入到softmax分类器，计算其属于第c类的概率，具体如下： 其中，W′j和b′j对应第j类的部分权重和偏置，W′c和b′c对应第c类的部分权重和偏置，C表示总的类别数；最后，测试样本xi的预测标签如下：labelxi＝argmaxPcxi7输出分类结果图：根据预测标签向量和测试样本的空间位置，绘制最后分类结果图。

全文数据：基于多特征与非负自动编码器的SAR图像分类方法技术领域本发明属于图像处理技术领域，特别是涉及地物分类领域中基于多特征与非负自动编码器的SAR图像分类方法，可用于SAR图像地物分类与目标识别。背景技术合成孔径雷达Syntheticapertureradar,SAR是一种主动成像传感器，具有全天候、全天时的数据获取能力，比传统的光学遥感技术具有明显优势。随着遥感技术的不断发展，SAR系统获取的图像分辨率越来越高，高分辨率SAR图像能够反映更加详细的地物信息，满足了许多实际应用的需求。SAR图像分类是SAR图像解译的重要内容，在军事侦察、资源探测、地理测绘等领域有着广泛的应用。由于SAR图像受到相干斑噪声的严重污染，导致缺乏有效的特征提取模型。一些基本的SAR图像特征被用于分类，如幅度、相位、密度等，但往往都不能得到很好的分类效果。SAR图像分类是分析获取的地物后向散射信号等特征进行地物分类的过程。目前的SAR图像分类按照是否需要人工的类别标签，可以分为有监督分类和无监督分类。在无监督分类方法中，不需要样本的类别标签，因此不需要人工参与，主要基于后向散射等信号实现分类。传统的无监督分类方法，如模糊C均值聚类、分水岭算法、迭代自组织分析算法，被用于SAR图像分类。在实际应用中，往往获得标签样本比较困难，这样无监督分类表现出一定优势。然而，由于SAR图像存在明显的相干斑噪声，无监督分类的准确度通常不是很理想。在有监督分类方法中，需要样本的类别标签，因此需要人工参与获得标签样本来训练分类器。不同的机器学习方法，如贝叶斯分类器、马尔可夫模型、支持向量机等，被用于SAR图像分类。例如，C.Tison等人2004年在IEEETransactionsonGeoscienceandRemoteSensing的第42卷第10期上发表的《AnewstatisticalmodelforMarkovianclassificationofurbanareasinhigh-resolutionSARimages》，提出了基于Fisher分布和对数矩估计的数学模型，并采用马尔可夫分割进行分类，获得了准确的高分辨率城市地区分类结果。A.Voisin等人2013年在IEEEGeoscienceandRemoteSensingLetters的第10卷第1期上发表的《ClassificationofveryhighresolutionSARimagesofurbanareasusingcopulasandtextureinahierarchicalmarkovrandomfieldmodel》，提出了基于上下文分层的监督贝叶斯分类方法，结合SAR图像振幅和结构信息建立统计模型，解决了城市区域的高分辨率SAR图像分类问题。近年来，深度学习方法在各种机器学习任务中取得了突破性的提升，获得了广泛应用，也逐渐应用于遥感图像分类任务中。例如，Z.Zhao等人2017年在PatternRecognition的第61卷上发表的《Discriminantdeepbeliefnetworkforhigh-resolutionSARimageclassification》，将集成学习与深度置信网络结合，学习高分辨率SAR图像的高层次特征，取得优异的分类效果。上述方法未考虑SAR图像相干斑噪声影响且未充分挖掘SAR图像的不同特征，未有效利用深度网络提高特征的区分性，从而影响了分类效果。发明内容本发明旨在克服上述现有技术的问题，提出一种基于多特征与非负自动编码器的SAR图像分类方法，进而提高SAR图像分类的准确度。本发明的技术方案：基于多特征与非负自动编码器的SAR图像分类方法，步骤如下：1基于灰度梯度共生矩阵的SAR图像空间域特征提取：1a输入一幅a×b的SAR图像，按照w×w的窗大小进行分块，得到图像块；1b基于灰度梯度共生矩阵，提取每个图像块的空间域特征；1c将每个图像块的空间域特征向量化；2基于二维Gabor变换的SAR图像变换域特征提取：2a对分块后得到的SAR图像的图像块，提取基于二维Gabor变换的变换域特征；2b对每个图像块的变换域特征向量化；3将空间域特征与变换域特征进行组合：将每个图像块的空间域特征向量和变换域特征向量组合成一个特征向量；4选取SAR图像块的训练样本集和测试样本集：从SAR图像所有带标签的图像块中，随机选取p％的图像块形成训练样本集，其余的图像块作为测试样本集；取每个图像块像素标签最多的类作为每个图像块的标签；5利用训练样本集对多层非负自动编码器与softmax分类器进行训练：5a先初始化三层非负自动编码器的权重和偏置参数，再利用训练样本集逐层预训练各层网络；5b预训练完各层非负自动编码器的权重和偏置参数后，第三层非负自动编码器的隐含输出作为优化后的样本特征；5c将优化后的训练样本特征与标签输入到softmax分类器，训练得到softmax分类器的权重和偏置参数；5d从softmax分类器、第三层非负自动编码器、第二层非负自动编码器、第一层非负自动编码器的自上而下的顺序，对整个网络的参数进行微调。6对测试样本集进行分类：6a将测试样本集输入到三层非负自动编码器中，第三层非负自动编码器的隐含输出作为测试样本优化后的样本特征；6b将优化后的测试样本特征输入到softmax分类器，得到预测标签；7输出分类结果图：根据预测标签向量和测试样本的空间位置，绘制最后的分类结果图。本发明与现有技术相比，主要具有如下的优点：第一，本发明提出了基于非负自动编码器的深度网络，不仅利用了深度网络优异的特征表示能力，获得高层次的特征；而且加入非负权重约束，提高编码表示的重建质量和稀疏性，能够提高特征的有效性。第二，本发明结合了基于灰度梯度共生矩阵的空间域特征和基于二维Gabor变换的变换域特征，获得了丰富的SAR图像块特征。采用图像块特征的提取方法，能够有效克服高分辨率SAR图像严重的相干斑噪声影响的问题，对提高分类精度起很大作用。附图说明图1为SAR图像分类的实现流程图；图2为本发明适用的SAR图像；图3为本发明适用的SAR图像的真实地物标记图；图4为利用现有自动编码器网络对SAR图像进行分类的分类结果图；图5为利用现有稀疏自动编码器网络对SAR图像进行分类的分类结果图；图6为利用本发明方法对SAR图像进行分类的分类结果图。具体实施方式下面结合具体实例和附图，对本发明作详细的阐述。根据图1，基于多特征与非负自动编码器的SAR图像分类方法，包括如下步骤：1基于灰度梯度共生矩阵的SAR图像空间域特征提取：1a输入一幅3580×2250的SAR图像，按照5×5的窗大小进行分块，得到322200个图像块；1b基于灰度梯度共生矩阵，提取每个图像块的15维的空间域特征，计算公式为：其中，Hij表示SAR图像块的灰度值为i且对应梯度图的梯度值为j的像素个数，表示归一化的Hij，Nh和Nt分别表示灰度级与梯度级；1c将每个图像块的空间域特征向量化，具体如下：T＝[T1,T2,T3,T4,…,T15]；2基于二维Gabor变换的SAR图像变换域特征提取：2a对分块后得到的SAR图像的图像块，提取基于二维Gabor变换的变换域特征，计算公式为：其中，Ix,y表示SAR图像，x和y表示图像像素点的位置，x0＝xcosθ+ysinθ，y0＝-xsinθ+ycosθ，w0表示Gabor滤波器的中心频率，取5个不同尺度，θ表示Gabor滤波器的方向，取8个方向σ表示高斯函数的标准差，表示卷积运算。这样提取了40维的Gabor特征，按照5×5的窗大小进行分块，并计算块的均值特征来代表图像块的40维空间域特征G。2b对每个图像块的变换域特征向量化，具体如下：G′＝[G1,G2,G3,G4,…,G40]；3将空间域特征与变换域特征进行组合：将每个图像块的空间域特征向量和变换域特征向量组合成一个特征向量，具体如下：f＝[T,G′]；4选取SAR图像块的训练样本集和测试样本集：从SAR图像所有带标签的图像块中，随机选取50％的图像块形成训练样本集，其余的图像块作为测试样本集；5利用训练样本集对多层非负自动编码器与softmax分类器进行训练：5a先初始化三层非负自动编码器的权重和偏置参数，再利用训练样本集逐层预训练各层网络。预训练第k层非负自动编码器的权重和偏置参数的目标方程为：其中，上式第一项为解码输出与编码前输入的平均重建误差，第二项为负权重的惩罚项，表示第k层网络的解码输出，表示第k层网络的输入同时也表示第k-1层网络的隐含输出，N表示训练样本数量，λ表示负权重惩罚因子。利用反向传播算法求解上式。5b预训练完各层非负自动编码器的权重和偏置参数后，第三层非负自动编码器的隐含输出作为优化后的样本特征；5c将优化后的训练样本特征与标签输入到softmax分类器，训练得到softmax分类器的权重和偏置参数；5d从softmax分类器、第三层非负自动编码器、第二层非负自动编码器、第一层非负自动编码器的自上而下的顺序，对整个网络的参数进行微调。反向微调整个网络权重和偏置参数的目标方程为：其中，上式第一项为所有编码层总的重建误差，第二项为各个编码层负权重的惩罚项，K表示网络层数，表示第k层网络的隐含输出。同样，利用反向传播算法求解上式，得到微调的各层非负自动编码器与softmax分类器的权重和偏置参数；6对测试样本集进行分类：6a将测试样本xi输入到三层非负自动编码器中，第三层非负自动编码器的隐含输出作为测试样本优化后的样本特征；6b将优化后的测试样本特征输入到softmax分类器，计算其属于第c类的概率，具体如下：其中，Wj′和bj′对应第j类的部分权重和偏置，Wc′和bc′对应第c类的部分权重和偏置，C表示总的类别数。最后，测试样本xi的预测标签如下：labelxi＝argmaxPcxi7输出分类结果图：把各个图像块按照空间位置，绘制最后3580×2250的分类结果图，其中每个图像块中各个像素的类别均等于该图像块分类后的类别。以下通过仿真实验，对本发明的技术效果进行说明：1、仿真条件与内容本发明的实验数据为2013年获取的挪威Lillestroem市的TerraSAR-X数据，如图2所示，图像大小为3580×2250，图3为该图像的真实地物标记图，共有5类不同地物，图4为利用现有自动编码器网络对SAR图像进行分类的分类结果图，图5为利用现有稀疏自动编码器网络对SAR图像进行分类的分类结果图，图6为利用本发明方法对SAR图像进行分类的分类结果图，表一为这三种技术的分类精度对比。仿真实验中，本发明和对比方法都是在MatlabR2017a中编程实现。2、仿真结果分析表一分类精度对比类别自动编码器网络稀疏自动编码器网络本发明河流84.9088.2491.69森林89.9092.0796.08草地85.6086.7190.01建筑85.6591.1297.15公路81.8580.4986.04总体精度86.7688.7892.49平均精度85.5887.7292.19Kappa系数0.80870.83810.8919由表一可知，本发明的方法比现有自动编码器网络、稀疏自动编码器网络都获得了更高的分类精度，证明本发明对高分辨率SAR图像分类的优异效果。从图4、图5和图6可知，本发明方法的分类结果图的误分类杂点更少，说明比现有自动编码器网络、稀疏自动编码器网络具有更好的相干斑噪声克服能力。总之，本发明的方法能够有效改善高分辨率SAR图像分类效果。

权利要求：1.一种基于多特征与非负自动编码器的SAR图像分类方法，步骤如下：1基于灰度梯度共生矩阵的SAR图像空间域特征提取：1a输入一幅a×b的SAR图像，按照w×w的窗大小进行分块，得到图像块；1b基于灰度梯度共生矩阵，提取每个图像块的空间域特征，计算公式为：其中，Hij表示SAR图像块的灰度值为i且对应梯度图的梯度值为j的像素个数，表示归一化的Hij，Nh和Nt分别表示灰度级与梯度级；1c将每个图像块的空间域特征向量化，具体如下：T＝[T1,T2,T3,T4,…,T15]；2基于二维Gabor变换的SAR图像变换域特征提取：2a对分块后得到的SAR图像的图像块，提取基于二维Gabor变换的变换域特征，计算公式为：其中，Ix,y表示SAR图像，x和y表示SAR图像像素点的位置，x0＝xcosθ+ysinθ，y0＝-xsinθ+ycosθ，w0表示Gabor滤波器的中心频率，取c个不同尺度，θ表示Gabor滤波器的方向，取d个方向，分别为σ表示高斯函数的标准差，表示卷积运算；提取c×d维的Gabor特征，按照w×w的窗大小进行分块，并计算块的均值特征来代表图像块的c×d维空间域特征G；2b对每个图像块的变换域特征向量化，具体如下：G′＝[G1,G2,G3,G4,…,Gc×d]；3将空间域特征与变换域特征进行组合：将每个图像块的空间域特征向量和变换域特征向量组合成一个特征向量，具体如下：f＝[T,G′]；4选取SAR图像块的训练样本集和测试样本集：从SAR图像所有带标签的图像块中，随机选取p％的图像块形成训练样本集，其余的图像块作为测试样本集；5利用训练样本集对多层非负自动编码器与softmax分类器进行训练：5a先初始化三层非负自动编码器的权重和偏置参数，再利用训练样本集逐层预训练各层网络；预训练第k层非负自动编码器的权重和偏置参数的目标方程为：其中，上式第一项为解码输出与编码前输入的平均重建误差，第二项为负权重的惩罚项，表示第k层网络的解码输出，表示第k层网络的输入同时也表示第k-1层网络的隐含输出，N表示训练样本数量，λ表示负权重惩罚因子；利用反向传播算法求解上式；5b预训练完各层非负自动编码器的权重和偏置参数后，第三层非负自动编码器的隐含输出作为优化后的样本特征；5c将优化后的训练样本特征与标签输入到softmax分类器，训练得到softmax分类器的权重和偏置参数；5d从softmax分类器、第三层非负自动编码器、第二层非负自动编码器、第一层非负自动编码器的自上而下的顺序，对整个网络的参数进行微调；反向微调整个网络权重和偏置参数的目标方程为：其中，上式第一项为所有编码层总的重建误差，第二项为各个编码层负权重的惩罚项，K表示网络层数，表示第k层网络的隐含输出；利用反向传播算法求解上式，得到微调的各层非负自动编码器与softmax分类器的权重和偏置参数；6对测试样本集进行分类：6a将测试样本xi输入到三层非负自动编码器中，第三层非负自动编码器的隐含输出作为测试样本优化后的样本特征；6b将优化后的测试样本特征输入到softmax分类器，计算其属于第c类的概率，具体如下：其中，W′j和b′j对应第j类的部分权重和偏置，W′c和b′c对应第c类的部分权重和偏置，C表示总的类别数；最后，测试样本xi的预测标签如下：labelxi＝argmaxPcxi7输出分类结果图：根据预测标签向量和测试样本的空间位置，绘制最后分类结果图。

百度查询： 大连理工大学;国网吉林省电力有限公司检修公司;国网吉林省电力有限公司;国家电网公司 基于多特征与非负自动编码器的SAR图像分类方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD