In [1]:

    

using Seismic, PyPlot


    

In [4]:

    

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


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100  650k  100  650k    0     0  1143k      0 --:--:-- --:--:-- --:--:-- 1143k

In [6]:

    

SegyToSeis("gom_cdp_nmo.su","gom_cdp_nmo",format="su",input_type="ieee",swap_bytes=true)


    

    




number of traces: 92
number of samples per trace: 1751

In [14]:

    

d,h,ext = SeisRead("gom_cdp_nmo");
dt = ext.d1;
nx = length(h);
dx = abs.(h[end].h - h[1].h)/nx;
println("dt = ",dt," dx = ",dx)


    

    




dt = 0.004 dx = 173.09782

In [17]:

    

SeisPlot(d,aspect=5000,dy=dt,dx=dx,title="GOM_NMO_CDPgather", xlabel="Offset", ylabel="Time", xunits="(ft)", yunits="(seconds)")


    

    













    
Out[17]:





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