PSN-L Email List Message
Subject: IRIS Wilber II
From: Bobhelenmcclure@.......
Date: Sat, 10 Feb 2007 20:11:15 EST
A memo to Larry Conklin and to others who might be interested.
Hi Larry,
You mentioned that you compare your data with PAL. I do also since they
are only 21km from me. Have you tried looking at BINY? They are closer to you.
WinQuake and Amaseis will open SAC binary files. You can download such
files from major events from almost any network station in the world easily from
IRIS at
_http://www.iris.washington.edu/_ (http://www.iris.washington.edu/)
Click on "Earthquakes" and "Wilber II" to get to the download page. BINY
belongs to the US network. Do not try to select an event using the picture
display. Go to the listings of events instead. After you have selected your
event and stations (you have to unselect almost everything), be sure you select
"SAC BINARY individual files" as the file format. Select BHZ, BHN, and BHE as
the channels. Give a user name and email address. Elect to have an email sent
when processing is complete. If you do not ask for many files, you will have
a high priority and will usually get your files prepared in a few minutes
after the request is entered.
I occasionally download events from PAL. My correspondence with their data
is remarkably close, once I have processed my data using my WQFilter
utility. Something is not right with the IIR filters in WinQuake. They have more
ringing and sharper cutoffs than they should, so I avoid using them. The FFT
filters are OK. Cochrane did not write the filter code himself. He lifted the
code from a book "Analog and Digital Filter Design using C" by Les Thede,
which is no longer in print.
The filters in Amaseis are OK. Both Dr. Jones and I worked out the
algorithms for Butterworth filters and ended up with identical results. I gave him
my algorithm for the period extending filter, and his filter performance
agrees with mine. I do have an optional back-and-forth filter pass capability,
which yields zero distortion filtering, equivalent to FFT filtering. Amaseis
does not implement that.
If you should want to convert PSN files to SAC binary format, WinQuake can
do that. However, there is a bug in WinQuake that prevents transfer of event
time. I am able to patch in a correction on the files I convert. I told
Cochrane what was wrong, but he has essentially retired from maintaining his
software, after a lifetime of good and dedicated work.
I am almost done with a heliplot program that filters and displays WinDaq
log files, including the one currently open, up to 27 hours before the
present hour or any lesser number of hours.
Regards,
Bob
A memo to Larry Conklin and to others who might be interested.
Hi Larry,
You mentioned that you compare your data with PAL. I do also sin=
ce=20
they are only 21km from me. Have you tried looking at BINY? They are closer=20=
to=20
you.
WinQuake and Amaseis will open SAC binary files. You can downloa=
d=20
such files from major events from almost any network station in the world ea=
sily=20
from IRIS at
Click on "Earthquakes" and "Wilber II" to get to the download=20
page. BINY belongs to the US network. Do not try to select an event usi=
ng=20
the picture display. Go to the listings of events instead. After you have=20
selected your event and stations (you have to unselect almost everything), b=
e=20
sure you select "SAC BINARY individual files" as the file format. Select BHZ=
,=20
BHN, and BHE as the channels. Give a user name and email address. Elect to h=
ave=20
an email sent when processing is complete. If you do not ask for many files,=
you=20
will have a high priority and will usually get your files prepared in a few=20
minutes after the request is entered.
I occasionally download events from PAL. My correspondence with=20
their data is remarkably close, once I have processed my data using my WQFil=
ter=20
utility. Something is not right with the IIR filters in WinQuake. They have=20=
more=20
ringing and sharper cutoffs than they should, so I avoid using them. The FFT=
=20
filters are OK. Cochrane did not write the filter code himself. He lifted th=
e=20
code from a book "Analog and Digital Filter Design using C" by Le=
s=20
Thede, which is no longer in print.
The filters in Amaseis are OK. Both Dr. Jones and I worked=20=
out=20
the algorithms for Butterworth filters and ended up with identical results.=20=
I=20
gave him my algorithm for the period extending filter, and his filter=20
performance agrees with mine. I do have an optional back-and-forth filt=
er=20
pass capability, which yields zero distortion filtering, equivalent to FFT=20
filtering. Amaseis does not implement that.
If you should want to convert PSN files to SAC binary forma=
t,=20
WinQuake can do that. However, there is a bug in WinQuake that prevents tran=
sfer=20
of event time. I am able to patch in a correction on the files I convert. I=20=
told=20
Cochrane what was wrong, but he has essentially retired=20
from maintaining his software, after a lifetime of good and dedicated=20
work.
I am almost done with a heliplot program that filters and displa=
ys=20
WinDaq log files, including the one currently open, up to 27 hours before th=
e=20
present hour or any lesser number of hours.
Regards,
Bob
[ Top ]
[ Back ]
[ Home Page ]