遥感地质蚀变异常信息提取实验

遥感地质蚀变异常信息提取实验年二月二十日 .....................................................3 ............4.......4 5.2.1利用主成分分析方法提取矿化蚀变信息...........................................................5 ETM9 ....................................................................................9 ........................................................................................10 6.7.1建立ROI(感兴趣区)......................................................................................11 ......................................................................................12 ................................................13 ......151实验目的通过学习基于地学软件的矿化蚀变异常信息提取的流程，了解利用主成分分析方法提取矿r2实验内容2.1熟悉遥感影像的辐射定标的方法与流程通过学习遥感影像的辐射定标的方法与流程，掌握用ENVI软件进行辐射定标的具体操作，2.2掌握遥感影像的波段合成、投影转换、影像裁剪的方法通过学习遥感影像的波段合成、投影转换、影像裁剪的方法与流程，掌握用ENVI软件进行这些操作的具体步骤，加深对遥感图像处理知识的理解。2.3掌握ETM+遥感影像的Flaash大气校正、掩膜的应用方法2.4掌握ETM+遥感影像羟基和铁染异常信息提取的方法与流程通过学习利用主成分分析方法对ETM+遥感影像进行羟基和铁染异常信息提取的方法与流3实验要求上课时，需认真听取实验老师的细致讲解；练习时，按照实验指导书的步骤，自己亲自动4实验条件三维图(等高线,imagemap,3dsurface)的软件，具有的强大插值功能和绘制图件能力，使它成为用来处理XYZ数据的首选软件，是地质工作者必备的专业成图软件。4.2遥感数据源：金川地区Landsat7ETM+遥感影像t5实验原理5.1蚀变异常提取的地质依据目前认为中等强度以上的蚀变带对于TM(ETM+)蚀变信息提取是十分有利的。有时尽象(如斑岩铜矿)。这便是以找矿(首先是大矿、富矿)为最终目的蚀变遥感异常提取的地质5.2蚀变异常提取的物理依据从HuntG.RC(1978)和他领导的实验室的研究成果以及阎积惠等(1995)依据矿物反射究中可以知道：主要造岩矿物在可见光——近红外光谱 (0.35-0.25um)并不产生具有鉴定意义的反射谱带，其光谱特征主要由岩石中为数不多的次要一是含铁(一是含铁(Fe2+，Fe3+)基团产生，含铁矿物主要有角闪石、赤铁矿、褐铁矿、针铁矿、磁铁矿、黄钾铁矾等，他们在TMl、TM4波段有强的吸收带，若岩石中含多量的Fe3+，而含Fe2+较少，这类岩石的主要吸收谱带位于TM4和TM1波段，反射波长相当于TM3波段的电磁波。若含大量二是含羟基(OH-)、水(H2O)或碳酸根(CO32-)基团产生，羟基的吸收谱带主要有二吸收谱带(1.9-2.55um)，较强的两个在2.35um和2.55um波长处(称为碳酸根谱带)，相对较利用主成分分析方法提取矿化蚀变信息变换的效。2.2铁染蚀变异常分析氧化物的特征光谱信息集中在了TM1—4波段，在TM4和TM1波段TM吸收而呈高反射。同时为了避免含羟基和碳酸根矿物的干扰，在特征向量(Eigenvectors)主主成分Band1Band3Band50.4153040.5278860.4966230.549750-0.204294-0.5514550.807515-0.045623-0.636068-0.097958-0.1858670.742479-0.6174170.6384660.258339-0.380024和TM3波段的加信息。根据铁染类蚀变矿物的波谱特征，包含这类蚀变信息的图像应该具有35.2.3含羟基类矿物和含CO2-矿物蚀变异常分析3类矿物(含羟基矿物)和含CO32-矿物的特征光谱信息集中在TM5和TM7波段，在TM7波段为特征吸收带，在TM5相对高反射。主成分变换的波段组合的选择也是以此为依据的，由于可见光波段对铁氧化物敏感，为了避免铁氧化物信息的干扰，故只选择了一个可见光波段参与特征向量(Eigenvectors)主主成分Band1Band5Band70.4250780.5101330.5650800.489650-0.1433960.809782-0.146893-0.549649-0.8701090.0487390.4575770.1765190.204102-0.2857310.670619-0.653429TMTM段的信息，且符6实验步骤6.1金川地区ETM+遥感影像辐射定标ETM+遥感影像中band1影像文件。然后，单击菜单栏中的BasicTools——preprocessing——OK所示。TMCalibrationInput对i框leb项，如图-所示。Landsat7ETM+遥感影像的头文件信息，依次填写下面红色标注处的信息，如图-所示。填写。d6.2金川地区ETM+遥感影像不同波段的合成ile图-6所示。图-打开波段合成工具6.3定义金川地区ETM+遥感影像的地理坐标ENVIETM并加载彩色合成影像；s图-打开定义地理坐标工具6.4对金川地区ETM+遥感影像的地理坐标进行投影转换g6.5对金川地区ETM+遥感影像进行裁剪首先，在ENVI4.6软件中打开波段合成后的金川地区ETM+遥感影像；在菜单栏中单击BasicTools——ResizeData(Spatial/Spectral)工具。图-0选择空间裁剪区域对话框6.6对金川地区ETM+遥感影像进行FLAASH大气校正6.6.1数据转换首先，在ENVI软件中打开裁剪后的金川地区ETM+遥感影像；再单击BasicTools——图-0FLAASH大气校正工具对话框图-1打开数据转换工具图-2输入数据转换文件窗口图-3转换参数设置窗口6.6.2编辑头文件信息然后，在AvailableBandsList窗口里选择刚才转换格式的文件，右击该文件，点击Edit图-4编辑头文件信息图-5编辑头文件中的中心波长信息图-6中心波长编辑对话框reprocessingCalibrationUtilitiesFLAASH图-7打开FLAASH大气校正工具K图-8FLAASH大气校正数据输入、影像选择及参数设置图图-9SensorTye数设置图并设置好保存的路径；然后，单击MultispectralSettings按钮，进行参数设置操作；之后，在s图-0打开掩膜应用工具A.在单击MultispectralSettings按钮，进行参数设置操作时；在新的弹出窗口中，安全按K6.7简易去除ETM+遥感影像的干扰信息6.7.1建立ROI(感兴趣区)图-4打开ROI工具图-5建立ROI区域图图-6保存ROI区域6.7.2建立掩膜图-7打开掩膜工具图-8定义掩膜对话框6.7.3应用掩膜plyMaskInputile6.7.4掩膜的反选图-3打开波段运算工具图-4波段运算对话框6.7.5掩膜反选后的应用图-6打开应用掩膜工具plyMaskInputile6.8主成分分析首先，打开经过掩膜之后的影像。然后，在主菜单栏中单击Transform，然后在下拉菜单图-8打开主成分分析工具接着，在弹出窗口的SelectInputFile中选择要进行主成分分析的影像，如本实验中的紧接着，在主菜单栏中点击BasicTools，在下拉菜单中单击Statistics，再选择View图-1打开查看统计结果工具图-2输入统计分析结果文件r图-3统计分析结果同理，可以用主成分分析的方法分析羟基蚀变异常信息；其中，对羟基蚀变异常进行主成6.9提取蚀变异常信息图-4第2主成分的灰度影像图-5打开交互式拉伸工具图-6交互式拉伸对话框图-7保存拉伸后的结果6.10用Surfer软件修饰铁染蚀变异常信息首先，加载前一步刚保存好的影像；再在主菜单栏中单击File，在下拉菜单中单击Save图-8将拉伸的结果保存为ASCII数据图-9设置ASCII参数图-61在surfer软件中打开图-62导入1345.grd数据的结果图图-63参数编辑对话框图-64修饰后的结果图-65地图参数设置对话框图-66打开属性工具图-67加载遥感地图的过程图-68加载遥感地图后的效果再在objectman