Using MSNoise 1.24 there are two python scripts that plot; they are
07plot_dtt.py and s07plot_dtt.py. I get errors for either.
Using s07plot_dtt.py yields attachment s07plot_dtt.png. There are two
problems:
1) I entered mov_stack values 1,3,5, and 10 days into the configurator but
the graphical output includes only 1, 3, and 5 days moving window graphs,
and the only graph with any data plotted onto it is the 5 days moving
window.
2) the graph x axes automatically include a much larger time period than
covered by the data plotted. This includes dates back to 2013 and forward
into 2016, while all the data I have is from 2014. I am unsure if this
problem is related to some of my settings in the configurator. My settings
that I can see might be related are: enddate=2100-01-01, ref_begin=-100,
ref_end=0, startdate=2014-05-04.
Using 07plot_dtt.py yields attachment image 07.plot_dtt.png the following:
>>> runfile('/home/devora/MSNoise-master/07.plot_dtt.py',
wdir=r'/home/devora/MSNoise-master')
UMD has deleted: database_tools, msnoise_table_def
loading 1 days
loading 3 days
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File
"/usr/lib/python2.7/site-packages/spyderlib/widgets/externalshell/sitecustomize.py",
line 540, in runfile
execfile(filename, namespace)
File "/home/devora/MSNoise-master/07.plot_dtt.py", line 128, in <module>
plt.fill_between(ALL.index,ALL[dttname]-ALL[errname],ALL[dttname]+ALL[errname],lw=1,color='red',zorder=-1,alpha=0.3)
File "/usr/lib64/python2.7/site-packages/matplotlib/pyplot.py", line 2757,
in fill_between
interpolate=interpolate, **kwargs)
File "/usr/lib64/python2.7/site-packages/matplotlib/axes.py", line 6988, in
fill_between
x = ma.masked_invalid(self.convert_xunits(x))
File "/usr/lib64/python2.7/site-packages/numpy/ma/core.py", line 2239, in
masked_invalid
condition = ~(np.isfinite(a))
TypeError: ufunc 'isfinite' not supported for the input types, and the
inputs could not be safely coerced to any supported types according to the
casting rule ''safe''
>>>
Has anyone else encountered these problems before? Does anyone have any
suggestions for what I might do to overcome them?
Thankyou
Josiah
Dear Thomas, dear MSNoise users,
No problem to fix, just a tip.
For people who are not familiar at all with MySQL (like me), here is how to
reset the correlation computation jobs from 'D'one to 'T'odo between a
certain range of dates. I needed to do that because the disk where I wrote
the correlations get full, but MSNoise kept trying to do the following job
even if the correlation couldn't be written. Actually, maybe there is a way
to let the job set to 'T' if the correlation is not really saved on a
disk...
Anyway, here is the command I used in phpMyAdmin, using the SQL tab, once
you selected your database:
UPDATE `job` SET `flag` = 'T' WHERE `day` BETWEEN 'date1' AND 'date2';
Then press Go button.
Before using it, I would like to have the advice from a MySQL expert: is it
the right way to do it? If yes, it could be worth to have this line in the
documentation somewhere.
Cheers,
Aurelien.
--
aurelien.mordret(a)gmail.com
Mobile: 0033 6 67 51 53 56
Dear all,
When I run Compute_cc the terminal show this message:
File "s03compute_cc.py", line 313, in <module>
X0 = s1.X
AttributeError: 'NoneType' object has no attribute 'X'
?In which way can I fix it?? The attribute X, Y and also Z exits in the table and I also populate it. So I don't understand this type of error. thank you very much. I hope to receive some answer..
cheers,
Carmelo
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,
I have just updated to 1.2.5 and tried to run *s07plot_dtt.py *(again on
both PCs) and I got this error message:
oscar@terminus:~/msnoise_sismo$ py s07plot_dtt.py
Traceback (most recent call last):
File "s07plot_dtt.py", line 13, in <module>
import matplotlib.pyplot as plt
File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/matplotlib/pyplot.py",
line 98, in <module>
_backend_mod, new_figure_manager, draw_if_interactive, _show =
pylab_setup()
File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/matplotlib/backends/__init__.py",
line 28, in pylab_setup
globals(),locals(),[backend_name],0)
File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/matplotlib/backends/backend_wxagg.py",
line 6, in <module>
import backend_wx # already uses wxversion.ensureMinimal('2.8')
File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/matplotlib/backends/backend_wx.py",
line 53, in <module>
raise ImportError(missingwxversion)
ImportError: Matplotlib backend_wx and backend_wxagg require wxversion,
which was not found.
Please help!!!!!
Oscar
Hi Achmad,
1/ Did you follow the steps of the workflow ?
2/ What is your archive folder structure / format ?
3/ for the filter problem, consider reading this thread of the ML :
http://news.gmane.org/gmane.science.geophysics.msnoise --> "Problem with
the filters tab"
Thomas
ps: your message was not forwarded to the ML because of attached images,
I have to find a way to fix this.
Le 09/10/2014 05:56, achmad gunar saadi a écrit :
> Good morning,
> I am Achmad Gunar Saadi, Gephisics student in Universitas Gadjah Mada.
> I wanna ask about menu in the configurator, why in the Filters and
> Networks tab doesn't appear any option?. And also about setting
> data_structure, How to fix this problems?. I have watched your video
> tutorial "Msnoise first install", but I still get those problems.
> Herewith I sent also the screenshot.
> Thank you very much.
>
> Yours sincerely,
>
> Achmad Gunar Saadi
> Universitas Gadjah Mada
> Inline image 1Inline image 2
Hello,
I am using MSNoise on two different PC's (Ubuntu 12.04 [AMD A6] and Ubuntu
14.04 [Intel i7]), each computer has the same version of MSNoise (1.2.4)
but differents sets of data. I compute all the steps without a problem but
when I try to compute *07.plot_dtt.py <http://07.plot_dtt.py>* the routine
fails. I got two different error messages:
Ubuntu 12.04:
*oscar@bayta:~/sweet_noise$ py 07.plot_dtt.py <http://07.plot_dtt.py> *
*/home/oscar/anaconda/lib/python2.7/site-packages/numpy/oldnumeric/__init__.py:11:
ModuleDeprecationWarning: The oldnumeric module will be dropped in Numpy
1.9*
* warnings.warn(_msg, ModuleDeprecationWarning)*
*/home/oscar/anaconda/lib/python2.7/site-packages/setuptools-2.2-py2.7.egg/pkg_resources.py:991:
UserWarning: /home/oscar/.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).*
*loading 1 days*
*Traceback (most recent call last):*
* File "07.plot_dtt.py <http://07.plot_dtt.py>", line 115, in <module>*
* py1_wmean, py1_wstd = get_wavgwstd(allbut)*
* File "07.plot_dtt.py <http://07.plot_dtt.py>", line 59, in get_wavgwstd*
* h = grouped.apply(wstd)*
* File
"/home/oscar/anaconda/lib/python2.7/site-packages/pandas/core/groupby.py",
line 424, in apply*
* return self._python_apply_general(f)*
* File
"/home/oscar/anaconda/lib/python2.7/site-packages/pandas/core/groupby.py",
line 427, in _python_apply_general*
* keys, values, mutated = self.grouper.apply(f, self.obj, self.axis)*
* File
"/home/oscar/anaconda/lib/python2.7/site-packages/pandas/core/groupby.py",
line 883, in apply*
* res = f(group)*
* File
"/home/oscar/anaconda/lib/python2.7/site-packages/pandas/core/groupby.py",
line 422, in f*
* return func(g, *args, **kwargs)*
* File "07.plot_dtt.py <http://07.plot_dtt.py>", line 52, in wstd*
* wstd = np.sqrt(np.sum(w * (d - wavg) ** 2) / ((N - 1) * np.sum(w) /
N))*
*ZeroDivisionError: float division by zero*
Ubuntu 14.04:
*oscar@terminus:~/msnoise_sismo$ py 07.plot_dtt.py <http://07.plot_dtt.py>*
*/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/setuptools-3.6-py2.7.egg/pkg_resources.py:1045:
UserWarning: /home/oscar/.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).*
*loading 1 days*
*Traceback (most recent call last):*
* File "07.plot_dtt.py <http://07.plot_dtt.py>", line 128, in <module>*
*
plt.fill_between(ALL.index,ALL[dttname]-ALL[errname],ALL[dttname]+ALL[errname],lw=1,color='red',zorder=-1,alpha=0.3)*
* File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/matplotlib/pyplot.py",
line 2757, in fill_between*
* interpolate=interpolate, **kwargs)*
* File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/matplotlib/axes.py",
line 6988, in fill_between*
* x = ma.masked_invalid(self.convert_xunits(x))*
* File
"/home/oscar/bin/python/anaconda/lib/python2.7/site-packages/numpy/ma/core.py",
line 2239, in masked_invalid*
* condition = ~(np.isfinite(a))*
*TypeError: ufunc 'isfinite' not supported for the input types, and the
inputs could not be safely coerced to any supported types according to the
casting rule ''safe''*
I am not a python expert and I haven't could solve these problems. If
anyone knows how to solve it please help me!
Thanks in advance!
Oscar
Hi Carmelo,
Tom suggested the issue may be linux related, and proposed some hard-coding
of the filters, but I have not had any problems on two linux installs (both
fedora and ubuntu) defining the filters. They are defined in the
configurator:
http://msnoise.org/doc/workflow.html#module-MSNoise.s001configurator
Do you see the filter tab on the left? This should open up boxes in the gui
where you can set the filters according to your desires. If your system
behaves differently, please describe in detail.
Cheers,
Kasper
On 4 October 2014 10:05, Sammarco, Carmelo <carmelosammarco(a)abdn.ac.uk>
wrote:
> Dear Thomas, Aurelien and MsNoise users,
>
> When I run s03compute_cc.py the terminal shows: 'NO FILTERS DEFINED,
> exiting'. If I go in the configurator and in the filter Tab the section is
> blank. In which way can I populate the filter tab?? Someone could explain
> me in which way can I fix the problem?? I am using linux and Mysql. Can I
> add the filters information on the database with phpmyadmin? And if yes in
> which way? I am a PhD student of the university of Aberdeen and I need to
> fix this problem as soon as possible. Please someone help me please. I need
> to do a lot of job so I can't rerun like Aurelien says:
>
> "IF you rerun s03compute_cc.py after this error, it will jump to the next
> job. It's a dirty quick fix and it does not help if the error happens every
> two days of data..."
>
>
> Thank you very much for the info!! Better than nothing.. I appreciate a
> lot but I am looking if possible for a permanent solution.
>
> I hope to recieve your replay.
>
> Cheers,
>
> Carmelo
>
>
>
>
>
>
> 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.
> _______________________________________________
> MSNoise mailing list
> MSNoise(a)mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
>
Tom, as you suggested, the problem may be bigger than file names/ networks.
Here's the output after changing to the NZ network name of our extra
station to match that of the rest of the stations. Both are mseed data
types:
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File
"/usr/lib/python2.7/site-packages/spyderlib/widgets/externalshell/sitecustomize.py",
line 540, in runfile
execfile(filename, namespace)
File "/home/devora/Test/s03compute_cc.py", line 214, in <module>
stream.merge(fill_value=0)
File "/usr/lib64/python2.7/site-packages/obspy/core/stream.py", line 1715,
in merge
self._mergeChecks()
File "/usr/lib64/python2.7/site-packages/obspy/core/stream.py", line 1660,
in _mergeChecks
raise Exception(msg)
Exception: Can't merge traces with same ids but differing data types!
On 3 October 2014 20:27, Thomas Lecocq <thomas.lecocq(a)oma.be> wrote:
> Dear Kasper,
>
> Le 02/10/2014 23:38, Kasper van Wijk a écrit :
> > Dear Tom et al.,
> >
> > I am working with a student to replicate the install of MSNOISE for the
> > Auckland Volcanic Field. Below are a few things we have ran into:
> >
> > 1. http://msnoise.org/doc/workflow.html#module-MSNoise.s001configurator
> > misses the explicit call to start the configurator: python
> > s001configurator.py
> >
> > (this is obvious to anyone who has used MSNOISE before, but stumped the
> > aforementioned new student)
> I've updated the documentation, thanks !
> >
> > 2. On CENTOS linux, I have an issue with sqlalchemy. The "fix" is as ugly
> > as it is mysterious:
> >
> >
> http://stackoverflow.com/questions/14600530/cant-import-from-sqlalchemy-ext…
> Well, this looks strange indeed, but good to know, I'll make a new
> "Troobleshooting" section in the doc !
> >
> > 3. We have a font problem (which I don't know how to fix), but starting
> the
> > configurator, we get the attached warning. Just a warning, though.
> Never seen this before... probably somehow specific to paranoid-CentOS
> ... :)
> >
> > 4. Is it correct that sqlite can only work on one thread?
> >
> > [devora@localhost MSNoise-master]$ python s01scan_archive.py --init -t 2
> > 2014-08-22 10:14:13,724 [INFO] *** Starting: Scan Archive ***
> > 2014-08-22 10:14:13,727 [INFO] Initializing (should be run only once)
> > 2014-08-22 10:14:13,727 [INFO] You can not work on 2 threads because
> SQLite
> > only supports 1 connection at a time
> > 2014-08-22 10:14:13,727 [INFO] Will work on 1 threads
>
> Well, in "theory", it could work on more than one thread, but there is a
> very big risk of problems, because SQLite can only accept one
> transaction at a time, so if two process try to insert, update (...) in
> the DB at the same time, you'll get a "Cannot insert, database locked"
> and the data will be lost...
>
> >
> > 5. And this is where are currently: we imported a station into our NZ
> > network with a different network name. This resulted in errors when
> pairing
> > the AK station with an NZ station. What should we do? After getting an
> > error, we tried renaming the AK data to NZ to be consistent, but we are
> > unable to cleanly perge the AK info in the database. Alternatively, can
> we
> > set up two networks in the configurator?
>
> Are you trying to run a second "project", independent of the AVF ? Or
> using a NZ station for the AVF monitoring ? I don't get this.
>
> Renaming station files is usually not enough, if you don't rename
> metadata within the files + update the database... I'd say it's usually
> not worth the effort, MSNoise should be able to process data from
> different networks... What is the error you get ?
>
> Best regards and thanks for using MSNoise :)
>
> Tom
> >
> > Cheers,
> > kasper
> >
> >
> > _______________________________________________
> > MSNoise mailing list
> > MSNoise(a)mailman-as.oma.be
> > http://mailman-as.oma.be/mailman/listinfo/msnoise
>
> _______________________________________________
> MSNoise mailing list
> MSNoise(a)mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
>
Dear all,
Someone can explain me in which way I can populate in the configurator the filters tab?? When I try to start cross correlation give me this error message: 'NO FILTERS DEFINED, exiting'. what it means? thank you very much and I hope to find a solution as soon as possible.
cheers,
Carmelo
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.