Adaptive Basis Function Method for the Detection of an Undersurface Magnetic Anomaly Target

: The orthogonal basis functions (OBFs) method is a prevailing choice for the detection of undersurface magnetic anomaly targets. However, it requires the detecting platform or target to move uniformly along a straight path. To circumvent the restrictions, a new adaptive basis functions (ABFs) approach is proposed in this article. It permits the detection platform to search for a possible target at different speeds along any course. The ABFs are constructed using the real-time data of the onboard triaxial ﬂuxgate, GPS module, and attitude gyro. Based on the pseudo-energy of an apparent target signal, the constant false alarm rate (CFAR) method is employed to judge whether a target is present. Moreover, by deﬁning the pixel as a relative possibility for a target at a geographic location, a magnetic anomaly target imaging scheme is introduced by displaying the pixels onto the searching area. On-site experimental data are utilized to demonstrate the proposed approach. Compared with the traditional OBFs method, the present ABFs approach can substantially improve the detection possibility and reduce false alarms.

A ferromagnetic object may be treated as a magnetic dipole if the separation from the object to the magnetic sensor is more than several times the target's dimension.Detection of a magnetic dipole target may involve two aspects: judgment of whether a magnetic dipole target is present and determination of the target location.
To judge the existence of a target from a piece of measured data, one may choose a signal-based or noise-based scheme.For a signal-based method, the pseudo signal energy is estimated first, and then the constant false alarm rate (CFAR) method is employed by dynamically comparing the signal energy with the referenced background noise.The orthogonal basis function (OBF) decomposition first introduced in [16] may be the most popular approach to estimate the signal energy, by which the scalar magnetic anomaly signal of a magnetic dipole target is expressed by three OBFs, and the signal energy is the sum of the squares of the three expansion coefficients.If a scalar gradient signal is available by using multiple scalar magnetometers, four OBFs can be used to express the gradient signal, and the pseudo signal energy is estimated using the four expansion coefficients [17].The expansion coefficients may be projected on an ellipsoid surface to reduce the impact of noise, and the modified coefficients are used to estimate the energy [18].Because the OBFs method hinges on matched filtering, it is very suitable to detect a known signal contaminated by Gaussian white noise.To enhance the detection capability under complex colored noise, a lot of techniques were proposed, including whitening filtering [19], optimized basis functions [20], adaptive filtering [21], Karhunen-Loeve expansion [22], etc.However, these signal-based methods assume that the detecting platform carrying the magnetic sensor moves uniformly along a straight path, so that the OBF expansion applies.If the detecting platform moves non-uniformly or along a hovering forward course, there is no target signal model available at the present stage.
Noise-based detection methods are more suitable when the target signal model cannot be established.By inspecting the statistical change of noise, whether a target signal is intruding may be judged.These kinds of schemes include the minimum entropy filtering (MEF) method [23], high-order crossing (HOC) method [24], fractal feature method [25], etc. Information entropy is calculated by estimating the probability density function of noise.The MEF approach leverages information entropy to detect a target signal, as it would experience a substantial decrease when a target signal is overlaid.The HOC method assesses the presence of a target signal by scrutinizing the disparity in HOC information between the measured data and background noise.Similarly, the difference in fractal characteristics between the target signal and noise is exploited to judge the presence of a target signal.However, these noise-based approaches are somewhat less robust and mainly operate by examining the differences between measured data and noise, rather than capitalizing on the features of target signals.In addition, some scholars have introduced detection methods based on nonlinear stochastic resonance systems to enhance the signal using noise energy [26,27].But the stochastic resonance method demands stringent system parameters and exhibits limited stability.Furthermore, a coherent filtering method may be employed to suppress the background magnetic noise if measured data from a reference magnetometer is available [28].Compared with signal-based methods, noise-based detection methods may have higher detection rates but higher false alarm rates, as well.
In recent years, machine learning methods have witnessed a remarkable upsurge in many research fields.In [29], the support vector machine (SVM) was applied to magnetic anomaly detection, which transformed the determination of signal presence into a binary classification conundrum.However, the SVM technique is susceptible to the pitfall of underfitting, which potentially diminishes its practical efficacy.Convolutional neural networks (CNNs) were attempted in [30][31][32], where both 1-D CNNs in tandem with signal features and 2-D CNNs based on time-frequency diagrams were designed and trained to recognize magnetic anomaly signals.A deep learning model via an attention mechanism for vector magnetic anomaly detection was introduced in [33].In [34], an end-to-end deep learning (DL) framework was introduced for magnetic anomaly detection (MAD) and denoising.In [35], a network with adaptive time-frequency feature expression to detect magnetic anomalies was proposed.Nevertheless, it should be noticed that machine learning (ML) methods exhibit a pronounced reliance on data and some degree of instability in terms of generalization capabilities.
As for the determination of the location and magnetic moment of a magnetic dipole target, there are six scalar unknowns to solve.Based on the single magnetic dipole model, a magnetic dipole tracking and locating approach utilizing magnetic gradient tensor data for direct inversion was initially proposed in [36].Both direct and iterative approaches have been reported in the literature.Direct methods are usually based on the availability of measured data from multiple magnetometers or magnetic gradiometers [37][38][39].The Euler deconvolution method [40] and phase-based filtering method [41] may also be categorized as direct schemes.It should be mentioned that direct methods are usually sensitive to measurement errors if high-order derivatives are required and calculated using the measured data, which results in unstable solutions.In contrast, iterative methods dispense with the calculations of derivatives and establish profound mathematical map-ping connecting the target location and magnetic moment with the magnetic anomaly fields.Directionally sensitive algorithms are usually employed in iterative methods, such as the Levenberg-Marquardt algorithm [42], Gauss-Newton algorithm [43], simulated annealing [44], and genetic algorithm [45].ML-based methods were also developed for the localization of undersurface magnetic targets [46,47].Nevertheless, these schemes are subject to the availability of sufficient training datasets.
The purpose of this article is to present an alternative approach for detection and localization of undersurface magnetic dipole targets, such as shipwrecks, using an aeromagnetic system.To meet the need to search for a target along an arbitrary course, a set of adaptive basis functions (ABFs) are derived to replace the orthonormal basis functions (OBFs) that are suitable only for straight courses.The target location and magnetic moment are determined by directly scanning the searching area and optimally matching the theoretical signal to the measured values.To show the possibility for a target localization in global geographic longitude and latitude coordinates, an imaging scheme is introduced via transformation of the geocentric coordinate system.
The remainder of this article is organized as follows: In Section 2, the signal model is established by using the ABFs.Data preprocessing, target detection, and location methods are described in Section 3. In Section 4, on-site experiment validation is presented.Some concluding remarks are given in Section 5.

Target Signal Modeling Using ABFs
We assume that there are four sets of sensors installed on an aeromagnetic detecting platform, including an optical-pumping magnetometer (OPM), a triaxial fluxgate, a GPS module, and an attitude gyro.The OPM measures the magnitude of a vector magnetic field, which has very high precision, as small as a few pT (pico-Tesla).The triaxial fluxgate measures the three components of a vector magnetic field, which are used to calculate the attitude of the platform with respect to the local geomagnetic field.The GPS module records the platform positions or flying course.The attitude gyro records the attitudes of the detecting platform with respect to the geographic coordinate system.
The magnetic anomaly target can be taken as a magnetic dipole if the minimum distance from the target to the flying path is greater than several times the target's dimension, which is well satisfied in this study.The magnetic anomaly target signal measured by the OPM can be written as where T(t) is the directional vector of the local geomagnetic field, M is the target magnetic moment, I is the identical matrix, R(t) = r(t) − r s is the radius vector from the target to the OPM, and R(t) = |r(t) − r s |.The coordinates r(t) are provided by the GPS module, M and r s are the unknowns to be solved by using a piece of data centered at r(t d ), as shown in Figure 1.The T is obtained by using the data of the fluxgate and gyro as where u(t) is the local geomagnetic directional vector with respect to the platform measured by the fluxgate, and M(t) is a transform matrix from the platform system to the geographic system.Specifically, in which ϕ(t), ϑ(t) and ψ(t) are the yaw, pitch, and roll angles of the platform with respect to the geographic coordinate system.The basis functions in (1) are where , and so forth.We call them adaptive basis functions (ABFs) because they are constructed using the real-time data of r(t) and T(t).
and so forth.We them adaptive basis functions (ABFs) because they are constructed using the real-ti data of ( ) t r and ˆ( ) t T .
are called the closest proximity approach (CPA) and charac istic time, respectively.Substituting ( 4)-( 6) into (3), we will obtain ) where , , ) are coefficients relevant to the depression angle s  but in pendent of the CPA 0 R , and their detailed expressions are omitted.It shows that i F linear combination of the Anderson-type function ).These th Anderson-type functions can be orthonormalized to become the three orthogonal ba functions (OBF) given by Ginzburg in [17], so that ( ) i F t , and then the target signal sig B , can be expressed by the three OBFs.
In this paper, the general signal model (1) will be adopted, which permits the dete ing platform to search for a possible target along any course and at different speeds If the platform moves uniformly at speed v along a straight line, refer to Figure 1; we will have where R 0 and τ = R 0 /v are called the closest proximity approach (CPA) and characteristic time, respectively.Substituting ( 4)-( 6) into (3), we will obtain where a i , b i , c i (i = 1, 2, 3) are coefficients relevant to the depression angle ϑ s but independent of the CPA R 0 , and their detailed expressions are omitted.It shows that F i is a linear combination of the Anderson-type function w n /(1 + w 2 ) 2.5 (n = 1, 2, 3).These three Anderson-type functions can be orthonormalized to become the three orthogonal basis functions (OBF) given by Ginzburg in [17], so that F i (t), and then the target signal B sig (t), can be expressed by the three OBFs.In this paper, the general signal model (1) will be adopted, which permits the detecting platform to search for a possible target along any course and at different speeds by using a GPS module and a flight gyro to record the path and attitude of the detecting platform, respectively.

Data Preprocessing
In general, the measured values by the OPM may be written as five parts: where B geo (t) is the local geomagnetic field, B pla (t) denotes the interference magnetic field generated by the detecting platform, B diu (t) reflects the diurnal magnetic variation, and B noi (t) accounts for other magnetic noises.
We fit the geomagnetic field using linear interpolation within the interval t 0 ≤ t ≤ t s : where α i (0 ≤ i ≤ 10) are undetermined coefficients, and In the above, x(t), y(t), and z(t) are the geographic coordinates of the detecting platform obtainable from the onboard GPS module.Please note that a constant term and a linear term are included in (9) because they are definitely not target signals and should be removed anyway.
The magnetic interference generated by the detecting platform is compensated by using a classical Tolles-Lawson (TL) model [48]: where the first term denotes the permanent interference due to the remanent magnetism of the aircraft material, the second term denotes the induction interference due to the magnetization of the aircraft body by the geomagnetic field, and the third term denotes the eddy interference due to the time-varying of the direction of the local geomagnetic field with respect to the platform system.In (11) h i (t) = dh i (t)/dt.The geomagnetic field and platform interference are the principal quantities that swamp the target signal and should be eliminated first.To this end, we use the least square method (LSM) with ridge regression to solve the coefficients α in (9) and a, b, c in (11), i.e., defining and minimizing the following function: In this equation, ( 9) and ( 11) are substituted into (12), making ∂ f /∂α = ∂ f /∂β = 0 to solve for α and β, where α and β contain the coefficients {α i } in (9) and a i b ij c ij in (11), respectively.
Once these coefficients are found, we will obtain the remaining magnetic field as follows: Because the presence of a target signal is a rare event, B rem (t) is largely the diurnal magnetic variation, which may be seen as a Gaussian-colored noise.So, a whitening filter is constructed as where w(t) is white noise.The filter coefficients are found by the Yule-Walker equation: where * represents the statistical average over time t.In actual operation, the B diu (t) in ( 16) and ( 17) is replaced with the B rem (t) of ( 13).The well-known Levinson-Durbin algorithm is used to solve (15).Finally, we obtain the sufficiently preprocessed data as which is regarded as a possible target signal superimposed on whiting background noise.

Target Signal Detection
First, we separate out the apparent target signal from (18) using a sliding window centered at the detecting instant t d , i.e., |t − t d | ≤ ∆T/2 by minimizing the functional where B sig (t) is the signal model of (1) based on the ABFs given in (3).As illustrated in Figure 1, the target is supposed to be located at The magnetic moment M is found by solving ∂ f (t d )/∂M i = 0 (i = 1, 2, 3).By scanning the horizontal parameter Y at a step ∆Y, i.e., Y = k∆Y, we obtain a group of apparent signals Next, we calculate the pseudo energy of each apparent signal and the averaged energy as where a total of K + 1 horizontal locations are checked.
Then, we adopt the well-known constant false alarm rate (CFAR) method to judge whether there is a target signal in the segment of data with the threshold value: α = − ln P CFAR (25) where P CFAR is the prescribed CFAR, say, if P CFAR = 0.1, α = 2.3.In (24), the length of protected elements centered at the detecting instant t d is taken to be ∆T, and the lengths of reference elements on each side are also taken to be ∆T.

Target Locating and Imaging
To display the relative possibility for a target residing at each location, a 2D image may be defined as Obviously, the larger the I(t d , k), the greater the possibility for a target to be present at the location given by (20).
Furthermore, we can convert the continuous strip image of ( 26) into a regional image on the global geographic coordinate system as t d steps forward.As shown in Figure 2, if the longitude, latitude, and height given by the onboard GPS are ϕ(t), ψ(t) and h(t), the position of the detecting platform in the geocentric coordinate system is [49] 1 − e 2 sin 2 θ where a = 6378, 137.0 m, b = 6356, 752.314245 m, and e = 1 − b 2 /a 2 ; the relationship between the geocentric latitude θ and geographic latitude ψ is Then, the local geographic coordinates x(t), y(t) and z(t) (x-axis to the north, y-axis to the east, and z-axis downward) within the window |t − In which In which ( ) According to ( 27)-( 30), the geographic longitude and latitude positions   As a result, the strip image of ( 26) is converted into a regional image labelled with the geographic longitude and latitude as Subsequent detection using image detection techniques may be attempted.

Experiment Validation
On-site experiments were conducted to validate the applicability of the proposed procedure.Figure 3 shows the magnetic detection airplane used for our experiments.As shown in Figure 4, a total of 10 artificial targets are assumed to be distributed in an area.The detection airplane is searching for the targets along a hovering forward course.The 10 targets are successfully found using the proposed procedure, as indicated in the figure.A 2D image of the region around the 9th target is placed on the side of the figure.Consequently, if the target is present at the platform coordinates {x s , y s , z s }, its geo- centric coordinates are According to ( 27)-( 30), the geographic longitude and latitude positions {ϕ s , ψ s } of the target can be obtained by 1 − e 2 (Z S /a) 2 As a result, the strip image of ( 26) is converted into a regional image labelled with the geographic longitude and latitude as Subsequent detection using image detection techniques may be attempted.

Experiment Validation
On-site experiments were conducted to validate the applicability of the proposed procedure.Figure 3 shows the magnetic detection airplane used for our experiments.As shown in Figure 4, a total of 10 artificial targets are assumed to be distributed in an area.The detection airplane is searching for the targets along a hovering forward course.The 10 targets are successfully found using the proposed procedure, as indicated in the figure.A 2D image of the region around the 9th target is placed on the side of the figure.
The raw data collected by the optical-pumping magnetometer (OPM) is shown in Figure 5a.The horizontal axis labels the sampling points with a sampling rate of 10 points per second, or 10 Hz.There are many peaks and troughs in the measured curve, which reflect the geological magnetic gradient along the course and magnetic interferences due to the jolting of the airplane.The 10 red segments in the curve indicate the places of the 10 target signals.Figure 5b shows the preprocessed data (cyan line) using the schemes described in Section 3.1, including geomagnetic fitting, interference compensating, and whiting filtering and, finally, a bandpass filtering of 0.04-0.6Hz.The black line shows the direct bandpass filtering of the raw data.The 10 true target signals are shown in Figure 5b, too.It is obvious that the SNR of preprocessed data is greater than the SNR of directly filtered data.The proposed preprocessed schemes have little impact on the target signals but reduce the magnetic noise/interference level.The raw data collected by the optical-pumping magnetometer (OPM) is shown in Figure 5a.The horizontal axis labels the sampling points with a sampling rate of 10 points per second, or 10 Hz.There are many peaks and troughs in the measured curve, which reflect the geological magnetic gradient along the course and magnetic interferences due to the jolting of the airplane.The 10 red segments in the curve indicate the places of the 10 target signals.Figure 5b shows the preprocessed data (cyan line) using the schemes described in Section 3.1, including geomagnetic fitting, interference compensating, and whiting filtering and, finally, a bandpass filtering of 0.04-0.6Hz.The black line shows the direct bandpass filtering of the raw data.The 10 true target signals are shown in Figure 5b, too.It is obvious that the SNR of preprocessed data is greater than the SNR of directly  The raw data collected by the optical-pumping magnetometer (OPM) is shown in Figure 5a.The horizontal axis labels the sampling points with a sampling rate of 10 points per second, or 10 Hz.There are many peaks and troughs in the measured curve, which reflect the geological magnetic gradient along the course and magnetic interferences due to the jolting of the airplane.The 10 red segments in the curve indicate the places of the 10 target signals.Figure 5b shows the preprocessed data (cyan line) using the schemes described in Section 3.1, including geomagnetic fitting, interference compensating, and whiting filtering and, finally, a bandpass filtering of 0.04-0.6Hz.The black line shows the direct bandpass filtering of the raw data.The 10 true target signals are shown in Figure 5b,  To further examine the detection performance, a series of Monte Carlo simulation experiments were conducted for different signal-to-noise rates (SNRs).In this paper, the SNR is defined as  Figure 5c shows the detection results using the proposed ABFs method described in Section 3.2.It can be seen that the 10 targets were successfully detected without missing or false alarm.As for comparison, the prevailing OBFs method is applied, and the detection results are shown in Figure 5d, from which it can be seen that six targets are successfully detected (green places), four real targets are missing (gray places), and three false targets are distinguished (lilac places).The error detections by the OBFs method may be attributed to the signal model that requires the search path within a detection window to be along a straight line, but, in fact, it is along a curved line.
In Table 1, the positioning errors are calculated, which are the distances between the "reported" true target locations and the detected locations given by the relative longitudes and latitudes.The average error is 156.55 m.The error may come from two aspects.One is the algorithm, and the other is the data or systematic errors.In terms of algorithms, we set the resolutions to be 70 m in the platform-moving direction and 100 m in the lateral direction, respectively.So, the algorithm error may be estimated to be √ 35 2 + 50 2 ≈ 61 m.As for the data, the on-board GPS has an error of about 30 m, and both the direction of the local geomagnetic field measured by the triaxial fluxgate and the attitude of the platform measured by the flight gyro have some degree of uncertainty.Therefore, the averaged error of 156.55 m is a little big but acceptable.Please note that the flying heights are about 400 m, and the flight speeds are about 70 m/s.The lateral scanning range is within 1300 m on both sides of the flight route.For reasons of confidentiality, the latitude and longitude have been processed.
Note that the magnetic moment M solved by ( 19) is just for extracting the target signal as (21).The value of M is discarded once the signal is constructed, because the magnetic moment is taken to be a random value in the generation of simulated target signals, and the simulated signals are scaled according to specified SNRs, so that the original magnetic moment has been lost/forgotten.Therefore, our detection task is to determine whether and where a magnetic dipole target is present; the size of the magnetic moment (amplitude and orientation) is not concerned in this study.
To further examine the detection performance, a series of Monte Carlo simulation experiments were conducted for different signal-to-noise rates (SNRs).In this paper, the SNR is defined as where The background or noise magnetic fields are real measured data.Every time, starting from a random place, a piece of data of 5000 points (or 500 s, because the sampling rate is 10 points/s) is taken from the measured sequence (the total length of a measured sequence in one day is about 180,000 points or about 5 h, and 6 days' data are used in the simulation experiments), and a total of 1000 sets of background magnetic field data (each is 5000 points in length) are obtained.In the meantime, 1000 sets of data for fluxgate, flight gyro, and GPS, respectively, are obtained.
The simulated signals are generated by using (1).The magnetic moment is taken to be M = M x x + M y ŷ + M z ẑ, where M x , M y and M z are empirically random values meeting −7.4 × 10 7 < M x < 6.3 × 10 7 , −3.0 × 10 7 < M y < 2.3 × 10 7 , and −7.6 × 10 7 < M z < 1.0 × 10 8 .The ABFs are constructed by using (3), in which T(t) = M(t)•u(t) and R(t) = r(t) − r s .M(t) is obtained by using the flight gyro data; u(t) is obtained by using the fluxgate data; r(t) is obtained by using the GPS data; and r s is the scanned target position of (20) with Y = k∆Y, k = −13, • • • , 13, and ∆Y = 100 m.A total of 1000 simulated signals are and then added to the 1000 sets of background noises according to the specified SNRs given by (36).Please note that the length of a simulated signal is 1000 points, while the length of a set of background noise is 5000 points.We randomly take a segment of 1000 points from the 5000 points (say from 233 to 1232) to calculate the ABFs, and the simulated signal is also added to this segment.
Finally, for each specified SNR, we obtain 1000 sets of simulated experimental data or "realizations".Then, the ABFs, OBFs, and MED methods are employed to carry out the detections.In the detection process, the ABFs are dynamically constructed using the real-time data of fluxgate, flight gyro, and GPS.As for the OBFs, the flight routes for the 1000 realizations are assumed to be straight lines, and only one parameter, i.e., the characteristic time τ = R 0 /v, is adjustable with 5 ≤ τ ≤ 16.In the simulated detections, only whether a target is present is of concern.
In Figure 6, comparisons of detection probability between the present ABFs method with the conventional OBFs method and the minimum entropy detection (MED) method [23] are given.Figure 6a,b are the results without and with the data preprocessing described in Section 3.1.The preprocessing is valuable for, say, raising the detection rate from 60% to 90% for 0 dB SNR.It is obvious that the present ABFs approach outperforms the typical signal-based OBFs method and noise-based MED method, especially for SNRs in the range of −5 dB to 0 dB.

Concluding Remarks
In this article, a new adaptive basis functions (ABFs) method is proposed to remedy the conventional orthogonal basis functions (OBFs) method.The OBFs method demands the detecting platform move uniformly along a straight course.For an arbitrary flying course and a changing flight speed, so far, there is no target signal model available to match the measured data.The present signal model based on ABFs fills the gap.Appropriate preprocessing schemes are suggested to raise the signal-to-noise ratio (SNR) and detection possibility, which include the fitting of geological magnetic gradients, compen- Please note that the present ABFs method is mainly for post-processing, because we need to collect some data to complete the data preprocessing as described in Section 3.1, which typically takes 500 s.Except for this delay, the detection can be carried out in real

Figure 1 .
Figure 1.Illustration of airborne magnetic anomaly target detection.If the platform moves uniformly at speed v along a straight line, refer to Figur we will have ˆˆ( ) [( )cos sin ] [( )sin cos ]

Figure 1 .
Figure 1.Illustration of airborne magnetic anomaly target detection.

Figure 2 .
Figure 2. Geocentric coordinate system vs. geographic coordinate system.Consequently, if the target is present at the platform coordinates   , , s s sx y z , its geo-

Figure 4 .
Figure 4. On-site experiment flight trajectory and positioning results.

Figure 4 .
Figure 4. On-site experiment flight trajectory and positioning results.

Figure 4 .
Figure 4. On-site experiment flight trajectory and positioning results.

Figure 5 .
Figure 5. Raw data, preprocessed data, and detection results.(a) Raw data.The 10 red segments in the curve indicate the locations of 10 target signals.(b) Preprocessed data.The black line shows the direct bandpass filtering of the raw data; the cyan line shows the preprocessed data using the schemes described in Section 3.1; the red line shows the direct bandpass filtering of the target signal.(c) Detection results using the present ABFs method.The black line represents the energy curve; the red line indicates the threshold.When the energy value exceeds the threshold, the existence of a target signal at that point is indicated (green stripes).(d) Detection results using the OBFs method.The gray stripes represent missed target signals; lilac stripes represent falsely detected target signals; and green stripes indicate correctly detected target signals.

Figure 5 .
Figure 5. Raw data, preprocessed data, and detection results.(a) Raw data.The 10 red segments in the curve indicate the locations of 10 target signals.(b) Preprocessed data.The black line shows the direct bandpass filtering of the raw data; the cyan line shows the preprocessed data using the schemes described in Section 3.1; the red line shows the direct bandpass filtering of the target signal.(c) Detection results using the present ABFs method.The black line represents the energy curve; the red line indicates the threshold.When the energy value exceeds the threshold, the existence of a target signal at that point is indicated (green stripes).(d) Detection results using the OBFs method.The gray stripes represent missed target signals; lilac stripes represent falsely detected target signals; and green stripes indicate correctly detected target signals.

Figure 6 .
Figure 6.Comparisons of the detection probability against SNR of the present ABFs method with the prevailing OBFs method and MED method.(a) Without data preprocessing.(b) With data preprocessing.

Figure 6 .
Figure 6.Comparisons of the detection probability against SNR of the present ABFs method with the prevailing OBFs method and MED method.(a) Without data preprocessing.(b) With data preprocessing.
Finally, two curves E ave (t d ) ∼ t d and E threshold (t d ) ∼ t d are plotted, and a target is said to be present at the instant t d if E ave (t d ) > E threshold (t d ).The lateral position of the target is taken to be at Y tar = k tar ∆Y if E(t d , k) achieves the maximum value when k = k tar .In the following experiments, we set ∆T = 1000 points or 100 s, t d = n∆t with ∆t = 10 points or 1 s, and ∆Y = 100 m.

Table 1 .
Detection and positioning statistics table.
[B mea (t d )] denotes a piece of measured background data of length ∆T centered at the instant t d , i.e., |t − t d | ≤ ∆T/2.A simple process of detrending and subtracting the mean value is carried out.[S(td )] denotes a simulated target signal that is added on [B mea (t d )] to form contaminated data."PPV"stands for the peak-to-peak value of the signal or background, and "bpf" stands for a bandpass filtering of 0.02-0.4Hz.The generations of background magnetic fields and simulated signals are described in detail, as follows: