[MSNoise] Stuck in Compute_cc step

[Sébastien Carniato]

Thomas,

It works with that config, but with that filter, I get once again a symmetric wave (like the one I put in my previous email)… The only change is the average amplitude, which is around -2.4e-11.
It does look like an autocorrelation. But why ?

Sébastien.
On 15 mai 2014, at 12:38, Thomas Lecocq <thomas.lecocq at oma.be> wrote:

> Sébastien,
> 
> yes sorry, the high & low vocabulary is not clear.
> 
> for the filters:
> 
> low : 1
> mwcs_low : 1
> high : 45
> mwcs_high : 45
> rms_treshold : 0
> mwcs_len : 10
> mwcs_step : 5
> 
> Should work
> 
> Thomas
> 
> Le 15/05/2014 11:59, Sébastien Carniato a écrit :
>> Thomas,
>> 
>> I changed my configurator which looks like this now :
>> 
>> 
>> startdate : 1970-01-01
>> enddate : 2100-01-01
>> analysis duration : 86400
>> cc_sampling_rate : 200
>> resampling_method : Decimate
>> decimation_factor : 2
>> preprocess_lowpass : 45
>> preprocess_highpass : 1
>> maxlag : 120
>> corr_duration : 1800
>> windsorizing : 3
>> crondays : -1
>> autocorr : N
>> keep_all : Y
>> keep_days : Y
>> ref_begin : -100
>> ref_end : 0
>> export_format : BOTH
>> sac_format : doublets
>> 
>> And my filter :
>> low : 45
>> mwcs_low : 45
>> high : 1
>> mwcs_high : 1
>> rms_treshold : 0
>> mwcs_len : 10
>> mwcs_step : 5
>> I’m not sure of getting what you meant by “same for the following”.
>> 
>> Now, I get the error :
>> 
>> 2014-05-15 11:52:01,333 [DEBUG] Processing pair: AN.LF1 vs AN.LF2
>> <msnoise_table_def.Station object at 0x4a94950>
>> <msnoise_table_def.Station object at 0x4a9f110>
>> Traceback (most recent call last):
>>   File "s03compute_cc.py", line 405, in <module>
>>     trame2h[i], Nfft, dt, low, high, plot=False)
>>   File "/home/stag01/Bureau/MSNoise-1.2.3/whiten.py", line 58, in whiten
>>     low = J[0] - Napod
>> IndexError: index out of bounds
>> 
>> Thanks.
>> Sébastien
>> 
>> On 15 mai 2014, at 11:22, Thomas Lecocq <thomas.lecocq at oma.be> wrote:
>> 
> 
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise