基于线性最小二乘支持向量机的高光谱图像端元选择方法技术

本发明专利技术提供的是一种基于线性最小二乘支持向量机的高光谱图像端元选择方法。步骤１．选取前Ｎ个像素点作为初始端元；步骤２．令目前所选端元中的第ｉ个为“１”类，其余Ｎ－１个为“０”类，首次执行ｉ＝１，建立相应的ＬＬＳＳＶＭ判别函数即距离测算函数；步骤３．依次计算每个像素的距离，如果某个像素的绝对距离大于１，则将该像素替换第ｉ个端元，置ｉ＝１，转入步骤２；步骤４．ｉ＝ｉ＋１，若ｉ＞Ｎ，转入步骤５，否则转入步骤２；步骤５．当前端元即为最终选择端元，结束。本发明专利技术为一种高光谱图像端元选择方法，采用ＬＬＳＶＭ作为主要工具来完成。无需降维预处理、复杂度低。

【技术实现步骤摘要】

本专利技术涉及一种高光谱图像的端元选择方法，特别是一种基于支持向量机(SVM) 的端元选择方法，属于遥感信息处理

技术介绍
高光谱图像的空间分辨率一般较低，这种情况导致了混合像元的广泛存在，即一 个像元可能是几种类别的混合。分析各类别成分在混合像元内所占的比例的技术称为光谱 解混，是高光谱数据分析的最基本、最重要内容之一。光谱解混实施的必要前提是要知道高 光谱数据中包含哪些地物类别，在此需求下提取各类别代表性纯光谱的技术称为光谱端元 选择，简称端元选择。在近十多年里，多种高光谱图像端元选择方法相继发展起来。N-FINDR 是基于N维空间谱凸多面体的搜索光谱端元的经典例子，因其具有全自动、无参数、选择效 果较好等优点而受到广泛欢迎。但该方法需要进行数据降维预处理，且包含大量的体积计 算，这也是它最为耗时的部分。并且，体积计算(即主要为行列式的计算)的复杂度将随着 所选择的光谱端元数目增大而呈现立方增长，从而导致算法运算速度大大降低。 目前已有一些典型文献提出了对N-FINDR算法的改进方案。WUCHA0-CHENG等采用 像元预选的方式来降低后续搜索的复杂性，是从侧面来降低算法计算量。CHOWDHURYA.等 利用顺次选择的方式来代替联合选择的方式，这种顺次选取的方式远离了 N-FINDR算法的 基本特征，像元一经选定便无法更新，光谱端元之间的相互依赖关系也无法得到最大满足。 TA0 XUETA0等提出的方法可以直接在原始数据空间上进行而免于降维预处理，因此选出的 光谱端元更具合理性，在理论上突破了 N-FINDR算法需要降维处理的传统模式，但该方法 也属于顺次选取。 另一方面，基于凸几何分析的端元选择方法容易受到野值点的影B向，而野值点在 高光谱图像中大量存在，现有文献对此并未提出相应的解决方案。
技术实现思路
本专利技术的目的在于提供一种一种无需降维预处理、低复杂度的基于线性最小二乘支持向量机的高光谱图像端元选择方法。 本专利技术的目的是这样实现的 步骤1.选取前N个像素点作为初始端元； 步骤2.令目前所选端元中的第i个为"1"类，其余N-l个为"0"类，首次执行i =l，建立相应的LLSSVM判别函数即距离测算函数 f i (s) = 〈w*， s》b* ; 步骤3.依次计算每个像素的距离，如果某个像素的绝对距离大于l，则将该像素 替换第i个端元，置i二l，步骤2; 步骤4. l，若i 〉N，转入步骤5，否则转入步骤2; 步骤5.当前端元即为最终选择端元，结束。 本专利技术还可以包括 1、在步骤1之前加入根据局域像元密度对野值点进行检测和去除过程，包括以每 个像元点为中心建立固定尺寸的邻域窗，计算方邻域窗内所包含的像元点数作为中心点的 孤立程度度量指标，孤立程度度量指标越大的点作为野值点被去除。 2、在步骤1之前加入利用光谱分量值直按进行像元预排序过程，包括a)由原始 数据空间的第一维到最后一维依次选出和排列对应于极大值和极小值坐标的点对；b)从 余下的数据空间中进行第一步操作；c)继续进行这样的过程，直至所有的数据点都被选 出、排列；待全部数据点排序完毕之后，排在最前面的若干特征数据点被选作初始端元，其 迭代更新过程也将按照排序结果依次进行。 3、在步骤1之前加入根据局域像元密度对野值点进行检测和去除过程和加入利 用光谱分量值直接进行像元预排序过程；所述根据局域像元密度对野值点进行检测和去除 过程包括以每个像元点为中心建立固定尺寸的邻域窗，计算方邻域窗内所包含的像元点 数作为中心点的孤立程度度量指标，孤立程度度量指标越大的点作为野值点被去除；所述 利用光谱分量值直接进行像元预排序过程包括a)由原始数据空间的第一维到最后一维 依次选出和排列对应于极大值和极小值坐标的点对；b)从余下的数据空间中进行第一步 操作；c)继续进行这样的过程，直至所有的数据点都被选出、排列；待全部数据点排序完毕 之后，排在最前面的若干特征数据点被选作初始端元，其迭代更新过程也将按照排序结果 依次进行。 4、在步骤1之前加入去除过程和加入利用光谱分量值直接进行像元预排序过程 和根据局域像元密度对野值点进行检测；所述利用光谱分量值直接进行像元预排序过程包 括a)由原始数据空间的第一维到最后一维依次选出和排列对应于极大值和极小值坐标 的点对；b)从余下的数据空间中进行第一步操作；c)继续进行这样的过程，直至所有的数 据点都被选出、排列；待全部数据点排序完毕之后，排在最前面的若干特征数据点被选作初 始端元，其迭代更新过程也将按照排序结果依次进行；所述根据局域像元密度对野值点进 行检测和去除过程包括以每个像元点为中心建立固定尺寸的邻域窗，计算方邻域窗内所 包含的像元点数作为中心点的孤立程度度量指标，孤立程度度量指标越大的点作为野值点 被去除。 本专利技术为一种高光谱图像端元选择方法，采用LLSVM作为主要工具来完成。本发 明在深入分析线性最小二乘支持向量机(LINEAR LEAST SQUARESUPPORT VECTOR MACHINE, 简记为LLSSVM)模型的基础上，提出建立免于降维预处理、免于体积计算、并对野值点干扰 具有鲁棒性的端元提取方法。无需降维预处理、复杂度低。 野值点预删除和像元预排序可以与基于LLSSVM进行端元选择算法结合使用，也 可以单独或联合用于其它端元选择算法中来提高其效果、效率。利用LLSSVM进行端元选择 过程可以推广到其他有关空间距离计算方面的应用。附图说明 图1为利用距离尺度替换体积尺度的原理图。 图2为2维合成数据及端元选择结果。 图3为真实高光谱图像在不同方法下的端元选择结果对比。其中图3a)为真实光谱，图3b)是未经野值点去除的端元选择结果，图3c)是经过野值点去除的端元选择结果。 图4的表1为1000点2维合成数据端元选择的运行时间/迭代次数比较。 图5的表2为10000点9维真实高光谱数据端元选择的运行时间/迭代次数比较。具体实施例方式下面结合附图举例对本专利技术做更详细地描述 本专利技术为一种高光谱图像端元选择方法，采用LLSVM作为主要理论工具来完成。 对野值点进行检测和去除基本过程如下以每个像元点为中心建立固定尺寸的邻域窗，通 过计算邻域窗内所包含的像元点数作为中心点的孤立程度度量指标来剔出野值点。像元预 排序基本过程如下利用投影统计的方法进行像元预排序，投影落在每个方向最远端次数 越多的点越被排在前面。端元选择基本过程如下根据LLSSVM具有相对距离测试功能，而 相对距离的比较和相应的体积比较具有一致性，建立基于LLSSVM的最大体积寻优迭代方 法。下面给出本专利技术的详细说明 1)根据局域像元密度对野值点进行检测和去除实施步骤如下采用邻域分析的 方法来确认并去除野值点。野值点通常以更加孤立的方式存在。这样，以每个像元点为中 心建立固定尺寸的邻域窗，通过计算邻域窗内所包含的像元点数作为中心点的孤立程度度 量指标。孤立程度度量指标越大的点作为野值点的可能性就越大而被去除。这里采用方邻 域(高维盒子)替换常用的圆邻域以降低计算量。 2)利用光谱分量值直接进行像元预排序实施步骤如下为了能够获得快速收敛 的重要端元迭代搜索，将每个像元点根据其潜在的纯度进行预先评价和排序。根据光谱端本文档来自技高网...

【技术保护点】
一种基于线性最小二乘支持向量机的高光谱图像端元选择方法，其特征是：步骤１．选取前Ｎ个像素点作为初始端元；步骤２．令目前所选端元中的第ｉ个为“１”类，其余Ｎ－１个为“０”类，首次执行ｉ＝１，建立相应的ＬＬＳＳＶＭ判别函数即距离测算函数：ｆ↓［ｉ］（ｓ）＝〈ｗ↑［＊］，ｓ〉＋ｂ↑［＊］；步骤３．依次计算每个像素的距离，如果某个像素的绝对距离大于１，则将该像素替换第ｉ个端元，置ｉ＝１，转入步骤２；步骤４．１，若ｉ＞Ｎ，转入步骤５，否则转入步骤２；步骤５．当前端元即为最终选择端元，结束。

【技术特征摘要】

【专利技术属性】
技术研发人员：王立国，张晶，邓禄群，赵春晖，乔玉龙，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：93[中国|哈尔滨]