CN108333129A - 一种蚀变信息获取方法及系统 - Google Patents

一种蚀变信息获取方法及系统 Download PDF

Info

Publication number
CN108333129A
CN108333129A CN201810145611.3A CN201810145611A CN108333129A CN 108333129 A CN108333129 A CN 108333129A CN 201810145611 A CN201810145611 A CN 201810145611A CN 108333129 A CN108333129 A CN 108333129A
Authority
CN
China
Prior art keywords
altered mineral
principal component
formula
mineral
band ratio
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810145611.3A
Other languages
English (en)
Other versions
CN108333129B (zh
Inventor
荆林海
丁海峰
唐韵玮
李慧
戈文艳
苗中杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Remote Sensing and Digital Earth of CAS
Original Assignee
Institute of Remote Sensing and Digital Earth of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Remote Sensing and Digital Earth of CAS filed Critical Institute of Remote Sensing and Digital Earth of CAS
Priority to CN201810145611.3A priority Critical patent/CN108333129B/zh
Publication of CN108333129A publication Critical patent/CN108333129A/zh
Application granted granted Critical
Publication of CN108333129B publication Critical patent/CN108333129B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/55—Specular reflectivity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793—Remote sensing

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Processing (AREA)

Abstract

本申请公开了一种蚀变信息获取方法及系统，其中，所述蚀变信息获取方法通过对传统的主成分分析法进行改进，根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别矿物信息的过程中分为三种情况，进行一次或多次主成分分析，以实现直接利用主成分分析法获得待识别蚀变信息的目的，避免了在深切割地形地区从遥感图像中提取蚀变信息可能产生的异常，以及在地表覆盖严重地区的遥感图像中难以有效提取蚀变矿物信息的问题。

Description

一种蚀变信息获取方法及系统
技术领域
本申请涉及遥感图像处理技术领域，更具体地说，涉及一种蚀变信息获取方法及系统。
背景技术
蚀变作用是指当流体或热液作用于岩石时，两者之间必然趋于达到化学平衡和发生物质交换作用，其结果导致岩石中部分或全部原矿物的消失和新矿物的形成，即原岩中某些物质部分或全部地带出和新物质带入的交代作用。蚀变信息对于矿物探测具有重大的意义，因此，通过提取遥感地图中的蚀变信息是矿物探测过程中的重要一环。
从遥感图像中提取蚀变信息的现有方法主要包括波段比值法和主成分分析法，例如Crosta技术。其中，在使用波段比值法提取深切割地区的遥感图像的蚀变信息时，可能会在波段比值过程中产生异常结果。这是因为深切割地区的海拔落差较大，山体坡度较陡，遥感图像中含有明显的地形阴影，如果遥感图像预处理过程中的大气校正不甚精确，遥感图像中阴影的某些波段的值有可能接近于0，从而在波段比值过程中产生异常高的结果，从而导致蚀变信息的提取产生较大的误差。在使用Crosta技术提取覆盖地区的遥感图像的蚀变矿物时，也难以达到良好的提取效果。这是因为在这些地区，矿物表面覆盖了大量的杂质(土、沙、砂石、植被等)。这些杂质遮盖了蚀变矿物的信息，使得蚀变矿物的信息在遥感图像中显示较弱，从而难以在主成分分析过程中将蚀变信息集中到单一的主成分中，进而导致蚀变矿物信息提取的效果较差。
发明内容
为解决上述技术问题，本申请提供了一种蚀变信息获取方法及系统，以解决在提取深切割地区的遥感图像的蚀变信息提取时可能产生的异常结果，以及在提取覆盖严重的地区的遥感图像的蚀变信息时，提取效果较差的问题。
为解决上述技术问题，本申请实施例提供了如下技术方案：
一种蚀变信息获取方法，用于遥感图像的蚀变矿物的信息提取，所述蚀变信息获取方法包括：
获取遥感图像，根据所述遥感图像确定待使用卫星数据的种类；
对确定的待使用卫星数据进行预处理；
根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式；
根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息；
所述根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息包括：
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行主成分分析，主成分分析结果中的第二主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bm进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bc进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的变信息，其中，n、m、i、c均为大于或等于1的正整数。
可选的，所述待使用卫星数据的种类为ASTER卫星数据或Landsat 8OLI 卫星数据或哨兵2号卫星数据或高分1号卫星数据或WorldView-3卫星数据。
可选的，所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或Fe3+铁染蚀变矿物或Fe2+铁染蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
可选的，所述根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式包括：
当所述待使用卫星数据的种类为ASTER卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为铁帽的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为Landsat 8OLI卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为WorldView-3卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe2+铁染蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
一种蚀变信息获取系统，用于遥感图像的蚀变矿物的信息提取，所述蚀变信息获取系统包括：
图像获取模块，用于获取遥感图像，根据所述遥感图像确定待使用卫星数据的种类；
预处理模块，用于对确定的待使用卫星数据进行预处理；
公式确定模块，用于根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式；
蚀变提取模块，用于根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息；
所述蚀变提取模块根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息具体用于，
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行主成分分析，主成分分析结果中的第二主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bm进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bc进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息，其中，n、m、i、c均为大于或等于1的正整数。
可选的，所述待使用卫星数据的种类为ASTER卫星数据或Landsat 8OLI 卫星数据或哨兵2号卫星数据或高分1号卫星数据或WorldView-3卫星数据。
可选的，所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或Fe3+铁染蚀变矿物或Fe2+铁染蚀变矿物或铁帽蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
可选的，所述公式确定模块具体用于，当所述待使用卫星数据的种类为 ASTER卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为 Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为铁帽蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为Landsat 8OLI卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为WorldView-3卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe2+铁染蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
从上述技术方案可以看出，本申请实施例提供了一种蚀变信息获取方法及系统，其中，所述蚀变信息获取方法通过对传统的主成分分析法进行改进，在对确定的待使用卫星数据进行预处理后，根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式，最后根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的过程中分为三种情况，进行一次或多次主成分分析，以实现直接利用主成分分析法获得待识别蚀变矿物的目的，避免了在深切割地区的遥感图像中进行波段除法运算，从而避免了部分阴影区会得到极大增强的可能，进而避免了对深切割地区的遥感图像进行蚀变信息提取时出现的高异常情况；而当遥感图像所对应的地区为覆盖严重的地区时，通过三种不同情况下的一次或多次主成分分析，可以将待识别蚀变矿物的信息集中到最终获得的第一主成分或第二主成分中，从而实现良好的蚀变信息的提取效果。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本申请的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其它的附图。
图1为本申请的一个实施例提供的一种蚀变信息获取方法的流程示意图；
图2为本申请的一个实施例提供的根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息的具体流程示意图；
图3为本申请的一个实施例提供的一种遥感图像的光谱波段设置情况的示意图；
图4为本申请的另一个实施例提供的一种遥感图像的光谱波段设置情况的示意图；
图5为本申请的又一个实施例提供的一种遥感图像的光谱波段设置情况的示意图；
图6(a)为龙巴其保研究区的哨兵2号图像的子区图像；
图6(b)为利用现有技术中的主成分分析法对图6(a)进行蚀变信息提取的结果图像；
图6(c)为利用现有技术中的波段比值法对图6(a)进行蚀变信息提取的结果图像；
图6(d)为利用本申请实施例提供的蚀变信息获取方法对图6(a)进行蚀变信息提取的结果图像；
图7(a)为龙巴其保研究区的哨兵2号图像的另一个子区图像；
图7(b)为利用现有技术中的主成分分析法对图7(a)进行蚀变信息提取的结果图像；
图7(c)为利用现有技术中的波段比值法对图7(a)进行蚀变信息提取的结果图像；
图7(d)为利用本申请实施例提供的蚀变信息获取方法对图7(a)进行蚀变信息提取的结果图像。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整的描述。显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
本申请实施例提供了一种蚀变信息获取方法，如图1所示，用于遥感图像的蚀变矿物的信息提取，所述蚀变信息获取方法包括：
S101：获取遥感图像，根据所述遥感图像确定待使用卫星数据的种类；
S102：对确定的待使用卫星数据进行预处理；
在确定了待使用卫星数据后的预处理过程以ASTER卫星数据为例，主要包括：辐射定标→大气校正→影像边框和串扰效应去除→几何纠正→正射纠正→去除云雪水体等干扰因素；这些步骤均在ENVI5.2和ArcGIS10.2中完成；每个步骤的具体过程已为本领域技术人员所熟知，本申请在此不做赘述。
S103：根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式；
S104：根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息；
参考图2，所述根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息包括：
S1041：当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行主成分分析，主成分分析结果中的第二主成分为所述待识别蚀变矿物的信息；
在S1041中的遥感图像的光谱排布情况参考图3，在图3中示出了光谱反射峰B1和光谱吸收峰B2。
S1042：当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰 Bm进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息；
在S1042中的遥感图像的光谱排布情况参考图4，在图4中示出了光谱反射峰B1、B3和光谱吸收峰B2。
S1043：当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰 Bc进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息，其中，n、m、i、c均为大于或等于1的正整数。
在S1043中的遥感图像的光谱排布情况参考图5，在图5中示出了光谱反射峰B1、B4和光谱吸收峰B2、B3。
所述蚀变信息获取方法通过对传统的主成分分析法进行改进，在对确定的待使用卫星数据进行预处理后，根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式，最后根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息的过程中分为三种情况，进行一次或多次主成分分析，以实现直接利用主成分分析法获得待识别蚀变矿物的信息的目的，避免了在深切割地区的遥感图像中进行除法运算，从而避免了部分阴影区会得到过分增强的可能，进而避免了对深切割地区的遥感图像进行蚀变信息提取时出现异常情况；而当遥感图像为覆盖严重的地区时，通过三种不同情况下的一次或多次主成分分析，可以将待识别蚀变矿物的信息集中到最终获得的第一主成分或第二主成分中，从而可以实现良好的蚀变信息提取效果。
在上述实施例的基础上，在本申请的一个实施例中，所述待使用卫星数据的种类为ASTER卫星数据或Landsat 8OLI卫星数据或哨兵2号卫星数据或高分1号卫星数据或WorldView-3卫星数据。
下面对其中几种卫星数据进行介绍：
Landsat 8OLI卫星数据包含11个波段，具体参考表1：
表1Landsat 8OLI数据参数表
高分1号卫星数据具体参数参考表2：
表2GF-1号系统参数表
ASTER卫星数据:ASTER数据在可见光-近红外(VNIR)范围有4个波段，空间分辨率为15m；在短波红外(SWIR)范围有6个波段，分辨率为30m；在热红外(TIR)范围有5个波段，分辨率为90m。
在上述实施例的基础上，在本申请的另一个实施例中，所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或Fe3+铁染蚀变矿物或Fe2+铁染蚀变矿物或明矾石蚀变矿物或铁帽或绿泥石和碳酸盐混合蚀变矿物。
在上述实施例的基础上，在本申请的一个具体实施例中，所述根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式包括：
当所述待使用卫星数据的种类为ASTER卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为铁帽的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为Landsat 8OLI卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为WorldView-3卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe2+铁染蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
为了检验本申请实施例提供的蚀变信息获取方法的实际效果，下面通过具体实施例来说明本申请实施例提供的蚀变信息获取方法的实施效果。
实施例1：
在本实施例中，遥感图像为龙巴其保研究区的哨兵2号图像(6593行× 5923列，13个波段)，通过实验对比，根据该遥感图像确定的待使用卫星数据为哨兵2号卫星数据；
图6(a)是原始遥感图像的子区(大小：1430行×930列)图像，图6 (b)为利用现有技术的主成分分析法提取的Fe3+铁染矿物图像，图6(c)为利用现有技术中的波段比值法提取的Fe3+铁染矿物图像，图6(d)为利用本申请实施例提供的蚀变信息获取方法提取的Fe3+铁染矿物图像。在进行主成分变换时，相关的参数统计只在无云雪覆盖的区域进行。在波段比值图像(图 6(c))中，中部突出的蚀变条带位于阳坡，是正地形，不会是阴影造成的虚假蚀变信息。从表3所示的主成分分析转换矩阵可以判断，第一主成分PC1 代表地表亮度，第二主成分PC2主要代表波段11与其它波段的反差，第三主成分PC3(图6(b))主要代表波段2和波段8的反差，第四主成分PC4主要代表了波段4和波段8的反差。4个主成分(PC)都与波段比值图像有明显差别，无一主成分代表波段2与4的反差，即主要包含Fe3+铁染信息。
表3波段2、4、8、11的主成分分析转换矩阵
在本实施例中，研究区如果属于深切割地区，地形阴影会比较明显。利用常规的波段比值来提取Fe3+铁染矿物信息时，部分阴影区会得到异常增强，并与真正的铁染信息较强的区域相混淆，严重误导铁染蚀变信息的提取。如果研究区内覆盖严重，铁染或羟基蚀变信息会被严重遮盖，在主成分分析中相关信息可能会分散在多个主成分中，不能集中于某一高阶主成分中。在此情况下，这类蚀变信息就难以用常规的主成分分析法来有效地提取。而用本申请实施例提供的蚀变信息提取方法得到的Fe3+铁染矿物信息和用波段比值得到的Fe3+蚀变信息(图6(c))相一致，中部突出的蚀变条带位于阳坡，是正地形，不会是阴影造成的虚假蚀变信息。
实施例2：
利用本申请实施例提供的蚀变信息的获取方法也可以有效地提取羟基蚀变信息。图7(a)是龙巴其保研究区的哨兵2号图像的子区(大小：1430行× 930列)图像，图7(b)为利用现有技术的主成分分析法提取的羟基矿物的图像，图7(c)为利用现有技术中的波段比值法提取的羟基矿物图像，图7(d) 为利用本申请实施例提供的蚀变信息获取方法提取的羟基矿物图像；从图中可以看出，3个羟基蚀变信息图像(图7(b)、图7(c)和图7(d))非常类似，都有效地提取出来图像中的羟基蚀变条带。
下面对本申请实施例提供的蚀变信息获取系统进行描述，下文描述的蚀变信息获取系统与上文描述的蚀变信息获取方法可相互对应参照。
一种蚀变信息获取系统，用于遥感图像的蚀变矿物的信息提取，所述蚀变信息获取系统包括：
图像获取模块，用于获取遥感图像，根据所述遥感图像确定待使用卫星数据的种类；
预处理模块，用于对确定的待使用卫星数据进行预处理；
公式确定模块，用于根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式；
蚀变提取模块，用于根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息；
所述蚀变提取模块根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息具体用于，
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行主成分分析，主成分分析结果中的第二主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bm进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bc进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息，其中，n、m、i、c均为大于或等于1的正整数。
可选的，所述待使用卫星数据的种类为ASTER卫星数据或Landsat 8OLI 卫星数据或哨兵2号卫星数据或高分1号卫星数据或WorldView-3卫星数据。
可选的，所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或Fe3+铁染蚀变矿物或Fe2+铁染蚀变矿物或铁帽蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
可选的，所述公式确定模块具体用于，当所述待使用卫星数据的种类为 ASTER卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为 Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为铁帽的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为Landsat 8OLI卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为WorldView-3卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe2+铁染蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
综上所述，本申请实施例提供了一种蚀变信息获取方法及系统，其中，所述蚀变信息获取方法通过对传统的主成分分析法进行改进，在对确定的待使用卫星数据进行预处理后，根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式，最后根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的过程中分为三种情况，进行一次或多次主成分分析，以实现直接利用主成分分析法获得待识别蚀变矿物的目的，避免了在深切割地区的遥感图像中进行波段除法运算，从而避免了部分阴影区会得到极大增强的可能，进而避免了对深切割地区的遥感图像进行蚀变信息提取时出现的高异常情况；而当遥感图像所对应的地区为覆盖严重的地区时，通过三种不同情况下的一次或多次主成分分析，可以将待识别蚀变矿物的信息集中到最终获得的第一主成分或第二主成分中，从而实现良好的蚀变信息的提取效果。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下，在其它实施例中实现。因此，本申请将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (8)

1.一种蚀变信息获取方法，其特征在于，用于遥感图像的蚀变矿物的信息提取，所述蚀变信息获取方法包括：
获取遥感图像，根据所述遥感图像确定待使用卫星数据的种类；
对确定的待使用卫星数据进行预处理；
根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式；
根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息；
所述根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息包括：
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行主成分分析，主成分分析结果中的第二主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bm进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bc进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的变信息，其中，n、m、i、c均为大于或等于1的正整数。
2.根据权利要求1所述的方法，其特征在于，所述待使用卫星数据的种类为ASTER卫星数据或Landsat 8 OLI卫星数据或哨兵2号卫星数据或高分1号卫星数据或WorldView-3卫星数据。
3.根据权利要求2所述的方法，其特征在于，所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或Fe3+铁染蚀变矿物或Fe2+铁染蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
4.根据权利要求3所述的方法，其特征在于，所述根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式包括：
当所述待使用卫星数据的种类为ASTER卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为铁帽的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为Landsat 8 OLI卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为WorldView-3卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe2+铁染蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
5.一种蚀变信息获取系统，其特征在于，用于遥感图像的蚀变矿物的信息提取，所述蚀变信息获取系统包括：
图像获取模块，用于获取遥感图像，根据所述遥感图像确定待使用卫星数据的种类；
预处理模块，用于对确定的待使用卫星数据进行预处理；
公式确定模块，用于根据待使用卫星数据，确定待识别蚀变矿物的波段比值提取公式；
蚀变提取模块，用于根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息；
所述蚀变提取模块根据确定的待识别蚀变矿物的波段比值提取公式，利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息具体用于，
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行主成分分析，主成分分析结果中的第二主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bm进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息；
当待识别蚀变矿物的波段比值提取公式为时，对光谱反射峰Bn和光谱吸收峰Bm进行第一次主成分分析，对光谱反射峰Bi和光谱吸收峰Bc进行第二次主成分分析，对第一次主成分分析结果中的第二主成分和第二次主成分分析结果中的第二主成分进行第三次主成分分析，第三次主成分分析结果中的第一主成分为所述待识别蚀变矿物的信息，其中，n、m、i、c均为大于或等于1的正整数。
6.根据权利要求5所述的系统，其特征在于，所述待使用卫星数据的种类为ASTER卫星数据或Landsat 8 OLI卫星数据或哨兵2号卫星数据或高分1号卫星数据或WorldView-3卫星数据。
7.根据权利要求6所述的系统，其特征在于，所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或Fe3+铁染蚀变矿物或Fe2+铁染蚀变矿物或铁帽蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
8.根据权利要求7所述的系统，其特征在于，所述公式确定模块具体用于，当所述待使用卫星数据的种类为ASTER卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为铁帽蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为Landsat 8 OLI卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为Fe2+铁染蚀变矿物的波段比值提取公式为当所述待使用卫星数据的种类为高分1号卫星数据时，Fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为WorldView-3卫星数据时，铝羟基蚀变矿物的波段比值提取公式为镁羟基蚀变矿物的波段比值提取公式为粘土蚀变矿物的波段比值为Fe2+铁染蚀变矿物的波段比值提取公式为明矾石蚀变矿物的波段比值提取公式为绿泥石和碳酸盐混合蚀变矿物的波段比值提取公式为
CN201810145611.3A 2018-02-12 2018-02-12 一种蚀变信息获取方法及系统 Active CN108333129B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810145611.3A CN108333129B (zh) 2018-02-12 2018-02-12 一种蚀变信息获取方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810145611.3A CN108333129B (zh) 2018-02-12 2018-02-12 一种蚀变信息获取方法及系统

Publications (2)

Publication Number Publication Date
CN108333129A true CN108333129A (zh) 2018-07-27
CN108333129B CN108333129B (zh) 2020-04-21

Family

ID=62929249

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810145611.3A Active CN108333129B (zh) 2018-02-12 2018-02-12 一种蚀变信息获取方法及系统

Country Status (1)

Country Link
CN (1) CN108333129B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110618106A (zh) * 2019-08-12 2019-12-27 中国地质调查局西安地质调查中心 一种基于近红外反射光谱的绿泥石矿物种类鉴定方法
CN113406041A (zh) * 2021-05-31 2021-09-17 核工业北京地质研究院 一种钠交代岩型铀矿关键蚀变矿物组合获取方法
CN114596430A (zh) * 2022-05-07 2022-06-07 交通运输通信信息集团有限公司 一种岩石蚀变信息提取方法、系统及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103984940A (zh) * 2014-06-03 2014-08-13 核工业北京地质研究院 一种基于高光谱数据识别赤铁矿化的方法
CN104050252A (zh) * 2014-06-12 2014-09-17 核工业北京地质研究院 一种高光谱遥感蚀变信息提取方法
CN104991287A (zh) * 2015-07-27 2015-10-21 中国地质调查局武汉地质调查中心 中南半岛成矿带典型区铜金矿遥感找矿模型建立方法
CN105068136A (zh) * 2015-07-27 2015-11-18 中国地质调查局武汉地质调查中心 基于多源信息的中南半岛示范区铜金矿潜力定位评价方法
CN105574621A (zh) * 2016-01-18 2016-05-11 中国地质科学院矿产资源研究所 基于遥感蚀变异常的斑岩铜矿预测系统及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103984940A (zh) * 2014-06-03 2014-08-13 核工业北京地质研究院 一种基于高光谱数据识别赤铁矿化的方法
CN104050252A (zh) * 2014-06-12 2014-09-17 核工业北京地质研究院 一种高光谱遥感蚀变信息提取方法
CN104991287A (zh) * 2015-07-27 2015-10-21 中国地质调查局武汉地质调查中心 中南半岛成矿带典型区铜金矿遥感找矿模型建立方法
CN105068136A (zh) * 2015-07-27 2015-11-18 中国地质调查局武汉地质调查中心 基于多源信息的中南半岛示范区铜金矿潜力定位评价方法
CN105574621A (zh) * 2016-01-18 2016-05-11 中国地质科学院矿产资源研究所 基于遥感蚀变异常的斑岩铜矿预测系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王彪: "遥感蚀变信息定量提取方法在成矿预测中的应用", 《大地构造与成矿学》 *
王烜: "基于主成分分析法对遥感蚀变异常提取", 《地质与资源》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110618106A (zh) * 2019-08-12 2019-12-27 中国地质调查局西安地质调查中心 一种基于近红外反射光谱的绿泥石矿物种类鉴定方法
CN110618106B (zh) * 2019-08-12 2022-04-01 中国地质调查局西安地质调查中心 一种基于近红外反射光谱的绿泥石矿物种类鉴定方法
CN113406041A (zh) * 2021-05-31 2021-09-17 核工业北京地质研究院 一种钠交代岩型铀矿关键蚀变矿物组合获取方法
CN114596430A (zh) * 2022-05-07 2022-06-07 交通运输通信信息集团有限公司 一种岩石蚀变信息提取方法、系统及介质

Also Published As

Publication number Publication date
CN108333129B (zh) 2020-04-21

Similar Documents

Publication Publication Date Title
Crowley et al. Airborne imaging spectrometer data of the Ruby Mountains, Montana: mineral discrimination using relative absorption band-depth images
Abdullah et al. Landsat ETM-7 for lineament mapping using automatic extraction technique in the SW part of Taiz Area, Yemen
Ismail et al. Satellite data classification accuracy assessment based from reference dataset
CN108333129A (zh) 一种蚀变信息获取方法及系统
CN110390267A (zh) 一种基于高分遥感影像的山地景观建筑提取方法和设备
Pasher et al. Multivariate forest structure modelling and mapping using high resolution airborne imagery and topographic information
CN105989322B (zh) 一种基于高分辨率遥感影像的多指标融合滑坡检测方法
Zhang et al. Spectral unmixing of normalized reflectance data for the deconvolution of lichen and rock mixtures
Trichon et al. Identifying spatial patterns in the tropical rain forest structure using hemispherical photographs
CN103761752A (zh) 一种极化合成孔径雷达图像的处理方法及装置
Nguyen Water body extraction from multi spectral image by spectral pattern analysis
CN103063311B (zh) 基于土壤指数的裸露基岩信息提取方法
CN101836233A (zh) 房屋变动判断方法、房屋变动判断程序、房屋变动判断用图像生成方法及房屋变动判断用图像
CN108171210A (zh) 一种植被覆盖区蚀变遥感异常信息的提取方法及系统
DE202014010866U1 (de) Spekularitätsbestimmung von Bildern
CN110533052A (zh) 一种协同遥感影像的航拍相片植被信息提取方法
Jing et al. Automated tree crown delineation from imagery based on morphological techniques
SA517390572B1 (ar) تحديد مادة حشو دعمي متحركة
CN106885772A (zh) 一种地质知识参与的etm遥感羟基蚀变矿物提取改进方法
Laake Structural mapping with spectral attributes
Harris et al. Modern isolated carbonate platforms: templates for quantifying facies attributes of hydrocarbon reservoirs
Tsendbazar Object based image analysis of geo-eye VHR data to model above ground carbon stock in Himalayan mid-hill forests, Nepal
Xu et al. Advanced exploratory data analysis for mapping regional canopy cover
JP2004061220A (ja) 地下水の探査方法
Ji et al. Subpixel measurement of mangrove canopy closure via spectral mixture analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant