) and I'll try to correct
it in the next release.
Brainstorming the issue:
- we could prevent MSNoise from adding (scanning) files that have a
sampling_rate lower than the `cc_sampling_rate`
- at `compute_cc` step, the data could be discarded.
This said, you should never have traces with different sampling rate
with the same seed id, e.g. "VG.LOK..HHZ" should always contain the same
sampling rate (at least during one day, I suppose). Actually, there
could be an issue in MSNoise now that I write these lines: if two
separate files are, e.g. BHZ and HHZ, MSNoise will try to read both
files and merge them. It's a baaaad idea !
I'll add this response to the issue mentionned above and find a solution
for next release !
Cheers
Thomas
Le 14/12/2015 09:03, Devy Kamil Syahbana a écrit :
Hello Thom,
Apa kabar?
MSNoise is great! I've succeeded to install it and happy to see the
results. But sometimes I got problems during the compute_cc process. It is
because some of the data have different sampling rate and/or gaps. These
force MSNoise stops. I could erase these "non-uniform data" manually in the
database and then reset CC and restart the compute_cc. But for a huge
amount of files, it is not convenient. Do you think is it possible to add
several lines in your MSNoise script so that when the data is not uniform,
MSNoise will skip it without stopping the process, but continue to compute
only the good data. Otherwise, I should fix my seismic data conversion. Or
do you have any better suggestions? Below is the error codes I sometime
found.
Traceback (most recent call last):
File "c:\anaconda\lib\runpy.py", line 162, in _run_module_as_main
"__main__", fname, loader, pkg_name)
File "c:\anaconda\lib\runpy.py", line 72, in _run_code
exec code in run_globals
File "C:\Anaconda\Scripts\msnoise.exe\__main__.py", line 9, in
<module>
File "c:\anaconda\lib\site-packages\msnoise\scripts\msnoise.py", line
393, in run
cli(obj={})
File "c:\anaconda\lib\site-packages\click-5.1-py2.7.egg\click\core.py",
line 700, in __call__
return self.main(*args, **kwargs)
File "c:\anaconda\lib\site-packages\click-5.1-py2.7.egg\click\core.py",
line 680, in main
rv = self.invoke(ctx)
File "c:\anaconda\lib\site-packages\click-5.1-py2.7.egg\click\core.py",
line 1027, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "c:\anaconda\lib\site-packages\click-5.1-py2.7.egg\click\core.py",
line 873, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "c:\anaconda\lib\site-packages\click-5.1-py2.7.egg\click\core.py",
line 508, in invoke
return callback(*args, **kwargs)
File "c:\anaconda\lib\site-packages\msnoise\scripts\msnoise.py", line
174, in compute_cc
main()
File "c:\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:\anaconda\lib\site-packages\msnoise\s03compute_cc.py", line 130,
in preprocess
stream[gap[0]] = stream[gap[0]].__add__(stream[gap[1]], method=0,
fill_value="interpolate")
File
"c:\anaconda\lib\site-packages\obspy-0.10.2-py2.7-win-amd64.egg\obspy\core\trace.py",
line 68
1, in __add__
raise TypeError("Sampling rate differs")
TypeError: Sampling rate differs
Cheers,
-----
Dr. Devy Kamil Syahbana
Head of Volcano Monitoring Section for Eastern Region of Indonesia
Ministry of Energy and Mineral Resources
Geological Agency
Center for Volcanology and Geological Hazard Mitigation
Jalan Diponegoro N°57
Bandung 40122
Indonesia
http://vsi.esdm.go.id
<https://www.avast.com/lp-safe-emailing?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
This
email has been sent from a virus-free computer protected by Avast.
www.avast.com
<https://www.avast.com/lp-safe-emailing?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>
<#DDB4FAA8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise