端元选择方法及操作

合集下载

改进的SGA端元选择的快速方法

改进的SGA端元选择的快速方法
Ab t a t s r c :SGA a e n a p p l ro e f risf Ⅱ a tmain a d h g e f in y p ro ma c .An i r v d S h s b e o u a n o t u u o to n ih. f c e c e r n e i f mp o e GA meh d wa r p s d a d a meho fc n tu tn y e pln sa o t d i h i h d me in ls c t o sp o o e n t d o o sr ci g a h p r a e wa d p e n t e h g - i nso a pa e.frt e 0 h pu p s fi r v n o s e lo t m ih i v le n u r b e v l me c lu ai n a d t us la s t o r o e o mp o i g a lw pe d a g r h wh c n o v s i n me a l ou a c l t n h e d o lw i o s e d o ac lto . Co lx v l me c mp rs n wa o v ne o smp e c mp rs n o h it n e fo a p i t p e fc l u ai n mp e o u o a io s c n e d t i l o a o ft e d sa c r m on i t y r ln n t i mpr v d meho o a h pe a e i h si p o e t d.Asa r s l.t e c mp e iy o a tS e u t h o lx t ffs GA sr d c d t ie ri p c i i e u e o l a n s a e d — n me so a i n in t l y.fo fr rc bi n o ii a GA.Ex e me t h w h tt e fs r m o me u c i rgn lS p r n s s o t a h a tSGA a he s me s lc in r s l i h st a e e t e u t o a rgna GA .whi he fr e nsg e ty f se h n t e lte n t r fe d mb rs l ci n s e d so i lS i l t o e m r r r al a t rt a h a tr i e ms o n me e ee to p e . u

三元选择表达式

三元选择表达式

三元选择表达式在计算机编程语言中,三元选择表达式(ternary operator)也被称为条件运算符(conditional operator),主要是用于给予条件的判断,从而返回不同的值。

即三元操作符会判断一个给定的条件,若条件为真,则返回一个值;若条件为假,则返回另一个值。

它通常的写法是"condition? value1 : value2",condition是用于判断的条件,value1和value2是两个分别对应的结果。

在一般的语言中,三元操作符都有明确定义的语法,大多数的开发语言都支持三元表达式,并提供了使用范例和语法规范。

常见的语言如JavaScript、PHP、Java、C# 等都有三元表达式。

三元选择表达式的语法规范如下:variableName = (expression) ? value1 : value2;在三元表达式中,条件表达式可以是任车数据类型,例如:整数、字符、布尔值等。

不过,一般来说,只有布尔值这种数据类型最常用于条件表达式中。

同时,value1 和 value2 也可以是任何类型数据,包括表达式。

在条件表达式为 true 时,所返回的值由 value1 所决定;在条件表达式为 false 时,所返回的值由 value2 所决定。

三元选择表达式的优点1. 增强了代码的可读性通过使用三元操作符,可以使代码更加简洁而且易于阅读。

使用三元表达式可以将多行的 if-else 语句简化为一行,易于编写和修改,也易于阅读和维护。

2. 缩短了代码的长度由于三元表达式将多行的 if-else 语句压缩至单行代码,所以能够大大缩减代码的长度,并在一定程度上减少了程序的复杂度。

3. 提高了代码的可维护性三元选择表达式作为一种简洁而高效的语句,可以大大提高代码的可读性和可维护性。

通过改善代码的结构和提高代码的可读性,可以使代码更加易于调试和修改。

三元选择表达式的缺点1. 可读性不强尽管三元选择表达式在简化代码长度和提高编写效率方面非常有用,但这种表达式并不总是可以使程序更加易于理解,特别是当条件和返回值非常复杂时。

dax中,多端表筛选一端表的计值规则

dax中,多端表筛选一端表的计值规则

dax中,多端表筛选一端表的计值规则【原创版】目录1.概述:DAX 中的多端表和一端表2.多端表筛选一端表的计值规则3.实际应用场景4.总结正文在数据分析领域,DAX(Data Analysis Expressions)是一种强大的数据处理和计算语言,被广泛应用于 Microsoft Power BI 和 Excel 等工具中。

在 DAX 中,我们经常会遇到多端表和一端表的概念。

为了更好地理解多端表筛选一端表的计值规则,我们先来了解一下这两种表的含义。

多端表是指一个表包含多个端,每个端可以看作是表的一个维度。

例如,一个销售表可以包含产品、时间、地区等多个端。

而一端表则是指一个表只包含一个端,例如,一个产品表只包含产品这个端。

