[MSNoise] Problems with s03compute_cc

Fri Feb 27 16:35:26 UTC 2015

Dear Thomas,
Indeed. I had an error in the filter values and resampling configuration.
Right now the problem is solved and the MSNoise software work perfectly.
Thank you very much
Juan

-----Mensaje original-----
De: msnoise-bounces at mailman-as.oma.be [mailto:msnoise-bounces at mailman-as.oma.be] En nombre de Thomas Lecocq
Enviado el: viernes, 20 de febrero de 2015 10:21
Para: msnoise at mailman-as.oma.be
Asunto: Re: [MSNoise] Problems with s03compute_cc

Juan

this problem is different from your previous one, here it seems you are trying to fit 10Hz data in a 20Hz prepared "white trace".

What is your data sampling_rate ? what is the configuration value for "cc_sampling_rate" ?

Thomas

Le 19/02/2015 17:17, Rueda Nuñez Juán José a écrit :
> Dear Thomas.
> In this moment the filter configuration I use :
>
> low : 0.01 or 0.1
> mwcs_low : 0.12
> high : 0.98
> mwcs_high : 1.0
> rms_treshold : 0
> mwcs_len : 10
> mwcs_step : 5
>
> And the result is:
>
> 2015-02-19 17:03:04,831 [DEBUG] ES.CGIN.Z Reading 1 Files
> 2015-02-19 17:03:06,029 [DEBUG] ES.CGIN.Z Merging Stream
> 2015-02-19 17:03:06,030 [DEBUG] ES.CGIN.Z Slicing Stream to 
> 2015-01-05T00:00:00.000000Z:2015-01-05T23:59:59.990000Z
> 2015-02-19 17:03:06,109 [DEBUG] ES.CGIN.Z Lowpass at 8.00 Hz
> 2015-02-19 17:03:08,506 [DEBUG] ES.CGIN.Z Highpass at 0.01 Hz
> 2015-02-19 17:03:09,214 [DEBUG] ES.CGIN.Z Decimate by a factor of 5
> 2015-02-19 17:03:09,263 [DEBUG] ES.CGUI.Z Reading 1 Files
> 2015-02-19 17:03:10,589 [DEBUG] ES.CGUI.Z Merging Stream
> 2015-02-19 17:03:10,590 [DEBUG] ES.CGUI.Z Slicing Stream to 
> 2015-01-05T00:00:00.000000Z:2015-01-05T23:59:59.990000Z
> 2015-02-19 17:03:10,590 [DEBUG] ES.CGUI.Z Lowpass at 8.00 Hz
> 2015-02-19 17:03:11,266 [DEBUG] ES.CGUI.Z Highpass at 0.01 Hz
> 2015-02-19 17:03:11,972 [DEBUG] ES.CGUI.Z Decimate by a factor of 5
> 2015-02-19 17:03:12,019 [DEBUG] ES.CHIE.Z Reading 1 Files
> 2015-02-19 17:03:13,379 [DEBUG] ES.CHIE.Z Merging Stream
> 2015-02-19 17:03:13,379 [DEBUG] ES.CHIE.Z Slicing Stream to 
> 2015-01-05T00:00:00.000000Z:2015-01-05T23:59:59.980000Z
> 2015-02-19 17:03:13,380 [DEBUG] ES.CHIE.Z Lowpass at 8.00 Hz
> 2015-02-19 17:03:13,717 [DEBUG] ES.CHIE.Z Highpass at 0.01 Hz
> 2015-02-19 17:03:14,066 [DEBUG] ES.CHIE.Z Decimate by a factor of 5 
> Traceback (most recent call last):
>    File "s03compute_cc.py", line 285, in <module>
>      tramef_Z[istation] = trame
> ValueError: could not broadcast input array from shape (864000) into 
> shape (1728000)
>
> Thanks.
> Juan
>
>
> -----Mensaje original-----
> De: msnoise-bounces at mailman-as.oma.be 
> [mailto:msnoise-bounces at mailman-as.oma.be] En nombre de Thomas Lecocq 
> Enviado el: jueves, 19 de febrero de 2015 9:19
> Para: msnoise at mailman-as.oma.be
> Asunto: Re: [MSNoise] Problems with s03compute_cc
>
> Hi,
>
> what are the frequency bounds for the configured filter(s) ?
>
> Thomas
>
> Le 18/02/2015 18:20, Rueda Nuñez Juán José a écrit :
>> Dear all,
>> I have a problem with the cross-correlation computation.
>> The output for the program is:
>>
>>
>> 2015-02-18 17:56:31,075 [INFO] Will compute ZZ
>> 2015-02-18 17:56:31,306 [INFO] New CC Job: 2015-01-01 (4 pairs with 5
>> stations)
>> 2015-02-18 17:56:31,408 [DEBUG] ES.CJUL.Z Reading 1 Files
>> 2015-02-18 17:56:33,837 [DEBUG] ES.CJUL.Z Merging Stream
>> 2015-02-18 17:56:33,839 [DEBUG] ES.CJUL.Z Slicing Stream to 
>> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
>> 2015-02-18 17:56:33,937 [DEBUG] ES.CJUL.Z Lowpass at 8.00 Hz
>> 2015-02-18 17:56:38,981 [DEBUG] ES.CJUL.Z Highpass at 0.01 Hz
>> 2015-02-18 17:56:39,705 [DEBUG] ES.CJUL.Z Decimate by a factor of 5
>> 2015-02-18 17:56:39,758 [DEBUG] ES.CORC.Z Reading 1 Files
>> 2015-02-18 17:56:42,276 [DEBUG] ES.CORC.Z Merging Stream
>> 2015-02-18 17:56:42,277 [DEBUG] ES.CORC.Z Slicing Stream to 
>> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
>> 2015-02-18 17:56:42,277 [DEBUG] ES.CORC.Z Lowpass at 8.00 Hz
>> 2015-02-18 17:56:42,924 [DEBUG] ES.CORC.Z Highpass at 0.01 Hz
>> 2015-02-18 17:56:43,625 [DEBUG] ES.CORC.Z Decimate by a factor of 5
>> 2015-02-18 17:56:43,670 [DEBUG] ES.CRST.Z Reading 1 Files
>> 2015-02-18 17:56:46,168 [DEBUG] ES.CRST.Z Merging Stream
>> 2015-02-18 17:56:46,168 [DEBUG] ES.CRST.Z Slicing Stream to 
>> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
>> 2015-02-18 17:56:46,169 [DEBUG] ES.CRST.Z Lowpass at 8.00 Hz
>> 2015-02-18 17:56:46,817 [DEBUG] ES.CRST.Z Highpass at 0.01 Hz
>> 2015-02-18 17:56:47,478 [DEBUG] ES.CRST.Z Decimate by a factor of 5
>> 2015-02-18 17:56:47,525 [DEBUG] ES.CTAB.Z Reading 1 Files
>> 2015-02-18 17:56:50,316 [DEBUG] ES.CTAB.Z Merging Stream
>> 2015-02-18 17:56:50,443 [DEBUG] ES.CTAB.Z Slicing Stream to 
>> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
>> 2015-02-18 17:56:50,445 [DEBUG] ES.CTAB.Z Lowpass at 8.00 Hz
>> 2015-02-18 17:56:51,115 [DEBUG] ES.CTAB.Z Highpass at 0.01 Hz
>> 2015-02-18 17:56:51,775 [DEBUG] ES.CTAB.Z Decimate by a factor of 5
>> 2015-02-18 17:56:51,815 [DEBUG] ES.CTIG.Z Reading 1 Files
>> 2015-02-18 17:56:54,611 [DEBUG] ES.CTIG.Z Merging Stream
>> 2015-02-18 17:56:55,737 [DEBUG] ES.CTIG.Z Slicing Stream to 
>> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
>> 2015-02-18 17:56:55,738 [DEBUG] ES.CTIG.Z Lowpass at 8.00 Hz
>> 2015-02-18 17:56:56,403 [DEBUG] ES.CTIG.Z Highpass at 0.01 Hz
>> 2015-02-18 17:56:57,074 [DEBUG] ES.CTIG.Z Decimate by a factor of 5
>> 2015-02-18 17:56:57,113 [DEBUG] Processing pair: ES.CJUL vs ES.CTIG 
>> Traceback (most recent call last):
>>     File "s03compute_cc.py", line 402, in <module>
>>       trame2h[i]*cp, Nfft, dt, low, high, plot=False)
>>     File "/home/juan/MSNoise-1.2.5/whiten.py", line 58, in whiten
>>       low = J[0] - Napod
>> IndexError: index 0 is out of bounds for axis 0 with size 0
>>
>>
>> Any solution?
>>
>> Thank you very much for your time.
>>
>> Juan Rueda
>> Instituto Geográfico Nacional
>> Spain
>> _______________________________________________
>> MSNoise mailing list
>> MSNoise at mailman-as.oma.be
>> http://mailman-as.oma.be/mailman/listinfo/msnoise
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise

_______________________________________________
MSNoise mailing list
MSNoise at mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise