In [46]:

    

# Managing noise in SeismicJulia
using Seismic, PyPlot


    

In [20]:

    

#Create linear events
d, ext = SeisLinearEvents();
dt = ext.d1;


    

In [41]:

    

#add noise
snr = 1;
dn = SeisAddNoise(d, snr, db=true);


    

In [42]:

    

# Now we can measure the SNR in bd 
SNR = MeasureSNR(d,dn, db=true)


    

    
Out[42]:





0.9999999999999992

In [43]:

    

# or amplitude
SNR = MeasureSNR(d,dn)


    

    
Out[43]:





1.1220184543019633

In [47]:

    

figure(1)
subplot(2, 1, 1)
SeisPlot(d, dy=dt, fignum=1)
subplot(2 ,1, 2)
SeisPlot(dn, dy=dt, fignum=1)


    

    













    
Out[47]:





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

In [ ]: