完整版遥感地质蚀变异常信息提取实验.docx

1、完整版遥感地质蚀变异常信息提取实验“基础地质学”创新性实验遥感地质蚀变异常信息提取实验“基础地质学”国家级教学实验示范中心二二年二月二十二日1 实验目的通过学习基于地学软件的矿化蚀变异常信息提取的流程，了解利用主成分分析方法提取矿化蚀变异常信息提取的原理，掌握ENVI4.6以及Surfer9软件的基本使用方法。2 实验内容2.1 熟悉遥感影像的辐射定标的方法与流程通过学习遥感影像的辐射定标的方法与流程，掌握用ENVI软件进行辐射定标的具体操作，加深对遥感图像处理知识的理解。2.2 掌握遥感影像的波段合成、投影转换、影像裁剪的方法通过学习遥感影像的波段合成、投影转换、影像裁剪的方法与流程，掌握用

2、ENVI软件进行这些操作的具体步骤，加深对遥感图像处理知识的理解。2.3 掌握ETM+遥感影像的Flaash大气校正、掩膜的应用方法通过学习ETM+遥感影像的Flaash大气校正，掩膜的建立与应用方法，熟练掌握ENVI软件进行这些操作的具体步骤，加深对大气校正知识的理解。2.4 掌握ETM+遥感影像羟基和铁染异常信息提取的方法与流程通过学习利用主成分分析方法对ETM+遥感影像进行羟基和铁染异常信息提取的方法与流程，熟练掌握ENVI软件进行这些操作的具体步骤，加深对矿化蚀变信息提取知识的理解。2.5 掌握ENVI与Surfer软件协同制图的方法。通过学习使用ENVI软件提取ETM+遥感影像矿化蚀

3、变信息与Surfer软件绘制矿化蚀变异常空间分布图，综合掌握ENVI与Surfer软件协同制图的方法，加深对所学知识的理解。3 实验要求上课时，需认真听取实验老师的细致讲解；练习时，按照实验指导书的步骤，自己亲自动手使用ENVI4.6软件提取矿化蚀变异常信息；以及在Surfer软件中制作矿化蚀变异常信息空间分布图。4 实验条件4.1 软件平台：ENVI4.6、Surfer9ENVI4.6软件为遥感图像处理软件，用于提取矿化蚀变异常信息；Surfer 9软件为一款以画三维图（等高线,image map, 3d surface）的软件，具有的强大插值功能和绘制图件能力，使它成为用来处理XYZ数据的

4、首选软件，是地质工作者必备的专业成图软件。4.2 遥感数据源：金川地区Landsat7 ETM+遥感影像Landsat7于1999年4月15日发射升空后，由于其优越的数据质量，以及与以前的Landsat系列卫星保持了在数据上的延续性，现在已成为中国遥感卫星地面站的主要产品之一，2003.5月出现故障退役。5 实验原理5.1 蚀变异常提取的地质依据目前认为中等强度以上的蚀变带对于TM（ETM+）蚀变信息提取是十分有利的。有时尽管有蚀变岩存在，但不一定有矿，然而围岩的剧烈而较大范围的蚀变常常与大矿及富矿石的生成互为隶属，大型特大型内生热液矿床一般均有强烈且较大范围的围岩蚀变，并且具有分带现象（如斑

5、岩铜矿）。这便是以找矿（首先是大矿、富矿）为最终目的蚀变遥感异常提取的地质依据。5.2 蚀变异常提取的物理依据从Hunt G.RC（1978）和他领导的实验室的研究成果以及阎积惠等（1995）依据矿物反射波谱特征吸收谱带特点的定性分类研究中可以知道：主要造岩矿物在可见光近红外光谱（0.35-0.25um）并不产生具有鉴定意义的反射谱带，其光谱特征主要由岩石中为数不多的次要矿物决定：一是含铁（Fe2+， Fe3+）基团产生，含铁矿物主要有角闪石、赤铁矿、褐铁矿、针铁矿、磁铁矿、黄钾铁矾等，他们在TMl、TM4波段有强的吸收带，若岩石中含多量的Fe3+，而含Fe2+较少，这类岩石的主要吸收谱带位于

6、TM4和TM1波段，反射波长相当于TM3波段的电磁波。若含大量的Fe2+、含Fe3+较少，则主要吸收谱带位于1波段，对于波长相当于TM2波段的电磁波有某种程度的反射。二是含羟基（OH-）、水（H2O）或碳酸根（CO32-）基团产生，羟基的吸收谱带主要有二处：2.2um，2.3um，由于OH-在2. 22.3um附近存在强吸收谷（称为羟基谱带），使得TM7产生低值，TM5产生高值，含羟基矿物大多为次级蚀变矿物，如高岭土、叶腊石、云母类矿物、绿泥石、绿帘石等，水在1.4um和1.9um处有特征吸收带。含碳酸根矿物主要有五个特征吸收谱带（1.9-2.55um），较强的两个在2.35um和2.55um

7、波长处（称为碳酸根谱带），相对较弱的在1.9um，2.0um，2.16um三处。5.2.1 利用主成分分析方法提取矿化蚀变信息遥感图像各波段间存有一定的相关性，为了减少相关性对分类的影响，常使用主分量分析法去相关。主分量分析基于变量之间的相互关系，在信息总量守恒的前提下，利用线性变换的方法来实现去相关性。由于所获各主分量之间不相关，故各主分量之间信息没有重复或冗余。故主分量分析这一基本性质在蚀变异常信息提取中被充分利用并很有成效。5.2.2 铁染蚀变异常分析针对各个波段波谱所反映的特征，选取TM1、TM3、TM4、TM5这4个波段作为组合波段做主成分分析。这是因为铁氧化物的特征光谱信息集中在了

8、TM14波段，在TM4和TM1波段有吸收峰，在TM3波段无特征吸收而呈高反射。同时为了避免含羟基和碳酸根矿物的干扰，在选取波段组合时舍弃了TM7波段。对TM1、TM3、TM4、TM5应用掩膜做主成分分析，统计分析如下：特征向量 （Eigenvectors）主成分Band 1Band 3Band 4Band 5PC10.4153040.5278860.4966230.549750PC2-0.204294-0.5514550.807515-0.045623PC3-0.636068-0.097958-0.1858670.742479PC4-0.6174170.6384660.258339-0.380

9、024从上表可以看出，PC1主要反映了TM3和TM5波段的信息；PC2主要反映了TM4波段的信息；PC3反映TM5波段的加信息和TM1波段的减信息；PC4反映了TM1波段的减信息和TM3波段的加信息。根据铁染类蚀变矿物的波谱特征，包含这类蚀变信息的图像应该具有TM3与TM1或TM3与TM4具有相反的贡献值且绝对值较大，从而选择PC2分量，故PC4分量中暗色调部分表征了铁染信息。5.2.3 含羟基类矿物和含CO32-矿物蚀变异常分析选取TM1、TM4、TM5、TM7这4个波段作为组合波段做主成分分析【4】。这是因为粘土类矿物（含羟基矿物）和含CO32-矿物的特征光谱信息集中在TM5和TM7波段，

10、在TM7波段为特征吸收带，在TM5相对高反射。主成分变换的波段组合的选择也是以此为依据的，由于可见光波段对铁氧化物敏感，为了避免铁氧化物信息的干扰，故只选择了一个可见光波段参与运算。对TM1、TM4、TM5、TM7应用掩膜做主成分分析，统计分析如下：特征向量 （Eigenvectors）主成分Band 1Band 4Band 5Band 7PC10.4250780.5101330.5650800.489650PC2-0.1433960.809782-0.146893-0.549649PC3-0.8701090.0487390.4575770.176519PC40.204102-0.285731

11、0.670619-0.653429从上表可以看出，PC1主要反映的是TM4和TM5波段的信息；PC2主要反映的是TM4波段的信息；PC3反映的是TM1波段的信息；PC4反映的是TM5和TM7波段的信息，且符号相反。根据含羟基类矿物和含CO32-矿物的波谱特征，PC4图像中的高亮度信息就是表征含羟基类矿物和含CO32-矿物的蚀变异常信息。6 实验步骤6.1 金川地区ETM+遥感影像辐射定标在ENVI 4.6软件中，单击菜单栏中的fileopen image file工具，打开金川地区Landsat7 ETM+遥感影像中band1影像文件。然后，单击菜单栏中的Basic Toolspreproce

12、ssingCalibration UtilitiesLandsat TM工具，为金川地区ETM+遥感影像中band1影像文件进行辐射定标，如图 61所示。图 61打开辐射定标工具对话框在TM Calibration Input File弹窗下的Select Input File栏中选择Landsat 7的band1文件，然后单击OK，如图 62所示。图 62 TM Calibration Input File对话框接来下，在TM Calibration Parameters弹窗中，单击Get Calibration Parameters from Web选项，如图 63所示。图 63 TM C

13、alibration Parameters对话框在TM Web Calibration Parameters弹窗中，选择Date/Path/Row/Band选项，根据金川地区Landsat7 ETM+遥感影像的头文件信息，依次填写下面红色标注处的信息，如图 64所示。注意：每一景影像的信息都不一样，一定要根据头文件信息填写。图 64 TM Web Calibration Parameters对话框填写完成后，在TM Calibration Parameters中依次填写下图红色标注部分的信息，如图 65所示。图 65 TM Calibration Parameters对话框单击OK，即完成了金

14、川地区Landsat7 ETM+遥感影像中band1影像文件的辐射定标操作；接下来，依次对band2、band3、band4、band5、band7波段进行辐射定标操作。6.2 金川地区ETM+遥感影像不同波段的合成首先，在ENVI 4.6软件中，单击菜单栏中的fileopen image file工具，打开经过辐射定标的金川地区ETM+遥感影像的band1、band2、band3、band4、band5、band7数据文件，如图 66所示。图 66 Available Band List 对话框然后，单击菜单栏中的Basic ToolsLayer Stacking工具；如图 67所示。图 6

15、7 打开波段合成工具在Layer Stacking Parameters弹窗中，单击Import File按钮，将经过辐射定标的金川地区ETM+遥感影像的band1、band2、band3、band4、band5、band7数据文件全部添加进来。之后，单击Reorder Files按钮，在Reorder Files弹窗中将顺序从上到下调整为：band1、band2、band3、band4、band5、band7。选择好保存的路径，单击OK即完成了金川地区ETM+遥感影像不同波段的波段合成操作。图 68 Layer Stacking Parameters对话框6.3 定义金川地区ETM+遥感影像

16、的地理坐标首先，在ENVI 4.6软件中打开金川地区ETM+影像（波段合成后），并加载彩色合成影像；在主菜单中单击Map中选择Customize Map Projections选项。 图 69 打开定义地理坐标工具然后，在新弹窗的Projection Name中输入投影的名字,如jinchuantouying；在False easting中输入500000.00；False northing中输入0.00；在Projection Type中选择Transverse Mercator；在Projection Datum中选择xian-80；在Latitude中分别输入0 ，0，0.00；在Lon

17、gitude中输入102，0，0.00。Scale factor中输入1。图 610 定义地理坐标对话框接来下，输入完以上数据后，在菜单栏中点击Projection，选择Add New Projection，即加载这个投影信息。图 611定义地理坐标对话框最后，点击File中Save Projections，保存这个投影信息。图 612定义地理坐标对话框6.4 对金川地区ETM+遥感影像的地理坐标进行投影转换首先，在主菜单栏中选择Convert Map Projection，进行投影转换。图 613打开投影转换工具其次，在新弹窗的Select Input File里选择要进行转换的影像，如甘肃

18、金川jubu.img；然后，单击OK。图 614 投影转换输出影像窗口接着，在新弹窗中Resampling里选择Bilinear方法；同时，选择好要保存的路径。然后，单击Change Proj。图 615 投影转换参数设置对话框再在新弹窗的Select New Projection中选择前面保存好的投影jinchuantouying；然后，单击OK。图 616 投影选择对话框下图为投影后的彩色影像。图 617 投影转换后的影像6.5 对金川地区ETM+遥感影像进行裁剪首先，在ENVI 4.6软件中打开波段合成后的金川地区ETM+遥感影像；在菜单栏中单击Basic ToolsResize Dat

19、a（Spatial/Spectral）工具。图 618 打开裁剪工具在Resize Data Input File弹窗中选择要进行裁剪的影像，单击Spatial Subset按钮；图 619输入裁剪数据对话框在Select Spatial Subset弹窗中，可以通过Image、Map、File、ROI/EVF四种方法进行裁剪，根据需要裁剪出相应的研究区域。最后，单击OK，即完成裁剪操作。图 620 选择空间裁剪区域对话框6.6 对金川地区ETM+遥感影像进行FLAASH大气校正6.6.1 数据转换首先，在ENVI软件中打开裁剪后的金川地区ETM+遥感影像；再单击Basic ToolsConv

20、ert Data（BSQ，BIL，BIP）工具，对需要进行FLAASH大气校正的数据进行数据转换操作；图 621 打开数据转换工具在Convert File Input File弹窗中，选择要进行数据转换的文件，再单击OK；图 622输入数据转换文件窗口再在Convert File Parameters弹窗中选择BIL格式，选择好存储路径后，单击OK，即完成了数据格式的转换。图 623 转换参数设置窗口6.6.2 编辑头文件信息然后，在Available Bands List窗口里选择刚才转换格式的文件，右击该文件，点击Edit Header选项。图 624 编辑头文件信息在弹出窗口中，单击E

21、dit Attributes，在下拉菜单中选择Wavelengths进行下一步编辑。图 625 编辑头文件中的中心波长信息在弹出的对话框中，依次将band1、band2、band3、band4、band5、band7波段的中心波长设置为: 0.4825、0.565、0.66、0.825、1.65、2.22；之后，单击OK。图 626 中心波长编辑对话框6.6.3 进行FLAASH大气校正最后，再单击菜单栏中Basic ToolsPreprocessingCalibration UtilitiesFLAASH工具，进行FLAASH大气校正；图 627 打开FLAASH大气校正工具在弹出的窗口中，

22、单击Input Radiance Image按钮，选择前步得到的影像数据后，单击OK；再在Radiance Scale Factors弹窗中，按下图参数进行设置，单击OK即可。图 628 FLAASH大气校正数据输入、影像选择及参数设置图之后，在Sensor Type中选择JNKNOWN-MSMultispectralLandsat TM7工具；图 629 Sensor Type参数设置图接着,根据金川地区ETM+遥感影像文件的头文件信息依次将下图红色标记部分的信息填完，并设置好保存的路径；然后，单击Multispectral Settings按钮，进行参数设置操作；之后，在单击Advance

23、d Settings按钮，进行参数设置操作。图 630 FLAASH大气校正工具对话框A.在单击Multispectral Settings按钮，进行参数设置操作时；在新的弹出窗口中，安全按照下面红色标注部分的参数进行设置；图 631 Multispectral Settings对话框按照红色标注信息继续设置其他的选项；之后，单击OK即可。图 632 Multispectral Settings对话框B.在单击Advanced Settings按钮，进行参数设置操作时，完全按照下面红色标注的信息进行填写。填写完后，单击OK。在图？中单击应用，即完成了FLAASH大气校正操作。图 633 Adv

24、anced Settings对话框6.7 简易去除ETM+遥感影像的干扰信息6.7.1 建立ROI（感兴趣区）首先，在IMAGE窗口中点击Overlay，选择Region of Interest选项。图 634 打开ROI工具然后，在弹出的#1 ROI Tool窗口中的Window里选择在Scroll中建立ROI。左图红色部分为建立的ROI。图 635 建立ROI区域图最后，将建立的ROI保存。图 636 保存ROI区域6.7.2 建立掩膜首先，在主菜单栏中单击Basic Tools，然后单击Masking，选择Build Mask，建立掩膜。图 637 打开掩膜工具其次，在新弹窗的菜单栏里点

