下载文档原格式
遥感蚀变信息提取方法嘿，朋友们！今天咱来聊聊遥感蚀变信息提取方法。
这可不是什么高深莫测的东西，就好像你找宝藏，得知道从哪儿开始挖，用啥工具一样。
遥感啊，就像是我们的千里眼，能从老远的地方就看到地面上的情况。
那蚀变呢，就是地面上一些特别的变化，就像你脸上长了颗痣一样明显。
那怎么把这些蚀变信息给提取出来呢？咱可以先从图像入手呀，就跟你看照片找不同一样。
仔细瞧瞧那些颜色、纹理啥的，说不定就能发现点蛛丝马迹。
这图像就像是一幅大地图，你得学会在上面找线索。
然后呢，可以利用一些专业的软件和算法。
这就好比你有了一把神奇的铲子，能帮你更准确、更快速地挖掘出那些蚀变信息。
这些软件和算法可厉害了，它们能把那些隐藏的信息都给揪出来。
你说这难不难？其实也没那么难啦！只要你有耐心，就像钓鱼一样，静静地等着鱼儿上钩。
而且啊，现在科技这么发达，工具这么多，还怕找不到那些蚀变信息吗？比如说，你可以想象一下，在一大片森林里找一朵特别的花，遥感就像是让你在空中俯瞰整个森林，然后软件和算法就是帮你把那朵花凸显出来的魔法。
再比如说，提取遥感蚀变信息就像是在一堆沙子里找金子，你得有好的方法和工具，还得有一双敏锐的眼睛。
还有啊，这可不是一个人能搞定的事儿，得大家一起合作。
就像一场足球比赛，每个人都有自己的位置和任务，大家齐心协力才能取得胜利。
遥感蚀变信息提取方法真的很有趣，也很有意义。
它能帮我们更好地了解地球，发现那些隐藏的宝藏和秘密。
所以啊，大家可别小瞧了它，要认真去学习、去探索。
总之，遥感蚀变信息提取方法就像是一把打开地球秘密之门的钥匙，只要我们掌握了这把钥匙，就能开启无限的可能。
让我们一起加油，去探索这个神奇的世界吧！。
蚀变岩石是在热液作用影响下，使矿物成分、化学成分、结构、构造等发生变化的岩石。
由于它们经常见于热液矿床的周围，因此被称为蚀变围岩，蚀变围岩是一种重要的找矿标志。
利用围岩蚀变现象作为找矿标志已有数百年历史，发现的大型金属、非金属矿床更是不胜枚举：北美、俄罗斯的大部分斑岩铜矿、我国的铜官山铜矿、犹他州的大铝矿、西澳大利亚的大型金矿、墨西哥的大铂矿、美国许多白钨矿、世界大多数锡矿、哈萨克斯坦的刚玉矿等，都属于以围岩蚀变作为找矿标志发现的矿床。
国内外遥感工作者，都在不断地设计、研制和总结对这种遥感信息的提取和识别技术。
矿化蚀变信息是找矿的一个重要标志，而这些对找矿有指导意义的矿化蚀变信息常常受其它地物信息的干扰，和受遥感图像的波谱分辨率和空间分辨率的制约，往往表现的很微弱。
因此，国内外学者也在不断尝试各种技术方法提取这种矿化蚀变弱信息。
本文总结了遥感蚀变信息提取的各类方法，及其在ENVI软件中的实现。
•原理遥感技术主要是建立在物体反射和发射电磁波的原理之上。
而地物波谱特性通常都是用地物反射辐射电磁波来描述。
由于地物反射发射电磁波的特性不同，其反射波谱曲线形态也有千差万别。
如植物的反射波谱曲线上，在绿光波段表现由于其叶绿素的存在表现为有一强反射峰，而在短波红外波段由于叶冠组织的相互作用表现为强反射峰，在红光波段则表现为强吸收谷。
遥感地质应用中，近矿围岩蚀变形成的蚀变岩石与其周围的正常岩石在矿物种类、结构、颜色等方面都有差异，这些差异导致了岩石反射光谱特征的差异，并且在某些特定的光谱波段形成了特定蚀变岩石的光谱异常。
光谱异常为用遥感图像的异常信息提取提供了理论依据。
•方法及实现依据矿化蚀变岩与围岩的波谱特征的差异，可采用图像增强处理方法获取矿化蚀变信息增强的图像变量，从而最终实现提取矿化蚀变信息的目的。
一般图像增强突出蚀变信息有以下几种方法。
(1)蚀变干扰信息剔除遥感数据包含地表的信息，遥感在地质方面的应用就是提取用户需要的信息，提取矿化蚀变信息的过程是计算影像中所有像素信息统计归类分析的过程，蚀变异常信息的提取对遥感图像的质量要求较高，因此首先要对遥感数据进行严格的筛选，干扰噪声小的数据，一般要求遥感数据的时相是植被发育较弱、冰雪覆盖少的季节，同时该时相的云覆盖量较少。
遥感构造蚀变异常信息提取及找矿预测——以老挝为例刘文兰;张微【摘要】遥感地质异常信息提取是遥感地质学的重要研究内容之一.以老挝为例,对研究区ETM+遥感图像采用几何纠正、图像镶嵌、彩色合成、空间变换、图像增强等技术手段,提取地质构造信息；利用掩模技术去除干扰信息,采用主成分分析和阈值分割技术进行铁染和羟基蚀变异常信息提取；结合区域地质、地球物理、地球化学等多元地学信息综合分析,开展研究区遥感找矿预测,圈定了5处遥感找矿远景区.研究成果为加快老挝地质找矿工作提供了基础数据和决策依据.%The extraction of remote sensing geological anomaly information is one of the important research methods in GARS. In this paper, taking Laos as an example and using such means as geometric rectification, image mosaicking, color composite, space transformation and image enhancement processing to the ETM + image in the study area, the authors extracted the geological structural information. The interference information was removed by mask technique, and the iron staining and hydroxyl alteration information was extracted by using principal component analysis and threshold segmentation techniques. Five perspective areas for mineral exploration were delineated on the basis of a synthetic analysis of the multi - geoscience information. The results can provide a reference and effective means for accelerating the geological prospecting work in Laos.【期刊名称】《国土资源遥感》【年(卷),期】2012(000)002【总页数】7页(P68-74)【关键词】断裂构造;蚀变异常;图像增强;找矿预测;老挝【作者】刘文兰;张微【作者单位】中国地质大学(北京)地球科学与资源学院,北京100083;中国国土资源航空物探遥感中心,北京100083【正文语种】中文【中图分类】TP75随着矿产资源勘查工作的不断深入，易于发现的露天矿、浅部矿和已知类型矿床日趋减少，找矿难度越来越大［1］。
专题四：基于crosta方法的遥感矿物蚀变信息提取1、专题概述蚀变岩石是在热液作用影响下，使矿物成分、化学成分、结构、构造等发生变化的岩石。
由于它们经常见于热液矿床的周围，因此被称为蚀变围岩，蚀变围岩是一种重要的找矿标志。
遥感地质应用中，近矿围岩蚀变形成的蚀变岩石与其周围的正常岩石在矿物种类、结构、颜色等方面都有差异，这些差异导致了岩石反射光谱特征的差异，并且在某些特定的光谱波段形成了特定蚀变岩石的光谱异常。
光谱异常为用遥感图像的异常信息提取提供了理论依据。
遥感蚀变异常信息提取的方法有多种，其中主成分分析法提取蚀变信息是相对最为广泛的。
主成分分析（PCA：Principal ponent Analysis）是基于信号二阶统计特性的分析方法，由于所获各主成分之间不相关，主成分之间信息没有重复或冗余。
多光谱遥感数据通过PCA 所获每一主成分常常代表一定的地质意义，且互不重复，即各主成分的地质意义有其独特性。
但是由于蚀变矿物形成的影像特征在遥感图像上往往表现得很微弱或不明显，甚至“淹没”在主体色调中。
以TM数据为例，通过Crosta方法说明其准则，通过TM1、TM3、TM4、TM5和TM1、TM4、TM5、TM7的波段组合分别进行主成分分析提取铁染蚀变和羟基蚀变信息。
由TM1、TM3、TM4、TM5做PCA处理，处理后的某个新的成分可能集中了铁染蚀变信息。
对代表铁染蚀变的主成分的判断准则是：TM3的系数应与TM1、TM4的系数相反。
由TM1、TM4、TM5、TM7作为输入波段进行主成分分析。
对代表羟基和碳酸根离子主成分的判断准则是：TM5系数应与TM7、TM4的系数符号相反，TM1一般与TM5系数符号相同。
依有关地物的波谱特征，羟基和碳酸根离子信息包含于符合这判断准则的主成分。
铁染蚀变和羟基蚀变存在于绝大多数成矿岩体中，提取这两种蚀变信息基本可以确定研究区成矿岩石的分布情况。
原理及其他方法参考博文：《遥感矿物蚀变信息提取方法及ENVI下实现》blog.sina../s/blog_764b1e9d0100so8i.html参考文献：《干旱区TM图像蚀变信息提取方法研究》、《ETM+(TM)蚀变遥感异常提取方法研究与应用--方法选择和技术流程》2、处理流程介绍主要包括：图像预处理、PCA分析、选择有效成分、异常切割、蚀变信息处理等步骤（图1）。
“基础地质学”创新性实验遥感地质蚀变异常信息提取实验“基础地质学”国家级教学实验示范中心二〇二二年四月二十八日目录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、Surfer9ENVI4.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）从上表可以看出，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波段，在TM7波段为特征吸收带，在TM5相对高反射。
主成分变换的波段组合的选择也是以此为依据的，由于可见光波段对铁氧化物敏感，为了避免铁氧化物信息的干扰，故只选择了一个可见光波段参与运算。
对TM1、TM4、TM5、TM7应用掩膜做主成分分析，统计分析如下：特征向量（Eigenvectors）TM4波段的信息；PC3反映的是TM1波段的信息；PC4反映的是TM5和TM7波段的信息，且符号相反。
根据含羟基类矿物和含CO32-矿物的波谱特征，PC4图像中的高亮度信息就是表征含羟基类矿物和含CO32-矿物的蚀变异常信息。
6 实验步骤6.1 金川地区ETM+遥感影像辐射定标在ENVI 4.6软件中，单击菜单栏中的file——open image file工具，打开金川地区Landsat7 ETM+遥感影像中band1影像文件。
然后，单击菜单栏中的Basic Tools——preprocessing——Calibration Utilities——Landsat TM工具，为金川地区ETM+遥感影像中band1影像文件进行辐射定标，如图 6-1所示。
图 6-1打开辐射定标工具对话框在TM Calibration Input File弹窗下的Select Input File栏中选择Landsat 7的band1文件，然后单击OK，如图 6-2所示。
图 6-2 TM Calibration Input File对话框接来下，在TM Calibration Parameters弹窗中，单击Get Calibration Parameters from Web选项，如图 6-3所示。
图 6-3 TM Calibration Parameters对话框在TM Web Calibration Parameters弹窗中，选择Date/Path/Row/Band选项，根据金川地区Landsat7 ETM+遥感影像的头文件信息，依次填写下面红色标注处的信息，如图 6-4所示。
注意：每一景影像的信息都不一样，一定要根据头文件信息填写。
图 6-4 TM Web Calibration Parameters对话框填写完成后，在TM Calibration Parameters中依次填写下图红色标注部分的信息，如图6-5所示。
图 6-5 TM Calibration Parameters对话框单击OK，即完成了金川地区Landsat7 ETM+遥感影像中band1影像文件的辐射定标操作；接下来，依次对band2、band3、band4、band5、band7波段进行辐射定标操作。
6.2 金川地区ETM+遥感影像不同波段的合成首先，在ENVI 4.6软件中，单击菜单栏中的file——open image file工具，打开经过辐射定标的金川地区ETM+遥感影像的band1、band2、band3、band4、band5、band7数据文件，如图 6-6所示。
