典型光谱矿物识别.docx
《典型光谱矿物识别.docx》由会员分享，可在线阅读，更多相关《典型光谱矿物识别.docx（18页珍藏版）》请在冰豆网上搜索。
典型光谱矿物识别
【ENVI入门系列】16.基本光谱分析
（2014-09-3017:
38:
25）
转载▼
标签：
杂谈
分类：
ENVI
版权声明：
本教程涉及到的数据仅供练习使用，禁止用于商业用途。
基本光谱分析
1. 概述
2. 详细操作步骤
2.1 标准波谱库与浏览
2.2 波谱库创建
2.3 高光谱地物识别
2.3.1 从标准波谱库选择端元进行地物识别
2.3.2 自定义端元进行地物识别
1.概述
光学遥感技术的发展经历了：
全色（黑白）—>彩色摄影—>多光谱扫描成像—>高光谱遥感四个历程。
高光谱分辨率遥感（HyperspectralRemoteSensing）用很窄（小于10nm）而连续的光谱通道对地物持续遥感成像的技术。
在可见光到短波红外波段其光谱分辨率高达纳米（nm）数量级，通常具有波段多的特点，光谱通道数多达数十甚至数百个以上，而且各光谱通道间往往是连续的，因此高光谱遥感又通常被称为成像光谱（ImagingSpectrometry）遥感。
相比多光谱传感器，如Landsat8OLI包括9个波段，光谱分辨率为106nm。
如下图为从多光谱和高光谱数据上获取的波谱曲线，更高波谱分辨率的图像可以用于识别物质，而相比多光谱图像，如TM只能用于区分物质。
图：
从多光谱和高光谱数据上获取的波谱曲线对比
本课程学习ENVI的高光谱基本分析功能，包括波谱库的浏览与创建、基本的高光谱物质识别流程。
2.详细操作步骤
2.1标准波谱库与浏览
ENVI自带多种标准波谱库，包括建立在JPL波谱库基础上的，从0.4~2.5μm三种不同粒径160种"纯"矿物的波谱。
美国USGS从0.4~2.5μm包括近500种典型的矿物和一些植被波谱。
来自JohnsHopkinsUniversity（JHU）的波谱包含0.4~14μm。
IGCP246波谱库有5部分组成，通过对26个优质样品用5个不同的波谱仪测量获得。
植被波谱库由ChrisElvidge提供，范围是0.4~2.5μm。
ENVI5.1波谱库中新增了2443种Aster的波谱文件，同时对应的波谱工具也有了很大的该进，用户可直观地看到每一种波谱库中的文件个数，以及更为方便的查看每一种波谱文件的波谱曲线。
ENVI的波谱库文件存放在HOME\ProgramFiles\Exelis\ENVI51\classic\spec_lib。
启动ENVI5.1，主菜单>Display>SpectralLibraryView，在对话框中显示的就是ENVI自带的波谱库文件；
图：
ENVI自带波谱库文件
（1）选择打开Veg_lib（99）中的几个植被波谱文件；
在vegetation波谱库中选择6种不同植被的波谱曲线，在下图可以看到起对应的波谱曲线，以及波谱文件的属性信息，包括常规信息和曲线信息。
图：
不同植被波谱曲线
（2）波谱曲线显示窗口中的功能；
在波谱曲线窗口中可以看到4个部分的内容：
∙
导入文件（Import）
∙
图：
导入数据方式
此功能可以导入两种格式的文件，包括ASCII及波谱库形式存在的文件。
∙
导出文件（Export）
∙
图：
导出数据方式
导出波谱文件的格式可以分类4大类：
∙
常见数据格式——ASCII与波谱库文件；
∙
∙
Image、PDF及Postscript输出格式；
∙
∙
复制波谱曲线；
∙
∙
直接打印曲线或在PowerPoint中展示；
∙
∙
选项工具（Options）
∙
图：
选项功能
选项工具中有三个功能：
∙
打开新的Plot窗口——自由拖拽收集的地物波谱；
∙
∙
波谱曲线上显示十字丝——一直保持十字丝可见，显示波谱显现十字丝节点含义；
∙
∙
添加波谱图例——不同颜色的波谱曲线代表什么样的地物，更直观方便；
∙
∙
波谱曲线X、Y轴代表含义
∙
X轴代表：
∙
Wavelength：
（默认显示）影像波长；
∙
∙
Index：
波段i，i代表影像具有i个波段；
∙
∙
Wavenumber:
波数，即1/wavelength，波数与波长成反比关系，波长越小，波数就越大。
∙
Y轴代表：
∙
DataValue：
（默认显示）影像原始值；
∙
∙
包络线去除ContinuumRemoved：
绘制数据与连续删除。
连续的是，套在光谱的顶部的凸包。
它分为原始数据值，以产生连续取出的值。
在连续使用中的绘制曲线中所显示的第一和最后一个数据点计算的，所以只对已缩放的图形，在连续的基础上所显示的数据来计算范围。
∙
∙
BinaryEncoding：
二进制编码，重新生成0与1的波频曲线。
∙
图：
导出PNG格式的波谱曲线
（3）波谱曲线属性显示窗口
∙
：
同一窗口中显示多个地类的波谱曲线不予重叠显示；
∙
∙
：
恢复原始数值范围曲线显示；
∙
∙
：
点击此功能显示或者隐藏PlotKey与曲线属性；
∙
∙
：
移除选中的一个曲线数据；
∙
∙
：
移除全部的曲线数据；
∙
∙
：
如果曲线节点有异常可以通过此工具进行编辑修订；
∙
2.2波谱库创建
ENVI可以从波谱源中构建波谱库，波谱来源包括：
ASCII文件，由ASD波谱仪获取的波谱文件，标准波谱库，感兴趣区/矢量区域平均波谱曲线，波谱剖面和曲线等。
下面介绍波谱库建立的操作步骤：
第一步：
输入波长范围
（1）在ENVI中，选择Toolbox/Spectral/SpectralLibraries/SpectralLibraryBuilder。
打开SpectralLibraryBuilder对话框。
（2）为波谱库选择波长范围和FWHM值，有三个选项：
∙
"DataFile"（ENVI图像文件）：
波长和FWHM值（若存在）从选择文件的头文件中读取
∙
∙
"ASCIIFile"：
波长值与FWHM值的列的文本文件
∙
∙
"FirstInputSpectrum"：
以第一次输入波谱曲线的波长信息为准
∙
选择"FirstInputSpectrum"，单击ok，打开SpectralLibraryBuilder面板。
图：
SpectralLibraryBuilder对话框
第二步：
波谱收集
在SpectralLibraryBuilder面板中，可以从各种数据源中收集波谱，见下表。
所有收集的波谱被自动重采样到选择的波长空间。
表2.1波谱收集方法说明
菜单命令
功能
文本文件（FromASCIIfile）
从包含波谱曲线x轴和y轴信息的文本文件，当选择好文本文件时候，需要在InputASCIIFile面板中为x轴和y轴选择文本文件中相应的列。
当选择fromASCIIfile（previoustemplate）时，自动按照前面设置导入波谱信息。
fromASDBinaryFiles
从ASD波谱仪中导入波谱曲线。
波谱文件将被自动重采样以匹配波谱库中的设置。
当ASD文件的范围与输入波长的范围不匹配，将会产生一个全0结果。
fromSpectralLibrary
从标准波谱库中导入波谱曲线。
fromROI/EVFfrominputfile
从ROI或者矢量EVF导入波谱曲线，这些ROI/EVF关联相应的图像，波谱就是ROI/EVF上每个要素对应图像上的平均波谱。
fromStatsfile
从统计文件中导入波谱曲线，统计文件的均值波谱将被导入。
FromPlotWindows
从Pot窗口中导入波谱曲线。
下面介绍从高光谱图像数据中收集波谱。
启动ENVI，打开高光谱数据CupriteReflectance.dat。
∙
收集图像上某个像元的波谱：
∙
（1）选择Display>Profile>Spectral，在SpectralProfile对话框中，显示当前鼠标点的剖面曲线。
找到要收集的像元，鼠标选中，该像元的波谱曲线显示
（2）回到SpectralLibraryBuilder面板中，选择Options>FromPlotWindows，将所收集的波谱选中导入。
（3）导入的波谱显示在列表中，在波谱名称（SpectrumName）字段对应的记录双击鼠标以修改波谱名称。
同样的方法可以修改颜色（Color）字段的信息。
图：
ImportfromPlotWindows面板
这种方法是从图像上获取单个点的波谱曲线，也可以获取某个区域的平均波谱曲线，如ROI文件或者矢量文件。
∙
收集ROI或矢量文件范围的平均波谱
∙
（1）在ENVI主菜单中，选择File>Open，打开"ROI\CupriteMineralROIs.xml"感兴趣样本文件。
（2）回到SpectralLibraryBuilder面板中，选择Options>fromROI/EVFfrominputfile，选择高光谱文件CupriteReflectance.dat作为波谱来源。
图：
SpectralLibraryBuilder面板
（3）选中某一类感兴趣区，如方解石Calcite，点击Plot，绘制该感兴趣区的平均光谱曲线。
图：
方解石Calcite的平均光谱曲线
第三步：
保存波谱库
（1）SpectralLibraryBuilder面板中，点击SelectAll，将样本全部选中；
（2）SpectralLibraryBuilder面板中，选择File>SaveSpectraAs>SpectralLibrary，打开OutputSpectralLibrary面板；
（3）在OutputSpectralLibrary面板中，可以输入以下参数：
∙
Z剖面范围（ZPlotRange）:
空白（Y轴的范围，根据波谱值自动调节）
∙
∙
X轴标题（XAxisTitle）：
波长
∙
∙
Y轴标题（YAxisTitle）：
反射率
∙
∙
反射率缩放系数（ReflectanceScaleFactor）：
空白
∙
∙
波长单位（WavelengthUnits）：
Nanometers
∙
∙
X值缩放系数（XScaleFactor）：
1
∙
∙
Y值缩放系数（YScaleFactor）：
1
∙
（4）选择输入路径及文件名，单击OK保存波谱库文件。
图：
保存波谱库参数设置面板
2.3高光谱地物识别
高光谱图像分类方法与传统的多光谱分类有本质的区别，从高光谱图像的每个像元均可以获取一条连续的波谱曲线，就可以考虑用已知的波谱曲线和图上每个像元获取的波谱曲线进行对比，理想情况下两条波谱曲线一样，就能说明这个像元是哪种物质。
我们把高光谱图像分类、物质识别、探测等称为波谱识别。
ENVI提供
升级会员 