[MSNoise] error in pws

孙杰 1010131334 at cugb.edu.cn
Thu Mar 30 02:06:04 UTC 2017


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


More information about the MSNoise mailing list