Hi,
Le 19/11/2015 07:41, Yoones Vaezi a écrit :
> Hi,
>
> Another problem I am facing is about the weights used for the MWCS
> program. The weights are supposed to be calculated using equation A5
> in Clarke et al. 2011 which is
> Inline image 2.
> But the following two lines in your Python code do not follow this
> equation:
>
> w = 1.0 / (1.0 / (coh[indRange]**2) - 1.0)
> w[coh[indRange] >= 0.99] = 1.0 / (1.0 / 0.9801 - 1.0)
>
>
> Could you please let me know where the values 0.99 and 0.9801 are
> coming from and what is the reason behind having these two lines in
> the code?
from the original code of Clarke, so, no, no idea precisely. I assume it
was to avoid infinite and zero division in the Fortran code.
Thomas
>
> Thank you very much for your time and consideration,
>
> Regards,
>
> Yoones
> __________________________________________
Hi,
Having looked at the MWCS.py program I have encountered a section which has
not been explained in the paper titled "MSNoise, a Python Package for
Monitoring Seismic Velocity Changes Using Ambient Seismic Noise" and I
could not find any justification why this has been done.
Within the program, before calculating the cross-spectrum between the
current and reference CCFs, the minimum values of each window is subtracted
and the result is divided by the maximum value. What is the reason for
doing this? In the paper it is just mentioned that the windows are
mean-adjusted and tapered.
I understand the reason behind detrending, mean adjustment, and the
tapering, but not the above-mentioned part of the code.
Thank you for your help,
Regards,
Yoones
______________________________
Dear MSNoise enthusiast,
I am currently processing data of Merapi volcano from 7 stations during
January 2013 - October 2015. When I type:
> msnoise compute_cc
I get:
2015-11-20 11:08:14 [INFO] *** Starting: Compute CC ***
> 2015-11-20 11:08:14 [INFO] Will compute ZZ
> 2015-11-20 11:08:14 [INFO] New CC Job: 2013-11-16 (1 pairs with 2 stations)
> Traceback (most recent call last):
> File "C:\Users\Dell\Anaconda\Scripts\msnoise-script.py", line 9, in
> <module>
> load_entry_point('msnoise==1.3.1', 'console_scripts', 'msnoise')()
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\msnoise\scripts\msnoise.py", line
> 393, in run
> cli(obj={})
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\click-4.0-py2.7.egg\click\core.py",
> line 664, in __call__
> return self.main(*args, **kwargs)
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\click-4.0-py2.7.egg\click\core.py",
> line 644, in main
> rv = self.invoke(ctx)
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\click-4.0-py2.7.egg\click\core.py",
> line 991, in invoke
> return _process_result(sub_ctx.command.invoke(sub_ctx))
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\click-4.0-py2.7.egg\click\core.py",
> line 837, in invoke
> return ctx.invoke(self.callback, **ctx.params)
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\click-4.0-py2.7.egg\click\core.py",
> line 464, in invoke
> return callback(*args, **kwargs)
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\msnoise\scripts\msnoise.py", line
> 174, in compute_cc
> main()
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\msnoise\s03compute_cc.py", line
> 271, in main
> basetime, tramef_Z = preprocess(db, stations, comps, goal_day, params,
> tramef_Z)
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\msnoise\s03compute_cc.py", line
> 148, in preprocess
> stream.merge(method=0, fill_value=0.0)
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\obspy-0.10.2-py2.7-win-amd64.egg\obspy\core\stream.py",
> line 1823, in merge
> interpolation_samples=interpolation_samples)
> File
> "C:\Users\Dell\Anaconda\lib\site-packages\obspy-0.10.2-py2.7-win-amd64.egg\obspy\core\trace.py",
> line 786, in __add__
> data = np.concatenate(data)
> MemoryError
Is this because the size of data is too large (~120 GB)?
Please help. Thank you.
Warm regards,
Mochammad Husni Rizal
B.Sc. (Geophysics) - Universitas Gadjah Mada
Phone : +62-857-151-22090
Dear msnoise community,
I have installed MSNoise and are really looking forward to use it. I
followed the documentation which is very helpful :). Unfortunately I
have a problem with scanning the archive. When I use
$ msnoise populate
I get:
>> Populating the Station table
Adding: ES BUENA
Adding: ES NATPA
That is correct. But the next step is then:
$ msnoise scan_archive
2015-11-17 18:09:51 [INFO] *** Starting: Scan Archive ***
2015-11-17 18:09:51 [INFO] Will work on 1 threads
2015-11-17 18:09:51 [INFO] Started:
/data/meivol/14-TEIDE/MINISEED/2014/ES/BUENA/HHZ.D
[INFO/Process-1] child process calling self.run()
[INFO/Process-1] process shutting down
[DEBUG/Process-1] running all "atexit" finalizers with priority >= 0
[DEBUG/Process-1] running the remaining "atexit" finalizers
[INFO/Process-1] process exiting with exitcode 0
2015-11-17 18:09:52 [INFO] Started:
/data/meivol/14-TEIDE/MINISEED/2014/ES/NATPA/HHZ.D
[INFO/Process-2] child process calling self.run()
[INFO/Process-2] process shutting down
[DEBUG/Process-2] running all "atexit" finalizers with priority >= 0
[DEBUG/Process-2] running the remaining "atexit" finalizers
[INFO/Process-2] process exiting with exitcode 0
2015-11-17 18:09:52 [INFO] *** Finished: Scan Archive ***
2015-11-17 18:09:52 [INFO] It took 1.00 seconds
[INFO/MainProcess] process shutting down
[DEBUG/MainProcess] running all "atexit" finalizers with priority >= 0
[DEBUG/MainProcess] running the remaining "atexit" finalizers
I attached the data_availability plot which shows that the data is
there... Is it maybe that I use the wrong data structure? I choose SDS
because my data looks like that:
/data/meivol/14-TEIDE/MINISEED/2014/ES/NATPA/HHZ.D/ES.NATPA..HHZ.D.2014.096
Thanks in advance,
Meike
--
Meike Volk
PhD student
School of Earth Sciences
University College Dublin
Hello Dylan,
Indeed, this looks like an incompatibility of types passed to resample
... You could maybre try to re-build the scikits-samplerate package
(uninstall and reinstall), that *could* work. Otherwise, either
downgrade numpy, or if you don't actually need the resampling (because
e.g. all your data is 100 Hz and you want to decimate to 20 Hz), you can
ignore it and use Decimate instead.
Could you open an issue on github for this ?
Thanks
Thomas
Le 13/11/2015 21:45, Dylan Mikesell a écrit :
> Hi MSNoise folks,
>
> I have a question about numpy. I installed cartopy today and during
> the install conda updated a bunch of other libraries, one of which was
> numpy. I had numpy 1.9.2 and updated to 1.10.0. Now when running
> 'compute_cc' I get the following error. I am not completely sure, but
> I believe it's related to samplerate and the input numpy data type.
> Has anyone ran into this error?
>
> Any thoughts on how to figure it out? I reinstalled scikits.samplrate,
> but that hasn't fixed the problem. And interestingly, the 'bugreport'
> says that samplerate is not found so perhaps this is the problem. But
> I haven't been able to figure it out yet and this error leads me to
> believe it's a problem with the numpy type.
> Thanks,
>
> Dylan
>
>
> [dmikesell@sonic deconData]$ msnoise compute_cc
>
> Traceback (most recent call last):
> File "/usr/local/anaconda/bin/msnoise", line 11, in <module>
> sys.exit(run())
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/msnoise/scripts/msnoise.py",
> line 393, in run
> cli(obj={})
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/click/core.py", line
> 700, in __call__
> return self.main(*args, **kwargs)
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/click/core.py", line
> 680, in main
> rv = self.invoke(ctx)
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/click/core.py", line
> 1027, in invoke
> return _process_result(sub_ctx.command.invoke(sub_ctx))
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/click/core.py", line
> 873, in invoke
> return ctx.invoke(self.callback, **ctx.params)
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/click/core.py", line
> 508, in invoke
> return callback(*args, **kwargs)
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/msnoise/scripts/msnoise.py",
> line 173, in compute_cc
> from ..s03compute_cc import main
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/msnoise/s03compute_cc.py",
> line 71, in <module>
> from scikits.samplerate import resample
> File
> "/usr/local/anaconda/lib/python2.7/site-packages/scikits.samplerate-0.3.3-py2.7-linux-x86_64.egg/scikits/samplerate/__init__.py",
> line 6, in <module>
> from _samplerate import resample, available_convertors,
> src_version_str, \
> File "numpy.pxd", line 43, in scikits.samplerate._samplerate
> (scikits/samplerate/_samplerate.c:3807)
> ValueError: numpy.dtype does not appear to be the correct type object
>
>
Dear Mr,
My name is Kike, currently I work in Chilean Volcano Observatory.
I am working with MSNoise 1.3.1 under Ubuntu 14.04.
After several steps:
Workflow
Installer
Configurator
Populate Station Table
Scan Archive
New Jobs
Compute Cross-Correlations -done and fine-
Stack (-r) -done-
Compute MWCS -done and fine-
Compute dt/t -done and fine-
...when I intended to plot, several graph have problems, especifically mwcs (wrong), distance(only color lines) and dvv (empty) : I attached the bug in file text and images.
Would I a help to you ?
I did the inscription yesterday and the processes was wrong.
Thanks so much for your orientation, KIKE
--
LUIS ENRIQUE FRANCO
Grupo de sismologia
Observatorio Volcanologico
Servicio de Geologia y Mineria
Temuco - Chile
Hi Angel,
Try to update traits and traitsui :
pip install traitsui==4.5.0
Thomas
Le 07/11/2015 23:12, msnoise-owner(a)mailman-as.oma.be a écrit :
> Dear msnoise community:
>
> I am installing msnoise in a 32-bits Ubuntu-Linux system without
> apparent errors.
>
> I have already created the database running the script msnoise
> install. But when I run the configurator script
>
> $ msnoise config
>
> I have an error like this:
>
> ImportError:
> /home/angel/anaconda/lib/python2.7/site-packages/traits/ctraits.so:
> undefined symbol: _PyTrash_thread_destroy_chain
>
> Can you help me with your experience?
>
> Thank you in advance
>
>
> Angel Figueroa-Soto
> Instituto de Investigaciones en Ciencias de la Tierra, UMSNH
> Research Institute for Earth Sciences, UMSNH
> Edif. "U" Ciudad Universitaria CP 58060
> Morelia, Michoacan. Mexico
> https://sites.google.com/site/angfsoto/
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.
Dear msnoise community:
I am installing msnoise in a 32-bits Ubuntu-Linux system without
apparent errors.
I have already created the database running the script msnoise install.
But when I run the configurator script
$ msnoise config
I have an error like this:
ImportError:
/home/angel/anaconda/lib/python2.7/site-packages/traits/ctraits.so:
undefined symbol: _PyTrash_thread_destroy_chain
Can you help me with your experience?
Thank you in advance
Angel Figueroa-Soto
Instituto de Investigaciones en Ciencias de la Tierra, UMSNH
Research Institute for Earth Sciences, UMSNH
Edif. "U" Ciudad Universitaria CP 58060
Morelia, Michoacan. Mexico
https://sites.google.com/site/angfsoto/
