[MSNoise] About s03compute_cc

[Thomas Lecocq]

Carmelo,

it depends on what you want to achieve, +with respect to Nyquist/Shannon 
etc...

if you want to cross-correlate at 1Hz, you'll have information up to 
0.5Hz, you have to preprocess-lowpass below that, and decimate by a 
factor of 125...

About the bug: I don't know, please fill in an Issue on the github, so I 
can have a look at it later.

Thomas

Le 09/02/2015 13:35, Sammarco, Carmelo a écrit :
> thanks Thomas for this info.
>
> 1) yes.. 125 Hz
> 2) the configured cc_sampling_rate is 1 HZ
>
> So what can you suggest me to run correctly the computations? What decimation factor and/or resample value I need to use?? May explain me in witch way I can understand what is the correct parameters?
>
> Is it possible resolve the bug or ignore it to finish the computation?
>
> thank you very much.
>
> -----------
> Mr Carmelo Sammarco MSc  BSc  FGS
> PhD Student in Geology and Petroleum Geology
>
> ________________________________________
> From: msnoise-bounces at mailman-as.oma.be <msnoise-bounces at mailman-as.oma.be> on behalf of Thomas Lecocq <thomas.lecocq at oma.be>
> Sent: 09 February 2015 11:57
> To: msnoise at mailman-as.oma.be
> Subject: Re: [MSNoise] About s03compute_cc
>
> Hi Carmelo,
>
> 1/ It looks like you are using 125 Hz data, is that right ?
> 2/ What is the configured "cc_sampling_rate" value ? I suppose it's 1 Hz.
>
> The first error is because you can't decimate by a factor of 5 to go
> from 125 to 1Hz.
> The second looks like an error in the resample method which gives back 2
> samples less than expected, and it's a bug.
>
> Thomas
>
> Le 09/02/2015 12:40, Sammarco, Carmelo a écrit :
>> Dear all,
>>
>>
>> I am running the CC process. The terminal show me this:
>>
>>
>> carmelo at COPS-023437:~/msnoise/FLEG$ python s03compute_cc.py
>> 2015-02-09 11:14:40,482 [INFO] *** Starting: Compute CC ***
>> 2015-02-09 11:14:40,509 [INFO] Will compute ZZ
>> 2015-02-09 11:14:40,572 [INFO] New CC Job: 2011-03-01 (1 pairs with 2 stations)
>> 2015-02-09 11:14:40,587 [DEBUG] FLEG.ACL2.Z Reading 24 Files
>> 2015-02-09 11:14:41,207 [DEBUG] FLEG.ACL2.Z Merging Stream
>> 2015-02-09 11:14:41,227 [DEBUG] FLEG.ACL2.Z Slicing Stream to 2011-03-01T00:00:00.000000Z:2011-03-01T23:59:59.992000Z
>> 2015-02-09 11:14:41,227 [DEBUG] FLEG.ACL2.Z Lowpass at 0.20 Hz
>> 2015-02-09 11:14:41,655 [DEBUG] FLEG.ACL2.Z Highpass at 0.05 Hz
>> 2015-02-09 11:14:42,082 [DEBUG] FLEG.ACL2.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 (2160000) into shape (86400)
>>
>>
>> -----------------------------------
>>
>>
>> carmelo at COPS-023437:~/msnoise/FLEG$ python s03compute_cc.py
>> 2015-02-09 11:13:07,728 [INFO] *** Starting: Compute CC ***
>> 2015-02-09 11:13:07,752 [INFO] Will compute ZZ
>> 2015-02-09 11:13:07,831 [INFO] New CC Job: 2011-03-01 (1 pairs with 2 stations)
>> 2015-02-09 11:13:07,849 [DEBUG] FLEG.ACL2.Z Reading 24 Files
>> 2015-02-09 11:13:08,471 [DEBUG] FLEG.ACL2.Z Merging Stream
>> 2015-02-09 11:13:08,492 [DEBUG] FLEG.ACL2.Z Slicing Stream to 2011-03-01T00:00:00.000000Z:2011-03-01T23:59:59.992000Z
>> 2015-02-09 11:13:08,492 [DEBUG] FLEG.ACL2.Z Lowpass at 0.20 Hz
>> 2015-02-09 11:13:08,921 [DEBUG] FLEG.ACL2.Z Highpass at 0.05 Hz
>> 2015-02-09 11:13:09,347 [DEBUG] FLEG.ACL2.Z Downsample to 1.0 Hz
>> 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 (86398) into shape (86400)
>>
>>
>>
>> The first I tried to do with decimate.. the second one with resample method.
>>
>>
>>
>> Someone know in witch way can I fix this error or give me some advice?
>>
>>
>> thank you very much for your time.
>>
>>
>> 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
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
>
>
> 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