25、击Options，选择Import ROIs。图 638 定义掩膜对话框最后，在弹出的窗口中，选择前面保存好的ROI；然后，单击OK建立掩膜。图 639Mask Definition Input ROIs对话框6.7.3 应用掩膜首先，在主菜单栏里单击Basic Tools，在下拉菜单里单击Masking，再选择Apply Mask。图 640 打开掩膜应用工具其次，在弹出的窗口中的Select Input File中选择投影后的影像，如本实验中的touying；然后再在Select Mask Band中选择mask-def文件；单击OK。图 641 Apply Mask Input File

26、对话框然后，在新弹窗中选择保存路径即可；单击OK。图 642 Apply Mask Parameters对话框6.7.4 掩膜的反选首先，在主菜单栏的Basic Tools中选择Band Math工具。图 643 打开波段运算工具其次，在Enter an expression中输入1-B1；再点击Add to List；单击OK即可。图 644 波段运算对话框最后，在Available Bands List中选择Mask Band，即定义了B1；然后选择好保存路径，再点击Ok。图 645 Variables to Bands Pairings对话框6.7.5 掩膜反选后的应用再在主菜单栏中Ba

27、sic Tools单击Masking，选择Apply Mask。图 646 打开应用掩膜工具其次，在弹出的窗口中的Select Input File中选择投影后的影像，如本实验中的touying；然后再在Select Mask Band中选择mask-fan文件；单击OK。图 647 Apply Mask Input File对话框最后，再在弹出的窗口中选择好保存路径。至此，应用掩膜的工作完成。6.8 主成分分析首先，打开经过掩膜之后的影像。然后，在主菜单栏中单击Transform，然后在下拉菜单Principal Components的Forward PC Rotation中选择Comput

28、e New Statistics and Rotate进行主成分分析工作。图 648打开主成分分析工具接着，在弹出窗口的Select Input File中选择要进行主成分分析的影像，如本实验中的Mask-apply2文件；再在Spectral Subset中选择1345这四个波段进行铁染蚀变异常的主成分分析。图 649Principal Component Input File对话框然后，在弹窗中选择好保存路径并命名，然后点击OK。图 650 Forward PC Parameters对话框紧接着，在主菜单栏中点击Basic Tools，在下拉菜单中单击Statistics，再选择View

29、Statistics File。图 651打开查看统计结果工具然后，单击打开前面刚存好的铁染信息主成分分析后的文件1345.sta。图 652 输入统计分析结果文件最后，对Eigenvector进行分析，确定铁染信息主要对应在第几主成分。图 653 统计分析结果同理，可以用主成分分析的方法分析羟基蚀变异常信息；其中，对羟基蚀变异常进行主成分分析所用到的波段为1457波段。注意：分析过程见：5.3利用主成分分析方法提取矿化蚀变信息6.9 提取蚀变异常信息通过分析特征向量后，得出铁染蚀变异常信息提取中1345文件的第2主成分主要包含铁染信息。首先，加载第2主成分的灰度影像。图 654 第2主成分的

30、灰度影像然后，在Image窗口里的菜单栏中单击Enhance，选择Interactive Stretching。图 655 打开交互式拉伸工具其次，通过调整弹出窗口中两根线条的位置进行交互式拉伸，突出暗色调中铁染蚀变异常信息。图 656 交互式拉伸对话框最后，再将交互式拉伸好的影像保存。图 657 保存拉伸后的结果6.10 用Surfer软件修饰铁染蚀变异常信息 首先，加载前一步刚保存好的影像；再在主菜单栏中单击File，在下拉菜单中单击Save File As，再选择ASC。图 658 将拉伸的结果保存为ASCII数据其次，在Output Style中选择X Y Z Format；然后，选择好路径；再单击OK。这样，就可以将铁染蚀变异常的信息导入到Surfer 软件中进行处理。图 659 设置ASCII参数接下来，在surfer软件中主菜单栏中单击Grid，选择Data，将上一步存储的含有铁染异常信息的1345.txt文件转换为1345.grd文件。图 660 在surfer中打开ASCII数据文件之后，点击New Contour Map工具，将上一步保存的1345.grd数据导入。图 661在surfer软件中打开下图为导入数据后的效果图。图 662 导入1345.grd数据的结果图然后，双击该图。单击Line，对线的颜色进行修改，同时可以单击其他工具