Thread: PROBLEM WITH SPECTROGRAM IN SAC 32 BITS

Started: 2015-09-10 19:25:39
Last activity: 2015-09-11 02:22:32
Topics: SAC Help
Silvana Spagnotto
2015-09-10 19:25:39
Hello:
I need plot a spectrogram, but this script only work in SAC 64 bits

sac > r 2005.*SAC
SAC> rmean
SAC> rtr
SAC> taper
SAC> interpolate delta 0.2
SAC> spectrogram window 200 slice 50 method mem order 400 ymin 0.0 ymax 0.1
cbar off noscaling

in 32 bits the error is

Window size: 999 Overlap: 749 FFT size: 2048
Spectrogram dimensions are 512 by 3 .
X Error of failed request: BadMatch (invalid parameter attributes)
Major opcode of failed request: 72 (X_PutImage)
Serial number of failed request: 69
Current serial number in output stream: 70



Attachment data

  • matt miller
    2015-09-10 21:59:49

    Silvana,

    The file 2005.102.03.46.00.0000.XJ.HUDS..HHN.D.SAC is only 300 seconds
    long, and has a high-frequency component. I would suggest reducing the
    values of window and slice in the spectrogram command, and increasing
    ymax.

    Also if you cut around the event then it looks better!

    qdp off
    cut 50 80
    r 2005.102.03.46.00.0000.XJ.HUDS..HHN.D.SAC
    rmean
    rtr
    taper
    spectrogram window 2 slice 0.2 method mem order 400 ymin 0 ymax 20 cbar
    off noscaling








    Hello:
    I need plot a spectrogram, but this script only work in SAC 64 bits

    sac > r 2005.*SAC
    SAC> rmean
    SAC> rtr
    SAC> taper
    SAC> interpolate delta 0.2
    SAC> spectrogram window 200 slice 50 method mem order 400 ymin 0.0 ymax
    0.1
    cbar off noscaling

    in 32 bits the error is

    Window size: 999 Overlap: 749 FFT size: 2048
    Spectrogram dimensions are 512 by 3 .
    X Error of failed request: BadMatch (invalid parameter attributes)
    Major opcode of failed request: 72 (X_PutImage)
    Serial number of failed request: 69
    Current serial number in output stream: 70



    Attachment data

    ----------------------
    SAC Help (http://ds.iris.edu/message-center/topic/sac-help/)

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



  • Milton Plasencia
    2015-09-11 02:22:32
    Silvana,

    Perhaps your trace have too much points… i do

    SAC> cut 55 90
    SAC> r ./2005.102.03.46.00.0000.XJ.HUDS..HHN.D.SAC
    SAC> synch
    SAC> p
    SAC> spectrogram
    Window size: 200 Overlap: 100 FFT size: 512
    Warning: ORDER begin reset to 1/2 length of correlation window.
    ORDER = 100
    Spectrogram dimensions are 512 by 34 .

    It works, i cut the trace between 55 to 90 s

    HTH,

    Milton





    On Sep 10, 2015, at 19:08, Silvana Spagnotto <pampa113<at>gmail.com> wrote:

    spectrogram window 200 slice 50 method mem order 400 ymin 0.0 ymax 0.1 cbar off noscaling


17:06:32 v.01697673