In [6]:

    

using PyPlot,Seismic


    

    




INFO: Loading help data...

In [2]:

    

download("http://certmapper.cr.usgs.gov/nersl/NPRA/seismic/1979/616_79/PROCESSED/616_79_PR.SGY","616_79_PR.SGY");


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100 11.3M  100 11.3M    0     0   762k      0  0:00:15  0:00:15 --:--:--  798k

In [29]:

    

SegyToSeis("616_79_PR.SGY","616_79_PR_ieee.seis",input_type="ieee");
SegyToSeis("616_79_PR.SGY","616_79_PR_ibm.seis",input_type="ibm");


    

    




number of traces: 1908
number of samples per trace: 1500
number of traces: 1908
number of samples per trace: 1500

In [16]:

    

d1,h,e = SeisRead("616_79_PR_ieee.seis");
d2,h,e = SeisRead("616_79_PR_ibm.seis");

println(maximum(d1[:]))
println(maximum(d2[:]))


    

    




213294.0
1.3471252e9

In [19]:

    

SeisPlot(d1,e,pclip=99.99,cmap="PuOr",wbox=9)


    

    













    
Out[19]:





PyObject <matplotlib.image.AxesImage object at 0x127f51510>

In [20]:

    

SeisPlot(d2,e,pclip=99.99,cmap="PuOr",wbox=9)


    

    













    
Out[20]:





PyObject <matplotlib.image.AxesImage object at 0x127d831d0>

In [3]:

    

d,h,e = SeisRead("616_79_PR.seis");
SeisPlot(d[1:500,:],e,cmap="PuOr",wbox=9);


    

In [4]:

    

SeisHeaderInfo("616_79_PR.seis")


    

    




Displaying information for 616_79_PR.seis (1908 traces):
       Key          Minimum          Maximum             Mean
=============================================================
  tracenum            1.000         1908.000          954.500
        o1            0.000            0.000            0.000
        n1         1500.000         1500.000         1500.000
        d1            0.004            0.004            0.004
        sx            0.000            0.000            0.000
        sy            0.000            0.000            0.000
        gx            0.000            0.000            0.000
        gy            0.000            0.000            0.000
        mx            0.000            0.000            0.000
        my            0.000            0.000            0.000
        hx            0.000            0.000            0.000
        hy            0.000            0.000            0.000
         h            0.000            0.000            0.000
        az            0.000            0.000            0.000
       ang            0.000            0.000            0.000
       isx            0.000            0.000            0.000
       isy            0.000            0.000            0.000
       igx            0.000            0.000            0.000
       igy            0.000            0.000            0.000
       imx          101.000         2008.000         1054.500
       imy            0.000            0.000            0.000
       ihx            0.000            0.000            0.000
       ihy            0.000            0.000            0.000
        ih            0.000            0.000            0.000
       iaz            0.000            0.000            0.000
      iang            0.000            0.000            0.000
     selev            0.000            0.000            0.000
     gelev            0.000            0.000            0.000
     sstat            0.000            0.000            0.000
     gstat            0.000            0.000            0.000
      trid            1.000            1.000            1.000

In [1]:

    

using Seismic
SeisWindow("616_79_PR.seis","616_79_PR_imx_501_1000.seis",key=["imx"],minval=[501],maxval=[1000])


    

    




type: SeisProcessHeaders: in typeassert, expected Symbol, got ASCIIString
while loading In[1], in expression starting on line 2

 in SeisProcessHeaders at /Users/astanton/Seismic/src/Utils/SeisProcessHeaders.jl:27
 in SeisWindowHeaders at /Users/astanton/Seismic/src/Utils/SeisWindowHeaders.jl:12
 in SeisWindow at /Users/astanton/Seismic/src/Utils/SeisWindow.jl:5

In [27]:

    

plot(d1[100:500,1000]/maximum(abs(d1[100:500,1000]))); hold; plot(d2[100:500,1000]/maximum(abs(d2[100:500,1000])))


    

    













    
Out[27]:





1-element Array{Any,1}:
 PyObject <matplotlib.lines.Line2D object at 0x130563490>

In [ ]: