pyrfu.pyrf.medfilt module#

pyrfu.pyrf.medfilt.medfilt(inp, n_pts: int = 11)[source]#

Applies a median filter over npts points to inp.

Parameters:

  • inp (xarray.DataArray) – Time series of the input variable.

  • n_pts (int, Optional) – Number of points of median filter.

Returns:

out – Time series of the median filtered input variable.

Return type:

xarray.DataArray

Examples

>>> import numpy
>>> from pyrfu import mms, pyrf

Time interval

>>> tint = ["2019-09-14T07:54:00.000", "2019-09-14T08:11:00.000"]

Spacecraft indices

>>> mms_list = numpy.arange(1,5)

Load magnetic field and electric field

>>> r_mms, b_mms = [[] * 4 for _ in range(2)]
>>> for mms_id in range(1, 5):
>>>         r_mms.append(mms.get_data("R_gse", tint, mms_id))
>>>         b_mms.append(mms.get_data("B_gse_fgm_srvy_l2", tint, mms_id))
>>>

Compute current density, etc

>>> j_xyz, _, b_xyz, _, _, _ = pyrf.c_4_j(r_mms, b_mms)

Get J sampling frequency

>>> fs = pyrf.calc_fs(j_xyz)

Median filter over 1s

>>> j_xyz = pyrf.medfilt(j_xyz,fs)