Landsat8OLI遥感影像的组合分析法对铁染蚀变的提取研究.docx
- 文档编号:7527075
- 上传时间:2023-01-24
- 格式:DOCX
- 页数:9
- 大小:480.96KB
Landsat8OLI遥感影像的组合分析法对铁染蚀变的提取研究.docx
《Landsat8OLI遥感影像的组合分析法对铁染蚀变的提取研究.docx》由会员分享,可在线阅读,更多相关《Landsat8OLI遥感影像的组合分析法对铁染蚀变的提取研究.docx(9页珍藏版)》请在冰豆网上搜索。
Landsat8OLI遥感影像的组合分析法对铁染蚀变的提取研究
8遥感影像的组合分析法对铁染蚀变的提取研究
以云南昭通毛坪铅锌矿为例
马威,温兴平,石振杰,沈攀
(1.昆明理工大学国土资源工程学院,昆明650093;2.云南省矿产资源预测评价工程实验室,昆明650093)
摘要:
本文以云南昭通毛坪铅锌矿为研究区域,利用最新的8影像,在详细解读云南昭通毛坪的8影像的基础上,以蚀变围岩的光谱异常为理论依据,结合研究区的地质资料,根据矿物的光谱特征,利用掩膜技术去除水体、阴影等干扰信息,采用主成分分析为主、比值法为辅的方法利用卫星数据提取研究区内铁染异常蚀变信息,并用影像的主成分分析作为比对,研究表明,利用组合分析法,也就是在主成分分析的基础上,采用比值法辅助,对8遥感影像上的铁染异常蚀变信息有较好的提取效果。
关键词:
影像;蚀变信息;提取方法;毛坪铅锌矿
引言
围岩蚀变(),又称围岩交代蚀变或主岩交代蚀变,是指近矿围岩在流体的作用下发生的化学和物理变化,从而引起围岩化学成分和结构构造的变化[1]。
蚀变围岩因其矿物成分、岩石结构与正常围岩有所不同,所以二者在遥感图像表现为不同的色调,在波谱维表现为截然不同的波谱曲线[2-3]。
这为利用遥感影像提取围岩蚀变信息提供了基础。
到目前为止,围岩蚀变信息的提取多利用等影像数据。
2013年2月,成功发射了8卫星,8既保持了7卫星的基本特征,又在波段数量、波段光谱范围和影响的辐射分辨率方面有了改进。
2013年,徐涵秋基于8的首幅影像,分析和研究了8遥感影像的新增特征及其生态环境意义[4]。
李旭文利用8卫星上运营性陆地成像仪获取的太湖流域遥感影像数据,开展了其在大气能见度、蓝藻水华及建设用地等生态环境监测领域中的应用研究[5]。
数据在地质遥感找矿方面的应用还没有文献报道。
此外,基于多光谱数据进行矿化蚀变信息提取的方法主要有主成分分析法、比值方法、光谱角法等[6-7]。
这些方法在大量的实践应用中均表现出了不同的问题[8],主成分分析法依据蚀变围岩与非蚀变围岩波谱差异提取矿化蚀变异常时,使得提取微弱的蚀变围岩信息淹没于背景噪声中,比值方法在波段相除时减弱源数据的信噪比,光谱角法因多光谱数据谱带过宽、数目偏少而导致结果精度偏低。
针对以上问题,本研究在详细解读云南昭通毛坪的影像的基础上,以蚀变围岩的光谱异常为理论依据,提取遥感蚀变信息。
结合研究区的地质资料,根据矿物的光谱特征,采用主成分分析为主、比值法为辅的方法利用卫星数据提取研究区内蚀变信息,并利用7影像的蚀变信息提取结果进行对比分析。
1研究区概况
毛坪位于云南省东北部的昭通地区,地形上属云岭高原与四川盆地的交接处,典型的山地地貌,地势南高北低,地形切割强烈,地理坐标是N26.80°-28.50°,E102.88°-105.30°[9]。
矿区出露的地层有上泥盆统宰格组(D3),下石炭统大塘组(C1d)、下石炭统摆佐组(C1b)、中石炭统威宁组(C2w),下二叠统梁山组(P1l)、下二叠统栖霞茅口组(P1),及第四系(Q)(图1)。
矿体赋存于上泥盆统宰格组第二段(D32)和下石炭统摆佐组上段(C1b2)及中石炭统威宁组(C2w)[10]。
主要有容矿层的白云岩的退色、结晶和白云石化、黄铁矿化、重晶石化、方解石化、硅化[9]。
各种不同的围岩蚀变作用大致具有一定分布规律:
即黄铁矿大多出现在矿体以内或近矿强烈蚀变的中—粗晶白云岩中,显见它和矿化关系十分密切;铁白云石和铁方解石化常伴随矿体出现,一般矿体都与铁白云石和铁方解石化有关,离矿体越近白云石结晶越粗大,方解石团块越发育。
2数据来源
近年来,遥感技术和地理信息技术飞速发展,2013年,美国地质调查局(及太空署发射了“陆地卫星数据连续性任务”卫星(卫星,发射后更名为8)。
8卫星上携带有两个主要载荷:
(陆地成像仪)和(热红外传感器)。
为了避免大气吸收特征,还在波段数量、波段光谱范围和影响辐射分辨率方面有了改进。
此外,还新增了两个波段,蓝色波
(1)段可应用于海岸带观测,短波红外段(9)包括水汽强吸收特征,可用于云检测[5]。
成像仪包括9个短波谱段(波段1~波段9),幅宽185,其中全色波段地面分辨率为15m,其他谱段地面分辨率为30m。
如表1所示。
表1卫星数据波段
1
序号
波段
空间分辨率(m)
1
0.433–0.453
30
2
0.450–0.515
30
3
0.525–0.600
30
4
5
6
7
8
9
0.630–0.680
0.845-0.885
01.560-1.650
2.100-2.300
0.500-0.680
1.360-1.390
30
30
30
30
15
30
成像波段在波谱范围上做出了很大的调整。
成像仪的波段5(0.845μm~0.885μm)去除了波谱范围中波长为0.825μm附近的水蒸气吸收段;波段7波谱范围变窄,这将有助于减少提取找矿信息时碳酸盐岩对羟基的干扰;另外,为了使得有植被的地表和没有植被的地表之间的反差更大,成像仪的全色波段(波段8)的波谱范围也很窄[5]。
3蚀变提取的的理论依据
遥感数据的事变信息是依据电磁波与出露的岩石和地质体发生作用,产生岩石和地质体的特征光谱来提取的[11-12]。
研究区黄铁矿化大多出现在矿体以内或近矿强烈蚀变的中-粗晶白云岩中,显见它和矿化关系十分密切铁离子矿物的提取具有一定意义。
含铁(3+、2+)矿物,铁矿物以次生氧化物为主,部分作为热液蚀变带的原生矿物,如常见的褐铁矿、针铁矿()、赤铁矿、黄钾铁矾等含大量3+,也有少量2+的纤铁矿等氧化物[62]。
0.45µm,0.55µm,0.85µm,0.90µm,0.94µm存在较强的吸收谷。
一般来讲,浅部黄铁矿会由于氧化作用形成褐铁矿化,也就是说,研究区内的铁染蚀变在遥感影像上读取为褐铁矿化[13]。
褐铁矿化实际上并不是一种单独矿物,而是由水针铁矿、针铁矿、水纤铁矿、纤铁矿、更富含铝的氢氧化物、水的氢氧化铁胶凝体、以及含水的泥质、氧化硅等常共同产出而形成;含铁量达30%~40%[14]。
其主要成分的光谱特征含量如下图,在可见光及近红外区中,岩石特征谱带形成占有优势。
在在0.6-0.9μm间产生强的吸收谱带,因而诊断性吸收峰位于0.9μm处。
通过野外地质调查和岩性地层剖面测量,在典型岩性段内建立地质观察点,采集具有代表性的岩石,保证样品有风化面和新鲜面,用定位采样点位置并拍摄野外采样点照片,同时记录岩石特征、地质地貌特征以及环境参数;然后对采样点进行野外光谱测试,每个采样点测试三条光谱曲线,以备后期实验室数据处理减小误差,最后,对每个观察点的光谱曲线取平均值,所得的波谱曲线图如图所示,反映在影像上,2与5波段上具有明显的吸收谷,在4波段具有相对较高的反射率
图1主要铁染蚀变矿物及对应波段
(根据波谱库)
.1
图2野外铁染蚀变矿物及对应波段
.2
4蚀变信息提取
4.1研究区干扰信息的去除
由于研究区地物光谱、遥感成像时间以及辐射分辨率的限制,传感器很难精确地记录复杂的地表信息,因而误差就不可避免地存在于遥感数据获取的过程中,误差的存在降低了遥感数据的质量,.所以,在数据处理、信息提取前,需要对原始数据进行预处理[15]。
受研究区的地形地貌及数据时相等原因,要采用掩膜[15-16]的方法去除研究区范围内的水体、植被等干扰元素的影响,,使蚀变提取的结果更加精确。
掩膜前
掩膜后
图3掩膜处理前后对比图
.3
4.2蚀变提取及结果
主成分分析法[17-19]很好地去除了多波段间的相关性,使得变换后各主分量之间相互独立,分别代表不同的地质意义,并且互不重复,矿物的蚀变信息集中于某一个主分量中。
主成分分析法虽然可以降低波段间的相关性和冗余性,但是干扰地物,以及乘性噪声(部分大气衰减和地形差异等)的存在,影响异常信息提取,因此采用比值法[20]作为辅助手段,提取蚀变信息。
与铁离子有关的蚀变矿物主要是黄铁矿,由于野外的风化剥蚀,在遥感影像上,其表现为褐铁矿化。
特征光谱主要集中在25波段,在4波段表现为相对高反射。
根据蚀变矿物的波普曲线及研究区地质概况,确定蚀变提取方法为主成分分析法,在主成分分析中,为避免含羟基和含碳酸根矿物的干扰,舍弃了7波段。
根据研究区蚀变矿物波普特征,65比值扩大了6与5的光谱反差,具有增强铁化蚀变图像的亮度反差功能。
用65代替6,即应用2、4、5、43作为输入波段组合进行主成分分析,以增强铁染蚀变信息。
得到主成分分析的特征向量矩阵(表1),
表2(2、4、5、4/3)主成分特征向量表
2(2,4,5,4/3)
主成分
2
4
5
4/3
1
-0.139400
-0.369006
-0.918913
-0.000156
2
-0.488834
-0.781377
0.387933
0.000318
3
-0.861167
0.503273
-0.071459
0.001117
4
-0.000785
0.000868
-0.00006
0.999999
从表中可以看出,3上5特征向量为负,6/5特征向量为正,因此对3图像上进行主成分分析提取蚀变异常,并利用同研究区的影像的蚀变信息提取结果进行对比,其结果如下:
图4铁染蚀变异常分布图
4
图5铁染蚀变异常分布图
(取值波段为1、3、4、5)
5+
(1,3,4,5)
图中的红色区域为铁染蚀变异常与对比的情况,通过比对可以看出,采取主成分分析法为主、比值法为辅的方法,能较好地提取铁染蚀变;采用主成分分析法,能较好地提取羟基蚀变。
5结论
(1)通过分析地质体波谱特征和遥感蚀变信息提取的地质依据,依据矿区矿化蚀变类型发育特征和各类典型蚀变矿物的波谱特征,提取研究区的铁染蚀变信息。
综合比较各蚀变信息提取方法,结合研究区实际情况采用主成分分析法为主,比值法为辅进行遥感蚀变信息提取。
其中,采用2、4、5、43波段组合进行主成分分析获取铁染蚀变信息,采用2,5,67波段组合进行羟基蚀变信息提取。
(2)主成分分析法、比值分析法是现今遥感蚀变信息提取中广泛应用的方法。
主成分分析法,需要将干扰地物去除才能取得较好的效果。
而根据研究区实际情况,对铁染蚀变信息而言,干扰地物的影响很难消除,因此在主成分分析的基础上,采用比值法辅助,很好地克服了使用单一方法会出现的问题。
从而达到较好的提取效果。
参考文献:
[1],.[J].,1999,5(3):
63—70.[陈淳福,梁有彬.内生金矿床围岩蚀变及其找矿意义[J]地质找矿论丛,1999,5(3):
63—70].
[2]EJM,HM.[J],2002,23:
22,4827-4852.
[3]AP,FCR.s
,[J].,2003,24(21):
4233-4244.
[4],,8[J],2013,33(11):
3249-3257.[徐涵秋,唐菲.新一代系列卫星:
8遥感影像新增特征及其生态环境意义[J].生态学报,2013,33(11):
3249-3257].
[5],,.8[J]20135(6):
1-5.[李旭文,牛志春,姜晟.8卫星遥感影像在生态环境监测中的应用研究[J].环境监控与预警20135(6):
1-5]
[6]F,F,H,al.T,,[J].,2008,29(10):
2851-2867.
[7][J]1998,36
(2)46-53.[张玉君,杨建民.基岩裸露区蚀变岩遥感信息的提取方法[J].国土资源遥感,1998,36
(2)46-53].
[8][J],2006,42(6):
71-76.[邹林,杨自安,朱谷昌等.多光谱遥感蚀变信息提取新方法研究[J].地质与勘探,2006,42(6):
71-76].
[9][M].:
2011:
80[崔银亮.滇东北铅锌银矿床遥感地质与成矿预测[M].北京:
地质出版社,2011:
80]
[10],,,,,,.[J],2014(10):
853-859.[余敏,温兴平,徐俊龙,晁江琴,杨炀,王军,易邦进.基于分形的遥感蚀变异常提取在毛坪铅锌矿中的应用[J].遥感技术与应用,2014(10):
853-859].
[11],,,,,,[J],2001,21
(2):
265-267.[韩润生,陈进,李元,马德云,赵德顺,马更生.云南会泽铅锌矿床构造控矿规律及其隐伏矿预测[J].矿物学报,2001,21
(2):
265-267].
[12]A().[J],2002,4(54):
30-37[张玉君,杨建民,陈薇.(TM)蚀变遥感异常提取方法研究与应用——地质依据和波谱前提.国土资源遥感,2002,4(54):
30-37]
[13],.[M].:
2004.[甘甫平,王润生.遥感岩矿信息提取基础与技术方法研究[M].北京:
地质出版社.2004].
[14],,.[J],2010,29(3):
66-68.[魏玉明,杨得光,罗琛.在找矿中不可忽视的褐铁矿化[J].吉林地质,2010,29(3):
66-68].
[15],.,,[J].,2012,29
(2):
63-65.[陈凤雨,顾尚义.盘县砂锅厂金矿大麦地——小尖山矿段褐铁矿化特征及地质意义[J].贵州大学学报,2012,29
(2):
63-65].
[16][J].,1994,9
(2):
84-88[马建文,张齐道.利用数据在多种环境因素干扰条件下填制蚀变岩的方法[J].地质找矿论丛,1994,9
(2):
84-88]
[17][J],2008,22
(2):
293-298[沈利霞,刘丽萍,苏新旭,祝艳,赵芳.不同植被覆盖率地区遥感矿化蚀变提取研究[J].现代地质,2008,22
(2):
293-298]
[18]():
[J],2003
(2):
44—49[张玉君,曾朝铭,陈薇()蚀变遥感异常提取方法研究与应用——方法选择和技术流程[J].国土资源遥感,2003
(2):
44—49].
[19]A.[J],2012
(2)[张楠楠,周可法,陈曦,李宏.基于的遥感蚀变信息提取方法对比研究[J]国土资源遥感2012
(2)]
[20].,[J],2010,34(3):
386—390.[张国荣,芦青山,费一清.数据在甘肃省肃北县黑刺沟一带蚀变遥感异常信息提取中的应用.大地构造与成矿学,2010,34(3):
386—390]
[21][J],1991,(3):
12-17.[赵元洪,张福祥,陈南峰.波段比值的主成分复合在热液蚀变信息提取中的应用[J].国土资源遥感,1991,(3):
12-17].
[22].[J],,2008,,(4):
92-95[段元彬,刘登忠,徐韬,徐争强,崔志强,赵培松.稻城地区遥感蚀变信息提取研究[J].国土资源遥感,2008,(4):
92-95]
a8
(1,,,650093,;
2,,650093,)
:
8,8,,,,,,,,+a,a,,,8.
:
;;;
作者简介:
马威(1989-),女,内蒙古赤峰人,博士研究生,研究方向:
遥感地质.
通信作者:
温兴平(1970-),男,山西兴县人,教授,博士生导师,研究方向:
遥感地质.
基金项目:
国家自然科学基金(41101343);云南省高技术产业发展项目计划
联系方式:
手机:
邮箱:
378486149
单位:
昆明理工大学,国土资源工程学院;
云南省矿产资源预测评价工程实验室
地址:
云南省昆明市一二一大街文昌路68号昆明理工大学莲华校区学生公寓1C309
邮编:
650093
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Landsat8OLI 遥感 影像 组合分析 铁染蚀变 提取 研究