ENVI高光谱数据处理流程

ENVI高光谱分析技术主要内容1、高光谱简介2、高光谱数据预处理3、物质制图与识别、探测4、植被分析1、高光谱遥感简介光学遥感技术的发展：全色（黑白）－－彩色摄影－－多光谱扫描成像－－高光谱遥感高光谱分辨率遥感（HyperspectralRemoteSensing）用很窄（10-2λ）而连续的光谱通道对地物持续遥感成像的技术。在可见光到短波红外波段其光谱分辨率高达纳米(nm)数量级，通常具有波段多的特点，光谱通道数多达数十甚至数百个以上，而且各光谱通道间往往是连续的，因此高光谱遥感又通常被称为成像光谱（ImagingSpectrometry）遥感。从光谱影像上获得光谱曲线高光谱图像空间成像的同时， 记录 混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载 下成百个连续光谱通道数据从每个像元均可提取一条连续的光谱曲线对高光谱图像的处理实质是对像元光谱曲线的定量化处理与分析高光谱成像技术成像光谱仪：与地面光谱辐射计相比，成像光谱仪不是在“点”上的光谱测量，而是在连续空间上进行光谱测量，因此它是光谱成像的；与传统多光谱遥感相比，其光谱通道不是离散而是连续的，因此从它的每个像元均能提取一条平滑而完整的光谱曲线。成像光谱仪系统介绍航空成像光谱仪系统国内系统：MAIS、OMIS-1、OMIS-2、PHI、WHI、LASIS国外系统：AIS、AVIRIS、TRWIS、GERIS、HYDICEAISA、DAIS、CASI、HYMAPAVIRISSpectralcoverage:VIStoNIR(400-2500nm)Spectralbands:224Spectralresolution:<10nmFOV:30°IFOV:1.0mradDigitization:12bitsHYMAPSpectralcoverage:VIS:400-800nm,15nmbands;NIR:881-1335nm,14nmbands;SWIR1:1400-1813nm,12nmbands;SWIR2:1950-2543nm,16nmbands;Spectralbands:126FOV:60°IFOV:2.5mrad(along_track)2.0mrad(across_track)Pixelsperline:512航天成像光谱仪系统——Hyperion/EO-1国家：美国时间：2000年11月卫星发射成功扫描带宽：7.5km,空间分辨率：30米，在0.4－2.5μm共有220波段:可见光－近红外(400-1000nm):60波段，短波红外(900-2500nm):160波段。环境与减灾小卫星星座（HJ-1B）2、高光谱数据预处理传感器定标大气校正传感器定标传感器定标是针对设备本身，建立传感器每个探测元件输出的数据量化值（DN）与它所对应像元内的实际地物的辐射亮度之间的定量关系（陈述彭等，1998）。辐射亮度（辐射率）单位可为：（μW）/（cm2*nm*sr）。ENVI提供针对特定传感器的定标，包括ASTER、AVHRR、MODIS、MSS、TM、IKONOS、QuickBird、WorldView等；通用方法，包括：平场域定标、对数残差、内部平均反射率法和经验线性；针对热红外数据，还提供大气校正工具、相对通道发射率、归一化发射率、Α残差等定标工具。为什么做大气纠正？太阳辐射通过大气以某种方式入射到物体表面然后再反射回传感器原始影像包含物体表面，大气，以及太阳的信息如果我们想要了解某一物体表面的光谱属性，我们必须将它的反射信息从大气和太阳的信息中分离出来。13IwouldliketotalkhereatthebeginningforacoupleofminutesaboutwhyAtmosphericCorrectionisoftenimportantforhyperspectralimageanalysis.Thebasicproblemwithimagesoftheearth’ssurfacethatusethesunasanilluminationsourceis...大气散射直接反射邻接反射14Sameinformationinpictures大气校正方法基于辐射传输模型LOWTRAN模型MORTRAN模型ATCOR模型6S模型基于统计学模型平场域定标对数残差内部平均反射率法经验线性基于简化辐射传输模型的黑暗像元法基于统计的不变目标法基于植被指数的大气阻抗植被指数法……ENVI大气校正模块ENVI的大气校正模块的模型为MODTRAN4+模型，它是由SpectralSciences,Inc.(SSI)和AirForceResearchLabs(AFRL)合作开发，ITTVIS进行整合和图形化。ENVI的大气校正模块能够对高光谱、多光谱影像进行校正。高光谱包括：HyMAP、AVIRIS、HYDICE、HYPERION、Probe-1,CASI、AISA等；多光谱包括：ASTER、AVHRR、IKONOS、IRS、Landsat、MODIS、SeaWiFS、SPOT、QuickBird等，以及航空（860nm-1135nm）数据。多光谱与高光谱的模型基础一样：MODTRAN4+。这个模块通过高光谱像素光谱上的特征来估计大气的属性，可以有效地去除水蒸气,气溶胶散射，漫反射的邻域效应。采用向导式操作 流程 快递问题件怎么处理流程河南自建厂房流程下载关于规范招聘需求审批流程制作流程表下载邮件下载流程设计 ，还包括快速大气校正功能。使用ENVI大气校正模块——输入文件准备数据是经过定标后的辐射亮度（辐射率）数据数据带有中心波长（wavelenth）值，如果是高光谱还必须有波段宽度（FWHM）,这两个参数都可以通过编辑头文件信息输入（EditHeader）。数据类型支持四种数据类型：浮点型（floating）、4-bytesignedintegers,2-bytesignedintegers,以及2-byteunsignedintegers。数据存储类型：ENVI 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 栅格格式文件，且是BIP或者BIL。波谱范围：flaash能够做的数据光谱范围是0.4－2500μm。使用ENVI大气校正模块——基本参数设置传感器基本信息设置使用ENVI大气校正模块——大气模型 ModelAtmosphere WaterVapor(stdatm-cm) WaterVapor(g/cm2) SurfaceAirTemperature Sub-ArcticWinter(SAW) 518 0.42 -16°Cor3°F Mid-LatitudeWinter(MLW) 1060 0.85 -1°Cor30°F U.S.Standard(US) 1762 1.42 15°Cor59° Sub-ArcticSummer(SAS) 2589 2.08 14°Cor57° Mid-LatitudeSummer(MLS) 3636 2.92 21°Cor70° Tropical(T) 5119 4.11 27°Cor80° Latitude(°N) Jan. March May July Sept. Nov. 80 SAW SAW SAW MLW MLW SAW 70 SAW SAW MLW MLW MLW SAW 60 MLW MLW MLW SAS SAS MLW 50 MLW MLW SAS SAS SAS SAS 40 SAS SAS SAS MLS MLS SAS 30 MLS MLS MLS T T MLS 20 T T T T T T 10 T T T T T T 0 T T T T T T -10 T T T T T T -20 T T T MLS MLS T -30 MLS MLS MLS MLS MLS MLS -40 SAS SAS SAS SAS SAS SAS -50 SAS SAS SAS MLW MLW SAS -60 MLW MLW MLW MLW MLW MLW -70 MLW MLW MLW MLW MLW MLW -80 MLW MLW MLW SAW MLW MLW水气去除设置WaterRetrieval水气去除设置，采用两种方式对水气进行去除：利用水气去除模型恢复影像中每个像元的水气量使用水气去除模型，数据必须具有15nm以上波谱分辨率，且至少覆盖以下波谱范围之一：1050-1210nm(对应1135nm)870-1020nm(对应940nm)770-870nm(对应820nm)单一的水气因数用于整体影像，默认是1，对于多光谱数据使用水气去除模型，可以在多光谱设置中手动设置水气波段气溶胶模型（AerosolModel）提供四种标准MODTRAN气溶胶模型Rural（乡村）、Urban（城市）、Maritime（海洋）、Tropospheric（对流层，能见度在40km以上）两种气溶胶去除方法2-Band（K-T）方法(类似模糊减少法)，如果没有找到适应的黑值（一般是阴影区或者水体），系统将采用能见度值来计算；所以即使选择了该选项也要给能见度。选择None，采用能见度值参与气溶胶去除，能见度值大约参考值参见表 天气条件 能见度 晴朗 40to100km 中等雾、阴霾 20to30km 厚雾、阴霾 15km或者更少高级设置光谱定义文件：内置AVIRIS、HYMAP、HYDICE、HYPERION、CASI、AISA。气溶胶厚度系数：用于技术邻域效应范围。一般值为1~2km。CO2混合比率：默认为390ppm。使用领域纠正。使用以前的MODTRAN模型计算结果。设置MODTRAN模型的光谱分辨率（推荐值5cm-1)。设置MODTRAN多散射模型。提供三种模型供选择Isaacs，DISORT和ScaledDISORT。默认是ScaledDISORT和streams为8，这种模型对于小于1000nm具有较高的精度；天顶角\方位角（针对非星下点传感器）。输出反射率缩放系数（OutputReflectanceScaleFactor）：为了降低结果储存空间，默认反射率乘于10000。输出结果表面反射率影像水气含量数据云图日志文件FLAASH大气校正模板文件3、物质制图与识别、探测波谱库波谱分析端元波谱收集高光谱制图与识别波谱库ENVI波谱库（安装 目录 工贸企业有限空间作业目录特种设备作业人员作业种类与目录特种设备作业人员目录1类医疗器械目录高值医用耗材参考目录 下spec_lib）JetPropulsionLaboratory0.4~2.5um160种纯矿物波谱美国地质调查局(USGS)0.4~2.5um500种质优矿物波谱和几个植被波谱Johnshopkinsuniversity0.4~14um矿物波谱IGCP264(项目)　26个质优样品应用波谱仪测量得到打开波谱库（spectral/spectrallibraries/…view)显示波谱曲线(点击）创建波谱库（spectral/spectrallibraries/…builder)波谱库的创建与浏览输入波长范围输入光谱从图像中获取外部文件（二进制）导入ASD波谱仪波谱库交互波谱库查看、编辑和分析波谱分割波谱重采样波谱分析工具SpectralAnalyst波谱分析首先需要打开一个波谱库，然后将未知波谱与波谱库中的波谱进行匹配处理，该工具运用波谱角分类,波谱特征拟和二进制编码技术,对一未知波谱与波谱库中要素的匹配进行排序,输出一个列表,按照波谱匹配的好坏依次排列,并纪录一个总体的得分.匹配时需要设置三种方法所占的权重,权重是任意的,最后输出一个总体得分,得分越高,表明匹配效果越好.光谱识别流程MNF变换重要作用用于判定图像内在的维数分离数据中的噪声减少计算量弥补了主成分分析在高光谱数据处理中的不足计算时需要输入的参数统计信息的图像范围　shiftdiffsubset噪声统计文件（可以用到另一副图像上做变换）MNF统计文件(反变换的时候要用)Mnf变换输出波段选择(根据特征值选择输出波段)利用MNF变换后的图像进行端元提取经过MNF变换后的图像波段之间的相关性有效地降低了，并且出现了多个拐点，这些拐点就是我们要找的端元—Endmembers。波谱曲线对比编辑端元生成端元地物平均波谱曲线利用波谱分析工具进行端元识别\识别每一类地物后,将端元导出成ROIPPI纯净像元指数生成的结果是一副灰度的影像,DN值越大表明像元越纯.作用及原理纯净像元指数法对图像中的像素点进行反复迭代，可以在多光谱或者高光谱影像中寻找最”纯”的像元.(通常基于MNF变换结果来进行)纯净像元指数可以将N维散点图迭代映射为一个随机单位向量来计算,每次映射的极值像元被记录下来,并且被标为极值的总数也被记录下来.按照多次映射每个像元被记录为极值像元的次数来决定该像元是否为纯净像元.计算时需要输入的参数进行迭代的次数设置域值的系数(极值像元的域值)数据二次采样(减少内存,但不能太小)处理时会出现纯净像元指数图X-像元迭代次数,Y-极值像元总数该图显示了PPI做为迭代次数的函数当所有极值像元都已经找到后,它应该接近于一条水平线处理结果:DN值越大表明像元越纯.PPI纯净像元指数提取ROI设置波段域值提取PPI_ROI通常只设最小值(值越大,像元越纯净)N维可视化器提取PPI_ROI后打开N维可视化器,选择MNF变换结果,系统默认导入PPI_ROIPPI_ROI端元提取波谱曲线对比编辑ROI生成地物平均波谱波谱分析,端元识别端元波谱收集端元波谱收集器ASCII文件ASD测量文件波谱库ROI/矢量文件统计文件剖面波谱图N维可视化分析器二维散点图SMACC工具所有监督分类与高光谱制图工具端元波谱收集的作用和多光谱中的训练样本的收集一样的。提供的途径非常多，比如可以从波谱库中获取，也可以通过波谱仪测量获得。34高光谱制图二进制编码波谱角分类线性波段预测（LS-Fit）线性波谱分离光谱信息散度匹配滤波混合调谐匹配滤波（MTMF）包络线去除光谱特征拟合多范围光谱特征拟合高光谱中的分类叫制图，Envi提供的高光谱制图方法很多。35高光谱制图练习从标准波谱库中获取端元波谱打开jpl_lib/jpl1.sli波谱库文件，获得以下物质的波谱ALUNITESO-4ABUDDINGTONITEFELDSTS-11ACALCITEC-3DCHABAZITETS-15AILLITEPS-11AKAOLINITEWELLORDEREDPS-1A选择SpectralAngleMapper制图方法波谱沙漏工具——目标识别向导式波谱分析工具影像亮度值定标为反射率最小噪声分离（MNF）纯净像元指数（PPI）N维散度分析选择终端单元地物制图（地物识别）ApparentReflectanceMNFPPIn-DIDMapDistributionandAbundance作用：将地物波谱与端元波谱进行对比来地物制图（识别地物）高光谱处理流程比较专业一些，这里提供了波谱沙漏工具，一步步引导你完成地图制图或者是地图识别。37SAM与BandMax结合进行目标探测基于BandMax向导的SAM目标探测引导你完成高光谱影像的目标探测。向导的BandMax部分能帮你找到最佳的波谱子集以区分背景和目标，并节省处理的时间。向导主要有以下几个步骤：选择输入/输出文件选择目标—选择目标波谱选择背景—选择需要抑制的背景信息利用BandMax计算有效波段—识别对SAM分析中有效的波段选择最大角阈值—定义SAM最大角检验制图成果—SAM分析以及成果检验。4、植被分析植被指数计算器植被分析植被指数计算器提供了6类27种植被指数的计算绿度　Greenness　光利用率　LightUseEfficiency氮干旱或炭衰减DryorSenescentCarbonStressPigments冠层水分含量CanopyWaterContent提供了生物物理学交叉检验通过植被指数计算器统一进行计算提供了每种植被指数的详细资料和计算公式能够根据影像信息自动显示可计算的植被指数没有一种其它的商业软件提供如此丰富的植被指数计算40植被对于很多领域都非常重要，ENVI提供了非常全的植被指数计算工具，它还可以根据您提供的数据源的情况，自动列出你能获得的植被指数，植被分析工具农作物胁迫AgriculturalStress易燃性分布分析FireFuels林木健康分析ForestHealth植被的抑制41还提供了植被分析工具，包括植被的胁迫性分析、着火性分析、林木健康情况和植被的抑制，其中植被的抑制分析对于地质上分析地质地貌作用非常大。农作物胁迫AgriculturalStress创建农作物胁迫的空间分布图判断出适合农作物生长的区域，可用于精确农业分析；绿度：来标识闲置农田、虚弱的植被、健康的作物；光利用效率：标识植被生长率；冠层氮含量：用于估计相关的氮等级；叶绿素：突出植被胁迫；冠层水分含量：标识水胁迫的层次。易燃性分布分析创建某一区域植被易燃性的空间分布图用于森林规划，也可用它来分析城郊混合区的火灾风险绿度：显示大多数的绿色指数；冠层水分含量：标识植被的冠层水分含量；干旱或炭衰竭：突出非生物指数；林木健康分析工具创建整个森林区域健康程度的空间分布图用于检测病虫害以及枯萎病的发生情况，也可以用于评估某地区的木材收获量绿度：表面绿色植被的分布；叶绿素：标识类胡萝卜素以及花青素的含量；冠层水分含量：标识水分含量；光使用效率：标识森林生长率；植被抑制从多光谱和高光谱减少或移除植被光谱信息利用红外与近红外波段方便地质或城市解译��影像文件最小噪声分离MNF数据维数判断是否从图像获得端元波谱计算纯净像元指数用户选定端元波谱N维可视化和端元选择是否输入用户选定端元波谱波谱识别结果否是