30 Mar
2017
30 Mar
'17
4:06 a.m.
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