Evaluation of Different Filtering Methods Devoted to Magnetometer Data Denoising

: In this article, we describe a performance comparison conducted between several digital filters intended to mitigate the intrinsic noise observed in magnetometers. The considered filters were used to smooth the control signals derived from the magnetometers, which were present in an autonomous forestry machine. Three moving average FIR filters, based on rectangular Bartlett and Hanning windows, and an exponential moving average IIR filter were selected and analyzed. The trade-off between the noise reduction factor and the latency of the proposed filters was also investigated, taking into account the crucial importance of latency on real-time applications and control algorithms. Thus, a maximum latency value was used in the filter design procedure instead of the usual filter order. The experimental results and simulations show that the linear decay moving average (LDMA) and the raised cosine moving average (RCMA) filters outperformed the simple moving average (SMA) and the exponential moving average (EMA) in terms of noise reduction, for a fixed latency value, allowing a more accurate heading angle calculation and position control mechanism for autonomous and unmanned ground vehicles (UGVs).


Introduction
The assessment of a robot's orientation with respect to a fixed reference is conducted using the heading, also called the azimuth, a vector representing the robot's direction.This parameter is essential to allow the robot to move along a specific trajectory.It can be obtained in several ways using accelerometers, which measure the acceleration of a given object in motion by calculating its direction, and gyroscopes, which are designed to measure and maintain the orientation of an object relative to a reference point [1].Another way involves using magnetometer sensors that can measure the intensity and direction of the local magnetic field, providing crucial information about the orientation of a robot relative to the Earth's magnetic field [2].Magnetometers have been widely used to estimate the orientation of vehicles, such as cars, and drones, and are fundamentally useful when combined with other sensors, such as gyroscopes, accelerometers, and Global Navigation Satellite System (GNSS) receivers.For example, in the work conducted by Jau-Hsiung Wang et al., an algorithm was proposed using neural networks to estimate the heading from a magnetometer sensor.The proposed algorithm does not require declination information or the estimation of magnetometer bias or the scale factor [3]. Xufei Cui et al. proposed an estimation algorithm integrating an INS/magnetometer system based on a Kalman filter, to combine inertial sensor and magnetometer information with attitude information.The

•
Hard iron: this category encompasses magnetic fields that are static or only slightly variable over time and whose origin is from ferromagnetic materials [7]; Soft iron: refers to the magnetic field generated internally by the device itself and whose origin is from ferromagnetic materials [7].
However, these types of interferences can be mitigated through a calibration process that adjusts a set of parameters to compensate for their effect.In the present work, the heading was obtained based on data from the magnetometer to indicate the azimuth of a forestry robot.It is assumed in this case that the magnetometer was previously calibrated, however it continues to be affected by noise that causes uncertainties in the correct measurement of the heading.
Finite impulse response (FIR) filters, designed using the DFT window method, and infinite impulse response (IIR) filters, based on analog filter prototypes, were the initial candidates to reduce the noise in the magnetometer's signals.The preliminary results, obtained from simulations, confirm that the filtered signal latency is too large, precluding the real-time operation of unmanned ground vehicles (UGVs) in forest environments.Thus, several moving average filters were then analyzed and implemented to achieve the abovementioned goals.
This article is organized as follows: Section 2 describes several similar related works, devoted to the minimization of the noise present in magnetometers.Section 3 presents the magnetometer used in this work, as well as a detailed description of the proposed digital filters, in order to allow its successful replication in other technical solutions.In Section 4, the results obtained are presented and discussed.Moreover, an innovative approach that considers the maximum delay value as the main criteria to compare the noise reduction capabilities of the filters is validated.Section 5 presents the main conclusions derived from the findings in this study.

Related Work
This section presents the state of the art related to the noise characteristics of magnetometers and its sources, as well as the digital filtering techniques used in previous works to minimize noise from magnetometer sensors [8], which complement methods designed to prevent the effects of electromagnetic interference (EMI) [9].

Magnetometer Noise
Conventional magnetometer sensors can be affected by various independent sources, such as environmental noise in urban and industrial scenarios.These sensors are specifically designed to operate in the nanotesla range (which refers to the unit of measurement for magnetic field intensity); however, they face challenges due to environmental magnetic noise.In this way, magnetometers exhibit distinct types of noise, including 1/f noise and Johnson-Nyquist noise [8].Flicker noise, also known as 1/f noise, is more pronounced at low frequencies and it can occur from various sources, such as temperature variations or even mechanical vibrations.On the other hand, the thermal agitation of charge carriers in a conducting material produces Johnson-Nyquist white noise, a type of random noise, which is uniform across all frequencies.When evaluating the performance of magnetometers, these types of noise should be taken into consideration as they can affect the sensitivity and accuracy of magnetometers across various frequency ranges [10].

Advancements in Noise Minimization in Magnetometer Sensors
In recent years, advancements have been observed in the field of magnetometer sensor noise minimization.This area of research aims to improve the accuracy and reliability of data from these sensors, which play a fundamental role in a wide variety of applications, such as autonomous localization and navigation systems.
For example, Dressler M. et al. managed to reduce the white noise of a fluxgate magnetometer by introducing a JFET input stage integrated directly into the magnetometer sensor head to achieve high common-mode rejection and low circuit loading, resulting in noise reduction contributions [11].In [12], the authors used a method to cancel noise transmitted by a BLDC motor, thus not affecting the magnetometer sensor readings.An inverse filter was applied in the time domain, demonstrating a positive performance.Porang T. et al. conducted a comparative study between two filtering techniques, namely the Kalman filter and the Butterworth low-pass filter, applied to readings from a GPS and magnetometer sensor.These sensors were coupled to a tank-shaped robot, and the Kalman filter demonstrated greater stability in terms of the magnetometer data readings [13].Ning Z. et al. proposed an adaptive method for reducing noise in magnetometer sensor readings.The authors designed an adaptive filtering algorithm with multiple fading memory factors to perform offline noise reduction on the magnetometer's output data.The results showed good filtering capacity of the proposed method, successfully reducing noise [14].
However, the current literature lacks methods for reducing noise originating from magnetometer sensors using digital filtering techniques with minimum delay and distortion, which motivated the undertaking of this study.

The UGV Platform
The UGV platform used in this study was a MDB LV600 PRO (MDB, Fossacesia, Italy) multipurpose remote-controlled platform designed to operate in outdoor environments, as illustrated in Figure 1.
Electronics 2024, 13, x FOR PEER REVIEW 3 of 22 in a conducting material produces Johnson-Nyquist white noise, a type of random noise, which is uniform across all frequencies.When evaluating the performance of magnetometers, these types of noise should be taken into consideration as they can affect the sensitivity and accuracy of magnetometers across various frequency ranges [10].

Advancements in Noise Minimization in Magnetometer Sensors
In recent years, advancements have been observed in the field of magnetometer sensor noise minimization.This area of research aims to improve the accuracy and reliability of data from these sensors, which play a fundamental role in a wide variety of applications, such as autonomous localization and navigation systems.
For example, Dressler M. et al. managed to reduce the white noise of a fluxgate magnetometer by introducing a JFET input stage integrated directly into the magnetometer sensor head to achieve high common-mode rejection and low circuit loading, resulting in noise reduction contributions [11].In [12], the authors used a method to cancel noise transmitted by a BLDC motor, thus not affecting the magnetometer sensor readings.An inverse filter was applied in the time domain, demonstrating a positive performance.Porang T. et al. conducted a comparative study between two filtering techniques, namely the Kalman filter and the Butterworth low-pass filter, applied to readings from a GPS and magnetometer sensor.These sensors were coupled to a tank-shaped robot, and the Kalman filter demonstrated greater stability in terms of the magnetometer data readings [13].Ning Z. et al. proposed an adaptive method for reducing noise in magnetometer sensor readings.The authors designed an adaptive filtering algorithm with multiple fading memory factors to perform offline noise reduction on the magnetometer's output data.The results showed good filtering capacity of the proposed method, successfully reducing noise [14].
However, the current literature lacks methods for reducing noise originating from magnetometer sensors using digital filtering techniques with minimum delay and distortion, which motivated the undertaking of this study.

The UGV Platform
The UGV platform used in this study was a MDB LV600 PRO (MDB, Fossacesia, Italy) multipurpose remote-controlled platform designed to operate in outdoor environments, as illustrated in Figure 1.The aim was to make this robotic platform autonomous by integrating a variety of multimodal sensors.The sensors were integrated in a sensor module called Sentry (as shown in Figure 1 and attached on the top of the platform), and were responsible for the The aim was to make this robotic platform autonomous by integrating a variety of multimodal sensors.The sensors were integrated in a sensor module called Sentry (as shown in Figure 1 and attached on the top of the platform), and were responsible for the machine's perception and control, enabling autonomous localization and navigation [15,16].
One of the fundamental parameters for the autonomous navigation of a robot is its orientation, which specifically consists of a direction vector indicating the alignment of the robot, with north coinciding with the Y-axis and east coinciding with the X-axis, projected onto a two-dimensional map.This orientation was determined by a Bosch BMM150 (Bosch Sensortec GmbH, Reutlingen, Germany) magnetometer present in the Sentry module.The magnetometer sensor measures the magnetic fields using the Hall effect [17].The sensor characteristics are shown in Table 1.

Parameters Values
Measuring Range ±1300 µT (X and Y axes), ±2500 The sensor has a semiconductor plate with an electric current flowing through it.When this semiconductor is exposed to an external magnetic field, a Lorentz force is applied, resulting in the movement of electrons to one side of the plate, generating a potential difference called the Hall voltage.This voltage is directly proportional to the intensity of the applied magnetic field, and the sensor measures the Hall voltage and converts it into a magnetic field reading [17].The BMM150 sensor contains three axes with high precision for detecting weak and strong magnetic fields.These characteristics contribute to the precise estimates of the robot's state.Its structure can be observed in Figure 2.
Electronics 2024, 13, x FOR PEER REVIEW 4 of 22 machine's perception and control, enabling autonomous localization and navigation [15,16].One of the fundamental parameters for the autonomous navigation of a robot is its orientation, which specifically consists of a direction vector indicating the alignment of the robot, with north coinciding with the Y-axis and east coinciding with the X-axis, projected onto a two-dimensional map.This orientation was determined by a Bosch BMM150 (Bosch Sensortec GmbH, Reutlingen, Germany) magnetometer present in the Sentry module.The magnetometer sensor measures the magnetic fields using the Hall effect [17].The sensor characteristics are shown in Table 1.The sensor has a semiconductor plate with an electric current flowing through it.When this semiconductor is exposed to an external magnetic field, a Lorentz force is applied, resulting in the movement of electrons to one side of the plate, generating a potential difference called the Hall voltage.This voltage is directly proportional to the intensity of the applied magnetic field, and the sensor measures the Hall voltage and converts it into a magnetic field reading [17].The BMM150 sensor contains three axes with high precision for detecting weak and strong magnetic fields.These characteristics contribute to the precise estimates of the robot's state.Its structure can be observed in Figure 2. In this sensor, Hall effect components are responsible for detecting the magnetic fields around it and converting this information into an electrical signal.Subsequently, there are signal conditioning circuits that amplify/process the Hall sensor signal to ensure a more reliable reading.Then, the analog signal is converted into a digital signal through an ADC, which provides compatibility with microcontrollers and other electronic devices.The BMM150 sensor includes communication interfaces, such as the inter-integrated circuit (I2C) or serial peripheral interface (SPI) [17].This sensor is coupled to a module called Duro Inertial (Figure 3) from Swift Navigation, which in turn is attached to the top of the Sentry.In this sensor, Hall effect components are responsible for detecting the magnetic fields around it and converting this information into an electrical signal.Subsequently, there are signal conditioning circuits that amplify/process the Hall sensor signal to ensure a more reliable reading.Then, the analog signal is converted into a digital signal through an ADC, which provides compatibility with microcontrollers and other electronic devices.The BMM150 sensor includes communication interfaces, such as the inter-integrated circuit (I2C) or serial peripheral interface (SPI) [17].This sensor is coupled to a module called Duro Inertial (Figure 3) from Swift Navigation, which in turn is attached to the top of the Sentry.The data from the magnetometer sensor is read using the Robot Operating System (ROS) framework, which communicates via the Ethernet with the Duro Inertial sensor, as depicted in Figure 4.The corresponding source code is available at [18].Through the raw data obtained from the magnetometer, it is possible to manipulate the data to achieve the heading.Therefore, the following sections will discuss the performance of several moving average filtering techniques, with the aim of reducing the noise present in the signals.

Moving Average Filters
Moving average filters are widely used in engineering and economics to identify trends, smooth time-series data that exhibit short-term fluctuations, or to filter signals that are corrupted by AWGN noise.The most common versions are implemented using FIR filters, with the output, at each instant, being the weighted sum of a fixed number of past values of the input signal.This process is thereafter repeated, moving the window in order to calculate the mentioned weighted sum using new time-series data that includes new input samples, while discarding the oldest ones.
In a simple moving average (SMA) filter [19], all the input data are equally weighted to produce an output value, which is obtained by multiplying the last N input samples by a 1/N factor and discarding those that fall outside this range.The SMA filter impulse response is, therefore, given by: where N is the filter length or number of samples taken in the averaging process.The data from the magnetometer sensor is read using the Robot Operating System (ROS) framework, which communicates via the Ethernet with the Duro Inertial sensor, as depicted in Figure 4.The corresponding source code is available at [18].The data from the magnetometer sensor is read using the Robot Operating System (ROS) framework, which communicates via the Ethernet with the Duro Inertial sensor, as depicted in Figure 4.The corresponding source code is available at [18].Through the raw data obtained from the magnetometer, it is possible to manipulate the data to achieve the heading.Therefore, the following sections will discuss the performance of several moving average filtering techniques, with the aim of reducing the noise present in the signals.

Moving Average Filters
Moving average filters are widely used in engineering and economics to identify trends, smooth time-series data that exhibit short-term fluctuations, or to filter signals that are corrupted by AWGN noise.The most common versions are implemented using FIR filters, with the output, at each instant, being the weighted sum of a fixed number of past values of the input signal.This process is thereafter repeated, moving the window in order to calculate the mentioned weighted sum using new time-series data that includes new input samples, while discarding the oldest ones.
In a simple moving average (SMA) filter [19], all the input data are equally weighted to produce an output value, which is obtained by multiplying the last N input samples by a 1/N factor and discarding those that fall outside this range.The SMA filter impulse response is, therefore, given by: where N is the filter length or number of samples taken in the averaging process.Through the raw data obtained from the magnetometer, it is possible to manipulate the data to achieve the heading.Therefore, the following sections will discuss the performance of several moving average filtering techniques, with the aim of reducing the noise present in the signals.

Moving Average Filters
Moving average filters are widely used in engineering and economics to identify trends, smooth time-series data that exhibit short-term fluctuations, or to filter signals that are corrupted by AWGN noise.The most common versions are implemented using FIR filters, with the output, at each instant, being the weighted sum of a fixed number of past values of the input signal.This process is thereafter repeated, moving the window in order to calculate the mentioned weighted sum using new time-series data that includes new input samples, while discarding the oldest ones.
In a simple moving average (SMA) filter [19], all the input data are equally weighted to produce an output value, which is obtained by multiplying the last N input samples by a 1/N factor and discarding those that fall outside this range.The SMA filter impulse response is, therefore, given by: where N is the filter length or number of samples taken in the averaging process.
The linear decay moving average (LDMA) filter, as the name indicates, is a moving average FIR filter whose weighting factors will decay linearly.The most recent input signal samples will have the highest weights and the oldest ones the lowest.The filter impulse response is similar to the Bartlett window, but only the decaying part is considered, thus h LDMA is no longer even symmetrical around a central peak value.The LDMA filter impulse response is defined using expression [20]: The raised cosine moving average (RCMA) filter, as before, emphasizes the most recent input samples, which have the maximum weight, and undervalues the oldest ones.The proposed filter is based on the Hanning filter; however, only the part in which the weight values are progressively decreasing, according to a raised cosine shape, is considered (see Figure A1 of Appendix A), resulting in the RCMA impulse response: Finally, an IIR filter named the exponential moving average (EMA) is considered to diminish the noise variance in the acquired signal.The EMA filter, first introduced in [21], has an impulse response in which the past input values are geometrically reduced accordingly to their oldness.The EMA filter impulse response is given by: The β parameter, called the smoothing coefficient, is a positive real value less than one, which directly reflects the time constant that is intended to be implemented by the filter.In order to hold a unit DC gain, the filter constants α and β must fulfill the equation α = 1 − β.Moreover, Refs.[20,21] defines the relationship between the smoothing geometric constant β and the SMA filter length N, which results in a group delay or lag value equal to (N − 1)/2, which is the one associated with an even symmetric or an asymmetric impulse response.
In order to design the most suitable filter, aiming to mitigate the inherent noise present in the magnetometer's signals, it is necessary to analyze the frequency responses of the proposed filters.
The analysis of the SMA filter is a mature subject [19,22]; its transfer function is given by: The zeros of H SMA (z) are located over the unit circle of the z plane, starting at z = 1, being separated by angles equal to 2π/N.A pole-zero cancellation occurs at z = 1, allowing the desired unitary DC gain.The location of the zeros on the unit circle leads to frequencies with zero gain and distinctive side lobes in the magnitude of the frequency response, as depicted in Figure 5a The LDMA filter transfer function is described by: Analyzing the frequency response of the LDMA filter, it is possible to verify the existence of two poles at z = 1 that will be canceled by two zeros of HLDMA(z); double polezero cancellation, being the remaining zeros now located inside the unit circle, resulting in a smoother frequency response at higher frequencies, as showed in Figure 5.
The proposed RCMA filter transfer function is obtained by applying the definition: developing the equation and dividing the summation in two, and thereafter applying the infinite geometric sum formula and the Euler identity, we get: The RCMA filter transfer function is now rewritten considering that  = cos   ⁄ : The filter transfer function HRCMA(z) is given by: The filter transfer function has three poles located on the unit circle that will cancel with three zeros, the remaining zeros are located inside the unit circle, resulting in a smooth magnitude response and a nonlinear phase response, which contrasts with the characteristics and the zero's location on the Hanning window (Figure A4 of Appendix A).
The EMA filter transfer function is given by: which corresponds to a first order IIR filter having a single pole at z = β, being its location fundamental to the noise mitigation and the definition of the time constant values.The LDMA filter transfer function is described by: Analyzing the frequency response of the LDMA filter, it is possible to verify the existence of two poles at z = 1 that will be canceled by two zeros of H LDMA (z); double polezero cancellation, being the remaining zeros now located inside the unit circle, resulting in a smoother frequency response at higher frequencies, as showed in Figure 5.
The proposed RCMA filter transfer function is obtained by applying the definition: developing the equation and dividing the summation in two, and thereafter applying the infinite geometric sum formula and the Euler identity, we get: The RCMA filter transfer function is now rewritten considering that γ = cos(π/N): The filter transfer function H RCMA (z) is given by: The filter transfer function has three poles located on the unit circle that will cancel with three zeros, the remaining zeros are located inside the unit circle, resulting in a smooth magnitude response and a nonlinear phase response, which contrasts with the characteristics and the zero's location on the Hanning window (Figure A4 of Appendix A).
The EMA filter transfer function is given by: which corresponds to a first order IIR filter having a single pole at z = β, being its location fundamental to the noise mitigation and the definition of the time constant values.

