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
Hi Thomas:
I define my own data structure in the config, but when I run the command "msnoise populate" it says:
(base) PS C:\Users\b2-2> msnoise populate
>> Populating the Station table
Can't parse the archive for format NET/STA/YEAR/NET.STA.YEAR.DAY.SAC !
trying to import local parser (should return a station dictionary)
Traceback (most recent call last):
File "c:\users\b2-2\anaconda3\lib\site-packages\msnoise\s002populate_station_table.py", line 127, in main
from custom import populate
ModuleNotFoundError: No module named 'custom'
No file named custom.py in the C:\Users\b2-2 folder
I tried to add a custom.py in the path saied above, it still wrong!
Hello Thomas,
I am Anthony Abi Nader, PhD student at the laboratory Chrono-environnement, Besançon, France. I am using MSNoise to do autocorrelation on a seismic station located at 23m beneath the surface, my thesis' objective is to develop an approach in order to monitor karstic aquifers combining seismology, hydrogeology, and climatology.
I have 2 questions to ask you concerning MSNoise, the first one is while plotting dvv I only have plots for the median (in orange) and nothing for the mean (in blue); why is that? and the second question is how can I extract the dvv data table used in the plotting so I can use it elsewhere and superpose it with other data for example?
Thank you for your help.
Best regards,
Anthony Abi Nader
Hei everyone,
As I try to run "msnoise p tomo iftan" I receive the error message below.
The GUI starts but it does not work.
I am using macOS Catalina, Python 3.7 via Anaconda and mysql. Any
suggestions are highly appreciated, struggling with this for quite some
time.
Thanks,
Digdem
********************************************
fmin=0.1, fmax=1, nfreqs=40
bmin=0.0022, bmax=0.025
ampMin=0.05
vgMin=0.5 vgMax=5
disp: no
diag: Period - Vgroup
out: enable (matrix)
max: 0.1
write amplitude
write phase
write FT
write TV
/opt/anaconda3/envs/msnoisetomo/lib/python3.7/site-packages/msnoise_tomo/ftan_call.py:26:
UserWarning: loadtxt: Empty input file: "write_TV.txt"
V = np.loadtxt('write_TV.txt')
/opt/anaconda3/envs/msnoisetomo/lib/python3.7/site-packages/msnoise_tomo/ftan_call.py:28:
UserWarning: loadtxt: Empty input file: "write_amp.txt"
amp=np.loadtxt('write_amp.txt')
Exception in Tkinter callback
Traceback (most recent call last):
File "/opt/anaconda3/envs/msnoisetomo/lib/python3.7/tkinter/__init__.py",
line 1705, in __call__
return self.func(*args)
File
"/opt/anaconda3/envs/msnoisetomo/lib/python3.7/site-packages/msnoise_tomo/iftan.py",
line 103, in process
ydata)
File
"/opt/anaconda3/envs/msnoisetomo/lib/python3.7/site-packages/msnoise_tomo/ftan_call.py",
line 31, in pickgroupdispcurv
index = np.unravel_index(np.argmax(amp), amp.shape)
File
"/opt/anaconda3/envs/msnoisetomo/lib/python3.7/site-packages/numpy/core/fromnumeric.py",
line 1037, in argmax
return _wrapfunc(a, 'argmax', axis=axis, out=out)
File
"/opt/anaconda3/envs/msnoisetomo/lib/python3.7/site-packages/numpy/core/fromnumeric.py",
line 51, in _wrapfunc
return getattr(obj, method)(*args, **kwds)
ValueError: attempt to get argmax of an empty sequence
Hey everyone,
I'm trying to use the MSNoise-Tomo plugin to compute dispersion curves
that I can then use for a simple inversion, but I'm having a hard time
to find "good" parameters. Could someone write me what the Bmin, Bmax
parameters do and what they're called? Or refer me to a paper/article,
where I can read a bit more about the subject.
In the papers referenced I could find something regarding beamforming,
is that maybe the B parameter?
I appreciate gladly any help and references :)
Best regards,
Jannik
Hey everyone,
while working on my master's thesis I have to plot a bunch of dvv curves
computed with the stretching method, so I gave it a go and created a
plugin to do that. There's a command which works similar to the 'msnoise
plot dvv' one, but I also included some other stuff.
Maybe it could be useful for some of you too.
The github link + documentation is:
https://github.com/janfer95/ms-stretch
<https://github.com/janfer95/ms-stretch>
I tried to test it as much as I could, but I wouldn't be surprised if
there are still some bugs, so write me if there seems to be one. :)
Also a big thank you to the MSNoise Community (and Thomas Lecocq) for
making my master's thesis a whole lot easier.
Greetings,
Jannik
Hello,
I am using MsNoise to analyze seismic noise in Moscow, Idaho. I am just wondering if there are output tables for the dv/v plots. I can plot the figures and they look great, but I was just wondering how I can find the results table from which the dv/v plot is created. Thanks for the help!
-Quinn Buzzard
Hello, I am getting the following error when running 'msnoise test' after
following the installation documentation.
This is on Windows 10, using Python 3.7 via Anaconda.
Any suggestions as to how to fix welcome, thanks!
======================================================================
ERROR: test_025_dtt (msnoise.test.tests.MSNoiseTests)
----------------------------------------------------------------------
Traceback (most recent call last):
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\series.py",
line 1000, in __setitem__
self._set_with_engine(key, value)
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\series.py",
line 1033, 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
"d:\applications\anaconda\envs\msnoise\lib\site-packages\msnoise\test\tests.py",
line 312, in test_025_dtt
main()
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\msnoise\s06compute_dtt.py",
line 231, in main
df['err'][tmp] = 1.0
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\series.py",
line 1026, in __setitem__
self._set_with(key, value)
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\series.py",
line 1061, in _set_with
self.loc[key] = value
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\indexing.py",
line 666, in __setitem__
indexer = self._get_setitem_indexer(key)
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\indexing.py",
line 617, in _get_setitem_indexer
return self._convert_to_indexer(key, axis=0, is_setter=True)
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\indexing.py",
line 1198, in _convert_to_indexer
return self._get_listlike_indexer(key, axis, raise_missing=True)[1]
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\indexing.py",
line 1254, in _get_listlike_indexer
self._validate_read_indexer(keyarr, indexer, axis,
raise_missing=raise_missing)
File
"d:\applications\anaconda\envs\msnoise\lib\site-packages\pandas\core\indexing.py",
line 1304, in _validate_read_indexer
raise KeyError(f"{not_found} 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, 46] not in index'
----------------------------------------------------------------------
Ran 52 tests in 128.839s
FAILED (errors=1)
Bugreport:
************* Computer Report *************
----------------+PYTHON+-------------------
Python:3.7.8 | packaged by conda-forge | (default, Jul 31 2020, 01:53:57)
[MSC v.1916 64 bit (AMD64)]
This script is at
d:\applications\anaconda\envs\msnoise\lib\site-packages\msnoise\bugreport.py
---------------+MODULES+-------------------
Required:
[X] setuptools: 49.6.0.post20201009
[X] numpy: 1.19.2
[X] scipy: 1.2.1
[X] pandas: 1.1.3
[X] matplotlib: 3.3.2
[X] sqlalchemy: 1.3.19
[X] obspy: 1.2.2
[X] click: 7.1.2
[X] pymysql: 0.10.1
[X] flask: 1.1.2
[X] flask_admin: 1.5.6
[X] markdown: 3.3
[X] wtforms: 2.3.2
[X] folium: 0.11.0
[X] jinja2: 2.11.2
Only necessary if you plan to build the doc locally:
[ ] sphinx: not found
[ ] sphinx_bootstrap_theme: not found
Graphical Backends: (at least one is required)
[ ] wx: not found
[ ] pyqt: not found
[ ] PyQt4: not found
[X] PyQt5: present (no version)
[ ] PySide: not found
Not required, just checking:
[X] json: 2.0.9
[ ] psutil: not found
[ ] reportlab: not found
[ ] configobj: not found
[X] pkg_resources: present (no version)
[ ] paramiko: not found
[X] ctypes: 1.1.0
[X] pyparsing: 2.4.7
[X] distutils: 3.7.8
[ ] IPython: not found
[ ] vtk: not found
[ ] enable: not found
[ ] traitsui: not found
[ ] traits: not found
[ ] scikits.samplerate: not found
--
Kevin Davidson
Dear MSNoise users,
Although my cross-correlation functions (CCF) are relatively
symmetrical, they have different positive and negative lag-time (left
and right side). When I used the MSNoise for the stretching method, my
results for different dtt_sides (both sides, right side, and left side)
are exactly the same, while they should have different results. So, do
you know how can I fix this problem?
Sincerely,
Navid
Dear msnoise community,
I come back to msnoise after finishing my Master’s and I am currently trying to study velocity variations on Merapi volcano. I have used this useful package a couple of years ago and see that some parts of it have changed. I was able to follow the “simple run” based on the doc, up to the command msnoise stack -r, where I get the following error:
msnoise stack -r
Lets STACK !
Traceback (most recent call last):
File "/Users/rizal/miniconda/bin/msnoise", line 8, in <module>
sys.exit(run())
File "/Users/rizal/miniconda/lib/python3.7/site-packages/msnoise/scripts/msnoise.py", line 1202, in run
cli(obj={})
File "/Users/rizal/miniconda/lib/python3.7/site-packages/click/core.py", line 829, in __call__
return self.main(*args, **kwargs)
File "/Users/rizal/miniconda/lib/python3.7/site-packages/click/core.py", line 782, in main
rv = self.invoke(ctx)
File "/Users/rizal/miniconda/lib/python3.7/site-packages/click/core.py", line 1259, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/Users/rizal/miniconda/lib/python3.7/site-packages/click/core.py", line 1066, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/Users/rizal/miniconda/lib/python3.7/site-packages/click/core.py", line 610, in invoke
return callback(*args, **kwargs)
File "/Users/rizal/miniconda/lib/python3.7/site-packages/click/decorators.py", line 21, in new_func
return f(get_current_context(), *args, **kwargs)
File "/Users/rizal/miniconda/lib/python3.7/site-packages/msnoise/scripts/msnoise.py", line 732, in stack
main('ref', loglevel=loglevel)
File "/Users/rizal/miniconda/lib/python3.7/site-packages/msnoise/s04stack.py", line 172, in main
start, end, datelist = build_ref_datelist(db)
File "/Users/rizal/miniconda/lib/python3.7/site-packages/msnoise/api.py", line 1495, in build_ref_datelist
DataAvailability.starttime).first().starttime.date()
AttributeError: 'NoneType' object has no attribute ‘starttime'
I assumed that the source came from not defining startdate, but I have actually set startdate
msnoise config get startdate
M startdate: 2019-08-13
I’ve resetted the reference stack as well by
msnoise reset STACK
but the error persists. For information, I initially set a database in mysql but then set a new database using sqlite since it looks simpler.
Where do I start to troubleshoot this error? Thank you very much for your guidance.
Warm regards,
Mochammad Husni Rizal
Master in Earth Science - IPGP/Université de Paris
