整理的高光谱数据的矿物填图技术ppt课件

基于高光谱数据的矿物填图 航天高光谱遥感器系统 LEWIS 1997 美国NASAOrbview 4 2001 发射失败 美国MODIS是EOS AM1卫星 1999年12月发射 和EOS PM1 2002年5月发射 上的主要探测仪器 中分辨率成像光谱仪 也是EOSTerra平台上唯一进行直接广播的对地观测仪器 通过MODIS可以获取0 4 14 m的36个波段的高光谱数据美国2000年7月发射的MightSat II卫星上搭载的傅立叶变换高光谱成像仪 FourierTransformHyperspectralImager FTHSI 是干涉成像光谱仪的成功典范 由于MightySatII出色的性能和成功运行 其研究组获得了美国空军研究实验室 AFRL 的 司令杯团体奖 和 空间运载工具董事会年度团体奖 2001年10月 欧空局 EuropeanSpaceAgency 成功发展了基于PROBA小卫星的紧凑型高分辨率成像光谱仪 CHRIS 并发射成功 CHRIS的地面采样间隔为18m 为14 14Km 在400nm 1050nm之间光谱分辨率为11nm 波段数19 62 本次数据简介 Hyperion传感器搭载于EO 1卫星平台 EO 1 EarthObserving 1 是美国NASA面向21世纪为接替LandSat 7而研制的新型地球观测卫星 于2000年11月发射升空 其卫星轨道参数与LandSat 7卫星的轨道参数接近 之所以设计相同轨道 目的是为了使EO 1和LandSat 7两颗星的图像每天至少有1 4景重叠 以便进行比对 Hyperion是第一个星载民用成像光谱仪 采用了谱像合一技术为使用者提供了242个波段 光谱范围覆盖355 2577nm Hyperion数据的三位光谱剖面图见图4 1 传感器空间分辨率30米 Hyperion产品命名规则 HyperionLevel1R数据文件命名具有一定的规则 EO1SPPPRRRYYYYDDDXXXML GGG VVEO1表示卫星的名字 S表示传感器 其H Hyperion A ALI PPP表示成像目标物所处的WRS系统轨道 RRR表示成像目标物所处的WRS系统的行 YYYY表示成像时的年份 DDD表示成像是处于该年中的第几天 第1个X表示Hyperion传感器的状态 第2个X表示ALI传感器的状态 第3个传感器表示LAC传感器的状态 其中0表示关 1表示开 M表示点模式 P pointedwithincurrentpath row K pointedoutsidecurrentpath row L表示一景的长度 F fullscene P partialscene GGG表示地面接收站 VV表示解译图号 EO1H1270322003301110PZ LGS 01 矿物的光谱特性 在0 4 1 3 m的光谱特性主要取决于矿物晶格结构中存在的铁 铜 镍 锰等过渡性金属元素的电子跃迁 1 3 2 5 m的光谱特性是由矿物组成中的碳酸根 CO32 羟基 OH 及可能存在的水分子 H2O 决定的 3 5 m的光谱特性是由Si O Al O等分子键的振动模式决定的 常见高光谱矿物填图流程 最大噪声分离 MNF 变换 成像光谱的光谱分辨率很高 波段多 数据海量 且波段之间的相关性很强 数据冗余度高 需要在端元选取和矿物识别之前对数据作减维 并弱化噪声 最常用的方法是 最大噪声组分变换 MNF相当于所有波段噪声方差都相等时的主成分分析 可分为两步实现 第一步 先将图像变换到一个新的坐标系统 使变换后图像噪声的协方差矩阵为单位阵 第二步 对变换后的图像进行主成分变换 像元纯度指数 PPI 分析 基于上述MNF变换 排序低的MNF波段被暂时忽略 仅选择高序次波段进一步处理 PPI设计指定光谱极值像元 对应为混合光谱端元 通过反复投影n维散点图到随机单元矢量来计算PPI 记录每次投影的极值像元 注记每个像元被标定为极值的总次数 PPI图像产生 其中 每个像元的DN值对应像元被记录为极值的次数 这些图像的直方图显示被PPI 击中 hit 的分布 从直方图中选择阈值 用于选择最纯的像元以保证被分析的像元数最小 这些像元被输入到分离特定光谱端元的交互式可视化算法中 N维可视化 n DimensionalVisualization 由于PPI算法本身确定的并不是最终的端元 而是从图像众多像元中选出包含有端元像元的较小像元子集 所以 一般要将PPI的处理结果输入到N维可视化工具中 选择出最终的端元像元 在N维可视化中 光谱可视为n维散点图中的一个点 n是波段数 对给定的像元 n维空间中组成n值的点的坐标是每个波段的光谱反射率 这些点在n维空间的分布可用于估计光谱端元数和它们的纯光谱特征 这种混合光谱的凸面几何特征 可用于确定端元光谱数并估算它们的光谱特征 经过PPI选择出的潜在端元光谱输入n维散点图中进行反复旋转以识别出纯端元 较好的端元通常会出现在n维散点图的顶点和拐角处 当一系列的端元点被确定后 就可以将其输入到图像中的感兴趣区 ROI 从图像中提取每个感兴趣区平均反射率光谱曲线作为成像光谱矿物填图的候选端元 矿物识别 高光谱矿物识别的基本原理是高光谱数据的重建光谱与矿物标准光谱或实测光谱的定量对比分析 由己知地物类型的反射波谱 通过波形或特征匹配比较来达到识别地物类型的目的 波谱匹配技术是矿物识别的关键技术 所谓波谱谱匹配 就是通过研究两个波谱曲线的相似度来判断地物的归属类别 常用的波谱匹配的算法有二值编码匹配 BE 波谱角度匹配 SAM 波谱特征拟合 SFF 选择合适的方法进行填图 常用的制图方法主要有 SAMMFMTMF光谱特征拟合线性光谱分解 SAM SAM SpectralAngleMapping 即光谱角制图技术 光谱角技术通过计算一个测试光谱 像元光谱 与参考光谱 实验室光谱等 之间的角度来估算两者之间的相似度 假定图像数据已转换为暗辐射或路径辐射偏差消除后的视反射率 光谱维数与波段数相等 MF MF Matchedfiltering 即匹配滤波 匹配滤波技术最大化突出目标 压缩背景 该技术减少了线性光谱分解中的端元光谱选择的复杂性问题 MTMF MTMF mixturetunedmatchedfiltering 即混合协调匹配滤波 MTMF模型融合线性光谱分解和匹配滤波的技术优势 避免单个方法的缺陷 匹配滤波技术具有突出目标 压制背景的作用 线性光谱分解具有单元相加和正值的特征 该技术方法能探测其它技术不能探测到的微细矿物成分 光谱特征拟合 光谱特征拟合又称之为波段填图算法 它是选择包含目标矿物特定吸收谱带的光谱区间 利用最小二乘拟合方法 比较像元光谱与目标光谱吸收特征的整体形态和吸收深度 线性光谱分解 如果每个光子与像元中的单个端元 endmember 相互作用 则混合模型是线性的 否则为非线性的 非线性混合模型最终要转化为线性模式计算 线性混合模型算法如下 DNb是指定像元在波段b的强度 Fi是端元i的丰度分量 DNib是端元i在波段b的强度 N是端元数 Eb是波段拟合误差 白云母族矿物 白云母族矿物的光谱诊断特征是在波长2200nm附近出现有Al OH的吸收带 随着白云母中Al含量的减少 Al OH的谱带位置逐渐向长波方向移动 从2190mn逐渐变化到2225mn 钠云母 富Al 含Na Al OH谱带位于2190mn 多硅白云母 贫Al 富Mg或Fe 位于2225nm 根据这一特征 不仅可以识别不同种类的白云母 还可以估计白云母中Al的金属含量 白云母的标准波谱曲线如图 绿帘石矿物 绿