Noise Reduction Analysis
The extraction of the desired signal from an acquired signal, which usually includes AWGN noise, is performed by means of an accurately designed filter H(z) that removes the unwanted noise leaving unaffected the desired signal features.Keeping in mind the peculiarities of the intended application, in this case the evaluation of the UGV heading based on coordinates taken from the magnetometer's sensors, the designed denoise filter must remove as much noise as possible, with the minimum delay and distortion, preserving a unitary gain at frequencies close to DC.
The noise reduction effectiveness of a digital linear filter is measured by the noisereduction factor (NRF), which quantifies the ratio of the output mean square noise values with respect to its input counterpart, defined by the following equations: The output and input mean square noise values of the noise components of the signals are defined by and v(n), being the output and input noise time sequences, respectively, with z(n) = h(n) * v(n), and also assuming linear filter usage that independently affects the desired signal and undesired noise that are superposed.Parseval's theorem is also used to evaluate the energy of the filtered noise signal in the time domain, which is equivalent to the energy evaluated in the frequency domain, integrating the square of the filter magnitude response.
The NRF expressions of the considered moving average filters, SMA, LDMA, RCMA, and EMA, shown in Table 2, were evaluated using (13).
Table 2. Noise-reduction factor expressions of the considered moving average filters.

Performance Metrics Expressions
Moving Average Filters SMA LDMA RCMA EMA Noise-Reduction Factor The noise-reduction factor of the proposed RCMA filter follows definition (13): By developing the binomial and dividing the summation into three parts, we obtain: From the individual sums, we obtain the NRF RCMA expression: The formulas were validated by simulations performed in MATLAB R2024a.In those simulations, it was considered that a deterministic signal (constant or sinusoidal) was corrupted by AWGN noise, with a particular variance value.The simulated NRF values were thereafter evaluated, dividing the filtered noise signal variance by the variance of the input noise for 10 8 length sequences.
From the results, shown in Figure 6, it is observed that the NRF value decreases with the increase in the filter length, improving the noise mitigation capabilities of the mentioned digital filters.For comparison purposes, an equivalent β parameter is used for the EMA filter that satisfies (5).Among the studied filters, the SMA filter is always the one that offers the best performance in terms of noise reduction for the same equivalent filter length.Keeping in mind that the planned solution will be implemented in the real-time trajectory control of an autonomous UGV, it is crucial that its position and direction values have been acquired and processed with minimum latency.The autonomous vehicle's heading values are obtained using the magnetometer's signals that are corrupted by noise, thus the noise minimizations filters must introduce minimum latency or lag.Thus, from now on, the performed analysis must consider the filter delay as the reference parameter instead of the filter order.Moreover, the literature has shown that a larger noise-reduction factor is always associated with an increase in the filter group delay or lag.
The application of linear-phase FIR filters, implemented by even or odd symmetrical windows, is a mature topic, resulting in a group delay value equal to  − 1 /2.The delay group is also independent of the frequency, being the derivative of the filter delay.That value is constant irrespective of the observed frequency.As presented above, the SMA filter has a symmetrical impulse response around a central point, satisfying the previous assumptions.The proposed LDMA and RCMA filter impulse responses are no longer even symmetrical, since they have been restricted to the descending part of the classical Bartlett and Hanning windows, respectively.
To evaluate the group delay (GD) at DC, the lag or latency of the proposed filters, we used a generic formula presented in [23].The group delay at DC is evaluated computing the following summation: where n is the index number and ℎ  is the filter impulse response.This expression is easily computed by simply multiplying, element by element, the filter impulse response vector, by a vector with linear increasing integer values (ramp signal).Appling (17) to the analyzed filters result in expressions, shown in Table 3, that can be used to evaluate the GD of the proposed filters based only on the filter length N. The GD expression for the RCMA filter approximates the exact result computed, as the filter length increases the error of that evaluation becomes negligible.Keeping in mind that the planned solution will be implemented in the real-time trajectory control of an autonomous UGV, it is crucial that its position and direction values have been acquired and processed with minimum latency.The autonomous vehicle's heading values are obtained using the magnetometer's signals that are corrupted by noise, thus the noise minimizations filters must introduce minimum latency or lag.Thus, from now on, the performed analysis must consider the filter delay as the reference parameter instead of the filter order.Moreover, the literature has shown that a larger noise-reduction factor is always associated with an increase in the filter group delay or lag.
The application of linear-phase FIR filters, implemented by even or odd symmetrical windows, is a mature topic, resulting in a group delay value equal to (N − 1)/2.The delay group is also independent of the frequency, being the derivative of the filter delay.That value is constant irrespective of the observed frequency.As presented above, the SMA filter has a symmetrical impulse response around a central point, satisfying the previous assumptions.The proposed LDMA and RCMA filter impulse responses are no longer even symmetrical, since they have been restricted to the descending part of the classical Bartlett and Hanning windows, respectively.
To evaluate the group delay (GD) at DC, the lag or latency of the proposed filters, we used a generic formula presented in [23].The group delay at DC is evaluated computing the following summation: where n is the index number and h(n) is the filter impulse response.This expression is easily computed by simply multiplying, element by element, the filter impulse response vector, by a vector with linear increasing integer values (ramp signal).Appling (17) to the analyzed filters result in expressions, shown in Table 3, that can be used to evaluate the GD of the proposed filters based only on the filter length N. The GD expression for the RCMA filter approximates the exact result computed, as the filter length increases the error of that evaluation becomes negligible.

Performance Metrics Expressions
Moving Average Filters From the results presented in Figure 7 it can be concluded that the GD value increases with the increase in the filter length.The RCMA filter produces the lowest latency, followed by the LDMA filter.Due to its even symmetry, the SMA filter produces a group delay value higher than the mentioned ones.It can also be concluded that the RCMA and LDMA filters have a lower latency value than that associated with the Bartlett and Hanning windows, for the same filter length.Figure 7 is a valuable resource in the assessment of the filter lengths that produce the same GD or latency value.As an example, for a target latency equal to 10 samples it is advisable to use: an SMA filter with N = 21; an LDMA filter with N = 31; an RCMA filter with N = 35; and an EMA filter with β equal to 10/11.
Electronics 2024, 13, x FOR PEER REVIEW 10 of 22 From the results presented in Figure 7 it can be concluded that the GD value increases with the increase in the filter length.The RCMA filter produces the lowest latency, followed by the LDMA filter.Due to its even symmetry, the SMA filter produces a group delay value higher than the mentioned ones.It can also be concluded that the RCMA and LDMA filters have a lower latency value than that associated with the Bartlett and Hanning windows, for the same filter length.Figure 7 is a valuable resource in the assessment of the filter lengths that produce the same GD or latency value.As an example, for a target latency equal to 10 samples it is advisable to use: an SMA filter with N = 21; an LDMA filter with N = 31; an RCMA filter with N = 35; and an EMA filter with β equal to 10/11.To validate the GD expressions, simulations were conducted to quantify the filtered signal delay regarding a sinusoidal input signal.Figure 8a,b presents the sinusoidal response of the considered filters, SMA, LDMA, RCMA and EMA, for lengths equal to 11 and 21, or equivalent values, respectively.From the obtained results it can be concluded that the RCMA filter presents the lowest latency at 119.38 ms, followed by the LDMA filter with a value of 133.33 ms and, finally, the SMA and EMA filters present a similar GD value of 200 ms.The same tendency is observed for filters with a length equal to 21.The filtered sinusoidal signal amplitudes are slightly different; this arises from the different gain values of the filter's magnitude frequency response at the signal frequency  = 0.05 rad/sample, as shown in Figure 5a,b, with the EMA filter magnitude being the lowest and the RCMA filter the highest one.To validate the GD expressions, simulations were conducted to quantify the filtered signal delay regarding a sinusoidal input signal.Figure 8a,b presents the sinusoidal response of the considered filters, SMA, LDMA, RCMA and EMA, for lengths equal to 11 and 21, or equivalent values, respectively.From the obtained results it can be concluded that the RCMA filter presents the lowest latency at 119.38 ms, followed by the LDMA filter with a value of 133.33 ms and, finally, the SMA and EMA filters present a similar GD value of 200 ms.The same tendency is observed for filters with a length equal to 21.The filtered sinusoidal signal amplitudes are slightly different; this arises from the different gain values of the filter's magnitude frequency response at the signal frequency ω = 0.05 rad/sample, as shown in Figure 5a,b, with the EMA filter magnitude being the lowest and the RCMA filter the highest one.
with a value of 133.33 ms and, finally, the SMA and EMA filters present a similar GD value of 200 ms.The same tendency is observed for filters with a length equal to 21.The filtered sinusoidal signal amplitudes are slightly different; this arises from the different gain values of the filter's magnitude frequency response at the signal frequency  = 0.05 rad/sample, as shown in Figure 5a,b, with the EMA filter magnitude being the lowest and the RCMA filter the highest one.The next stage of our analysis consists of a comparative study on the performance of moving average filters, assuming that the filters are of equal length.From the obtained results, as shown in Table 4, we can conclude that the SMA filtering technique outperforms the LDMA and RCMA techniques.The NRF value of the SMA filter is the lowest one among the considered filters; on the other hand, the latency or group delay value of the SMA filter is the largest.Although the EMA filter presents a similar result with respect to the SMA filter, its step response exhibits an undesirable behavior, underperforming in comparison to the other analyzed filters.In the conducted analysis, the magnetometer signals have been taken at a sampling rate of 25 Hz, i.e., a sampling time interval of 40 ms.Filters processing time delay, as already known, have a noticeable degrading effect on control system performance, namely those that involve close loop control [24].Thus, for the considered UGV position control mechanism, two maximum latency values were defined for the conducted analysis and simulations: 200 ms and 400 ms.For the abovementioned sampling time of 40 ms, the latency values correspond to 5 and 10 samples, respectively.Thus, knowing that the SMA filter latency is equal to (N SMA − 1)/2, then these values correspond to filter lengths equal to N SMA = 11 (200 ms latency) and N SMA = 21 (400 ms latency).
The next step is the LDMA filter N LDMA value assessment, which results in the same group delay value with respect to a particular SMA filter.Rearranging the expressions presented in Table 3, it is possible to evaluate the N LDMA , based on the SMA filter N SMA value, using (18).
The obtained N LDMA value is now used to evaluate the NRF of the LDMA filter, which produces the same latency as an SMA filter with length N SMA .From the obtained expression and results, we conclude that this innovative solution always outperforms the SMA filter, ensuring a given latency target value.
An SMA filter of length 11 reduces the noise by a factor of 0.09091 (1/11) and produces a latency equal to five samples, those values are compared with an LDMA filter of length 16, which produces the same latency and reduces the noise by a factor of 0.08088 (11/136), as can be seen in Table 5.As explained before, the SNR is inversely proportional to the NRF, thus the proposed LDMA filter usage outperforms the SMA filter by 0.50753 dB.For longer filter lengths that value converges to 0.51153 dB.Moreover, the observed trend in the results is independent of the latency value, as can be seen for the LDMA filter of length 31, which outperforms the SMA filter of length 21, both filters having a group delay of 10 samples.Despite the lack of a closed formula to assess the group delay of the RCMA filter, it is possible to observe from the results presented in Table 5 that an RCMA filter with N = 17 will reduce the noise by a factor of 0.08488, with a latency value of 190.60 ms.This result is better than the one obtained for the SMA filter, but slightly worse than the corresponding one attained for the LDMA filter.If a 1.2% increase in the latency value to 202.48 ms is allowed, then an N = 18 filter with an NRF value of 0.08033 can be used, which performs better than the equivalent LDMA filter of length 16.The same trend is observed in the filter's performance, if a 400 ms group delay was considered as a design constraint value.
The obtained results can also be verified in the frequency domain.Figure 9 presents the frequency response of the filters for an equal group delay, with it being possible to observe the lowest magnitude of the LDMA and RCMA filters' frequency responses for frequencies higher than that associated with the first null of the SMA filter's response.That observation explains the improved performance in regard to the noise reduction goal of those filters with different lengths, while maintaining the latency constant.A non-significant advancement was obtained from the EMA filter usage.
the frequency response of the filters for an equal group delay, with it being possible to observe the lowest magnitude of the LDMA and RCMA filters' frequency responses for frequencies higher than that associated with the first null of the SMA filter's response.That observation explains the improved performance in regard to the noise reduction goal of those filters with different lengths, while maintaining the latency constant.A non-significant advancement was obtained from the EMA filter usage.Now, to validate the new proposed method, a transient analysis of the considered filters is performed [25].Figure 10a shows the step response of the filters, assuming that Now, to validate the new proposed method, a transient analysis of the considered filters is performed [25].Figure 10a shows the step response of the filters, assuming that they have a length of 11.From the results obtained, it is possible to conclude that the LDMA and RCMA filters respond faster than the SMA filter, this is consistent with the lower GD values observed; however, all three filters, SMA, LDMA and RCMA, reach the final stage value after 11 samples only, due to the inherent DC unit normalization [20,25].The EMA filter step response is slower and takes a significantly extended time period to reach the final value of the step.In that analysis, a β = 5/6 value was considered for the EMA filter, which is equivalent to an SMA filter with N = 11, according to (5).  Figure 10b presents the step response of the filters, with different length values, which produce the same GD or latency value.Once again, all the filters respond faster than the SMA filter until a particular point, around 80% of its final value.Thereafter, the LDMA and RCMA filter responses are slower, since the weights associated with the oldest values of their impulsive response are significantly smaller, contributing less to the output value.Moreover, only after 16 and 18 times, the sampling interval, respectively, did the mentioned LDMA and RCMA filters reach the exact step value.

Results
This section presents the experimental results obtained from the new proposed approach, which considers latency as the main design criterion for filters devoted to reduc- Figure 10b presents the step response of the filters, with different length values, which produce the same GD or latency value.Once again, all the filters respond faster than the SMA filter until a particular point, around 80% of its final value.Thereafter, the LDMA and RCMA filter responses are slower, since the weights associated with the oldest values of their impulsive response are significantly smaller, contributing less to the output value.
Moreover, only after 16 and 18 times, the sampling interval, respectively, did the mentioned LDMA and RCMA filters reach the exact step value.

Results
This section presents the experimental results obtained from the new proposed approach, which considers latency as the main design criterion for filters devoted to reducing noise in magnetometers.The experiment involves collecting data in static conditions, by attaching a Duro Inertial sensor to a wall at the top of a building.
Since the position of the device does not change, the result obtained is a reference signal superimposed on the noise (the X, Y, and Z signals acquired from the magnetometer are corrupted by noise).The following sections in this paper are dedicated to the statistical analysis of the raw and filtered signals, with the aim of determining which filtering technique and filtering order provides the best performance in terms of noise reduction.

Data Collection System
As mentioned above, data collection from the magnetometer sensor was conducted under stationary conditions.The Duro Inertial sensor was placed on a platform and then fixed to a wall at the top of a building forming part of the Association for the Development of Industrial Aerodynamics, located on the campus of the Faculty of Science and Technology at the University of Coimbra, Portugal, as shown in Figure 11.The magnetometer's signals were then acquired at a sampling rate of 25 Hz.Those signals were filtered, in real time, using the proposed moving average filters.The filters were implemented on a processing unit equipped with an Intel ® Core™ i7-5500U CPU (Intel Corporation, Santa Clara, CA, USA.), operating at 2.4 GHz, and with 4 GB of DDR3L SO-DIMM RAM (Kingston Technology, Valley, CA, USA).The filter source code was developed in Python, using the NumPy library.The implemented solution uses floatingpoint arithmetic.
After collecting the raw and filtered signals, statistical analysis was carried out which included assessing the mean, standard deviation (SD), and mean absolute error (MAE) of the raw and filtered time-series data.
The MAE is a metric that calculates the arithmetic mean of the absolute differences between the values predicted by a model and the actual values of a given data set.The lower the MAE value, the better the model's performance, indicating that the predictions are close to the actual values.In summary, the MAE is defined as follows [26]: where, as follows: • N is the number of samples; • y is the observed value for each sample; •  is the value predicted by the model for each sample.
The MAE metric was selected, among others, because it offers the better resistance The magnetometer's signals were then acquired at a sampling rate of 25 Hz.Those signals were filtered, in real time, using the proposed moving average filters.The filters were implemented on a processing unit equipped with an Intel ® Core™ i7-5500U CPU (Intel Corporation, Santa Clara, CA, USA), operating at 2.4 GHz, and with 4 GB of DDR3L SO-DIMM RAM (Kingston Technology, Valley, CA, USA).The filter source code was developed in Python, using the NumPy library.The implemented solution uses floatingpoint arithmetic.
After collecting the raw and filtered signals, statistical analysis was carried out which included assessing the mean, standard deviation (SD), and mean absolute error (MAE) of the raw and filtered time-series data.
The MAE is a metric that calculates the arithmetic mean of the absolute differences between the values predicted by a model and the actual values of a given data set.The lower the MAE value, the better the model's performance, indicating that the predictions are close to the actual values.In summary, the MAE is defined as follows [26]: where, as follows: • N is the number of samples; • y is the observed value for each sample; • ŷi is the value predicted by the model for each sample.
The MAE metric was selected, among others, because it offers the better resistance against distortions introduced by outliers present in the signal time series, treating all errors equally.

Experimental Results
From now on, our analysis is devoted to the evaluation of the heading direction, i.e., the orientation of the robot relative to the Earth's magnetic field using the Duro Inertial, in a fixed position, based on the magnetometer's signals.
Figure 12 illustrates the three raw data signals acquired from the magnetometer corresponding to the X, Y, and Z axes.By analyzing the acquired data, it is possible to observe the quantized noise superimposed on the signals that will be processed.From the data, it is also possible to state that the highest magnitude values are associated with the X and Z axes, with the Y-axis measured values being residual.Moreover, from the statistical analysis of the raw data obtained from the magnetometer, it is possible to conclude that the amplitude distribution of the sampled signals related to the magnetic field in the X, Y, and Z axes, under static conditions, roughly follows a Gaussian distribution.The raw data acquired from the Duro Inertial magnetometer associated with the Xaxis was then filtered using two SMA filters with lengths equal to 11 and 21, and two LDMA filters with lengths equal to 16 and 31, the results of which are shown in Figure 13a and Figure 13b, respectively.The results show that increasing the length of the filters always results in a greater capacity to reduce noise and signal fluctuations.For filters with the same latency, the LDMA technique performs better than the SMA technique.
(a) (b) The raw data acquired from the Duro Inertial magnetometer associated with the X-axis was then filtered using two SMA filters with lengths equal to 11 and 21, and two LDMA filters with lengths equal to 16 and 31, the results of which are shown in Figure 13a and Figure 13b, respectively.The results show that increasing the length of the filters always results in a greater capacity to reduce noise and signal fluctuations.For filters with the same latency, the LDMA technique performs better than the SMA technique.
The overall noise reduction performance of the filters was also evaluated based on descriptive statistics of the data acquired by the magnetometer before and after filtering.From the results presented in Table 6, it can be concluded that all proposed filters significantly reduce the noise in the magnetometer signals, a conclusion based on the reduction in the standard deviation value observed in the filtered signals.
The raw data acquired from the Duro Inertial magnetometer associated with the Xaxis was then filtered using two SMA filters with lengths equal to 11 and 21, and two LDMA filters with lengths equal to 16 and 31, the results of which are shown in Figure 13a and Figure 13b, respectively.The results show that increasing the length of the filters always results in a greater capacity to reduce noise and signal fluctuations.For filters with the same latency, the LDMA technique performs better than the SMA technique.The overall noise reduction performance of the filters was also evaluated based on descriptive statistics of the data acquired by the magnetometer before and after filtering.From the results presented in Table 6, it can be concluded that all proposed filters significantly reduce the noise in the magnetometer signals, a conclusion based on the reduction in the standard deviation value observed in the filtered signals.
Among them, the LDMA filter provides the greatest reduction in the value of the standard deviation on the X and Z axes, closely followed by the RCMA filter.The above assumption that the raw signals follow a Gaussian distribution makes it possible to calculate the noise reduction factor using the variance values of the filtered   2 and raw signals    2 , the square values of the standard deviations presented.Considering the Zaxis and an LDMA filter with N = 16, the approximate NRF is equal to 0.2476 or 0.514 1.033 ⁄ , which is higher than the theoretical value of 0.081.Among them, the LDMA filter provides the greatest reduction in the value of the standard deviation on the X and Z axes, closely followed by the RCMA filter.The above assumption that the raw signals follow a Gaussian distribution makes it possible to calculate the noise reduction factor using the variance values of the filtered σ 2 f iltered and raw signals σ 2  raw , the square values of the standard deviations presented.Considering the Z-axis and an LDMA filter with N = 16, the approximate NRF is equal to 0.2476 or 0.514 2 /1.033 2 , which is higher than the theoretical value of 0.081.The results presented in Table 6 are now compared with other filtering methods available in the literature [14] to reduce the magnetometer noise.
Table 7 shows the NRF values calculated from the standard deviation values of the filtered σ f iltered and raw signals σ raw , using Equation (13).For the wavelet denoising, EKF, and adaptive filter techniques, the NRF values were calculated based on the data available in Table V of [14]; for the moving average filters studied in this article we used the data shown in Table 6.
From the results presented in Table 7, we conclude that the LDMA filter with N = 31 and the RCMA filter with N = 35 outperform the wavelet denoising technique, behave slightly better than the extended Kalman filter (EKF), and perform worse than the adaptive filter, with larger NRF values.Note that [14] uses the RM1300 sensor, while we have used the BMM150 sensor.In addition, compliance with the latency limit is not guaranteed for the analyzed adaptive filter.The final step in our analysis is devoted to the heading signal.Firstly, we examine the raw heading values evaluated from the noisy magnetometer signals, which also exhibit a quantized nature.The raw unprocessed heading data is depicted in the upper graph in Figure 14, with the filtered version shown in the lower graph; once again, an SMA filter of length 11 was used to smooth the original heading signal.Table 8 presents the descriptive statistics of the heading signal before and after filtering, using different moving average filters and lengths.The mean, standard deviation, and MAE values of those signals are used to access the filtering technique performance.The results exhibit a notable reduction in the filtered heading signal standard deviation, with respect to the raw heading data for all the considered filters.As an example, considering a latency limit of five samples (or 200 ms), the standard deviation σ values were equal to 0.8787°, 0.8303°, 0.8278°, and 0.8695° from an SMA filter with N = 11; an LDMA filter with N = 16; an RCMA filter with N = 18; and an EMA filter with β = 5/6, respectively.Those values must be compared with the raw heading data standard deviation equal to 1.8007°.For filters with a longer length, the noise reduction capabilities increase further.
The results show that the RCMA filter has the best noise attenuation capacity for a latency value of 5 or 10 samples and that the LDMA filter presents similar performance.Table 8 presents the descriptive statistics of the heading signal before and after filtering, using different moving average filters and lengths.The mean, standard deviation, and MAE values of those signals are used to access the filtering technique performance.The results exhibit a notable reduction in the filtered heading signal standard deviation, with respect to the raw heading data for all the considered filters.As an example, considering a latency limit of five samples (or 200 ms), the standard deviation σ values were equal to 0.8787 • , 0.8303 • , 0.8278 • , and 0.8695 • from an SMA filter with N = 11; an LDMA filter with N = 16; an RCMA filter with N = 18; and an EMA filter with β = 5/6, respectively.Those values must be compared with the raw heading data standard deviation equal to 1.8007 • .For filters with a longer length, the noise reduction capabilities increase further.
The results show that the RCMA filter has the best noise attenuation capacity for a latency value of 5 or 10 samples and that the LDMA filter presents similar performance.
Regarding heading noise reduction, we claim that for a latency threshold of 10 samples, all the proposed filtering techniques, SMA of length 21, LDMA of length 31, RCMA of length 35, and EMA with β equal to 10/11, outperform the adaptive filtering presented in [14].However, if that latency limit reduces to five samples, the adaptive filtering performs better.The noise-reduction factor for the adaptive filtering presented in [14] was evaluated by dividing the variance of the filtered heading signal by the variance of the raw heading signal; the result obtained is equal to σ 2 f iltered /σ 2 raw = 0.1985.The variance values were obtained by squaring the standard deviation values of the heading signals before and after filtering, which were equal to σ f iltered = 0.5387 • and σ raw = 1.2089 • , respectively.Those values were obtained from Table VII of [14].
A different approach was proposed and analyzed, which involved filtering the heading signal instead of filtering the magnetometer signals, used to calculate the heading.This approach, which swaps the order in which the filtering and heading calculations were performed, results in a slightly lower-quality performance.

Conclusions and Future Work
This article presents a performance evaluation of different digital filtering techniques to minimize the noise observed in magnetometers.These sensors are used in diverse applications, such as mobile robotics and autonomous UGVs, where they provide crucial orientation information needed for navigation.However, these sensors are sensitive and can be affected by several types of noise, including flicker noise and Johnson-Nyquist noise.
Four moving average filters, SMA, LDMA, RCMA, and EMA, were thereafter selected and analyzed to smooth the magnetometer's data, which was corrupted with different types of noise.The RCMA filter is an innovative technique, proposed by the authors, which only considers the descending part of the Hanning window impulse response.A comparative analysis of these filtering techniques is available in Appendix A.
The results indicate that the SMA filtering technique outperforms the LDMA and RCMA techniques in terms of noise reduction, if filters of equal length are considered; although the SMA filter presents a higher latency, when compared with the LDMA and RCMA proposals.
Considering the impact of latency on autonomous navigation systems, a maximum delay value was considered in the filter design procedure, instead of the usual filter order value.Under the same latency assumption, the LDMA and RCMA filters outperformed the SMA and the EMA filters, with a similar transient step response.
All the techniques were then implemented using a microprocessor and an experimental procedure was conducted, using a Duro Inertial sensor.From the experimental results obtained by placing the Duro Inertial in a fixed position, it can be seen that the noise in the magnetometer is significantly reduced by filtering.The sensor heading values calculated from the magnetometer's filtered data also show a reduction in the standard deviation values, which is consistent with a reduction in noise.The RCMA filter stands out, with the lowest standard deviation and MAE values of all the techniques analyzed, independently of the latency value: 200 ms or 400 ms.The proposed approach, which defines the allowed latency and then selects the technique that minimizes noise mitigation, offers significant The DC normalization performed on the Hanning window modifies the gain of the transfer function only, while maintaining all its original characteristics, namely: pole-zero patterns, magnitude and phase responses, and group delay and impulse response shape.The example used to illustrate the differences between the aforementioned techniques considers an RCMA filter and a normalized Hanning window, both with N equal to 11.The impulse responses of the above filters are shown in Figure A1.The group delay of the Hanning window, shown in Figure A3b, is constant for all frequencies, because the filter has a linear phase.The group delay of the RCMA filter decreases significantly as the frequency increases, with the highest value occurring at DC, as illustrated in Figure A3a.The group delay of the Hanning window, shown in Figure A3b, is constant for all frequencies, because the filter has a linear phase.The group delay of the RCMA filter decreases significantly as the frequency increases, with the highest value occurring at DC, as illustrated in Figure A3a.The group delay of the Hanning window, shown in Figure A3b, is constant for all frequencies, because the filter has a linear phase.The group delay of the RCMA filter decreases significantly as the frequency increases, with the highest value occurring at DC, as illustrated in Figure A3a.

Figure 4 .
Figure 4. Interface to collect data from magnetometer sensor.

Figure 4 .
Figure 4. Interface to collect data from magnetometer sensor.

Figure 4 .
Figure 4. Interface to collect data from magnetometer sensor.

Electronics 2024 , 22 Figure 6 .
Figure 6.Noise-reduction factor value as a function of the filter length.

Figure 6 .
Figure 6.Noise-reduction factor value as a function of the filter length.

Figure 7 .
Figure 7. Group delay value as a function of the filter length.

Figure 7 .
Figure 7. Group delay value as a function of the filter length.

Figure 8 .
Figure 8. Sinusoidal response of moving average filters: (a) group delay for filters when N = 11; (b) group delay for filters when N = 21.
Electronics 2024, 13, x FOR PEER REVIEW 13 of 22 they have a length of 11.From the results obtained, it is possible to conclude that the LDMA and RCMA filters respond faster than the SMA filter, this is consistent with the lower GD values observed; however, all three filters, SMA, LDMA and RCMA, reach the final stage value after 11 samples only, due to the inherent DC unit normalization [20,25].The EMA filter step response is slower and takes a significantly extended time period to reach the final value of the step.In that analysis, a β = 5/6 value was considered for the EMA filter, which is equivalent to an SMA filter with N = 11, according to (5).

Figure 10 .
Figure 10.Transient step response of moving average filters: (a) filters with N = 11; (b) filters with lengths that correspond to an equal group delay.

Figure 10 .
Figure 10.Transient step response of moving average filters: (a) filters with N = 11; (b) filters with lengths that correspond to an equal group delay.

22 Figure 11 .
Figure 11.Duro Inertial sensor on a platform, used for data collection.

Figure 11 .
Figure 11.Duro Inertial sensor on a platform, used for data collection.

22 Figure 12 .
Figure 12.Unfiltered magnetometer signals acquired in the X, Y, and Z axes.

Figure 12 .
Figure 12.Unfiltered magnetometer signals acquired in the X, Y, and Z axes.

Figure 13 .
Figure 13.Magnetometer X-axis raw and filtered signals: (a) SMA filters with N equal to 11 and 21; (b) LDMA filters with N equal to 16 and 31.

Figure 13 .
Figure 13.Magnetometer X-axis raw and filtered signals: (a) SMA filters with N equal to 11 and 21; (b) LDMA filters with N equal to 16 and 31.

Electronics 2024 , 22 Figure 14 .
Figure 14.Heading signals: raw data and filtered data using an SMA filter.

Figure 14 .
Figure 14.Heading signals: raw data and filtered data using an SMA filter.

Figure
Figure A2 illustrates the magnitude and phase responses of the RCMA filter and Hanning window for N equal to 11.

Figure
Figure A4a,b depicts the pole-zero patterns of the considered RCMA filter and normalized Hanning window, respectively.

Figure
Figure A4a,b depicts the pole-zero patterns of the considered RCMA filter and normalized Hanning window, respectively.

FigureFigure A4 .
Figure A4a,b depicts the pole-zero patterns of the considered RCMA filter and normalized Hanning window, respectively.Electronics 2024, 13, x FOR PEER REVIEW 21 of 22

Table 3 .
Group delay expressions of the considered filters.

Table 3 .
Group delay expressions of the considered filters.

Table 4 .
Performance metrics for equal length moving average filters.

Table 5 .
Moving average filters performance metrics for target group delay values.

Table 6 .
Descriptive statistics of the BMM 150 magnetometer data before and after filtering.

Table 7 .
[14]arison between several filtering methods[14]and the studied moving average filters, using the signal noise variance reduction metric.

Table 8 .
Descriptive statistics of the heading signal before and after filtering.

Table 8 .
Descriptive statistics of the heading signal before and after filtering.