[MSNoise] MSNoise 1.5.1

Corentin Caudron corentin.caudron at gmail.com
Fri Jan 26 18:52:57 UTC 2018


Yihaaaa

On 26 Jan 2018 7:19 PM, "Thomas Lecocq" <thomas.lecocq at seismology.be> wrote:

Dear noise fans,

I've released a bugfix (msnoise 1.5.1) today to correct:

* a nasty issue when computing autocorrelations with multiple filters
configured, resulting in a damaged spectrum (see
https://github.com/ROBelgium/MSNoise/issues/105)

* reverted the usage of next_fast_len in favor of the previous nextpox2 for
the MWCS FFT computations. "next_fast_len" is an optimised value for the
FFT library (faster computation), but usually results in a smaller number
than nextpow2. Then, as currently the smoothing of the spectrum (and
cross-spectrum) is done in number of samples (not number of samples/Hz) AND
because of the small number of points in the FFT, the resulting MWCS values
were highly noisy. The smoother provides some stability and for sure,
somehow enhances the mean coherence measured. From tests I've done, the
final results of the analysis (at the compute_dtt step) is identical, but
the noise (scatter) level was much higher. This reverts the behaviour to
what was obtained for the original version.

* a small bug that prevented the "-i" parameter to be passed to msnoise
stack.

to update:

$ pip uninstall msnoise

$ pip install msnoise

and you're done! All you need to do is then to reset the DTT jobs and rerun
compute_mwcs & compute_dtt.

----

MSNoise 1.6, the next release, will most probably appear around April 2018.
It will include :

* rewritten job logic (new job types, 1 for each step)

* rewritten compute_cc that is muuuuuuuuch faster for most jobs

* support for other spectral whitening methods (brutal "all to 1.0",
division by PSD)

* support for nomalization by the mean 3C PSD energy, to allow
post-processing analysis (particle motion, or rotation to RR and TT after
computing cross-components)

* rewritten stack, compute_mwcs and compute_dtt to make them parallelisable

* a tutorial on how to use msnoise on a 2000 core HPC using PBS


Best regards from Brussels !

Cheers

Thomas

_______________________________________________
MSNoise mailing list
MSNoise at mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise


More information about the MSNoise mailing list