In [1]:

    

using PyPlot,Seismic


    

    




INFO: Loading help data...

In [2]:

    

download("http://s3.amazonaws.com/teapot/npr3_gathers.sgy","npr3_gathers.sgy")


    

    




  % Total    % Received % Xferd  Average Speed   Time    Time     Time  Current
                                 Dload  Upload   Total   Spent    Left  Speed
100 5824M  100 5824M    0     0  2283k      0  0:43:31  0:43:31 --:--:-- 2236k






    
Out[2]:





"npr3_gathers.sgy"

In [3]:

    

SegyToSeis("npr3_gathers.sgy","npr3_gathers.seis");


    

    




number of traces: 723991
number of samples per trace: 2049

In [4]:

    

SeisGeometry("npr3_gathers.seis",ang=90,omx=788937,omy=938845,dmx=110,dmy=110,oh=0,oaz=0,dh=420,daz=45)


    

In [3]:

    

h = SeisReadHeaders("npr3_gathers.seis")
im1 = SeisPlotCoordinates(h,style="sxsygxgy")
im2 = SeisPlotCoordinates(h,style="fold",cmap="jet",vmin=0,vmax=50,aspect="auto",xlabel="imx",ylabel="imy")


    

    













    













    













    
Out[3]:





PyObject <matplotlib.image.AxesImage object at 0x115c6b850>

In [ ]: