Hi Sun Jie,
It's again a numpy-too-recent problem (and our mistake to leave non-int
values where they should be). I've now corrected this in the development
version, in the meantime just replace:
line 982, in stack
coh = np.convolve(sp.signal.boxcar(timegate_samples)/timegate_samples, coh,
'same')
by
line 982, in stack
coh = np.convolve(sp.signal.boxcar(int(timegate_samples))/timegate_samples, coh,
'same')
and it should work.
Thomas
Le 30/03/2017 à 04:06, 孙杰 a écrit :
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
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise