In [1]:

    

using Seismic, SeismicPlotting


    

    




INFO: Loading help data...

In [2]:

    

download("http://seismic.physics.ualberta.ca/data/sigmoid.rsf","sigmoid.rsf");
download("http://seismic.physics.ualberta.ca/data/sigmoid.rsf@","sigmoid.rsf@");
SegyToSeis("sigmoid.rsf","sigmoid",["format"=>"rsf"]);


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
151   151  151   151    0     0  51659      0 --:--:-- --:--:-- --:--:--  147k
  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100   97k  100   97k    0     0  16.0M      0 --:--:-- --:--:-- --:--:-- 31.7M






    




nx=100
nt=250

In [3]:

    

d,h = SeisRead("sigmoid");
plotpar = ["style"=>"wiggles",
           "vmin"=>-1e-3,"vmax"=>1e-3,
           "aspect"=>"auto",
           "xlabel"=>"X","xunits"=>"","ox"=>0,"dx"=>1,
           "ylabel"=>"Time","yunits"=>"(seconds)","oy"=>0,"dy"=>h[1].d1,
           "cmap"=>"gray",
           "title"=>"sigmoid.rsf"];
SeisPlot(d[:,:],plotpar);


    

    






Figure(PyObject <matplotlib.figure.Figure object at 0x120e99410>)






    

In [ ]: