[MSNoise] CC problem

[Thomas Lecocq]

Hey Carmelo,

It's a problem with the Resample (scikits-samplerate) routine, you 
should use Decimate if possible (125->25 = factor 5 indeed)

Thomas

Le 15/05/2015 18:16, Sammarco, Carmelo a écrit :
> Dear Thomas and msnoise users,
>
>
> I have the following problem after running the CC process:
>
>
> carmelo at COPS-023437:~/Desktop/FLEGREI/msnoise$ msnoise compute_cc
> 2015-05-15 17:00:43 [INFO] *** Starting: Compute CC ***
> 2015-05-15 17:00:43 [INFO] Will compute ZZ
> 2015-05-15 17:00:43 [INFO] New CC Job: 2011-01-01 (28 pairs with 8 stations)
> Traceback (most recent call last):
>    File "/home/carmelo/anaconda/bin/msnoise", line 9, in <module>
>      load_entry_point('msnoise==1.3.1', 'console_scripts', 'msnoise')()
>    File "/home/carmelo/anaconda/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 393, in run
>      cli(obj={})
>    File "/home/carmelo/anaconda/lib/python2.7/site-packages/click/core.py", line 610, in __call__
>      return self.main(*args, **kwargs)
>    File "/home/carmelo/anaconda/lib/python2.7/site-packages/click/core.py", line 590, in main
>      rv = self.invoke(ctx)
>    File "/home/carmelo/anaconda/lib/python2.7/site-packages/click/core.py", line 936, in invoke
>      return _process_result(sub_ctx.command.invoke(sub_ctx))
>    File "/home/carmelo/anaconda/lib/python2.7/site-packages/click/core.py", line 782, in invoke
>      return ctx.invoke(self.callback, **ctx.params)
>    File "/home/carmelo/anaconda/lib/python2.7/site-packages/click/core.py", line 416, in invoke
>      return callback(*args, **kwargs)
>    File "/home/carmelo/anaconda/lib/python2.7/site-packages/msnoise/scripts/msnoise.py", line 174, in compute_cc
>      main()
>    File "/home/carmelo/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/carmelo/anaconda/lib/python2.7/site-packages/msnoise/s03compute_cc.py", line 191, in preprocess
>      tramef_Z[istation] = trace.data
> ValueError: could not broadcast input array from shape (2159998) into shape (2160000)
> carmelo at COPS-023437:~/Desktop/FLEGREI/msnoise$
>
>
> My data have a sample rate of 125 Hz and in the configurator I choose the resampling method = 25 Hz
>
> I try also with 20, 125 but I have always this problem
>
>
> The same  with decimation=5..
>
>
> some suggestions?? Every help is always welcome..
>
>
> Thank you very much!!
>
>
> cheers,
>
>
>
> -----------
> Mr Carmelo Sammarco MSc  BSc  FGS
> PhD Student in Geology and Petroleum Geology
>
>
> The University of Aberdeen is a charity registered in Scotland, No SC013683.
> Tha Oilthigh Obar Dheathain na charthannas clàraichte ann an Alba, Àir. SC013683.
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise