In [1]:

    

using PyPlot,Seismic


    

In [7]:

    

# I lost link from USGS, data are now in our local rep. 

download("http://seismic.physics.ualberta.ca/data/616_79_PR.SGY","616_79_PR.SGY");
SegyToSeis("616_79_PR.SGY","616_79_PR.seis");


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
 95 11.3M   95 10.8M    0     0  11.1M      0  0:00:01 --:--:--  0:00:01 11.1M





    




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






    




100 11.3M  100 11.3M    0     0  11.1M      0  0:00:01  0:00:01 --:--:-- 11.1M

In [6]:

    

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


    

In [ ]: