发明名称 一种SAM加权KEST高光谱异常检测算法
摘要 本发明公开了一种SAM加权KEST高光谱异常检测算法(SKEST)。首先,推导SKEST算法;其次采用双矩形窗口,对高光谱图像中每个像元计算其SKEST值,进行阈值分割,检测出异常点。SKEST算法在基于核的特征空间分离变换算法(KEST)基础上,采用光谱角匹配(SAM)测度对高维特征空间检测点邻域差异相关矩阵(DCOR)中每个样本引入权重因子,各样本权重因子取决于该样本光谱向量与检测窗口数据中心向量夹角。从而抑制检测窗口中的异常数据,突出主成分数据的贡献,使得DCOR矩阵能够更好地描述目标、背景数据分布差异;同时SAM对光谱能量是鲁棒的,结合径向基核函数,SKEST算法既考虑了信号的光谱能量差异,也考虑了信号的光谱曲线形状差异,因而更符合高光谱数据特性。
申请公布号 CN102663752A 申请公布日期 2012.09.12
申请号 CN201210104037.X 申请日期 2012.04.11
申请人 南京理工大学 发明人 柏连发;张毅;陈钱;顾国华;韩静;岳江;王博;徐杭威;祁伟;金左轮
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 1.一种SAM加权KEST高光谱异常检测算法,其特征在于:在基于核的特征空间分离变换KEST算法基础上,对高维特征空间DCOR矩阵中的每个样本引入权重因子,从而抑制检测窗口中的异常数据,具体实现步骤如下:步骤(1),输入高光谱图像数据;步骤(2),根据图像中异常目标大小设置双矩形窗尺寸;步骤(3),逐行扫描各个像素点,统计其在双矩形窗邻域内的目标光谱矩阵<img file="957265DEST_PATH_IMAGE001.GIF" wi="20" he="18" />和背景光谱矩阵<img file="347927DEST_PATH_IMAGE002.GIF" wi="16" he="18" />,以及目标样本数目<img file="333200DEST_PATH_IMAGE003.GIF" wi="21" he="25" />、背景样本数目<img file="174117DEST_PATH_IMAGE004.GIF" wi="24" he="25" />;其中高光谱数据波段数为<img file="697371DEST_PATH_IMAGE005.GIF" wi="17" he="18" />,定义<img file="891723DEST_PATH_IMAGE006.GIF" wi="45" he="25" />矩阵<img file="731503DEST_PATH_IMAGE007.GIF" wi="114" he="26" />表示目标数据,样本<img file="743322DEST_PATH_IMAGE008.GIF" wi="17" he="25" />为<img file="35763DEST_PATH_IMAGE005.GIF" wi="17" he="18" />维光谱向量<img file="846855DEST_PATH_IMAGE009.GIF" wi="121" he="28" />,<img file="947666DEST_PATH_IMAGE010.GIF" wi="84" he="25" />,<img file="599228DEST_PATH_IMAGE003.GIF" wi="21" he="25" />为目标样本个数;<img file="910123DEST_PATH_IMAGE011.GIF" wi="46" he="25" />矩阵<img file="961125DEST_PATH_IMAGE012.GIF" wi="113" he="26" />表示背景数据,样本<img file="509918DEST_PATH_IMAGE013.GIF" wi="20" he="26" />为<img file="801222DEST_PATH_IMAGE005.GIF" wi="17" he="18" />维光谱向量<img file="599413DEST_PATH_IMAGE014.GIF" wi="133" he="28" />,<img file="532734DEST_PATH_IMAGE015.GIF" wi="89" he="25" />,<img file="90361DEST_PATH_IMAGE004.GIF" wi="24" he="25" />为背景样本个数;步骤(4),根据各像素点的<img file="552567DEST_PATH_IMAGE001.GIF" wi="20" he="18" />、<img file="572475DEST_PATH_IMAGE002.GIF" wi="16" he="18" />,由式【1】计算其SKEST值;循环计算高光谱图像中所有像素点的SKEST值,得到SKEST图像;<img file="965279DEST_PATH_IMAGE016.GIF" wi="468" he="32" />【1】其中<img file="488665DEST_PATH_IMAGE017.GIF" wi="42" he="26" />即为SKEST值,<img file="184088DEST_PATH_IMAGE018.GIF" wi="13" he="14" />为检测点光谱向量,<img file="425714DEST_PATH_IMAGE019.GIF" wi="72" he="21" />,核函数<img file="169679DEST_PATH_IMAGE020.GIF" wi="202" he="30" />,<img file="704827DEST_PATH_IMAGE021.GIF" wi="17" he="16" />为常量参数;<img file="305573DEST_PATH_IMAGE022.GIF" wi="509" he="41" />【2】<img file="34495DEST_PATH_IMAGE023.GIF" wi="194" he="49" />【3】<img file="582151DEST_PATH_IMAGE024.GIF" wi="210" he="36" />为各像素光谱向量权重因子:<img file="80128DEST_PATH_IMAGE025.GIF" wi="331" he="190" />【4】<img file="241988DEST_PATH_IMAGE026.GIF" wi="46" he="25" />为<img file="458206DEST_PATH_IMAGE027.GIF" wi="168" he="90" />的特征向量<img file="606290DEST_PATH_IMAGE028.GIF" wi="273" he="130" />【5】步骤(5),对SKEST图进行阈值分割;设置阈值,逐行扫描各个像素点,若其SKEST值大于阈值即为异常点。
地址 210094 江苏省南京市孝陵卫200号