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.lec...@seismology.be>
> 发送时间: 2016-09-20 18:30:05 (星期二)
> 收件人: msnoise@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, hsy...@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@mailman-as.oma.be
> > http://mailman-as.oma.be/mailman/listinfo/msnoise
> 
> _______________________________________________
> MSNoise mailing list
> MSNoise@mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
_______________________________________________
MSNoise mailing list
MSNoise@mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise

Reply via email to