发明名称 高维空间定向投影端元提取方法
摘要 本发明涉及一种高维空间定向投影端元提取方法,包括以下步骤:1)输入图像和端元数量;2)初始化投影平面x轴对应的基向量;3)计算投影平面y轴对应的基向量;4)计算每个像元在投影平面内的坐标,按坐标将投影点显示在可视化窗口中;5)在可视化窗口中勾选出新的端元;6)根据新端元计算投影平面x轴对应的基向量;7)判断端元数量是否达到要求;8)输出所有端元:输出所有的<img file="dsa00000483607400011.GIF" wi="124" he="58" />作为端元。本发明的有益效果为:对投影后的数据进行定量分析,将最可能是下一个端元的数据点进行标记,指导操作人员选择,也允许操作人员根据实际情况自由选择;在不改变顶点选择原理的基础上降低了选择的难度,使端元提取变得更加直观、定量化、自动化。
申请公布号 CN102184400A 申请公布日期 2011.09.14
申请号 CN201110107797.1 申请日期 2011.04.28
申请人 中国科学院对地观测与数字地球科学中心 发明人 张兵;高连如;孙旭;吴远峰;张文娟;申茜
分类号 G06K9/46(2006.01)I 主分类号 G06K9/46(2006.01)I
代理机构 北京纽乐康知识产权代理事务所 11210 代理人 唐忠庆
主权项 1.一种高维空间定向投影端元提取方法,其特征在于,包括以下步骤:1)输入图像和端元数量:输入包含n个像元的L波段高光谱图像像元集合<img file="FSA00000483607600011.GIF" wi="102" he="51" />和端元数量m;将循环计数器j赋值为1,x方向投影向量w<sub>1</sub>赋值为∈<sub>1</sub>,∈<sub>1</sub>是第1个元素为1其余元素为0的L维向量,L是遥感图像波段数;2)初始化投影平面x轴对应的基向量:如果j=1,将w<sub>2</sub>赋值为∈<sub>2</sub>;如果j不等于1,将w<sub>2</sub>赋值为u<sub>j-1</sub>;3)计算投影平面y轴对应的基向量:将|(w<sub>1</sub>,r<sub>i</sub>)|赋值给x<sub>i</sub>,将(w<sub>2</sub>,r<sub>i</sub>)赋值给y<sub>i</sub>;其中(x<sub>i</sub>,y<sub>i</sub>)为本次循环像元r<sub>i</sub>在投影平面上的坐标,而{w<sub>1</sub>,w<sub>2</sub>}为这个投影平面的一组基底;4)计算每个像元在投影平面内的坐标,按坐标将投影点显示在可视化窗口中;具体步骤如下:在显示窗口中建立平面直角坐标系(白色背景、黑色坐标轴),为每个像元r<sub>i</sub>按照其在步骤3)中求得的坐标(x<sub>i</sub>,y<sub>i</sub>),在坐标系中绘制一个蓝色点;然后在显示窗口中将会看到平面上的n个点;把这些点中最右侧的一个,也就是x<sub>i</sub>最大的那个点用红色圆圈圈出,x<sub>i</sub>坐标记为(x<sub>p</sub>,y<sub>p</sub>);则(x<sub>p</sub>,y<sub>p</sub>)对应的像元r<sub>p</sub>是下一个端元;5)在可视化窗口中勾选出新的端元:在显示窗口中绘制一条封闭的曲线,表示选中曲线中包围的蓝色点;将这些蓝色的点对应的像元(L维向量)取平均值,作为下一个端元e<sub>j</sub>;6)根据新端元计算投影平面x轴对应的基向量:如果j等于1,将u<sub>j</sub>赋值为e<sub>j</sub>;如果j不等于1,将u<sub>j</sub>赋值为<img file="FSA00000483607600012.GIF" wi="379" he="131" />7)判断端元数量是否达到要求:将j+1赋值给j,把<img file="FSA00000483607600013.GIF" wi="357" he="128" />赋值给w<sub>1</sub>;如果j≤m,则转到步骤3);否则,转下一步骤;8)输出所有端元:输出所有的<img file="FSA00000483607600014.GIF" wi="122" he="57" />作为端元。
地址 100094 北京市海淀区邓庄南路9号