[MSNoise] Error when plot dtt figure

hsy119 at mail.ustc.edu.cn hsy119 at mail.ustc.edu.cn
Thu Sep 22 02:28:24 UTC 2016


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 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