多元数据分析与遥感矿化蚀变信息提取模型.docxVIP

多元数据分析与遥感矿化蚀变信息提取模型 第1期，总第67期 2006年3月15日 国土资源遥感 REMOTESENSINGFORLANDamp;RESOURCES No. 1,2006 M81H.,2006 多元数据分析与遥感矿化蚀变信息提取模型 吴德文?一,朱谷昌，张远飞,袁继明 （1 ?中国地质大学，北京100083;2.有色金属矿产地质调查中心，北京100814;3 .桂林 矿产地 质研究院，桂林541004） 摘要:蚀变矿物的光谱行为造成了蚀变岩的光谱特征?利用典型岩石的多元数据分 析，可以对蚀变矿物成份和光 谱数据Z间的关系作出定性和定量评价，进而建立遥感矿化蚀变信息提取模型?笔 者利用青海某金多金属矿化区 典型岩（矿）石样品的野外实测波谱数据和化学分析数据，进行了基于多元数据分 析的遥感矿化蚀变信息提取模型 研究，建立的比值组合线性回归模型较之单一比值方法具有更好的应用效果. 关键词:多元数据分析;遥感信息提取模型 中图分类号:TP79文献标识码?A文章编号:1001 —070X（2006）01 —0022 一 o4 0引言 图像比值组合是一种十分有用的蚀变信息增强 处理方法，但由于蚀变岩受其自身反射和辐射强度 以及环境条件影响，常规増强蚀变信息的比值因子 （如TM3/TM1增强铁化信息,TM5/TM7增强泥化信 息）并非在所有地区都有效.换言之,对于不同的地 区.由于蚀变岩石及”背景”地物的光谱特征存在差 异，在提取矿化蚀变信息时,需要根据岩石（矿物）的 光谱特征分析来调整比值因子或波段组合;另外，不 同地区的蚀变岩石与”背景”地物在光谱空间的”聚 类结构”特征也不同，所以，矿化蚀变信息提取的方 法及遥感信息模型的建立也不存在固定的模式.事 实上，利用对岩石光谱数据进行多元数据分析，可以 研究和揭示不同地区的蚀变岩的光谱个性特征，进 而建立遥感矿化蚀变信息提取模型. 1多元数据获取 选择青海茫崖柴水沟一采石沟金多金属矿化区 作为研究区,采用的多元数据包括37个典型岩（矿） 石样品的野外实测波谱数据，微量Au分析和岩石组 分分析结果. 为了同TM图像数据进行类比,对每个样品的 波谱数据在光谱范围450?520nm,520?600nm, 630?690nm,760?900nm,1550?1750nm 以及 2080?2350nm 中分別获取与 TM1,TM2,TM3,TM4, TM5和TM7波段图像对应的光谱值,而后求取比值 3/TM1 ,TM3/TM2,TM3/TM4,TM5/TM7,TM5/ TM3,TM5/TM4,TM7/TM3 及 TM7/TM4.把岩矿石 的组份分析结果（包括SiO:,Al:0?CaO,Mgo, Fe203,FeO,0,Na20,Au,0+Na20+A1203,CaO +Mgo及Fe:0,/FeO等）和微量金分析结果也参与 到多元数据统计分析中来. 2多元数据的数理统计分析 对所有数值变量进行相关矩阵分析和R一型聚 类分析（图1）,可以得出如下结果： R 一型聚类分成两个大类时,变量分成两大 组,一组为 CaO,CaO+Mgo,Mgo,FeO,Na20,A1203, 0+Na20+A1203,Si02,0 和 TM3/TM4,另一组 为 TM5/TM3,TM5/TM4,TM7/TM4,TM2,TM3,TM4, TM1 ,TM5,TM7,TM3/TM 1 ,Fe203,Fe203/FeO,TM5/ TM7和A”?这说明:①与泥化，碳酸盐化和硅化相 关的变量同与铁化有关的变量，刚好与两大变量聚 类相对应;②光谱测试变量只有比值TM3/TM4与 泥化，碳酸盐化和硅化蚀变有关，其它波段和比值均 与铁化相关，这反映了样品的光谱变化主要来铁化 蚀变;③本地区的泥化，碳酸盐化和硅化3种蚀变相 互关系密切;④铁化蚀变与A”有较好的相关性. 收稿日期:2005—05—23;修订日期:2005—06—27 基金项目:”十?五”国家科技攻关计划项目”西部高寒山区铜金矿产勘查评价技术 与综合示范研究”(编号-2001BA609A 一 04). 第1期吴德文，等:多元数据分析与遥感矿化蚀变信息提取模型?23? 变量序号+0…?? 5 !D????.j: . C 口 015 CaO+MgO23 MgO16 Fe018 NaO20 A12OJ14 a.f22 f0j!3 019 TM3,TM4,0 TM5,TM3? TM5,TM4,2 TM7,TM49 Af22 3 44 TM1 J5 76 TM3,TMU Fe20J17 FejOtFeO24 TMs 九 M78 A21 图124个变量的R 一型聚类树状图 (K—Na-Al 表示 K20+Na20+J4|c203) 当R 一型聚类细分成5类时，则有[CaO, Ca