高光谱数据的制图方法简介
第五讲高光谱数据分析

第五讲高光谱数据分析高光谱技术可提供空间域信息和光谱域信息,即“谱图合一”,且由图像数据反演出的像元光谱曲线可以与实验室所测的同类地物光谱曲线相类比,因此可以用于鉴别物质,比如鉴别矿物、岩石的类型,区分环境中各种污染物的成分以及农作物、森林的种类等。
一、提取波谱剖面廓线打开文件cup95eff.int,并RGB彩色显示band 183,band 193,band 207。
在主影像菜单栏中选择Tools/Profiles/Z profiles(Spectrum),打开并显示波谱曲线。
二、采集波谱曲线在Spectral Profile窗口中,选择Options->Collect Spectra,采集绘图窗口中的波谱曲线。
相应的,要将波谱曲线采集到另一个绘图窗口中,先打开一个新的绘图窗口,然后将Spectral Profile窗口中的波谱曲线保存到新的绘图窗口中。
具体步骤如下:1.从绘图窗口的菜单栏中选择Options->New Window:blank,打开一个新的绘图窗口。
2.在先前的绘制窗口中,点击鼠标右键,选择Plot Key,将波谱曲线的名字显示在绘图窗口的右边。
3.在第一条波谱曲线的名字上,点击并按住鼠标左键不放,将波谱曲线的名字拖到新的绘图窗口中,然后松开鼠标左键。
4.在主影像窗口或缩放窗口中移动当前光标像素定位器,从影像中选择一条新的波谱曲线。
重复上面点击拖拽的过程,在新绘图窗口中建立一系列的波谱曲线。
要改变不同波谱曲线的颜色和线形,选择新绘图窗口中的Edit-Data Parameters.每一条波谱曲线的名字/位置都将在Data Parameters对话框中列出。
三、动画显示数据在先前的灰阶影像显示的主影像窗口中,选择Tools-Animation生成动画显示。
弹出的Animation Input Parameters对话框中列出了可用波段列表中的所有波段。
从所有波段中选择一个子集来生成动画。
高光谱ENVI使用方法简介

高光谱制图—FLAASH大气校正
FLAASH是目前精度最高的大气辐射校正模型, 使用了 MODTRAN 4+ 辐射传输模型的代码,基 于像素级的校正 FLAASH可对Landsat, SPOT, AVHRR, ASTER, MODIS, MERIS, AATSR, IRS等多光谱、高光谱 数据、航空影像及自定义格式的高光谱影像进行 快速大气校正分析。能有效消除大气和光照等因 素对地物反射的影响,获得地物较为准确的反射 率和辐射率、地表温度等真实物理模型参数
高光谱制图—FLAASH大气校正(5)
如果要自动保存前面所输入的FLAASH参 数 如果需要生成相关诊断文件(如通道定义 文件等)
高光谱影像地理坐标定位
空间遥感平台在传感器采集数据的同时也精确地 记录了自身的几何信息,使用这些几何信息如星 历、姿态数据以及传感器探元与成像数据上像元 间的几何关系等,可以计算出影像上每一个像元 所对应的经纬度,其结果将作为影像数据的辅助 地理信息一并打包发布给用户。利用这些详细的 输入几何信息(Input Geometry)使得影像不需 要选择大量地面控制点就可以进行几何精纠正, 即ENVI所谓的地理坐标定位Georeference)。
比较N维散点图和二维散点图 利用N维散点图进行端元选取,理解使用菜 单Class Controls的使用 N维可视化仪同光谱剖面的链接,使用鼠 标中键来进行光谱曲线的绘制 光谱分析与N维可视化仪连接起来
高光谱影像分析-光谱切面
光谱切面包括水平切面、垂直切面和任意 方向切面。 切面是一幅ENVI影像,沿水平方向的切面, 样本数等于光谱波段数,行数等于采样数; 沿垂直方向的切面,样本数等于行数;对 于任意方向的切面,样本数等于沿ROI折 线的像元总数
ENVI高光谱数据分析操作手册

感兴趣区和掩膜的选择和使用可具体情况具体分析,运行一项或两项均可。
北京卓立汉光仪器有限公司
4. 滤波
打开图像,FilterConvolutions and Morphology。在Convolutions and Morphology Tools 中,选择 Convolutions滤波类型(高通滤波 器、低通滤波 器、拉普拉斯算子、方向滤波器、高斯高通滤波器、高斯低通滤波器、中值滤波 器、Sobel、Roberts、自定义卷积核)。
2.3.2.3. 保存波谱库
北京卓立汉光仪器有限公司 在Spectral Libraries Resampling Parameters对话框中,为Resample Wavelength To选择匹配源,一般选择图像文件为参考。 输出重采样波谱库.sli
北京卓立汉光仪器有限公司
3. 感兴趣区和掩膜
3.1. 感兴趣区(ROI)
Display 窗 口 中 , Overlay → Region of Interest , 在 ROI 对 话 框 中 , 单 击 ROI_Type→Polygon. 绘制窗口中,选择Image,绘制一个多边形,右键结束,可根据需要多绘制 几个。
主菜单→Basic Tools→Subset Data via ROIs,选择裁剪图像。 在Saptial Subset via ROIs Parameters中,设置参数。 Select Input ROIs,选择绘制的ROI。 Mask Pixel Outside of ROIs选择yes。
4.1. 设置参数
Kernel Size(卷积核大小):奇数。 Image Add Back(加回值):将原始图像中的部分加回到卷积滤波结果图像中, Editable Kernel(卷积核中各项的值)。
高光谱数据的制图方法简介

高光谱数据的制图方法简介ENVI软件在Spectral菜单中提供许多波谱制图方法,包括:二进制编码、波谱角制图、线性波段预测(LS-Fit)、线性波谱分匹配滤波、混合调制匹配滤波、包络线去除,以及波谱特征拟合等。
本文主要介绍几种高光谱数据处理的过程操作。
1.二进制编码二进制编码分类技术根据波段值落在波谱均值的下方或上方,将数据和端元分别编码为0和1。
在编码过程中,使用一个高级的(exclusive)OR函数,用于将需要编码的数据波谱与参照波谱相比较,从而生成一幅分类图像。
选择菜单栏Spectral—Mapping Methods—Binary Encoding。
在打开的窗口设置参数如下:图1-1 二进制编码分类参数设置注意:“OutputRuleImages”切换按钮被设置为“No”,规则图像将不被保存。
分类结束后,规则图像将出现在可用波段列表中,可以在任何显示窗口中显示(或链接/覆盖),并可以使用ENVI的像元位置/值功能进行查询。
结果显示如图1-2:图1-2 原影像图(左)与二进制编码分类结果图(右)2. 波谱角分类波谱角分类(SAM)是一种基于自身的波谱分类方法,这种方法将图像波谱与参照波谱在N-维空间进行匹配。
SAM用到的参照端元波谱可以来自于ASCII文件、波谱库、统计文件或直接从图像中抽取(如ROI均值波谱),本实验中用的是ROI均值波谱。
SAM把端元波谱(被认为是一个N维向量,N维波段数)和像元向量放在n维空间中进行角度比较。
较小的角度代表象元与参照波谱匹配紧密。
这一技术用于数据定标时,对照度和反照率的影响并不敏感。
选择菜单栏Spectral—Mapping Methods—Spectral Angle Mapper。
设置参数如图2-1,波谱角分类结果,如图2-2:图2-1 波谱角分类参数设置图2-2 波谱角分类结果影像3.LS-Fit(线性波段预测)LS-Fit使用一个最小方框(least squares)拟合技术来进行线性波段预测。
高光谱图像处理技术的使用方法与技巧

高光谱图像处理技术的使用方法与技巧高光谱图像处理技术是一种在特定波长范围内连续获取大量光谱信息的技术。
它不仅可以提供丰富的光谱数据,还能提供高分辨率的空间信息,因此在许多领域都有广泛的应用。
本文将介绍高光谱图像处理技术的使用方法与一些常用的技巧。
首先,高光谱图像的处理流程主要包括预处理、特征提取和分类三个步骤。
预处理是为了去除图像中的噪声和杂质,使得后续的特征提取和分类工作更加准确。
常见的预处理方法包括影像校正、光谱校正和噪声去除等。
影像校正是为了解决图像中的光照不均匀问题,常用的方法有常规平滑和直方图匹配等。
常规平滑方法可以通过滤波算法去除图像中的噪声和高频分量,提高图像的可视性。
而直方图匹配则可以通过调整图像的亮度和对比度,使得不同图像之间的光照条件保持一致。
光谱校正是为了解决不同设备采集的高光谱数据存在光谱偏移的问题。
通常可以通过使用已知光谱的参考物质进行校正,如大气校正和地物光谱响应校正等。
大气校正可以去除大气对光谱数据的影响,使得数据更加准确。
地物光谱响应校正则是为了减少不同地物对光谱数据的影响。
噪声去除是为了去除因设备等原因造成的图像噪声,提高图像的质量。
常见的噪声去除方法包括中值滤波、高斯滤波和小波分析等。
中值滤波是一种基于排序统计的方法,通过对图像中的像素排序并取中值来去除噪声。
高斯滤波则是一种常用的线性滤波方法,通过将像素的值与周围像素的值按照一定的权重进行加权平均,得到滤波后的像素值。
小波分析是一种基于频域的方法,通过对图像进行频域分解和重构来去除噪声。
接下来是特征提取。
高光谱图像的特征提取是为了从原始数据中提取出与目标信息相关的特征。
常用的特征提取方法包括光谱特征提取、纹理特征提取和形状特征提取等。
光谱特征提取是通过对高光谱数据进行光谱统计分析来获得与目标信息相关的参数。
常用的统计参数包括均值、方差、偏度和峰度等。
这些参数可以反映出光谱数据在不同波段上的分布情况。
纹理特征提取是为了从高光谱图像中提取出纹理信息。
专题六:高光谱数据介绍

光谱库
• ENVI中 的相关应用-标准波谱库: *.sli ,*.hdr
高光谱图像
• 对地物进行多波段成像所得到的一组二 维图像,每个波段对应的一个二维图像。 • 高光谱图像与自然图像的区别在于多了 一维光谱信息。
• 图像立方体 Spectral>Build 3D Cube
光谱空间及光谱角
• 光谱曲线图是区分不同地物的 主要方法 • 区分大量光谱时,需要用光谱 空间来表述。 • 以n=2为例,光谱向量(右图) • 多维光谱空间在可视化绘图是 困难的,数学构建上是可能的。 • 光谱间的相似性可以通过光谱 向量间的角度来判断——光谱 角。
光谱端元
将相关性很小的图像波段,如PCA、MNF 的前两个波段,作为X,Y构成二维散点图。 在理想状态下,根据线性混合模型数学描 述,三角形顶端为纯净像元。
在实际选择中,往往选择凸出部分,再获 取这个区域相应的平均波谱。
主要流程
• • • • 查看高光谱图像 打开常见图谱库 端元波谱提取(MNF) 高光谱分类—波谱角(SAM)
专题六:高光谱遥感 hyperspectral remote sensing
遥感的发展趋势
• 平台、传感器——数据
– 高(空间)分辨率 – 高光谱 – 高时间、高辐射 – 遥感反演 – 面向对象 – 光谱端元
• 信息处理方法
“天地一体化”
• 应用方面
– 实用化、商业化、国际化、一体化
背景
• 随着对地观测技术的迅速发展,图像的光谱 分辨率、空间分辨率和时间分辨率有了较显 著的提高,高光谱和高空间分辨率图像得到 了越来越多的应用。 • 精细的观测地物——地物的识别,地物的成 分信息
• 高光谱图像由成像光谱仪产生 • 测谱学和遥感成像技术的融合。
高光谱数据可视化python实现

高光谱数据可视化python实现在遥感领域,高光谱数据处理是一项重要的任务。
高光谱数据能够提供丰富的光谱信息,但由于数据的复杂性,我们需要使用适当的工具来处理和可视化这些数据。
本文将介绍如何使用Python语言来实现高光谱数据的可视化。
1. 准备工作在开始之前,我们需要准备一些高光谱数据。
这些数据可以来自卫星、无人机或其他测量设备。
高光谱数据通常包含数百个波段,每个波段都代表了不同的光谱信息。
我们可以使用Python中的numpy库来处理这些数据,并使用matplotlib库进行可视化。
2. 数据加载与处理首先,我们需要将高光谱数据加载到Python环境中。
可以使用numpy库的loadtxt函数来导入数据文件。
假设我们的数据文件是以逗号分隔的文本文件,每一行代表一个像素点,每个像素点的光谱信息以逗号分隔。
加载数据的代码如下:```import numpy as npdata = np.loadtxt('data.csv', delimiter=',')```加载完数据后,我们可以使用numpy库的各种方法来对数据进行处理。
例如,可以计算每个波段的平均值、最大值、最小值等统计指标。
通过这些统计指标,我们可以更好地了解数据的分布情况,为后续的可视化提供参考。
3. 数据可视化接下来,我们可以使用matplotlib库进行高光谱数据的可视化。
matplotlib是一个功能强大的绘图库,可以绘制各种类型的图形,包括折线图、散点图、柱状图等。
首先,我们可以使用matplotlib的plot函数来绘制每个波段的光谱曲线。
代码如下:```import matplotlib.pyplot as plt# 绘制光谱曲线for i in range(data.shape[1]):plt.plot(data[:, i])plt.xlabel('Wavelength')plt.ylabel('Reflectance')plt.title('Spectral Curve')plt.show()```上述代码中,我们使用for循环遍历数据的每个波段,然后使用plot函数将每个波段的光谱曲线绘制出来。
高光谱数据分析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,动态显示各波段图像,能很快的分辨出包含信息量较多的波段。
高光谱数据格式

高光谱数据格式摘要:一、高光谱数据简介二、高光谱数据格式概述1.数据结构2.数据存储方式3.数据处理与分析方法三、常见高光谱数据格式介绍1.HDF5格式2.ENVI格式3.ASCII格式4.其他格式四、高光谱数据在实际应用中的案例与优势五、总结与展望正文:一、高光谱数据简介高光谱数据是一种特殊的遥感数据,它通过获取连续波段的光谱信息,为地表目标物识别、环境监测、资源勘探等领域提供了强大的技术支持。
高光谱数据具有光谱分辨率高、波段数量多、数据量大等特点,可以帮助我们从不同角度和深度挖掘地表特征信息。
二、高光谱数据格式概述1.数据结构高光谱数据通常包括两部分:一是光谱数据,即波段强度信息;二是与之相关的元数据,如波段名称、波段宽度、光谱采样间隔等。
2.数据存储方式高光谱数据的存储方式有多种,如HDF5、ENVI、ASCII等格式。
这些格式在存储数据时,既有单一波段的文件,也有多波段的文件。
3.数据处理与分析方法高光谱数据处理与分析方法包括预处理、特征提取、分类和聚类等。
预处理主要包括去除噪声、辐射校正、大气校正等;特征提取是从光谱数据中提取有用信息,如光谱指数、连续波段组合等;分类和聚类是对高光谱数据进行地物识别和分类。
三、常见高光谱数据格式介绍1.HDF5格式HDF5(Hierarchical Data Format 5)是一种高效、可扩展的文件格式,适用于存储大量数据。
HDF5文件具有灵活的结构,可以存储多波段的高光谱数据,同时支持复杂的数据类型和复杂的数组结构。
2.ENVI格式ENVI(Environmental Visualization Infrastructure)是一种专为遥感数据设计的数据格式。
ENVI文件结构清晰,易于阅读和分析,支持多波段高光谱数据存储。
此外,ENVI还提供了一系列图像处理和分析工具,便于高光谱数据的处理。
3.ASCII格式ASCII格式是一种简单的文本格式,适用于存储单波段的高光谱数据。
高光谱数据可视化python实现 -回复

高光谱数据可视化python实现-回复标题:基于Python的高光谱数据可视化实现导言:高光谱数据是一种含有大量连续波段信息的数据,它具有广泛的应用价值,如农业、环境监测、地质勘探等领域。
可视化这些数据能够帮助我们更好地理解和分析数据,从而发现隐藏在其中的规律和信息。
本文将介绍如何使用Python实现高光谱数据的可视化,通过一步一步的详细实践,让读者熟悉高光谱数据可视化的实现过程。
第一步:获取高光谱数据首先,我们需要获取一组高光谱数据。
可以通过公开的数据集或者自行采集所需数据。
数据集通常包含多个波段,并以像素为单位存储。
这里我们以一个包含8个波段的高光谱遥感图像数据集为例。
第二步:导入所需库和数据集在Python中,我们可以使用多个库来处理和可视化高光谱数据。
常用的库包括NumPy、Pandas、Matplotlib和Seaborn。
首先,我们要确保这些库已经安装在本地环境中。
然后,我们可以使用以下代码导入数据集并加载所需库:pythonimport numpy as npimport pandas as pdimport matplotlib.pyplot as pltimport seaborn as sns# 导入高光谱数据集dataset = pd.read_csv('spectral_data.csv')第三步:数据探索与预处理在进行数据可视化前,我们需要对数据进行探索和预处理。
这包括查看数据的基本信息,检查数据是否存在缺失值或异常值,并对数据进行必要的处理。
python# 查看前几行数据print(dataset.head())# 查看数据形状print(dataset.shape)# 检查是否存在缺失值print(dataset.isnull().sum())# 检查数据的统计摘要print(dataset.describe())根据探索结果,我们可以根据需要对数据进行清洗、处理缺失值或异常值等操作。
高光谱图像数据处理技术研究

高光谱图像数据处理技术研究
随着科技的不断发展,高光谱图像技术已经被广泛应用于许多领域,例如矿产勘探、农业、环境监测和陆地利用等,这种技术可以为人们提供丰富的信息和数据,帮助人们更好地了解和利用地球资源。
高光谱图像技术可以获取到比普通彩色图像更多的光谱信息,每个像素点可以对应到可见光谱到近红外谱段之间的多个波长,同时也拥有更高的空间分辨率。
因此,将高光谱图像数据处理起来是一项十分复杂的任务,需要大量的计算和分析。
在高光谱图像数据处理中,主要分为预处理、特征提取和分类识别三个部分。
其中,预处理主要是去除图像中的噪声、平滑图像以及校正图像等。
特征提取则是选取有用的光谱、空间和纹理特征,以便后续的分类识别。
而分类识别则是将图像分成不同的类别或者判断图像中是否存在某种物体等。
各种高光谱图像处理方法的发展,使得高光谱图像数据的处理能力不断提高,具有更广泛的应用前景。
以下是常用的高光谱图像处理方法:
1、传统的多光谱图像处理方法,是通过对光谱带的光谱反射率进行分析和处理,来实现图像的处理。
缺点是分辨率较低,不适合处理复杂图像。
2、基于小波变换的高光谱图像处理方法,通过将高光谱图像转换到小波空间中处理,可以实现更好的空间分辨率和高光谱特征提取能力。
3、基于深度学习的高光谱图像处理方法,通过构建深度学习模型,可以自动提取光谱、空间和纹理特征,同时也能够很好地处理高光谱图像中的噪声,取得了很好的效果。
综上所述,高光谱图像技术在数据处理方面具有广泛的应用前景,同时也需要我们不断研究和探索,来发现更多的应用场景和解决方案。
高光谱图像处理算法的使用教程与特征提取

高光谱图像处理算法的使用教程与特征提取近年来,随着高光谱技术的快速发展,高光谱图像处理算法在各个领域得到了广泛应用,尤其在地球观测、农业、遥感和医学影像等领域发挥着重要作用。
本文旨在介绍高光谱图像处理算法的使用教程,并重点讨论其中的特征提取方法。
高光谱图像是一种具有连续的光谱信息的多光谱图像,它在不同的波段中可以提供更加详细的信息,因此可以更准确地描述和识别目标。
在进行高光谱图像处理之前,首先需要了解图像处理的基本概念和技术。
一、高光谱图像处理基础知识1. 高光谱图像的表示与获取高光谱图像可以通过散射光谱仪或成像光谱仪获取。
它通常由一系列的波段组成,每个波段对应着特定的频率范围。
具体表示方法可以是一个三维数据集,其中两个维度代表图像的行和列,第三个维度表示不同的波段。
2. 高光谱图像的预处理在进行特征提取之前,通常需要对高光谱图像进行预处理,以消除噪声、纠正伪迹和增强图像的对比度等。
预处理的具体方法包括平滑滤波、去噪、边缘增强等。
3. 高光谱图像的分割与分类图像分割是指将图像分为具有相似特征的不同区域,而图像分类是指将像素或区域分别分配给不同的类别。
分割和分类是高光谱图像处理中常用的方法,它们可以帮助我们理解图像中的对象和区域。
二、高光谱图像处理算法的使用教程1. 特征提取特征提取是高光谱图像处理的核心任务之一,它旨在从高光谱数据中提取出能够最好地描述目标的特征。
特征可以是单波段的或多波段的,常用的特征包括光谱特征、空间特征和统计特征等。
光谱特征是指从每个像素的波段值中提取的特征,包括平均值、最大值、最小值、方差等。
空间特征是指在空间维度上从图像中提取的特征,常用的方法包括纹理特征、边缘特征等。
统计特征是指从整个图像或图像区域中提取的统计信息,如均值、标准差、直方图等。
2. 高光谱图像处理工具与库为了方便进行高光谱图像处理,有许多开源的工具和库可供使用。
其中,常用的包括ENVI、IDL、MATLAB等。
机载高光谱技术的特点 高光谱数据 点云建模 三维制图

