I believe if you have different xmls from different years, but there all in
the same inventory directory, it should attach the correct one
automatically. The code loops through the inventory director reading in
each xml response file and appending it to a list of station responses
(actually a list of obspy inventories). It then attaches this to current
stream of seismograms, with obspy's attach_response handling everything.
if response_format == "inventory":
firstinv = True
inventory = None
for file in files:
try:
inv = read_inventory(file)
if firstinv:
inventory = inv
firstinv = False
else:
inventory += inv
except:
traceback.print_exc()
pass
if inventory:
stream.attach_response(inventory)
stream.remove_response(output='VEL',
pre_filt=response_prefilt)
On Thu, Oct 6, 2016 at 12:26 PM, VICIC BLAZ [PHD0900016] <
BLAZ.VICIC(a)phd.units.it> wrote:
Hi,
I will use this correspondence for my question, since you already
mentioned removal of the response.
Do i have to change the code or something if I have different stationxmls
for different years, so they are read correctly, or does the code manage
this and i just put them into the folder?
Thanks
________________________________
From: msnoise-bounces(a)mailman-as.oma.be <msnoise-bounces(a)mailman-as.oma.be>
on behalf of Flinders, Ashton <aflinders(a)usgs.gov>
Sent: Thursday, October 6, 2016 7:38:40 PM
To: Python Package for Monitoring Seismic Velocity Changes using
AmbientSeismic Noise
Subject: Re: [MSNoise] MSNoise Digest, Vol 31, Issue 1
Thanks Esteban,
I have looked through the documentation extensively, as well as the SRL
paper, and a bit of grepping on the source code (e.g. it took me a while to
figure out why in the world I couldnt remove instrument responses using a
RESP or PZ file, since it is listed on the documentation page, until I
grepped the source code and saw those were not implemented.)
1) Yes. It is more efficient to compute the cross-correlation functions or
Green’s functions using 30-min chunks. Although you can change this feature
if you will. Segments are re-build afterwards.
This is only true if the daily segments are always preloaded into memory
(IO vs CC time). The SRL paper makes it sounds like this is the case, but I
want to make sure that this is always the case, even for processing many
stations/days or long segments. I also wasnt sure if these segmented CC
were used for any statistics (e.g. uncertainty in daily CC values).
2) You need to define the filters in order to compute the
cross-correlations and extract the coherent part of the ambient noise for a
given frequency range. MWCS is used to compute the moving window cross
spectral analysis after the cross-correlation procedure.
Yes, I understand that (the way its been written) filter values should be
defined before the CC (which really should only be the whiteneing band or
preprocessing), however I dont understand why in the web gui the whiten and
MWCS filters are put within the same group. This make me predefine a MWCS
filter band, along with my whiten band, even though the MWCS doesnt happed
for another two processing steps. This leads to confusion, for example;
What happens if I compute the CC's with a defined filter (as written, a
MWCS and whiten filter defined together) and then want to compute different
MWCS's later on when I actually get to the MWCS step. Can I just define a
new filter, with new mwcs_low mwcs_high values? Since the web gui forces me
to define a filter high low value (whiten, I assume) does it not use these?
Do I have to go back and recalculate the cross correlations everytime I
want to add an MWCS filter (which I wouldnt expect I have too)?
Similarly, what happens if I want to delete a set of processed MWCS or
stacks? I cant just delete the folder, since their are still pointers in
the sql database. Just reprocessing doesnt seem to work either.
thanks again
-ashton
On Thu, Oct 6, 2016 at 10:17 AM, Esteban Chaves <echfisica(a)gmail.com>
wrote:
Hi Ashton,
As an active user:
I highly recommend you to check the documentation:
http://msnoise.org/doc/
workflow.html#module-msnoise.s03compute_cc
<http://msnoise.org/doc/
workflow.html#module-msnoise.s03compute_cc>
I think the pre-processing and processing steps are well documented in
the
web site.
1) Yes. It is more efficient to compute the cross-correlation functions
or
Green’s functions using 30-min chunks. Although
you can change this
feature
if you will. Segments are re-build afterwards.
2) You need to define the filters in order to compute the
cross-correlations and extract the coherent part of the ambient noise
for a
given frequency range. MWCS is used to compute
the moving window cross
spectral analysis after the cross-correlation procedure.
3) dunno. May be Thomas can give you more info.
-Esteban J. Chaves
On Oct 6, 2016, at 5:00 AM,
msnoise-request(a)mailman-as.oma.be wrote:
Send MSNoise mailing list submissions to
msnoise(a)mailman-as.oma.be
To subscribe or unsubscribe via the World Wide Web, visit
http://mailman-as.oma.be/mailman/listinfo/msnoise
or, via email, send a message with subject or body 'help' to
msnoise-request(a)mailman-as.oma.be
You can reach the person managing the list at
msnoise-owner(a)mailman-as.oma.be
When replying, please edit your Subject line so it is more specific
than "Re: Contents of MSNoise digest..."
Today's Topics:
1. new user (Flinders, Ashton)
----------------------------------------------------------------------
Message: 1
Date: Wed, 5 Oct 2016 13:10:53 -0700
From: "Flinders, Ashton" <aflinders(a)usgs.gov>
To: <msnoise(a)mailman-as.oma.be>
Subject: [MSNoise] new user
Message-ID:
<CAPbWHEO7J_9dyMdKJ3SQkpGy9hG7N1M-COwgzGTpkVmy85ig5g(a)mail.gmailail.
com>
> Content-Type: text/plain; charset="UTF-8"
>
> Hi all, long-time obspy user new MSnoise user.
>
> I just had a couple quick questions I was hoping to get clarification
on;
From the SRL article;
"Once all
waveforms are loaded in memory, the computation is done by
iteration of the station pairs, on the different components to compute
and
the on the different filters for each defined
window (30 minutes slices
by
default)."
1) To be clear, the two time-series from a station-station pair are
divided
into segments (30 minutes long, no overlap, by
default), and the CC is
done
> on these segments? Are these segments then rebuilt into the daily CCF?
Is
this
segmentation purely for efficieny reasons, or is there something
else
going on?
2) Why does msnoise compute_cc, not work if there is no MCWS filter
defined?
3) Is there an intelligent way to search the mailing list archive?
Thanks for the help!
-ashton
------------------------------
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
End of MSNoise Digest, Vol 31, Issue 1
**************************************
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
--
Ashton F. Flinders, Ph.D
U.S. Geological Survey
345 Middlefield Road
Menlo Park, CA 94025
(650) 329-5050
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
--
Ashton F. Flinders, Ph.D
U.S. Geological Survey
345 Middlefield Road
Menlo Park, CA 94025
(650) 329-5050