澳大利亚Carpentaria锰矿蚀变信息提取研究.docVIP

澳大利亚Carpentaria锰矿蚀变信息提取研究 　　摘要：根据锰矿相关蚀变矿物光谱特征的分析以及在分析和研究ASTER多光谱卫星遥感数据的基础上，利用ERDAS空间建模的方法，根据锰矿相关蚀变矿物光谱特征的分析对澳大利亚Carpentaria锰矿先做出了数据预处理以及提取蚀变异常信息，主要提取了铁染、羟基、碳酸盐类的蚀变信息 关键词：ASTER；蚀变信息；主成分分析；锰矿 中图分类号：TB 文献标识码：A doi：10.19311/j.cnki2016.25.092 0引言 以澳大利亚Carpentaria地区为例，利用ASTER对该区的矿化蚀变信息提取的研究。蚀变信息提取方法的前提条件：（1）该区的蚀变信息以及与蚀变信息对立的干扰信息的光谱特征；（2）提取的蚀变信息的最优波段组合。先分析该区的成矿条件和蚀变矿物光谱特征，再用主成分分析的方法提取出矿化蚀变信息 1研究背景 该研究区的锰矿有沉积型、热液型、表生型三种锰矿床，软锰矿，硬锰矿及菱锰矿等锰矿。浅层的锰矿床中富含大量羟基离子和三价铁离子是由于在锰矿中存在大量铁化和泥化，使锰矿化异常在一定程度上可由羟基、铁染异常表征。根据菱锰矿的性质，锰矿化异常可以依据碳酸根离子来判断 2ASTER数据简介 ASTER获得的影像数据有14个波段，集中分布于可见光-近红外、短波红外、热红外3个光谱区。主要信息如表2所述 3ASTER数据的蚀变信息波谱依据 3.1羟基离子波谱特征 羟基离子（-OH）在8波段呈现强烈的吸收峰现象，在4波段呈反射峰。利用ASTER数据主成分分析提取羟基蚀变物主要对ASTER1-ASTER（N）（1，3，4，6波段主成分分析方法提取Al-OH基团信息；1，3，4，8波段主成分分析方法提取Mg-OH基团信息），波段的选择取决于研究区蚀变矿物的波谱特征。在波段8为反射峰，系数取正；Band4为吸收谷，取负 3.2Fe3+波谱特征 Fe3+的反射率在ASTER 1-ASTER4波段是递增的状态，第4波段的反射率远大于3波段的反射率，从而在3波段形成一个相对较弱的吸收峰，在4波段有明显的反射峰。1波段和3波段为吸收谷 铁染蚀变矿物在ASTER1、3波段处是吸收谷，在ASTER2、4波段处是反射峰，由此得到对代表铁染物主成分的判断准则是：选取ASTER1、2、3、4等4个波段构成与改该主成分的特征向量。其中Band2系数应与Band1及Band3的系数符号相反，与Band4的系数符号相同，因为铁染物对应的1、3波段处是波谷，因此取负；同时Band2、4取正 3.3CO32-波谱特征 离子在波段8有明显的吸收峰，在波段5有一个次级吸收峰 4蚀变信息提取流程 ASTER数据的蚀变信息提取的基本流程如图1 5ASTER蚀变信息的提取 5.1ASTER数据校正及波段提取 检查并处理ASTER数据，再根据需要选取相应的波段：铁染（ASTER1-4波段）、羟基（ASTER1、3、4、6波段）、碳酸盐化（ASTER1、3、5、8波段） 5.2掩膜去干扰并进行主成分分析 5.2.1制作mask文件 植被在ASTER数据光谱特征上选band4/band3为判断依据：当只有植被因素时，其反射特征呈band4低于band3亮度值应超过band4亮度值，相反，还有其他因素的作用时，band4亮度值才会超过band3亮度值 水体在ASTER近红外波段吸收较强，在第4波段附近吸收率近100%； 在可见光波段的灰度比较高，因此可以选用波段band4/band1　　此步骤采用ERDAS建模进行。建模如图2 运算函数为：EITHER 1 IF （ $n1_prdat015_1234（4） / （$n1_prdat015_1234（1） + 0000000001） = 056 ） OR 0 OTHERWISE 5.2.2掩膜主成分分析 该模型，主要是对数据进行波段选择，然后用生成的mask文件进行掩膜，最后进行主成分分析，经过主成分分析即可得到蚀变异常信息 5.3铁染矿物遥感异常的提取 用主成分分析ASTER数据的1-4波段来提取铁染矿物信息。Fe2+/Fe3+异常分量表征特征为ASTER2、4波段贡献与3波段贡献相反 通过对表中特征向量的观察，可得第3主分量（PCA3）为蚀变异常分量 5.4OH-矿物遥感异常的提取 通过主成分分析1、3、4、6波段，可以提取绝大多数OH-基团矿物信息。OH-物异常分量表征特征应为在ASTER1、3、4波段贡献与6波段贡献相反，并在6波段具有高载荷 通过对表中特征向