analysis.make_more_incomplete

seismostats.analysis.make_more_incomplete(magnitudes: ndarray, times: array, delta_t: timedelta64 = np.timedelta64(60, 's'), return_idx: bool = False) → tuple[ndarray, ndarray]

Returns filtered magnitudes and times. Filters the magnitudes and times in the following way: If an earthquake is smaller than the previous one and less than delta_t away, the earthquake is removed.

Source:

E. Lippiello and G. Petrillo. Journal of Geophysical Research: Solid Earth, 129(2):e2023JB027849, 2024.

Parameters:

• magnitudes – Array of magnitudes, ordered in time (first entry is the earliest earthquake).

• times – Array of datetime objects of occurrence of each earthquake.

• delta_t – Time window in seconds to filter out events. Default is 60 seconds.

• return_idx – If True, the indices of the events that were kept are also returned.

Returns:

• magnitudes – Array of filtered magnitudes.

• times – Array of filtered datetime objects.

• idx – Array of indices of the events that were kept.