专利名称：高维空间定向投影端元提取方法
技术领域：
本发明涉及高维空间定向投影图像分析领域，尤其涉及一种高维空间定向投影端元提取方法。
背景技术：
高光谱图像进行端元提取是依据线性光谱混合模型的几何学解释，认为混合像元分布在以端元为顶点的高维空间单形体中，而单形体在低维空间的投影一定是凸集且投影凸集的顶点一定是单形体顶点的投影。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)计算每个像元在投影平面内的坐标，按坐标将投影点显示在可视化窗口中；具体步骤如下在显示窗口中建立平面直角坐标系(白色背景、黑色坐标轴)，为每个像元巧按照其在步骤幻中求得的坐标(Xi，yi)，在坐标系中绘制一个蓝色点；然后在显示窗口中将会看到平面上的η个点；把这些点中最右侧的一个，也就是Xi最大的那个点用红色圆圈圈出，Xi坐标记为(ν yp)；则(ν yp)对应的像元rp是下一个端元；5)在可视化窗口中勾选出新的端元在显示窗口中绘制一条封闭的曲线，表示选中曲线中包围的蓝色点；将这些蓝色的点对应的像元(L维向量)取平均值，作为下一个端元6)根据新端元计算投影平面χ轴对应的基向量如果j等于1，将赋值为。；如果j不等于1，将Uj赋值为
全文摘要
本发明涉及一种高维空间定向投影端元提取方法，包括以下步骤1)输入图像和端元数量；2)初始化投影平面x轴对应的基向量；3)计算投影平面y轴对应的基向量；4)计算每个像元在投影平面内的坐标，按坐标将投影点显示在可视化窗口中；5)在可视化窗口中勾选出新的端元；6)根据新端元计算投影平面x轴对应的基向量；7)判断端元数量是否达到要求；8)输出所有端元输出所有的作为端元。本发明的有益效果为对投影后的数据进行定量分析，将最可能是下一个端元的数据点进行标记，指导操作人员选择，也允许操作人员根据实际情况自由选择；在不改变顶点选择原理的基础上降低了选择的难度，使端元提取变得更加直观、定量化、自动化。
文档编号G06K9/46GK102184400SQ20111010779
公开日2011年9月14日 申请日期2011年4月28日 优先权日2011年4月28日
发明者吴远峰, 孙旭, 张兵, 张文娟, 申茜, 高连如 申请人:中国科学院对地观测与数字地球科学中心