In [1]:

    

using PyPlot,Seismic


    

In [2]:

    

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,ext = SeisRead("data_with_noise");
plotpar = Dict(:style=>"wiggles",
           :vmin=>-5,:vmax=>5,
           :xlabel=>"X",: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
100  109k  100  109k    0     0   222k      0 --:--:-- --:--:-- --:--:--  223k






    













    




number of traces: 50
number of samples per trace: 501

In [3]:

    

download("http://seismic.physics.ualberta.ca/data/syn_cmp.su","syn_cmp.su");
SegyToSeis("syn_cmp.su","syn_cmp",format="su");
d,h,ext = SeisRead("syn_cmp");
plotpar = Dict(:style=>"wiggles",
           :vmin=>-2,:vmax=>2,
           :xlabel=>"X",: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
 15 88464   15 14073    0     0  66077      0  0:00:01 --:--:--  0:00:01 66070





    













    




number of traces: 38
number of samples per trace: 522






    




100 88464  100 88464    0     0   283k      0 --:--:-- --:--:-- --:--:--  283k

In [17]:

    

download("http://seismic.physics.ualberta.ca/data/small_stack.su","small_stack.su");
SegyToSeis("small_stack.su","small_stack",format="su");
d,h,ext = 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,e_wav = SeisRead("wavelet_for_small_stack");
nt_wav = size(w,1);
plotpar = Dict(:xlabel=>"X",:dx=>1,
           :ylabel=>"Time",:yunits=>"(seconds)",:oy=>0,:dy=>h[1].d1,
           :cmap=>"PuOr",:fignum=>1,
           :title=>"small_stack.su");

subplot(122);
wlong = hcat(w[:]', zeros(nt-nt_wav)');
ax = gca();
ax[:invert_yaxis]();  

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

subplot(121);
SeisPlot(d;pclip=200,plotpar...);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100 20800  100 20800    0     0  43939      0 --:--:-- --:--:-- --:--:-- 43974






    




number of traces: 20
number of samples per trace: 200






    




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





    













    




number of traces: 1
number of samples per trace: 45






    




100   420  100   420    0     0   2179      0 --:--:-- --:--:-- --:--:--  2187

In [16]:

    

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,ext = SeisRead("gom_cdp_nmo");
nx = length(h);
dx = abs.(h[end].h - h[1].h)/nx;
plotpar = Dict(           :aspect=>5000,
           :xlabel=>"Offset",:xunits=>"(ft)",: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
100  650k  100  650k    0     0  1238k      0 --:--:-- --:--:-- --:--:-- 1237k






    













    




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

In [28]:

    

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",swap_bytes=false,input_type="ieee");
d,h,ext = SeisRead("min_phase_wavelet");
nt = size(d,1);
dt = h[1].d1;
plot(d,collect(0:nt-1)*dt,color="k");
title("min_phase_wavelet.su");
ax = gca();
ax[:invert_yaxis]();
ylabel("Time (s)");
axis([-20,20,(nt-1)*dt,0.])


    

    




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





    













    




number of traces: 1
number of samples per trace: 35






    




100   380  100   380    0     0   1730      0 --:--:-- --:--:-- --:--:--  1735






    
Out[28]:





4-element Array{Float64,1}:
 -20.0  
  20.0  
   0.136
   0.0  

In [ ]: