ENVI高光谱分析技术

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,ENVI高光谱分析技术,主要内容,1、高光谱简介,2、高光谱数据预处理,3、物质制图与识别、探测,4、植被分析,1、高光谱遥感简介,光学遥感技术的发展：,全色（黑白）彩色摄影多光谱扫描成像高光谱遥感,高光谱分辨率遥感（,HyperspectralRemote Sensing,）,用很窄（,10-2,）而连续的光谱通道对地物持续遥感成像的技术。在可见光到短波红外波段其光谱分辨率高达纳米,(nm),数量级，通常具有波段多的特点，光谱通道数多达数十甚至数百个以上，而且各光谱通道间往往是连续的，因此高光谱遥感又通常被称为成像光谱（,Imaging Spectrometry,）遥感。,从光谱影像上获得光谱曲线,高光谱图像,空间成像的同时，记录下成百个连续光谱通道数据,从每个像元均可提取一条连续的光谱曲线,对高光谱图像的处理实质是对像元光谱曲线的定量化处理与分析,高光谱成像技术,成像光谱仪：,与地面光谱辐射计相比，成像光谱仪不是在“点”上的光谱测量，而是在连续空间上进行光谱测量，因此它是光谱成像的；,与传统多光谱遥感相比，其光谱通道不是离散而是连续的，因此从它的每个像元均能提取一条平滑而完整的光谱曲线。,成像光谱仪系统介绍,航空成像光谱仪系统,国内系统：MAIS、OMIS-1、OMIS-2、PHI、WHI、LASIS,国外系统：AIS、AVIRIS、TRWIS、GERIS、HYDICEAISA、DAIS、CASI、HYMAP,AVIRIS,Spectral coverage:,VIS to NIR (400-2500nm),Spectral bands: 224,Spectral resolution: 10nm,FOV: 30,IFOV: 1.0 mrad,Digitization:12 bits,HYMAP,Spectral coverage:,VIS:400-800nm,15nm bands;,NIR:881-1335nm, 14nm bands;,SWIR1:1400-1813nm, 12nm bands;,SWIR2:1950-2543nm, 16nm bands;,Spectral bands: 126,FOV: 60,IFOV: 2.5 mrad(along_track),2.0 mrad(across_track),Pixels per line: 512,航天成像光谱仪系统Hyperion/EO-1,国家：美国,时间：2000年11月卫星发射成功,扫描带宽：7.5km, 空间分辨率：30米，,在0.42.5m共有220波段:,可见光近红外(400-1000nm): 60波段，,短波红外(900-2500nm): 160波段。,环境与减灾小卫星星座（HJ-1B）,2、高光谱数据预处理,传感器定标,大气校正,传感器定标,传感器定标是针对设备本身，建立传感器每个探测元件输出的数据量化值（,DN,）与它所对应像元内的实际地物的辐射亮度之间的定量关系（陈述彭等，,1998,）。辐射亮度（辐射率）单位可为：（,W,）,/,（,cm,2,*nm*sr,）。,ENVI,提供针对特定传感器的定标，包括,ASTER,、,AVHRR,、,MODIS,、,MSS,、,TM,、,IKONOS,、,QuickBird,、,WorldView,等；通用方法，包括：平场域定标、对数残差、内部平均反射率法和经验线性；针对热红外数据，还提供大气校正工具、相对通道发射率、归一化发射率、,残差等定标工具。,为什么做大气纠正？,太阳辐射通过大气以某种方式入射到物体表面然后再反射回传感器,原始影像包含物体表面，大气，以及太阳的信息,如果我们想要了解某一物体表面的光谱属性，我们必须将它的反射信息从大气和太阳的信息中分离出来。,大气散射,直接反射,邻接反射,大气校正方法,基于辐射传输模型,LOWTRAN,模型,MORTRAN,模型,ATCOR,模型,6S,模型,基于统计学模型,平场域定标,对数残差,内部平均反射率法,经验线性,基于简化辐射传输模型的黑暗像元法,基于统计的不变目标法,基于植被指数的大气阻抗植被指数法,ENVI大气校正模块,ENVI,的大气校正模块的模型为,MODTRAN 4+,模型，它是由,Spectral Sciences, Inc. (SSI),和,Air Force Research Labs (AFRL),合作开发，,ITT VIS,进行整合和图形化。,ENVI,的大气校正模块能够对高光谱、多光谱影像进行校正。,高光谱包括：,HyMAP,、,AVIRIS,、,HYDICE,、,HYPERION,、,Probe-1, CASI,、,AISA,等；,多光谱包括：,ASTER,、,AVHRR,、,IKONOS,、,IRS,、,Landsat,、,MODIS,、,SeaWiFS,、,SPOT,、,QuickBird,等，以及航空（,860nm-1135nm,）数据。,多光谱与高光谱的模型基础一样：,MODTRAN 4+,。这个模块通过高光谱像素光谱上的特征来估计大气的属性，可以有效地去除水蒸气,气溶胶散射，漫反射的邻域效应。采用向导式操作流程，还包括快速大气校正功能。,使用,ENVI,大气校正模块,输入文件准备,数据是经过定标后的辐射亮度（辐射率）数据,数据带有中心波长（wavelenth）值，如果是高光谱还必须有波段宽度（FWHM）,这两个参数都可以通过编辑头文件信息输入（Edit Header）。,数据类型,支持四种数据类型：浮点型（floating）、4-byte signed integers, 2-byte signed integers,以及 2-byte unsigned integers。,数据存储类型： ENVI标准栅格格式文件，且是BIP或者BIL。,波谱范围：flaash能够做的数据光谱范围是0.42500m。,使用,ENVI,大气校正模块,基本参数设置,传感器基本信息设置,使用ENVI大气校正模块大气模型,Model Atmosphere,Water Vapor,(std atm-cm),Water Vapor (g/cm2),Surface Air Temperature,Sub-Arctic Winter (SAW),518,0.42,-16 C or 3 F,Mid-Latitude Winter (MLW),1060,0.85,-1 C or 30 F,U.S. Standard (US),1762,1.42,15 C or 59 ,Sub-Arctic Summer (SAS),2589,2.08,14 C or 57 ,Mid-Latitude Summer (MLS),3636,2.92,21 C or 70 ,Tropical (T),5119,4.11,27 C or 80 ,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,水气去除设置Water Retrieval,水气去除设置，采用两种方式对水气进行去除：,利用水气去除模型恢复影像中每个像元的水气量,使用水气去除模型，数据必须具有,15nm,以上波谱分辨率，且至少覆盖以下波谱范围之一：,1050-1210 nm (,对应,1135 nm),870-1020 nm (,对应,940 nm),770-870 nm (,对应,820 nm),单一的水气因数用于整体影像，默认是,1,，,对于多光谱数据使用水气去除模型，可以在多光谱设置中手动设置水气波段,气溶胶模型（Aerosol Model）,提供四种标准,MODTRAN,气溶胶模型,Rural,（乡村）、,Urban,（城市）、,Maritime,（海洋）、,Tropospheric,（对流层，能见度在,40km,以上）,两种气溶胶去除方法,2-Band,（,K-T,）方法,(,类似模糊减少法,),，如果没有找到适应的黑值（一般是阴影区或者水体），系统将采用能见度值来计算；所以即使选择了该选项也要给能见度。,选择,None,，采用能见度值参与气溶胶去除，能见度值大约参考值参见表,天气条件,能见度,晴朗,40 to 100 km,中等雾、阴霾,20 to 30 km,厚雾、阴霾,15 km,或者更少,高级设置,光谱定义文件：内置,AVIRIS,、,HYMAP,、,HYDICE,、,HYPERION,、,CASI,、,AISA,。,气溶胶厚度系数：用于技术邻域效应范围。一般值为,12km,。,CO2,混合比率：默认为,390ppm,。,使用领域纠正。,使用以前的,MODTRAN,模型计算结果。,设置,MODTRAN,模型的光谱分辨率（推荐值,5 cm-1),。,设置,MODTRAN,多散射模型。,提供三种模型供选择,Isaacs,，,DISORT,和,Scaled DISORT,。默认是,Scaled DISORT,和,streams,为,8,，这种模型对于小于,1000nm,具有较高的精度；,天顶角,方位角（针对非星下点传感器）。,输出反射率缩放系数（,Output Reflectance Scale Factor,）：为了降低结果储存空间，默认反射率乘于,10000,。,输出结果,表面反射率影像,水气含量数据,云图,日志文件,FLAASH大气校正模板文件,3、物质制图与识别、探测,波谱库,波谱分析,端元波谱收集,高光谱制图与识别,波谱库,ENVI,波谱库 （安装目录下,spec_lib,）,Jet Propulsion Laboratory 0.42.5um 160,种纯矿物波谱,美国地质调查局,(USGS) 0.42.5um 500,种质优矿物波谱,和几个植被波谱,Johns hopkins university,0.414um,矿物波谱,IGCP264 (,项目,),26,个质优样品应用波谱仪测量得到,打开波谱库（,spectral/spectral libraries/view),显示波谱曲线,(,点击）,创建波谱库（,spectral/spectral libraries/builder),波谱库的创建与浏览,输入波长范围,输入光谱,从图像中获取,外部文件（二进制）导入,ASD,波谱仪,波谱库交互,波谱库查看、编辑和分析,波谱分割,波谱重采样,波谱分析工具,Spectral Analyst,波谱分析首先需要打开一个波谱库，然后将未知波谱与波谱库中的波谱进行匹配处理，该工具运用,波谱角分类,波谱特征拟和二进制编码技术,对一未知波谱与波谱库中要素的匹配进行排序,输出一个列表,按照波谱匹配的好坏依次排列,并纪录一个总体的得分,.,匹配时需要设置三种方法所占的权重,权重是任意的,最后输出一个总体得分,得分越高,表明匹配效果越好,.,光谱识别流程,MNF变换,重要作用,用于判定图像内在的维数,分离数据中的噪声,减少计算量,弥补了主成分分析在高光谱数据处理中的不足,计算时需要输入的参数,统计信息的图像范围,shift diff subset,噪声统计文件（可以用到另一副图像上做变换）,MNF,统计文件,(,反变换的时候要用,),Mnf,变换输出波段选择,(,根据特征值选择输出波段,),利用,MNF,变换后的图像进行端元提取,经过,MNF,变换后的图像波段之间的相关性有效地降低了，并且出现了多个拐点，这些拐点就是我们要找的端元,Endmembers,。,波谱曲线对比编辑端元,生成端元地物平均波谱曲线,利用波谱分析工具进行端元识别,识别每一类地物后,将端元导出成,ROI,PPI,纯净像元指数,生成的结果是一副灰度的影像,DN,值越大表明像元越纯,.,作用及原理,纯净像元指数法对图像中的像素点进行反复迭代，可以在多光谱或者高光谱影像中寻找最,”,纯,”,的像元,.(,通常基于,MNF,变换结果来进行,),纯净像元指数可以将,N,维散点图迭代映射为一个随机单位向量来计算,每次映射的极值像元被记录下来,并且被标为极值的总数也被记录下来,.,按照多次映射每个像元被记录为极值像元的次数来决定该像元是否为纯净像元,.,计算时需要输入的参数,进行迭代的次数,设置域值的系数,(,极值像元的域值,),数据二次采样,(,减少内存,但不能太小,),处理时会出现纯净像元指数图,X-,像元迭代次数,Y-,极值像元总数,该图显示了,PPI,做为迭代次数的函数,当所有极值像元都已经找到后,它应该接近于一条水平线,处理结果,: DN,值越大表明像元越纯,.,PPI,纯净像元指数提取,ROI,设置波段域值提取,PPI_ROI,通常只设最小值,(,值越大,像元越纯净,),N维可视化器,提取,PPI_ROI,后打开,N,维可视化器,选择,MNF,变换结果,系统默认导入,PPI_ROI,PPI_ROI,端元提取,波谱曲线对比编辑,ROI,生成地物平均波谱,波谱分析,端元识别,端元波谱收集,端元波谱收集器,ASCII,文件,ASD,测量文件,波谱库,ROI/,矢量文件,统计文件,剖面波谱图,N,维可视化分析器,二维散点图,SMACC,工具,所有监督分类与高光谱制图工具,高光谱制图,二进制编码,波谱角分类,线性波段预测（LS-Fit）,线性波谱分离,光谱信息散度,匹配滤波,混合调谐匹配滤波（MTMF）,包络线去除,光谱特征拟合,多范围光谱特征拟合,高光谱制图练习,从标准波谱库中获取端元波谱,打开,jpl_lib/jpl1.sli,波谱库文件，获得以下物质的波谱,ALUNITE SO-4A,BUDDINGTONITE FELDS TS-11A,CALCITE C-3D,CHABAZITE TS-15A,ILLITE PS-11A,KAOLINITE WELL ORDERED PS-1A,选择,Spectral Angle Mapper,制图方法,波谱沙漏工具目标识别,向导式波谱分析工具,影像亮度值定标为反射率,最小噪声分离（,MNF,）,纯净像元指数（,PPI,）,N,维散度分析,选择终端单元,地物制图（地物识别）,Apparent Reflectance,MNF,PPI,n-D,ID,Map Distribution,and Abundance,作用：将地物波谱与端元波谱进行对比来地物制图（识别地物）,SAM与BandMax结合进行目标探测,基于,BandMax,向导的,SAM,目标探测引导你完成高光谱影像的目标探测。向导的,BandMax,部分能帮你找到最佳的波谱子集以区分背景和目标，并节省处理的时间。,向导主要有以下几个步骤：,选择输入,/,输出文件,选择目标,选择目标波谱,选择背景,选择需要抑制的背景信息,利用,BandMax,计算有效波段,识别对,SAM,分析中有效的波段,选择最大角阈值,定义,SAM,最大角,检验制图成果, SAM,分析以及成果检验。,4、植被分析,植被指数计算器,植被分析,植被指数计算器,提供了,6,类,27,种植被指数的计算,绿度,Greenness,光利用率,Light Use Efficiency,氮,干旱或炭衰减,Dry or Senescent Carbon,Stress Pigments,冠层水分含量,Canopy Water Content,提供了生物物理学交叉检验,通过植被指数计算器统一进行计算,提供了每种植被指数的详细资料和计算公式,能够根据影像信息自动显示可计算的植被指数,没有一种其它的商业软件提供如此丰富的植被指数计算,植被分析工具,农作物胁迫,Agricultural,Stress,易燃性分布分析,Fire,Fuels,林木,健康分析,Forest,Health,植被的,抑制,农作物胁迫Agricultural Stress,创建农作物胁迫的空间分布图,判断出适合农作物生长的区域，可用于精确农业分析；,绿度：来标识闲置农田、虚弱的植被、健康的作物；,光利用效率：标识植被生长率；,冠层氮含量：用于估计相关的氮等级；,叶绿素：突出植被胁迫；,冠层水分含量：标识水胁迫的层次。,易燃性分布分析,创建某一区域植被易燃性的空间分布图,用于森林规划，也可用它来分析城郊混合区的火灾风险,绿度：显示大多数的绿色指数；,冠层水分含量：标识植被的冠层水分含量；,干旱或炭衰竭：突出非生物指数；,林木健康分析工具,创建整个森林区域健康程度的空间分布图,用于检测病虫害以及枯萎病的发生情况，也可以用于评估某地区的木材收获量,绿度：表面绿色植被的分布；,叶绿素：标识类胡萝卜素以及花青素的含量；,冠层水分含量：标识水分含量；,光使用效率：标识森林生长率；,植被抑制,从多光谱和高光谱减少或移除植被光谱信息,利用红外与近红外波段,方便地质或城市解译,