Robert Pickle
Robert Pickle
I think this is a great idea. I also once dreamt up an idea of something like `stream,inv = client.get(...)` which would just return both at once (ALL servers offer...
I think the fundamental issue is that we're still using evalresp which doesn't work with polynomials. I suppose previously these instruments were just written as a simple PAZ type response....
FWIW tried to AI (**warning!!!**) a full example stationXML based on the published FDSN 1.2 xsd.. had to fix a few items but looks pretty good, or at least decent...
> > FWIW tried to AI (**warning!!!**) a full example stationXML based on the published FDSN 1.2 xsd.. had to fix a few items but looks pretty good, or at...
this one looks super useful, will try to run some tests in the next week
Tested this on a 16 station multi-year network (7200 sec windows), worked great and it nearly halved the processing time. Nice!
here's the function in q: ```python def _generic_reader(pathname_or_url=None, callback_func=None, **kwargs): # convert pathlib.Path objects to str for compatibility. if isinstance(pathname_or_url, PurePath): pathname_or_url = str(pathname_or_url) if not isinstance(pathname_or_url, str): # not...
i can, just waiting for the board to clear a bit bytestring only-- good question, i assume in this context they will only arrive via XML but let me reconsider...
FWIW I think the xml you attached may not be completely valid FDSNXML. I am able to read it, but it's missing most of the proper stages in the file...
my apologies i think i misread what type of stage i was looking at! according to this: https://docs.obspy.org/_modules/obspy/core/inventory/response.html#ResponseListResponseStage the interpolation is `scipy.interpolate.InterpolatedUnivariateSpline` with a default k=3 ``` f = np.array([...