[MSNoise] Error when plot dtt figure
hsy119 at mail.ustc.edu.cn
hsy119 at mail.ustc.edu.cn
Tue Sep 20 11:13:15 UTC 2016
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 at seismology.be>
> 发送时间: 2016-09-20 18:30:05 (星期二)
> 收件人: msnoise at 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 at 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 at mailman-as.oma.be
> > http://mailman-as.oma.be/mailman/listinfo/msnoise
>
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
More information about the MSNoise
mailing list