Problems with obspy arclink

coralie_aubert1 · August 23, 2013, 6:21am

Hi all,

I use obspy arclink library to download miniseed data from different networks. Until last week, I did not get any problems but since several days I get errors whereas I did not modify my program.

Here is my code. In the example below, I want to download data from the CH networf, AIGLE station.

megies · August 23, 2013, 12:16pm

Hi Coralie,

I can not say what exactly the problem with the routing is but if you
definitively did not change your ObsPy installation or your code then
there might be a problem with the routing table at webdc.eu (you could
ask at eida@gfz-potsdam.de).

If you need a quick fix you can just directly fetch the data from the
Swiss arclink node:

clientarclink = obspy.arclink.Client(host="eida.ethz.ch", port=18001)

best,
Tobias

P.S.: The building the timestamps does not work for me and it looks very
weird. Why not just UTCDateTime("2012-200T000000.0") or
UTCDateTime("%i-%03iT000000.0" % (year, doy)) if you need to put in
integers dynamically?

Amin_Khalil · August 23, 2013, 7:03pm

Dear Colleagues,
I have problems in reading data from nanometrics in both Y-format and their special seed archives. Is there any possibilities or library to read these data with ObsPy.

Thank you and best regards.

Amin

??? 23/08/2013 08:21 ?, ??? coralie aubert:

card2 - 1.jpg

barsch · August 23, 2013, 9:37pm

Dear Amin,

ObsPy has currently *no* support for Nanometrics Y-File Format -
however I looked into the format description at

and IMHO it seems fairly easy to implement a reader if someone could
provide us with small sample files and a e.g. MSEED export of the data
in order to write proper tests.

I will give it a try if you (or anyone) is willing to provide such
files

Cheers,
Robert

Dear Colleagues, I have problems in reading data from nanometrics
in both Y-format and their special seed archives. Is there any
possibilities or library to read these data with ObsPy.

Thank you and best regards.

Amin

على 23/08/2013 08:21 ص, كتب coralie aubert:

Hi all,

I use obspy arclink library to download miniseed data from
different networks. Until last week, I did not get any problems
but since several days I get errors whereas I did not modify my
program.

Here is my code. In the example below, I want to download data
from the CH networf, AIGLE station.
################################### from obspy.core import
UTCDateTime, Stream import obspy.arclink

clientarclink = obspy.arclink.Client() stT1 = str(2012) + 200 +
"T000000.0" t1 = UTCDateTime(stT1) stT2 = str(2012) + 200 +
"T235959.999999" t2 = UTCDateTime(stT2)

av = clientarclink.getInventory("CH", "AIGLE", '*', "HHZ", t1,
t2, restricted=False, permanent=None) st =
clientarclink.getWaveform("CH", "AIGLE",
av.keys()[0].split('.')[2], "HHZ", t1, t2)
#####################################

Here is the error I get after getInventory or getWaveform:

Traceback (most recent call last): File "<stdin>", line 1, in
<module> File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",
line 394, in getWaveform compressed=compressed, route=route)
File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",
line 514, in saveWaveform data = self._fetch(rtype, rdata,
route=route) File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",
line 200, in _fetch raise ArcLinkException(msg %
(request_data[2], request_data[3]))
obspy.arclink.client.ArcLinkException: Could not find route to
CH.AIGLE

I know this data exists because I already downloaded it. Has
something been modified in the code?

Thanks,

Regards, Coralie -- Coralie Aubert

Institut des Sciences de la Terre (ISTerre)

Tél. : 04 76 63 52 27 Fax : 04 76 63 52 52

web : isterre.fr/

_______________________________________________ obspy-users
mailing list obspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

_______________________________________________ obspy-users mailing
list obspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

- --
Dr. Robert Barsch

EGU Office Munich
Luisenstr. 37
80333 Munich
Germany

Phone: +49-89-21806549
Fax: +49-89-218017855
eMail: barsch@egu.eu

barsch · August 23, 2013, 9:44pm

