Thread: Question on filtering signal

Started: 2021-04-28 22:18:45
Last activity: 2021-05-12 20:21:13
Topics: SAC Help
Alessio Lodato
2021-04-28 22:18:45
Dear,

i need to filter a signal derived from a seismometer (output in
displacement) from the transfer function of the instrument.

What is the procedure to be followed in SAC?

Thank you in advance.

Regards,



AL


  • Milton Plasencia
    2021-05-11 21:41:59
    Ciao Alessio,

    Sac have a Spectral Analysis Module, below a list of commands,
    you can read the help for the command that you want to apply.
    For example:
    sac> help bandpass

    HTH,

    Cheers,

    m.


    fft: Performs a discrete Fourier transform.
    keepam: Keep only the amplitude component of the spectral file.
    ifft: Performs an inverse discrete Fourier transform.
    plotsp: Plots spectral data in several different formats.
    writesp: Writes spectral files to disk as "normal" data files.
    readsp: Reads spectral files written by @writesp@.
    lowpass: Applies an IIR lowpass filter.
    highpass: Applies an IIR highpass filter.
    bandpass: Applies an IIR bandpass filter.
    bandrej: Applies an IIR bandreject filter.
    wiener: Applies an adaptive Wiener filter.
    fir: Applies a finite-impulse-response filter.
    hanning: Applies a "hanning" window to each data file.
    unwrap: Computes amplitude and unwrapped phase.
    khronhite: Applies a Khronhite filter to the data.
    benioff: Applies a Benioff filter to the data.
    divomega: Performs integration in the frequency domain.
    mulomega: Performs differentiation in the frequency domain.
    hilbert: Applies a Hilbert transform.
    convolve: Computes the convolution of waveforms.
    correlate: Computes the auto- and cross- correlation functions.
    envelope: Computes the envelope function using a Hilbert transform.


    ************************************
    Milton P. PLASENCIA LINARES
    Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS
    Borgo Grotta Gigante 42/C
    (34010) Sgonico - Trieste - Italia
    Tel: +39 040 2140 156 (Udine)
    Tel: +39 040 2140 256 (Trieste)
    Cell.: +39 331 6481 935

    E-mail: mplasencia<at>inogs.it

    Antarctic Seismographic Argentinean Italian Network (ASAIN)
    **************************************


    On Tue, 11 May 2021 at 21:21, Alessio Lodato (via IRIS) <
    sac-help-bounce<at>lists.ds.iris.edu> wrote:

    Dear,

    i need to filter a signal derived from a seismometer (output in
    displacement) from the transfer function of the instrument.

    What is the procedure to be followed in SAC?

    Thank you in advance.

    Regards,



    AL

    ----------------------
    SAC Help
    Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
    Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu

    Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
    Update subscription preferences at http://ds.iris.edu/account/profile/


    • Peter Schmidt
      2021-05-11 21:54:46
      Hi Alessio


      If I understand your question correctly you want to deconvolve the instrument response, if so then have a look at the transfer command. Sorry have not used it myself but try:

      sac> help transfer

      for further insight


      \p


      On 2021-05-11 21:43, Milton Plasencia (via IRIS) wrote:
      Ciao Alessio,

      Sac have a Spectral Analysis Module, below a list of commands,
      you can read the help for the command that you want to apply.
      For example:
      sac> help bandpass

      HTH,

      Cheers,

      m.


      fft: Performs a discrete Fourier transform.
      keepam: Keep only the amplitude component of the spectral file.
      ifft: Performs an inverse discrete Fourier transform.
      plotsp: Plots spectral data in several different formats.
      writesp: Writes spectral files to disk as "normal" data files.
      readsp: Reads spectral files written by @writesp@.
      lowpass: Applies an IIR lowpass filter.
      highpass: Applies an IIR highpass filter.
      bandpass: Applies an IIR bandpass filter.
      bandrej: Applies an IIR bandreject filter.
      wiener: Applies an adaptive Wiener filter.
      fir: Applies a finite-impulse-response filter.
      hanning: Applies a "hanning" window to each data file.
      unwrap: Computes amplitude and unwrapped phase.
      khronhite: Applies a Khronhite filter to the data.
      benioff: Applies a Benioff filter to the data.
      divomega: Performs integration in the frequency domain.
      mulomega: Performs differentiation in the frequency domain.
      hilbert: Applies a Hilbert transform.
      convolve: Computes the convolution of waveforms.
      correlate: Computes the auto- and cross- correlation functions.
      envelope: Computes the envelope function using a Hilbert transform.


      ************************************
      Milton P. PLASENCIA LINARES
      Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS
      Borgo Grotta Gigante 42/C
      (34010) Sgonico - Trieste - Italia
      Tel: +39 040 2140 156 (Udine)
      Tel: +39 040 2140 256 (Trieste)
      Cell.: +39 331 6481 935

      E-mail: mplasencia<at>inogs.it<mplasencia<at>inogs.it>

      Antarctic Seismographic Argentinean Italian Network (ASAIN)
      **************************************


      On Tue, 11 May 2021 at 21:21, Alessio Lodato (via IRIS) <sac-help-bounce<at>lists.ds.iris.edu<sac-help-bounce<at>lists.ds.iris.edu>> wrote:
      Dear,
      i need to filter a signal derived from a seismometer (output in displacement) from the transfer function of the instrument.
      What is the procedure to be followed in SAC?
      Thank you in advance.
      Regards,

      AL

      ----------------------
      SAC Help
      Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu<sac-help-unsubscribe<at>lists.ds.iris.edu>

      Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
      Update subscription preferences at http://ds.iris.edu/account/profile/




      ----------------------
      SAC Help
      Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu<sac-help-unsubscribe<at>lists.ds.iris.edu>

      Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
      Update subscription preferences at http://ds.iris.edu/account/profile/


      --
      ********************************************************************************
      Peter Schmidt Tel: +46-18-4717104
      Swedish National Seismological Network (SNSN) Mobile: +46-73-3190975
      Dept. of Earth Sciences:geophysics e-mail: peter.schmidt<at>geo.uu.se<peter.schmidt<at>geo.uu.se>
      Uppsala University
      Villavagen 16
      SE-75236 Uppsala
      ********************************************************************************









      När du har kontakt med oss på Uppsala universitet med e-post så innebär det att vi behandlar dina personuppgifter. För att läsa mer om hur vi gör det kan du läsa här: http://www.uu.se/om-uu/dataskydd-personuppgifter/

      E-mailing Uppsala University means that we will process your personal data. For more information on how this is performed, please read here: http://www.uu.se/en/about-uu/data-protection-policy

      • Alessio Lodato
        2021-05-12 13:59:21
        Thank you very much,

        i have another question, i want to use SAC on Windows10, how it is possible to do it?

        Best regards,



        AL



        Da: sac-help-bounce<at>lists.ds.iris.edu [sac-help-bounce<at>lists.ds.iris.edu]
        Inviato: martedì 11 maggio 2021 21:55
        A: SAC Help <sac-help<at>lists.ds.iris.edu>
        Oggetto: Re: [IRIS][sac-help] Question on filtering signal



        Hi Alessio



        If I understand your question correctly you want to deconvolve the instrument response, if so then have a look at the transfer command. Sorry have not used it myself but try:

        sac> help transfer

        for further insight



        \p



        On 2021-05-11 21:43, Milton Plasencia (via IRIS) wrote:

        Ciao Alessio,



        Sac have a Spectral Analysis Module, below a list of commands,

        you can read the help for the command that you want to apply.

        For example:

        sac> help bandpass



        HTH,



        Cheers,



        m.





        fft: Performs a discrete Fourier transform.

        keepam: Keep only the amplitude component of the spectral file.

        ifft: Performs an inverse discrete Fourier transform.

        plotsp: Plots spectral data in several different formats.

        writesp: Writes spectral files to disk as "normal" data files.

        readsp: Reads spectral files written by @writesp@.

        lowpass: Applies an IIR lowpass filter.

        highpass: Applies an IIR highpass filter.

        bandpass: Applies an IIR bandpass filter.

        bandrej: Applies an IIR bandreject filter.

        wiener: Applies an adaptive Wiener filter.

        fir: Applies a finite-impulse-response filter.

        hanning: Applies a "hanning" window to each data file.

        unwrap: Computes amplitude and unwrapped phase.

        khronhite: Applies a Khronhite filter to the data.

        benioff: Applies a Benioff filter to the data.

        divomega: Performs integration in the frequency domain.

        mulomega: Performs differentiation in the frequency domain.

        hilbert: Applies a Hilbert transform.

        convolve: Computes the convolution of waveforms.

        correlate: Computes the auto- and cross- correlation functions.

        envelope: Computes the envelope function using a Hilbert transform.






        ************************************

        Milton P. PLASENCIA LINARES

        Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS

        Borgo Grotta Gigante 42/C

        (34010) Sgonico - Trieste - Italia

        Tel: +39 040 2140 156 (Udine)

        Tel: +39 040 2140 256 (Trieste)

        Cell.: +39 331 6481 935



        E-mail: mplasencia<at>inogs.it <mplasencia<at>inogs.it>



        Antarctic Seismographic Argentinean Italian Network (ASAIN)

        **************************************





        On Tue, 11 May 2021 at 21:21, Alessio Lodato (via IRIS) <sac-help-bounce<at>lists.ds.iris.edu <sac-help-bounce<at>lists.ds.iris.edu> > wrote:

        Dear,

        i need to filter a signal derived from a seismometer (output in displacement) from the transfer function of the instrument.

        What is the procedure to be followed in SAC?

        Thank you in advance.

        Regards,



        AL


        ----------------------
        SAC Help
        Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu <sac-help-unsubscribe<at>lists.ds.iris.edu>

        Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
        Update subscription preferences at http://ds.iris.edu/account/profile/






        ----------------------
        SAC Help
        Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu <sac-help-unsubscribe<at>lists.ds.iris.edu>

        Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
        Update subscription preferences at http://ds.iris.edu/account/profile/

        --
        ********************************************************************************
        Peter Schmidt Tel: +46-18-4717104
        Swedish National Seismological Network (SNSN) Mobile: +46-73-3190975
        Dept. of Earth Sciences:geophysics e-mail: peter.schmidt<at>geo.uu.se <peter.schmidt<at>geo.uu.se>
        Uppsala University
        Villavagen 16
        SE-75236 Uppsala
        ********************************************************************************









        När du har kontakt med oss på Uppsala universitet med e-post så innebär det att vi behandlar dina personuppgifter. För att läsa mer om hur vi gör det kan du läsa här: http://www.uu.se/om-uu/dataskydd-personuppgifter/

        E-mailing Uppsala University means that we will process your personal data. For more information on how this is performed, please read here: http://www.uu.se/en/about-uu/data-protection-policy


        • Milton Plasencia
          2021-05-12 17:09:53
          Ciao,
          Sac running on OSX and Linux.

          Perhaps a virtual PC with linux in your Windows box can do the trick.

          M.
          ************************************
          Milton P. PLASENCIA LINARES
          Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS
          Borgo Grotta Gigante 42/C
          (34010) Sgonico - Trieste - Italia
          Tel: +39 040 2140 156 (Udine)
          Tel: +39 040 2140 256 (Trieste)
          Cell.: +39 331 6481 935

          E-mail: mplasencia<at>inogs.it

          Antarctic Seismographic Argentinean Italian Network (ASAIN)
          **************************************


          On Wed, 12 May 2021 at 14:05, Alessio Lodato (via IRIS) <
          sac-help-bounce<at>lists.ds.iris.edu> wrote:

          Thank you very much,

          i have another question, i want to use SAC on Windows10, how it is
          possible to do it?

          Best regards,



          AL



          *Da:* sac-help-bounce<at>lists.ds.iris.edu [mailto:
          sac-help-bounce<at>lists.ds.iris.edu]
          *Inviato:* martedì 11 maggio 2021 21:55
          *A:* SAC Help <sac-help<at>lists.ds.iris.edu>
          *Oggetto:* Re: [IRIS][sac-help] Question on filtering signal



          Hi Alessio



          If I understand your question correctly you want to deconvolve the
          instrument response, if so then have a look at the transfer command. Sorry
          have not used it myself but try:

          sac> help transfer

          for further insight



          \p



          On 2021-05-11 21:43, Milton Plasencia (via IRIS) wrote:

          Ciao Alessio,



          Sac have a Spectral Analysis Module, below a list of commands,

          you can read the help for the command that you want to apply.

          For example:

          sac> help bandpass



          HTH,



          Cheers,



          m.





          fft: Performs a discrete Fourier transform.

          keepam: Keep only the amplitude component of the spectral file.

          ifft: Performs an inverse discrete Fourier transform.

          plotsp: Plots spectral data in several different formats.

          writesp: Writes spectral files to disk as "normal" data files.

          readsp: Reads spectral files written by @writesp@.

          lowpass: Applies an IIR lowpass filter.

          highpass: Applies an IIR highpass filter.

          bandpass: Applies an IIR bandpass filter.

          bandrej: Applies an IIR bandreject filter.

          wiener: Applies an adaptive Wiener filter.

          fir: Applies a finite-impulse-response filter.

          hanning: Applies a "hanning" window to each data file.

          unwrap: Computes amplitude and unwrapped phase.

          khronhite: Applies a Khronhite filter to the data.

          benioff: Applies a Benioff filter to the data.

          divomega: Performs integration in the frequency domain.

          mulomega: Performs differentiation in the frequency domain.

          hilbert: Applies a Hilbert transform.

          convolve: Computes the convolution of waveforms.

          correlate: Computes the auto- and cross- correlation functions.

          envelope: Computes the envelope function using a Hilbert transform.




          ************************************

          Milton P. PLASENCIA LINARES

          Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS

          Borgo Grotta Gigante 42/C

          (34010) Sgonico - Trieste - Italia

          Tel: +39 040 2140 156 (Udine)

          Tel: +39 040 2140 256 (Trieste)

          Cell.: +39 331 6481 935



          E-mail: mplasencia<at>inogs.it



          Antarctic Seismographic Argentinean Italian Network (ASAIN)

          **************************************





          On Tue, 11 May 2021 at 21:21, Alessio Lodato (via IRIS) <
          sac-help-bounce<at>lists.ds.iris.edu> wrote:

          Dear,

          i need to filter a signal derived from a seismometer (output in
          displacement) from the transfer function of the instrument.

          What is the procedure to be followed in SAC?

          Thank you in advance.

          Regards,



          AL


          ----------------------
          SAC Help
          Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
          Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu

          Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
          Update subscription preferences at http://ds.iris.edu/account/profile/





          ----------------------

          SAC Help

          Topic home: http://ds.iris.edu/message-center/topic/sac-help/ | Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu



          Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)

          Update subscription preferences at http://ds.iris.edu/account/profile/

          --

          ********************************************************************************

          Peter Schmidt Tel: +46-18-4717104

          Swedish National Seismological Network (SNSN) Mobile: +46-73-3190975

          Dept. of Earth Sciences:geophysics e-mail: peter.schmidt<at>geo.uu.se

          Uppsala University

          Villavagen 16

          SE-75236 Uppsala

          ********************************************************************************









          När du har kontakt med oss på Uppsala universitet med e-post så innebär
          det att vi behandlar dina personuppgifter. För att läsa mer om hur vi gör
          det kan du läsa här: http://www.uu.se/om-uu/dataskydd-personuppgifter/

          E-mailing Uppsala University means that we will process your personal
          data. For more information on how this is performed, please read here:
          http://www.uu.se/en/about-uu/data-protection-policy

          ----------------------
          SAC Help
          Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
          Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu

          Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
          Update subscription preferences at http://ds.iris.edu/account/profile/


  • Milton Plasencia
    2021-05-12 19:55:37
    Also give a look on:

    http://www.cygwin.com
    it installs a Linux on your windows box
    m.
    ************************************
    Milton P. PLASENCIA LINARES
    Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS
    Borgo Grotta Gigante 42/C
    (34010) Sgonico - Trieste - Italia
    Tel: +39 040 2140 156 (Udine)
    Tel: +39 040 2140 256 (Trieste)
    Cell.: +39 331 6481 935

    E-mail: mplasencia<at>inogs.it

    Antarctic Seismographic Argentinean Italian Network (ASAIN)
    **************************************


    On Tue, 11 May 2021 at 21:21, Alessio Lodato (via IRIS) <
    sac-help-bounce<at>lists.ds.iris.edu> wrote:

    Dear,

    i need to filter a signal derived from a seismometer (output in
    displacement) from the transfer function of the instrument.

    What is the procedure to be followed in SAC?

    Thank you in advance.

    Regards,



    AL

    ----------------------
    SAC Help
    Topic home: http://ds.iris.edu/message-center/topic/sac-help/ |
    Unsubscribe: sac-help-unsubscribe<at>lists.ds.iris.edu

    Sent from the IRIS Message Center (http://ds.iris.edu/message-center/)
    Update subscription preferences at http://ds.iris.edu/account/profile/


20:37:04 v.01697673