can you please give a step by step instructions regarding how to input data
into the database. I have setup everything (phpmyadmin, mysql database,
msnoise web interface) . I have data from stations in ASCII format and are
in separate files of one hour time intervals. Data files have an extension
of seisan and seis_saf. Also i would like to know how to add stations and
populate station list.
Dear Thomas and msnoise users,
I'm using msnoise for the first time and so far I have been able to compute dtt without any trouble. However, when I'm trying to plot "mwcs" and "dtt" I have the following errors.
I already check the msnoise bugreport -m and it seems I have everything in order, so I really don't know what seems to be the problem.
I appreciate in advance for any help.
Kind regards,
Jose Castillo
josecastillo@josecastillo-550P5C-550P7C:~/Desktop/Project_msnoise/MS_Noise$ msnoise plot mwcs YH_DA07 YH_DB07
7.79420817017
New Data for YH_DA07_YH_DB07-ZZ-1-1
47
/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/plots/mwcs.py:146: FutureWarning: how in .resample() is deprecated
the new syntax is .resample(...).mean()
alldt = alldt.resample('D', how='mean')
/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/plots/mwcs.py:147: FutureWarning: how in .resample() is deprecated
the new syntax is .resample(...).mean()
allcoh = allcoh.resample('D', how='mean')
Traceback (most recent call last):
File "/home/josecastillo/anaconda3/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 716, in __call__
return self.main(*args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 696, in main
rv = self.invoke(ctx)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 889, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 534, in invoke
return callback(*args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 502, in mwcs
main(sta1, sta2, filterid, comp, mov_stack, show, outfile)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/plots/mwcs.py", line 159, in main
plt.imshow(alldt.T, extent=xextent,aspect="auto",interpolation='none',origin='lower',cmap=cmap_center_point_adjust(cm.seismic,crange,0))
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/plots/mwcs.py", line 80, in cmap_center_point_adjust
abs(center - range[0]) / abs(range[1] - range[0]))
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/plots/mwcs.py", line 68, in cmap_center_adjust
return cmap_powerlaw_adjust(cmap, a)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/plots/mwcs.py", line 54, in cmap_powerlaw_adjust
cdict[key].sort()
AttributeError: 'map' object has no attribute 'sort'
josecastillo@josecastillo-550P5C-550P7C:~/Desktop/Project_msnoise/MS_Noise$ msnoise plot dtt YH_DA07 YH_DB07 2012-10-11
MWCS/01/001_DAYS/ZZ/YH_DA07_YH_DB07/2012-10-11.txt
Empty DataFrame
Columns: [Date, A, EA, EM, EM0, M, M0, Pairs]
Index: []
Traceback (most recent call last):
File "/home/josecastillo/anaconda3/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 716, in __call__
return self.main(*args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 696, in main
rv = self.invoke(ctx)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 889, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/core.py", line 534, in invoke
return callback(*args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 559, in dtt
main(sta1, sta2, filterid, comp, day, mov_stack, show, outfile)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/msnoise/plots/dtt.py", line 91, in main
plt.plot(xlineM0, ylineM0, 'r', label='M0')
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/matplotlib/pyplot.py", line 3154, in plot
ret = ax.plot(*args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/matplotlib/__init__.py", line 1812, in inner
return func(ax, *args, **kwargs)
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_axes.py", line 1424, in plot
for line in self._get_lines(*args, **kwargs):
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 386, in _grab_next_args
for seg in self._plot_args(remaining, kwargs):
File "/home/josecastillo/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 374, in _plot_args
seg = func(x[:, j % ncx], y[:, j % ncy], kw, kwargs)
ZeroDivisionError: integer division or modulo by zero
Let's Bug Report MSNoise !
************* Computer Report *************
----------------+PYTHON+-------------------
Python:3.5.1 |Anaconda 4.0.0 (64-bit)| (default, Dec 7 2015, 11:16:01)
[GCC 4.4.7 20120313 (Red Hat 4.4.7-1)]
---------------+MODULES+-------------------
Required:
[X] setuptools: 20.3
[X] click: 6.6
[X] numpy: 1.10.4
[X] scipy: 0.17.0
[X] pandas: 0.18.0
[X] matplotlib: 1.5.1
[X] statsmodels: 0.6.1
[X] sqlalchemy: 1.0.12
[X] traitsui: 5.1.0
[X] traits: 4.5.0
[X] scikits.samplerate: present (no version)
[X] obspy: 1.0.1
[X] flask: 0.10.1
[X] flask.ext.admin: 0.10.1
[X] flask_admin: 1.4.2
[X] wtforms: 2.1
[X] bokeh: 0.11.1
Only necessary if you plan to build the doc locally:
[X] sphinx: 1.3.5
[X] jinja2: 2.8
Graphical Backends: (at least one is required)
[ ] wx: not found (install via None)
[X] PyQt4: present (no version)
[ ] PySide: not found (install via None)
Not required, just checking:
[X] json: 2.0.9
[X] psutil: 4.1.0
[ ] reportlab: not found (install via None)
[X] configobj: 5.0.6
[X] pkg_resources: present (no version)
[ ] paramiko: not found (install via None)
[X] ctypes: 1.1.0
[X] pyparsing: 2.0.3
[X] distutils: 3.5.1
[X] IPython: 4.1.2
[ ] vtk: not found (install via None)
[ ] enable: not found (install via None)
Hi all,
This is an annoucement for a small bugfix release: MSNoise 1.4.1
contains corrections for:
* Typos
* Hardcode version to prevent pypi errors
* dynamic loading of CDN or local jquery
* fixed raw_input vs input (Py2 vs Py3) at "installation step"
Nothing 'core' business, so these "bugs" will not have affected your
results using MSNoise!
By the way: here are two recent really cool papers using MSNoise !
* Mordret, Aurélien, T. Dylan Mikesell, Christopher Harig, Bradley P.
Lipovsky, et Germán A. Prieto. « Monitoring Southwest Greenland’s Ice
Sheet Melt with Ambient Seismic Noise ». /Science Advances/ 2, n^o 5 (1
mai 2016): e1501538. doi:10.1126/sciadv.1501538.
* Wu, Chunquan, Andrew Delorey, Florent Brenguier, Celine Hadziioannou,
Eric G Daub, et Paul Johnson. « Constraining Depth Range of S-Wave
Velocity Decrease after Large Earthquakes near Parkfield, California. »
/Geophysical Research Letters/, 1 mai 2016, 2016GL069145.
doi:10.1002/2016GL069145.
Cheers !
Thomas
--
Dr. Thomas Lecocq
Geologist - Seismologist
Seismology - Gravimetry
Royal Observatory of Belgium
*
* * * * *
* * * *
---------
http://www.seismology.behttp://twitter.com/#!/Seismologie_behttps://www.facebook.com/seismologie.be
Hi Thomas,
As far as I understand you haven't changed anything since then in the code
concerning the pre-processing when computing autocorrelations. Right?
Greetings from Athens,
Christos
From: Thomas Lecocq <thomas.lecocq@...
<http://gmane.org/get-address.php?address=thomas.lecocq%2d7ewb5K7nbqnl%2fIpH…>
>
Subject: Re: Autocorrelations in MSNoise
<http://news.gmane.org/find-root.php?message_id=54DB931F.5080705%40seismolog…>
Newsgroups: gmane.science.geophysics.msnoise
<http://news.gmane.org/gmane.science.geophysics.msnoise>
Date: 2015-02-11 17:36:31 GMT (1 year, 14 weeks, 6 days, 11 hours and 16
minutes ago)
Hello Aurélien,
Le 11/02/2015 18:30, Aurélien Mordret a écrit :
> Hi Thomas,
>
> What is the pre-processing of the noise when computing autocorrelations? Is
> it the same as for the cross-correlations?
Yeah, currently it's the same, WHICH IS NOT GOOD ! (that's why it's
written somewhere NOT TO USE AUTOCORR), but that might not be trivial to
find...
Your suggestion makes absolutely sense ! Are you used to github etc ? If
yes, please fork and PR an updated version of whiten ! (including a
autocorr=False flag in the def, maybe?)!
Regards from Brussels
Thomas
>
> If yes, would it be better to not fully whiten the data before
> autocorrelation as the only information comes from the amplitude spectrum
> (and not the phase)?
>
> I suggest to modify a bit the whiten.py function to keep the structure of
> the whole code similar. Following Gorbatov et al. (2013), instead of
> putting the amplitude spectrum =1 (with tapper) I would use a normalisation
> with a water-level in the form
>
> A(freq) = A(freq)/max(A(freq),c*max(A(freq)))
>
> With A the amplitude spectrum and c the water-level factor (which is
> typically on the order of 1%).
>
> It should not change the results for the cross-correlations and improve the
> autocorrelations.
>
> Does it make sense or did I missed something?
>
> Cheers,
>
> Aurelien
>
> REF:
>
> Gorbatov, A., Saygin, E., & Kennett, B. L. N. (2013). Crustal properties
> from seismic station autocorrelograms. *Geophysical Journal International*,
> **192**(2), 861-870.
>
