msnoise compute_cc problem - MSNoise - KSB-ORB Mailing Lists

3 Nov 2015

Hi,
I've installed MSNoise 1.3.1 on a Ubuntu 32bit computer, but when I try to compute
the cross- correlation functions it does it well for the first 3 years of data (from 1
January 2011 to 31 December 2013), but after that it reports this error:

2015-11-02 18:45:57 [INFO] *** Starting: Compute CC ***

/home/juan/anaconda/lib/python2.7/site-packages/setuptools-5.8-py2.7.egg/pkg_resources.py:1032:
UserWarning: /home/juan/.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).

2015-11-02 18:45:58 [INFO] Will compute ZZ ZT

2015-11-02 18:45:58 [INFO] New CC Job: 2014-01-15 (6 pairs with 4 stations)

Traceback (most recent call last):

  File "/home/juan/anaconda/bin/msnoise", line 9, in <module>

    load_entry_point('msnoise==1.3.1', 'console_scripts',
'msnoise')()

  File
"/home/juan/anaconda/lib/python2.7/site-packages/msnoise-1.3.1-py2.7.egg/msnoise/scripts/msnoise.py",
line 393, in run

    cli(obj={})

  File "/home/juan/anaconda/lib/python2.7/site-packages/click/core.py", line
700, in __call__

    return self.main(*args, **kwargs)

  File "/home/juan/anaconda/lib/python2.7/site-packages/click/core.py", line
680, in main

    rv = self.invoke(ctx)

  File "/home/juan/anaconda/lib/python2.7/site-packages/click/core.py", line
1027, in invoke

    return _process_result(sub_ctx.command.invoke(sub_ctx))

  File "/home/juan/anaconda/lib/python2.7/site-packages/click/core.py", line
873, in invoke

    return ctx.invoke(self.callback, **ctx.params)

  File "/home/juan/anaconda/lib/python2.7/site-packages/click/core.py", line
508, in invoke

    return callback(*args, **kwargs)

  File
"/home/juan/anaconda/lib/python2.7/site-packages/msnoise-1.3.1-py2.7.egg/msnoise/scripts/msnoise.py",
line 174, in compute_cc

    main()

  File
"/home/juan/anaconda/lib/python2.7/site-packages/msnoise-1.3.1-py2.7.egg/msnoise/s03compute_cc.py",
line 266, in main

    basetime, tramef_Z, tramef_E, tramef_N = preprocess(db, stations, comps, goal_day,
params, tramef_Z, tramef_E, tramef_N)

  File
"/home/juan/anaconda/lib/python2.7/site-packages/msnoise-1.3.1-py2.7.egg/msnoise/s03compute_cc.py",
line 144, in preprocess

    trace.detrend(type="linear")

  File
"/home/juan/anaconda/lib/python2.7/site-packages/obspy-0.9.2-py2.7-linux-i686.egg/obspy/core/util/decorator.py",
line 249, in new_func

    return func(*args, **kwargs)

  File
"/home/juan/anaconda/lib/python2.7/site-packages/obspy-0.9.2-py2.7-linux-i686.egg/obspy/core/util/decorator.py",
line 232, in new_func

    return func(*args, **kwargs)

  File
"/home/juan/anaconda/lib/python2.7/site-packages/obspy-0.9.2-py2.7-linux-i686.egg/obspy/core/trace.py",
line 1732, in detrend

    self.data = func(self.data, **options)

  File
"/home/juan/anaconda/lib/python2.7/site-packages/scipy/signal/signaltools.py",
line 1560, in detrend

    coef, resids, rank, s = linalg.lstsq(A, newdata[sl])

  File "/home/juan/anaconda/lib/python2.7/site-packages/scipy/linalg/basic.py",
line 550, in lstsq

    resids = np.sum(np.abs(x[n:])**2, axis=0)

MemoryError

I've made more attempts but now it takes 4000 seconds per day and sometimes it is
interrupted by reporting the same error above.

Some of the parameters used are: cc_sampling_rate: 20 Hz, ZZ correlation only, resampling
method: resample, filter parameters (Low: 0.01, Mwcs low: 0.12, Mwcs high: 0.98, high: 1,
Rms threshold: 0.0, Mwcs wlen:10 and Mwcs step: 5)

What can I do? Thank you very much in advance.
Laura.