Next Article in Journal
An Optimization Method for Non-IID Federated Learning Based on Deep Reinforcement Learning
Next Article in Special Issue
Advantages of Bistable Microwires in Digital Signal Processing
Previous Article in Journal
Guided Wave Characteristic Research and Probabilistic Crack Evaluation in Complex Multi-Layer Stringer Splice Joint Structure
Previous Article in Special Issue
Use of Different Types of Magnetic Field Sensors in Diagnosing the State of Ferromagnetic Elements Based on Residual Magnetic Field Measurements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on a Time Difference Processing Method for RTD-Fluxgate Data Based on the Combination of the Mahalanobis Distance and Group Covariance

College of Computer Science and Technology, Beihua University, No. 3999 East Binjiang Road, Jilin 132013, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(22), 9223; https://doi.org/10.3390/s23229223
Submission received: 19 October 2023 / Revised: 12 November 2023 / Accepted: 14 November 2023 / Published: 16 November 2023
(This article belongs to the Special Issue Advances and Applications of Magnetic Sensors)

Abstract

:
During the measurement of magnetic fields, Residence Time Difference (RTD)-fluxgate sensors suffer from abnormal time difference jumps due to the random interference of magnetic core noise and environmental noise, which results in gross errors. This situation restricts the improvement of sensor accuracy and stability. In order to solve the above problems efficiently, a time difference gross error processing method based on the combination of the Mahalanobis distance (MD) and group covariance is presented in this paper, and the processing effects of different methods are compared and analyzed. The results of the simulation and experiment indicate that the proposed method is more advantageous in identifying the gross error in time difference. The signal-to-noise ratio for the time difference is improved by about 34 times, while the fluctuation of the Negative Magnetic Saturation Time (NMST) ΔTNMST is reduced by 95.402%, which significantly reduces the fluctuation of time difference and effectively improves the accuracy and stability of the sensor.

1. Introduction

In recent years, Unmanned Aerial Vehicle (UAV) aeromagnetic measurement as a new aeromagnetic technology has become a hot research topic [1]. Due to the characteristics of high sensitivity, small size, low power consumption, etc., the fluxgate sensor contributes to the trend of portable measurement with UAV as a carrier [2]. In particular, the RTD-fluxgate can avoid odd harmonic interference and allows for complete symmetry of the two-axis structure of the sensing unit in terms of the measurement principle [3,4,5]. In order to better meet the requirements for UAV aeromagnetic measurement, further improving the accuracy and stability of an RTD-fluxgate has become a key research topic [6,7,8].
When the RTD-fluxgate is working, the magnetic core produces random magnetic noise under the action of the excitation magnetic field [9,10]. The detection circuit and time difference during transmission are easily affected by environmental noise. The above situation causes abnormal time difference jumps, reducing the stability of the RTD-fluxgate [11,12]. Therefore, it is necessary to estimate and deal with the gross error in time difference. The Laida criterion can be used if the processing method adopted includes discriminating and eliminating gross errors. However, this method is only suitable for cases with small measurement times and low processing accuracy [13]. When the data contain only a single outlier, the Grubbs criterion can be used more easily and effectively, but it cannot deal with a large number of outliers [14,15]. When the comparatively large amount of data does not follow a normal distribution, using the Chauville criterion reduces the accuracy of the measurement [16]. This shows that the scope of application for different criteria is diverse. The direct elimination of gross errors can lead to important data information being removed. If we use the moving average filtering method, the random fluctuation caused by noise can be effectively suppressed, but the selection of window parameters directly affects the smoothing effect on the data. When larger window parameters are selected, the determinacy of the high-frequency variation is weakened on average. However, when smaller window parameters are selected, the low-frequency random fluctuation is not averaged [17]. At the same time, the amount of data is reduced due to the added windows processing. It can be seen that the methods above have limitations when dealing with the gross error in time difference.
Because the MD is sensitive to data outliers, it is suitable as a statistic for measuring outliers [18,19]. The MD of the time difference is regarded as the standard, so the gross error in the time difference can be accurately estimated. In addition, covariance is the measure of the degree of deviation from its mean. Therefore, in order to eliminate the gross errors in real-time dynamic magnetic field measurement using RTD-fluxgate, this article proposes a method for processing time difference gross errors by combining MD with group covariance as weights cooperating with mean values. This method can reduce the time difference fluctuation caused by random noise, improve the stability of RTD-fluxgate, and better process the time difference dynamically in real time.

2. Structure and Working Principle of RTD-Fluxgate

The RTD-fluxgate sensor is mainly composed of a sensing unit and a detection circuit, as is shown in Figure 1a. When it is working, the periodic alternating excitation current is input to both ends of the excitation coil of the sensing unit in order to generate an excitation magnetic field, where the magnetic core is magnetized to a bidirectional supersaturated state. As is shown in Figure 1b, after the excitation magnetic field is modulated by the hysteresis loop of the magnetic core, a pulse shape induction signal is generated through the induction coil. We can obtain the target magnetic field Hx by determining the time difference between the positive and negative pulses of the induction signal relative to the magnetic saturation [20,21].
In order to minimize the influence of noise on the detection of and reduce the uncertainty in the time difference, we use the excitation signal and the output pulse signal, combined, to find the Negative Magnetic Saturation Time (NMST) ∆TNMST, which is the quantity used to measure the target magnetic field Hx [22,23,24]. The schematic diagram for the NMST reading technology is shown in Figure 2. When the magnetic core becomes saturated, an output pulse signal is generated on the induction coil; the transition time when the excitation signal amplitude is zero is used as the reference time, tT; and when the applied magnetic field exceeds the coercive field, −Hc, the induced voltage produces a negative pulse at tP.
When the trapezoidal excitation field He(t) is used, assuming the total time of the slope within a cycle is t , the stable state time is t . Te is the excitation period, and f is the excitation frequency, so t + t  = Te. The slope of the trapezoidal excitation is ±a, and the stable amplitude is ±Hm. The expression is as follows:
H e ( t ) = a t               N T e t 4 < t < N T e + t 4 H m             N T e + t 4 < t < N T e + T e 2 t 4 a t             N T e + T e 2 t 4 < t < N T e + T e 2 + t 4 H m           N T e + T e 2 + t 4 < t < N T e + T e t 4
The relationship between the excitation field He(t) and the target field Hx is as follows:
t T : α t T = 0
t P : H x α ( t P T e 2 ) = H c
Deduced using Equations (2) and (3):
t T = 0
t P = ( H c + H x ) α + T e 2
The time difference ∆TNMST between tP and tT is defined by the Negative Magnetic Saturation Time, which is given by the following:
Δ T NMST = t P t T
The expression of the relationship between ∆TNMST and Hx is as follows:
Δ T NMST = t P t T = ( H c + H x ) 4 H m f + T e 2
The sensitivity SNMST expression found using the NMST reading strategy is as follows:
S NMST = Δ T NMST H x = 1 α = T e 4 H m = 1 4 H m f

3. Discriminating and Processing Methods for Gross Error in Time Difference

The detection circuit in an RTD-fluxgate sensor adopts certain measures to suppress fluctuation and utilizes the NMST readout strategy to reduce the influence of induced signals. However, due to random external interference and other comprehensive factors from the environment, the original time difference still has certain fluctuation, and even some individual values fluctuate greatly, which causes a decrease in sensor stability. Therefore, it is necessary to estimate the gross error accurately and replace it with the effective value after time difference processing. Due to the sensitivity of the MD to anomalous values, the MD is regarded as a test statistic for the time difference; it can accurately estimate whether the time difference is interfered by noise or not. The covariance of the time difference can measure the deviation degree from its mean. The MD of the time difference data combined with the covariance of the corresponding group as the weight cooperate with the mean value replaces the gross time difference and reduces the impact of random fluctuation.

3.1. Time Difference Gross Error Discrimination Method Based on the Mahalanobis Distance

The Mahalanobis distance is the closest distance between a single sample and the “center of gravity” of a sample set [25]; this can help calculate the similarity of two unknown sample sets and overcome the defect of Euclidean distance which lacks the overall sample to affect the distance. Considering the relationship between various features, the MD is sensitive to abnormal values [26,27,28]. For a sample set xi(=1,… n.), the Dm expression of the MD between two samples, xi and xj, is expressed as follows:
D i = ( x i x j ) T · Σ 1 · ( x i x j )
Because of the sensitivity of the MD to anomalous values, it is suitable to be used as a statistic to measure anomalous values. The MD from the time difference to the total data is used as the standard to estimate the gross error. First, the time difference ∆TNMST is read, and then the Mahalanobis distance Di of each time difference ∆TNMST(i) to the overall data is computed. Then, the time difference ∆TNMST(i) whose MD is less than the threshold d is chosen to form a new time difference sequence N1.
If the average value D i ¯ is taken as the threshold d, the selected threshold value d results in a large amount of the gross time difference being ignored, and it cannot be discriminated. If the minimum value of the time difference min (Di) is selected, the effective time difference is discriminated. Because the time difference’s Mahalanobis distance approximately obeys the chi-squared distribution with the degree of freedom P [27], the chi-squared distribution value X P 2 ( ) can be selected as the threshold value d. When the value Di of the time difference ΔTNMST(i) is larger than X P 2 ( ) , the ΔTNMST(i) is considered to be the gross error. The discriminant expression is shown in Formula (10).
D i > X P 2 ( )
In the equation above, the parameter represents the confidence level and X P 2 ( ) represents the chi-squared distribution value. In this article, the confidence probability P = 97.5% is selected. There are two variables in the formula; in the corresponding table with α = 0.975, n = 2, the value X P 2 ( ) is 0.051. The process of time difference gross error discrimination is shown in Figure 3.

3.2. Time Difference Gross Error Processing Algorithm Combined with the Mahalanobis Distance and Group Covariance

As seen above, the MD can be used as a statistic for outliers and it can estimate the gross error in the time difference accurately. While the data contain useful information, they need to be processed. As described in Equation (8), the MD can be understood as the degree of difference between two random variables which obey the same distribution and have a covariance matrix of Σ (Σ represents the overall covariance matrix of the sample).
In order to better reflect the variation in the original time difference, this paper divides the time difference data into groups to obtain the covariance of each group. The MD for the time difference data combined with the covariance of the corresponding group as the weight cooperates with the mean value, with which gross time difference is replaced. The above processing can reduce random fluctuation and is suitable for the real-time dynamic processing of time difference data. Therefore, a time difference gross error processing method based on the combination of the MD and group covariance is presented in the paper. The distance between two samples (xixj) in Equation (9) is replaced by the distance from the ΔTNMST(i) to the mean value Δ T NMST ( i ) ¯ , which can be expressed as follows:
D i = ( Δ T NMST ( i ) Δ T ¯ ) T · C 1 · ( Δ T NMST ( i ) Δ T ¯ )
From the above equation, the distance between the time difference ∆TNMST(i) and the mean value Δ T NMST ( i ) ¯ is C · D i . In this paper, the time difference ∆TNMST(i) is divided into n groups, for which the covariance Cn of each group can be calculated. Therefore, the Di of ∆TNMST(i) is multiplied by the corresponding covariance Cn of the number n group of time difference data as the weight. The number i time difference ∆TNMST(i) processed by this algorithm is Fi, shown in Equation (12):
F i = Δ T NMST ( i ) ± ( C n · D i )   ( I = 1 ,   ,   n )
The flow chart of the time difference gross error processing method proposed in this paper is shown in Figure 4. First, the mean value of the time difference data is obtained, and the data are divided into n groups to calculate the covariance Cn of each group. Then, Di of ∆TNMST(i), for calculating time difference, is cyclically traversed to determine whether Di is less than the threshold d or not.
If Di < d, the time difference ∆TNMST(i) is normal, and this value is stored in the array N1 to continue traversing the time difference ∆TNMST(i+1); if Di ≥ d, the time difference ∆TNMST(i) is considered to be a gross error and needs to be processed. When the time difference ∆TNMST(i) is greater than the mean value of the data Δ T NMST ¯ , the ∆TNMST(i) is subtracted from the square root of the covariance Cn of the corresponding group of time difference data, multiplied by the time difference ∆TNMST(i). When the time difference ∆TNMST(i) is smaller than the mean value Δ T NMST ¯ of the data, ∆TNMST(i) is added to the square root of the covariance Cn of the corresponding group of time difference data, multiplied by the time difference ∆TNMST(i).
In this process, when each time difference gross error processing finishes, it is stored in the array N2 until all the gross error processing is completed. The time difference data Fi in N2 and the normal time difference array N1 form a new time difference array N.
The above processing can reduce random fluctuation and is suitable for the real-time dynamic processing of time difference data. Meanwhile, the algorithm is fast, simple, and easy to implement, which enables it to effectively meet the real-time dynamic processing requirements and improves the stability of the RTD-fluxgate system.

4. Simulation Data Analysis

In order to verify the effectiveness of the algorithm proposed in this paper, the Matlab Simulink random signal module is used to simulate the time difference in RTD-fluxgate; three groups of time difference data are obtained, with 1000 time differences in each group, and with different mean values and variance. The Laida criterion, moving average filtering, and the algorithm proposed in the paper are, respectively, used for processing. The statistical values resulting from the different processing methods are compared, as shown in Table 1.
As seen in Table 1, the mean values after processing using the three methods remain basically unchanged. The deviation and the fluctuation in the time difference are reduced by 2.034% and 1.797%, respectively, after processing using the Laida criterion. Through the use of moving average filtering, the variance of time difference is reduced by 98.428% and the fluctuation decreases by at least 88.963%. However, after using the algorithm proposed in this paper, the variance of the time difference is reduced by 99.392% and the fluctuation decreases by at least 90.896%. It can be concluded that the method of combining the MD and group covariance has a better processing effect.
The traditional signal-to-noise ratio reflects the ratio between the energy of the signal and the noise, but it can also evaluate the effect of noise processing. In order to better evaluate the processing effect of the time difference gross error, we consider the working principle of RTD-fluxgate and the particularity of the physical quantity of the magnetic core negative saturation time ΔTNMST. In this paper, the ratio of the core negative saturation time mean value Δ T NMST ¯ and the noise-induced uncertainty δTNMST) are defined as the time difference signal-to-noise ratio ΔTSNR, as is shown in Equation (13), and this parameter is used to quantitatively evaluate the processing effect of the time difference coarse error.
Δ T SNR = Δ T NMST ¯ δ ( Δ T NMST )
When the observation time is Tob and the number of cycles collected is N, then
T ob = N T e = N ( Δ T PMST + Δ T NMST )
In this formula, ΔTPMST and ΔTNMST represent the residence time of the magnetic core in the positive and negative saturation state, respectively. Since the residence time distribution of the positive and negative saturated states is the same and uncorrelated, and the standard deviation is approximately equal, then
σ Δ T NMST σ Δ T PMST
The uncertainty δTNMST) of the magnetic core negative saturation time ΔTNMST is shown in Equation (16):
δ ( Δ T NMST ) = σ Δ T NMST N
In this formula, σΔTNMST is the standard deviation of the ΔTNMST distribution. By substituting Equation (16) into Equation (13), the expression for the time difference signal-to-noise ratio ΔTSNR can be obtained as follows:
Δ T SNR = Δ T NMST ¯ δ ( Δ T NMST ) = Δ T NMST ¯ σ Δ T NMST N
The first set of simulated time difference data is selected, and the ΔTSNR processing using the three different methods is shown in Table 2. It can be seen from the table that the ΔTSNR obtained using the method based on the combination of the MD and group covariance is obviously higher than that obtained using the other two methods, so this method is suitable for processing the coarse value of the time difference.

5. Experiment and Preliminary Results

Due to the existence of geomagnetic fields and external interference from magnetic fields, the test was greatly influenced. In order to improve the stability and accuracy of the test, a Helmholtz coil was placed in the middle area of the electromagnetic shielding cylinder made of a five-layer permalloy in an electromagnetic shielding room of the National Geophysical Exploration Instrument Engineering and Technology Research Center of Jilin University. The RTD-fluxgate sensor was made by the Key laboratory of geophysical exploration equipment and it was placed in the uniform region of the Helmholtz coil; the measured magnetic field Hx, which was parallel to the core axis, was applied to the sensing unit.
The precision current sources of KEITHLEY 6221 were utilized in this experiment to drive the excitation coil of the RTD-fluxgate sensor to generate an excitation magnetic field. The excitation magnetic field waveform adopted a trapezoid-wave magnetic field with high time difference stability. When the frequency f = 30 Hz and the current I = 60 mA, the induction voltage produced by the induction coil pass through, in turn, the preamplifier circuit Pre-Amp, the secondary amplifier circuit Sec-Amp, and the shaping circuit. The shaping circuit consists of an addition circuit, a bias voltage generating circuit, and a hysteresis comparison circuit. When the induction voltage passes through the hysteresis comparison circuit, if using a threshold slightly lower than the peak value of the induction voltage, the shaping circuit obtains a rectangular wave with a duty cycle which varies with the information of the measured magnetic field Hx. This rectangular signal is input to the CH1 channel of the FPGA logic signal processor. The excitation current source is adjusted to generate a synchronous trigger pulse. When the excitation voltage amplitude is zero, this is set as a trigger point, and the synchronous trigger pulse is input to the CH2 channel of the FPGA. It uses two-channel signals to count the time difference, in which the counting frequency fc of FPGA is 100 MHz. The number of time points N is converted into the time difference ΔTNMST and transmitted to STM32 for storage. The relevant instrument connection and test environment are shown in Figure 5.
By changing the current at both ends of the Helmholtz coil, Hx is gradually moved through the range of −52,500~+52,500 nT and the least squares method is used to fit it linearly. When the sum of squared deviations from each data point to the fitting curve is minimized, the relationship between ΔTNMST and Hx is obtained as follows:
ΔTNMST = 0.008039 × Hx + 1.612 × 104
As in the equation above, the sensitivity of the sensor is SNMST = 0.008039 μs/nT. The relationship between the output ΔTNMST and the measured magnetic field Hx, as well as the linear deviations, are shown in Figure 6.
Figure 6a shows the fitting curve of ΔTNMST. According to the linear relationship, the sensitivity of the sensor SNMST is 0.008039 μs/nT. Figure 6b shows the linear deviations of ΔTNMST. When Hx = +52,500 nT, ΔTNMST is 16,517.17 μs and the maximum residual value between the measured results and the fitting curve is 8.87 μs. Therefore, the linear error (relative error) of the sensor is γ < ±0.054%. It can be seen from the figure that the RTD-fluxgate system possesses good linearity in the whole range of measurement.
To verify the effectiveness of the proposed algorithm in this paper, when Hx = +50,000 nT and the observation time t = 1 h, the ΔTNMST are divided into 30 groups and averaged. Then, the ΔTNMST are processed using the Laida criterion, moving average filtering, and the algorithm proposed in this paper, respectively. Because the observation time is longer, the data of the time difference ΔTNMST are larger. We only present the ΔTNMST processed using different methods for 60 s among 1 h, as shown in Figure 7.
As is illustrated in Figure 7, after being processed using the proposed algorithm, the stability of ΔTNMST is obviously higher than that of the other two methods. The statistical value of ΔTNMST processed using different methods is shown in Table 3. According to the table, with the mean of ΔTNMST unchanged, the time difference signal-to-noise ΔTSNR obtained using the method of combining the MD and group covariance is 47.776, which is about 34 times higher than that without any processing and is significantly higher than the other two methods. The standard deviation is reduced from 0.340 μs to 0.010 μs, which effectively decreases the time difference fluctuation caused by noise. The histogram of the ΔTNMST fluctuation distribution after being processed using the different methods is further calculated, as shown in Figure 8.
It can be seen from the above figure that the distribution of the ΔTNMST fluctuation is approximately a normal distribution, and the percentage of the ΔTNMST fluctuation greater than 1 μs is 0.360%. In Figure 8a, the fluctuation in ΔTNMST is 2.414 μs without processing. In Figure 8b, it is 1.995 μs after processing using the Laida criterion. In Figure 8c, it is 0.144 μs after processing using moving average filtering. In Figure 8d, it is 0.111 μs after processing using the combination of the MD and group covariance. This comparison indicates that the fluctuation in the ΔTNMST after using the proposed method is reduced by 95.402% compared with that without any processing, which significantly reduces the fluctuation in the ΔTNMST compared with other methods, so the stability of the RTD-fluxgate system can be improved.

6. Conclusions

On the basis of describing the structure and working principle of RTD-fluxgate, it is known that the time difference detection is affected by external interference and environmental factors, which causes coarse error in the time difference data and reduces the stability of the sensor. In this paper, a processing method is proposed to deal with the gross error in time difference by combining the MD with group covariance as weights cooperating with mean values. The results of our simulation and experiment indicate that the proposed method is more advantageous in identifying the time difference gross error: the signal-to-noise ratio of the time difference is improved by 34 times and the fluctuation in the ΔTNMST is reduced by 95.402%. The processing effect is better than that of other methods, and the influence of random interference on time difference detection is reduced. Meanwhile, the algorithm is fast, simple, and easy to implement, which can effectively meet the real-time dynamic processing requirements and improve the stability of RTD-fluxgate.

Author Contributions

Conceptualization, N.P.; methodology, N.P. and D.W.; software, Y.Y. and R.W.; investigation, N.P., D.W., Y.Y. and R.W.; writing—original draft preparation, N.P. and Y.Y.; writing—review and editing, N.P. and R.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 42004153), Natural Science Foundation of Jilin Province (Grant No. YDZJ202201ZYTS489), Science and Technology Research Project of Education Department of Jilin Province (Grant No. JJKH20220051KJ), and Science and Technology Bureau Outstanding Young Talent Project of Jilin City (Grant No. 20190104170).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank the Key Laboratory of Geo-Exploration Instrumentation (Jilin University) Ministry of Education for their additional assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, S.; Hu, X.Y.; Guo, N.; Cai, H.Z.; Zhang, H.L.; Li, Y.T. Overview of Aeromagnetic Measurement Technology of UAV. J. Wuhan Univ. (Inf. Sci. Ed.) 2023, 48, 823–840. [Google Scholar]
  2. Li, Z.P.; Gao, S.; Wang, X.B. Research on Aeromagnetic Measurement of Rotorcraft UAV in Special Area. J. Geophys. 2018, 61, 3825–3834. [Google Scholar]
  3. Li, J.J.; Zhang, X.; Shi, J.Q.; Heidari, H.D.; Wang, Y.Z. Performance Degradation Effect Countermeasures in Residence Times Difference (RTD) Fluxgate Magnetic Sensors. IEEE Sens. J. 2019, 19, 11819–11827. [Google Scholar] [CrossRef]
  4. Ando, B.; Baglio, S.; Bulsara, A.; Sacco, V. Residence times difference fluxgate. Measurement 2005, 38, 89–112. [Google Scholar] [CrossRef]
  5. Nikitin, A.; Stocks, N.G.; Bulsara, A.R. Bistable sensors based on broken symmetry phenomena: The residence time difference vs. the second harmonic method. Eur. Phys. J. Spec. Top. 2013, 222, 2583–2593. [Google Scholar] [CrossRef]
  6. Trigona, C.; Sinatra, V.; Ando, B.; Baglio, S.; Bulsara, A.R.; Mostile, G.; Nicoletti, A.; Zappia, M. Measurements of iron compound content in the brain using a flexible core fluxgate magnetometer at room temperature. IEEE Instrum. Meas. Mag. 2018, 67, 971–980. [Google Scholar] [CrossRef]
  7. Andò, B.; Baglio, S.; Sacco, V.; In, V. Noise effects in RTD-fluxgate. In Proceedings of the SENSORS, 2005 IEEE, Irvine, CA, USA, 30 October–3 November 2005; pp. 935–938. [Google Scholar]
  8. Carlo, T.; Valentina, S.; Bruno, A.; Salvatore, B.; Ratan, B.A. Flexible microwire residence times difference fluxgate magnetometer. IEEE Instrum. Meas. Mag. 2017, 66, 559–568. [Google Scholar]
  9. Pang, N.; Wang, Y.Z.; Cheng, D.F. Research on sensitive unit Barkhausen noise processing of residence times difference fluxgate. Chin. J. Sci. Instrum. 2015, 36, 2594–2601. [Google Scholar]
  10. Farrell, G.P.; Hill, E.W. The limit of fluxgate sensitivity due to Barkhausen noise for single layer and bi-layer permalloy thin film cores. IEEE Trans. Magn. 2002, 31, 4050–4052. [Google Scholar] [CrossRef]
  11. Chen, S.Y.; Wang, Y.Z.; Li, J.J.; Piao, H.; Lin, J. Sensitivity Model for Residence Times Difference Fluxgate Magnetometers Near Zero Magnetic Field. IEEE Sens. J. 2019, 20, 868–875. [Google Scholar] [CrossRef]
  12. Pang, N.; Cheng, D.F.; Wang, Y.Z.; Wu, S.J.; Liu, F.M.; Li, F. Study on Sensitivity and Stability of Time Difference Fluxgate. Instrum. Tech. Sens. 2016, 3, 17–20+27. [Google Scholar]
  13. Wang, X.K.; Zhu, X.W. Three Statistical Methods of Discriminating and Eliminating Outliers in Metrological Verification. China Test Conf. 2018, 44, 41–44. [Google Scholar]
  14. Zhang, B.M.; Jiang, Q.Y.; Mo, D.; Xiao, L.Y. A New Method to Deal with Coarse Errors in Electromagnetic Exploration Data. J. Geophys. 2015, 58, 2087–2102. [Google Scholar]
  15. Hao, X.E.; Shen, D.W.; Zhang, Y.B.; Huang, F.F. A Method of Fitting Calibration Data of Bore Pressure Tester Based on Grubbs Criterion. Sens. Microsyst. 2022, 41, 61–64+76. [Google Scholar]
  16. Wang, B.; Li, Z.Y.; Sun, Y.; Wang, K.X.; Pei, Y. Abnormal Power Data Cleaning and Prediction Application of PV Based on Least Square Filter Chauvelle Criterion. Kunming Univ. Sci. Technol. 2021, 46, 13. [Google Scholar]
  17. Tang, Y.F.; Zhang, H.; Xu, G.T.; Liu, J.J. Implementation of The Attitude Angle Measurement Technology Based on Weighted Moving Average. Chin. J. Sci. Instrum. 2012, 33, 1844–1849. [Google Scholar]
  18. Xu, B.; Wang, Q.D.; Li, S.X.; Wang, L.Z. Collaborative Localization Method Based on Mahalanobis Distance Combined with Adaptive Filtering. Chin. J. Inertia Technol. 2021, 29, 8. [Google Scholar]
  19. Wang, K.Y.; Wang, H.W.; Zhao, Q.; Wang, S.S. A Modified Mahalanobis Distance Discrimination Method. J. Beijing Univ. 2022, 48, 824–830. [Google Scholar]
  20. Andò, B.; Baglio, S.; Caruso, V.; Sacco, V.; Bulsara, A. Multilayer based technology to build RTD fluxgate magnetometer. IEEE Instrum. Meas. Mag. 2006, 65, 509–514. [Google Scholar]
  21. Wang, Y.; Wu, S.; Lu, H. Output performance optimization for RTD fluxgate sensor based on dynamic permeability. Sci. China Inf. Sci. 2016, 59, 112213. [Google Scholar] [CrossRef]
  22. Pang, N.; Cheng, D.F.; Wang, Y.Z. A High Stability Time Difference Readout Technique of RTD-fluxgate. Sensors 2017, 17, 2325. [Google Scholar] [CrossRef] [PubMed]
  23. Andò, B.; Ascia, A.; Baglio, B.; Bulsara, A.R.; Neff, J.D.; In, V. Towards an optimal readout of a residence times difference (RTD) fluxgate magnetometer. Sens. Actuators A 2008, 142, 73–79. [Google Scholar] [CrossRef]
  24. Chen, S.Y.; Wang, Y.Z.; Lin, J. A SFTD algorithm for optimizing the performance of the readout strategy of residence time difference fluxgate. Sensors 2018, 18, 3985. [Google Scholar] [CrossRef] [PubMed]
  25. Zu, Z.W.; Li, Q. Mahalanobis Distance Fuzzy Clustering Algorithm Based on Particle Swarm Optimization. J. Chongqing Univ. Posts Telecommun. 2019, 2, 279–284. [Google Scholar]
  26. Daniel, I.; Eduardo, G.; Jesús, S. Incipient Wear Detection of Welding Gun Secondary Circuit by Virtual Resistance Sensor Using Mahalanobis Distance. Sensors 2023, 23, 894. [Google Scholar]
  27. Nguyen, V.D.; Trung, N.L.; Abed, M.K. Robust Subspace Tracking Algorithms using Fast Adaptive Mahalanobis Distance. Signal Process. 2022, 195, 108402. [Google Scholar] [CrossRef]
  28. Mei, J.Y.; Liu, M.Z.; Wang, Y.F.; Gao, H.J. Learning a Mahalanobis Distance Based Dynamic Time Warping Measure for Multivariate Time Series Classification. IEEE Trans. Cybern. 2016, 46, 1363–1374. [Google Scholar] [CrossRef]
Figure 1. Schematic diagrams of the RTD-fluxgate structure and its working principle. (a) The structure of an RTD-fluxgate sensor; (b) A schematic diagram of the RTD-fluxgate system’s working principle.
Figure 1. Schematic diagrams of the RTD-fluxgate structure and its working principle. (a) The structure of an RTD-fluxgate sensor; (b) A schematic diagram of the RTD-fluxgate system’s working principle.
Sensors 23 09223 g001
Figure 2. Diagram of the NMST readout strategy.
Figure 2. Diagram of the NMST readout strategy.
Sensors 23 09223 g002
Figure 3. The flow chart for time difference gross error discrimination based on the MD.
Figure 3. The flow chart for time difference gross error discrimination based on the MD.
Sensors 23 09223 g003
Figure 4. The flow chart of the processing algorithm.
Figure 4. The flow chart of the processing algorithm.
Sensors 23 09223 g004
Figure 5. Instrument connection and test environment.
Figure 5. Instrument connection and test environment.
Sensors 23 09223 g005
Figure 6. Fitting linearity deviations of ΔTNMST with I= 60 mA and f = 30 Hz. (a) Fitting curve of ΔTNMST; (b) deviations of ΔTNMST.
Figure 6. Fitting linearity deviations of ΔTNMST with I= 60 mA and f = 30 Hz. (a) Fitting curve of ΔTNMST; (b) deviations of ΔTNMST.
Sensors 23 09223 g006
Figure 7. Comparison diagram of ΔTNMST processed using different methods.
Figure 7. Comparison diagram of ΔTNMST processed using different methods.
Sensors 23 09223 g007
Figure 8. Histograms of ΔTNMST fluctuation distribution using different methods. (a) Fluctuation in ΔTNMST without processing; (b) fluctuation in ΔTNMST after processing using the Laida criterion; (c) fluctuation in ΔTNMST after processing using moving average filtering; (d) fluctuation in ΔTNMST after processing using the combination of the MD and covariance.
Figure 8. Histograms of ΔTNMST fluctuation distribution using different methods. (a) Fluctuation in ΔTNMST without processing; (b) fluctuation in ΔTNMST after processing using the Laida criterion; (c) fluctuation in ΔTNMST after processing using moving average filtering; (d) fluctuation in ΔTNMST after processing using the combination of the MD and covariance.
Sensors 23 09223 g008
Table 1. The statistical values from the different processing algorithms.
Table 1. The statistical values from the different processing algorithms.
First Set of DataSecond Set of DataThird Set of Data
Mean ValueVariance
(≈)
FluctuationMean ValueVariance
(≈)
FluctuationMean ValueVariance
(≈)
Fluctuation
original data804.54530.91132.280500.00010.17819.881200.00020.57328.042
Laida criterion804.52429.93831.700499.8579.97118.023200.08319.76525.435
moving average804.5690.3243.299499.9050.1602.194200.0350.2222.591
algorithm in paper804.5460.1262.543499.9100.0371.771200.0220.1252.553
Table 2. The ΔTSNR of the simulated time difference using different processing methods.
Table 2. The ΔTSNR of the simulated time difference using different processing methods.
Original DataLaida CriterionMoving Average Algorithm in Paper
Δ T NMST ¯ (μs)804.545804.524804.569804.546
σΔTNMST (μs)5.6825.4720.5690.355
ΔTSNR4.4774.64944.71471.667
Table 3. Statistical values of ΔTNMST using different processing methods.
Table 3. Statistical values of ΔTNMST using different processing methods.
Mean Value (μs)Fluctuation (μs)Standard Deviation (μs)ΔTSNR
original data16,516.1692.4140.3401.349
Laida criterion16,516.1681.9950.3351.370
moving average filtering16,516.1900.1440.03214.246
algorithm in this paper16,516.1680.1110.01047.776
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pang, N.; Wang, D.; Yang, Y.; Wang, R. Research on a Time Difference Processing Method for RTD-Fluxgate Data Based on the Combination of the Mahalanobis Distance and Group Covariance. Sensors 2023, 23, 9223. https://doi.org/10.3390/s23229223

AMA Style

Pang N, Wang D, Yang Y, Wang R. Research on a Time Difference Processing Method for RTD-Fluxgate Data Based on the Combination of the Mahalanobis Distance and Group Covariance. Sensors. 2023; 23(22):9223. https://doi.org/10.3390/s23229223

Chicago/Turabian Style

Pang, Na, Dan Wang, Yuhan Yang, and Rui Wang. 2023. "Research on a Time Difference Processing Method for RTD-Fluxgate Data Based on the Combination of the Mahalanobis Distance and Group Covariance" Sensors 23, no. 22: 9223. https://doi.org/10.3390/s23229223

APA Style

Pang, N., Wang, D., Yang, Y., & Wang, R. (2023). Research on a Time Difference Processing Method for RTD-Fluxgate Data Based on the Combination of the Mahalanobis Distance and Group Covariance. Sensors, 23(22), 9223. https://doi.org/10.3390/s23229223

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop