矿化蚀变信息提取的TM_ETM_遥感影像模式_张廷斌

  • 格式:pdf
  • 大小:173.58 KB
  • 文档页数:5

矿化蚀变信息提取的TM /ETM +遥感影像模式张廷斌①,②,唐菊兴③,黄丁发①(①西南交通大学土木工程学院,成都610031;②成都理工大学地球科学学院,成都610059;③中国地质科学院矿产资源研究所,北京100037)收稿日期:2008-06-10 修订日期:2008-07-24基金项目:国家科技支撑计划项目(2006BAB01A01);四川省教育厅自然科学青年项目(08ZB085);成都理工大学青年科学基金项目(2008Q J10)。

作者简介:张廷斌(1978年~),男,山东沂水人,讲师,西南交通大学在读博士研究生,主要从事3S 技术教学与研究工作。

E -m ail :zhangtb @cdut .edu .cn摘要:在遥感找矿中,矿化蚀变信息的提取是重要的工作内容之一。

以TM /ET M +为数据源进行矿化蚀变信息的提取,经历了一个由简单到综合的发展历程。

为更好地开展蚀变信息的提取工作,本文在总结前人研究的基础上,将各种蚀变信息的提取方法归结为两个模式,即“经验模式”和“准定量模式”,并对两种模式下的各种方法结合实例进行了介绍,最后对两种模式的关系作了初步的探讨。

关键词:蚀变信息;提取;TM ;ET M +;遥感;模式中图分类号:T P75 文献标识码:A 文章编号:1000-3177(2009)102-0047-051 引 言矿化蚀变信息是重要的找矿标志,最常见的蚀变有硅化、绢云母化、绿泥石化、褐铁矿化、黄铁矿化、云英岩化、矽卡岩化、白云岩化及碳酸盐化等。

由于蚀变岩和非蚀变岩在矿物成分和结构上的差异,形成了遥感影像上的光谱差异,用遥感技术提取蚀变信息是遥感找矿的方法之一。

三十多年来,国内外学者尝试了许多方法,取得了丰硕的成果。

本文在研究前人工作的基础上,总结了TM /ETM +遥感影像蚀变信息的提取模式。

2 蚀变信息提取模式概述早在1977年Row an 等[1]采用MSS4/5、MSS5/6和MSS5/7比值合成图像对次生粘土、氧化硅及褐铁矿区特征的热液蚀变区进行识别填图;H unt 等[2](1978年)利用羟基在可见光到近红外的波谱响应特征来确定特定的矿物以圈定蚀变岩石;Ambrams (1984年)利用TM 图像波段比值+主成份变换方法排除了植被覆盖干扰提取了蚀变信息;C roksta 等[3](1989年)提出选择主成份分析;Lo ug hlin [4]指出,有目的地仅对一定的波段组合进行主成份分析可将特定的信息聚集到单一的主成份中;Robert W .M 等[5](1996年)在新墨西哥南部Blue G reek 盆地内热液蚀变填图使用TM 图像数据成功地圈定了粘土、云母含羟基矿物的矿化蚀变区;Timothy K M 等(2002年)利用TM 5/TM7、TM5/TM 1、(TM 5/TM4)×(TM 3/TM 4)进行假彩色合成,在阿拉伯—努比亚地盾的干旱气候下提取金矿化蚀变信息。

国内地质矿产遥感信息提取技术的应用研究始于上个世纪80年代。

刘燕君[6](1983年)利用M SS 图像采用波段比值方法在北京密云水库周围的丘陵区变质铁矿实验区,成功地把矿化点及含矿围岩内强烈体而部位增强出来;赵元洪等[7](1991年)在浙江新昌地区采用“TM 4/3+主成份变换+密度分割+消除噪声”压抑植被干扰达到提取蚀变信息的目的;何国金等[8](1994年)提出了“微差信息处理”这一直接提取金矿化信息的方法;1994年马建文等[9]提出了采用“图像掩膜+波段选择主成份变换(SPCA )+人机交互与PC 组分数据二维散点图掩膜分类”方法,有效地减少了多种环境因素干扰下环境与目标之间光谱信息混合,目标信息得到增强;陈赶良等[10](1996年)通过对岩矿的光谱特征、地球化学因子、TM 波段比值因子与地区化学因子的相关性等研究,探索从室内光谱数据中提取与金矿化密切相关的蚀变信息;1997年马建文[11]提出了“TM 掩膜+主成份变换+分类识别”提取矿化弱信息的方法;1999年张远飞等[12]利用“多元数据分析+比值+主成份变换+掩膜+分类(分割)”的方法在新疆、内蒙古及江西、云南成功地提取了金矿化蚀变信息;以及刘素红[13](2000年)、王晓鹏[14] (2003年)、张玉君[15](2003年)、杨波[16]等均进行了相关的研究。

利用TM/ETM+遥感图像提取蚀变信息的方法,经历了一个由简单到繁杂、由单一到综合的发展历程。

为更好地开展蚀变信息的提取工作,有必要对这些方法进行系统的归纳和总结。

在此,笔者将这些提取方法归结为两种模式,即“经验模式”和“准定量模式”。

用“模式”化的思想来理解和选择蚀变信息提取方法,同时也可以用“模式”化的思想来丰富和发展这一方法。

3 经验模式20世纪70年代,H unt和他领导的实验室对300个粒状矿物的测定结果显示,在可见光—近红外光谱区(0.325μm~2.5μm),主要造岩矿物中各主要化学成分(Si、Al、M g和O)并不产生具有鉴定意义的反射谱带;在岩石反射谱带中占据主导地位的是岩石中为数不多的次要矿物(含铁矿物及蚀变矿物)中的Fe2+、Fe3+、OH-、CO2-3等离子或离子基团,形成了反射谱的特征吸收谷[17~19]。

蚀变岩石因含有这些离子(或离子基团)而在TM/ETM+相应波段产生特征光谱。

经验模式就是利用这些典型地物的(矿物)光谱特征(曲线),根据地物在不同波段上的反射、吸收特征,在已有地质资料的基础上,通过遥感图像的信息增强技术,反复比较试验而总结出来的一种模式。

这种模式缺乏精确的工作区地物波谱数据,或者说是在未进行工作区波谱测试的条件下进行的一种蚀变信息提取工作。

这种模式提取的效果因工作区的地物环境和遥感图像的质量以及蚀变的大小而产生比较大的差异。

但这种经验模式的提取工作成本低、操作简单,在一些基岩裸露区或植被低覆盖区可以定性的评判蚀变的范围和种类,往往是一种行之有效的方法。

3.1 波段比值法波段比值处理是将一个光谱波段中的灰度值与另一个波段图像中对应像元灰度值相除,比值的结果反映了地物波谱曲线变化的斜率,从而增强了地物波谱特征的微小差异。

波段比值法是早期遥感蚀变信息提取的主要方法之一,如赵元洪[20](1991年)、平仲良等[21](1993年)利用TM波段比值法提取蚀变信息;陈赶良等[22](1996年)利用岩矿的成分因子与波段比值进行相关性分析,选择研究区最优的波段比值进行蚀变信息的提取。

识别热液蚀变常用的波段比值有:TM3/1,用于识别褐铁矿;TM5/4,用于区分有植被和无植被覆盖的土壤和岩石,区分云母及黄钾铁矾,明矾石及石膏、方解石及粘土这三类矿物,识别褐铁矿化; TM5/7,识别含羟基矿物,水合硫酸盐和碳酸盐; TM7/4,区分云母、石膏与明矾石;TM3/4,识别植被和区分褐铁矿化岩石[23]。

3.2 主成份分析法主成份分析也称为K-L(Karhunen-Lo eve)变换,它是在对多波段图像进行特征统计基础上的多维正交线性变换。

其本质是对多光谱空间的坐标系按一定规律进行旋转,从而达到去除相关、进行特征提取和数据压缩的目的。

主成份分析法也是遥感蚀变信息提取中应用较早的方法之一,且到目前为止其应用也比较广泛,并日渐成熟。

通常进行主成份分析是把一幅图像的所有波段一起处理,但也可以基于地物的光谱特征把所有波段先进行分组,然后进行主成份分析。

Cros-ta和M C.M.Mo ore提出选择主成份分析方法(SP-CA)[3],选择有地质意义的波段进行主成份分析。

Loughlin(1991年)的研究表明,有目的地对一定的波段组合进行主成份分析可将特定的信息聚集到单一的主成份中[4],提出了面向特征的主成份分析。

国内如赵元洪[20]、马建文[24]、张满郎[25]、张玉君[26]等都对主成份分析的方法进行了改进和丰富。

3.3 综合提取法蚀变信息相对于遥感图像丰富的信息而言是一种弱信息,且常被背景信息所掩盖,通常在蚀变信息的提取过程中,单纯用一种或几种方法很难取得满意的效果。

综合提取法就是各种蚀变信息提取方法的综合应用,一般是在TM原始波段、比值波段、和主成份分析后的主分量中进行的二次波段比值或二次主成份分析。

这种方法有机地综合了原始波段、比值波段和主分量波段的信息,可以定向的提取感兴趣信息,同时压抑或去除掉干扰信息。

在实践中,这种方法应用较多,国内许多学者也取得了一些值得借鉴的经验。

如:张玉君等[26]利用TM1+TM2、TM4/3、TM5、TM7波段组合代替TM1、TM4、TM5、TM7组合进行主成份分析,提取羟基蚀变信息;赵元洪[20]和刘志杰[27]利用TM5/7和TM4/3进行主成份分析,在提取羟基蚀变的同时压抑了植被信息;朱小鸽[28]利用不同原始波段组合以及主成份分析波段、原始波段和比值波段的组合进行主成份分析,提取各种相关地质信息;陈松龄等[29]利用TM1-TM5,TM7波段主成份分量的比值及波段比值的加减运算识别围岩蚀变信息;孟新[30]利用TM7(R)、TM5(G)、TM4(B)进行H IS变换后的饱和度(S)分量与TM5/1波段比值进行主成分分析提取褐铁矿化信息;还有其他研究者如张远飞[12]、张满郎[25]、张晋开[31]等,均取得了良好的效果。

4 准定量模式波谱测试研究是遥感图像蚀变信息提取的物理基础,研究区地面波谱测试工作对揭示研究区遥感图像色调异常的成像机理和地质成因,对定向的选取蚀变信息提取方法和提取结果的解释等方面都具有重要的指导意义。

基于地面波谱测试结果开展的蚀变信息提取工作,因理论依据充分、目标明确、针对性强而带有“定量”的色彩。

由于光谱分辨率和空间分辨率的限制,利用TM/ETM+图像一般仅能识别和区分铁染异常(铁的氧化物和氢氧化物)和羟基异常(含羟基矿物、水合硫酸盐及碳酸盐等)两大类,因此称基于研究区典型地物实测光谱数据进行的蚀变信息提取方法为“准定量模式”。

4.1 基于视反射率的蚀变信息提取由于卫星传感器所记录的数据是图像的相对亮度(DN),并非辐射亮度值L(反射的太阳辐射能量),因此需将各TM波段的图像相对亮度值(DN)转换为辐射亮度值L,再将辐射亮度值转换为TM 图像的视反射率。

闫积惠等[32]讨论了TM图像的相对亮度值与岩石反射波谱的关系,指出了图像相对亮度值与岩石反射波谱的相似性;王海平等[33]在西支沟地区以地面实测波谱数据为基础,建立了典型矿区TM单波段图像视反射率与地面波谱反射率的相关模式,用于矿化蚀变信息的反演识别研究。

4.2 基于多元数据分析的蚀变信息提取该模式主要以研究区典型地物的光谱分析或其与岩石地球化学分析相结合为基础,采用多元数据分析的方法来考查TM图像及其变换图像(波段比值、主成份分析的向量等)、岩石地球化学数据和蚀变目标间的相关关系,然后采用多元统计回归的方法建立多元变量的回归方程进行蚀变信息的准定量提取。