[MSNoise] msnoise compute_cc problems

[Thomas Lecocq]

Hi Silvio,

I have noticed this error too... I haven't had time to check it, but it 
happens when we try to merge two traces that need merging...

it's definitively a bug :-)

Tom

Le 2/04/2015 17:07, De Angelis, Silvio a écrit :
> Hi everyone,
>
> I am testing release 1.3.1 on Ubuntu 64bit (12.04). Things went pretty smooth for a first attempt with some CCF being calculated. Except every now and again the msnoise compute_cc stops working delivering the error message below, apparently linked to sampling rate issues?
>
> Here’s a link to a SDS data structure with data (two stations/two days) that generated the error:
>
> https://copy.com/blzIarTs4QrAllNr
>
> I have kept most options in the config default, just for testing purposes, except start and end date to match my data. I also added a filter in the config window (0.1 (low)- 0.3 (mwcs low)- 0.9 (mwcs high) -1.1 (high) Hz, 100 (RMS threshold), 1 (window), 0.5 (step)).
>
> msnoise compute_cc did not work unless a filter was defined in the configuration. I had 27 days worth of CCF's computed that I could plot, thus, I suppose the general config was ok? The computed CCF’s were what i expected (matched another workflow I use for ambient noise).
>
> Any hints?  The data files look ok as far as I can tell but there may be something I am missing?
>
> I am now testing with a different dataset to see whether problems are with data files.
>
> Thank you.
> Ciao,
> Silvio
> Silvio De Angelis
>
> Lecturer in Geophysics
> University of Liverpool
> Jane Herdman Building
> 4 Brownlow Street,
> Liverpool, Merseyside, L69 3GP
> United Kingdom
> email: silvioda at liverpool.ac.uk<mailto:silvioda at liverpool.ac.uk>
>
>
>
> 2015-04-02 15:22:48 [INFO] *** Starting: Compute CC ***
> 2015-04-02 15:22:48 [INFO] Will compute ZZ
> 2015-04-02 15:22:48 [INFO] New CC Job: 2006-01-25 (1 pairs with 2 stations)
> 2015-04-02 15:22:55 [INFO] Pre-Whitening Traces
> 2015-04-02 15:22:56 [INFO] Processing CC
> 2015-04-02 15:22:56 [INFO] Job Finished. It took 7.45 seconds
> 2015-04-02 15:22:56 [INFO] New CC Job: 2006-01-26 (1 pairs with 2 stations)
> 2015-04-02 15:23:03 [INFO] Pre-Whitening Traces
> 2015-04-02 15:23:03 [INFO] Processing CC
> 2015-04-02 15:23:03 [INFO] Job Finished. It took 7.37 seconds
> 2015-04-02 15:23:03 [INFO] New CC Job: 2006-01-27 (1 pairs with 2 stations)
> 2015-04-02 15:23:10 [INFO] Pre-Whitening Traces
> 2015-04-02 15:23:11 [INFO] Processing CC
> 2015-04-02 15:23:11 [INFO] Job Finished. It took 7.30 seconds
> 2015-04-02 15:23:11 [INFO] New CC Job: 2006-01-28 (1 pairs with 2 stations)
> Traceback (most recent call last):
>    File "/home/silvio/anaconda/bin/msnoise", line 9, in <module>
>      load_entry_point('msnoise==1.3.1', 'console_scripts', 'msnoise')()
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 393, in run
>      cli(obj={})
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/click-4.0-py2.7.egg/click/core.py", line 664, in __call__
>      return self.main(*args, **kwargs)
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/click-4.0-py2.7.egg/click/core.py", line 644, in main
>      rv = self.invoke(ctx)
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/click-4.0-py2.7.egg/click/core.py", line 991, in invoke
>      return _process_result(sub_ctx.command.invoke(sub_ctx))
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/click-4.0-py2.7.egg/click/core.py", line 837, in invoke
>      return ctx.invoke(self.callback, **ctx.params)
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/click-4.0-py2.7.egg/click/core.py", line 464, in invoke
>      return callback(*args, **kwargs)
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 174, in compute_cc
>      main()
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/msnoise/s03compute_cc.py", line 271, in main
>      basetime, tramef_Z = preprocess(db, stations, comps, goal_day, params, tramef_Z)
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/msnoise/s03compute_cc.py", line 130, in preprocess
>      stream[gap[0]] = stream[gap[0]].__add__(stream[gap[1]], method=0, fill_value="interpolate")
>    File "/home/silvio/anaconda/lib/python2.7/site-packages/obspy-0.10.1-py2.7-linux-x86_64.egg/obspy/core/trace.py", line 681, in __add__
>      raise TypeError("Sampling rate differs")
> TypeError: Sampling rate differs
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise