[MSNoise] problems with different sampling rate and gap duration

[Thomas Lecocq]

Hello Devy,

Kabar baik !

Well, yes, it's (sadly) a known issue 
(https://github.com/ROBelgium/MSNoise/issues/51) 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 at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise