本发明属于异常数据提取,特别是涉及一种基于时空联合的震前热异常提取方法。
背景技术:
1、地震发生前普遍存在的红外辐射异常现象,是当前评估区域发震危险性的重要参数之一,然而,并非所有的地表红外异常都与构造活动或地震有关,如何排除非构造因素对地表热红外辐射的影响,从强噪声背景中提取微弱信号,是当前利用热红外遥感技术研究构造活动的难点。
2、国内外针对震前热异常的研究最初是应用卫星红外图像进行地震异常分析,主要用目视解译法获得遥感图像增温区域的演变趋势,同时结合已知的地震地质知识,经过对震例的研究总结出卫星红外亮温增温异常与地震时、空、强三要素的关系,这种方法无法消除很多非构造因素的影响,主观性较强。后来逐渐发展了基于数理统计与分析的方法和基于背景场的方法。基于数理统计的方法主要包括基于差值分析的异常提取方法和基于信号分析的异常提取方法,基于差值分析的方法虽然可以在一定程度上突出地震前异常的时间空间分布,但其仅适用于地震热异常与断裂带密切相关或者震前震后亮温有明显变化的异常提取工作中。而基于信号分析的异常提取方法可以对热红外信号进行处理,将时间序列数据分解成相互独立的频率分量,因此这种方法提取出的异常信息的物理意义比其他方法更为清晰,但这种方法对空间域的异常信息分析较少,并且需要具备一定的地理和地质学科的相关知识,其机理解释也需要加强。而基于背景场分析的热异常检测方法是通过建立基于多年数据的区域背景场并与实际观测数据进行差值分析提取震前热异常,例如空间域rst算法可以消除气象影响,提取出空间域震前热异常信息,但却忽略了地形地貌的影响。
3、综上所述,现有的震前热异常的提取方法虽然原理简单易于实现,但该类方法没有做到时空联合,即同时考虑到时间域异常信息和空间域异常信息,并同时做到剔除大气活动、人类活动、地表覆盖等非构造因素影响,以提取构造相关性更强的震前热异常。
技术实现思路
1、本发明的目的是提供一种基于时空联合的震前热异常提取方法,以解决上述现有技术存在的问题。
2、为实现上述目的,本发明提供了一种基于时空联合的震前热异常提取方法,包括以下步骤:
3、获取待检测地区的地表温度数据,基于小波变换法提取原地温度;
4、基于所述原地温度,采用涡度算法获取原地涡度;
5、基于所述原地涡度,采用历年同期rst算法获取震前热异常信息。
6、可选地,基于小波变换法提取原地温度的过程包括:
7、基于一维小波变换对所述地表温度数据进行分解,剔除太阳活动、地理位置和地表覆盖的影响,获取分解数据;
8、基于二维小波变换对所述分解数据进行处理,剔除非构造因素的影响,获取所述原地温度。
9、可选地,基于一维小波变换对所述地表温度数据进行分解的方法包括:
10、
11、其中,是子小波,是一维母小波,表示一维时间尺度,是平移因子,是尺度因子。
12、可选地,基于二维小波变换对所述分解数据进行处理的方法包括:
13、
14、其中,表示二维空间信息,表示位置信息,为二维母小波,为部分子小波,是尺度因子,、为平移因子。
15、可选地,采用涡度算法获取原地涡度的过程包括:
16、获取待计算位置的中心像素,采用涡度算法计算所述中心像素相对于四邻域像素的大小差异,基于所述大小差异获取热异常信号的位置信息。
17、可选地,采用涡度算法计算所述中心像素相对于四邻域像素的大小差异的方法包括:
18、
19、其中,为待计算位置的中心像素,为中心像素相对于四邻域像素的大小差异。
20、可选地,采用历年同期rst算法获取震前热异常信息的方法为:
21、
22、其中,为历年同期一系列有效值的均值,为处时的地表温度值;
23、
24、其中, r表示震前热异常信息,代表位置处,时刻的遥感图像的dn值,和是处的均值和标准差。
25、本发明的技术效果为:
26、本发明提出的震前热异常提取方法兼顾了空间域的异常信息和时间域的异常信息,并剔除了大气活动、人类活动和地表覆盖等非构造因素的影响,提取出的热异常既有更强的构造相关性和稳定性,并且在时间和空间上都具有一定的规律性,提高了震前热异常检测的精度,便于进一步的分析与研究。
1.一种基于时空联合的震前热异常提取方法,其特征在于,包括以下步骤:
2.根据权利要求1所述的基于时空联合的震前热异常提取方法,其特征在于,基于小波变换法提取原地温度的过程包括:
3.根据权利要求2所述的基于时空联合的震前热异常提取方法,其特征在于,基于一维小波变换对所述地表温度数据进行分解的方法包括:
4.根据权利要求2所述的基于时空联合的震前热异常提取方法,其特征在于,基于二维小波变换对所述分解数据进行处理的方法包括:
5.根据权利要求1所述的基于时空联合的震前热异常提取方法,其特征在于,采用涡度算法获取原地涡度的过程包括:
6.根据权利要求5所述的基于时空联合的震前热异常提取方法,其特征在于,采用涡度算法计算所述中心像素相对于四邻域像素的大小差异的方法包括:
7.根据权利要求1所述的基于时空联合的震前热异常提取方法,其特征在于,采用历年同期rst算法获取震前热异常信息的方法为: