pyrfu.pyrf.medfilt module#

pyrfu.pyrf.medfilt.medfilt(inp: DataArray, kernel_size: int | None = None) → DataArray[source]#

Applies a median filter over npts points to inp.

Parameters:

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

  • kernel_size (int, Optional) – Number of points of median filter. Default is a 3-point median filter.

Returns:

Time series of the median filtered input variable.

Return type:

DataArray

Raises:

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)