一种利用WorldView-3遥感数据提取铁染信息的方法技术

本发明专利技术属于遥感地质技术领域，具体涉及一种利用WorldView

【技术实现步骤摘要】
一种利用WorldView
‑
3遥感数据提取铁染信息的方法


[0001]本专利技术属于遥感地质
，具体涉及一种利用WorldView
‑
3遥感数据提取铁染信息的方法。

技术介绍

[0002]铁染是指岩石遭受中、低温热液作用或者氧化作用而发生赤铁矿化、褐铁矿化等蚀变，使岩石发生铁质浸染而变成红色。铁染可以为寻找某些矿产提供重要线索，比如中、低温热液成因的铜、铅、锌、金、铀等矿产，有铁帽标志的各种硫化物矿床，以及由氧化还原作用导致部分成矿元素富集的矿床等。因此，如果能知道某地区的铁染信息分布，对于寻找上述矿产具有重要的指导意义。遥感技术在提取矿物蚀变信息方面能发挥重要作用，已有学者利用LandsatTM、ASTER、Hyperion、GF
‑
1、WorldView
‑
2和WorldView
‑
3等遥感数据，采用主成分分析、波段比值运算、光谱角匹配、掩膜分析等方法提取了铁染等蚀变信息，但主要应用于地表岩石地层出露较好的热液型矿床，对于发育在盆地内部的沉积型矿床，目前应用的相对偏少，主要是因为盆地内多覆盖严重，利用遥感技术提取相关蚀变信息的难度较大。在热液型铀矿床中，上世纪末就被地质工作者发现在矿床附近多发育由赤铁矿化或针铁矿化引起的红化蚀变，不少学者利用遥感技术在热液型铀矿附近提取了铁染等相关蚀变信息，为下一步铀矿勘查提供了重要线索。在砂岩型铀矿附近，通常也发育因氧化作用引起的褐铁矿化和赤铁矿化，但多蚀变较弱，并且由于地表覆盖严重，只在局部出露。根据砂岩型铀成矿理论，这些铁染信息对寻找砂岩型铀矿具有重要的指示意义。利用遥感技术提取蚀变信息多应用于岩石地层出露较好的热液成因矿床区，在出露较差的盆地内应用较少，也缺乏有效的方法。因此，利用遥感技术提取盆地内或者盆山结合带地区的铁染信息的意义也十分重大。

技术实现思路