在实际操作中,我们常常需要对多端表进行筛选,只保留其中一个端,并根据这个端进行计值。

例如,我们需要统计某个时间段内各个产品的销售总额。

在这种情况下,我们需要筛选出时间端,并根据这个端计算销售额。

DAX 中提供了丰富的函数和操作符来实现多端表筛选一端表的计值规则。

以下是一些常用的方法:1.使用"SELECT"函数:我们可以使用"SELECT"函数来选择需要的端。

例如,在上面的例子中,我们可以使用"SELECT(时间,产品)"来选择时间和产品这两个端。

2.使用"FILTER"函数:"FILTER"函数可以帮助我们筛选出需要的端。

例如,我们可以使用"FILTER(时间,时间 >= 开始时间 && 时间 <= 结束时间)"来筛选出指定时间段内的时间端。

3.使用"CALCULATE"函数:"CALCULATE"函数可以在筛选出的端上执行计算。

例如,我们可以使用"CALCULATE(销售额 = 销售额 * 数量)"来计算销售额。

端元名词解释

端元名词解释

端元名词解释
端元是指在计算机程序语言中用于标识程序中某个特定位置的
符号。

端元可以用于分隔程序中不同部分的标识符,例如在 C 语言中,分号就是一个常用的端元。

端元也可以用于表示程序中的操作,例如在Python 语言中,冒号就是用于标识代码块的端元。

端元通常是单个字符,但也可能是多个字符组成的符号。

端元在程序语言中扮演着非常重要的角色,因为它们可以帮助解析器理解程序的结构和含义。

例如,在C 语言中,分号是一个常用的端元,它用于在语句末尾标识一条语句的结束。

因此,在编写C 程序时,每条语句都必须以分号结尾。

同样,在Python 语言中,冒号是用于标识代码块的端元,因此在编写Python 程序时,必须在每个代码块的开头使用冒号。

一种高光谱遥感影像端元自动提取方法

一种高光谱遥感影像端元自动提取方法

An Automatic Endmember Extraction Algorithm from Hyperspectral Image
WANG Xiao2 ling ¹ , º , DU Pei 2 jun ¹ , TAN Kun ¹ , CAO Wen ¹ ( ¹ China Univer sity of Mining and T echnology , key labor a tory f or land Envir onment a nd Disaster Monitor ing of SBS M, Xuz hou 221116 ; º China univer sity of Mining and Techology , J ia ngsu key l abor ator y of Resour ces a nd Envir onment I nf orma tion Engineeing, Xuzhou 221116 ) Abstr act: An automat ic endmember extraction algorit hm is pr oposed based on unsuper vised classification, pixel pur ity in2 dex, liner spectra l mixing model and simplex of convex geometr y concepts. This pr oposed algorit hm can avoid the effects of un2 cert ainty, heavy wor kload and other shortcomings of existing artif icial sampling procedures. The approach is experimented by an example of OMIS hyperspectral image, and the experimental result indicates that the algorithm is effective and has high degree of automation. Key words: liner spectral mixtur e model; convex simplex; endmember

高光谱数据分析ENVI操作手册

高光谱数据分析ENVI操作手册

高光谱数据分析ENVI操作手册1.常见参数选择主菜单→File→Preferences●用户自定义文件(User Defined Files)图形颜色文件,颜色表文件,ENVI的菜单文件,地图投影文件等。

需重启ENVI ●默认文件目录(Default Directories)默认数据目录,临时文件目录,默认输出文件目录,ENVI补丁文件、光谱库文件、备用头文件目录等,需重启ENVI。

●显示设置(Display Default)可以设置三窗口中各个分窗口的显示大小,窗口显示式样等。

其中可以设置数据显示拉伸方式(Display Default Stretch),默认为2%线性拉伸。

●其他设置(Miscollaneous)制图单位(Page Unit),默认为英寸(Inches),可设置为厘米(Centimeters)还有缓冲大小(cache size),可以设置为物理内存的50-75%左右。

Image Tile Size不能超过4M。

2.显示图像及其波谱2.1.打开文件●主菜单,Open Image File→文件名.raw。

●或Window→Available Bands List→File →Open Image File→文件名.raw。

2.2.显示图像●显示单波段灰度级图像:Gray color,选择的波段一般是图像显示最清晰的波段。

●显示伪彩色图像:RGB color,选择具有明显吸收谷、强烈反射作用和所含信息量较大的波段作为彩色合成RGB波段。

●显示真彩色图像:波段列表(Available Bands List)中,右键→Load TrueColor 。

