In [8]:

    

using PyPlot,Seismic


    

In [9]:

    

download("http://seismic.physics.ualberta.ca/data/data_with_noise.su","data_with_noise.su");
SegyToSeis("data_with_noise.su","data_with_noise",["format"=>"su"]);
d,h = SeisRead("data_with_noise");
plotpar = ["style"=>"wiggles",
           "vmin"=>-5,"vmax"=>5,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"gray",
           "title"=>"data_with_noise.su"];
SeisPlot(d[:,:],plotpar);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0





    




nx=50





    













    




nt=501






    




100  109k  100  109k    0     0  17.0M      0 --:--:-- --:--:-- --:--:-- 26.7M

In [10]:

    

download("http://seismic.physics.ualberta.ca/data/syn_cmp.su","syn_cmp.su");
SegyToSeis("syn_cmp.su","syn_cmp",["format"=>"su"]);
d,h = SeisRead("syn_cmp");
plotpar = ["style"=>"wiggles",
           "vmin"=>-2,"vmax"=>2,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"gray",
           "title"=>"syn_cmp.su"];
SeisPlot(d[:,:],plotpar);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0





    




nx=38





    













    




nt=522






    




100 88464  100 88464    0     0  14.7M      0 --:--:-- --:--:-- --:--:-- 28.1M

In [11]:

    

download("http://seismic.physics.ualberta.ca/data/syn_cmp_mult.su","syn_cmp_mult.su");
SegyToSeis("syn_cmp_mult.su","syn_cmp_mult",["format"=>"su"]);
d,h = SeisRead("syn_cmp_mult");
plotpar = ["style"=>"overlay",
           "vmin"=>-1,"vmax"=>1,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"seismic",
           "title"=>"syn_cmp_mult.su"];
SeisPlot(d[:,:],plotpar);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0





    




nx=49





    













    




nt=1001






    




100  203k  100  203k    0     0  20.1M      0 --:--:-- --:--:-- --:--:-- 28.3M

In [12]:

    

download("http://seismic.physics.ualberta.ca/data/small_stack.su","small_stack.su");
SegyToSeis("small_stack.su","small_stack",["format"=>"su"]);
d,h = SeisRead("small_stack");
nt = size(d,1);
dt = h[1].d1;
download("http://seismic.physics.ualberta.ca/data/wavelet_for_small_stack.su","wavelet_for_small_stack.su");
SegyToSeis("wavelet_for_small_stack.su","wavelet_for_small_stack",["format"=>"su"]);
w,h_wav = SeisRead("wavelet_for_small_stack");
nt_wav = size(w,1);
plotpar = ["style"=>"color",
           "vmin"=>-1,"vmax"=>1,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>dt,
    "cmap"=>"PuOr","fignum"=>1,
           "title"=>"small_stack.su"];

subplot(122);
w = [w[:]' zeros(nt-nt_wav)'];
plot(w[:],[0:nt-1]*dt);
title("wavelet_for_small_stack.su");
suptitle("Stack and Wavelet");
ax = gca();
ax[:invert_yaxis](); 

subplot(121);
SeisPlot(d[:,:],plotpar);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100 20800  100 20800    0     0  4631k      0 --:--:-- --:--:-- --:--:-- 19.8M






    




nx=20





    













    




nt=200
nx=1
nt=45






    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
105   420  105   420    0     0   140k      0 --:--:-- --:--:-- --:--:--  410k

In [13]:

    

download("http://seismic.physics.ualberta.ca/data/gom_cdp_nmo.su","gom_cdp_nmo.su");
SegyToSeis("gom_cdp_nmo.su","gom_cdp_nmo",["format"=>"su"]);
d,h = SeisRead("gom_cdp_nmo");
nx = length(h);
dx = abs(h[end].h - h[1].h)/nx;
plotpar = ["style"=>"color",
           "vmin"=>-5,"vmax"=>5,
           "aspect"=>5000,
           "xlabel"=>"Offset","xunits"=>"(ft)","ox"=>0,"dx"=>dx,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"gray",
           "title"=>"gom_cdp_nmo.su"];
SeisPlot(d[:,:],plotpar);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
  0     0    0     0    0     0      0      0 --:--:-- --:--:-- --:--:--     0





    




nx=92





    













    




nt=1751






    




100  650k  100  650k    0     0  32.4M      0 --:--:-- --:--:-- --:--:-- 37.3M

In [14]:

    

download("http://seismic.physics.ualberta.ca/data/min_phase_wavelet.su","min_phase_wavelet.su");
SegyToSeis("min_phase_wavelet.su","min_phase_wavelet",["format"=>"su"]);
d,h = SeisRead("min_phase_wavelet");
nt = size(d,1);
dt = h[1].d1;

plot(d[:],[0:nt-1]*dt,color="k");
title("min_phase_wavelet.su");
ax = gca();
ax[:invert_yaxis]();
axis([-20,20,(nt-1)*dt,0.])
ylabel("Time (s)");


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
126   380  126   380    0     0   120k      0 --:--:-- --:--:-- --:--:--  371k






    




nx=1
nt=35






    

In [ ]: