[obspy-users] obspy-users Digest, Vol 42, Issue 10

gonzalo antonio fernandez marañon tonino.gafm at gmail.com
Thu May 15 16:28:45 CEST 2014


Dear Lion, Bernhard, Tobias

Thank you very much for your time and help, I will apply what you told me
many thanks have a wonderful day!!!

Best Regards.

Gonzalo


2014-05-15 6:00 GMT-04:00 <obspy-users-request at lists.swapbytes.de>:

> Send obspy-users mailing list submissions to
>         obspy-users at lists.swapbytes.de
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         http://lists.swapbytes.de/mailman/listinfo/obspy-users
> or, via email, send a message with subject or body 'help' to
>         obspy-users-request at lists.swapbytes.de
>
> You can reach the person managing the list at
>         obspy-users-owner at lists.swapbytes.de
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of obspy-users digest..."
>
>
> Today's Topics:
>
>    1. Probabilistic Power Spectral Densities Questions
>       (gonzalo antonio fernandez mara?on)
>    2. Re: Probabilistic Power Spectral Densities Questions
>       (Lion Krischer)
>    3. Re: Probabilistic Power Spectral Densities Questions
>       (Bernhard Morgenstern)
>    4. Re: Probabilistic Power Spectral Densities Questions
>       (Tobias Megies)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 14 May 2014 22:49:32 -0400
> From: gonzalo antonio fernandez mara?on  <tonino.gafm at gmail.com>
> To: users at obspy.org
> Subject: [obspy-users] Probabilistic Power Spectral Densities
>         Questions
> Message-ID:
>         <
> CAE897hgzMWgxHhUQf0-UHqE0pyCTFX_EYTpyj2r-6jQsACsLMg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear Obspy forum.
>
> I am glad to learn Obspy, sorry for the following questions and thank you
> for your time.
>
> I am trying to get the PPSD following the web page example but doing some
> modification, let me show you:
>
>
> st=read('F:\\+0n1n0\\JICA\\strong
> motion\\Acel02_Umsa\\AC02.HHZ.12-0000.mseed')
>
> parser=Parser('F:\\+0n1n0\\JICA\\strong
> motion\\Acel02_Umsa\\AC02_z_ok_dataless.seed')
>
> paz=parser.getPAZ(st)
>
>
>
> Then I got the following error:
>
> SEEDParserException: No channel found with the given SEED id: 1
> Trace(s) in Stream:
> 12.AC02.04.HHZ | 2013-08-24T00:00:03.300000Z -
> 2013-08-24T01:00:02.025000Z | 40.0 Hz, 143950 samples
>
> my dataless file has the same informmation that the waveform, let me
> show you the waveform stats
>
>
> Stats({'network': '12', '_format': 'MSEED', 'mseed':
> AttribDict({'record_length': 512, 'encoding': 'STEIM1', 'filesize':
> 303616L, 'dataquality': 'D', 'number_of_records': 593L, 'byteorder':
> '>'}), 'sampling_rate': 40.0, 'station': 'AC02', 'location': '04',
> 'starttime': UTCDateTime(2013, 8, 24, 0, 0, 3, 300000), 'delta':
> 0.025, 'calib': 1.0, 'npts': 143950, 'endtime': UTCDateTime(2013, 8,
> 24, 1, 0, 2, 25000), 'channel': 'HHZ'})
>
>
> Can you help me or give me some tips to keep going ahead with obspy.
>
> Thank you for your time, best regards.
>
> Gonzalo
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> http://lists.swapbytes.de/pipermail/obspy-users/attachments/20140514/738cbeb4/attachment-0001.html
> >
>
> ------------------------------
>
> Message: 2
> Date: Thu, 15 May 2014 08:08:19 +0200
> From: Lion Krischer <krischer at geophysik.uni-muenchen.de>
> To: users at obspy.org
> Subject: Re: [obspy-users] Probabilistic Power Spectral Densities
>         Questions
> Message-ID:
>         <BB89B507-F586-41F8-95CA-2ABFE4629BC2 at geophysik.uni-muenchen.de>
> Content-Type: text/plain; charset=iso-8859-1
>
> Hi Gonzalo,
>
> the parser.getPAZ() method takes a string denoting the SEED id, not a
> stream object.
>
>
> http://docs.obspy.org/packages/autogen/obspy.xseed.parser.Parser.getPAZ.html
>
> if you do
>
> >>> paz = parser.getPAZ(st[0].id, st[0].stats.starttime)
>
> it should work. This will extract the poles and zeros for the first trace
> in your stream objects. You should always also specify the time as
> instrument responses might change with time.
>
> Cheers!
>
> Lion
>
>
> On 15 May 2014, at 04:49, gonzalo antonio fernandez mara?on <
> tonino.gafm at gmail.com> wrote:
>
> > Dear Obspy forum.
> >
> > I am glad to learn Obspy, sorry for the following questions and thank
> you for your time.
> >
> > I am trying to get the PPSD following the web page example but doing
> some modification, let me show you:
> >
> >
> > st=read('F:\\+0n1n0\\JICA\\strong
> motion\\Acel02_Umsa\\AC02.HHZ.12-0000.mseed')
> >
> > parser=Parser('F:\\+0n1n0\\JICA\\strong
> motion\\Acel02_Umsa\\AC02_z_ok_dataless.seed')
> >
> > paz=parser.getPAZ(st)
> >
> >
> > Then I got the following error:
> >
> > SEEDParserException: No channel found with the given SEED id: 1 Trace(s)
> in Stream:
> > 12.AC02.04.HHZ | 2013-08-24T00:00:03.300000Z -
> 2013-08-24T01:00:02.025000Z | 40.0 Hz, 143950 samples
> >
> > my dataless file has the same informmation that the waveform, let me
> show you the waveform stats
> >
> >
> > Stats({'network': '12', '_format': 'MSEED', 'mseed':
> AttribDict({'record_length': 512, 'encoding': 'STEIM1', 'filesize':
> 303616L, 'dataquality': 'D', 'number_of_records': 593L, 'byteorder': '>'}),
> 'sampling_rate': 40.0, 'station': 'AC02', 'location': '04', 'starttime':
> UTCDateTime(2013, 8, 24, 0, 0, 3, 300000), 'delta': 0.025, 'calib': 1.0,
> 'npts': 143950, 'endtime': UTCDateTime(2013, 8, 24, 1, 0, 2, 25000),
> 'channel': 'HHZ'})
> >
> >
> > Can you help me or give me some tips to keep going ahead with obspy.
> > Thank you for your time, best regards.
> >
> > Gonzalo
> >
> >
> > _______________________________________________
> > obspy-users mailing list
> > obspy-users at lists.swapbytes.de
> > http://lists.swapbytes.de/mailman/listinfo/obspy-users
>
>
>
> ------------------------------
>
> Message: 3
> Date: Thu, 15 May 2014 10:59:53 +0200
> From: Bernhard Morgenstern <bmorgenstern at geophysik.uni-muenchen.de>
> To: users at obspy.org
> Subject: Re: [obspy-users] Probabilistic Power Spectral Densities
>         Questions
> Message-ID: <53748209.1020701 at geophysik.uni-muenchen.de>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Hi Gonzalo!
>
> You can simply provide the parser object to the PPSD at instantiation:
>
> ppsd = PPSD(st[0].stats, parser=parser)
>
> You can then add() new data to the PPSD and the correct response
> information will be selected from the parser automatically.
>
> Hope that helps,
> Bernhard
>
>
> On 15/05/14 04:49, gonzalo antonio fernandez mara?on wrote:
> > Dear Obspy forum.
> >
> > I am glad to learn Obspy, sorry for the following questions and thank you
> > for your time.
> >
> > I am trying to get the PPSD following the web page example but doing some
> > modification, let me show you:
> >
> >
> > st=read('F:\\+0n1n0\\JICA\\strong
> > motion\\Acel02_Umsa\\AC02.HHZ.12-0000.mseed')
> >
> > parser=Parser('F:\\+0n1n0\\JICA\\strong
> > motion\\Acel02_Umsa\\AC02_z_ok_dataless.seed')
> >
> > paz=parser.getPAZ(st)
> >
> >
> >
> > Then I got the following error:
> >
> > SEEDParserException: No channel found with the given SEED id: 1
> > Trace(s) in Stream:
> > 12.AC02.04.HHZ | 2013-08-24T00:00:03.300000Z -
> > 2013-08-24T01:00:02.025000Z | 40.0 Hz, 143950 samples
> >
> > my dataless file has the same informmation that the waveform, let me
> > show you the waveform stats
> >
> >
> > Stats({'network': '12', '_format': 'MSEED', 'mseed':
> > AttribDict({'record_length': 512, 'encoding': 'STEIM1', 'filesize':
> > 303616L, 'dataquality': 'D', 'number_of_records': 593L, 'byteorder':
> > '>'}), 'sampling_rate': 40.0, 'station': 'AC02', 'location': '04',
> > 'starttime': UTCDateTime(2013, 8, 24, 0, 0, 3, 300000), 'delta':
> > 0.025, 'calib': 1.0, 'npts': 143950, 'endtime': UTCDateTime(2013, 8,
> > 24, 1, 0, 2, 25000), 'channel': 'HHZ'})
> >
> >
> > Can you help me or give me some tips to keep going ahead with obspy.
> >
> > Thank you for your time, best regards.
> >
> > Gonzalo
> >
> >
> >
> > _______________________________________________
> > obspy-users mailing list
> > obspy-users at lists.swapbytes.de
> > http://lists.swapbytes.de/mailman/listinfo/obspy-users
> >
>
>
>
> ------------------------------
>
> Message: 4
> Date: Thu, 15 May 2014 11:16:37 +0200
> From: Tobias Megies <megies at geophysik.uni-muenchen.de>
> To: users at obspy.org
> Subject: Re: [obspy-users] Probabilistic Power Spectral Densities
>         Questions
> Message-ID: <537485F5.4070101 at geophysik.uni-muenchen.de>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> To follow up on this, as the PPSD's docs page states, it is much safer
> to use a Parser object instead of static PAZ dictionary. I would
> strongly advise you go that way.
>
> best,
> Tobias
>
>
> On 15/05/14 10:59, Bernhard Morgenstern wrote:
> > Hi Gonzalo!
> >
> > You can simply provide the parser object to the PPSD at instantiation:
> >
> > ppsd = PPSD(st[0].stats, parser=parser)
> >
> > You can then add() new data to the PPSD and the correct response
> > information will be selected from the parser automatically.
> >
> > Hope that helps,
> > Bernhard
> >
> >
> > On 15/05/14 04:49, gonzalo antonio fernandez mara?on wrote:
> >> Dear Obspy forum.
> >>
> >> I am glad to learn Obspy, sorry for the following questions and thank
> you
> >> for your time.
> >>
> >> I am trying to get the PPSD following the web page example but doing
> some
> >> modification, let me show you:
> >>
> >>
> >> st=read('F:\\+0n1n0\\JICA\\strong
> >> motion\\Acel02_Umsa\\AC02.HHZ.12-0000.mseed')
> >>
> >> parser=Parser('F:\\+0n1n0\\JICA\\strong
> >> motion\\Acel02_Umsa\\AC02_z_ok_dataless.seed')
> >>
> >> paz=parser.getPAZ(st)
> >>
> >>
> >>
> >> Then I got the following error:
> >>
> >> SEEDParserException: No channel found with the given SEED id: 1
> >> Trace(s) in Stream:
> >> 12.AC02.04.HHZ | 2013-08-24T00:00:03.300000Z -
> >> 2013-08-24T01:00:02.025000Z | 40.0 Hz, 143950 samples
> >>
> >> my dataless file has the same informmation that the waveform, let me
> >> show you the waveform stats
> >>
> >>
> >> Stats({'network': '12', '_format': 'MSEED', 'mseed':
> >> AttribDict({'record_length': 512, 'encoding': 'STEIM1', 'filesize':
> >> 303616L, 'dataquality': 'D', 'number_of_records': 593L, 'byteorder':
> >> '>'}), 'sampling_rate': 40.0, 'station': 'AC02', 'location': '04',
> >> 'starttime': UTCDateTime(2013, 8, 24, 0, 0, 3, 300000), 'delta':
> >> 0.025, 'calib': 1.0, 'npts': 143950, 'endtime': UTCDateTime(2013, 8,
> >> 24, 1, 0, 2, 25000), 'channel': 'HHZ'})
> >>
> >>
> >> Can you help me or give me some tips to keep going ahead with obspy.
> >>
> >> Thank you for your time, best regards.
> >>
> >> Gonzalo
> >>
> >>
> >>
> >> _______________________________________________
> >> obspy-users mailing list
> >> obspy-users at lists.swapbytes.de
> >> http://lists.swapbytes.de/mailman/listinfo/obspy-users
> >>
> >
> > _______________________________________________
> > obspy-users mailing list
> > obspy-users at lists.swapbytes.de
> > http://lists.swapbytes.de/mailman/listinfo/obspy-users
>
> --
> Dipl.-Geophys. Tobias Megies
>
> Geophysikalisches Observatorium
> Ludwigsh?he 8
> 82256 F?rstenfeldbruck
>
> Ludwig-Maximilians-Universit?t
> Department f?r Geo- und Umweltwissenschaften
> Sektion Geophysik
> Theresienstrasse 41/IV
> 80333 M?nchen
>
> Tel: +49 (0) 89 2180-73981
>       +49 (0) 89 2180-4326
> Mail: tobias.megies at geophysik.uni-muenchen.de
>
>
> ------------------------------
>
> _______________________________________________
> obspy-users mailing list
> obspy-users at lists.swapbytes.de
> http://lists.swapbytes.de/mailman/listinfo/obspy-users
>
> End of obspy-users Digest, Vol 42, Issue 10
> *******************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.swapbytes.de/pipermail/obspy-users/attachments/20140515/0c46b069/attachment-0001.html>


More information about the obspy-users mailing list