一种利用worldview-3遥感数据提取铁染信息的方法
技术领域
1.本发明属于遥感地质技术领域，具体涉及一种利用worldview-3遥感数据提取铁染信息的方法。
背景技术：
2.铁染是指岩石遭受中、低温热液作用或者氧化作用而发生赤铁矿化、褐铁矿化等蚀变，使岩石发生铁质浸染而变成红色。铁染可以为寻找某些矿产提供重要线索，比如中、低温热液成因的铜、铅、锌、金、铀等矿产，有铁帽标志的各种硫化物矿床，以及由氧化还原作用导致部分成矿元素富集的矿床等。因此，如果能知道某地区的铁染信息分布，对于寻找上述矿产具有重要的指导意义。遥感技术在提取矿物蚀变信息方面能发挥重要作用，已有学者利用landsattm、aster、hyperion、gf-1、worldview-2和worldview-3等遥感数据，采用主成分分析、波段比值运算、光谱角匹配、掩膜分析等方法提取了铁染等蚀变信息，但主要应用于地表岩石地层出露较好的热液型矿床，对于发育在盆地内部的沉积型矿床，目前应用的相对偏少，主要是因为盆地内多覆盖严重，利用遥感技术提取相关蚀变信息的难度较大。在热液型铀矿床中，上世纪末就被地质工作者发现在矿床附近多发育由赤铁矿化或针铁矿化引起的红化蚀变，不少学者利用遥感技术在热液型铀矿附近提取了铁染等相关蚀变信息，为下一步铀矿勘查提供了重要线索。在砂岩型铀矿附近，通常也发育因氧化作用引起的褐铁矿化和赤铁矿化，但多蚀变较弱，并且由于地表覆盖严重，只在局部出露。根据砂岩型铀成矿理论，这些铁染信息对寻找砂岩型铀矿具有重要的指示意义。利用遥感技术提取蚀变信息多应用于岩石地层出露较好的热液成因矿床区，在出露较差的盆地内应用较少，也缺乏有效的方法。因此，利用遥感技术提取盆地内或者盆山结合带地区的铁染信息的意义也十分重大。
技术实现要素：
3.本发明的目的在于：提供一种利用worldview-3遥感数据提取铁染信息的方法，该方法能够在盆山结合地带和盆地内部提取出地表弱的铁染信息，从而为寻找砂岩型铀矿或其它沉积型矿床提供线索。
4.实现本发明目的的技术方案：
5.一种利用worldviews-3遥感数据提取铁染信息的方法，该方法具体包括如下步骤：
6.步骤1、确定研究区或工作区范围，获取相应范围的wv-3数据；
7.步骤2、对获取的wv-3数据辐射校正；
8.步骤3、对辐射校正后的wv-3数据进行波段合成；
9.步骤4、在envi软件中打开波段合成后的文件；
10.步骤5、选择vnir-3，vnir-7，swir-1和swir-3四个波段进行主成分分析；
11.步骤6、选择步骤5分析完成后的文件，根据统计的特征向量矩阵选择主成分分析
波段；
12.步骤7、计算平均值和标准方差；
13.步骤8、根据计算的平均值和标准方差进行阈值分割和信息提取。
14.所述的步骤1中获取的wv-3数据为16个波段，包括8个可见-近红外波段和8个短波红外波段，一个全色波段。
15.所述可见-近红外波段的空间分辨率为1.2m，短波红外波段的空间分辨率为 3.7m，全色波段的空间分辨率0.3m。
16.所述的步骤2中的辐射校正包括辐射定标和flassh大气校正。
17.所述的步骤3中的波段合成，是将短波红外波段进行重采样，统一为可见光-近红外波段一致的空间分辨率。
18.所述步骤5中选择vnir-3、vnir-7、swir-1和vnir-7四个波段进行主成分分析的原因是vnir-3波段贡献系数与vnir-7波段的贡献系数符号相反， swir-1波段的贡献系数与swir-3波段的贡献系数符号也相反，vnir-7波段位置是三价铁的特征吸收峰位置。
19.所述的步骤6中获得特征向量矩阵后，首先看各波段的值正负是否与光谱吸收特征一致，光谱在相应波段吸收，其特征值为负，反射则为正；其次，看特征值的绝对值大小，绝对值越大，说明光谱信息较丰富，更能反映所提取的矿物信息；以上述两个特点来选择合适的波段，然后进行阈值分割。
20.所述的步骤8中阈值分割是根据步骤7中计算获得的平均值mean和标准方差stdev，以mean+2stdev，mean+2.5stdev和mean+3stdev为阈值分割界线，将提取的铁染信息分为三级、二级和一级，其中三级代表弱铁染信息，二级代表中等铁染信息，一级为强铁染信息。
21.本发明的有益技术效果在于：本发明提供的一种利用worldview-3遥感数据提取铁染信息的方法，选择wv-3(16波段)遥感数据的vnir-3、vnir-7、 swir-1和vnir-7四个波段，利用主成分分析方法，通过阈值分割和多级别信息提取，可以准确地获取研究区相关蚀变信息的分布和强弱。该方法不仅可以应用于地层岩石出露较好的裸露区，也可以应用于覆盖较为严重的盆山结合部位和盆地内部，其主要技术效果在于解决了盆地覆盖区难以提取一些重要地质弱信息的问题。本方法可以较为准确地提取出盆山结合带和盆地内部弱的铁染信息，为寻找盆地内沉积型矿床(如砂岩型铀矿)提供重要线索。
具体实施方式
22.下面结合实施例对本发明提供的一种利用worldview-3遥感数据提取铁染信息的方法作进一步详细说明。
23.本发明所提供的一种利用worldview-3遥感数据提取铁染信息的方法，该方法具体包括如下步骤：
24.步骤1、确定研究区或工作区范围，获取相应范围的wv-3数据
25.本次选择内蒙古乌拉特后旗本巴图地区作为试验区，获得该地区 worldview-3(wv-3)数据，面积约300km2。获取的wv-3数据为16波段，包括8个可见-近红外波段(空间分辨率1.2m)和8个短波红外波段(空间分辨率3.7m)，另外还有一个全色波段(空间分辨率0.3m)。
26.步骤2、辐射校正
27.对获取的wv-3数据进行辐射校正，主要包括辐射定标和flassh大气校正；
28.步骤3、波段合成
29.对辐射校正后的wv-3数据进行波段合成，将短波红外波段(swir)进行重采样，统一为可见光-近红外波段(vnir)一致的空间分辨率，即1.2m，然后另存为16个波段的文件f。
30.步骤4、打开波段合成文件
31.在envi软件中，打开步骤3另存的文件f，加载处理后的影像文件；
32.步骤5、主成分分析
33.选择加载的文件f，根据赤铁矿、褐铁矿等三价铁矿物的光谱特征，选择 vnir-3，vnir-7，swir-1和swir-3四个波段进行主成分分析，然后另存为文件p；
34.步骤6、主成分分析波段选择
35.选择步骤5完成的文件p，获得特征向量矩阵。
36.在获得特征向量矩阵后，首先看各波段的值正负是否与光谱吸收特征一致，光谱在相应波段吸收，其特征值为负，反射则为正。其次，看特征值的绝对值大小，绝对值越大，说明光谱信息较丰富，更能反映所提取的矿物信息。以上述两个特点来选择合适的波段，然后进行阈值分割，本次选择第四波段作为主成分进行分析，因为band-3和band-7符号相反，band-9和band-11符号相反，且绝对值相对较大。
37.步骤7、计算平均值和标准方差
38.加载步骤6所选的波段影像，然后选择快速统计，计算出主成分分析波段的最小值(min＝-1287.759766)、最大值(max＝1287.936401)、平均值(mean＝0) 和标准方差(stdev＝54.616399)。
39.步骤8、阈值分割和信息提取
40.阈值分割方法：根据步骤7中计算获得的平均值(mean)和标准方差(stdev)，以mean+2stdev＝109.2，mean+2.5stdev＝136.5和mean+3stdev＝163.8为阈值分割界线，将提取的铁染信息分为三级(阈值范围109.2～136.5)、二级(阈值范围 136.5～163.8)和一级(136.8～514.2875)，，其中三级代表弱铁染信息，二级代表中等铁染信息，一级为强铁染信息。
41.采用本发明的提取铁染信息的方法，解决了盆地覆盖区难以提取一些重要地质弱信息的问题，能够在盆地内部和盆山结合地区提取地层中和地表的弱铁染信息，不仅可以直观地获取铁染信息的分布，而且可以区分蚀变信息的强弱，从而帮助地质人员总结和分析地质规律，并为盆地内下一步矿产勘查提供方向。
42.上面结合实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。