In [5]:

    

using PyPlot,Seismic


    

In [6]:

    

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   955k      0  0:00:12  0:00:12 --:--:





    




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






    




--  855k

In [7]:

    

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


    

In [ ]: