Hi Thomas,
The filter number 3 is of 40 until 20 seconds and for the correlations are Preprocess_highpass is 0.01 and lowpass 8.0, therefore I want ask why the correlations have a duration of 240 seconds?
Edilson Fernando Salazar Monroy
Master in Earht Sciences
National Autonomus University of Mexico
Dear All,
After much attempt to install the MSnoise software on my Centos 7 Linux
system, i feel i should ask for some guidance.
[root@localhost ojo]# sudo yum install libsamplerate0 libsamplerate0-dev
Loaded plugins: fastestmirror, langpacks, priorities
Loading mirror speeds from cached hostfile
* base: centos.ustc.edu.cn
* epel: mirrors.ustc.edu.cn
* extras: centos.ustc.edu.cn
* rpmforge: mirror-fpt-telecom.fpt.net
* updates: centos.ustc.edu.cn
No package libsamplerate0 available.
No package libsamplerate0-dev available.
[root@localhost ojo]# easy_install scikits.samplerate
Searching for scikits.samplerate
Reading https://pypi.python.org/simple/scikits.samplerate/
Best match: scikits.samplerate 0.3.3
Downloading
https://pypi.python.org/packages/source/s/scikits.samplerate/scikits.sample…
Processing scikits.samplerate-0.3.3.tar.gz
Writing /tmp/easy_install-aF6CLy/scikits.samplerate-0.3.3/setup.cfg
Running scikits.samplerate-0.3.3/setup.py -q bdist_egg --dist-dir
/tmp/easy_install-aF6CLy/scikits.samplerate-0.3.3/egg-dist-tmp-L4szBy
error: SRC (http://www.mega-nerd.com/SRC/) library not found. Directories
to search
for the libraries can be specified in the site.cfg file, in section
[samplerate].
Kindly advice me. Thank you.
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
Dear all,
When I try to run some script now appear this:
/home/carmelo/anaconda/lib/python2.7/site-packages/setuptools-5.8-py2.7.egg/pkg_resources.py:1032: UserWarning: /home/carmelo/.python-eggs is writable by group/others and vulnerable to attack when used with get_resource_filename. Consider a more secure location (set with .set_extraction_path or the PYTHON_EGG_CACHE environment variable).
may someone help me please?? I have no idea what happen. Before i run the script without problem..
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.
Hello,
I am using 20sps data and everything is fine until I reach *s06compute_dtt.py.
*I changed the mwcs_len parameter to 30 and the *mwcs_step* to 10. When I
do that I get the next error message whan I run s06compute_dtt.py:
*2015-01-28 11:30:17,546 [DEBUG] Found 731 updated days*
*2015-01-28 11:30:17,556 [INFO] Loading mov=20 days for filter=03*
*2015-01-28 11:30:17,556 [DEBUG] Processing 2005-01-01*
*Traceback (most recent call last):*
* File "s06compute_dtt.py", line 187, in <module>*
* dtArray = np.vstack((dtArray, df['dt']))*
* File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/numpy/core/shape_base.py",
line 228, in vstack*
* return _nx.concatenate([atleast_2d(_m) for _m in tup], 0)*
*ValueError: all the input array dimensions except for the concatenation
axis must match exactly*
Does anyone know how to fix it? Thanks in advance!
Oscar
Hi Thomas,
I didn't send this to the entire mail list because they are comments on the
master (not release) of msnoise. Thanks for all the work you have been
doing on msnoise!
In api.py, module check_and_phase_shift:
Any modifications to trace result in trace.data containing floats. If
trace.data originally contained ints this presents an issue if there was
"No Offset" because ints can't later be merged with floats. My quick
workaround was:
Dear MSnoisers:
In step s07, I have a problem with the resampling in pandas. The result is
that I lose the errorbars on my dv/v's.
py1_wmean, py1_wstd = get_wavgwstd(allbut)
py1_wmean = py1_wmean.resample('D', how='mean')
py1_wstd = py1_wstd.resample('D', how='mean').fillna(0.0)
data = detrend(py1_wmean)
print 'data= ', data
[kasper@localhost MSNoise-master]$ python s07plot_dtt.py
data= [ nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan nan nan nan nan nan
nan nan nan nan]
If I don't resample, life is good:
py1_wmean, py1_wstd = get_wavgwstd(allbut)
#py1_wmean = py1_wmean.resample('D', how='mean')
#py1_wstd = py1_wstd.resample('D', how='mean').fillna(0.0)
data = detrend(py1_wmean)
print 'data= ', data
[kasper@localhost MSNoise-master]$ python s07plot_dtt.py
data= [ 1.05332051e-02 1.24279614e-02 -4.84688864e-01 1.11580621e-02
8.28039210e-03 1.04757327e-02 1.07993328e-02 -5.96110125e-04
7.23556870e-03 9.96330072e-03 9.54712307e-03 1.10273139e-02
1.06786539e-02 1.11155446e-02 9.51223128e-03 1.05238797e-02
1.00996973e-02 9.37587307e-03 7.18620810e-03 8.97052425e-03
9.08720458e-03 5.78305838e-03 8.33841882e-03 9.67250220e-03
8.61533043e-03 9.87386678e-03 8.22291954e-03 1.00831938e-02
9.36887461e-03 6.17043697e-03 8.22567034e-03 8.65890679e-03
8.22813898e-03 7.76614092e-03 7.84201987e-03 6.62494706e-03
8.02592632e-03 7.57269605e-03 8.68610949e-03 8.13714761e-03
7.87601472e-03 7.59459241e-03 7.61133449e-03 7.23192283e-03
Cheers,
kasper
