Dear all,
rotation from ZNE to LQT components, looks working in my computer, but after all, I can not see the any changes in my wave amplitudes when compare the original ZNE and new LQT. Am I doing something wrong ?
here is my code and sample data,could you give a try if it is possible ?
best wishes,
Metin
2009015_ACPY_17_49.BHE.sac (18.1 KB)
2009015_ACPY_17_49.BHN.sac (18.1 KB)
2009015_ACPY_17_49.BHZ.sac (18.1 KB)
Your waveform data has strange channel codes, certainly not adhering to
seed standard (as implied by the file names):
st = read('2009015_ACPY_17_49*')
print st
.ACPY..E-W | 2009-01-15T18:11:19.650024Z - 2009-01-15T18:13:11.625024Z
40.0 Hz, 4480 samples
.ACPY..N-S | 2009-01-15T18:11:19.650024Z - 2009-01-15T18:13:11.625024Z
40.0 Hz, 4480 samples
.ACPY..VERT | 2009-01-15T18:11:19.650024Z - 2009-01-15T18:13:11.625024Z
40.0 Hz, 4480 samples
Therefore the rotate method can't find the right traces to use. You can
either fix the channel codes in the traces manually (by setting e.g.
st[0].stats.channel) or fall back to the low level routine from
obspy.signal that you were trying to use in your first post.
best,
Tobias