In [1]:

    

using PyPlot,Seismic


    

    




INFO: Loading help data...

In [2]:

    

download("http://certmapper.cr.usgs.gov/nersl/NPRA/seismic/1979/616_79/PROCESSED/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
100 11.3M  100 11.3M    0     0  1460k      0  0:00:07  0:00:07 --:--:-- 1625k






    




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

In [3]:

    

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