Hi, Thomas
I'm not sure if you have got my reply email,so I write again. If I ask too much, sorry for that......Please let me know...
I use "msnoise plot dtt net.station1 net.station2 2015-01-01 -f 01 -m 10" to plot figure, but there's an error:
Traceback (most recent call last):
File "/public/software/anaconda3/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 716, in __call__
return self.main(*args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 696, in main
rv = self.invoke(ctx)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 889, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 534, in invoke
return callback(*args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/plots/dtt.py", line 91, in main
plt.plot(xlineM0, ylineM0, 'r', label='M0')
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/pyplot.py", line 3154, in plot
ret = ax.plot(*args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/__init__.py", line 1812, in inner
return func(ax, *args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_axes.py", line 1424, in plot
for line in self._get_lines(*args, **kwargs):
File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 360, in _plot_args
y = _check_1d(tup[-1])
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/cbook.py", line 2207, in _check_1d
return np.atleast_1d(x)
File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/shape_base.py", line 50, in atleast_1d
ary = asanyarray(ary)
File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/numeric.py", line 533, in asanyarray
return array(a, dtype, copy=False, order=order, subok=True)
File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/core/series.py", line 583, in __getitem__
result = self.index.get_value(self, key)
File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/indexes/base.py", line 1980, in get_value
tz=getattr(series.dtype, 'tz', None))
File "pandas/index.pyx", line 103, in pandas.index.IndexEngine.get_value (pandas/index.c:3332)
File "pandas/index.pyx", line 111, in pandas.index.IndexEngine.get_value (pandas/index.c:3035)
File "pandas/index.pyx", line 159, in pandas.index.IndexEngine.get_loc (pandas/index.c:4018)
File "pandas/hashtable.pyx", line 303, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6610)
File "pandas/hashtable.pyx", line 309, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6554)
KeyError: 0
It's strange that when I want to plot dtt figure of THE station pair(that is G2.1592_G2.1647),it works. But for the rest pairs, I got an error.
No error when I compute ccf,mwcs,dtt.
So what should I do?
Thank you for help.
(I misspell "rest" in my first email, very sorry for that)
Shuye
> ----------
> Hi Shuye,
>
> I don't understand what you mean with "the reset pairs" ? What is the
> command you used to start this plot ?
>
> Thomas
>
>
> On 20/09/2016 03:10, hsy119(a)mail.ustc.edu.cn wrote:
> > Hi,
> >
> > Sorry to bother again.
> > I'm trying to plot dtt figure but I meet an error. It's strange because there's no problem when I plot a certain station pair, whatever the data is. But for the reset pairs, there's a bug as follows.
> >
> >
> > Date A EA EM EM0 M M0 \
> > 1 2015-09-01 0.010862 0.010306 0.000708 0.000445 -0.000198 0.000383
> >
> > Pairs
> > 1 G2_1592_G2_1650
> > Traceback (most recent call last):
> > File "/public/software/anaconda3/bin/msnoise", line 11, in <module>
> > sys.exit(run())
> > File "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 614, in run
> > cli(obj={})
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 716, in __call__
> > return self.main(*args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 696, in main
> > rv = self.invoke(ctx)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
> > return _process_result(sub_ctx.command.invoke(sub_ctx))
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
> > return _process_result(sub_ctx.command.invoke(sub_ctx))
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 889, in invoke
> > return ctx.invoke(self.callback, **ctx.params)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 534, in invoke
> > return callback(*args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/decorators.py", line 17, in new_func
> > return f(get_current_context(), *args, **kwargs)
> > File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/plots/dtt.py", line 91, in main
> > plt.plot(xlineM0, ylineM0, 'r', label='M0')
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/pyplot.py", line 3154, in plot
> > ret = ax.plot(*args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/__init__.py", line 1812, in inner
> > return func(ax, *args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_axes.py", line 1424, in plot
> > for line in self._get_lines(*args, **kwargs):
> > File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 360, in _plot_args
> > y = _check_1d(tup[-1])
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/cbook.py", line 2207, in _check_1d
> > return np.atleast_1d(x)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/shape_base.py", line 50, in atleast_1d
> > ary = asanyarray(ary)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/numeric.py", line 533, in asanyarray
> > return array(a, dtype, copy=False, order=order, subok=True)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/core/series.py", line 583, in __getitem__
> > result = self.index.get_value(self, key)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/indexes/base.py", line 1980, in get_value
> > tz=getattr(series.dtype, 'tz', None))
> > File "pandas/index.pyx", line 103, in pandas.index.IndexEngine.get_value (pandas/index.c:3332)
> > File "pandas/index.pyx", line 111, in pandas.index.IndexEngine.get_value (pandas/index.c:3035)
> > File "pandas/index.pyx", line 159, in pandas.index.IndexEngine.get_loc (pandas/index.c:4018)
> > File "pandas/hashtable.pyx", line 303, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6610)
> > File "pandas/hashtable.pyx", line 309, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6554)
> >
> > There's no error when I compute cc, mwcs, dtt,only some warning.
> > So how can I fix this?
> > Thanks for your help
> >
> > Shuye
> >
> > _______________________________________________
> > 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
Dear all
I would like to try to use msnoise. I installed it on Ubuntu 16.04, and as far as I know, everything was installed as planed.
Now, I open msnoise admin and try to manually add stations in the browser. I set NET, STA, X, Y, ALT, COORD and click save. When I do that, I get a lot of errors:
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 2000, in __call__
return self.wsgi_app(environ, start_response)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 1991, in wsgi_app
response = self.make_response(self.handle_exception(e))
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 1567, in handle_exception
reraise(exc_type, exc_value, tb)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 1988, in wsgi_app
response = self.full_dispatch_request()
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 1641, in full_dispatch_request
rv = self.handle_user_exception(e)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 1544, in handle_user_exception
reraise(exc_type, exc_value, tb)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 1639, in full_dispatch_request
rv = self.dispatch_request()
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask/app.py", line 1625, in dispatch_request
return self.view_functions[rule.endpoint](**req.view_args)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask_admin/base.py", line 69, in inner
return self._run_view(f, *args, **kwargs)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask_admin/base.py", line 368, in _run_view
return fn(self, *args, **kwargs)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask_admin/model/base.py", line 1920, in create_view
model = self.create_model(form)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask_admin/contrib/sqla/view.py", line 1028, in create_model
if not self.handle_view_exception(ex):
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask_admin/contrib/sqla/view.py", line 1011, in handle_view_exception
return super(ModelView, self).handle_view_exception(exc)
*
File "/home/blaz/miniconda2/lib/python2.7/site-packages/flask_admin/contrib/sqla/view.py", line 1022, in create_model
model = self.model()
TypeError: __init__() takes exactly 9 arguments (1 given)
What am I doing wrong?
Thanks,
Blaz
Hi,Thomas
Well, I use "msnoise plot dtt net.station1 net.station2 2015-01-01 -f 01 -m 10" to plot.
I mean, only when I want to plot a certain station pair,the command works(actually the pair is G2.1592_G2.1647, G2 is net)
But when I try to plot other station pairs, for example, G2.1592_G2.1650, there's an error. Just like this.
Traceback (most recent call last):
File "/public/software/anaconda3/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 716, in __call__
return self.main(*args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 696, in main
rv = self.invoke(ctx)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 889, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 534, in invoke
return callback(*args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/plots/dtt.py", line 91, in main
plt.plot(xlineM0, ylineM0, 'r', label='M0')
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/pyplot.py", line 3154, in plot
ret = ax.plot(*args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/__init__.py", line 1812, in inner
return func(ax, *args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_axes.py", line 1424, in plot
for line in self._get_lines(*args, **kwargs):
File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 360, in _plot_args
y = _check_1d(tup[-1])
File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/cbook.py", line 2207, in _check_1d
return np.atleast_1d(x)
File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/shape_base.py", line 50, in atleast_1d
ary = asanyarray(ary)
File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/numeric.py", line 533, in asanyarray
return array(a, dtype, copy=False, order=order, subok=True)
File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/core/series.py", line 583, in __getitem__
result = self.index.get_value(self, key)
File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/indexes/base.py", line 1980, in get_value
tz=getattr(series.dtype, 'tz', None))
File "pandas/index.pyx", line 103, in pandas.index.IndexEngine.get_value (pandas/index.c:3332)
File "pandas/index.pyx", line 111, in pandas.index.IndexEngine.get_value (pandas/index.c:3035)
File "pandas/index.pyx", line 159, in pandas.index.IndexEngine.get_loc (pandas/index.c:4018)
File "pandas/hashtable.pyx", line 303, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6610)
File "pandas/hashtable.pyx", line 309, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6554)
KeyError: 0
There's 78 pairs in my net, only 1 pair can plot dtt figure. I've tried 20 pairs, it doesn't work, probably the same for the rest ones.
Sorry for my poor expression in English...(I want to say "rest", but "reset" instead)
Shuye
> -----原始邮件-----
> 发件人: "Thomas Lecocq" <Thomas.Lecocq(a)seismology.be>
> 发送时间: 2016-09-20 18:30:05 (星期二)
> 收件人: msnoise(a)mailman-as.oma.be
> 抄送:
> 主题: Re: [MSNoise] Error when plot dtt figure
>
> Hi Shuye,
>
> I don't understand what you mean with "the reset pairs" ? What is the
> command you used to start this plot ?
>
> Thomas
>
>
> On 20/09/2016 03:10, hsy119(a)mail.ustc.edu.cn wrote:
> > Hi,
> >
> > Sorry to bother again.
> > I'm trying to plot dtt figure but I meet an error. It's strange because there's no problem when I plot a certain station pair, whatever the data is. But for the reset pairs, there's a bug as follows.
> >
> >
> > Date A EA EM EM0 M M0 \
> > 1 2015-09-01 0.010862 0.010306 0.000708 0.000445 -0.000198 0.000383
> >
> > Pairs
> > 1 G2_1592_G2_1650
> > Traceback (most recent call last):
> > File "/public/software/anaconda3/bin/msnoise", line 11, in <module>
> > sys.exit(run())
> > File "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py", line 614, in run
> > cli(obj={})
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 716, in __call__
> > return self.main(*args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 696, in main
> > rv = self.invoke(ctx)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
> > return _process_result(sub_ctx.command.invoke(sub_ctx))
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 1060, in invoke
> > return _process_result(sub_ctx.command.invoke(sub_ctx))
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 889, in invoke
> > return ctx.invoke(self.callback, **ctx.params)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py", line 534, in invoke
> > return callback(*args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/click/decorators.py", line 17, in new_func
> > return f(get_current_context(), *args, **kwargs)
> > File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/plots/dtt.py", line 91, in main
> > plt.plot(xlineM0, ylineM0, 'r', label='M0')
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/pyplot.py", line 3154, in plot
> > ret = ax.plot(*args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/__init__.py", line 1812, in inner
> > return func(ax, *args, **kwargs)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_axes.py", line 1424, in plot
> > for line in self._get_lines(*args, **kwargs):
> > File "/public/software/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 "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 360, in _plot_args
> > y = _check_1d(tup[-1])
> > File "/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/cbook.py", line 2207, in _check_1d
> > return np.atleast_1d(x)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/shape_base.py", line 50, in atleast_1d
> > ary = asanyarray(ary)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/numeric.py", line 533, in asanyarray
> > return array(a, dtype, copy=False, order=order, subok=True)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/core/series.py", line 583, in __getitem__
> > result = self.index.get_value(self, key)
> > File "/public/software/anaconda3/lib/python3.5/site-packages/pandas/indexes/base.py", line 1980, in get_value
> > tz=getattr(series.dtype, 'tz', None))
> > File "pandas/index.pyx", line 103, in pandas.index.IndexEngine.get_value (pandas/index.c:3332)
> > File "pandas/index.pyx", line 111, in pandas.index.IndexEngine.get_value (pandas/index.c:3035)
> > File "pandas/index.pyx", line 159, in pandas.index.IndexEngine.get_loc (pandas/index.c:4018)
> > File "pandas/hashtable.pyx", line 303, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6610)
> > File "pandas/hashtable.pyx", line 309, in pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6554)
> >
> > There's no error when I compute cc, mwcs, dtt,only some warning.
> > So how can I fix this?
> > Thanks for your help
> >
> > Shuye
> >
> > _______________________________________________
> > 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
Hello,
I periodically receive this error when running (very large) cross correlation jobs:
“pymysql.err.OperationalError 2006, "MySQL server has gone away (error(32, 'Broken pipe'))")”
I have modified wait_timeout to 24 hours and increased the max_allowed_packet variable in the mySQL server (the two suggested sources of the error I’ve found through Google,) but this error is still occurring.
Any recommendations?
Thanks,
Chrysm
Hello,
I am trying to look for reference for msnoise, but when I opened
http://msnoise.org/doc it shows "Page not found". Is there a problem with
the msnoise website?
Warm regards,
Mochammad Husni Rizal
Technician of geophysics equipments for DOMERAPI -- Research Engineer
Phone : (+62) 857-151-220-90
Hi Shuye,
there are some possible causes of this:
1/ the start & end dates are not properly defined --> edit config
2/ the ref start & end dates are not properly defined --> edit config
3/ the CC jobs have been done more than 1 day ago --> run "msnoise stack
-r -m -i 9999 " to include all computations,
Let us know how it goes,
Thomas
Le 06/09/2016 à 12:08, hsy119 a écrit :
> Hi,
>
> I've do the cross-correlation jobs, and now stacking, but there's a
> problem. I have data in my data folder from 2015-07-01 to 2015-12-01
> of 3 stations, and 100+ cross-correlation functions for each pair have
> been computed(stored in 001_DAYS I guess)
> But when I stacking, only a few days' stacking have been done. The
> folder 010_DAYS/station1_station2/ only have tens of .mseed files.
> What has happend?
>
> Everything is right when I only populate 2 stations. Stacking is all
> right and 100+ files is in 010_DAYS/station1_station2/.
>
> MSNoise is very nice tool, even for the beginner of python like me.
> But I'm stucking in this step.
> I'm appreciated for your help
>
> Shuye
>
