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