oh I forgot about the SEED volumes - have a look at an old thread we
had on our mailing list

http://lists.swapbytes.de/archives/obspy-users/2011-July/000544.html

especially

http://lists.swapbytes.de/archives/obspy-users/2011-July/000549.html

the attachment there is a plain python source file - maybe this helps

Cheers,
Robert

Dear Colleagues, I have problems in reading data from nanometrics
in both Y-format and their special seed archives. Is there any
possibilities or library to read these data with ObsPy.

Thank you and best regards.

Amin

??? 23/08/2013 08:21 ?, ??? coralie aubert:

Hi all,

I use obspy arclink library to download miniseed data from
different networks. Until last week, I did not get any problems
but since several days I get errors whereas I did not modify my
program.

Here is my code. In the example below, I want to download data
from the CH networf, AIGLE station.
################################### from obspy.core import
UTCDateTime, Stream import obspy.arclink

clientarclink = obspy.arclink.Client() stT1 = str(2012) + 200 +
"T000000.0" t1 = UTCDateTime(stT1) stT2 = str(2012) + 200 +
"T235959.999999" t2 = UTCDateTime(stT2)

av = clientarclink.getInventory("CH", "AIGLE", '*', "HHZ", t1,
t2, restricted=False, permanent=None) st =
clientarclink.getWaveform("CH", "AIGLE",
av.keys()[0].split('.')[2], "HHZ", t1, t2)
#####################################

Here is the error I get after getInventory or getWaveform:

Traceback (most recent call last): File "<stdin>", line 1, in
<module> File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",
line 394, in getWaveform compressed=compressed, route=route)
File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",
line 514, in saveWaveform data = self._fetch(rtype, rdata,
route=route) File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",
line 200, in _fetch raise ArcLinkException(msg %
(request_data[2], request_data[3]))
obspy.arclink.client.ArcLinkException: Could not find route to
CH.AIGLE

I know this data exists because I already downloaded it. Has
something been modified in the code?

Thanks,

Regards, Coralie -- Coralie Aubert

Institut des Sciences de la Terre (ISTerre)

T�l. : 04 76 63 52 27 Fax : 04 76 63 52 52

web : isterre.fr/

_______________________________________________ obspy-users
mailing list obspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

_______________________________________________ obspy-users mailing
list obspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

- --
Dr. Robert Barsch

EGU Office Munich
Luisenstr. 37
80333 Munich
Germany

Phone: +49-89-21806549
Fax: +49-89-218017855
eMail: barsch@egu.eu

Amin_Khalil · August 24, 2013, 4:41am

Hi Robert,
Please find attached a sample Y-format file.

Thank you and best regards

Amin

على 23/08/2013 11:37 م, كتب Robert Barsch:

card2 - 1.jpg

YAYT_BHZ_20021223.124800 (71.2 KB)

barsch · August 24, 2013, 6:34pm

Hi Amin,

have a look at

github.com

obspy/obspy/blob/e267dc37401a7562af3592a2f1652b2ce31680ca/obspy/y/core.py

# -*- coding: utf-8 -*-
"""
Y bindings to ObsPy core module.

:copyright:
    The ObsPy Development Team (devs@obspy.org)
:license:
    GNU Lesser General Public License, Version 3
    (http://www.gnu.org/copyleft/lesser.html)
"""

from obspy import Stream
from obspy.core.trace import Trace
from obspy.core.utcdatetime import UTCDateTime
from obspy.core.util import AttribDict
from struct import unpack
import numpy as np


def __parseTag(fh):

This file has been truncated. show original

- - I guess it will take a while until it will be released - however you
may just take the script and run readY(filename) to retrieve a ObsPy
stream object. Please report any bugs at GitHub.

Also it would be very helpful to have a copy of your submitted file in
another format (MSEED, SAC, etc.) in order to compare the results. A Y
file written on a big endian system would complete the test suite.

Cheers,
Robert

Hi Robert, Please find attached a sample Y-format file.

Thank you and best regards

Amin

على 23/08/2013 11:37 م, كتب Robert Barsch: Dear Amin,

ObsPy has currently *no* support for Nanometrics Y-File Format -
however I looked into the format description at
http://www.ipgp.fr/~arnaudl/NanoCD/documentation/Nanometrics/Software/DataFormats/NmxDataFormats_RefGuide_14602R4.pdf

and IMHO it seems fairly easy to implement a reader if someone could

provide us with small sample files and a e.g. MSEED export of the
data in order to write proper tests.

I will give it a try if you (or anyone) is willing to provide such
files

Cheers, Robert

Dear Colleagues, I have problems in reading data from
nanometrics in both Y-format and their special seed archives.
Is there any possibilities or library to read these data with
ObsPy.

Thank you and best regards.

Amin

على 23/08/2013 08:21 ص, كتب coralie aubert:

Hi all,

I use obspy arclink library to download miniseed data from
different networks. Until last week, I did not get any
problems but since several days I get errors whereas I did
not modify my program.

Here is my code. In the example below, I want to download
data from the CH networf, AIGLE station.
################################### from obspy.core import
UTCDateTime, Stream import obspy.arclink

clientarclink = obspy.arclink.Client() stT1 = str(2012) +
200 + "T000000.0" t1 = UTCDateTime(stT1) stT2 = str(2012) +
200 + "T235959.999999" t2 = UTCDateTime(stT2)

av = clientarclink.getInventory("CH", "AIGLE", '*', "HHZ",
t1, t2, restricted=False, permanent=None) st =
clientarclink.getWaveform("CH", "AIGLE",
av.keys()[0].split('.')[2], "HHZ", t1, t2)
#####################################

Here is the error I get after getInventory or getWaveform:

Traceback (most recent call last): File "<stdin>", line 1,
in <module> File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",

line 394, in getWaveform compressed=compressed, route=route)

File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",

line 514, in saveWaveform data = self._fetch(rtype, rdata,

route=route) File
"/usr/local/lib/python2.7/site-packages/obspy-0.8.3-py2.7-linux-x86_64.egg/obspy/arclink/client.py",

line 200, in _fetch raise ArcLinkException(msg %

(request_data[2], request_data[3]))
obspy.arclink.client.ArcLinkException: Could not find route
to CH.AIGLE

I know this data exists because I already downloaded it.
Has something been modified in the code?

Thanks,

Regards, Coralie -- Coralie Aubert

Institut des Sciences de la Terre (ISTerre)

Tél. : 04 76 63 52 27 Fax : 04 76 63 52 52

web : isterre.fr/

_______________________________________________
obspy-users mailing listobspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

_______________________________________________ obspy-users
mailing listobspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

_______________________________________________ obspy-users
mailing list obspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

_______________________________________________ obspy-users mailing
list obspy-users@lists.swapbytes.de
http://lists.swapbytes.de/mailman/listinfo/obspy-users

- --
Dr. Robert Barsch

EGU Office Munich
Luisenstr. 37
80333 Munich
Germany

Phone: +49-89-21806549
Fax: +49-89-218017855
eMail: barsch@egu.eu

megies · August 24, 2013, 8:04pm

Hi obspy.arclink users,

the arclink node at GFZ which is used by default (i.e. when calling
Client() without arguments) has updated its ArcLink software. The update
introduces possible problems with the routing in between arclink nodes.

This will probably be fixed during a week or so at the server side but
if you experience problems fetching data in the next days or if you need
to be on the safe side, it might be advisable to use one of the other
entry nodes to the EIDA ArcLink network during the next week or so (see
the listing here:
http://docs.obspy.org/master/packages/autogen/obspy.arclink.client.Client.html?highlight=eida).

best,
Tobias

megies · September 6, 2013, 2:01pm

Hi obspy.arclink users,

the issue was fixed on the server side by the ArcLink protocol maintainers.

best,
Tobias