In [1]:

    

using PyPlot,Seismic


    

In [2]:

    

# Download from SEG repository at AWS..

# https://wiki.seg.org/wiki/Teapot_dome_3D_survey

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  6384k      0  0:15:34  0:15:34 --:--:-- 4579k0:06  0:43:00 2513k10k  0:03:49 8708k






    
Out[2]:





"npr3_gathers.sgy"

In [4]:

    

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


    

    




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

In [5]:

    

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


    

In [6]:

    

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")


    

    













    













    




/Users/msacchi/anaconda/lib/python2.7/site-packages/matplotlib/collections.py:590: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison
  if self._edgecolors == str('face'):






    
Out[6]:





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

In [ ]: