本发明专利技术一种基于高阶统计量的高光谱图像目标检测方法，它有四大步骤：步骤一：计算机在ＭＡＴＬＡＢ?Ｒ２００８ｂ环境下读取高光谱图像数据；步骤二：计算机对数据进行预处理，即对数据进行去均值化和白化处理；步骤三：构造检测滤波器在对目标光谱增益为１的约束下最小化输出数据的高阶统计量，求解检测滤波器最佳权向量；步骤四：设定合适的阈值，获取检测结果图像。本发明专利技术克服了现有技术的不足，充分利用数据的高阶统计量，取得了较好的检测效果，特别是能够提高低虚警率情况下的检测概率，它在高光谱遥感图像目标检测技术领域里具有实用价值和广阔的应用前景。
【技术实现步骤摘要】
本专利技术涉及，属于高光谱遥感图像目标检测
技术介绍
随着高光谱成像技术在过去三十年的迅猛发展，高光谱图像(Hyperspectral Images)的分析和处理方法成为当前国内外遥感图像领域的研究热点之一。高光谱图像和传统遥感图像相比，突出特点是光谱分辨率高，可获得观测对象的几十个甚至上百个光谱波段的图像信息。成像光谱系统获得的连续波段宽度一般都在10nm以内，因此这种数据能以足够的光谱分辨率区分出那些具有诊断性光谱特征的地表物质，这也是高光谱技术用于探测甚至识别地物的物质基础(耿修瑞.高光谱遥感图像目标探测与分类技术研究.北京中国科学院遥感应用研究所，2005)。高光谱图像数据是一个数据立方体，其一层图像对应一个光谱波段；每个像素点位置对应一条光谱曲线(贺霖，潘泉等.高光谱图像目标检测研究进展.电子学报，2009，37(9)2016-2024)。因此，高光谱图像具有图谱合一的特性，即同时包含被观测物体的空间信息和光谱信息。 正是由于高光谱图像的高光谱分辨率，使得许多原来利用普通光学图像信息不能解决的问题，可以通过分析高光谱数据得到解决。高光谱图像目标检测是利用已知的目标光谱信息在高光谱图像中对感兴趣的目标进行检测、定位的技术(刘颖.高光谱图像目标识别新技术研究.哈尔滨哈尔滨工业大学，2007)。高光谱图像目标检测技术在军事和民用领域中都有重要的应用价值。在军事领域可用于对飞机、坦克等军事目标进行检测、定位。在民用领域也可用于房屋等目标的检测。 目前存在的高光谱图像目标检测方法，主要是基于统计方法。其中经典的有匹配滤波器(MF，Matched Filter)，能量约束最小化(CEM，Constrained Energy Minimization)，自适应匹配滤波(AMF，Adaptive Matched Filter)，自适应一致性估计(ACE，Adaptive Coherence Estimator)等。MF假设目标像素和背景像素分别服从均值不同但协方差相同的高斯分布，然后利用Fisher线性判别准则构造检测算子。CEM在对目标光谱增益为1的约束下，最小化输出数据的能量，求解出检测算子。AMF假设背景服从零均值高斯分布，包含背景的像素也服从高斯分布，协方差和背景的协方差相同，但均值由目标子空间决定。ACE假设背景服从零均值高斯分布，包含背景的像素也服从高斯分布，协方差和背景的协方差结构相同但方差不同，均值由目标子空间决定。AMF和ACE都是基于广义似然比检验的方法，根据像素服从的概率分布的假设进行广义似然比检验，得出检测结果。 目前存在的这些目标检测方法主要利用数据的二阶统计量进行计算，主要涉及到了数据的协方差阵或相关阵，而没有利用数据的高阶统计量。而数据的高阶统计量往往包含着数据更重要的特征信息，对检测来说是非常有用的。在复杂背景下，目标所占的像素往往很少，高阶统计量能够更好的描述目标的统计特征。本专利技术针对上述情况，提出一种基于高阶统计量的目标检测方法，充分利用了数据的高阶统计信息，取得了较好的效果。
技术实现思路
1、目的本专利技术的目的是提供，它充分利用数据的高阶统计量，取得了较好的检测效果，特别是能够提高低虚警率情况下的检测概率，这点对于实际应用是十分有利的。 2、技术方案本专利技术是通过以下技术方案实现的 本专利技术，它包括如下步骤 步骤一用计算机读取数据。计算机在MATLAB R2008b环境下读取高光谱图像数据，数据来源于成像光谱仪采集到的遥感图像，得到数据立方体。高光谱图像数据应去除被水汽吸收的波段和信噪比较低的波段。目标光谱来自已知的光谱库或对目标所在像素的光谱曲线取平均值得到。假设高光谱图像共有M个波段，像素的光谱曲线用向量形式表示x0＝T，x0i为该像素第i个波段的值。已知的目标光谱也用向量形式表示s0＝T，s0i为目标光谱第i个波段的值。 步骤二数据预处理。用计算机读取数据后需要对数据进行预处理，即对数据进行去均值化和白化。 (1)去均值化。对高光谱图像去均值化，每个像素的光谱曲线都减去所有像素光谱曲线的平均值，去均值后使得整个高光谱图像数据的均值为零。另外，已知的目标光谱也要减去均值。 去均值化可通过下式进行 x′＝x0-E{x0} 其中E{}表示取数学期望，可用图像各像素光谱曲线的平均值近似代替，即 其中x0i表示第i个像素的光谱曲线，N等于图像的像素个数。相应的已知的目标光谱也需要进行相应的处理，减去图像各像素光谱曲线的期望即 s′＝s0-E{x0} (2)数据白化。对图像进行去均值化以后，还需要对高光谱图像数据进行白化。白化是去相关的过程，白化之后像素的光谱曲线各分量之间是不相关的，即像素光谱曲线的协方差矩阵是单位阵。另外，目标光谱也要进行相应的处理，即用白化矩阵乘以目标光谱。白化可通过下式进行 x＝Γ-1/2x′ 其中Γ为数据的协方差矩阵，因为已进行了去均值化，协方差矩阵等于自相关阵，即 xi′表示去均值化后的第i个像素的光谱曲线，T表示转置。相应的目标光谱也需要进行相应的处理即 s＝Γ-1/2s′ 数据经过去均值化和白化处理后均值为零，协方差矩阵是单位阵即 E{x}＝0 E{xxT}＝I 步骤三求解检测滤波器最佳权向量。数据进行去均值化和白化后，接下来的检测过程可以看作是一个滤波的过程。检测滤波器可以写为 y＝wTx x是每个像素的光谱曲线，y是滤波器的输出，w＝T是滤波器的权向量。接下来就是要在对目标光谱增益为1的约束下，最小化输出数据的高阶统计量，求出最佳权向量。这样，最佳权向量w的求解问题可以写为 其中，s为经过预处理后的目标光谱，y＝wTx，G为高阶统计量函数。本专利技术选取了四种高阶统计量，即G(y)＝y3，G(y)＝y4，G(y)＝tanh(y)，G(y)＝sign(y)。 求解最佳权向量本质上是条件极值的求解问题，可用拉格朗日乘子法转化为无条件极值的问题，即求下列函数的极值 J(w)＝E{G(wTx)}-λ(wTs-1) 等价于求解下列方程组 可从方程组解得λ＝E{(wTx)G′(wTx)}。本专利技术采用梯度下降法求解函数J(w)的极值，从而求解出最佳权向量w。具体步骤如下 (1)、初始化w。w的初始值可以随机给定，然后进行归一化，即 w+＝rand(w) w＝w+/‖w+‖ (2)、迭代w。利用梯度下降法对w进行迭代，迭代规则为 λ＝E{(wTx)G′(wTx)} w+＝w-μ▽J＝w-μ{E{G′(wTx)x}-λs} w＝w+/‖w+‖ 其中，μ为步长，本专利技术取μ＝10-4。 (3)、停止迭代条件。当两次迭代的权向量w在相邻两次迭代的过程中变化不大时停止迭代，本专利技术中停止条件为 ‖w-wold‖＜tol wold表示上一次迭代中w的值。本专利技术中取tol＝10-4。 步骤四获取检测结果图像。在求解出最佳权向量后，将每个像素的光谱曲线通过检测滤波器，得出输出数据。然后设定合适的阈值，滤波器输出大于阈值的判定该像素存在目标，小于阈值的判定该像素不存在目标。将输出大于阈值所对应的像素本文档来自技高网...
【技术保护点】
一种基于高阶统计量的高光谱图像目标检测方法，其特征在于：该方法具体步骤如下：步骤一：用计算机读取数据；计算机在ＭＡＴＬＡＢＲ２００８ｂ环境下读取高光谱图像数据，数据来源于成像光谱仪采集到的遥感图像，得到数据立方体；高光谱图像数据应去除被水汽吸收的波段和信噪比较低的波段，目标光谱来自已知的光谱库或对目标所在像素的光谱曲线取平均值得到；假设高光谱图像共有Ｍ个波段，像素的光谱曲线用向量形式表示ｘ↓［０］＝［ｘ↓［０１］，ｘ↓［０２］，．．．，ｘ↓［０Ｍ］］↑［Ｔ］，ｘ↓［０ｉ］为该ｈ（ｙ），Ｇ（ｙ）＝ｓｉｇｎ（ｙ）；求解最佳权向量本质上是条件极值的求解问题，可用拉格朗日乘子法转化为无条件极值的问题，即求下列函数的极值：Ｊ（ｗ）＝Ｅ｛Ｇ（ｗ↑［Ｔ］ｘ）｝－λ（ｗ↑［Ｔ］ｓ－１）等价于求解下列方程组：＊＊＊可从方程组解得λ＝Ｅ｛（ｗ↑［Ｔ］ｘ）Ｇ′（ｗ↑［Ｔ］ｘ）｝；本专利技术采用梯度下降法求解函数Ｊ（ｗ）的极值，从而求解出最佳权向量ｗ；其具体步骤如下：（１）、初始化ｗ；ｗ的初始值可以随机给定，然后进行归一化，即：ｗ↑［＋］＝ｒａｎｄ（ｗ）ｗ＝ｗ↑［＋］／‖ｗ↑［＋］‖（２）、迭代ｗ；利用梯度下降法对ｗ进行迭代，迭代规则为：λ＝Ｅ｛（ｗ↑［Ｔ］ｘ）Ｇ′（ｗ↑［Ｔ］ｘ）｝ｗ↑［＋］＝ｗ－μ▽↓［Ｊ］＝ｗ－μ｛Ｅ｛Ｇ′（ｗ↑［Ｔ］ｘ）ｘ｝－λｓ｝ｗ＝ｗ↑［＋］／‖ｗ↑［＋］‖其中，μ为步长，本专利技术取μ＝１０↑［－４］；（３）、停止迭代条件；当两次迭代的权向量ｗ在相邻两次迭代的过程中变化不大时停止迭代，本专利技术中停止条件为：‖ｗ－ｗ↓［ｏｌｄ］‖＜ｔｏｌｗ↓［ｏｌｄ］表示上一次迭代中ｗ的值，本专利技术中取ｔｏｌ＝１０↑［－４］；步骤四：获取检测结果图像；在求解出最佳权向量后，将每个像素的光谱曲线通过检测滤波器，得出输出数据；然后设定合适的阈值，滤波器输出大于阈值的判定该像素存在目标，小于阈值的判定该像素不存在目标；将输出大于阈值所对应的像素灰度值设为２５５，输出小于阈值所对应的像素灰度值设为０，就得到了检测结果的二值图像；在二值图像中目标对应的区域为白色，非目标对应的区域为黑色，从而完成了对目标的检测和定位。像素第ｉ个波段的值，已知的目标光谱也用向量形式表示ｓ↓［０］＝［ｓ↓［０１］，ｓ↓［０２］，．．．，ｓ↓［０Ｍ］］↑［Ｔ］，ｓ↓［０ｉ］为目标光谱第ｉ个波段的值；步骤二：数据预处理；用计算机读取数据后需要对数据进行预处理，即对数据进...
【技术特征摘要】
【专利技术属性】
技术研发人员：史振威，杨硕，姜志国，赵卫，王扬，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：11[中国|北京]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。