In [1]:

    

using PyPlot,Seismic


    

    




INFO: Precompiling module PyPlot.

In [20]:

    

#Create linear events
d,ext = SeisLinearEvents(p1 = [-.001, 0.0015],tau=[1., 1./3.],dx1=5); 

#Randomly decimate
dec = SeisDecimate(d);


    

In [22]:

    

param = Dict(:Niter=>100,:fmax=>60,:padt=>2,:padx=>2,:dt=>0.004)
dpocs = SeisPOCS(dec;param...);


    

In [21]:

    

figure(1, figsize=(10, 5))
subplot(121)
SeisPlot(dec,cmap="seismic",fignum=1,pclip=200,fignum=1)
subplot(122)
SeisPlot(dec,plot_type="FK", cmap="seismic", dy=0.004,fignum=1,hbox=5,pclip=200,fignum=1)


    

    













    
Out[21]:





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

In [23]:

    

figure(2,figsize=(10, 5))
subplot(121)
SeisPlot(dpocs[:,:,1,1,1],cmap="seismic",fignum=2,pclip=200)
subplot(122)
SeisPlot(dpocs[:,:,1,1,1],plot_type="FK",cmap="seismic",dy=0.004,fignum=2,pclip=200)


    

    













    
Out[23]:





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