●图像保存:Display窗口,File→Save Image As→Image File,选择输出格式、路径和名称,OK。

●动画显示:Display窗口,Tools→Animation,动态显示各波段图像,能很快的分辨出包含信息量较多的波段。

混合像元分解及其应用(ENVI操作文本)讲解

混合像元分解及其应用(ENVI操作文本)讲解

一基于PPI的端元提取借助纯净像元指数(PPI)和n维可视化工具用于端元波谱收集。

第一步、获取纯净像元这个步骤是在MNF变换的结果上计算纯净像元指数(PPI),之后选择阈值范围从PPI图像上获得感兴趣区,感兴趣区包含的像元就是比较纯净的像元。

(1)打开高光谱数据。

(2)在ENVI主菜单中,选择Spectral ->MNF Rotation- > Forward MNF -> Estimate Noise Statistics From Dat a。

在标准ENVI文件选择对话框中,选择高光谱图像文件。

(3)打开Forward MNF Transform Parameters面板,选择MNF输出路径及文件名,单击OK执行MNF变换。

(4)在波段列表中输出MNF影像及特征曲线值。

从图中可以看出,大约20个波段以后的MNF的特征值很小(5)MNF变换后,在ENVI主菜单中,选择 Spectral-> Pixel Purity Index->[FAST] New Output Band。

在打开的Pixel Purity Index Input File对话框中,选择MNF变换结果,单击Spectral Subset按钮,选择前面10个波段(MNF后面波段基本为噪声),单击OK。

(6)在Display窗口中显示PPI结果。

选择Overlay->Region of Interest,在ROI Tool 面板中,选择Options->Band Threshold to ROI,选择PPI图像作为输入波段,单击OK,打开Band Threshold to ROI 面板(图14.19)。

Min Thresh Value:10,Max Thresh Value:空(PPI图像最大值),其他默认设置,单击OK计算感兴趣区,得到的感兴趣区显示在Display窗口中。

第二步、构建n维可视化窗口(1)在ENVI主菜单中,选择Spectral ->n-Dimensional Visualizer,在n-D Visualizer Input File 对话框中选择MNF变换结果,单击OK。

混合像元分解及其应用ENVI操作文本

混合像元分解及其应用ENVI操作文本

一基于PPI的端元提取借助纯净像元指数(PPI)和n维可视化工具用于端元波谱收集。

第一步、获取纯净像元这个步骤是在MNF变换的结果上计算纯净像元指数(PPI),之后选择阈值范围从PPI图像上获得感兴趣区,感兴趣区包含的像元就是比较纯净的像元。

(1)打开高光谱数据。

(2)在ENVI主菜单中,选择Spectral ->MNF Rotation- > Forward MNF -> Estimate Noise Statistics From Dat a。

在标准ENVI文件选择对话框中,选择高光谱图像文件。

(3)打开Forward MNF Transform Parameters面板,选择MNF输出路径及文件名,单击OK执行MNF变换。

(4)在波段列表中输出MNF影像及特征曲线值。

从图中可以看出,大约20个波段以后的MNF的特征值很小(5)MNF变换后,在ENVI主菜单中,选择 Spectral-> Pixel Purity Index->[FAST] New Output Band。

在打开的Pixel Purity Index Input File对话框中,选择MNF变换结果,单击Spectral Subset按钮,选择前面10个波段(MNF后面波段基本为噪声),单击OK。

(6)在Display窗口中显示PPI结果。

选择Overlay->Region of Interest,在ROI Tool 面板中,选择Options->Band Threshold to ROI,选择PPI图像作为输入波段,单击OK,打开Band Threshold to ROI 面板(图14.19)。

Min Thresh Value:10,Max Thresh Value:空(PPI图像最大值),其他默认设置,单击OK计算感兴趣区,得到的感兴趣区显示在Display窗口中。

第二步、构建n维可视化窗口(1)在ENVI主菜单中,选择Spectral ->n-Dimensional Visualizer,在n-D Visualizer Input File 对话框中选择MNF变换结果,单击OK。

端元提取方法

端元提取方法

端元提取是一种数学计算方法,用于从一组数据中提取一个特定的数值。

端元提取通常用于解决一些实际问题,例如从图像中提取特征、从音频信号中提取特征等。

端元提取的方法有很多种,其中最常用的方法是通过特征提取算法来提取端元。

特征提取算法是一种数学变换方法,通过将数据转化为新的形式来提取特征。

常见的特征提取算法包括傅里叶变换、小波变换、能量谱等。

在进行端元提取时,首先需要将数据输入到一个合适的特征提取算法中,然后根据算法的要求进行计算和分析,最后得到所需的端元。

不同的端元提取方法适用于不同的数据类型和问题,因此需要根据具体的情况选择合适的方法。

envi操作(2)

envi操作(2)

八、定义感兴趣区及分类主窗口—Classification:监督分类,非监督分类,决策树分类。

2)监督分类:按照分类以前自定义的样本进行分类。

2 样本选择:主影像窗口—Tools—Region Of Interest—ROI Tool 调出感兴趣区工具窗口进行样本选择,可以进行样本编辑(名称,颜色,填充方式等),样本选择越精确,分类结果越好。

2 选择分类方式:分类方式包括平行六面体法、最短距离法、马氏距离法、最大似然法、波谱角分类以及二进制编码法等,选择合适的分类方式。

2 引入影像—确定分类范围和波段—选择样本—给定阈值—确定存储路径和文件名—OK。

平行六面体法:用一条简单的判定规则对多波谱数据进行分类。

判定边界在图像数据空间中,形成了一个N维平行六面体。

平行六面体的维数由来自每一种选择的分类的平均值的标准差的阈值确定。

如果像元值位于N个被分类波段的低阈值与高阈值之间,则它归属于这一类。

如果像元值落在多个类里,那么ENVI将这一像元归到最后一个匹配的类里。

没有落在平行六面体的任何一类里的区域被称为无类别的。

最短距离法:用到每一个终端单元的均值矢量,计算每一个未知像元到每一类均值矢量的欧几里德距离。

所有像元都被归为最近的一类,除非限定了标准差和距离的极限(这时,会出现一些像元因不满足选择的标准,而成为“无类别”)。

马氏距离法:是一个方向灵敏的距离分类器,分类时用到了统计。

它与最大似然分类有些类似,但是假定所有类的协方差相等,所以是一种较快的方法。

所有像元都被归到最临近的ROI类,除非用户限定了一个距离阈值(这时,如果一些像元不在阈值内,就会被划为无类别)。

最大似然法:假定每个波段每一类统计呈均匀分布,并计算给定像元属于一特定类别的可能性。

除非选择一个可能性阈值,所有像元都将参与分类。

每一个像元被归到可能性最大的那一类里。

波谱角分类:(SAM)是一个基于自身的波谱分类,它是用n维角度将像元与参照波谱匹配。

前端数字区间选择方法

前端数字区间选择方法

前端数字区间选择方法
前端数字区间的选择方法,一般需要按照以下步骤进行:
1. 确定需求:首先需要明确数字区间的目的和用途,例如是为了满足特定的业务需求、提高用户体验、优化性能等。

2. 确定数据类型:根据需求,选择合适的数据类型。

例如,如果需要存储整数,可以选择int类型;如果需要存储浮点数,可以选择float或decimal
类型。

3. 确定范围:根据需求和数据类型,确定数字区间的范围。

例如,如果需要存储0到100的整数,可以选择int类型并设置数字区间的范围为0-100。

4. 考虑性能:在选择数字区间时,需要考虑性能因素。

如果数字区间的范围非常大,可能会导致查询效率降低或占用大量存储空间。

因此,需要根据实际情况进行权衡和优化。

5. 考虑可扩展性:在选择数字区间时,还需要考虑未来的扩展性。

如果未来数字区间的范围可能会发生变化,需要选择一种易于扩展的方案。

6. 测试和验证:在选择数字区间后,需要进行测试和验证,确保其满足需求并且性能良好。

总之,前端数字区间的选择方法需要根据实际情况进行综合考虑,权衡利弊,选择最适合的方案。

农作物种植面积遥感信息提取分析

农作物种植面积遥感信息提取分析

农作物种植面积遥感信息提取分析韦春月,欧阳义(桂林市临桂区农业农村局,广西桂林541100)摘要:农作物种植面积是我国最为基础的农业情报数据,也是我国粮食安全最为重要的参考指标之一,借助遥感信息技术对其种植面积进行信息提取,是我国当前对农作物种植面积进行信息提取的主要方法。

因此,以广西壮族自治区桂林市临桂区为例,对农作物种植面积遥感信息提取进行分析与研究。

首先对其临桂区区域概况进行分析,其次对其技术路线进行简析,最后对其农作物种植面积遥感信息提取进行分析,以供参考与借鉴。

关键词:农作物;种植面积;遥感信息;提取中图分类号:S127文献标识码:A文章编号:1005-7897(2023)06-0181-030引言美国国家航空航天局在20世纪70年代开展了遥感提取实验,以此来对美国的农作物基本需求进行快速分析,从而切实提升美国对其农业资源的管理质量及效率。

我国科研学者于20世纪70年代末期开始提倡遥感估产,其在多年的发展与完善过程中,其遥感估产的精准度得到了全面的提升,并借助专业的设备以及先进的技术,对我国农作物种植面积进行了全面的分析,以此来切实保证我国粮食安全,满足广大人民的基本粮食需求。

1临桂区区域概况临桂区地处我国广西壮族自治区东北部,桂林市的西南方,临桂区东北部,比邻桂林市老城区,南部与阳朔县相接壤,西部与融安县以及永福县相交。

临桂区的地势西高东低,有着“桂都首邑”的别称,临桂区区域总面积约为2193km2,下辖11个乡镇。

2016年总人口数量为50.86万人,农业人口数量为43.9万人,其实际耕地面积约为47.14万亩,其中包含28.47万亩的水田,农民人均年收入为14487元。

临桂区地处亚热带与温带的交界处,辖区内四季分明,光伏及热能较为充足,降雨量充沛,年降雨量约为1862mm,年平均气温为18~19益,全年无霜期较长,约为270~300d,全年平均湿度为76%,临桂区的土壤性质多为砂岩、页岩以及花岗岩共同构成的酸性红壤以及黄壤土,其土层较为深厚,土质肥沃,土质中含有丰富的腐殖质,其土壤pH为5.5~6.5,十分适合种农作物,除此之外,还较为适合种植反季节蔬菜、特色水果、中药材以及经济作物[1]。

端元提取算法matlab程序

端元提取算法matlab程序

端元提取算法matlab程序
端元提取算法matlab程序是一种基于光谱数据的数据处理方法,旨在从复杂的光谱数据中提取出反映物质本质特征的端元。

此程序采用多元统计学的方法,通过对光谱数据的降维和聚类处理,可以有效地提取出端元,并可用于化学、生物、医学等领域的数据分析与研究。

该程序具有以下特点:
1. 精度高,能够准确提取出端元数据。

2. 速度快,可处理大规模数据,提高数据处理效率。

3. 易于使用,操作简单,适合初学者使用。

4. 可扩展性强,可以根据不同领域和数据需求进行扩展和优化。

该程序的使用方法简单,只需输入原始光谱数据,选择端元提取算法和参数设置,即可得到提取出的端元数据。

该程序可用于科研工作和实际应用中,为各行各业的数据分析和研究提供了有效的工具和方法。

- 1 -。

高光谱数据处理基本操作规范

高光谱数据处理基本操作规范

高光谱分辨率遥感用很窄(10-2l)而连续的光谱通道对地物持续遥感成像的技术。

在可见光到短波红外波段其光谱分辨率高达纳米(nm)数量级,通常具有波段多的特点,光谱通道数多达数十甚至数百个以上,而且各光谱通道间往往是连续的,每个像元均可提取一条连续的光谱曲线,因此高光谱遥感又通常被称为成像光谱(ImagingSpectrometry)遥感。

行输出。

三、波谱库1、标准波谱库软件自带多种标准波谱库,单击波谱名称可以显示波谱信息。

2、自定义波谱库ENVI提供自定义波谱库功能,允许基于不同的波谱来源创建波谱库,波谱来源包括收集任意点波谱、ASCII文件、由ASD波谱仪获取的波谱文件、感兴趣区均值、波谱破面和曲线等等。

3、波谱库交互浏览波谱库浏览器提供很多的交互功能,包括设置波谱曲线的显示样式、添加注记、优化显示曲线等四、端元波谱提取端元的物理意义是指图像中具有相对固定光谱的特征地物类型,它实际上代(1)MNF变换了主成分分析在高光谱数据处理中的不足。

(2)计算纯净像元指数PPIPPI生成的结果是一副灰度的影像,DN值越大表明像元越纯。

作用及原理:纯净像元指数法对图像中的像素点进行反复迭代,可以在多光谱或者高光谱影像中寻找最“纯”的像元。

(通常基于MNF变换结果来进行)纯净像元指数可以将N维散点图映射为一个随机单位向量来计算,每次映射的极值像元被记录下来,并且被标为极值的总数也被记录下来。

按照多次映射每个像元被记录为极值像元的次数来决定该像元是否为纯净像元。

