Hi,
Since I have problems with storage memory on my machine, I was wondering if I can modify the path of the cross-correlation output in the script, and save the folder into another directory. Then modify the line of the next script to recall the cross-correlation file.
Thank you for your help!
Mari
Dear MSNoise developers.
First of all a big thank you for this wonderful piece of software!
Currently I am trying to compute a couple of CCF's from data acquired
via free-fall ocean bottom seismometers with hydrophones.
To facilitate the cc of the hydrophones (channel=H) I extented the list
of components in s03compute_no_rotation.py source-code file. Actually,
the hydrophones are more broadband than my geophone channels.
The computation works fine and when I dump the config.csv file the
paramter 'components_to_compute=HH' is in agreement with the folders and
files in ./STACKS.
However, when I try to plot the ccf's and spectra the plotting routines
wrongly assign the components as 'ZZ'.
I guess, have to expand another list components, but I was unable to
figure it out.
Thanks a lot for your help!
Best,
Florian
Hello everyone!
Hope everyone is doing well! I ran into an error while trying to run
compute_mwcs. I am not quite sure what the error fully means and was
wondering if anyone could help me with the error or point me in the right
direction to troubleshoot and resolve the error. Thanks!
[image: Screen Shot 2021-11-23 at 8.43.25 AM.png]
Best Wishes,
Michelle Lee
----------------------------------------------------------------------------
Ph.D. Candidate
Columbia University: Lamont-Doherty Earth Observatory
Department of Earth and Environmental Science
Marine Geology and Geophysics
c: (206)397-2488
alt. email: mkl2158(a)columbia.edu
Hello,
Hope everyone is doing well. I was wondering if there is any way to only
run a compute_mwcs and compute_dtt on just a single station pair. I tried
putting the station from "use" to "unuse" in the admin interface but that
didn't work. Thanks in advance.
Best Wishes,
Michelle Lee
----------------------------------------------------------------------------
Ph.D. Candidate
Columbia University: Lamont-Doherty Earth Observatory
Department of Earth and Environmental Science
Marine Geology and Geophysics
c: (206)397-2488
alt. email: mkl2158(a)columbia.edu
Hello,
Hope everyone is doing well. I am trying to plot dvv for individual station
pairs and am running into the problem of it coming up as an empty dataframe
(see attached screenshots), which in turn means nothing gets plotted for
the individual station pair and only the median and mean dvv gets plotted
(see attachment)
. I was wondering if anyone else had run into this issue before and knows
what to do or if there is something that I might have done wrong to cause
this issue. Thanks in advance! Have a great day![image: Screen Shot
2021-10-11 at 2.21.28 PM.png]
[image: plottingindvstatpair.png]
Best Wishes,
Michelle Lee
----------------------------------------------------------------------------
Graduate Student
Columbia University: Lamont-Doherty Earth Observatory
Department of Earth and Environmental Science
Marine Geology and Geophysics
c: (206)397-2488
alt. email: mkl2158(a)columbia.edu
Hello
I’m trying to install MSNoise, but keep getting unintelligible (to me) errors when running msnoise test:
*****
Traceback (most recent call last):
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/base.py", line 1276, in _execute_context
self.dialect.do_execute(
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/default.py", line 593, in do_execute
cursor.execute(statement, parameters)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/cursors.py", line 148, in execute
result = self._query(query)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/cursors.py", line 310, in _query
conn.query(q)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 548, in query
self._affected_rows = self._read_query_result(unbuffered=unbuffered)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 775, in _read_query_result
result.read()
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 1156, in read
first_packet = self.connection._read_packet()
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 725, in _read_packet
packet.raise_for_error()
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/protocol.py", line 221, in raise_for_error
err.raise_mysql_exception(self._data)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/err.py", line 143, in raise_mysql_exception
raise errorclass(errno, errval)
pymysql.err.ProgrammingError: (1146, "Table 'msnoise.config' doesn't exist")
The above exception was the direct cause of the following exception:
Traceback (most recent call last):
File "/opt/anaconda3/bin/msnoise", line 5, in <module>
from msnoise.scripts.msnoise import run
File "/opt/anaconda3/lib/python3.8/site-packages/msnoise/scripts/msnoise.py", line 1183, in <module>
plugins = get_config(db, "plugins")
File "/opt/anaconda3/lib/python3.8/site-packages/msnoise/api.py", line 215, in get_config
config = session.query(table).filter(table.name == name).first()
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/orm/query.py", line 3429, in first
ret = list(self[0:1])
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/orm/query.py", line 3203, in __getitem__
return list(res)
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/orm/query.py", line 3535, in __iter__
return self._execute_and_instances(context)
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/orm/query.py", line 3560, in _execute_and_instances
result = conn.execute(querycontext.statement, self._params)
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/base.py", line 1011, in execute
return meth(self, multiparams, params)
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/sql/elements.py", line 298, in _execute_on_connection
return connection._execute_clauseelement(self, multiparams, params)
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/base.py", line 1124, in _execute_clauseelement
ret = self._execute_context(
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/base.py", line 1316, in _execute_context
self._handle_dbapi_exception(
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/base.py", line 1510, in _handle_dbapi_exception
util.raise_(
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/util/compat.py", line 182, in raise_
raise exception
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/base.py", line 1276, in _execute_context
self.dialect.do_execute(
File "/opt/anaconda3/lib/python3.8/site-packages/sqlalchemy/engine/default.py", line 593, in do_execute
cursor.execute(statement, parameters)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/cursors.py", line 148, in execute
result = self._query(query)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/cursors.py", line 310, in _query
conn.query(q)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 548, in query
self._affected_rows = self._read_query_result(unbuffered=unbuffered)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 775, in _read_query_result
result.read()
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 1156, in read
first_packet = self.connection._read_packet()
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/connections.py", line 725, in _read_packet
packet.raise_for_error()
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/protocol.py", line 221, in raise_for_error
err.raise_mysql_exception(self._data)
File "/opt/anaconda3/lib/python3.8/site-packages/pymysql/err.py", line 143, in raise_mysql_exception
raise errorclass(errno, errval)
sqlalchemy.exc.ProgrammingError: (pymysql.err.ProgrammingError) (1146, "Table 'msnoise.config' doesn't exist")
[SQL: SELECT config.name AS config_name, config.value AS config_value
FROM config
WHERE config.name = %(name_1)s
LIMIT %(param_1)s]
[parameters: {'name_1': 'plugins', 'param_1': 1}]
(Background on this error at: http://sqlalche.me/e/13/f405)
*****
I’ve got all the required Python packages, MySQLServer and Workbench versions 8.0.18 (compatible with my macOS Mojave 10.14.6). Local database and user are set-up.
Any ideas?
Thanks,
Dave
David Healy [he/him]
School of Geosciences
Meston Walk
University of Aberdeen
Aberdeen AB24 3UE
My working day may not be your working day. Please do not feel obliged to reply to this e-mail outside of your normal working hours.
Publications: https://scholar.google.com/citations?hl=en&user=ug9DnVoAAAAJ
Software: http://www.fracturedrockphysics.com/software.html
Home page: http://www.fracturedrockphysics.com/index.html
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.
Dear Thomas,
When I run MSnoise 1.6 with msnoise test command, the error is as follows.
I'd like to kown what the problem is.
thank you very much, if you could reply to me.
======================================================================
ERROR: test_025_dtt (msnoise.test.tests.MSNoiseTests)
----------------------------------------------------------------------
Traceback (most recent call last):
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 972, in __setitem__
self._set_with_engine(key, value)
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 1005, in _set_with_engine
loc = self.index._engine.get_loc(key)
File "pandas\_libs\index.pyx", line 70, in pandas._libs.index.IndexEngine.get_loc
File "pandas\_libs\index.pyx", line 75, in pandas._libs.index.IndexEngine.get_loc
TypeError: '[ 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 23 31 32 33 34 35 36 37 38
39 40 41 42 43 44 45 46]' is an invalid key
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "c:\anaconda3\lib\site-packages\msnoise\test\tests.py", line 312, in test_025_dtt
main()
File "c:\anaconda3\lib\site-packages\msnoise\s06compute_dtt.py", line 231, in main
df['err'][tmp] = 1.0
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 998, in __setitem__
self._set_with(key, value)
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 1031, in _set_with
self._set_labels(key, value)
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 1042, in _set_labels
raise KeyError(f"{key[mask]} not in index")
KeyError: '[ 1 2 3 4 6 7 8 9 11 12 13 14 23 31 32 33 34 36 37 38 39 41 42 43 44\n 46] not in index'
----------------------------------------------------------------------
Ran 52 tests in 147.687s
祝 好
江苏省地震监测中心 王 俊
------------------ Original ------------------
From: "msnoise" <msnoise-request(a)mailman-as.oma.be>gt;;
Date: Wed, Apr 21, 2021 08:00 PM
To: "msnoise"<msnoise(a)mailman-as.oma.be>gt;;
Subject: MSNoise Digest, Vol 75, Issue 3
Send MSNoise mailing list submissions to
msnoise(a)mailman-as.oma.be
To subscribe or unsubscribe via the World Wide Web, visit
http://mailman-as.oma.be/mailman/listinfo/msnoise
or, via email, send a message with subject or body 'help' to
msnoise-request(a)mailman-as.oma.be
You can reach the person managing the list at
msnoise-owner(a)mailman-as.oma.be
When replying, please edit your Subject line so it is more specific
than "Re: Contents of MSNoise digest..."
Today's Topics:
1. Re: Issue-regarding Grid.dat (Dr. Arun Singh)
----------------------------------------------------------------------
Message: 1
Date: Tue, 20 Apr 2021 21:25:36 +0530 (IST)
From: "Dr. Arun Singh" <arun(a)gg.iitkgp.ac.in>
To: msnoise <msnoise(a)mailman-as.oma.be>
Subject: Re: [MSNoise] Issue-regarding Grid.dat
Message-ID:
<1913349746.354648.1618934136707.JavaMail.zimbra(a)gg.iitkgp.ac.in>
Content-Type: text/plain; charset=utf-8
Dear Thomas,
I think, I did a mistake. Thank you for pointing it out.
arun
----- Original Message -----
From: "Dr. Arun Singh" <arun(a)gg.iitkgp.ac.in>
To: "msnoise" <msnoise(a)mailman-as.oma.be>
Sent: Tuesday, April 20, 2021 6:44:42 PM
Subject: Re: [MSNoise] Issue-regarding Grid.dat
Hi Thomas,
Yes I think they are correct? Can you help from where is it taking these values? I can cross check them. If it is
----- Original Message -----
From: "Thomas Lecocq" <Thomas.Lecocq(a)seismology.be>
To: "msnoise" <msnoise(a)mailman-as.oma.be>
Sent: Tuesday, April 20, 2021 4:25:04 PM
Subject: Re: [MSNoise] Issue-regarding Grid.dat
Hi,
Are you sure your X & Y values are OK in the database ? The code
https://github.com/ThomasLecocq/msnoise-tomo/blob/master/msnoise_tomo/prepa…
is correct: element 2 is Latitude (Y) and element 3 is Longitude (X)
Thomas
Le 20/04/2021 ? 12:04, Dr. Arun Singh a ?crit?:
> Dear All,
> Can any one of you help me out.
>
> prepare_tomo.py is writing "Grid.dat", However, it seems it is writing wrong minlat and minlon. In place of longitude it seems it is writing latitude. Station coordinates are correct in STACoord.dat. How X and
>
> minlon = df[3].min() - xstep*1.5
> maxlon = df[3].max() + xstep*1.5
> minlat = df[2].min() - ystep*1.5
> maxlat = df[2].max() + ystep*1.5
> f.write("%f %f\n"%(minlon, maxlon)) # It is writing latitude
> f.write("%f %f\n"%(minlat, maxlat)) # It is writing longitude
> f.write("%f %f\n"%(xstep, ystep))
> f.close
>
> Somehow in the final tomographic plots (result_*), I am getting a a mirror image of expected image. (Latitude in x-axis and Longitude in Y-axis). If I change Grid.dat manually (by changing lat and lon) before running ANSWT.py, it is not working well.
>
> Regards,
>
> arun
> _______________________________________________
> 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
------------------------------
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
End of MSNoise Digest, Vol 75, Issue 3
**************************************
Dear Thomas,
When I run MSnoise 1.6 with msnoise test command, the error is as follows.
I'd like to kown what the problem is.
thank you very much, if you could reply to me.
======================================================================
ERROR: test_025_dtt (msnoise.test.tests.MSNoiseTests)
----------------------------------------------------------------------
Traceback (most recent call last):
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 972, in __setitem__
self._set_with_engine(key, value)
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 1005, in _set_with_engine
loc = self.index._engine.get_loc(key)
File "pandas\_libs\index.pyx", line 70, in pandas._libs.index.IndexEngine.get_loc
File "pandas\_libs\index.pyx", line 75, in pandas._libs.index.IndexEngine.get_loc
TypeError: '[ 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 23 31 32 33 34 35 36 37 38
39 40 41 42 43 44 45 46]' is an invalid key
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "c:\anaconda3\lib\site-packages\msnoise\test\tests.py", line 312, in test_025_dtt
main()
File "c:\anaconda3\lib\site-packages\msnoise\s06compute_dtt.py", line 231, in main
df['err'][tmp] = 1.0
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 998, in __setitem__
self._set_with(key, value)
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 1031, in _set_with
self._set_labels(key, value)
File "c:\anaconda3\lib\site-packages\pandas\core\series.py", line 1042, in _set_labels
raise KeyError(f"{key[mask]} not in index")
KeyError: '[ 1 2 3 4 6 7 8 9 11 12 13 14 23 31 32 33 34 36 37 38 39 41 42 43 44\n 46] not in index'
----------------------------------------------------------------------
Ran 52 tests in 147.687s
祝 好
江苏省地震监测中心 王 俊
------------------ Original ------------------
From: "msnoise" <msnoise-request(a)mailman-as.oma.be>gt;;
Date: Wed, Apr 21, 2021 08:00 PM
To: "msnoise"<msnoise(a)mailman-as.oma.be>gt;;
Subject: MSNoise Digest, Vol 75, Issue 3
Send MSNoise mailing list submissions to
msnoise(a)mailman-as.oma.be
To subscribe or unsubscribe via the World Wide Web, visit
http://mailman-as.oma.be/mailman/listinfo/msnoise
or, via email, send a message with subject or body 'help' to
msnoise-request(a)mailman-as.oma.be
You can reach the person managing the list at
msnoise-owner(a)mailman-as.oma.be
When replying, please edit your Subject line so it is more specific
than "Re: Contents of MSNoise digest..."
Today's Topics:
1. Re: Issue-regarding Grid.dat (Dr. Arun Singh)
----------------------------------------------------------------------
Message: 1
Date: Tue, 20 Apr 2021 21:25:36 +0530 (IST)
From: "Dr. Arun Singh" <arun(a)gg.iitkgp.ac.in>
To: msnoise <msnoise(a)mailman-as.oma.be>
Subject: Re: [MSNoise] Issue-regarding Grid.dat
Message-ID:
<1913349746.354648.1618934136707.JavaMail.zimbra(a)gg.iitkgp.ac.in>
Content-Type: text/plain; charset=utf-8
Dear Thomas,
I think, I did a mistake. Thank you for pointing it out.
arun
----- Original Message -----
From: "Dr. Arun Singh" <arun(a)gg.iitkgp.ac.in>
To: "msnoise" <msnoise(a)mailman-as.oma.be>
Sent: Tuesday, April 20, 2021 6:44:42 PM
Subject: Re: [MSNoise] Issue-regarding Grid.dat
Hi Thomas,
Yes I think they are correct? Can you help from where is it taking these values? I can cross check them. If it is
----- Original Message -----
From: "Thomas Lecocq" <Thomas.Lecocq(a)seismology.be>
To: "msnoise" <msnoise(a)mailman-as.oma.be>
Sent: Tuesday, April 20, 2021 4:25:04 PM
Subject: Re: [MSNoise] Issue-regarding Grid.dat
Hi,
Are you sure your X & Y values are OK in the database ? The code
https://github.com/ThomasLecocq/msnoise-tomo/blob/master/msnoise_tomo/prepa…
is correct: element 2 is Latitude (Y) and element 3 is Longitude (X)
Thomas
Le 20/04/2021 ? 12:04, Dr. Arun Singh a ?crit?:
> Dear All,
> Can any one of you help me out.
>
> prepare_tomo.py is writing "Grid.dat", However, it seems it is writing wrong minlat and minlon. In place of longitude it seems it is writing latitude. Station coordinates are correct in STACoord.dat. How X and
>
> minlon = df[3].min() - xstep*1.5
> maxlon = df[3].max() + xstep*1.5
> minlat = df[2].min() - ystep*1.5
> maxlat = df[2].max() + ystep*1.5
> f.write("%f %f\n"%(minlon, maxlon)) # It is writing latitude
> f.write("%f %f\n"%(minlat, maxlat)) # It is writing longitude
> f.write("%f %f\n"%(xstep, ystep))
> f.close
>
> Somehow in the final tomographic plots (result_*), I am getting a a mirror image of expected image. (Latitude in x-axis and Longitude in Y-axis). If I change Grid.dat manually (by changing lat and lon) before running ANSWT.py, it is not working well.
>
> Regards,
>
> arun
> _______________________________________________
> 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
------------------------------
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
End of MSNoise Digest, Vol 75, Issue 3
**************************************
Good-day,
I am currently using msnoise to work on the recent St. Vincent eruption.
However, when I try to compute_cc, I am unable to get it to remove
instrument response. I am using python version 3.7 and MySQL version 8.0.
I have ensured that the dataless files are matching the data structure as
the mseed files and that the dataless files are located as instructed. I
have also set a direct path to the files in the configuration but I keep
getting the same error.
If I set the instrument response to N, all the other stages work perfectly
with no errors.
I have attached the errors I got as well as my msnoise configuration page.
Thanks in advance for any assistance and please let me know if you require
further information.
Much thanks,
Yadava Persad
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…>
Virus-free.
www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…>
<#m_-3219373662840979284_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>