class CsvCatalogueParserTestCase(unittest.TestCase):
    """ Unit tests for the csv Catalogue Parser Class"""
    BASE_DATA_PATH = os.path.join(os.path.dirname(__file__), 'data')

    def setUp(self):
        self.gardner_knopoff_window = GardnerKnopoffWindow()
        self.gruenthal_window = GruenthalWindow()
        self.uhrhammer_window = UhrhammerWindow()

    def test_gardner_knopoff_window(self):
        """
        Test the Gardner and Knopoff Distance-Time window
        """
        mag = np.array([5.0, 6.6])
        sw_space, sw_time = self.gardner_knopoff_window.calc(mag)
        self.assertAlmostEqual(sw_space[0], 39.99447, places=5)
        self.assertAlmostEqual(sw_space[1], 63.10736, places=5)
        self.assertAlmostEqual(sw_time[0], 143.71430 / 364.75, places=5)
        self.assertAlmostEqual(sw_time[1], 891.45618 / 364.75, places=5)

    def test_gruenthal_window(self):
        """
        Test the Gruenthal Distance-Time window
        """
        mag = np.array([5.0, 6.6])
        sw_space, sw_time = self.gruenthal_window.calc(mag)
        #self.assertAlmostEqual(sw_space[0], 39.99447, places=5)
        #self.assertAlmostEqual(sw_time[0], 143.71430/364.75, places=5)
        #self.assertAlmostEqual(sw_space[1], 63.10736, places=5)
        #self.assertAlmostEqual(sw_time[1], 891.45618/364.75, places=5)

    def test_uhrhammer_window(self):
        mag = np.array([5.0, 6.6])
        sw_space, sw_time = self.uhrhammer_window.calc(mag)
class CsvCatalogueParserTestCase(unittest.TestCase):
    """ Unit tests for the csv Catalogue Parser Class"""

    BASE_DATA_PATH = os.path.join(os.path.dirname(__file__), "data")

    def setUp(self):
        self.gardner_knopoff_window = GardnerKnopoffWindow()
        self.gruenthal_window = GruenthalWindow()
        self.uhrhammer_window = UhrhammerWindow()

    def test_gardner_knopoff_window(self):
        """
        Test the Gardner and Knopoff Distance-Time window
        """
        mag = np.array([5.0, 6.6])
        sw_space, sw_time = self.gardner_knopoff_window.calc(mag)
        self.assertAlmostEqual(sw_space[0], 39.99447, places=5)
        self.assertAlmostEqual(sw_space[1], 63.10736, places=5)
        self.assertAlmostEqual(sw_time[0], 143.71430 / 364.75, places=5)
        self.assertAlmostEqual(sw_time[1], 891.45618 / 364.75, places=5)

    def test_gruenthal_window(self):
        """
        Test the Gruenthal Distance-Time window
        """
        mag = np.array([5.0, 6.6])
        sw_space, sw_time = self.gruenthal_window.calc(mag)
        # self.assertAlmostEqual(sw_space[0], 39.99447, places=5)
        # self.assertAlmostEqual(sw_time[0], 143.71430/364.75, places=5)
        # self.assertAlmostEqual(sw_space[1], 63.10736, places=5)
        # self.assertAlmostEqual(sw_time[1], 891.45618/364.75, places=5)

    def test_uhrhammer_window(self):
        mag = np.array([5.0, 6.6])
        sw_space, sw_time = self.uhrhammer_window.calc(mag)
Exemple #3
0
                  'figname': output_base+"_pp_decluster.png",
                  }
        declustered_catalogue = decluster_catalogue(catalogue, config)
        plot_rate(catalogue, cumulative=True, color='#FF9D84', new_figure=False, overlay=True, linewidth=1.5, label="Gardner-Knopoff/Gardner-Knopoff", alpha=0.8)

        
        catalogue = read_catalog(input_catalogue_file)
        valid_months = catalogue.data['day'] != 0
        catalogue.select_catalogue_events(valid_months)
        valid_years = catalogue.data['year'] >= 1900
        catalogue.select_catalogue_events(valid_years)
        config = {
                  'decluster_method': 'gardner_knopoff', # afteran | gardner_knopoff
                  #'time_distance_window': GardnerKnopoffWindow(),
                  #'time_distance_window': GruenthalWindow(),
                  'time_distance_window': UhrhammerWindow(),
                  # afteran
                  'time_window': 60.,
                  # gardner_knopoff
                  'fs_time_prop': 1.0,
                  # write declustered
                  'filename': None,
                  #'filename': output_base+"_pp_decluster.csv",
                  # plot cluster
                  'plot':False,
                  'figname': output_base+"_pp_decluster.png",
                  }
        declustered_catalogue = decluster_catalogue(catalogue, config)
        plot_rate(catalogue, cumulative=True, color='#FF5460', new_figure=False, overlay=True, linewidth=1.5, label="Gardner-Knopoff/Uhrhammer", alpha=0.8)

        
 def setUp(self):
     self.gardner_knopoff_window = GardnerKnopoffWindow()
     self.gruenthal_window = GruenthalWindow()
     self.uhrhammer_window = UhrhammerWindow()
 def setUp(self):
     self.gardner_knopoff_window = GardnerKnopoffWindow()
     self.gruenthal_window = GruenthalWindow()
     self.uhrhammer_window = UhrhammerWindow()