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(a)mailman-as.oma.be [mailto:msnoise-bounces@mailman-as.oma.be] En nombre
de Thomas Lecocq
Enviado el: viernes, 20 de febrero de 2015 10:21
Para: msnoise(a)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(a)mailman-as.oma.be
[mailto:msnoise-bounces@mailman-as.oma.be] En nombre de Thomas Lecocq
Enviado el: jueves, 19 de febrero de 2015 9:19
Para: msnoise(a)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(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be