[0003]本专利技术的目的在于：提供一种利用WorldView
‑
3遥感数据提取铁染信息的方法，该方法能够在盆山结合地带和盆地内部提取出地表弱的铁染信息，从而为寻找砂岩型铀矿或其它沉积型矿床提供线索。
[0004]实现本专利技术目的的技术方案：
[0005]一种利用WorldViews
‑
3遥感数据提取铁染信息的方法，该方法具体包括如下步骤：
[0006]步骤1、确定研究区或工作区范围，获取相应范围的WV
‑
3数据；
[0007]步骤2、对获取的WV
‑
3数据辐射校正；
[0008]步骤3、对辐射校正后的WV
‑
3数据进行波段合成；
[0009]步骤4、在ENVI软件中打开波段合成后的文件；
[0010]步骤5、选择VNIR
‑
3，VNIR
‑
7，SWIR
‑
1和SWIR
‑
3四个波段进行主成分分析；
[0011]步骤6、选择步骤5分析完成后的文件，根据统计的特征向量矩阵选择主成分分析
波段；
[0012]步骤7、计算平均值和标准方差；
[0013]步骤8、根据计算的平均值和标准方差进行阈值分割和信息提取。
[0014]所述的步骤1中获取的WV
‑
3数据为16个波段，包括8个可见
‑
近红外波段和8个短波红外波段，一个全色波段。
[0015]所述可见
‑
近红外波段的空间分辨率为1.2m，短波红外波段的空间分辨率为 3.7m，全色波段的空间分辨率0.3m。
[0016]所述的步骤2中的辐射校正包括辐射定标和Flassh大气校正。
[0017]所述的步骤3中的波段合成，是将短波红外波段进行重采样，统一为可见光
‑
近红外波段一致的空间分辨率。
[0018]所述步骤5中选择VNIR
‑
3、VNIR
‑
7、SWIR
‑
1和VNIR
‑
7四个波段进行主成分分析的原因是VNIR
‑
3波段贡献系数与VNIR
‑
7波段的贡献系数符号相反， SWIR
‑
1波段的贡献系数与SWIR
‑
3波段的贡献系数符号也相反，VNIR
‑
7波段位置是三价铁的特征吸收峰位置。
[0019]所述的步骤6中获得特征向量矩阵后，首先看各波段的值正负是否与光谱吸收特征一致，光谱在相应波段吸收，其特征值为负，反射则为正；其次，看特征值的绝对值大小，绝对值越大，说明光谱信息较丰富，更能反映所提取的矿物信息；以上述两个特点来选择合适的波段，然后进行阈值分割。
[0020]所述的步骤8中阈值分割是根据步骤7中计算获得的平均值Mean和标准方差Stdev，以Mean+2Stdev，Mean+2.5Stdev和Mean+3Stdev为阈值分割界线，将提取的铁染信息分为三级、二级和一级，其中三级代表弱铁染信息，二级代表中等铁染信息，一级为强铁染信息。
[0021]本专利技术的有益技术效果在于：本专利技术提供的一种利用WorldView
‑
3遥感数据提取铁染信息的方法，选择WV
‑
3(16波段)遥感数据的VNIR
‑
3、VNIR
‑
7、 SWIR
‑
1和VNIR
‑
7四个波段，利用主成分分析方法，通过阈值分割和多级别信息提取，可以准确地获取研究区相关蚀变信息的分布和强弱。该方法不仅可以应用于地层岩石出露较好的裸露区，也可以应用于覆盖较为严重的盆山结合部位和盆地内部，其主要技术效果在于解决了盆地覆盖区难以提取一些重要地质弱信息的问题。本方法可以较为准确地提取出盆山结合带和盆地内部弱的铁染信息，为寻找盆地内沉积型矿床(如砂岩型铀矿)提供重要线索。
具体实施方式
[0022]下面结合实施例对本专利技术提供的一种利用WorldView
‑
3遥感数据提取铁染信息的方法作进一步详细说明。
[0023]本专利技术所提供的一种利用WorldView
‑
3遥感数据提取铁染信息的方法，该方法具体包括如下步骤：
[0024]步骤1、确定研究区或工作区范围，获取相应范围的WV
‑
3数据
[0025]本次选择内蒙古乌拉特后旗本巴图地区作为试验区，获得该地区 WorldView
‑
3(WV
‑
3)数据，面积约300km2。获取的WV
‑
3数据为16波段，包括8个可见
‑
近红外波段(空间分辨率1.2m)和8个短波红外波段(空间分辨率3.7m)，另外还有一个全色波段(空间分辨率0.3m)。
[0026]步骤2、辐射校正
[0027]对本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种利用WorldViews
‑
3遥感数据提取铁染信息的方法，其特征在于，该方法具体包括如下步骤：步骤(1)、确定研究区或工作区范围，获取相应范围的WV
‑
3数据；步骤(2)、对获取的WV
‑
3数据辐射校正；步骤(3)、对辐射校正后的WV
‑
3数据进行波段合成；步骤(4)、在ENVI软件中打开波段合成后的文件；步骤(5)、选择VNIR
‑
3，VNIR
‑
7，SWIR
‑
1和SWIR
‑
3四个波段进行主成分分析；步骤(6)、选择步骤(5)分析完成后的文件，根据统计的特征向量矩阵选择主成分分析波段；步骤(7)、计算平均值和标准方差；步骤(8)、根据计算的平均值和标准方差进行阈值分割和信息提取。2.根据权利要求1所述的一种利用WorldView
‑
3遥感数据提取铁染信息的方法，其特征在于：所述的步骤(1)中获取的WV
‑
3数据为16个波段，包括8个可见
‑
近红外波段和8个短波红外波段，一个全色波段。3.根据权利要求2所述的一种利用WorldView
‑
3遥感数据提取铁染信息的方法，其特征在于：所述可见
‑
近红外波段的空间分辨率为1.2m，短波红外波段的空间分辨率为3.7m，全色波段的空间分辨率0.3m。4.根据权利要求2所述的一种利用WorldView
‑
3遥感数据提取铁染信息的方法，其特征在于：所述的步骤(2)中的辐射校正包括辐射定标和Flassh大气校正。5.根据权利要求4所述的一种利用WorldView
‑
3遥...

【专利技术属性】
技术研发人员：童勤龙，叶发旺，易敏，伍炜超，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：