Hi,
Thanks,
I fix it, it works.
SAC is nice because of the metadata coming with the data. For a correlation
file, the possibility to have the event coordinates as well as the station
coordinates in the header is very convenient because you can fill the event
coordinates with STA1 coordinates and station coordinates with STA2
coordinates. So distance, azimuth and back-azimuth are easily (and
automatically) computed.
And I developed a code to automatically compute group velocity dispersion
curves from CC, which reads the distance between the stations in the
header, without the need of other metadata file. ;)
Cheers,
Aurelien
2015-03-31 10:15 GMT-04:00 Thomas Lecocq <thomas.lecocq(a)seismology.be>be>:
Hi Aurélien !
Le 31/03/2015 16:06, Aurélien Mordret a écrit :
Hi Thomas,
Great job for the 1.3.1 MSNoise version!
Thanks!
I have however encountered a problem when computing the MWCS: the function
only takes Miniseed files as input, which is no compatible with the fact
that you can save the CC as SAC files.
Aha ! Indeed... Why would one still work with SAC anyway :-)
Until I provide a bugfix, hack the code at two places:
https://github.com/ROBelgium/MSNoise/blob/master/msnoise/
s05compute_mwcs.py#L109
and
https://github.com/ROBelgium/MSNoise/blob/master/msnoise/
s05compute_mwcs.py#L115
Let me know,
Tom
Here is the error below.
Cheers,
Aurelien
Traceback (most recent call last):
File "/usr/local/bin/msnoise", line 9, in <module>
load_entry_point('msnoise==1.3.1', 'console_scripts',
'msnoise')()
File "/usr/local/lib/python2.7/dist-packages/msnoise/scripts/
msnoise.py",
line 393, in run
cli(obj={})
File "/usr/local/lib/python2.7/dist-packages/click/core.py", line
610, in
__call__
return self.main(*args, **kwargs)
File "/usr/local/lib/python2.7/dist-packages/click/core.py", line
590, in
main
rv = self.invoke(ctx)
File "/usr/local/lib/python2.7/dist-packages/click/core.py", line
936, in
invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/usr/local/lib/python2.7/dist-packages/click/core.py", line
782, in
invoke
return ctx.invoke(self.callback, **ctx.params)
File "/usr/local/lib/python2.7/dist-packages/click/core.py", line
416, in
invoke
return callback(*args, **kwargs)
File "/usr/local/lib/python2.7/dist-packages/msnoise/scripts/
msnoise.py",
line 199, in compute_mwcs
main()
File "/usr/local/lib/python2.7/dist-packages/msnoise/
s05compute_mwcs.py",
line 110, in main
ref = read(rf)[0].data
File
"/usr/local/lib/python2.7/dist-packages/obspy/core/util/decorator.py",
line
371, in new_func
return func(*args, **kwargs)
File "/usr/local/lib/python2.7/dist-packages/obspy/core/stream.py",
line
224, in read
raise IOError(2, "No such file or directory", pathname)
IOError: [Errno 2] No such file or directory:
'STACKS/01/REF/ZZ/DK_ANGG_DK_ANGG.MSEED'
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
--
==========================
Postdoctoral fellow
Department of Earth, Atmospheric and Planetary Sciences
Massachusetts Institute of Technology
54-526
Cambridge, MA 02139-4307
e-mail: mordret(a)mit.edu
web: