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