Hi, Thomas
I'm Sun Jie and i am a python beginner.Following your suggestion, I have solved the 'pws compute_cc' problem ,many thanks!
However, i have another question when i run 'msnoise stack -r -m ' I 'config' the stack_method=pws ,pws power=2 , pws_timegate=10.0,autocorr=Y,ZZ=Y....,
There is following error when I 'stack'. How can i solve this problem? many thanks!
$ msnoise stack -r -m
rusult in
Lets STACK !
2017-03-30 16:46:58 [DEBUG] Starting the ref stack
[]
2017-03-30 16:46:58 [DEBUG] Processing YP_NE49:YP_NE49-ZZ-1
2017-03-30 16:46:58 [DEBUG] Found 278 updated days
2017-03-30 16:46:58 [DEBUG] New Data for YP_NE49:YP_NE49-ZZ-1
/home/sj/anaconda2/lib/python2.7/site-packages/obspy/io/sac/util.py:370: UserWarning: Old header has invalid reftime.
warnings.warn(msg)
2017-03-30 16:47:01 [DEBUG] Processing YP_NE49:YP_NE49-RR-1
2017-03-30 16:47:01 [DEBUG] Found 278 updated days
2017-03-30 16:47:01 [DEBUG] New Data for YP_NE49:YP_NE49-RR-1
Traceback (most recent call last):
File "/home/sj/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 722, in __call__
return self.main(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 535, in invoke
return callback(*args, **kwargs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 301, in stack
main('ref', interval)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/s04stack.py", line 209, in main
db, sta1, sta2, filterid, components, datelist, format=format)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line 1020, in get_results
corr = stack(session, stack_data)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line 978, in stack
coh = 1. / data.shape[0] * np.abs(phasestack)
ZeroDivisionError: float division by zero
Hi all
I am trying to set up MSNoise in my new windows PC(windows 10 Pro) and
have a problem when I attempt to install scikits.samplerate by a command
below:
pip install scikits.samplerate-0.3.3-cp27-none-win_amd64.whl
If I type this command, I had some error message(below) and could not
install scikits.samplerate appropriately. Does anyone have the same
problem before? I appreciate if you let me know how to avoid this error.
I am sorry for dammy question.
Thanks in advance,
Run
##############################################
Processing
c:\windows\system32\scikits.samplerate-0.3.3-cp27-none-win_amd64.whl
Exception:
Traceback (most recent call last):
File "C:\ProgramData\Anaconda2\lib\site-packages\pip\basecommand.py",
line 215, in main
status = self.run(options, args)
File
"C:\ProgramData\Anaconda2\lib\site-packages\pip\commands\install.py",
line 335, in run
wb.build(autobuilding=True)
File "C:\ProgramData\Anaconda2\lib\site-packages\pip\wheel.py", line
749, in build
self.requirement_set.prepare_files(self.finder)
File "C:\ProgramData\Anaconda2\lib\site-packages\pip\req\req_set.py",
line 380, in prepare_files
ignore_dependencies=self.ignore_dependencies))
File "C:\ProgramData\Anaconda2\lib\site-packages\pip\req\req_set.py",
line 620, in _prepare_file
session=self.session, hashes=hashes)
File "C:\ProgramData\Anaconda2\lib\site-packages\pip\download.py",
line 809, in unpack_url
unpack_file_url(link, location, download_dir, hashes=hashes)
File "C:\ProgramData\Anaconda2\lib\site-packages\pip\download.py",
line 715, in unpack_file_url
unpack_file(from_path, location, content_type, link)
File
"C:\ProgramData\Anaconda2\lib\site-packages\pip\utils\__init__.py", line
599, in unpack_file
flatten=not filename.endswith('.whl')
File
"C:\ProgramData\Anaconda2\lib\site-packages\pip\utils\__init__.py", line
482, in unzip_file
zipfp = open(filename, 'rb')
IOError: [Errno 2] No such file or directory:
'C:\\Windows\\system32\\scikits.samplerate-0.3.3-cp27-none-win_amd64.whl'
Hi, Thomas
I'm Sun Jie , very sorry to sent you email twice due to some problems in my first email system.
Following your suggestion, I have solved the autocorrelation plotting problem ,many thanks!
However, i have another question about pws stacking. I 'config' the stack_method=pws ,pws power=2 , pws_timegate=10.0,autocorr=Y,ZZ=Y....,
There is following error when I 'compute_cc'. How can i solve this problem? many thanks!
$ msnoise compute_cc
result in
2017-03-30 09:38:15 [INFO] *** Starting: Compute CC ***
2017-03-30 09:38:15 [INFO] Will compute ZZ RR TT
2017-03-30 09:38:16 [INFO] New CC Job: 2010-03-11 (1 pairs with 1 stations)
2017-03-30 09:38:16 [DEBUG] YP.NE49.Z Reading 1 Files
2017-03-30 09:38:16 [DEBUG] Checking sample alignment
2017-03-30 09:38:16 [DEBUG] Checking Gaps
2017-03-30 09:38:16 [DEBUG] YP.NE49.Z Slicing Stream to 2010-03-11T00:00:00.000000Z:2010-03-11T23:59:59.800000Z
2017-03-30 09:38:16 [DEBUG] YP.NE49.Z Highpass at 0.01 Hz
2017-03-30 09:38:16 [INFO] Processing pair: YP.NE49 vs YP.NE49
2017-03-30 09:38:16 [DEBUG] Azimuth=0.0
Process Process-1:
Traceback (most recent call last):
File "/home/sj/anaconda2/lib/python2.7/multiprocessing/process.py", line 258, in _bootstrap
self.run()
File "/home/sj/anaconda2/lib/python2.7/multiprocessing/process.py", line 114, in run
self._target(*self._args, **self._kwargs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/s03compute_cc.py", line 690, in main
corr = stack(db, corrs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line 982, in stack
coh = np.convolve(sp.signal.boxcar(timegate_samples)/timegate_samples, coh, 'same')
File "/home/sj/anaconda2/lib/python2.7/site-packages/scipy/signal/windows.py", line 59, in boxcar
return np.ones(M, float)
File "/home/sj/anaconda2/lib/python2.7/site-packages/numpy/core/numeric.py", line 192, in ones
a = empty(shape, dtype, order)
TypeError: 'float' object cannot be interpreted as an index
Hi,Thomas and everyone ,
I compute the autocorr ZZ and autocorr=Y in the config,and it worked,and generate a foder named STACKS which include 01/005_days/ZZ/YP_NE7D_YP_NE7D/*SAC.
but run
$msnoise plot interferogram YP.NE7D YP.NE7D -f01 -cZZ -m5
result in nothing,no error .
many thanks!
From: Thomas Lecocq<mailto:thomas.lecocq@oma.be>
Date: 2017-03-28 23:01
To: msnoise<mailto:msnoise@mailman-as.oma.be>
Subject: Re: [MSNoise] error in autocorrelation plotting
Hi,
did you compute the autocorr ZZ (autocorr=Y in the config) and it worked
OK during compute_cc ?
Thomas
Le 28/03/2017 à 16:52, sun jie a écrit :
> Hi,Thomas and everyone
>
> First ,the problems of interferogram ,ccftime... plotting has been solved,thank you very much!
> Second,the interferogram,ccftime... plotting only can plot the cross-correlation data.The autocorrelation data cannot be plotted,
> run
> $ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
> result in nothing,
> no error,no anything .
>
> how can i solve this problem? many thanks.
>
>
> _______________________________________________
> MSNoise mailing list
> MSNoise(a)mailman-as.oma.be<mailto:MSNoise@mailman-as.oma.be>
> http://mailman-as.oma.be/mailman/listinfo/msnoise
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be<mailto:MSNoise@mailman-as.oma.be>
http://mailman-as.oma.be/mailman/listinfo/msnoise
Hi,Thomas and everyone ,
I��m Sun Jie ,i changed the email adress to this.
I compute the autocorr ZZ and autocorr=Y in the config,and it worked,and generate a foder named STACKS which include 01/005_days/ZZ/YP_NE7D_YP_NE7D/*SAC.
but run
$msnoise plot interferogram YP.NE7D YP.NE7D -f01 -cZZ -m5
result in nothing,no error .
many thanks!
From: Thomas Lecocq
Date: 2017-03-28 23:01
To: msnoise
Subject: Re: [MSNoise] error in autocorrelation plotting
Hi,
did you compute the autocorr ZZ (autocorr=Y in the config) and it worked
OK during compute_cc ?
Thomas
Le 28/03/2017 �� 16:52, sun jie a ��crit :
> Hi,Thomas and everyone
>
> First ,the problems of interferogram ,ccftime... plotting has been solved,thank you very much!
> Second,the interferogram,ccftime... plotting only can plot the cross-correlation data.The autocorrelation data cannot be plotted,
> run
> $ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
> result in nothing,
> no error,no anything .
>
> how can i solve this problem? many thanks.
>
>
> _______________________________________________
> 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
�&N&j��r�*�謚�h�&�Y�Lj��Ӟx㭴״^���=�aN&j��r�*�謚�h�&�\��)ڕ:&�z"�隊Y���(���
�b���u:&�z"�隊Y���(���E� ���~�&�I��ɚ�Z.����u�j:�f������gz�:�N{�m7o)��Y���(���>�-~,p��DH�ωݺ���N~��z"�隊Y���(���Y�v�ƫ�M{Ӟx㯴�ML$^q�y��jek'��!j���[���v�}|�}������&j)h��/�w�>�-~,p��DH�ω����
�~��z"�隊Y���(���Y�v�ƫ�M{Ӿx㟴�M$�yǢ��_����mz�M���u��y�m��gr���r�(m�x�^x�m6� ����ƭy� �����Kh���ʙZ����HZ��gz�������D�L�л��u�M;��+���Ǧj)fjv��f�yg����m5�N��~�m4D�N�����zf��f�j�&i��y��zw ��:u����u^[ז�N]
�ޙ���Mt�}w߇.�����Ů���������ʓ�������ʬz+&�Z ɷ�z����^<��:i��q�wƻk�]�W�s���謚�h�&�
�^Y�v�ƫ�M{Ӿx㟴�M������2��V��E���0����E�[�+��t�$����M5�T�w^���w�9�C0E^�Ȩ�]���$��Νy�;{�_ז�巟ӗB���j)u�]5�]��ˠm�nrp��ק�<�zױ��ej)܅��z���+��
Hi,Thomas and everyone
First ,the problems of interferogram ,ccftime... plotting has been solved,thank you very much!
Second,the interferogram,ccftime... plotting only can plot the cross-correlation data.The autocorrelation data cannot be plotted,
run
$ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
result in nothing,
no error,no anything .
how can i solve this problem? many thanks.
Hi,Thomas and everyonw.
when i run $ msnoise plot data_availability ,result in ok;
And run
$ msnoise plot station_map
$ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
$ msnoise plot ccftime...
$ msnoise plot mwcs...
$ msnoise plot dtt...
result in the same problems such as $ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
But run $ msnoise plot dvv...,result in ok
The error report is described below.
how can I solve this problem? many thanks
Ubuntu 16.04 LTS 64bit msnoise 1.4
$ msnoise plot station_map
result in
Traceback (most recent call last):
File "/home/sj/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 722, in __call__
return self.main(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 535, in invoke
return callback(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 535, in station_map
from ..plots.station_map import main
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/plots/station_map.py", line 36, in <module>
from mpl_toolkits.basemap import Basemap
ImportError: No module named basemap
run
$ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
result in
Traceback (most recent call last):
File "/home/sj/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 722, in __call__
return self.main(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 535, in invoke
return callback(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 456, in interferogram
from ..plots.interferogram import main
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/plots/interferogram.py", line 21, in <module>
from scipy.stats.stats import nanmean
ImportError: cannot import name nanmean
Hi,Thomas and everyonw.
when i run $ msnoise plot data_availability ,result in ok;
And run
$ msnoise plot station_map
$ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
$ msnoise plot ccftime...
$ msnoise plot mwcs...
$ msnoise plot dtt...
result in the same problems such as $ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
But run $ msnoise plot dvv...,result in ok
The error report is described below.
how can I solve this problem? many thanks
Ubuntu 16.04 LTS 64bit msnoise 1.4
$ msnoise plot station_map
result in
Traceback (most recent call last):
File "/home/sj/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 722, in __call__
return self.main(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 535, in invoke
return callback(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 535, in station_map
from ..plots.station_map import main
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/plots/station_map.py", line 36, in <module>
from mpl_toolkits.basemap import Basemap
ImportError: No module named basemap
run
$ msnoise plot interferogram YP.NE12 YP.NE12 -f01 -cZZ -m5
result in
Traceback (most recent call last):
File "/home/sj/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 614, in run
cli(obj={})
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 722, in __call__
return self.main(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 535, in invoke
return callback(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/decorators.py", line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 456, in interferogram
from ..plots.interferogram import main
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/plots/interferogram.py", line 21, in <module>
from scipy.stats.stats import nanmean
ImportError: cannot import name nanmean
Hi, Thomas
Sorry to reply so late.
I tried to uninstall numpy and use conda to install again, it seems no problem.
But I tried again to run compute_cc command, some error occurs again.
The bug report shows:
Process Process-1:
Traceback (most recent call last):
File "/public/software/anaconda3/lib/python3.5/multiprocessing/process.py", line 249, in _bootstrap
self.run()
File "/public/software/anaconda3/lib/python3.5/multiprocessing/process.py", line 93, in run
self._target(*self._args, **self._kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/s03compute_cc.py", line 653, in main
trame2h[i]*cp, Nfft, dt, low, high, plot=False)
File "/public/software/anaconda3/lib/python3.5/site-packages/msnoise/whiten.py", line 39, in whiten
freqVec = scipy.fftpack.fftfreq(Nfft,d=delta)[:Nfft/2]
TypeError: slice indices must be integers or None or have an __index__ method
Maybe it's because I set the wrong filter? But I do not know how to correct it.
Best Wishes,
Shuye Huang
> -----原始邮件-----
> 发件人: "Thomas Lecocq" <thomas.lecocq(a)oma.be>
> 发送时间: 2017-03-20 16:08:31 (星期一)
> 收件人: msnoise(a)mailman-as.oma.be
> 抄送:
> 主题: [SPAM] Re: [MSNoise] Error when computing_cc
>
> Hi Shuye,
>
> Never seen this error, it is indeed something wrong with numpy. It looks
> like numpy and its dependencies haven't been installed the same way (or
> installed at all). you could remove numpy (pip uninstall numpy) and,
> assuming you are using anaconda, use "conda install numpy"
>
> Thomas
>
>
> Le 20/03/2017 à 08:17, hsy119(a)mail.ustc.edu.cn a écrit :
> > Hi, Thomas and users
> >
> > I just got an error when I tried to run compute_cc command. It shows:
> >
> > [hsy119@yhjs3 msnoise_my]$ msnoise compute_cc
> > 2017-03-20 15:18:30 [INFO] *** Starting: Compute CC ***
> > 2017-03-20 15:18:30 [INFO] Will compute ZZ RR TT TR RT ZR RZ TZ ZT
> > 2017-03-20 15:18:30 [INFO] New CC Job: 2015-01-10 (3 pairs with 3 stations)
> > 2017-03-20 15:18:31 [DEBUG] G2.XJ16.Z Reading 1 Files
> > 2017-03-20 15:18:31 [DEBUG] Checking sample alignment
> > 2017-03-20 15:18:31 [DEBUG] Checking Gaps
> > Intel MKL FATAL ERROR: Cannot load libmkl_def.so.
> >
> > Is there something wrong with numpy?
> >
> > Best Regards,
> > Shuye Huang
> > _______________________________________________
> > 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