(3)端元波谱收集n维可视化工具-选取样本像元-生成地物平均波谱五、波谱识别和图像分类4、高光谱和多光谱实质上的差别就是,高光谱的波段较多,谱带较窄(比如hyperion有242个波段,带宽10nm);5、多光谱相对波段较少;6、高光谱遥感就是比多光谱遥感的光谱分辨率更高,但是光谱分辨率高的同时空间分辨率会降低。

优化端元提取方法的高光谱混合像元分解

优化端元提取方法的高光谱混合像元分解

优化端元提取方法的高光谱混合像元分解高光谱混合像元分解是一种常用的遥感数据分析方法,可以用于提取地物信息和监测环境变化。

在实际应用中,为了提高分解结果的准确性和可靠性,需要进行端元提取的优化。

端元提取是指从高光谱数据中选择代表地物的像元进行分解。

传统的端元提取方法主要基于经验或人工选择,存在以下问题:首先,传统方法需要人工选择代表地物样本进行端元提取,这种方式受主观因素干扰较大,容易引入误差。

为了减少主观因素的干扰,可以使用统计学方法来进行自动化的端元提取。

常用的统计学方法有聚类分析、主成分分析和最大似然分类等。

其次,传统方法在进行端元提取时通常只考虑了光谱信息,而忽略了空间信息。

然而,地物的空间分布特征对端元提取和混合像元分解结果的准确性和可靠性有重要影响。

因此,应该将空间信息考虑进来,可以利用地物边界信息和多源遥感数据进行融合,以提高端元提取的准确性。

此外,高光谱混合像元分解还需要考虑混合像元的数量和选择。

传统方法通常假设混合像元是由两个或三个端元组成的,但实际情况往往更为复杂,混合像元可能由多个端元组成。

因此,可以利用自适应光谱混合方法,对混合像元数量进行估计,并选择最优的混合像元组合来进行分解。

最后,在进行端元提取和混合像元分解时还应考虑光谱响应的非线性和光谱混叠的影响。

非线性效应会导致混合像元分解结果的偏差,光谱混叠则会造成端元提取的困难。

因此,可以利用非线性光谱混合像元分解方法和反混叠技术,来克服这些问题,提高分解结果的准确性。

综上所述,优化端元提取方法的高光谱混合像元分解可以使用统计学方法进行自动化的端元提取,同时考虑空间信息和光谱非线性效应等因素。

通过合理选择混合像元数量和采用反混叠技术,可以提高分解结果的准确性和可靠性,从而更好地应用于地物信息提取和环境监测等领域。

植被覆盖度反演

植被覆盖度反演

题目:植被遥感监测具有很长的历史,植被遥感中植被覆盖反演是主要内容之一,线性分解模型是混合像元分解法的一种,根据混合像元分解及线性分解模型的原理,利用红碱淖2016年6月17日Landsat8 OLI数据,求取研究区植被覆盖度(写出详细计算过程和步骤)。

一、操作思路:端元代表影像中的纯净像元,求取植被覆盖度即求取端元的丰度,因此采用混合像元分解及线性分解模型进行混合像元分解,得到端元以及相应的丰度影像,即为植被覆盖度影像。

二、操作步骤:1.影像预处理由于操作时间的限制和硬件设备的不允许,将影像选取一定的区域进行裁剪,感兴趣区域为roi1文件,得到裁剪后影像2016new文件,导入影像,进行影像的预处理,包括辐射定标和大气校正步骤,辐射定标采用Radiometric Calibration工具,大气校正采用FLAASH工具,分别得到辐射定标结果2016_rad1和2016_ref文件。

变换采用MNF变换工具,可将数据波段进行“降维”,提取出有用信息集中的波段,去除噪声信息集中的波段,利用Forward MNF Estimate Noise Statistic 将数据进行MNF变换,得到MNF变换结果2016_MNF和噪声文件MNF_Noise和统计文件MNF_Statistic文件。

变换纯净像元指数法指像元被标记为纯净像元的次数,可以将混合像元进行分解,有效的提取出端元。

由于操作时间的限制和硬件设备的不允许,将MNF变换后影像与原影像选取一定的同样区域进行裁剪,感兴趣区域选择ROI2文件,得到裁剪结果为2016_MNF_sub和2016new_sub文件。

在MNF变换后发现有用信息集中在1、2、3波段,因此利用Pixel PurityIndex 工具,将MNF变换后影像选择1、2、3波段进行端元提取,阈值设为,操作的结果为2016_PPI文件。

维可视化利用N维可视化工具可以将端元更好的显示,有利于更加直观的确定端元,在N维可视化窗口中,选中1、2、3波段,在显示窗口中将较为集中的区域定为端元,进行类(class)的划定,一共确定三类,利用mean all工具将三类端元的波谱显示出来,并保存为波谱库2016_sli文件。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

