Landsat8OLI遥感影像的组合分析法对铁染蚀变的提取研究报告

1、-Landsat8OLI遥感影像的组合分析法对铁染蚀变的提取研究-以毛坪铅锌矿为例马威，温兴平，石振杰，攀(1. 理工大学国土资源工程学院， 650093；2.省矿产资源预测评价工程实验室, 650093)摘要：本文以毛坪铅锌矿为研究区域，利用最新的Landsat8 OLI影像，在详细解读毛坪的Landsat8 OLI影像的根底上，以蚀变围岩的光谱异常为理论依据，结合研究区的地质资料，根据矿物的光谱特征，利用掩膜技术去除水体、阴影等干扰信息，采用主成分分析为主、比值法为辅的方法利用OLI卫星数据提取研究区铁染异常蚀变信息，并用ETM+影像的主成分分析作为比对，研究说明，利用组合分析法，也就是在

2、主成分分析的根底上，采用比值法辅助，对Landsat8 OLI遥感影像上的铁染异常蚀变信息有较好的提取效果。关键词：OLI影像;蚀变信息;提取方法;毛坪铅锌矿. z-引言围岩蚀变wall-rock Alteration,又称围岩交代蚀变或主岩交代蚀变，是指近矿围岩在流体的作用下发生的化学和物理变化，从而引起围岩化学成分和构造构造的变化1。蚀变围岩因其矿物成分、岩石构造与正常围岩有所不同, 所以二者在遥感图像表现为不同的色调, 在波谱维表现为截然不同的波谱曲线 2-3。这为利用遥感影像提取围岩蚀变信息提供了根底。到目前为止，围岩蚀变信息的提取多利用ETM+,ASTER等影像数据。2021年2 月

3、，NASA 成功发射了Landsat 8 卫星，Landsat 8既保持了Landsat 7 卫星的根本特征，又在波段数量、波段光谱围和影响的辐射分辨率方面有了改进。2021年，*涵秋基于Landsat 8 的首幅影像，分析和研究了Landsat 8 遥感影像的新增特征及其生态环境意义4。旭文利用Landsat 8 卫星上运营性陆地成像仪OLI 获取的太湖流域遥感影像数据，开展了其在大气能见度、蓝藻水华及建立用地等生态环境监测领域中的应用研究5。OLI 数据在地质遥感找矿方面的应用还没有文献报道。此外，基于多光谱数据进展矿化蚀变信息提取的方法主要有主成分分析法、比值方法、光谱角法等6-7。这些

4、方法在大量的实践应用中均表现出了不同的问题 8，主成分分析法依据蚀变围岩与非蚀变围岩波谱差异提取矿化蚀变异常时, 使得提取微弱的蚀变围岩信息淹没于背景噪声中, 比值方法在波段相除时减弱源数据的信噪比,光谱角法因多光谱数据谱带过宽、数目偏少而导致结果精度偏低。针对以上问题，本研究在详细解读毛坪的OLI影像的根底上，以蚀变围岩的光谱异常为理论依据，提取遥感蚀变信息。结合研究区的地质资料，根据矿物的光谱特征，采用主成分分析为主、比值法为辅的方法利用OLI卫星数据提取研究区蚀变信息，并利用Landsat 7ETM+影像的蚀变信息提取结果进展比照分析。1研究区概况毛坪位于省东北部的地区，地形上属云岭高原

5、与盆地的交接处，典型的山地地貌，地势南高北低，地形切割强烈，地理坐标是N26.80-28.50，E102.88-105.309。矿区出露的地层有上泥盆统宰格组D3zg，下石炭统大塘组C1d、下石炭统摆佐组C1b、中石炭统威宁组C2w，下二叠统梁山组P1l、下二叠统栖霞茅口组P1q+m，及第四系Q图1。矿体赋存于上泥盆统宰格组第二段D3zg2和下石炭统摆佐组上段C1b2及中石炭统威宁组C2w10。主要有容矿层的白云岩的退色、结晶和白云石化、黄铁矿化、重晶石化、方解石化、硅化9。各种不同的围岩蚀变作用大致具有一定分布规律：即黄铁矿大多出现在矿体以或近矿强烈蚀变的中粗晶白云岩中，显见它和矿化关系十分

6、密切；铁白云石和铁方解石化常伴随矿体出现，一般矿体都与铁白云石和铁方解石化有关，离矿体越近白云石结晶越粗大，方解石团块越发育。2数据来源近年来，遥感技术和地理信息技术飞速开展，2021年，美国地质调查局USGS及太空署NASA发射了“陆地卫星数据连续性任务卫星Landsat Data Continuity Mission卫星，发射后更名为Landsat-8。Landsat-8卫星上携带有两个主要载荷：OLI(陆地成像仪)和TIRS(热红外传感器)。为了防止大气吸收特征，还在波段数量、波段光谱围和影响辐射分辨率方面有了改进。此外，还新增了两个波段，蓝色波OLI1段可应用于海岸带观测，短波红外段O

7、LI9包括水汽强吸收特征，可用于云检测5。OLI成像仪包括9个短波谱段波段1波段9，幅宽185km，其中全色波段地面分辨率为15m，其他谱段地面分辨率为30m。如表1所示。表1OLI卫星数据波段Table1 OLI-band satellite data序号波段空间分辨率m10.4330.453 3020.4500.515 3030.5250.600 304567890.6300.680 30 30 30 301530 OLI成像波段在波谱围上做出了很大的调整。OLI成像仪的波段50.845m0.885m去除了波谱围中波长为0.825附近的水蒸气吸收段；波段波谱围变窄，这将有助于减少提取找矿信

8、息时碳酸盐岩对羟基的干扰；另外，为了使得有植被的地表和没有植被的地表之间的反差更大，OLI成像仪的全色波段波段8的波谱围也很窄5。3蚀变提取的的理论依据遥感数据的事变信息是依据电磁波与出露的岩石和地质体发生作用，产生岩石和地质体的特征光谱来提取的11-12。研究区黄铁矿化大多出现在矿体以或近矿强烈蚀变的中-粗晶白云岩中，显见它和矿化关系十分密切铁离子矿物的提取具有一定意义。含铁Fe3+、Fe2+矿物，铁矿物以次生氧化物为主，局部作为热液蚀变带的原生矿物，如常见的褐铁矿、针铁矿Goethite、赤铁矿、黄钾铁矾等含大量Fe3+，也有少量Fe2+的纤铁矿等氧化物62。0.45m，0.55m，0.8

9、5m，0.90m，0.94m存在较强的吸收谷。一般来讲，浅部黄铁矿会由于氧化作用形成褐铁矿化，也就是说，研究区的铁染蚀变在遥感影像上读取为褐铁矿化13。褐铁矿化实际上并不是一种单独矿物，而是由水针铁矿、针铁矿、水纤铁矿、纤铁矿、更富含铝的氢氧化物、水的氢氧化铁胶凝体、以及含水的泥质、氧化硅等常共同产出而形成；含铁量达30% 40%14。其主要成分的光谱特征含量如以下图，在可见光及近红外区中，岩石特征谱带形成占有优势。在在m间产生强的吸收谱带，因而诊断性吸收峰位于0.9m处。通过野外地质调查和岩性地层剖面测量，在典型岩性段建立地质观察点，采集具有代表性的岩石，保证样品有风化面和新鲜面，用GPS定

10、位采样点位置并拍摄野外采样点照片，同时记录岩石特征、地质地貌特征以及环境参数；然后对采样点进展野外光谱测试，每个采样点测试三条光谱曲线，以备后期实验室数据处理减小误差，最后，对每个观察点的光谱曲线取平均值，所得的波谱曲线图如以下图，反映在OLI影像上，OLI2与OLI5波段上具有明显的吸收谷，在OLI4波段具有相对较高的反射率图1 主要铁染蚀变矿物及对应波段根据USGS-MIN波谱库Fig.1The main alteration minerals iron staining and the corresponding band图2野外铁染蚀变矿物及对应波段Fig.2 Field iron s

11、tained alteration minerals and the corresponding band4蚀变信息提取4.1研究区干扰信息的去除由于研究区地物光谱、遥感成像时间以及辐射分辨率的限制, 传感器很难准确地记录复杂的地表信息, 因而误差就不可防止地存在于遥感数据获取的过程中, 误差的存在降低了遥感数据的质量,.所以, 在数据处理、信息提取前, 需要对原始数据进展预处理15。受研究区的地形地貌及数据时相等原因, 要采用掩膜15-16的方法去除研究区围的水体、植被等干扰元素的影响,，使蚀变提取的结果更加准确。掩膜前Before Mask掩膜后After Mask图3掩膜处理前后比照图F

12、ig.3Mask before and after parison chart4.2蚀变提取及结果主成分分析法17-19很好地去除了多波段间的相关性, 使得变换后各主分量之间相互独立, 分别代表不同的地质意义, 并且互不重复, 矿物的蚀变信息集中于*一个主分量中。主成分分析法虽然可以降低波段间的相关性和冗余性, 但是干扰地物, 以及乘性噪声( 局部大气衰减和地形差异等)的存在, 影响异常信息提取, 因此采用比值法20作为辅助手段，提取蚀变信息。与铁离子有关的蚀变矿物主要是黄铁矿，由于野外的风化剥蚀，在遥感影像上，其表现为褐铁矿化。特征光谱主要集中在OLI2-OLI5波段，在OLI4波段表现为相

13、对高反射。根据蚀变矿物的波普曲线及研究区地质概况，确定蚀变提取方法为主成分分析法，在主成分分析中，为防止含羟基和含碳酸根矿物的干扰，舍弃了OLI7波段。根据研究区蚀变矿物波普特征，OLI6/OLI5比值扩大了OLI6与OLI5的光谱反差, 具有增强铁化蚀变图像的亮度反差功能。用OLI6/OLI5代替OLI6, 即应用OLI2、OLI4、OLI5、OLI4/OLI3作为输入波段组合进展主成分分析,以增强铁染蚀变信息。得到主成分分析的特征向量矩阵(表1)，表2OLI2、4、5、4/3主成分特征向量表Table 2 OLI (2,4,5,4 / 3) main ponent feature vect

14、or主成分OLI2OLI4OLI5OLI4/3PC1-0.139400-0.369006-0.918913-0.000156PC2-0.488834-0.7813770.3879330.000318PC3-0.8611670.503273-0.0714590.001117PC4-0.0007850.000868-0.000060.999999从表中可以看出，PC3上OLI5特征向量为负，OLI6/5特征向量为正，因此对PC3图像上进展主成分分析提取蚀变异常，并利用同研究区的ETM+影像的蚀变信息提取结果进展比照，其结果如下:图4OLI铁染蚀变异常分布图Figure 4 OLI iron sta

15、ined alteration anomaly maps图5 ETM+铁染蚀变异常分布图取值波段为1、3、4、5Figure 5 ETM + iron stained alteration anomaly maps (Value band of 1,3,4,5)图中的红色区域为铁染蚀变异常与ETM+比照的情况，通过比对可以看出，采取主成分分析法为主、比值法为辅的方法，能较好地提取铁染蚀变；采用主成分分析法，能较好地提取羟基蚀变。5结论1通过分析地质体波谱特征和遥感蚀变信息提取的地质依据，依据矿区矿化蚀变类型发育特征和各类典型蚀变矿物的波谱特征，提取研究区的铁染蚀变信息。综合比较各蚀变信息提取方

16、法，结合研究区实际情况采用主成分分析法为主，比值法为辅进展遥感蚀变信息提取。其中，采用OLI2、OLI4、OLI5、OLI4/OLI3波段组合进展主成分分析获取铁染蚀变信息，采用OLI2，OLI5，OLI6,OLI7波段组合进展羟基蚀变信息提取。2主成分分析法、比值分析法是现今遥感蚀变信息提取中广泛应用的方法。主成分分析法, 需要将干扰地物去除才能取得较好的效果。而根据研究区实际情况，对铁染蚀变信息而言，干扰地物的影响很难消除，因此在主成分分析的根底上，采用比值法辅助，很好地抑制了使用单一方法会出现的问题。从而到达较好的提取效果。参考文献：1Chen Chun-fu, Liang You-bi

17、n. Wall Rock Alterations of Endogenic Gold Deposit and the Significance to Gold E*plorationJ.Contributions to Geology and Mineral Resources Research,1999，5(3)：6370.淳福,梁有彬.生金矿床围岩蚀变及其找矿意义J地质找矿论丛，1999，5(3)：63702Carranza E J M, H ale M. Mineral Imaging with Landsat Thermatic Mapper Data for Hydrothermal

18、 Alteration Mapping in Heavily Vegetated TerraneJ.International Journal of Remote Sensing, 2002, 23: 22, 4827-4852.3 Crosta A P, Desouza F C R. Targeting Key Alteration Minerals in Epithermal Deposit s inPatagonia, Argentina, Using AS ER Imagery and Principal ponent Analysis J. International Journal

19、 of Remote Sensing, 2003, 24(21): 4233-4244. 4*u Han-qiu, Tang Fei, Analysis of new characteristics of the first Landsat 8 image and their eco-environmental significanceJ.ActaEcologicaSinica,2021,33 ( 11) : 32493257*涵秋,唐菲.新一代Landsat 系列卫星: Landsat8 遥感影像新增特征及其生态环境意义J.生态学报,2021,33 ( 11) : 324932575Li *

20、u-wen,NiuZhi-chun, Jiang Sheng. Study on the Usage of Landsat 8 Satellite Remote Sensing Image in Ecological Environment MonitoringJ.Environmental Monitoring and Forewarning20215(6):15.旭文,牛志春,晟Landsat8卫星OLI遥感影像在生态环境监测中的应用研究J.环境监控与预警20215(6):156 Moore F, Rastmanesh F, Asadi H, et a l. Mapping Mineral

21、ogical Alteration Using Principalponent Analysis and Matched Filter Processing in the T ak ab Area, Northwest Iran, from ASTER Data J . International Journal of Remote Sensing, 2021, 29(10): 2851-2867.7 Zhang Yu-jun Yang Jian-min The Method of Abstracting Remote Sensing Information of Alterated Rock

22、s in the Uncovered Bedrocks AreaJ Remote Sensing for Land and Resources1998, 36 (2) 46-53.玉君,建民.基岩裸露区蚀变岩遥感信息的提取方法J.国土资源遥感,1998,36(2)46-53.8 Zou Lin Yang Zi-an Zhu Gu-chang Wu De-wen *iang Ai-qin New Method of Alteration Information E*traction from Multi-spectrum Remote Sensing DateJ Geology and E*pl

23、oration,2006,42(6):71-76.邹林,自安,朱谷昌等.多光谱遥感蚀变信息提取新方法研究J.地质与勘探,2006, 42(6) : 71-76. 9Cui Yin-liang.Remote sensing geology and metallogenic prediction for lead-zinc deposit in the northeast of Yunnan ProvinceM. Beijing: Geological Publishing House,2021:80 银亮.滇东北铅锌银矿床遥感地质与成矿预测M.:地质,2021:8010Yu Min, Wen *

24、ing-ping, *u Jun-long, Chao Jiang-qin,Yang Yang, Wang Jun, Yi Bang-jin. Application of E*traction of Remote Sensing Alteration Anomalies based on Fractal Theory in Mao-ping Lead-zinc DepositJ Remote Sensing Technology and Application,2021(10):853-859.余敏,温兴平,*俊龙,晁江琴，炀，王军，易邦进.基于分形的遥感蚀变异常提取在毛坪铅锌矿中的应用J.

25、 遥感技术与应用,2021(10):853-859.11 Han Run-sheng,Chen Jin, Li Yuan,Ma De-yun, Zhao De-shun, Ma Geng-shengOre-controlling Tectonics and Prognosis of Concealed ores in Hui-zePb-Zn Deposit, YunnanJ,ActaMineralogicaSinica,2001,21(2):265-267.润生,进,元,马德云，德顺，马更生.会泽铅锌矿床构造控矿规律及其隐伏矿预测J.矿物学报,2001,21(2):265-267.12 Zha

26、ng Yu-jun Yang Jian-min Chen Wei A study of the Method for E*traction of Alteration Anomalies from the ETM+ (TM) Date and its Application:Geologic Basis and Spectral Precondition.J Remote Sensing for Land and Resources, 2002,4 (54 ): 30-37玉君,建民,薇. ETM (T M)蚀变遥感异常提取方法研究与应用地质依据和波谱前提. 国土资源遥感, 2002,4(54

27、):30 -3713Gan Fu-ping,Wang Run-sheng. Remote sensing rock information e*traction technology on the basis and methodsM. Beijing: Geological Publishing House, 2004.甘甫平,王润生.遥感岩矿信息提取根底与技术方法研究M.:地质.200414 Wei Yu-ming, Yang De-guang, Luo Chen.Can not be ignored in the prospecting limoniteJ.Jilin Geology,2

28、021,29(3):66 -68.玉明,得光,罗琛.在找矿中不可无视的褐铁矿化J.地质,2021,29(3):66 -68.15Chen Feng-yu, Gu Shang-yi. Limonite Characteristics and Its Significanse of Damaidi-*iaoJianshan Portion in Shaguochang Gold Deposit, Pan*ian County, Southwestern GuizhouProvinceJ. Journal of Guizhou University,2021,29(2):63-65. 凤雨,顾尚义.

29、盘县砂锅厂金矿大麦地小尖山矿段褐铁矿化特征及地质意义J.大学学报,2021,29(2):63-65.16Ma Jianwen Zhang QIdao The Method of Mapping Hydrothermally Altered Rocks under Various Environmental Conditions using TM DateJ.Contributions to Geology and Mineral Resources Research, 1994, 9 (2): 84-88马建文, 齐道.利用TM数据在多种环境因素干扰条件下填制蚀变岩的方法J .地质找矿论丛,

30、1994, 9 (2): 84-8817Shen Li-*ia Liu Li-ping Su *in-*u Zhu Yan Zhao Fang Study on Identification of Altered Wall-rock in Areas of Different Vegetation Coverages Based on Remote SensingJ.Geoscience, 2021,22(2): 293-298利霞,丽萍,新旭,祝艳,芳.不同植被覆盖率地区遥感矿化蚀变提取研究J. 现代地质, 2021,22(2): 293-29818 Zhang Yu-jun Zeng Ch

31、ao-ming Chen-wei The Method for E*traction of Alteration Anomalies from the ETM+ (TM) Date and its Application : Method Selection and Technological Flow ChartJ Remote Sensing for Land and Resources,2003 (2)：4449玉君,曾朝铭,薇.ETM+(TM)蚀变遥感异常提取方法研究与应用方法选择和技术流程J.国土资源遥感, 2003 (2)：444919 Zhang Nan-nan Zhou Ke-

32、fa Chen *i Li Hong A parative Study of E*traction Methods for Alteration Information Based on ETM+.J Remote Sensing for Land and Resources, 2021(2)楠楠,周可法,曦,宏.基于ETM+的遥感蚀变信息提取方法比照研究J国土资源遥感2021(2)20. Zhang Guo-rong Liu qing-shanFei Yi-qing Application of ETM+ Data to Alteration Remote Sensing Anomaly E

33、*traction in the Heicigou Area of Su-bei County, GansuJ Remote Sensing for Land and Resources,2021,34(3)：386390国荣,芦青山,费一清. ETM+数据在省肃北县黑刺沟一带蚀变遥感异常信息提取中的应用构造与成矿学,2021,34(3)：38639021 Zhao Yuan-hong Zhang Fu-*iang Chen Nan-feng The Application of Principal ponent Integration of Band Ratios to E*tracting

34、 Hydrothermal Alteration InformationJ Remote Sensing for Land and Resources, 1991, (3):12 -17.元洪,福祥,南峰.波段比值的主成分复合在热液蚀变信息提取中的应用J.国土资源遥感, 1991, (3):12 -17.22Duan Yuan-bin Liu Deng-zhong *u Tao *u Zheng-qiang Cui Zhi-qiang Zhao Pei-song The E*traction of Remotely Sensed Alteration Information from Dao-

35、cheng Area.J Remote Sensing for Land and Resources, ,2021, ,(4):92-95段元彬，登忠，*韬，*争强，志强，培松. 稻城地区遥感蚀变信息提取研究J. 国土资源遥感,2021, (4):92-95. z-E*traction of a bination of remote sensing image analysis Landsat8 OLI of iron stained alteration- To Zhao-tong hair Ping Pb-Zn CaseMa Wei Wen *ing-ping Shi Zhen-jie Shen Pan1.Faculty of Land Resource Engineering, Kunming University of Science and Technology, Kunming, 650093, China;2.Mineral Resources Prediction and Evaluation Engineering Laboratory of Yunnan Province, Kunming, 650093, ChinaAb