基于脊波和深度卷积网络的高光谱图像分类方法
【技术领域】
[0001] 本发明属于图像处理技术领域，特别涉及一种高光谱图像分类方法，可用于气象 监测、环境监测、土地利用、城市规划以及防灾减灾。
【背景技术】
[0002] 高光谱分辨率遥感是指利用很多很窄的电磁波波段从感兴趣的物体中获取有关 数据。它的最大特点就是在获得目标图像二维空间景象信息的同时，还可以获得高分辨率 的一维表征其物理属性的光谱信息，即具有"图谱合一"的特性，是代表遥感最新成就的新 型技术之一。高光谱遥感与常规遥感数据的主要区别在于高光谱遥感是窄波段成像，可以 获得连续的光谱信息，探测到常规遥感中不可探测的物质。因此，当一个宽波段的系统只能 大致区分不同物质种类时，高光谱传感器却可以为物质的详细鉴定和更准确估计它的丰富 程度提供了潜在的可能。
[0003] 高光谱数据处理的一个主要内容就是地物目标分类。分类是一种描述地物目标或 种类的分析技术，其主要任务是对数据体的每个像素点赋予一个类别标记以产生专题地图 的一种过程，它是人们从遥感影像上提取有用信息的重要途径之一。分类后产生的专题地 图可以清晰地反映出地物的空间分布，便于人们从中认识和发现其规律，使高光谱遥感图 像具有真正的使用价值并有效的投入到实际应用中。
[0004] 传统的图像分类方法是目视解译，目视技术利用了人类优秀的思维能力来定性评 价图像中的空间模式。这种方法存在一定的缺点，它需要图像目视判读者具有丰富的地学 知识和目视判读经验，并且劳动强度大，需要花费大量的时间。另外光谱特性不是都可以用 目视解译的方法全面评定的。为了提高分类的质量和效率，从上世纪70年代起，人们开始 重视由计算机自动获取遥感图像中专题信息的方法研究。当时主要是利用传统的统计模式 识别方法进行遥感计算机解译，分类精度不能够令人满意，随着遥感影像不断的发展变化， 对分类算法也在不断的提出新的要求，因此改进现有分类算法，寻找新的方法一直是遥感 应用研究中的热点之一。
[0005] 高光谱图像的特征提取是高光谱图像分类的一个重要部分，对分类精度有很大的 影响。目前，市场上使用的高光谱图像特征提取方法主要有基于光谱信息的特征提取方法， 基于空间信息的特征提取方法，和结合空间信息和光谱信息的特征提取方法。
[0006] 在基于光谱信息的特征提取方法中，高光谱数据中的每个像素在光谱空间中表现 为一条波谱响应曲线。不同的地物具有不同的波谱反射率和吸收特性；相同的地物在不 同的波段也具有不同的反射波谱率，在遥感数据中表现为不同的辐射强度。因此不同的地 物的波谱曲线形态不同；同一地物的波谱曲线也并不是平坦的而是起伏变化的，常常有多 个峰点与谷点。基于光谱信息的分类方法就是利用不同地物的波谱曲线对地物进行分类， 常用的特征提取算法有稀疏PCA(PrincipalComponentsAnalysis)，ICA(Independent ComponentAnalysis)以及LDA(LatentDirichletAllocation)等。但是地物光谱响应会 受到很多因素的影响，如太阳照度、大气透明度和风速等，而这些因素通常都是很难准确测 量的，所以实际测量得到的地物光谱响应曲线与实际的曲线可能有很大的差异。这种情况 下，光谱空间描述方式不能很好的适应高光谱数据的分析，所以基于光谱信息的特征提取 方法常常得到的分类精度不够理想。
[0007] 基于空间信息的特征提取方法是只利用高光谱图像的空间信息进行分类，典型的 方法如基于方差的特征提取方法、基于灰度共生矩阵的特征提取方法、以及小波分析的特 征提取方法。该类方法是一种人工经验的特征提取方法，需要预先知道图像的特征，再对应 的选择合适方法，所以该类方法需要较好的先验知识才能达到较好的分类效果。
[0008] 为此，很多研究者提出了结合空间和光谱信息的特征提取方法，借助于高光谱图 像的光谱和空间信息来提高分类精度，如IFRF(ImageFusionandRecursiveFiltering)、 EPF(Edge-PreservingFilters)和NMFL(NonlinearandLinearMultipleFeature Learning)方法。该类方法虽然在一定程度上克服了仅仅使用光谱信息或空间信息造成的 地物错分问题问题，但是仍然需要较多的先验知识才能够获得较好的分类精度。
[0009] 神经网络是一类有效的提取空谱特征的方法，也是一种主动的特征学习的方法， 不需要对图像有先验知识，典型的神经网络如BP神经网络、小波神经网络、和脊波神经网 络。但是这些都是浅层的神经网络，都只包含3层，为了能够更好的挖掘图像更深层的特 征，深度神经网络的模型被提出，典型的深度神经网络有自编码深度网络、受限玻尔兹曼机 深度网络、深度卷积网络等。由于深度卷积网络是一个真正的二维神经网络，对于二维的图 像来说，深度卷积网络能够更好的表示图像的特征。然而传统的深度卷积网络滤波器的初 始化都是采用随机初始化，或者高斯函数初始化，好的初始化对于网络的性能和逼近速率 有很大的影响，而这些传统的初始化方法很难达到一个理想的效果。
【发明内容】
[0010] 本发明的目的在于针对上述已有技术的不足，提出一种基于脊波和深度卷积网络 的高光谱图像分类方法，以解决现有技术在高光谱图像分类问题中难以学习到有效的分类 特征，以及传统的深度卷积网络难以达到更小的计算复杂度的问题，提高光谱图像分类的 精度和速度。
[0011] 为实现上述目的，本发明的实现步骤如下：
[0012] 1)输入一幅高光谱图像和该图像中地物信息的类标，从该高光谱图像中选择 10%的像素作为训练样本；
[0013] 2)提取训练样本的光谱信息：沿着高光谱图像光谱的维度，提取每一个训练样本 的光谱信息，组成光谱向量f]，j= 1，...，J，J是训练样本的个数；
[0014] 3)对高光谱图像进行降维处理，保留前4个主分量，组成降维后的图像；
[0015] 4)提取训练样本的空间信息：以每个训练样本为中心，在降维后的图像中每一维 上，选取7X7大小的窗口，得到该样本在该维上的空间信息g/./=l4 ；
[0016] 5)将每一个训练样本的空间信息gf和光谱向量组成一个正方形的训练样本图 像块，并对该图像进行归一化，得到归一化后的训练样本正方形图像块F];
[0017] 6)构造一个5层的深度卷积网络，并用归一化后的训练样本正方形图像块FM乍为 该卷积网络的输入，对该网络进行训练，得到训练好的网络；
[0018] 7)将高光谱图像中剩余90%的像素作为测试样本，提取每一个样本的光谱向量 r，以及空间向量n1，i= 1，. . .，4,组成一个正方形的测试样本图像块，并对该图像块进行 归一化，得到归一化后的测试样本正方形图像块P，q= 1，...，Q，Q是测试样本的个数；
[0019] 8)将归一化后的测试样本正方形图像块P输入到步骤6)训练好的网络中，根据 网络中训练好的参数值，得到该样本的类标值，完成分类。
[0020] 本发明与现有技术相比，具有如下效果：
[0021] (a)本发明使用了空间信息和光谱信息结合的高光谱图像分类方法，克服了传统 的高光谱图像分类方法中仅仅采用光谱信息，忽略了有效利用空间信息的问题，提高了分 类精度。
[0022] (b)本发明在传统深度卷积神经网络的框架下，对网络中卷积层滤波器的初始化 采用脊波函数，克服了传统滤波器初始化方法难以有效、快速的逼近高光谱图像的空间和 光谱信息的问题，提高了分类的速度。
【附图说明】
[0023]图1是本发明的基于脊波和深度卷积网络的高光谱图像分类方法流程框图；
[0024] 图2是本发明实验使用图像及其真实地物分类图；
[0025] 图3是用本发明和现有技术对图2的分类结果对比图；
[0026] 图4是用本发明和现有方法对图2的误差下降对比图。
【具体实施方式】
[0027] 以下参照附图，对本发明的技术方案和效果做进一步详细描述。
[0028]参照图1，本发明的实现步骤如下：
[0029] 步骤1，输入图像。
[0030] 输入一幅高光谱图像，如图所示，其中2(a)为输入的高光谱图像，图2(b)为2(a) 对应的类标图像，从2(a)中选择10%的像素作为训练样本。
[0031] 步骤2,提取训练样本的光谱信息。
[0032] 设步骤1中输入的高光谱图像的光谱维数为V，对每一个训练样本，提取该样本每 一维的光谱值，组成光谱向量f°,j= 1，. . .，J，J是训练样本的个数，光谱向量的维数为 V。
[0033] 步骤3,对高光谱图像降维。
[0034] 对图像进行降维的方法有稀疏PCA，LDA，PCA，ICA等方法，本发明使用PCA方法对 步骤1中输入的高光谱图像进行降维处理，具体步骤如下：
[0035] 3a)求出步骤1中输入的高光谱图像的协方差矩阵；
[0036] 3b)根据协方差矩阵得到步骤1中输入的高光谱图像的分量图；
[0037] 3c)最后保留前4个能量值最大的分量图，将这4个分量图组成降维后的图像。
[0038] 步骤4，提取训练样本的空间信息。
[0039] 在步骤2中降维后的图像中的每一个分量图上，以每个训练样本为中心，选取 7X7大小的窗口，得到该样本在分量图上的空间信息g/V= 1，.…4。
[0040] 步骤5,将空间信息与光谱信息相结合。
[0041] 5a)将4个7X7大小的代表图像空间信息的图像块拼成一个大小为
[0042] 5b)将光谱向量f]重新排列成一个大小为H*1 2的光谱信息长方形图像块F，H为 长方形图像块的长，12为长方形图像块的宽，H*12=V，其中V为光谱向量的维数；
[0043] 5c)在空间信息正方形图像块G中随机选取一个大小为（H-14)*(H-12)的图像块 B，根据该图像块B和空间信息正方形图像块G及光谱信息长方形图像块F，构建大小为H*H
[0044] 5d)对得到的训练样本正方形图像块A进行归一化，得到归一化的训练样本正方 形图像块F]，j= 1，. . .，J，J是训练样本的个数。
[0045] 步骤6,构造5层深度卷积网络并对其进行训练。
[0046] 6a)构造5层深度卷积网络：其中第1层为输入层，第2层和第3层为卷积层，第 4层为全链接层，第5层为softmax分类器；该第1层的输入为归一化的训练样本正方形图 像块F];该第2层包含L