Dear community,
About 1 year after the last major release (MSNoise 1.4
<http://msnoise.org/doc/releasenotes/msnoise-1.4.html>) we are proud to
announce the new MSNoise 1.5
<http://msnoise.org/doc/releasenotes/msnoise-1.5.html#>. It is a *major*
release, with a massive amount of work since the last one: in GitHub
numbers
<https://github.com/ROBelgium/MSNoise/graphs/contributors?from=2016-06-02&to…>
, it’s over 120 commits and over 2500 lines of code and documentation
changed or added!
MSNoise 1.5 introduces a series of *new features* :
* We have started to move core math functions to ObsPy, currently the
only one ready is linear_regression, a function I wrote to remove
the dependency to |statsmodels|, required to move mwcs to ObsPy later.
* The preprocessing routine has been isolated, rewritten and
optimized. It is now a standalone script, callable by plugins. It
returns a Stream object with all the data needed for the analysis.
* This change in preprocessing was done mostly to allow
cross-component, auto- correlation and cross-correlation, with or
without rotation, to be done with the same code. CC, SC and AC are
now supported in MSNoise with proper whitening (possible to disable
spectral whitening for specific cases).
* This documentation is now available in PDF
<http://msnoise.org/doc/MSNoise.pdf> too (easier for offline usage)
and it also includes a new tutorial for setting up the MySQL server
and Workbench.
* Last but not least: MSNoise is “tested” automatically on Linux
(thanks to TravisCI) & Windows (thanks to Appveyor), for Python
versions 2.7 and 3.5. With MSNoise 1.5 we also added the MacOSX
tests on TravisCI. With these tests, we can guarantee MSNoise works
on different platforms and Anaconda (or miniconda) python versions.
This version has benefited from outputs/ideas/pull requests/questions
from several users/friends (listed alphabetically):
* Raphael De Plaen
* Clare Donaldson
* Robert Green
* Aurelien Mordret
* Lukas Preiswerk
* The participants to the NERC MSNoise Liverpool Workshop in January 2017
* all others (don’t be mad :-) )
Thanks to all for using MSNoise, and please, let us know why/how you use
it (and please cite it!)!
To date, we found/are aware of 25 publications using MSNoise ! That’s
the best validation of our project ever ! See the full list on the
MSNoise website <http://www.msnoise.org/they-cite-msnoise/>.
/Thomas, Corentin and others/
------------------------------------------------------------------------
PS: if you use MSNoise for your research and prepare publications,
*please consider citing it*:
*Lecocq, T., C. Caudron, et F. Brenguier (2014)*, MSNoise, a Python
Package for Monitoring Seismic Velocity Changes Using Ambient Seismic
Noise, /Seismological Research Letters/, 85(3), 715‑726,
doi:10.1785/0220130073.
--
Dr. Thomas Lecocq
Geologist - Seismologist
Seismology - Gravimetry
Royal Observatory of Belgium
*
* * * * *
* * * *
---------
http://www.seismology.behttp://msnoise.orghttp://twitter.com/#!/Seismologie_behttps://www.facebook.com/seismologie.be
Hi,Thomas
i have already installed msnosie 1.4.how can i solve the autocorrlation problem? many thanks!
> -----Original Messages-----
> From: msnoise-request(a)mailman-as.oma.be
> Sent Time: Saturday, April 15, 2017
> To: msnoise(a)mailman-as.oma.be
> Cc:
> Subject: MSNoise Digest, Vol 36, Issue 9
>
> 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. Re: autocorrelation resulting problem (Thomas Lecocq)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Sat, 15 Apr 2017 11:33:29 +0200
> From: Thomas Lecocq <Thomas.Lecocq(a)seismology.be>
> To: msnoise(a)mailman-as.oma.be
> Subject: Re: [MSNoise] autocorrelation resulting problem
> Message-ID: <7cd0b509-461c-37a9-4596-16207a0474fb(a)seismology.be>
> Content-Type: text/plain; charset=UTF-8; format=flowed
>
> Hi,
>
> You are using the master version or release 1.4?
>
> Thomas
>
>
> On 15/04/2017 03:46, ?? wrote:
> > Hi,Thomas
> > I'm Sun Jie . I already I 'config' the 'autocorr=Y','stack_method=linear','mov_stack=10','keep_days=Y',and etc. And there in no error in the working process .But in the 'ccftime plot',the waveforms of every '1day-stack'is exactly same,the '10days-stack'is so on.I cannot find the reason why the autocorrelation waveforms is the same.How can i solve this problem? Many thanks.
> >
> >
> >
> > _______________________________________________
> > 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
>
>
> End of MSNoise Digest, Vol 36, Issue 9
> **************************************
Hi,Thomas
I'm Sun Jie . I already I 'config' the 'autocorr=Y','stack_method=linear','mov_stack=10','keep_days=Y',and etc. And there in no error in the working process .But in the 'ccftime plot',the waveforms of every '1day-stack'is exactly same,the '10days-stack'is so on.I cannot find the reason why the autocorrelation waveforms is the same.How can i solve this problem? Many thanks.
Dear Thomas and users,
It is a while that i noticed that the final stacking or mov windows don't include in the header the distance info and a no standard Beginning/ending value of the trace. I mean if i have maxlag =120.0 it means that when i plot the race i should have B=-120.0 and E=120.0 so to have the 0 in the middle. I wrote two script to solve this ( it is not an issue but have it done should be better). i attached the two script and i hope that in some way they can be implemented as python code or implemented as plugin. It will be very cool! I have not time to do it but my c-shell script can be a start point.
I remember that in an old version of msnoise i changed the python code to include the distance in the header file...and i give it personally to Thomas.. so if possible may you include it in the next msnoise version? if not possible because i understand that this version is a lot different than the msnoise beginning (when still it was not a python package) please if you are able to realise some plugin just let me know.. it will be very nice do this operation in automatic also using plugin. [😊]
Have a lovely day everyone.
cheers
-----------
Mr Carmelo Sammarco MSc BSc FGS
PhD Student in Geology and Petroleum Geology
The University of Aberdeen is a charity registered in Scotland, No SC013683.
Tha Oilthigh Obar Dheathain na charthannas clàraichte ann an Alba, Àir. SC013683.
Hi Thomas
I am C.J Liao, and i'm installing the Msnoise 1.4 version in Ubuntu system with Quick, i have a question when i run 'msnoise stack -r -m ' I 'config' the stack_method=pws ,pws power=2 , pws_timegate=10.0,autocorr=N,RR=Y,ZZ=Y....,
There is following error when I 'stack'. How can i solve this problem? many thanks!
$ msnoise stack -r -m
rusult in
Lets STACK !
2017-04-10 16:16:00 [DEBUG] Starting the ref stack
[]
2017-04-10 16:16:00 [DEBUG] Processing YP_NE12:YP_NE48-ZZ-1
2017-04-10 16:16:00 [DEBUG] Found 295 updated days
2017-04-10 16:16:00 [DEBUG] New Data for YP_NE12:YP_NE48-ZZ-1
/home/l/anaconda2/lib/python2.7/site-packages/obspy/io/sac/util.py:370: UserWarning: Old header has invalid reftime.
warnings.warn(msg)
2017-04-10 16:16:04 [DEBUG] Processing YP_NE12:YP_NE48-RR-1
2017-04-10 16:16:04 [DEBUG] Found 295 updated days
2017-04-10 16:16:04 [DEBUG] New Data for YP_NE12:YP_NE48-RR-1
Traceback (most recent call last):
File "/home/l/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/l/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/l/.local/lib/python2.7/site-packages/click/core.py", line 722, in __call__
return self.main(*args, **kwargs)
File "/home/l/.local/lib/python2.7/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/home/l/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/l/.local/lib/python2.7/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/l/.local/lib/python2.7/site-packages/click/core.py", line 535, in invoke
return callback(*args, **kwargs)
File "/home/l/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 301, in stack
main('ref', interval)
File "/home/l/anaconda2/lib/python2.7/site-packages/msnoise/s04stack.py", line 209, in main
db, sta1, sta2, filterid, components, datelist, format=format)
File "/home/l/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line 1020, in get_results
corr = stack(session, stack_data)
File "/home/l/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line 978, in stack
coh = 1. / data.shape[0] * np.abs(phasestack)
ZeroDivisionError: float division by zero
Hello,
I am a new using.
I am at first step. Please let me know how can I pupplate my satation
information and link to data directory.
Thanks
--
Regards,
Yawar Hussain
Ph.D. Student, Geotechnical Engineering
University of Brasilia, Brazil.
Hi, Thomas
I'm Sun Jie and i am a python beginner.Following your suggestion, I have solved the 'pws compute_cc' problem ,many thanks!
However, i have another question when i run 'msnoise stack -r -m ' I 'config' the stack_method=pws ,pws power=2 , pws_timegate=10.0,autocorr=Y,ZZ=Y....,
There is following error when I 'stack'. How can i solve this problem? many thanks!
$ msnoise stack -r -m
rusult in
Lets STACK !
2017-03-30 16:46:58 [DEBUG] Starting the ref stack
[]
2017-03-30 16:46:58 [DEBUG] Processing YP_NE49:YP_NE49-ZZ-1
2017-03-30 16:46:58 [DEBUG] Found 278 updated days
2017-03-30 16:46:58 [DEBUG] New Data for YP_NE49:YP_NE49-ZZ-1
/home/sj/anaconda2/lib/python2.7/site-packages/obspy/io/sac/util.py:370: UserWarning: Old header has invalid reftime.
warnings.warn(msg)
2017-03-30 16:47:01 [DEBUG] Processing YP_NE49:YP_NE49-RR-1
2017-03-30 16:47:01 [DEBUG] Found 278 updated days
2017-03-30 16:47:01 [DEBUG] New Data for YP_NE49:YP_NE49-RR-1
Traceback (most recent call last):
File "/home/sj/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 722, in __call__
return self.main(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 535, in invoke
return callback(*args, **kwargs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 301, in stack
main('ref', interval)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/s04stack.py", line 209, in main
db, sta1, sta2, filterid, components, datelist, format=format)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line 1020, in get_results
corr = stack(session, stack_data)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line 978, in stack
coh = 1. / data.shape[0] * np.abs(phasestack)
ZeroDivisionError: float division by zero
