CleanRotator example code

from obspy.core.stream import read
from obspy.core.inventory import read_inventory
from tiskitpy import CleanRotator, SpectralDensity

# Read the data and metadata, and calculate a CleanRotator
stream = read('data/XS.S11D.LH.2016.12.11.mseed', 'MSEED')
inv = read_inventory('data/XS.S11_decimated.station.xml', 'STATIONXML')
cr = CleanRotator(stream)

[INFO] Didn't find local EQ file '20161210-20161211_MM5.85_eqcat.qml', reading from USGS online catalog...
[INFO] Done
[INFO] writing catalog to "20161210-20161211_MM5.85_eqcat.qml"
[INFO] CleanRotator: angle, azimuth, var_red =  0.09,  212.3, 0.75

# Apply the CleanRotator to the data
stream_rotated = cr.apply(stream)

# Combine the rotated and unrotated data, for comparison
z_compare = stream_rotated.select(channel='*Z') + stream.select(channel='*Z')

# Print and plot the combined data
print(z_compare)
z_compare.plot()

2 Trace(s) in Stream:
XS.S11D.-ROT.LHZ | 2016-12-10T23:59:59.992583Z - 2016-12-11T23:59:59.992583Z | 1.0 Hz, 86401 samples
XS.S11D..LHZ     | 2016-12-10T23:59:59.992583Z - 2016-12-11T23:59:59.992583Z | 1.0 Hz, 86401 samples

# Create and plot a SpectralDensity object for the rotated and unrotated data
sd_compare = SpectralDensity.from_stream(z_compare, inv=inv)
sd_compare.plot(overlay=True)

[INFO] z_threshold=3, rejected 2% of windows (2/84)