1. xarray traces accessor

class sisicepy.traces.traces.traces(xarray_obj)

This is a classe to work on seismic traces using xarray.

__init__(xarray_obj)

Constructor for traces.

Parameters:

xarray_obj (xr.DataArray)

__weakref__

list of weak references to the object

1.1. Loading data

sisicepy.traces.io.datacube_loader.load_datacube(adr, comp=[0], force_freq=None)

Loader for DIGOS DataCube

Parameters:

  • adr (str) – Path to the DataCube file

  • comp (list of int) – List of componant to load

Returns:

xarray Dataset containing the loaded data and metadata

Return type:

xarray.Dataset

sisicepy.traces.io.datacube_loader.load_mfdatacube(adr_list, **kwargs)

Loader for DIGOS DataCube

Parameters:

adr_list (list)

sisicepy.traces.io.datacube_loader.open_day_datacube(path, day, **kwargs)

Loader for DIGOS DataCube for a given day

Parameters:

day (str) – format %Y-%m-%d

1.2. Recorders

sisicepy.traces.recorder.datacube.bob_S100(self)

Apply Digos BreakOutBox for S100 the divide by 4 the response

sisicepy.traces.recorder.datacube.count2volt_datacube(self, Vpp=4.096, nb_bin=16777216)

Convert numeric value to Voltage for DataCube

Parameters:

  • Vpp (float) – peak to peak value, default: 4.096 (unit V)

  • nb_bin (int) – number of integer for the numerisation, default: 2**24

1.3. Captors

sisicepy.traces.captor.digos.geophone_ZPK()

Value for Poles and Zeros function of the Digos Geophone

Zeros: [0,0] Poles :[-19.78+20.20i,-19.78-20.20i] K=27.7 V/(m.s^-1)

sisicepy.traces.captor.geobit.C100_ZPK()

Value for Poles and Zeros function of the Geobit C100

Zeros: [0+0i,0+0i] Poles :[-15.8336-23.4251i,-15.8336+23.4251i] K = 28.8 V/(m.s^-1)

sisicepy.traces.captor.geobit.S100_ZPK()

Value for Poles and Zeros function of the Geobit S100

Zeros: [0+0i,0+0i,0+0i] Poles :[-615+0i,-0.406-0.606i,-0.406+0.606i,-1.226313+0j] K = 1500*616 V/(m.s^-1)

1.4. Signal processing

sisicepy.traces.signal.freq_filter.poles_zeros_filter(self, ZPG, waterlevel)

Pole and zeros filter (https://en.wikipedia.org/wiki/Pole%E2%80%93zero_plot)

Parameters:

  • self (xr.DataArray) – signal to filter

  • ZPG (list) – [zeros,poles,gains]

  • waterlevel (float)

sisicepy.traces.signal.freq_filter.trapezoid(self, band, delta)

Trapeziod band filter: .. line-block:

........Lower band...........Higher band........
........band[0]...............band[1]...........
...........|------------------|.................
..........||..................||................
.........|||..................|||...............
........||||..................||||..............
.......|||||..................|||||.............
......||||||..................||||||............
_____|||||||..................|||||||_______....
.......delta..................delta.............
Parameters:

  • self (xr.DataArray) – signal to filter with trapeziode band

  • band (list) – lower band and higher band ex: [4.5,90] Hz

  • delta (float) – delta for the begining and the end of the trapez

sisicepy.traces.signal.utils.remove_mean(self)

Remove mean from da signal

Parameters:

self (xr.DataArray) – signal

sisicepy.traces.signal.utils.remove_trend(self)

Remove linear trend from da signal

Parameters:

self (xr.DataArray) – signal

sisicepy.traces.signal.processing.PSD(self, **kwargs)

Wrapper of scipy.signal.welch function for Power Spectral Density. :param self: :type self: xr.DataArray

sisicepy.traces.signal.processing.spectrogram(self, win, t_start, dB=True, **kwargs)

Wrapper of scipy.signal.welch to performed PSD spectrogram

Parameters:

  • self (xr.DataArray)

  • win (int) – time length in seconds

  • t_start (str) – data for the start

  • dB (bool) – value in dB

2. xarray spectrogram accessor

class sisicepy.spectrogram.spectrogram.spectrogram(xarray_obj)

This is a classe to work on seismic spectrogram using xarray.

__init__(xarray_obj)

Constructor for traces.

Parameters:

xarray_obj (xr.DataArray)

__weakref__

list of weak references to the object

2.1. Processing

sisicepy.spectrogram.processing.to_energy(self, low_freq, high_freq)

Compute the energy between the given range of frequencies

Parameters:

  • low_freq (float)

  • high_freq (float)