In [2]:

    

using Seismic, PyPlot


    

In [2]:

    

param2 = Dict(:style=>"mxmyhaz", :min_imx=>10,:max_imx=>100, :min_imy=>35, :max_imy=>45,
               :min_ih=>1, :max_ih=>6, :min_iaz=>0, :max_iaz=>7);

param3 = Dict( :ix1_WL=>30,:ix1_WO=>5);


    

In [ ]:

    

patch_out,npatch=SeisPatch(file_bin,"patch"; param2...,param3... );


    

In [ ]:

    

SeisUnPatch(patch_out,file_final; param2...,param3...,nt=251);


    

In [3]:

    

#Plot for QC
db,hb,eb=SeisRead("section_bin");


    

In [5]:

    

SeisPlot(db[:,:,5,4,5],xlabel="Trace",ylabel="Time(s)", dy=0.002,oy=0.3,title="Binned data")


    

    













    
Out[5]:





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

In [2]:

    

df,hf,ef=SeisRead("section_final");


    

In [4]:

    

SeisPlot(df[:,:,5,4,5],xlabel="Trace",ylabel="Time(s)",title="Unpatched result",dy=0.002,oy=0.3)


    

    













    
Out[4]:





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

In [3]:

    

dp1,h1,e1=SeisRead("patch_1");
dp2,h2,e2=SeisRead("patch_2");
dp3,h3,e3=SeisRead("patch_3");
dp4,h4,e4=SeisRead("patch_4");


    

In [4]:

    

SeisPlot(hcat(dp1[:,:,5,4,5],dp2[:,6:end,5,4,5],dp3[:,6:end,5,4,5],dp4[:,6:end,5,4,5]),xlabel="Trace",ylabel="Time(s)", dy=0.002,oy=0.3,title="Overlapping patches")


    

    













    
Out[4]:





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

In [ ]: