基于crosta方法的遥感矿物蚀变信息提取操作文档.doc

《基于crosta方法的遥感矿物蚀变信息提取操作文档.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、......专题四：基于crosta方法的遥感矿物蚀变信息提取1、专题概述蚀变岩石是在热液作用影响下，使矿物成分、化学成分、结构、构造等发生变化的岩石。由于它们经常见于热液矿床的周围，因此被称为蚀变围岩，蚀变围岩是一种重要的找矿标志。遥感地质应用中，近矿围岩蚀变形成的蚀变岩石与其周围的正常岩石在矿物种类、结构、颜色等方面都有差异，这些差异导致了岩石反射光谱特征的差异，并且在某些特定的光谱波段形成了特定蚀变岩石的光谱异常。光谱异常为用遥感图像的异常信息提取提供了理论依据。遥感蚀变异常信息提取的方法有多种，其中主成分分析

2、法提取蚀变信息是相对最为广泛的。主成分分析（PCA：PrincipalComponentAnalysis）是基于信号二阶统计特性的分析方法，由于所获各主成分之间不相关，主成分之间信息没有重复或冗余。多光谱遥感数据通过PCA所获每一主成分常常代表一定的地质意义，且互不重复，即各主成分的地质意义有其独特性。但是由于蚀变矿物形成的影像特征在遥感图像上往往表现得很微弱或不明显，甚至“淹没”在主体色调中。以TM数据为例，通过Crosta方法说明其准则，通过TM1、TM3、TM4、TM5和TM1、TM4、TM5、TM7的波段组合

3、分别进行主成分分析提取铁染蚀变和羟基蚀变信息。由TM1、TM3、TM4、TM5做PCA处理，处理后的某个新的成分可能集中了铁染蚀变信息。对代表铁染蚀变的主成分的判断准则是：TM3的系数应与TM1、TM4的系数相反。由TM1、TM4、TM5、TM7作为输入波段进行主成分分析。对代表羟基和碳酸根离子主成分的判断准则是：TM5系数应与TM7、TM4的系数符号相反，TM1一般与TM5...专业........系数符号相同。依有关地物的波谱特征，羟基和碳酸根离子信息包含于符合这判断准则的主成分。铁染蚀变和羟基蚀变存在于绝大多数

4、成矿岩体中，提取这两种蚀变信息基本可以确定研究区成矿岩石的分布情况。原理及其他方法参考博文：《遥感矿物蚀变信息提取方法及ENVI下实现》blog.sina..cn/s/blog_764b1e9d0100so8i.html参考文献：《干旱区TM图像蚀变信息提取方法研究》、《ETM+(TM)蚀变遥感异常提取方法研究与应用--方法选择和技术流程》2、处理流程介绍主要包括：图像预处理、PCA分析、选择有效成分、异常切割、蚀变信息处理等步骤（图1）。原始数据经过一系列的图像预处理，包括几何纠正，大气校正，掩膜去背景等；再选择波

5、段进行主成分分析，根据上述的原则进行成分的选择；将选择的成分用异常切割标准处理，突出蚀变信息；并与已知的矿点资料进行叠合分析，验证提取的蚀变结果。在做异常切割时，利用(X+kó)确定异常下限和划分异常强度等级。X是某一成分的统计均值代表区域背景，ó是该成分的标准差。k取值一般是（1～3），有了这一标准，切割异常时可以减少主观任意性，并使操作较为规化。经PCA处理后得到的各成分是互不相关的，理论上可以看成是单一的地物。但是由于大气，传感器、异物同谱等影响，各成分并不只是单一的地物，所以在提取异常时利用阈值切割法。参考《

6、ETM+(TM)蚀变遥感异常提取方法研究与应用--方法选择和技术流程》...专业........TM/ETM+原始数据数据读取辐射定标几何校正大气校正研究区裁剪基准影像蚀变异常等级分布图结果验证1345PCA分析4铁染异常成分判断数据预处理铁染蚀变信息提取规则裁剪1457PCA分析羟基异常成分判断异常等级划分异常等级划分羟基蚀变信息提取图专4-1流程图1、详细的处理过程本专题的数据存放在“18-基于crost方法遥感蚀变信息提取”文件夹。...专业........3.1数据预处理第一步：数据读取和定标（1）打开tm原

7、始影像数据选择主菜单->File->OpenExternalFile->Landsat->GeoTIFFwithMetadata，打开“XXXXLT249IKR01L5145028__MTL.txt”文件。可以看到ENVI自动进行了波段合成。（2）TM数据辐射定标选择主菜单->BasicTools->Preprocessing->CalibrationUtilities->LandsatCalibration，选择含有多波段的L5145028__MTL.txt，弹出如下对话框，图2，选择保存路径，点击OK。图专4

8、-2辐射定标参数设置对话框...专业........第二步：FLAASH大气校正（3）辐射量度单位转换选择主菜单->BasicTools->BandMath，在对话框中输入转换公式，见下图3：图专4-3Bandmath输入单位转换公式输入转换公式以后，点击AddtoList，再点击OK。弹出如下对话框（图4），点击‘MapVariabletoI