高维空间定向投影端元提取方法技术

本发明专利技术涉及一种高维空间定向投影端元提取方法，包括以下步骤：1)输入图像和端元数量；2)初始化投影平面x轴对应的基向量；3)计算投影平面y轴对应的基向量；4)计算每个像元在投影平面内的坐标，按坐标将投影点显示在可视化窗口中；5)在可视化窗口中勾选出新的端元；6)根据新端元计算投影平面x轴对应的基向量；7)判断端元数量是否达到要求；8)输出所有端元：输出所有的作为端元。本发明专利技术的有益效果为：对投影后的数据进行定量分析，将最可能是下一个端元的数据点进行标记，指导操作人员选择，也允许操作人员根据实际情况自由选择；在不改变顶点选择原理的基础上降低了选择的难度，使端元提取变得更加直观、定量化、自动化。

【技术实现步骤摘要】

本专利技术涉及高维空间定向投影图像分析领域，尤其涉及一种。
技术介绍
高光谱图像进行端元提取是依据线性光谱混合模型的几何学解释，认为混合像元分布在以端元为顶点的高维空间单形体中，而单形体在低维空间的投影一定是凸集且投影凸集的顶点一定是单形体顶点的投影。N维可视化技术将高光谱图像的所有像元点从高维特征空间投影到二维平面上，并允许用户人工选取投影点中可能是顶点的那些点，通过不断地随机生成投影平面，可以改变像元点的投影位置，避免漏选或错选。但是N维可视化技术存在以下问题1、N维可视化中投影平面随机产生，导致顶点选择缺乏定量的评价依据； 2、主要靠操作人员人工判读，随意性比较大，难以提高精度；3、对操作人员理论水平要求比较高，难以向具有实际应用需求但不具备大量专业知识的使用人员普及。
技术实现思路
本专利技术的目的是提供一种，本方法是一种具有导向性的半自动端元提取方法，提供一种可以评价一个投影点是否为原单形体顶点的定量指标，并通过可视化的手段向用户进行展示，同时根据这个定量指标自动判断端元并在可视化界面中对用户作出提示，以克服现有N维可视化技术存在的上述不足。为达到上述目的，本专利技术的技术方案提供一种， 所述方法包括以下步骤1)输入图像和端元数量输入包含η个像元的L波段高光谱图像像元集合忉}；=1和端元数量m;将循环计数器j赋值为1，χ方向投影向量W1赋值为e 工是第1个元素为 1其余元素为0的L维向量，L是遥感图像波段数；2)初始化投影平面χ轴对应的基向量如果j = 1，将％赋值为e 2 ；如果j不等于1，将《2赋值为I^1 ；3)计算投影平面y轴对应的基向量将I (W15Ti) I赋值给Xi，将O^ri)赋值给； 其中(Xyyi)为本次循环像元α在投影平面上的坐标，而{wi，w2}为这个投影平面的一组基底；4)计算每个像元在投影平面内的坐标，按坐标将投影点显示在可视化窗口中；具体步骤如下在显示窗口中建立平面直角坐标系(白色背景、黑色坐标轴)，为每个像元A按照其在步骤幻中求得的坐标(Xi，yi)，在坐标系中绘制一个蓝色点；然后在显示窗口中将会看到平面上的η个点；把这些点中最右侧的一个，也就是Xi最大的那个点用红色圆圈圈出，Xi 坐标记为(xp，yp)；则(χρ, yp)对应的像元rp是下一个端元；5)在可视化窗口中勾选出新的端元在显示窗口中绘制一条封闭的曲线，表示选中曲线中包围的蓝色点；将这些蓝色的点对应的像元(L维向量)取平均值，作为下一个端元eJ ；6)根据新端元计算投影平面χ轴对应的基向量如果j等于1，将赋值为h ；如/-ι (e u )果j不等于1，将Uj赋值为ey·；7)判断端元数量是否达到要求将j+Ι赋值给九把。-赋值给 ；如果j ( m，则转步骤3)；否则，转下一步骤；8)输出所有端元输出所有的；^·。作为端元。本专利技术的有益效果为1)本专利技术基于凸面单形体的结构特点，在投影平面的选择上使用与已有端元正交的平面，保证各个数据点在投影平面上的特征与已有端元特征的差异最大化，优于N维可视化方法中的随机选择投影平面；2)本专利技术对投影后的数据进行定量分析，将最可能是下一个端元的数据点进行标记，指导操作人员选择，同时也允许操作人员根据实际情况自由选择，优于N维可视化方法中的完全自由选择；3)本专利技术在不改变顶点选择原理的基础上降低了选择的难度，使得端元提取变得更加直观、定量化、自动化。附图说明下面根据附图对本专利技术作进一步详细说明。图1是本专利技术实施例所述的的流程示意图。 具体实施例方式如图1所示，本专利技术实施例所述的一种，所述方法包括以下步骤1)输入图像和端元数量输入包含η个像元的L波段高光谱图像像元集合kK^P 端元数量m;将循环计数器j赋值为1，χ方向投影向量W1赋值为e 工是第1个元素为 1其余元素为0的L维向量，L是遥感图像波段数；2)初始化投影平面χ轴对应的基向量如果j = 1，将％赋值为e 2 ；如果j不等于1，将《2赋值为I^1 ；3)计算投影平面y轴对应的基向量将I (W15Ti) I赋值给Xi，将O^ri)赋值给； 其中(Xyyi)为本次循环像元α在投影平面上的坐标，而{wi，w2}为这个投影平面的一组基底；4)计算每个像元在投影平面内的坐标，按坐标将投影点显示在可视化窗口中；具体步骤如下在显示窗口中建立平面直角坐标系(白色背景、黑色坐标轴)，为每个像元A按照其在步骤幻中求得的坐标(Xi，yi)，在坐标系中绘制一个蓝色点；然后在显示窗口中将会看到平面上的η个点；把这些点中最右侧的一个，也就是Xi最大的那个点用红色圆圈圈出，Xi坐标记为(xp，yp)；则(χρ, yp)对应的像元rp是下一个端元；5)在可视化窗口中勾选出新的端元在显示窗口中绘制一条封闭的曲线，表示选中曲线中包围的蓝色点；将这些蓝色的点对应的像元(L维向量)取平均值，作为下一个端元eJ ；6)根据新端元计算投影平面χ轴对应的基向量如果j等于1，将赋值为。；如本文档来自技高网...

【技术保护点】
１．一种高维空间定向投影端元提取方法，其特征在于，包括以下步骤：１）输入图像和端元数量：输入包含ｎ个像元的Ｌ波段高光谱图像像元集合和端元数量ｍ；将循环计数器ｊ赋值为１，ｘ方向投影向量ｗ１赋值为∈１，∈１是第１个元素为１其余元素为０的Ｌ维向量，Ｌ是遥感图像波段数；２）初始化投影平面ｘ轴对应的基向量：如果ｊ＝１，将ｗ２赋值为∈２；如果ｊ不等于１，将ｗ２赋值为ｕｊ－１；３）计算投影平面ｙ轴对应的基向量：将｜（ｗ１，ｒｉ）｜赋值给ｘｉ，将（ｗ２，ｒｉ）赋值给ｙｉ；其中（ｘｉ，ｙｉ）为本次循环像元ｒｉ在投影平面上的坐标，而｛ｗ１，ｗ２｝为这个投影平面的一组基底；４）计算每个像元在投影平面内的坐标，按坐标将投影点显示在可视化窗口中；具体步骤如下：在显示窗口中建立平面直角坐标系（白色背景、黑色坐标轴），为每个像元ｒｉ按照其在步骤３）中求得的坐标（ｘｉ，ｙｉ），在坐标系中绘制一个蓝色点；然后在显示窗口中将会看到平面上的ｎ个点；把这些点中最右侧的一个，也就是ｘｉ最大的那个点用红色圆圈圈出，ｘｉ坐标记为（ｘｐ，ｙｐ）；则（ｘｐ，ｙｐ）对应的像元ｒｐ是下一个端元；５）在可视化窗口中勾选出新的端元：在显示窗口中绘制一条封闭的曲线，表示选中曲线中包围的蓝色点；将这些蓝色的点对应的像元（Ｌ维向量）取平均值，作为下一个端元ｅｊ；６）根据新端元计算投影平面ｘ轴对应的基向量：如果ｊ等于１，将ｕｊ赋值为ｅｊ；如果ｊ不等于１，将ｕｊ赋值为７）判断端元数量是否达到要求：将ｊ＋１赋值给ｊ，把赋值给ｗ１；如果ｊ≤ｍ，则转到步骤３）；否则，转下一步骤；８）输出所有端元：输出所有的作为端元。...

【技术特征摘要】
1. 一种高维空间定向投影端元提取方法，其特征在于，包括以下步骤1)输入图像和端元数量输入包含η个像元的L波段高光谱图像像元集合^KL1和端元数量m;将循环计数器j赋值为1，χ方向投影向量W1赋值为e 工是第1个元素为1其余元素为0的L维向量，L是遥感图像波段数；2)初始化投影平面χ轴对应的基向量如果j= 1，将W2赋值为e 2 ；如果j不等于1， 将W2赋值为I^1 ；3)计算投影平面y轴对应的基向量将I(wi; r,) I赋值给Xi，将(w2，r,)赋值给Ji ；其中Ui, Yi)为本次循环像元A在投影平面上的坐标，而Iw1, w2}为这个投影平面的一组基底；4)计算每个像元在投影平面内的坐...

【专利技术属性】
技术研发人员：张兵，高连如，孙旭，吴远峰，张文娟，申茜，
申请(专利权)人：中国科学院对地观测与数字地球科学中心，
类型：发明
国别省市：11