Dear all,
I am running the CC process. The terminal show me this:
carmelo@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@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.
Hi Fernando,
Can you tell me what is filter 03 ? and at which frequency do you
calculate the CCFs ?
An usual reason for failure at this point is that you don't have enough
points in the calculation to compute the weighted linear regression...
Thomas
Le 8/02/2015 17:44, Fernando Salazar Monroy a écrit :i
> Greetings for all,
>
>
> I'm trying estimate the s05compute_mwcs.py for all my data base, when
> I run this code, this appeared:
>
>
> 2015-02-08 10:33:39,545 [DEBUG] Processing MWCS for:
> IG_ARIG_IG_PLIG.ZZ.03 - 2011-12-29 - 10 days
> Traceback (most recent call last):
> File "s05compute_mwcs.py", line 118, in <module>
> cur, ref, f.mwcs_low, f.mwcs_high, goal_sampling_rate, -maxlag,
> f.mwcs_wlen, f.mwcs_step)
> File "/media/fernando/Seagate Expansion
> Drive/MSNOISE/GGAP_SSN/MWCS.py", line 151, in mwcs
> res = sm.regression.linear_model.WLS(phi, v, w**2).fit()
> File
> "/usr/local/lib/python2.7/dist-packages/statsmodels-0.5.0-py2.7-linux-x86_64.egg/statsmodels/regression/linear_model.py",
> line 381, in __init__
> weights=weights, hasconst=hasconst)
> File
> "/usr/local/lib/python2.7/dist-packages/statsmodels-0.5.0-py2.7-linux-x86_64.egg/statsmodels/regression/linear_model.py",
> line 79, in __init__
> super(RegressionModel, self).__init__(endog, exog, **kwargs)
> File
> "/usr/local/lib/python2.7/dist-packages/statsmodels-0.5.0-py2.7-linux-x86_64.egg/statsmodels/base/model.py",
> line 136, in __init__
> super(LikelihoodModel, self).__init__(endog, exog, **kwargs)
> File
> "/usr/local/lib/python2.7/dist-packages/statsmodels-0.5.0-py2.7-linux-x86_64.egg/statsmodels/base/model.py",
> line 52, in __init__
> self.data = handle_data(endog, exog, missing, hasconst, **kwargs)
> File
> "/usr/local/lib/python2.7/dist-packages/statsmodels-0.5.0-py2.7-linux-x86_64.egg/statsmodels/base/data.py",
> line 397, in handle_data
> return klass(endog, exog=exog, missing=missing, hasconst=hasconst,
> **kwargs)
> File
> "/usr/local/lib/python2.7/dist-packages/statsmodels-0.5.0-py2.7-linux-x86_64.egg/statsmodels/base/data.py",
> line 78, in __init__
> self._check_integrity()
> File
> "/usr/local/lib/python2.7/dist-packages/statsmodels-0.5.0-py2.7-linux-x86_64.egg/statsmodels/base/data.py",
> line 246, in _check_integrity
> if len(self.exog) != len(self.endog):
> TypeError: len() of unsized object
>
>
> I haven't modified anything inside of s05compute_mwcs.py, I don't
> understand the error.
>
>
> I hope that somebody help me with this.
>
>
>
> Edilson Fernando Salazar Monroy
> Master in Earth Sciences
> National Autonomus University of Mexico
Dear Thomas Lecocq,
I am writing to find out if i can use this application to compute empirical
green functions by performing cross correlation of ambient noise data
already acquired for a year.
If it is possible, i will appreciate some guidance.
Many thanks.
Ojo, Adebayo
I am writing to find out if i can use this application to compute empirical
green functions by performing cross correlation of ambient noise data
already acquired for a year.
If it is possible, i will appreciate some guidance.
Many thanks.
Ojo, Adebayo
