watex.methods.em.MT.get_ss_correction_factors#

MT.get_ss_correction_factors(station, r=100.0, nfreq=21, skipfreq=7, tol=0.12, bounds_error=True)[source]#

Compute the static shift correction factor from a station using a spatial median filter.

This will find those station within the given radius (meters). Then it will find the median static shift for the x and y modes and remove it, given that it is larger than the shift tolerance away from 1.

Parameters:

station (str, int,) – The index of station to compute the static shift factors. If the station name is passed as string object, it should include the position number.
r (float, default=1000.) – radius to look for nearby stations, in meters.
nfreq (int, default=21) – number of frequencies calculate the median static shift. This is assuming the first frequency is the highest frequency. Cause usually highest frequencies are sampling a 1D earth.
skipfreq (int, default=5) – number of frequencies to skip from the highest frequency. Sometimes the highest frequencies are not reliable due to noise or low signal in the AMT deadband. This allows you to skip those frequencies.
bounds_error (bool, default=True) –
Check whether the frequency for interpolation is within the frequency range. Otherwise, raises an error.

New in version 0.2.8.
tol (float, default=0.12) – Tolerance on the median static shift correction. If the data is noisy the correction factor can be biased away from 1. Therefore the shift_tol is used to stop that bias. If 1-tol < correction < 1+tol then the correction factor is set to 1

Returns:

(sx_x, ss_y) – static shift corrections factor for x and y modes

Return type:

(float, float)

Examples

>>> import watex
>>> from watex.methods import ZC
>>> edi_sample = watex.fetch_data ('edis', samples =17 ,
                                   return_data =True )
>>> zo = ZC ().fit(edi_sample).get_ss_correction_factors (station =0 )
Out[16]: (1.5522030221266643, 0.742682340427651)