Source code for openquake.hmtk.seismicity.declusterer.distance_time_windows
# -*- coding: utf-8 -*-# vim: tabstop=4 shiftwidth=4 softtabstop=4## LICENSE## Copyright (C) 2010-2023 GEM Foundation, G. Weatherill, M. Pagani,# D. Monelli.## The Hazard Modeller's Toolkit is free software: you can redistribute# it and/or modify it under the terms of the GNU Affero General Public# License as published by the Free Software Foundation, either version# 3 of the License, or (at your option) any later version.## You should have received a copy of the GNU Affero General Public License# along with OpenQuake. If not, see <http://www.gnu.org/licenses/>## DISCLAIMER## The software Hazard Modeller's Toolkit (openquake.hmtk) provided herein# is released as a prototype implementation on behalf of# scientists and engineers working within the GEM Foundation (Global# Earthquake Model).## It is distributed for the purpose of open collaboration and in the# hope that it will be useful to the scientific, engineering, disaster# risk and software design communities.## The software is NOT distributed as part of GEM’s OpenQuake suite# (https://www.globalquakemodel.org/tools-products) and must be considered as a# separate entity. The software provided herein is designed and implemented# by scientific staff. It is not developed to the design standards, nor# subject to same level of critical review by professional software# developers, as GEM’s OpenQuake software suite.## Feedback and contribution to the software is welcome, and can be# directed to the hazard scientific staff of the GEM Model Facility# (hazard@globalquakemodel.org).## The Hazard Modeller's Toolkit (openquake.hmtk) is therefore distributed WITHOUT# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License# for more details.## The GEM Foundation, and the authors of the software, assume no# liability for use of the software.""""""importabcimportnumpyasnpfromopenquake.baselib.generalimportCallableDictTIME_DISTANCE_WINDOW_FUNCTIONS=CallableDict()DAYS=364.75
[docs]deftime_window_cutoff(sw_time,time_cutoff):""" Allows for cutting the declustering time window at a specific time, outside of which an event of any magnitude is no longer identified as a cluster """sw_time=np.array([(time_cutoff/DAYS)ifx>(time_cutoff/DAYS)elsexforxinsw_time])returnsw_time
[docs]classBaseDistanceTimeWindow(object):""" Defines the space and time windows, within which an event is identified as a cluster. """
[docs]@abc.abstractmethoddefcalc(self,magnitude,time_cutoff=None):""" Allows to calculate distance and time windows (sw_space, sw_time) see reference: `Van Stiphout et al (2011)`. :param magnitude: magnitude :type magnitude: numpy.ndarray :param time_cutoff: time window cutoff in days (optional) :type time_cutoff: int :returns: distance and time windows :rtype: numpy.ndarray """return
[docs]@TIME_DISTANCE_WINDOW_FUNCTIONS.add("GardnerKnopoff")classGardnerKnopoffWindow(BaseDistanceTimeWindow):""" Gardner Knopoff method for calculating distance and time windows """
[docs]@TIME_DISTANCE_WINDOW_FUNCTIONS.add("Gruenthal")classGruenthalWindow(BaseDistanceTimeWindow):""" Gruenthal method for calculating distance and time windows """
[docs]@TIME_DISTANCE_WINDOW_FUNCTIONS.add("UrhammerWindow")classUhrhammerWindow(BaseDistanceTimeWindow):""" Uhrhammer method for calculating distance and time windows """