基于区域光谱梯度特征对比的高光谱图像显著性检测方法
【技术领域】
[0001] 本发明涉及一种高光谱图像显著性检测方法，特别涉及一种基于区域光谱梯度特 征对比的高光谱图像显著性检测方法。
【背景技术】
[0002] 高光谱图像以精细的光谱分辨率采集了场景中的物体在不同波长上的反射值，波 长范围为可见光波段的高光谱图像可采集能为人的视觉系统感知但普通真彩图像(RGB三 色图像)却无法表达的信息。因此，利用高光谱图像来解决传统机器视觉领域内的研究问题 已成为一大热点。一些真彩图像显著性检测算法已被成功推广到波长范围为可见光的高光 谱图像，并取得了一定的效果。然而，这些方法通常将高光谱图像转换成真彩图像，再利用 真彩图像上的显著性检测算法进行检测，并未充分利用高光谱图像中丰富的光谱信息。近 年来，有研究表明利用高光谱数据中的光谱信息能够有效地检测场景中的显著性区域，而 无需将高光谱数据转换成真彩图像。
[0003] Jie Liang等人在文献"Salient object detection in hyperspectal imagery, IEEE International Conference on Image Processing(ICIP)，2013，pp.2393_2397"中 利用光谱向量的欧式距离和余弦相似度来进行显著性检测。具体是通过在输入数据上构建 尺度空间，并通过欧氏距离和余弦相似度分别计算像元及其邻域之间的差异。然而，该方法 以单个像元为显著性检测单元，并直接利用光谱向量进行相关计算；因此，其检测结果容易 受到光照等因素的影响，对图像中亮度高的区域过于敏感，检测准确率低。
【发明内容】
[0004] 为了克服现有高光谱图像显著性检测方法准确率低的不足，本发明提供一种基于 区域光谱梯度特征对比的高光谱图像显著性检测方法。该方法首先提取光谱梯度特征，并 在所得数据上对图像数据进行过分割，获得一组超像素;再使用聚类方法将在光谱梯度特 征上相似的超像素标记成同一类别，从而形成具有不同光谱梯度特征的区域。对每一区域， 利用区域对比方法，计算其与相邻区域的光谱梯度特征对比值作为其显著性响应值。在哈 佛(Harvard)大学和曼彻斯特(Manchester)大学提供的数据集上的实验结果表明，在查全 率为70%时，本发明方法查准率为82%，并且在查全率从0到100 %变化的过程中，本方法查 准率以超过80 %的概率高于【背景技术】方法。
[0005] 本发明解决其技术问题所采用的技术方案是:一种基于区域光谱梯度特征对比的 高光谱图像显著性检测方法，其特点是包括以下步骤：
[0006] 步骤一、对于输入的高光谱图像数据Xn，在光谱维上求梯度得到光谱梯度图像Gn = {gl，g2,…，gn}，其中第i个光谱梯度向量gi表示为
[0007]
[0008] 式中，Λλ为相邻波段的波长间隔。
[0009 ]步骤二、对光谱梯度图像Gn中的每一个梯度向量，在其末尾标记空间维坐标信息， 生成梯度扩展向I
。所有扩展向量成数据G\。使用简 单线性迭代聚类对G' n进行过分割获得一组超像素，分割过程如下：
[0010] (a)在空间维上以S为步长形成mXm2的网格，每个网格作为一个初始超像素，取每 个网格中心像素的扩展向量为该超像素的初始中心向量Ck = g' k，ke {1，2，…，nu Xm2}，赋 予网格内的所有扩展向量与Ck同样的类别标记；
[0011] (b)对于每个中心Ck，计算其与落在以Ck为中心大小为2SX2S的区域所对应的扩展 向量f」间的距离D( j，k)，计算式如下
[0012]
(1)
[0013]其中g、（l)和Ck(l)表示对应向量中的第1个分量，α为调节光谱梯度距离和空间距 离的系数；
[0014] (c)若D(j，k)小于g、与其当前归属的中心Cr之间的距离D(j，r)，则将g、的类别标 记为属于中心否则，保持g、的类别标记不变；
[0015] (d)通过步骤(b)、步骤(c)得到每个扩展向量归属的新类别标记，根据新的类 别标记计算当前每个中心&对应的新的超像素中心Ck;
[0016] (e)计算每组Ck和C'k之间的差异，并由此求出迭代误
[0017] (f)若err小于给定阈值，则结束计算；否则，更新当前每个中心Ck为k，并返回步 骤(b) 〇
[0018]得到一个由超像素中心向量所组成的集合C = {Ck| k=l，2,…,ι?ιΧπμ}，及每个梯 度向量的超像素类别标记Li。
[0019] 步骤三、对获得的这组超像素使用均值漂移算法进行聚类，将在光谱梯度上较相 似的超像素标记为同一类，从而形成若干区域。聚类时，第j个超像素中心与第k个超像素中 心之间的光谱梯度距离定义为
，它们之间的空间距离定义 为
2，Cj(l)和Ck(l)表示对应向量中的第1个分量。
[0020] 令k = l，具体聚类过程如下：
[0021 ] (g)从超像素中心向量集合C中取出第k个中心向量&;
[0022] (h)从 C 中选择(/ ={Cj|ds(j，k)〈TsAdc(j，k)〈Tc}，j = l，2,…，miXm2。其中 Ts 和 Tc 分别为给定的梯度向量距离阈值和空间距离阈值；
[0023] ⑴计算C'中各中心向量的均值匕，并计算其与&的差异err= ;
[0024] (j)令Ci^C/k，若err大于给定阈值，贝lj返回步骤(h);否则，转步骤(k);
[0025] (k)设当前聚类中心集合为Θ，遍历Θ并选择第一个满足光谱梯度距离小于!^/2 且空间距离小于Tc/2的中心Θ;
[0026] (1)若步骤(k)中没有满足条件的中心，则Θ = Θυ{&}，并将Ck标记为新的一类;否 则，Θ = (0+Ck)/2，将Ck标记为Θ所对应的类别。
[0027] (m)若C中所有成员已被遍历，则结束计算;否则，k = k+1，转步骤(g)。
[0028]将C中的超像素聚成若干类，再结合每个梯度向量的超像素类别标记。确定每个 梯度向量最终的类别标记L2;由此，Gn在空间维上被分割成若干区域。
[0029]步骤四、根据算每一类别中心所对应的扩展向量^，_]_6{1，2，一少}，其中~为 总的类数，每一个ζ」对应着一个空间区域心。对于区域心，其基于区域光谱特征对比的显著 性响应值为：
[0030](2) Kk
[0031] 其 4
为光谱梯度距离，ζ」（1)和Gk(l)表示对应向量中的第1个分量，ω (Rj)为Rj中的像元数，σ为用 于调节空间距离权重的系数。
[0032] 本发明的有益效果是:该方法首先提取光谱梯度特征，并在所得数据上对图像数 据进行过分割，获得一组超像素;再使用聚类方法将在光谱梯度特征上相似的超像素标记 成同一类别，从而形成具有不同光谱梯度特征的区域。对每一区域，利用区域对比方法，计 算其与相邻区域的光谱梯度特征对比值作为其显著性响应值。在哈佛(Harvard)大学和曼 彻斯特(Manchester)大学提供的数据集上的实验结果表明，在查全率为70%时，本发明方 法查准率为82 %，并且在查全率从0到100 %变化的过程中，本方法查准率以超过80 %的概 率高于【背景技术】方法。
[0033]以下结合【具体实施方式】详细说明本发明。
【具体实施方式】
[0034] 本发明基于区域光谱梯度特征对比的高光谱图像显著性检测方法具体步骤如下：
[0035] 高光谱图像是由空间维和光谱维构成的一个三维立方体结构，空间维反映了成像 场景中不同位置像元在某一光谱波段上的反射值;而光谱维反映像元在不同波段上的反射 值。因而高光谱图像可视为P幅2D图像(波段图像），其中p为波段数，若将每一个波段图像拉 成一个行向量，贝一幅高光谱图像可以表示成一pXn的二维矩阵X n= {xi，X2, . . .，xn}，其中 η为图像中像元总数。乂"的每一行对应某一波段上的所有像素，其每一列则是某一像素对应 的光谱向量。
[0036] 1、提取光谱梯度特征。
[0037] 对于输入的高光谱图像数据Χη，在光谱维上求梯度得到光谱梯度图像Gn={gl， 沿，"_，811}，其中第1个光谱梯度向量81可以表示：
其中Λλ为相邻波段的波长间隔。
[0038] 2、过分割获得超像素。
[0039] 对光谱梯度图像6"中的每一个梯度向量，在其末尾标记空间维坐标信息，生成梯 度扩展向量g = 彳e {1,2,…，《丨。所有扩展向量g'满成数据G'n。使用简单线性 迭代聚类（Simple Linear Iterative Clustering，SLIC)对G'n进行过分割获得一组超像 素，具体分割过程如下：
[0040] (a)在空间维上以S为步长形成mXm2的网格