机载高光谱技术的特点高光谱技术是对地观测的高尖端技术之一,在地学领域具有广泛的应用前景。
随着我国经济的高速增长,城市化水平的不断提高,资源枯竭、环境污染与破坏等问题也日趋严重,在这种条件下,高光谱技术将发挥重要的作用,并将带来显著的经济与社会效益。
目前,该技术已开始应用到地质与矿产资源调查、环境污染与生态检测和防治、土地资源和水资源利用与管理等领域。
高光谱技术的发展已促使地学研究的范围、尺度、内容和研究方法产生革命性的变化,其与GIS、 GPS技术的结合,为地学研究源源不断地提供高精度定位、高频度、多频谱不同级次的宏观影像,极大地拓宽了人类的视野和视觉能力,使人类真正能够从总体上把地球作为一个统一的系统加以分析,在不同级次上揭示地球各个圈层的相互联系和相互作用。
归纳起来,目前高光谱技术应用发展存在五大特点。
1)高光谱分辨率与高空间分辨率美国科学家通过对AVIRIS数据所进行的一系列研究,包括低空和高空试验飞行后认为,20m/pixel 的空间分辨率能够对岩石矿物进行一定程度的识别,但是对更为精细信息(如植株信息)的探测却较为困难。
同时,高光谱仪幅宽普遍较窄(含星载数据),以致重复飞行造成飞行成本较高。
因此随着高光谱应用技术的不断深入,以及对地物精微特征探测的需求,传感器更加趋向于朝高光谱分辨率、更高空间分辨率,以及更大幅宽方向发展。
20世纪90年代至21世纪初,高光谱技术已经由试验研究逐步走向商业化运营,其全球研究所需空间分辨率大致锁定在30m/pixel。
而军事和商业运营空间分辨率趋向几米甚至更小,数据获取幅宽也极大地提高,光谱分辨率也呈线性增长。
利用高光谱仪所获取的海量数据,可以提供更为精细的地表信息,实现采用宏观手段对地表微观特征进行识别与研究,使人们可以在不同领域分别得到相关的、丰富的地物内在信息,从而从物质的内在本质进行宏观把控、可持续管理等;这样也促使地学研究将空间尺度的信息与地学机理有机结合,实现遥感地学应用的综合与模型化。
hyperion高光谱数据矿物填图流程

hyperion高光谱数据矿物填图1.hyperion数据的简介EO-1是美国航空航天局为接替Landsat7而研制的新型地球观测卫星,于2000年11月21日发射升空,EO-1卫星轨道与Landsat7基本相同,为太阳同步轨道,轨道高度为705km,倾角98.7°,与Landsat7相隔50km,比Landsat7差1min过赤道。
EO-1卫星具有在当前轨道从各个方向获取相邻轨道重叠区域的影像的功能。
EO-1上搭载了3种传感器,即高级陆地成像仪ALI(Ad-vancedLandImager)、大气校正仪AC(Atmospheric Correc-tor)和高光谱成像光谱仪Hyperion。
Hyperion传感器以推扫的方式获取可见光、近红外(VNIR,400-1000nm)和短波红外(SWIR,900-2500nm)光谱数据。
Hyperion产品分两级:Level0和Level1,L0是原始数据,仅用来生成L1产品。
Hyperion是以推扫方式获取可见光-近红外(VNIR,400-1 000nm)和短波红外(SWIR,900-2 500nm)光谱数据。
Hyperion产品分两级:Level0和Level1,L0是原始数据,仅用来生产成L1产品。
L1产品有242个波段,1-70为可见近红外波段(VNIR),71-242为短波红外波段(SWIR),其中198个波段经过辐射定标处理,定标的波段分别为VNIR 8~57,SWIR77-224。
表1 Hyperion L1产品的特性说明2.hyperion数据的预处理流程图1 hyperion数据预处理流程图2.1.ENVI补丁下数据波段组合:Hyperion工具的用途:设计Hyperion工具组件是为了方便hyperion数据的使用。
其最基本的功能是把L1R HDF和L1G/L1T GeoTIFF数据集转换成包含波长、最大半波宽和坏波段信息的ENVI格式文件。
ENVI高光谱数据处理流程

高级设置
• 光谱定义文件:内置AVIRIS、HYMAP、HYDICE、HYPERION、 CASI、AISA。
• 气溶胶厚度系数:用于技术邻域效应范围。一般值为1~2km。 • CO2混合比率:默认为390ppm。 • 使用领域纠正。 • 使用以前的MODTRAN模型计算结果。 • 设置MODTRAN模型的光谱分辨率(推荐值5 cm-1)。 • 设置MODTRAN多散射模型。 • 提供三种模型供选择Isaacs,DISORT和 Scaled DISORT。默认是
和几个植被波谱
– Johns hopkins university 0.4~14um 矿物波谱
– IGCP264 (项目) 到
26个质优样品应用波谱仪测量得
• 打开波谱库(spectral/spectral libraries/…view)
• 显示波谱曲线(点击)
• 创建波谱库(spectral/spectral libraries/…builder)
• 两种气溶胶去除方法
– 2-Band(K-T)方法(类似模糊减少法),如果没有找到适应的黑 值(一般是阴影区或者水体),系统将采用能见度值来计算;所 以即使选择了该选项也要给能见度。
– 选择None,采用能见度值参与气溶胶去除,能见度值大约参考值 参见表
天气条件 晴朗 中等雾、阴霾 厚雾、阴霾
能见度
3、物质制图与识别、探测
•波谱库 •波谱分析 •端元波谱收集 •高光谱制图与识别
波谱库
• ENVI波谱库 (安装目录下spec_lib)
– Jet Propulsion Laboratory 0.4~2.5um 160种纯矿物波谱
处理高光谱影像和lidar数据的典型算法

处理高光谱影像和LiDAR数据的典型算法随着遥感技术的不断发展,高光谱影像和LiDAR数据已经成为遥感领域中的重要数据源。
这些数据具有丰富的信息量,可以用于地表覆盖分类、三维地形建模、环境监测等多个领域。
为了充分利用这些数据,研究人员开发了许多处理高光谱影像和LiDAR数据的典型算法。
本文将分别介绍处理高光谱影像和LiDAR数据的一些典型算法,并分析它们的特点和应用。
一、处理高光谱影像的典型算法1. 高光谱影像预处理算法高光谱影像预处理是处理高光谱影像数据的第一步,其主要目的是去除噪声、增强图像、降维等。
在高光谱影像预处理中,常用的算法包括主成分分析(PCA)、小波变换、均值滤波等。
这些算法能够有效地提取高光谱影像中的有效信息,为后续的分析和处理打下良好的基础。
2. 高光谱数据分类算法高光谱数据分类是遥感领域中的重要应用之一,其目的是根据高光谱数据对地物进行分类。
经典的高光谱数据分类方法包括支持向量机(SVM)、最大似然分类、随机森林等。
这些算法能够有效地对高光谱数据进行分类,实现对地物的识别和分类。
3. 高光谱数据融合算法高光谱数据融合是将高光谱数据与其他多源数据进行融合,获得更加丰富和准确的信息。
常用的高光谱数据融合算法包括基于像素的融合、基于特征的融合、基于决策的融合等。
这些算法能够有效地融合多源数据,为遥感应用提供更加丰富的信息源。
二、处理LiDAR数据的典型算法1. LiDAR数据滤波算法LiDAR数据滤波是处理LiDAR数据的关键步骤,其主要目的是去除噪声、提取地物信息等。
常用的LiDAR数据滤波算法包括统计滤波、高斯滤波、小波变换等。
这些算法能够有效地去除LiDAR数据中的噪声,提取出地物的准确信息。
2. LiDAR数据特征提取算法LiDAR数据特征提取是从原始LiDAR数据中提取出有用的地物特征信息,为后续的分析和应用提供支持。
常用的LiDAR数据特征提取算法包括高度提取、点云拟合、特征点提取等。
ENVI高光谱数据处理流程

ENVI高光谱数据处理流程一、显示图像波谱1.打开文件:主菜单中,File→Open Image File→文件名.raw或者Window→Available Bands List→File →Open Image File→文件。
2.显示真彩色图像:波段列表(Available Bands Lis)中,右键→Load TrueColor。
3.*设置像素大小:主窗口(Display)中,右键→Pixel Locator。
4.绘制波谱:主窗口中,右键→Z Profile(Spectrum)。
5.收集任意点波谱:Spectral Profile中,Options→Collect Spectra,点击图像任6.光谱平滑:Spectral Profile中,Options→Set Z Profile Avg Window,将window7.部分光谱:主菜单→Basic Tools→Resize Data(Spatial/Spectral)→Spectral Subset,选择需要的光谱波段。
生成新的文件,右键→Load True Color to<new>。
显示新图像。
8.关闭所有文件:File→Close All Files。
二、标准波谱库主菜单→Spectral→Spectral Libraries→Spectral Library Viewer→安装文件夹下,ITT\IDL\IDL80\products\envi48\spec_lib。
共有usgs_min、veg _lib、jpl_lib、jhu_lib四个标准波谱库。
在Spectral Library Viewer中,单击波谱名称,自动显示波谱。
三、自定义波谱库1.输入波长范围:在菜单中,Spectral Spectral Library→Spectral Library Builder2.波谱收集:以从影像数据中收集波谱为例:a)打开高光谱图像,收集任意点波谱。
高光谱数据分析ENVI操作手册

4.1. 设置参数
Kernel Size(卷积核大小):奇数。 Image Add Back(加回值):将原始图像中的部分加回到卷积滤波结果图像中, Editable Kernel(卷积核中各项的值)。
有助于保持图像的空间连续性。
滤波前
滤波后
北京卓立汉光仪器有限公司
5. 主成分分析列出各波段以及相应的百分比,可自主选择主成分波段。“No”系统会计 算特征值和显示供选择的输出波段。
5.2. 协方差矩阵、特征向量矩阵的统计
主菜单,Basic ToolsStatisticsView Statistics File,打开主成分分析中得到 的统计文件,可以得到各个波段的基本统计值、协方差矩阵、相关系数矩阵和特 征向量矩阵。 当协方差矩阵数据量较大时,不能直接在统计文件中显示,这时可通过输出 ASCII文件并导入到excel中来查看协方差矩阵和特征向量矩阵。 波长、 反射率和协方差矩阵、特征向量矩阵的数据分析可采用其他数值统计 分析软件进行。
2.2. 添加注记
在Spectral Library Plots窗口中,Option→Annotate Plot,手动添加注记,如文 Annotation窗口中,Object选择注记类型后,在Spectral Library Plots窗口中左 在Spectral Library Plots窗口中,右键→Plot Key,添加注记,名称和颜色在
选择Memory或在Enter Output Filename输入文件名生成新的文件。 右键→Load True Color to<new>,显示新图像。
北京卓立汉光仪器有限公司
1.6. 光谱数据输出
光谱曲线窗口中,File→Save Plot As→ASCII,在Output Plots to ASCII File文 件中,Selsct Plot To Output选中需要输出曲线的点,输出路径和名称,OK。
高光谱图像技术

光谱图像简介根据传感器的光谱分辨率对光谱成像技术进行分类, 光谱成像技术一般可分成3类。
光谱图像是指在特定波长范围内由一系列波长处的光学图像组成的图像块。
(1)多光谱成像———光谱分辨率在Δλ/λ= 0.1数量级,这样的传感器在可见光和近红外区域一般只有几个波段。
(2) 高光谱成像———光谱分辨率在Δλ/λ= 0.01数量级,这样的传感器在可见光和近红外区域有几十到数百个波段,光谱分辨率可达nm级。
(3) 超光谱成像———光谱分辨率在Δλ/λ= 0.001数量级,这样的传感器在可见光和近红外区域可达数千个波段。
高光谱图像基本原理高光谱图像技术是由高光谱遥感成像技术发展起来的一项技术。
高光谱图像是在特定波长范围内由一系列波长处的光学图像组成的三维图像块。
图为高光谱图像三维数据块的示意图。
图1中,x和y表示二维平面像素信息坐标轴,第三维(λ轴)是波长信息坐标轴。
从中可以看出,高光谱图像既具有某个特定波长λi下的图像信息,并且针对xy平面内某个特定像素,又具有不同波长下的光谱信息。
图1 高光谱图像数据块高光谱图像技术优势高光谱图像集样品的图像信息与光谱信息于一身。
由于光谱信息能充分反映样品内部的物理结构、化学成分,内部结构的差异可以通过特定波长下的光谱值来表现,在每个特定波长下,xy平面内每个像素点的灰度值又与其在该波长下的光谱值之间一一对应。
图像信息可以反映样品的形状、缺陷等外部品质特征,由于不同成分对光谱吸收亦不同的影响,在某个特定波长下图像对某个缺陷会有较显著的反映。
这些特点决定了高光谱图像技术在农产品内外部品质的检测方面的独特优势。
基于高光谱图像的无损检测系统基于高光谱图像的无损检测系统主要包括硬件平台(高光谱图像的获取) 和软件数据处理(高光谱图像数据分析) 两部分。
硬件平台主要由光源、分光部件、CCD、图像采集系统和计算机所组成。
根据分光部件的不同,其硬件平台又有2种不同的组建方式,即基于滤光片的高光谱图像获取系统和基于成像光谱仪的高光谱图像获取系统。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
高光谱数据的制图方法简介ENVI软件在Spectral菜单中提供许多波谱制图方法,包括:二进制编码、波谱角制图、线性波段预测(LS-Fit)、线性波谱分匹配滤波、混合调制匹配滤波、包络线去除,以及波谱特征拟合等。
本文主要介绍几种高光谱数据处理的过程操作。
1.二进制编码二进制编码分类技术根据波段值落在波谱均值的下方或上方,将数据和端元分别编码为0和1。
在编码过程中,使用一个高级的(exclusive)OR函数,用于将需要编码的数据波谱与参照波谱相比较,从而生成一幅分类图像。
选择菜单栏Spectral—Mapping Methods—Binary Encoding。
在打开的窗口设置参数如下:图1-1 二进制编码分类参数设置注意:“OutputRuleImages”切换按钮被设置为“No”,规则图像将不被保存。
分类结束后,规则图像将出现在可用波段列表中,可以在任何显示窗口中显示(或链接/覆盖),并可以使用ENVI的像元位置/值功能进行查询。
结果显示如图1-2:图1-2 原影像图(左)与二进制编码分类结果图(右)2. 波谱角分类波谱角分类(SAM)是一种基于自身的波谱分类方法,这种方法将图像波谱与参照波谱在N-维空间进行匹配。
SAM用到的参照端元波谱可以来自于ASCII文件、波谱库、统计文件或直接从图像中抽取(如ROI均值波谱),本实验中用的是ROI均值波谱。
SAM把端元波谱(被认为是一个N维向量,N维波段数)和像元向量放在n维空间中进行角度比较。
较小的角度代表象元与参照波谱匹配紧密。
这一技术用于数据定标时,对照度和反照率的影响并不敏感。
选择菜单栏Spectral—Mapping Methods—Spectral Angle Mapper。
设置参数如图2-1,波谱角分类结果,如图2-2:图2-1 波谱角分类参数设置图2-2 波谱角分类结果影像3.LS-Fit(线性波段预测)LS-Fit使用一个最小方框(least squares)拟合技术来进行线性波段预测。
它可以用于在数据集中找出异常波谱响应区。
LS-Fit先计算出输入数据的协方差,用它对所选的波段进行预测,预测值等于所用预测波段的线性组合再加上一个偏移值,在该过程中计算了实际波段和模拟波段之间的残差,并输出为一幅图像。
残差大的像元(无论正负)表示出现了不可预测的迹象(即一个吸收波段)。
模拟波段的图像也被包括在输出中。
被预测的波段可以根据现有统计或新的统计进行计算。
3.1用新的统计信息进行预测选择Spectral--Mapping Methods--LS-Fit(Linear Band Prediction)—Predict with new Statistics。
选择输入文件,根据需要选取空间子集。
将出现LS-Fit Parameters对话框。
设置参数如图3-1:图3-1 用新的统计信息进行预测的参数设置输出包含两个波段:模拟波段和残差图像。
残差图像中,值较大的像元,无论正负,都表明了所在位置的实际波段和模拟波段的差异。
预测结果如图3-2:图3-2 用新统计信息进行预测的模拟波段(左)和残差图像(右)3.2根据现有统计信息进行预测使用Predict from Existing Stats选项可以根据现有的统计文件对波段进行预测。
注意:只有在统计文件中包含协方差值的波段才能作为预测波段或模拟波段。
选择Spectral--Mapping Methods--LS-Fit(Linear Band Prediction)--Predict from Existing Stats。
选择输入文件,根据需要选取空间子集,选择与输入的数据文件相对应的统计文件,这里用到3.1中生成的统计文件。
在出现LS-Fit Parameters对话框,只有统计计算中用到的波段才可以被选为预测波段或模拟波段。
设置参数如图3-3:图3-3用现有的统计信息进行预测的参数设置输出包含的两个波段:模拟波段和残差图像。
残差图像中,值较大的像元,无论正负,都表明了所在位置的实际波段和模拟波段的差异。
预测结果如图3-4:图3-4 根据现有统计信息进行预测的模拟波段(左)和残差图像(右)由图3-2和图3-4得:两种预测方法结果类似,区别不大。
4.匹配滤波使用Matched Filtering选项可以进行局部分离,用于确定用户自定义端元的权重。
该方法并不需要图像中所有的端元都是已知的。
这项技术使已知端元的响应最大化,并抑制了未知背景合成的反应,因此“匹配”了已知信号。
它根据图像要素对波谱库或图像端元波谱的匹配程度,无需对图像中所有端元波谱进行了解,就可以快速探测出特定要素。
选择菜单栏Spectral--Mapping Methods--Matched Filtering。
当出现Matched Filter Input File selection 对话框时,选择输入文件,或根据需要选取任意空间和波谱子集或掩膜。
点击“OK”,将出现Endmember Collection:Matched Filter 对话框,选择相应的波谱,如图4-1图4-1 匹配滤波对应参数和端元选择匹配滤波的结果将以一系列灰阶图像的形式出现,每幅图像对应一个选择的端元。
浮点型结果提供了像元与参照波谱相对匹配程度的估计方法(1.0表示完全匹配),以及亚像元的权重。
图4-2 匹配滤波生成的一系列灰阶影像图5.混合调制匹配滤波使用Mixture Tuned Matched Filtering TM(MTMF )TM选项可以进行匹配滤波,混合调制匹配滤波需要输入经过MNF变换的文件或其它的isotropic数据,单位变化噪声。
MNF Rotation(Minimum Noise Fraction)工具用于判定图像数据内在的维数(即波段数),分离数据中的噪声,减少随后处理中的计算需求量。
Forward MNF变换用于估计第一次旋转中所用的噪声统计。
该选项包括:从输入的数据中估计噪声;运用以前计算的噪声统计;使用与数据集相关的“暗色图像”(dark Image)进行噪声统计。
本文应用“从输入的数据中估计噪声”的方法,选择菜单栏Transforms--MNF Rotation--Forward MNF--Estimate Noise Statistics From Data,选择要输入的文件以及子集,在出现的对话框保存统计文件:图5-1 Forward MNF Transform对话框保存设置ENVI处理完毕后,MNF波段将被导入到可用波段列表中,并显示MNF EigenValues 图表窗口如图5-2。
输出的MNF波段数等于选择输出的波段数,本实验选取了前50个波段。
特征值较大(大于1)的波段包含数据,特征值接近于1的波段包含噪声。
显示可用波段列表中的特征图像(MNF波段),并与MNF特征值图表相比较,可以判定出哪些波段包含数据,哪些波段主要包含噪声。
在随后的数据处理中,所选取的MNF波段的子集应该仅包括那些图像在空间上连续显示且特征值大于MNF图表中曲线陡坡转折处的特征值的波段。
图5-2 MNF EigenValues 图表窗口选择Spectral--Mapping Methods--Mixture Tuned Matched Filtering。
当出现Mixture Tuned Matched Filter Input File selection 对话框时,根据需要选择所需的MNF文件或其空间和波谱子集。
点击“OK”。
将出现Endmember Collection:Mixture Tuned Matched Filter 对话框。
输入要被匹配的波谱(MNF空间)。
当所需波谱选择完毕后,点击“Apply”。
将出现Mixture Tuned Matched Filter Parameters 对话框。
用箭头切换按钮选择“Compute New Covariance Stats”,选择输出到“Memory”或“File”。
在“Output Data Type”菜单中,选择输出数据类型:浮点型,点击“OK”,开始处理。
图5-3 混合调制匹配滤波结果图6.去除包络线包络线去除是将反射波谱标准化的一种方法,它允许从通用的基线对每个吸收特征进行比较。
包络线去除是一个在波谱顶部的凸起的外壳拟合,它用直线段连接局部的波谱最大值。
注意:使用不同的波谱子集将得到不同的结果,因此应当抽取包含吸收特征的感兴趣区作为子集。
通过将包络线区分为图像中每个像元的实际波谱而把它消除。
在最终图像中的包络线和匹配波谱处,波谱等于1.0,出现吸收特征的区域波谱小于1.0。
也可以对数据文件或绘图窗口中的单个波谱进行包络线去除。
选择菜单栏Spectral--Mapping Methods--Continuum Removal。
当出现Continuum Removal Input File对话框时,选择输入文件,若需要,选取任意子集或掩模。
为得到最好的结果,应当抽取包含吸收特征的感兴趣区作为子集。
点击“OK”。
当出现Continuum Removal Parameters对话框时,选择输出到“Memory”或“File”。
点击“OK”,开始处理。
结果如图7-1:图7-1 去包络线结果图(右)与原影像图(左)对比7.波谱特征拟合Spectral Feature Fitting TM(SFF TM)是一种基于吸收特征的方法,使用最小方块技术将图像波谱的拟合比选择的参照波谱。
在对数据集进行包络线去除以后,参照波谱被缩放,从而与图像波谱相匹配。
该方法将为每一个参照波谱输出一幅比例图像,它可以对与要素权重相关的吸收特征的强度进行度量。
波谱特征拟合在一个小方块范围内,在每个选择的波长处对图像和参照波谱进行比较,并为每个参照波谱评定平方根误差。
选择菜单栏Spectral--Mapping Methods--Spectral Feature Fitting。
当出现Spectral Feature Fitting Input File对话框时,选择经过包络线去除的输入文件,若需要,选取任意子集或掩膜。
为得到最好的结果,应当抽取包含吸收特征的感兴趣区作为子集。
在出现的Endmember Collection: Feature Fitting对话框中输入所需的参照波谱。
点击“Apply”。
当出现Spectral Feature Fitting Parameters对话框时,用箭头切换按钮选择“Output separate Scale and RMS Images”或“Output Combined (Scale/RMS) Image”,如图7-1,并选择输出位置。