本申请涉及遥感图像处理技术领域,更具体地说,涉及一种蚀变信息获取方法及系统。
背景技术:
蚀变作用是指当流体或热液作用于岩石时,两者之间必然趋于达到化学平衡和发生物质交换作用,其结果导致岩石中部分或全部原矿物的消失和新矿物的形成,即原岩中某些物质部分或全部地带出和新物质带入的交代作用。蚀变信息对于矿物探测具有重大的意义,因此,通过提取遥感地图中的蚀变信息是矿物探测过程中的重要一环。
从遥感图像中提取蚀变信息的现有方法主要包括波段比值法和主成分分析法,例如crosta技术。其中,在使用波段比值法提取深切割地区的遥感图像的蚀变信息时,可能会在波段比值过程中产生异常结果。这是因为深切割地区的海拔落差较大,山体坡度较陡,遥感图像中含有明显的地形阴影,如果遥感图像预处理过程中的大气校正不甚精确,遥感图像中阴影的某些波段的值有可能接近于0,从而在波段比值过程中产生异常高的结果,从而导致蚀变信息的提取产生较大的误差。在使用crosta技术提取覆盖地区的遥感图像的蚀变矿物时,也难以达到良好的提取效果。这是因为在这些地区,矿物表面覆盖了大量的杂质(土、沙、砂石、植被等)。这些杂质遮盖了蚀变矿物的信息,使得蚀变矿物的信息在遥感图像中显示较弱,从而难以在主成分分析过程中将蚀变信息集中到单一的主成分中,进而导致蚀变矿物信息提取的效果较差。
技术实现要素:
为解决上述技术问题,本申请提供了一种蚀变信息获取方法及系统,以解决在提取深切割地区的遥感图像的蚀变信息提取时可能产生的异常结果,以及在提取覆盖严重的地区的遥感图像的蚀变信息时,提取效果较差的问题。
为解决上述技术问题,本申请实施例提供了如下技术方案:
一种蚀变信息获取方法,用于遥感图像的蚀变矿物的信息提取,所述蚀变信息获取方法包括:
获取遥感图像,根据所述遥感图像确定待使用卫星数据的种类;
对确定的待使用卫星数据进行预处理;
根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式;
根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息;
所述根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息包括:
当待识别蚀变矿物的波段比值提取公式为
当待识别蚀变矿物的波段比值提取公式为
当待识别蚀变矿物的波段比值提取公式为
可选的,所述待使用卫星数据的种类为aster卫星数据或landsat8oli卫星数据或哨兵2号卫星数据或高分1号卫星数据或worldview-3卫星数据。
可选的,所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或fe3+铁染蚀变矿物或fe2+铁染蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
可选的,所述根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式包括:
当所述待使用卫星数据的种类为aster卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为landsat8oli卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为worldview-3卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
一种蚀变信息获取系统,用于遥感图像的蚀变矿物的信息提取,所述蚀变信息获取系统包括:
图像获取模块,用于获取遥感图像,根据所述遥感图像确定待使用卫星数据的种类;
预处理模块,用于对确定的待使用卫星数据进行预处理;
公式确定模块,用于根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式;
蚀变提取模块,用于根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息;
所述蚀变提取模块根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息具体用于,
当待识别蚀变矿物的波段比值提取公式为
当待识别蚀变矿物的波段比值提取公式为
当待识别蚀变矿物的波段比值提取公式为
可选的,所述待使用卫星数据的种类为aster卫星数据或landsat8oli卫星数据或哨兵2号卫星数据或高分1号卫星数据或worldview-3卫星数据。
可选的,所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或fe3+铁染蚀变矿物或fe2+铁染蚀变矿物或铁帽蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
可选的,所述公式确定模块具体用于,当所述待使用卫星数据的种类为aster卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为landsat8oli卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为worldview-3卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
从上述技术方案可以看出,本申请实施例提供了一种蚀变信息获取方法及系统,其中,所述蚀变信息获取方法通过对传统的主成分分析法进行改进,在对确定的待使用卫星数据进行预处理后,根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式,最后根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的过程中分为三种情况,进行一次或多次主成分分析,以实现直接利用主成分分析法获得待识别蚀变矿物的目的,避免了在深切割地区的遥感图像中进行波段除法运算,从而避免了部分阴影区会得到极大增强的可能,进而避免了对深切割地区的遥感图像进行蚀变信息提取时出现的高异常情况;而当遥感图像所对应的地区为覆盖严重的地区时,通过三种不同情况下的一次或多次主成分分析,可以将待识别蚀变矿物的信息集中到最终获得的第一主成分或第二主成分中,从而实现良好的蚀变信息的提取效果。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本申请的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其它的附图。
图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:当待识别蚀变矿物的波段比值提取公式为
在s1041中的遥感图像的光谱排布情况参考图3,在图3中示出了光谱反射峰b1和光谱吸收峰b2。
s1042:当待识别蚀变矿物的波段比值提取公式为
在s1042中的遥感图像的光谱排布情况参考图4,在图4中示出了光谱反射峰b1、b3和光谱吸收峰b2。
s1043:当待识别蚀变矿物的波段比值提取公式为
在s1043中的遥感图像的光谱排布情况参考图5,在图5中示出了光谱反射峰b1、b4和光谱吸收峰b2、b3。
所述蚀变信息获取方法通过对传统的主成分分析法进行改进,在对确定的待使用卫星数据进行预处理后,根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式,最后根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息的过程中分为三种情况,进行一次或多次主成分分析,以实现直接利用主成分分析法获得待识别蚀变矿物的信息的目的,避免了在深切割地区的遥感图像中进行除法运算,从而避免了部分阴影区会得到过分增强的可能,进而避免了对深切割地区的遥感图像进行蚀变信息提取时出现异常情况;而当遥感图像为覆盖严重的地区时,通过三种不同情况下的一次或多次主成分分析,可以将待识别蚀变矿物的信息集中到最终获得的第一主成分或第二主成分中,从而可以实现良好的蚀变信息提取效果。
在上述实施例的基础上,在本申请的一个实施例中,所述待使用卫星数据的种类为aster卫星数据或landsat8oli卫星数据或哨兵2号卫星数据或高分1号卫星数据或worldview-3卫星数据。
下面对其中几种卫星数据进行介绍:
landsat8oli卫星数据包含11个波段,具体参考表1:
表1landsat8oli数据参数表
高分1号卫星数据具体参数参考表2:
表2gf-1号系统参数表
aster卫星数据:aster数据在可见光-近红外(vnir)范围有4个波段,空间分辨率为15m;在短波红外(swir)范围有6个波段,分辨率为30m;在热红外(tir)范围有5个波段,分辨率为90m。
在上述实施例的基础上,在本申请的另一个实施例中,所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或fe3+铁染蚀变矿物或fe2+铁染蚀变矿物或明矾石蚀变矿物或铁帽或绿泥石和碳酸盐混合蚀变矿物。
在上述实施例的基础上,在本申请的一个具体实施例中,所述根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式包括:
当所述待使用卫星数据的种类为aster卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为landsat8oli卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为worldview-3卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
为了检验本申请实施例提供的蚀变信息获取方法的实际效果,下面通过具体实施例来说明本申请实施例提供的蚀变信息获取方法的实施效果。
实施例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))非常类似,都有效地提取出来图像中的羟基蚀变条带。
下面对本申请实施例提供的蚀变信息获取系统进行描述,下文描述的蚀变信息获取系统与上文描述的蚀变信息获取方法可相互对应参照。
一种蚀变信息获取系统,用于遥感图像的蚀变矿物的信息提取,所述蚀变信息获取系统包括:
图像获取模块,用于获取遥感图像,根据所述遥感图像确定待使用卫星数据的种类;
预处理模块,用于对确定的待使用卫星数据进行预处理;
公式确定模块,用于根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式;
蚀变提取模块,用于根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息;
所述蚀变提取模块根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的信息具体用于,
当待识别蚀变矿物的波段比值提取公式为
当待识别蚀变矿物的波段比值提取公式为
当待识别蚀变矿物的波段比值提取公式为
可选的,所述待使用卫星数据的种类为aster卫星数据或landsat8oli卫星数据或哨兵2号卫星数据或高分1号卫星数据或worldview-3卫星数据。
可选的,所述待识别蚀变矿物为铝羟基蚀变矿物或镁羟基蚀变矿物或铁羟基蚀变矿物或粘土蚀变矿物或fe3+铁染蚀变矿物或fe2+铁染蚀变矿物或铁帽蚀变矿物或明矾石蚀变矿物或绿泥石和碳酸盐混合蚀变矿物。
可选的,所述公式确定模块具体用于,当所述待使用卫星数据的种类为aster卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为landsat8oli卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为哨兵2号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为高分1号卫星数据时,fe3+铁染蚀变矿物的波段比值提取公式为
当所述待使用卫星数据的种类为worldview-3卫星数据时,铝羟基蚀变矿物的波段比值提取公式为
综上所述,本申请实施例提供了一种蚀变信息获取方法及系统,其中,所述蚀变信息获取方法通过对传统的主成分分析法进行改进,在对确定的待使用卫星数据进行预处理后,根据待使用卫星数据,确定待识别蚀变矿物的波段比值提取公式,最后根据确定的待识别蚀变矿物的波段比值提取公式,利用主成分分析法确定所述遥感图像中的待识别蚀变矿物的过程中分为三种情况,进行一次或多次主成分分析,以实现直接利用主成分分析法获得待识别蚀变矿物的目的,避免了在深切割地区的遥感图像中进行波段除法运算,从而避免了部分阴影区会得到极大增强的可能,进而避免了对深切割地区的遥感图像进行蚀变信息提取时出现的高异常情况;而当遥感图像所对应的地区为覆盖严重的地区时,通过三种不同情况下的一次或多次主成分分析,可以将待识别蚀变矿物的信息集中到最终获得的第一主成分或第二主成分中,从而实现良好的蚀变信息的提取效果。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。