选取合适的端元是成功的混合像元分解的关键。

端元选取包括确定端元数量以及端元的光谱。

理论上,只要端元数量m小于等于b+1(b表示波段数),线性方程组就可以求解。

然而实际上由于端元波段间的相关性,选取过多的端元会导致分解结果更大的误差。

端元光谱的确定有两种方式:(1) 使用光谱仪在地面或实验室测量到的“参考端元”;(2) 在遥感图像上得到的“图像端元”。

方法(1)一般从标准波谱库选择,方法(2)直接从图像上寻找端元可选择的方法有:从二维散点图中基于几何顶点的端元提取,借助纯净像元指数(Pixel Purity Index——PPI)和n维可视化工具用于端元波谱收集,基于连续最大角凸锥(Sequential Maximum Angle Convex Cone——简称SMACC)的端元自动提取。

下面介绍几种端元选择的方法。

1基于几何顶点的端元提取将相关性很小的图像波段,如PCA、IC、MNF等变换结果的前面两个波段,作为X、Y轴构成二维散点图。

在理想情况下,散点图是三角形状,根据线性混合模型数学描述,纯净端元几何位置分布在三角形的三个顶点,而三角形内部的点则是这三个顶点的线性组合,也就是混合像元,如图所示。

根据这个原理,我们可以在二维散点图上选择端元波谱。

在实际的端元选择过程中,往往选择散点图周围凸出部分区域,后获取这个区域相应原图上的平均波谱作为端元波谱。

图散点图上的纯净像元与混合像元下面以MNF变换后的第一、第二波段作为X、Y轴构建二维散点图,如图所示。

图 Scatter Plot窗口2基于PPI的端元提取借助纯净像元指数(PPI)和n维可视化工具用于端元波谱收集,下面详细介绍操作步骤。

第一步、获取纯净像元这个步骤是在MNF变换的结果上计算纯净像元指数(PPI),之后选择阈值范围从PPI图像上获得感兴趣区,感兴趣区包含的像元就是比较纯净的像元。

(1)打开高光谱数据。

(2)在ENVI主菜单中,选择Spectral ->MNF Rotation- > Forward MNF -> Estimate Noise Statistics From Data。

在标准ENVI文件选择对话框中,选择高光谱图像文件。

打开Forward MNF Transform Parameters面板,选择MNF输出路径及文件名,单击OK执行MNF变换。

(3)在ENVI主菜单中,选择Spectral-> Pixel Purity Index->[FAST] New Output Band。

在打开的Pixel Purity Index Input File对话框中,选择MNF变换结果,单击Spectral Subset按钮,选择前面10个波段(MNF后面波段基本为噪声),单击OK。

(4)在Pixel Purity Index Parameters面板中,设置Threshold Factor:3,其他参数默认,选择输出路径及文件名,单击OK执行PPI计算。

(5)在Display窗口中显示PPI结果。

选择Overlay->Region of Interest,在ROI Tool 面板中,选择Options->Band Threshold to ROI,选择PPI图像作为输入波段,单击OK,打开Band Threshold to ROI面板(图)。

Min Thresh Value:10,Max Thresh Value:空(PPI图像最大值),其他默认设置,单击OK计算感兴趣区,得到的感兴趣区显示在Display窗口中。

图 Band Threshold to ROI面板第二步、构建n维可视化窗口(1)在ENVI主菜单中,选择Spectral ->n-Dimensional Visualizer,在n-D Visualizer Input File对话框中选择MNF变换结果,单击OK。

(2)在n-D Controls面板中,选择1、2、3、4、5波段,构建5维的散点图。

第三步:选择端元波谱(1)在n-D Controls面板中,设置适当的速度(Speed),单击Start 按钮,在n-D Visualizer窗口中的点云随机旋转,当在n-D Visualizer窗口中的点云有部分聚集在一块时,单击Stop按钮。

(2)在n-D Visualizer窗口中,用鼠标左键勾画“白点”集中区域,选择的点被标示颜色。

(3)在n-D Controls面板中,选择Class->Items 1:20->White(用于删除点),单击Start按钮,当看到有部分选择的点云分散时候,单击Stop按钮,在n-D Visualizer窗口中选择分散的点,自动会将选择的点删除。

借助<-,->,New按钮可以一帧帧从不同视角浏览以辅助删除分散点。

(4)在n-D Visualizer窗口中,单击右键选择New Class快捷菜单,重复(1)~(3)选择其他“白点”集中区域。

图 n-D Visualizer窗口中的端元第四步、输出端元波谱(1)在n-D Controls面板中,选择Options->Mean All,在Input File Associated with n-D Scatter Plot对话框中选择原图像,单击OK。

(2)获取的平均波谱曲线绘制在n_D Mean绘图窗口中。

(3)参考“波谱分析工具”章节,识别每条波谱曲线对应的地物类型。

(4)在n_D Mean绘图窗口中,选择File->Save Plot As->Spectral Library(或者ASCII),将端元波谱保存为波谱库文件或者文本文件。

3基于SMACC的端元提取连续最大角凸锥(Sequential Maximum Angle Convex Cone )简称SMACC。

SMACC方法可从图像中提取端元波谱以及丰度图像(abundance Image)。

它提供了更快,更自动化的方法来获取端元波谱,但是它的结果近似程度较高,精度较低。

SMACC方法是基于凸锥模型(也称为残余最小化)借助约束条件识别图像端元波谱。

采用极点来确定凸锥,并以此定义第一个端元波谱;然后,在现有锥体中应用一个具有约束条件的斜投影生成下一个端元波谱;继续增加锥体生成新的端元波谱。

重复这个过程直至生成的凸锥中包括了已有的终端单元(满足一定的容差),或者直至满足了指定的端元波谱类别个数。

通俗的解释,SMACC方法首先找到图像中最亮的像元,然后找到和最亮的像元差别最大的像元;继续再找到与前两种像素差别最大的像素。

重复该方法直至SMACC找到一个在前面查找像素过程已经找到的像素,或者端元波谱数量已经满足。

SMACC方法找到的像素波谱转成波谱库文件格式的端元波谱。

下面以一个高光谱数据为例,详细介绍这个工具的操作过程。

(1)在ENVI主菜单中,选择File->Open Image File,打开高光谱数据文件。

(2)在ENVI主菜单中,选择Spectral ->SMACC Endmember Extraction,在Select Input Image对话框中选择高光谱数据文件,单击OK打开SMACC Endmember Extraction Parameters面板(图)。

(3)在SMACC Endmember Extraction Parameters面板中,需要填写以下参数:端元波谱提取数量(Number of Endmembers):15误差容限值(RMS Error Tolerance):0默认值0表示只有达到Number of Endmembers 参数指定的终端个数,SMACC 才会结束。

如果指定一个RMS误差,那么达到这个RMS误差的话,SMACC就会结束,不管是否获取指定数量的端元波谱。

反射率数据推荐使用,辐射亮度值数据推荐使用1。

但是要注意反射率数据常常扩大了倍数,比如扩大了10000倍,这个时候RMS Error Tolerance参数设置应该为10000x1%=100。

选择分离端元波谱的约束条件(unmixing constraint For Endmember Abundances):Positivity Only:把每个波长的正值端元波谱作为约束条件。

这个选项常用于反射率数据,因为负反射率值没有物理意义,Sum to Unity or Less:等于或者小于每个像素计算得到每种物质的组分之和作为约束条件。

当想从反射率数据中获取物质的物理意义和丰度图像的阴影图时候,可以选择这个约束条件,结果中会单独生成一个丰度阴影图像(Shadow Abundance)。

Sum to Unity:等于每个像素计算得到每种物质的组分之和作为约束条件。

当零端元波谱没有物理意义或者想获得暗端元波谱可以选择这个约束条件,这个约束条件推荐用于辐射亮度数据和热辐射数据。

合并相似端元波谱(Coalesce Redundant Endmembers):该选项是基于波谱角制图方法把阈值(在SAM Coalesce Value对话框中定义的值)内的所有端元波谱合并为一个端元波谱。

如果想要区分波谱比较相似的地物,不要选择该选项。

输出结果文件Endmember Location ROIs:该输出包括从终端单元波谱结果中产生的像元感兴趣区文件,这个输出文件是可选的。

Abundance Image:输出丰度图像,该输出文件将包括阴影图像和终端单元聚集图像。

该输出图像是可选的。

Select Output Spectral Library Enter Output Filename:该输出文件中包括提取出的终端单元的波谱库信息。

这个是必先项目。

(4)单击OK,执行SMACC过程。

图 SMACC Endmember Extraction Parameters面板获取的端元波谱以ENVI波谱库文件形式保存,设置的端元数为15,由于设置合并相似端元波谱选项,实际获得6种端元波谱,借助Spectral Analyst功能识别获得的端元波谱。

同时还可以得到每种端元波谱的丰度图像。

相关文档
最新文档