ETM+和ASTER数据在岩体提取中的对比研究.docVIP

ETM+和ASTER数据在岩体提取中的对比研究 　　【摘 要】遥感岩体信息的提取是地质调查和找矿的一个重要技术手段。本文选择位于西昆仑西段的某铜矿区作为研究区，选用ASTER和ETM+影像为数据源，利用这两种影像来分别提取研究区内的花岗岩体，对两者提取的结果进行对比分析。结果表明，相较于ETM+数据，ASTER数据在岩体信息提取方面具优势性，岩体的成功提取，为区域地质调查和地质找矿都能提供良好的科学依据。 　　【关键词】ASTER;ETM+;主成分分析;提取;岩体 　　1 引言 　　成矿过程中，岩体和围岩之间的蚀变带往往是矿物富集的地带。因此，在地质找矿过程中，圈定岩体的界线对成功找矿将会奠定一个良好基础。现今找矿主要集中在我国西部条件艰苦的无人区，若能借助遥感作为辅助手段进行找矿，在节省人力物力的同时，也使找矿范围更具针对性。本文分别利用ASTER和ETM+两种影像对西昆仑西段的某铜矿区的花岗岩体进行提取，旨在对比分析两种影像对该岩体的识别能力，圈定岩体与围岩之间的蚀变带，为外业找矿明确目标和范围。 　　2 ASTER和ETM+数据特征分析及预处理 　　2.1 ETM+和ASTER数据特征 　　ETM+是NASA于1999年7月15日成功发射的美国陆地卫星Landsat携带的对地观测传感器，是一台8波段的多光谱扫描辐射计，工作于可见光、近红外、短波红外和热红外波段[1]。 　　ASTER是搭载于Terra卫星上的多光谱成像仪，它涵盖了可见光到近红外14个波段，包括3个15m空间分辨的可见光和近红外波段、6个30m空间分辨率的短波红外波段以及5个90m空间分辨率的热红外波段。 　　2.2 数据预处理 　　在提取岩体之前先进行影像预处理，预处理包括几何校正、裁剪、辐射定标、大气校正和掩膜等处理。预处理最为关键的环节在于大气校正，大气校正（本文采用FLAASH大气校正模型）是为消除大气和光照等因素对地物反射的影响，获得地物真实反射率和辐射率[2、3]。 　　3 岩体的提取及对比分析 　　多光谱遥感影像提取地物信息的方法主要有波段比值法、主成分分析法、图像掩模、图像分割、光谱角填图等几种方法。综合多方面因素考虑及其多次试验结果分析，主成分分析法较为理想。 　　3.1 ASTER数据的岩体提取 　　利用ENVI软件的统计分析工具分别对ASTER影像进行统计分析，选用信息量丰富且相互间数据冗余量小的4个波段进行主成分分析。经分析，ASTER的1、4、6和7四个波段满足条件，选为主成分分析的组合波段。 　　利用CROSTA主成分分析法通过用ASTER 1、4、6和7进行主成分分析，得到其主成分变化特征向量矩阵如表1所示。 　　通过分析表1中的特征向量矩阵，PC1和PC2虽然包含80%以上的地物信息，但是在这些信息中，几乎每种地物的信息都蕴含其中，没具有特殊性;PC3和PC4包含的信息虽少，但它们往往蕴含较特殊的地物信息，因此对岩体矿物的识别有特殊的意义。表中分析得知，主成分PC3特征向量绝对值大的是b4和b7，其绝对值分别为0.725760和0.625661，并且两个波段的特征向量符号相反，加上b1和b6载荷因子较小，说明PC3的信息量大部分来自b4和b7这两个波段。特征向量符号为负的b7表示减信息，为正的表示加信息，经过求反PC3，可以增强对岩体的信息。因此，可用PC3作为提取岩体的信息增强图像。 　　得到可作为提取岩体的信息增强图像PC3后，经过3*3的低值滤波处理后，再经过密度分割处理，最后得到提取后的图像如图2所示（绿色表示岩体）。 　　表1 ASTER 1467波段主成分变化特征向量矩阵 　　主成分 　　Band 1 　　Band 4 　　Band 5 　　Band 6 　　PC1 　　-0.328850 　　-0.657373 　　-0.677393 　　-0.029285 　　PC2 　　-0.728943 　　-0.279768 　　0.624344 　　0.023814 　　PC3 　　-0.600061 　　0.699408 　　-0.385386 　　-0.047256 　　PC4 　　-0.020665 　　0.020500 　　-0.053014 　　0.998169 　　表2 ETM+ 1457波段主成分变化特征向量矩阵 　　主成分 　　Band 1 　　Band 4 　　Band 5 　　Band 6 　　PC1 　　-0.328850 　　-0.657373 　　-0.677393 　　-0.029285 　　PC2 　　-0.728943 　　-0.279768 　　0.624