Skip to content

Clipping test gives false negative (algorithm false to detect clipping)

Some traces are passing the clipping detection test when they should not be. I assume the histogram algorithm should detect this clipping.

Example data: event.json KO.KIZT.xml KO.KIZT..HHE.mseed KO.KIZT..HHN.mseed KO.KIZT..HHZ.mseed

Code to run clipping test for example data:

#!/usr/bin/env python3

import pathlib
import obspy
import json

from gmprocess.core.stationstream import StationStream
from gmprocess.waveform_processing.clipping.clipping_check import check_clipping
from gmprocess.utils.event import get_event_object


EVENT = "us6000jlqa"
STATION = "KO.KIZT"

traces = obspy.read(f"{STATION}*.mseed")
inventory = obspy.read_inventory(f"{STATION}.xml")
station_stream = StationStream(traces, inventory)

with open("event.json") as fin:
    ev_info = json.load(fin)
event = get_event_object(ev_info)

import pdb; pdb.set_trace()
check_clipping(station_stream, event)