教程-ENVI遥感地质蚀变异常信息提取过程(全)

“基础地质学”创新性实验遥感地质蚀变异常信息提取实验

“基础地质学”国家级教学实验示范中心

二〇二二年四月二十八日

目录

1实验目的 (3)

2实验内容 (3)

2.1熟悉遥感影像的辐射定标的方法与流程 (3)

2.2掌握遥感影像的波段合成、投影转换、影像裁剪的方法 (3)

2.3掌握ETM+遥感影像的Flaash大气校正、掩膜的应用方法 (3)

2.4掌握ETM+遥感影像羟基和铁染异常信息提取的方法与流程 (3)

2.5掌握ENVI与Surfer软件协同制图的方法。 (3)

3实验要求 (3)

4实验条件 (4)

4.1软件平台：ENVI4.6、Surfer9 (4)

4.2遥感数据源：金川地区Landsat7 ETM+遥感影像 (4)

5实验原理 (4)

5.1蚀变异常提取的地质依据 (4)

5.2蚀变异常提取的物理依据 (4)

5.2.1利用主成分分析方法提取矿化蚀变信息 (5)

5.2.2铁染蚀变异常分析 (5)

5.2.3含羟基类矿物和含CO32-矿物蚀变异常分析 (6)

6实验步骤 (7)

6.1金川地区ETM+遥感影像辐射定标 (7)

6.2金川地区ETM+遥感影像不同波段的合成 (9)

6.3定义金川地区ETM+遥感影像的地理坐标 (10)

6.4对金川地区ETM+遥感影像的地理坐标进行投影转换 (12)

6.5对金川地区ETM+遥感影像进行裁剪 (15)

6.6对金川地区ETM+遥感影像进行FLAASH大气校正 (17)

6.6.1数据转换 (17)

6.6.2编辑头文件信息 (18)

6.6.3进行FLAASH大气校正 (19)

6.7简易去除ETM+遥感影像的干扰信息 (23)

6.7.1建立ROI（感兴趣区） (23)

6.7.2建立掩膜 (24)

6.7.3应用掩膜 (25)

6.7.4掩膜的反选 (27)

6.7.5掩膜反选后的应用 (28)

6.8主成分分析 (29)

6.9提取蚀变异常信息 (32)

6.10用Surfer软件修饰铁染蚀变异常信息 (34)

6.11总结 (40)

7课后练习及作业 (41)

1 实验目的

通过学习基于地学软件的矿化蚀变异常信息提取的流程，了解利用主成分分析方法提取矿化蚀变异常信息提取的原理，掌握ENVI4.6以及Surfer9软件的基本使用方法。

2 实验内容

2.1 熟悉遥感影像的辐射定标的方法与流程

通过学习遥感影像的辐射定标的方法与流程，掌握用ENVI软件进行辐射定标的具体操作，加深对遥感图像处理知识的理解。

2.2 掌握遥感影像的波段合成、投影转换、影像裁剪的方法

通过学习遥感影像的波段合成、投影转换、影像裁剪的方法与流程，掌握用ENVI软件进行这些操作的具体步骤，加深对遥感图像处理知识的理解。

2.3 掌握ETM+遥感影像的Flaash大气校正、掩膜的应用方法

通过学习ETM+遥感影像的Flaash大气校正，掩膜的建立与应用方法，熟练掌握ENVI软件进行这些操作的具体步骤，加深对大气校正知识的理解。

2.4 掌握ETM+遥感影像羟基和铁染异常信息提取的方法与流程

通过学习利用主成分分析方法对ETM+遥感影像进行羟基和铁染异常信息提取的方法与流程，熟练掌握ENVI软件进行这些操作的具体步骤，加深对矿化蚀变信息提取知识的理解。

2.5 掌握ENVI与Surfer软件协同制图的方法。

通过学习使用ENVI软件提取ETM+遥感影像矿化蚀变信息与Surfer软件绘制矿化蚀变异常空间分布图，综合掌握ENVI与Surfer软件协同制图的方法，加深对所学知识的理解。

3 实验要求

上课时，需认真听取实验老师的细致讲解；练习时，按照实验指导书的步骤，自己亲自动

手使用ENVI4.6软件提取矿化蚀变异常信息；以及在Surfer软件中制作矿化蚀变异常信息空间分布图。

4 实验条件

4.1 软件平台：ENVI4.6、Surfer9

ENVI4.6软件为遥感图像处理软件，用于提取矿化蚀变异常信息；Surfer 9软件为一款以画三维图（等高线,image map, 3d surface）的软件，具有的强大插值功能和绘制图件能力，使它成为用来处理XYZ数据的首选软件，是地质工作者必备的专业成图软件。

4.2 遥感数据源：金川地区Landsat7 ETM+遥感影像

Landsat7于1999年4月15日发射升空后，由于其优越的数据质量，以及与以前的Landsat系列卫星保持了在数据上的延续性，现在已成为中国遥感卫星地面站的主要产品之一，2003.5月出现故障退役。

5 实验原理

5.1 蚀变异常提取的地质依据

目前认为中等强度以上的蚀变带对于TM（ETM+）蚀变信息提取是十分有利的。有时尽管有蚀变岩存在，但不一定有矿，然而围岩的剧烈而较大范围的蚀变常常与大矿及富矿石的生成互为隶属，大型特大型内生热液矿床一般均有强烈且较大范围的围岩蚀变，并且具有分带现象（如斑岩铜矿）。这便是以找

矿（首先是大矿、富矿）为最终

目的蚀变遥感异常提取的地质依

据。

5.2 蚀变异常提取

的物理依据

从Hunt G.RC（1978）和他

领导的实验室的研究成果以及阎

积惠等（1995）依据矿物反射波谱特征吸收谱带特点的定性分类研究中可以知道：主要造岩矿物在可见光——近红外光谱（0.35-0.25um）并不产生具有鉴定意义的反射谱带，其光谱特征主要由岩石中为数不多的次要矿物决定：

一是含铁（Fe2+， Fe3+）基团产生，含铁矿物主要有角闪石、赤铁矿、褐铁矿、针铁矿、磁铁矿、黄钾铁矾等，他们在TMl、TM4波段有强的吸收带，若岩石中含多量的Fe3+，而含

Fe2+较少，这类岩石的主要吸收谱带位于TM4和TM1波段，反射波长相当于TM3波段的电磁波。若含大量的Fe2+、含Fe3+较少，则主要吸收谱带位于1波段，对于波长相当于TM2波段的电磁波有某种程度的反射。

二是含羟基（OH-）、水（H2O）或碳酸根（CO32-）基团产生，羟基的吸收谱带主要有二处：2.2um，2.3um，由于OH-在2. 2—2.3um附近存在强吸收谷（称为羟基谱带），使得TM7产生低值，TM5产生高值，含羟基矿物大多为次级蚀变矿物，如高岭土、叶腊石、云母类矿物、绿泥石、绿帘石等，水在1.4um和1.9um处有特征吸收带。含碳酸根矿物主要有五个特征吸收谱带（1.9-2.55um），较强的两个在2.35um和2.55um波长处（称为碳酸根谱带），相对较弱的在1.9um，2.0um，2.16um三处。

5.2.1 利用主成分分析方法提取矿化蚀变信息

遥感图像各波段间存有一定的相关性，为了减少相关性对分类的影响，常使用主分量分析法去相关。主分量分析基于变量之间的相互关系，在信息总量守恒的前提下，利用线性变换的方法来实现去相关性。由于所获各主分量之间不相关，故各主分量之间信息没有重复或冗余。故主分量分析这一基本性质在蚀变异常信息提取中被充分利用并很有成效。

5.2.2 铁染蚀变异常分析

针对各个波段波谱所反映的特征，选取TM1、TM3、TM4、TM5这4个波段作为组合波段做主成分分析。这是因为铁氧化物的特征光谱信息集中在了TM1—4波段，在TM4和TM1波段有吸收峰，在TM3波段无特征吸收而呈高反射。同时为了避免含羟基和碳酸根矿物的干扰，在选取波段组合时舍弃了TM7波段。

对TM1、TM3、TM4、TM5应用掩膜做主成分分析，统计分析如下：

特征向量（Eigenvectors）