In [3]:

    

using Seismic


    

In [ ]:

    

using PyPlot,PyCall

include("plot_demo.jl")
include("plot_demo2.jl")


    

In [2]:

    

download("http://seismic.physics.ualberta.ca/data/prestack_section.su","section.su");


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100 65.8M  100 65.8M    0     0  27.0M      0  0:00:02  0:00:02 --:--:-- 27.0M

In [ ]:

    

SegyToSeis("section.su","section",format="su", input_type="ieee");


    

In [3]:

    

d,h,ext=SeisRead("section");


    

In [6]:

    

SeisPlotCoordinates(h,title="Acquisition geometry", xlabel="X(m)",ylabel="Y(m)")


    

    













    
Out[6]:





PyObject <matplotlib.collections.PathCollection object at 0x7ff3c8b07e90>

In [9]:

    

SeisPlot(d[:,1500:1800],xlabel="Trace",ylabel="Time(s)", dy=0.002,oy=0.3,title="Prestack gathers",ox=1500,wbox=10)


    

    













    
Out[9]:





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

In [1]:

    

param1 = Dict( :dmx=>15, :dmy=>15, :dh=>30, :daz=>45 );
 
param2 = Dict(:style=>"mxmyhaz", :min_imx=>10,:max_imx=>100, :min_imy=>35, :max_imy=>45,
               :min_ih=>1, :max_ih=>6, :min_iaz=>0, :max_iaz=>7);


    

In [4]:

    

SeisGeometry("section"; param1...)


    

In [5]:

    

SeisBinHeaders("section","section_bin"; param1..., param2...);


    

    




nt= 251
The final binned cube will have an approximate size of 0.048 Gb

In [6]:

    

SeisBinData("section","section_bin"; param1..., param2...);


    

In [8]:

    

db,hb,eb=SeisRead("section_bin");


    

In [10]:

    

plot_offsets(hb)


    

    













    
Out[10]:





PyObject <matplotlib.collections.PathCollection object at 0x7f6893777650>

In [ ]: