[obspy-users] Dip and azimuth information using arclink client

Januka Attanayake jattanayake at yahoo.com
Fri May 31 07:40:51 CEST 2019


Hi All, I am following this tutorial of yours to calculate PPSDs (https://docs.obspy.org/tutorial/code_snippets/probabilistic_power_spectral_density.html) and getting the following error. How do I fix it? Thanks! 
print(ppsd.times[:2])

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-17-cbfe003c07f0> in <module>
----> 1 print(ppsd.times[:2])

AttributeError: 'PPSD' object has no attribute 'times'


Januka Attanayake
Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake

School of Earth Sciences | McCoy Bldg. 200 
University of Melbourne | Parkville 3010 VIC 
Australia
 

    On Tuesday, April 16, 2019, 10:35:32 PM GMT+10, Lion Krischer <lion.krischer at erdw.ethz.ch> wrote:  
 
 Hi Coralie,

sorry for the late reply.

You can get dataless SEED files from ArcLink with:

>>> client.save_response(filename="dataless.seed", network="BW",
station="MANZ", location="", channel="EHZ", starttime=t, endtime=t+10)

You can then read the resulting file with ObsPy and access the channel
level information:

>>> obspy.read_inventory("dataless.seed")[0][0][0]
Channel 'EHZ', Location ''
    Time range: 2006-01-01T00:00:00.000000Z - 2015-10-21T12:00:00.000000Z
    Latitude: 49.99, Longitude: 12.11, Elevation: 635.0 m, Local Depth:
0.0 m
    Azimuth: 0.00 degrees from north, clockwise
    Dip: -90.00 degrees down from horizontal
    Sampling Rate: 200.00 Hz
    Sensor (Description): STS-2 (None)
    Response information available

Hope it helps,

Lion


On 12.03.19 11:44, Coralie Aubert wrote:
> Hello,
>
> Does anybody know how to get the dip and azimuth information when using
> arclink client from obspy?
>
> I do not find this information in the attached response and I do not
> know how to get it.
>
> Here is my program:
>
> from obspy.clients.arclink.client import Client as arclink_client
>
> client = arclink_client(user=UserID)
> st =
> client.get_waveforms(net,station,locid,comp,T1,T2,attach_response=True)
>
> Regards,
> Coralie
>
> -- 
> Coralie Aubert
> Institut des Sciences de la Terre
> Maison des Géosciences
> 1381 rue de la piscine
> 38400 Saint Martin d-Hères
> FRANCE
> Tel : 04.76.63.51.14
> Tel : 06 14 83 07 05
>
> _______________________________________________
> obspy-users mailing list
> obspy-users at lists.swapbytes.de
> https://lists.swapbytes.de/mailman/listinfo/obspy-users

_______________________________________________
obspy-users mailing list
obspy-users at lists.swapbytes.de
https://lists.swapbytes.de/mailman/listinfo/obspy-users  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.swapbytes.de/pipermail/obspy-users/attachments/20190531/3bf874c1/attachment-0001.html>


More information about the obspy-users mailing list