Robust Pedestrian Dead Reckoning Integrating Magnetic Field Signals and Digital Terrestrial Multimedia Broadcasting Signals

Currently, many positioning technologies complementary to Global Navigation Satellite System (GNSS) are providing ubiquitous positioning services, especially the coupling positioning of Pedestrian Dead Reckoning (PDR) and other signals. Magnetic field signals are stable and ubiquitous, while Digital Terrestrial Multimedia Broadcasting (DTMB) signals have strong penetration and stable transmission over a large range. To improve the positioning performance of PDR, this paper proposes a robust PDR integrating magnetic field signals and DTMB signals. In our study, the Spiking Neural Network (SNN) is first used to learn the magnetic field signals of the environment, and then the learning model is used to detect the magnetic field landmarks. At the same time, the DTMB signals are collected by the self-developed signal receiver, and then the carrier phase ranging of the DTMB signals is realized. Finally, robust pedestrian positioning is achieved by integrating position information from magnetic field landmarks and ranging information from DTMB signals through Extended Kalman Filter (EKF). We have conducted indoor and outdoor field tests to verify the proposed method, and the outdoor field test results showed that the positioning error cumulative distribution of the proposed method reaches 2.84 m at a 68% probability level, while that of the PDR only reaches 8.77 m. The proposed method has been validated to be effective and has good positioning performance, providing an alternative solution for seamless indoor and outdoor positioning.


Introduction
The signal of the Global Navigation Satellite System (GNSS) will be seriously attenuated and blocked in some environments (such as urban canyons, indoors, and underground) [1].At present, many positioning methods complementary to GNSS have been developed to provide ubiquitous location-based services (LBS).Compared to positioning methods that require the deployment of positioning facilities in advance, such as WiFi [2], Ultra-Wide Band (UWB) [3], Long Term Evolution (LTE) [4], 5G [5], audio [6], etc. Pedestrian Dead Reckoning (PDR) is an independent positioning method.However, due to its sensor characteristics and algorithm design, PDR faces serious challenges of position drift [7].
Magnetic field signals are ubiquitous, and if the environment remains unchanged, they can remain stable for several years [8,9]; it does not require the deployment of base stations and is low-cost.In addition, the magnetic field signals are not affected by multipath effects due to it being a physical field existing on the earth.So, magnetic field signals are effective sources of positioning that can assist PDR.Most of the positioning methods based on magnetic field signals are first to build the magnetic field signal database of the environment in advance and then realize the positioning using the magnetic field signal matching algorithm.Currently, the commonly used matching algorithms include Dynamic Time Warping (DTW) [10], Convolutional Neural Network (CNN) [11], Recurrent Neural Network (RNN) [12], etc. DTW is the most widely used magnetic field signal matching method with small computation and stable positioning performance.CNN has a strong ability for feature extraction of magnetic field signals, while RNN has a strong learning ability for sequence magnetic field signals; however, these existing methods take the magnetic field signal value directly as the input.As the third-generation neural network, Spiking Neural Network (SNN) [13] has the characteristics of event-driven and low-power consumption, and has attracted widespread attention.The neuron model in the traditional neural networks is based on numerical calculation, while SNN transmits information by action potential sequence rather than numerical value, which is a neuromorphological model with biological significance.However, SNN is difficult to train due to its complex neuron dynamics and indistinguishable spike operation.With the development of SNN, research and applications based on SNN are also increasing.Gutig et al. [14] proposed an SNN for classification, called Templeton.This method defined the Leaky Integrateand-Fire (LIF) neuron model and gave the learning mechanism based on neuron potential.Fang et al. [15] proposed a learning algorithm based on back-propagation using spiking neurons with learnable membrane parameters, called Parameter Leaky Integrate-and-Fire (PLIF) spike neurons.This method better represents the heterogeneity of neurons, thus enhancing the performance of SNN.Kim et al. [16] presented the first spiked-based object detection model, called Spiking-YOLO.The method introduced channel-wise normalization and signed neurons with imbalanced thresholds, both of which provide fast and accurate information transmission for deep SNN, making its detection accuracy up to 98%.Thus far, SNN is mainly used for image classification and object detection.Inspired that the magnetic field sequence signals can be mapped to the image, this paper first proposes using SNN to detect magnetic field landmarks to assist PDR.Nevertheless, magnetic field landmarks do not have global uniqueness because similar magnetic field features can be found in different places.
Digital Terrestrial Multimedia Broadcasting (DTMB) signals have low frequency, and stable transmission over a large range in cities, which is a good signal of opportunity (SoP) for positioning [17,18].The DTMB signals show the following obvious advantages.(1) The DTMB signal transmission is continuous and stable.The signal transmitter is fixed in the DTMB network, so it is not affected by the Doppler frequency shift, which will lead to easier signal processing and the possibility of integration over a longer period time [19].(2) The DTMB system has a strong anti-multipath interference ability.The Time Domain Synchronous Orthogonal Frequency Division Multiplexing (TDS-OFDM) is used in DTMB mode.TDS-OFDM can process Pseudo-random Noise (PN) sequences of multiple OFDM frames at the same time so that the delay length of anti-multipath interference is not limited by the length of the guard interval [20].(3) The infrastructure of DTMB has been deployed in major cities in China, which makes it easier to obtain as a standard operating procedure, and no more infrastructure investment is required for receiving equipment.Moreover, there have been many studies on fusing positioning based on DTMB signals.For example, Dai et al. [21] reported a positioning method that adopts the Time Domain Synchronous Orthogonal Frequency Division Multiplexing (TDS-OFDM) signals in singlefrequency networks, and the simulation positioning accuracy is within 1 m.Cong et al. [22] presented a positioning method using the GPS, DTMB, and Frequency-Modulated (FM) signals, which achieved outdoor positioning accuracy of 10.56 m.Jiao et al. [23] used DTMB signals to aid GNSS positioning.They first proposed a ranging method based on DTMB signals and then used the DTMB signals ranging information and GNSS pseudorange measurement to obtain fusing positioning, and the mean position errors are 6.23 m.Liu et al. [24] proposed an enhanced PDR based on DTMB signals, and the positioning error cumulative distribution reaches 2.36 m at a 68% probability level.All these works above have demonstrated the positioning potential of the DTMB signals.
Therefore, to further enhance the positioning performance of the PDR, this paper proposes a robust PDR method integrating magnetic field signals and DTMB signals.This paper uses an SNN for magnetic field landmark detection.To avoid the interference of highly similar landmarks, the DTMB ranging signal is introduced as a reference.Firstly, compare the distance between the candidate magnetic field landmarks and the DTMB base station with the DTMB signal ranging results, and then select the landmarks with a difference less than the threshold as the optimal landmarks.Finally, the position information from the optimal landmark and the ranging information from the DTMB signal are fused with the PDR to improve positioning performance.The contributions of the paper are as follows.(1) In this paper, a magnetic field landmark detection method based on SNN is proposed.The magnetic field signals are encoded based on the discharge characteristics of biological neurons, and then SNN is used to learn the encoded magnetic field signals.
Finally, the learning model is used to detect the magnetic field landmark.(2) We selected optimal landmarks for candidate magnetic field landmarks based on DTMB ranging information.To achieve robust PDR, the position information from the optimal landmarks and the ranging information from the DTMB signal are fused with the PDR using the Extended Kalman Filter (EKF) to achieve the fusion positioning.(3) To validate the proposed method, we tested it indoors and outdoors in the real world, and the outdoor test results showed that the positioning error cumulative distribution of the proposed method reaches 2.84 m at a 68% probability level.The proposed method will provide an alternative solution for seamless indoor and outdoor positioning.
The paper is organized as follows.Section 2 introduces the magnetic field landmark detection method based on SNN in detail.Section 3 describes the ranging method based on DTMB signals.Section 4 gives a fusion positioning method based on magnetic field signal landmarks, DTMB ranging signals, and PDR.Section 5 describes the tests and results.The conclusion is summarized in Section 6.The framework of this paper is shown in Figure 1.

Magnetic Field Landmark Detection Based on SNN
The detection of magnetic field landmarks first needs to build the magnetic field landmark database of the environment in advance, then use the SNN to learn the magnetic field signals, and, finally, detect the magnetic field landmarks based on the learning model.

Construction of Magnetic Field Landmark Database
To save the cost of establishing the magnetic field landmark database, this paper dynamically collects magnetic field signals.After collecting the magnetic field signals, this paper adopts the Butterworth filter to eliminate the noise.When the tester holding the smartphone walks at a constant speed, the triaxial components of the magnetic field signal have obvious differences at the same position due to the pose change of the smartphone.So, this paper uses the modulus of the triaxial component of the magnetic field signal to build the database.The calculation formula for magnetic field signal modulus is as follows, where B m indicates the module of triaxial magnetic field signals and B x , B y , and B z indicate the x-, y-, and z-axis magnetic field signals.The accurate construction of the magnetic field landmark database is an important condition to ensure the accuracy of magnetic field matching and positioning.Magnetic field signal landmarks can be divided into intrinsic landmarks and organic landmarks based on whether they contain visible ferromagnetic substances.Among them, inherent landmarks refer to visible objects containing ferromagnetic material, including fixed structures in buildings such as elevators, escalators, stairs, water dispensers, iron doors, etc.; organic landmarks refer to non-visible objects containing ferromagnetic material.Although the specific structure of the landmark cannot be seen, the magnetic signal will produce long-term stable distortion features when passing through the landmark.The construction of a landmark database includes three steps, extracting magnetic field feature sequences, mapping magnetic field images, and expanding the magnetic field image database.
The specific steps for extracting magnetic field feature sequences are as follows, (1) detecting peaks and valleys in magnetic field signal sequences; (2) using the detected peaks/valleys as the starting point, detect the first point after/before the peak that starts to increase (or the first point after/before the valley that starts to decrease) from the left and right directions, respectively; (3) use the first increase/decrease point discovered on both sides as the window radius, and extend the window radius in the left and right directions from the starting point (peaks/valleys) to extract feature sequences; and (4) calculate the change rate of magnetic field signal intensity within the window and set a threshold.If the change rate exceeds the threshold, the sequence will be recorded as a landmark.Otherwise, return to (2) until the data traversal is completed.
The specific steps for magnetic field image mapping are as follows.
(1) Interpolating the magnetic field signal sequence to increase sampling points to avoid sparse pixel points on the image; (2) subtract the minimum value from the magnetic field signal strength, and then multiply it by a scale parameter (in this paper, it is 100); (3) the difference between the maximum and minimum values of the magnetic field signal strength is mapped as columns of the image, and the sampling point index of the magnetic field signal are rows of the image, which are assigned one by one to complete the image mapping; and (4) annotate the position of the mapped image and store it as a landmark.
After mapping the magnetic field image, the changing trend of the magnetic field signal sequence can be extracted, which weakens the impact of different smartphones during the landmark detection process (different smartphones often have the same changing trend, but the overall magnetic field signal strength will have a deviation [8]).At the same time, converting magnetic field signals into images can effectively perform two-dimensional spatial calculations.Considering the mismatch in walking speed during data collection, using magnetic field signal images for detection may result in a decrease in accuracy.Therefore, this paper expands the database by stretching magnetic field images.This paper considers the differences in walking speed between different tests when setting the stretching factor.After empirical testing, the stretching factor is set to 0.5∼1.5.

SNN Architecture
In this paper, SNN is used to learn and detect magnetic field landmarks.The SNN architecture is shown in Figure 2.After mapping the magnetic field image, we use the phase encoding method to convert the image into spike signals that conform to the characteristics of biological neurons and input them into SNN for model training.In the process of magnetic field landmark detection, the neural network module is directly called to output the position of the matching landmark.

Neuron Encoding Model
SNN simulates biological neurons receiving spike signal sequences as inputs.To design an efficient SNN architecture, it is necessary to adopt appropriate spike encoding methods to encode magnetic field images into spike signal sequences.This paper uses phase encoding to encode magnetic field images into spike signal sequences of neurons.The phase encoding formula is shown as where A is the amplitude factor, B is the period factor, t is the current time of the spike, φ i represents the phase of the i th pixel, and it is calculated as follows where φ 0 is the initial phase, and ∆φ i is the phase difference between two adjacent pixel points, calculated as follows where N RF represents the number of pixels in the time window.

Neuron Dynamics Model
The neuron dynamics model is the basic computing unit of the neural network.LIF is a typical neuron dynamics model.It can capture the key dynamics of the input and output characteristics of biological neural networks, so it has been widely used in neuron computing.In this paper, the LIF neuron dynamics model is used as the basic computing unit of the SNN.For the nervous system, the information transmission of neurons can be expressed as the transmission of membrane potential from presynaptic neurons to postsynaptic neurons.The membrane potential of a postsynaptic neuron is determined by the membrane potential of all presynaptic neurons connected to the neuron.The membrane potential of postsynaptic neurons is calculated as follows [14] where η is the number of input spike; i is the index of the input spike; t represents the current time; t i is the release time of the i th input spike; w i represents the weight of the i th input synapse; V rest is the resting potential of the current neuron; and K(t − t i ) is the kernel function representing the delay effect.The specific calculation formula is as follows: Among them, τ indicates the delay parameter of membrane integration current; τ s is the delay parameter of synaptic current; and V 0 is the normalization factor.
According to the ( 5) and ( 6), the membrane potential intensity changing trend can be described as that after the peak value is generated, the membrane potential intensity rises from V rest to V th , and then gradually decays and tends to V rest .When the potential of the neuron reaches V th , the spike will be activated.However, the transmission of the spike is not only related to its release time but also related to the strength of the synaptic connection.When the synaptic weight w i is positive, the presynaptic potential has a longterm enhancement effect; when w i is negative, the presynaptic potential will produce a long-term inhibitory effect.

Neuron Learning Mechanism
The cost function of neurons in the learning process is defined as where V th is the spike release threshold and t max is the maximum release time.The derivative of the cost function over weight is: The updated rules for synaptic transmission weights are as follows: Among them, parameters λ represent the learning rate and describe the magnitude of synaptic weight adjustment.

Ranging Based on DTMB Signals
DTMB is a Chinese standard [25], and together with the Advanced Television System Committee (ATSC) of the United States and the Digital Video Broadcasting (DVB) of Europe become the main digital broadcasting system.In this paper, the DTMB signals are used for ranging estimation.

DTMB System Description
The DTMB system uses OFDM modulation to achieve robust transmission in multipath scenarios.The DTMB system supports single-carrier and multi-carrier modulation.This study focuses on multi-carrier modulation.The DTMB signals are transmitted in streams and organized in frames.The DTMB transmitting system completes the conversion from the input data stream to the DTMB channel transmission signal, and the transmitting principle of the DTMB system is shown in Figure 3.The input data stream is encoded by scrambler (randomization), forward error correction coding, constellation mapping from bitstream to symbol stream, and interweaved to form a data block.After the data block and system information are multiplexed, the frame body is formed by frame body data processing.The frame body and the corresponding frame head (PN sequence) are multiplexed into a signal frame and converted into a baseband output signal with bandwidth within 8 MHz after baseband post-processing.Finally, the signal is converted into an RF signal (UHF and VHF) by orthogonal upconversion.The frame structure of the DTMB system is shown in Figure 4, which is a four-layer structure.The basic unit of the frame structure is the signal frame, which is composed of the frame head and the frame body.The frame head is a PN sequence with a preamble and a postamble.The PN sequence includes three different length modes, 420, 595, and 945, at the sampling rate of 7.56 MHz.In our study, the length mode of the PN sequence is 945.The frame body is an Inverse Discrete Fourier Transform (IDFT) block with 3780 subcarriers, 36 of which are used to transmit Transmission Parameter Signaling (TPS) to provide demodulation and to decode information at the receiver.The superframe consists of 200 signal frames.In a superframe, each PN sequence is different from the other 199, but it is repeated in every superframe.The minute frame is defined as 480 superframes that last 1 min.The top layer of the frame structure is called the calendar day frame, starting from 00:00:00 and ending at 24:00:00.The signal structure is periodic and keeps synchronization with natural time.

Ranging Estimation
The DTMB signal system adopts Time-Domain Synchronous Orthogonal Frequency Division Multiplexing (TDS-OFDM) mode.Assuming X(k) represents the k th subcarrier in the OFDM system, and k = 0, ..., N − 1.Then, the transmitted l th OFDM sample is described as [26]: After the Fast Fourier Transform (FFT), the k th subcarrier of the m th signal frame received in the frequency domain is expressed as [27]: where N s is the total sample of a signal frame, N g is the length of the guard interval, α ∆T/T is the normalized Sampling Frequency Offset (SFO), T is the sampling period, ∆ f NT is the normalized Carrier Frequency Offset (CFO), C(k) is the channel frequency response, w m (k)∼N(0, σ 2 ) is the Gaussian noise, A N (k) is the attenuation factor and I m (k) represents the Inter-Carrier Interference (ICI): The data modulated on the TPS can remain unchanged for a long time, such as for dozens of days.Consequently, the CFO of the m th signal frame ∆f m can be obtained as [24]: The CFO is mainly caused by oscillator mismatch between the transmitter and receiver and the Doppler effect.It can be defined as where f d m is caused by the Doppler effect, f s m is caused by oscillator mismatch between receiver and transmitter, and q m is the Gaussian noise.Generally, f s m is stable, so suppose it remains unchanged during a test.Thus, the f s m can be obtained by holding still for a few seconds after the test starts.The Doppler speed ∆v dtmb m+1,m between the (m+1) th signal frame to the transmitter and the m th signal frame to the transmitter is expressed as [28] where c is the light speed and f c is the carrier frequency.Finally, the range difference ∆r dtmb m between the m th signal frame to the transmitter and the first signal frame to the transmitter is expressed as:

Hybrid PDR Integrating Magnetic Field Signals and DTMB Signals
Magnetic field landmarks cannot provide unique detection results, because similar magnetic field features can be found in different places.To eliminate magnetic field landmarks with significant anomalies, this paper proposes an optimal landmark selection method based on DTMB ranging information, and then the optimal landmarks are used for fusion positioning.

Optimal Landmark Selection
Assume that P mag = P mag 1 , P mag 2 , . . ., P mag n is the position coordinates of the candidate landmarks obtained based on magnetic field signal matching, where n represents the number of landmarks.Then, compare the distance between the candidate magnetic field landmark and the DTMB base station with the DTMB signal ranging results, and then select the landmarks with a difference less than the threshold parameter as the optimal landmarks.The calculation formula for the optimal landmark P mag is as follows where δ r is a distance difference threshold, P base is the position coordinate of the DTMB base station, and ∆r dtmb is the ranging result from DTMB signals.

Hybrid Positioning Model
In this paper, we use the Attitude and Heading Reference System (AHRS) [29,30] to obtain the raw PDR trajectory.To achieve robust PDR, we use the position information from the optimal landmark and the ranging information from the single base station DTMB signals as the observation, and then use the EKF to achieve the fusion positioning.The state model of the fusion framework in this paper is defined as: The state vector is defined as: where t represents the current time, x is the position of the eastern direction, y is the position of the northern direction, θ is the heading change, θ is the heading, and v is the motion velocity.The state equation is defined as follows where ∆t is the time interval between the time t − 1 and t, w t = {w 1 , w 2 , ..., w 5 } t is the noise vector of the state equation.The state transition matrix F t is written as: The process noise matrix W t [24,31] is expressed as where a = cosθ t−1 and b = sinθ t−1 .The observation vector Z t is defined as: Among them, x mag and y mag are the positions from optimal magnetic field landmark, which is detailed in Sections 2 and 4.1.∆r dtmb is from the DTMB signals, which is detailed in Section 3. The heading θ pdr and speed v pdr are estimated by AHRS based on Micro-Electro-Mechanical Systems (MEMS) sensors.The observation model is defined as: The observation equation can be expressed as: The r t represents the range from the pedestrian to the DTMB transmitter and it is expressed as where x base , y base , and z base are the position coordinates of the DTMB signal transmitter, respectively; z t is the height of the pedestrian movement plane, which is set to 0. The observation Jacobian matrix H t is shown as: The observation noise is δ t ∼N(0, R t ) and the observation covariance matrix R t is:

Tests and Results
To validate the proposed method, we conducted field tests in both indoor and outdoor environments.Among them, there is no DTMB information for indoor testing, mainly testing the feasibility of magnetic field landmarks.

Indoor Tests
We collected indoor test data at the Xinghu Building of Wuhan University.The floor map of the indoor data collection scenario is shown in Figure 5.The collected mobile device is the Huawei P30.To achieve mutual verification of multiple sets of tests, the testers walked along the testing route three times.We use three tests for mutual verification to test the accuracy of landmark detection.During testing, one of the tests is used as the test set, and the other two tests are used as the magnetic field landmark database for SNN model matching.Then, the intersection of the matching results is extracted as the final result of this test.To verify the performance of SNN, we use the DTW [32] method as a comparative test.The matching results for each test set of the two methods are shown in Figure 10.Among them, Figure 10a shows the results of the DTW method, and it can be seen that a total of nine landmark intersections were matched in test 1, while eight landmark intersections were matched in both test 2 and test 3. Figure 10b shows the results of the SNN method, and it can be seen that a total of seven landmark intersections were matched in test 1, while six landmark intersections were matched in both test 2 and test 3. Comparing the matching results of the two methods, it can be concluded that SNN matching results are more accurate than DTW overmatching.This is important for magnetic field matching, as similar magnetic field sequences are often in different positions.Overmatching will lead to landmark recognition errors and significant position deviations.After further selection, we selected six intersection landmarks from three tests based on the matching results of the SNN method, as the final landmarks.Meanwhile, we also extracted the turning points of the motion trajectory and used them as the calibration points for the final PDR trajectory, along with the magnetic field landmarks.The PDR trajectory calibration results are shown in Figure 11.It can be seen that after landmark calibration, the PDR trajectories of the three tests can almost overlap, which is in line with the actual motion state.

Outdoor Tests
We organized a real-world outdoor test to verify the proposed method.The test address is located in Lianhua Lake Park, Wuhan, China.During the test, the tester followed the planned route from the starting point to the endpoint, then returned to the starting point, walking for two laps.The trajectory is about 116 × 2 m, or 232 m. Figure 12 is a top view of the remote sensing of the real-world outdoor test.To adapt to the SNN structure, the magnetic field image is converted into a spike sequence, as shown in Figure 16.We use the landmarks in the first lap as the database to train the SNN model, and then perform matching tests on the landmarks in the second lap.The matching results are shown in Figure 17, indicating that five magnetic field landmarks have been matched.

The Ranging Results from DTMB Signals
The ranging results based on DTMB signals are shown in Figure 18.In this paper, the ground position reference is first obtained through U-blox-M8U, and then obtained through RTK processing using data from the Wuhan IGS station (jfng) as the reference station.The distances from the ground position reference to the DTMB signal transmitter are calculated as a ranging reference, and then the ranging errors of the DTMB signal are counted.After statistics, the median error of DTMB ranging is 1.79 m and the maximum error is 5.92 m.

The Fusing Positioning Results
The outdoor field test results for the three positioning systems are shown in Figure 19, including ground reference, PDR, and the proposed method.It can be seen from Figure 19 that the trajectory direction calculated by the PDR deviates significantly compared to the ground reference, especially for the second lap of trajectory, where the heading of the PDR drifted severely.This is mainly due to the unstable attitude of the smartphone during the test.It is worth noting that after the fusion of the magnetic field landmarks and the ranging results from DTMB signals, the proposed method greatly corrects the PDR positioning trajectory.Then, magnetic field landmarks were optimized based on the DTMB ranging information to eliminate the interference of similar magnetic field landmarks in positioning.Finally, EKF is used to fuse the position information from the optimal magnetic field landmarks and the ranging information from DTMB signals.The proposed method has been validated to be effective in both indoor and outdoor environments with positioning accuracy at the meter level.It will provide an alternative solution for seamless indoor and outdoor positioning.In future work, we will explore the indoor fusion positioning based on DTMB signals and integration positioning with other signals, such as WiFi, UWB, LTE, 5G, and audio.

Figure 1 .
Figure 1.The framework of this paper.

Figure 3 .
Figure 3.The transmitting principle of the DTMB system.

Figure 4 .
Figure 4.The frame structure in DTMB system.

Figure 5 .
Figure 5. Indoor floor map.After data collection is completed, the magnetic field sequence signals of the three tests are smoothed and feature sequences are extracted to obtain a landmark database for

Figure 6 .
Figure 6.Indoor magnetic field data processing results.Using the smoothed triaxial modulus for indoor landmark feature extraction, the extraction results are shown in Figure 7.Among them, nine landmarks were extracted in test 1, and eight features were extracted in test 2 and test 3.The results of image mapping from landmark feature curves are shown in Figure8, which shows that the curve features of the magnetic field signals are well preserved in the form of images.To adapt to the SNN structure, the magnetic field image is converted into a spike sequence, as shown in Figure9.

Figure 7 .
Figure 7.The indoor landmark feature extraction results.

Figure 8 .
Figure 8.The indoor landmark image mapping results.

Figure 9 .
Figure 9.The spike sequences of the indoor landmark images.

Figure 10 .
Figure 10.Landmark detection results of three indoor tests.

Figure 11 .
Figure 11.The positioning trajectory of three indoor tests.

Figure 12 .
Figure 12.The remote sensing of the real-world outdoor test.5.2.1.The Magnetic Field Landmark Detection ResultsThe original magnetic field signals and the smoothed magnetic field signals of the outdoor collected route are shown in Figure13.Using the smoothed triaxial modulus for landmark feature extraction, the extraction results are shown in Figure14, where numbers 1-8 represent the landmark features detected in the first lap and numbers 9-16 represent the landmark features detected in the second lap.The results of image mapping from landmark feature curves are shown in Figure15, which shows that the curve features of the magnetic field signals are well preserved in the form of images.

Figure 13 .
Figure 13.Magnetic field data processing results.

Figure 14 .
Figure 14.The outdoor landmark feature extraction results.

Figure 15 .
Figure 15.The image mapping results of the outdoor landmark feature.

Figure 16 .
Figure 16.The spike sequences of the outdoor landmark images.

Figure 17 .
Figure 17.Landmark detection results of the outdoor field test.

Figure 18 .
Figure 18.The DTMB signals ranging results of the outdoor field test.

Figure 19 .
Figure 19.The positioning trajectory of the outdoor field test.The cumulative distribution function (CDF) of the PDR and the proposed method are shown in Figure 20.The CDF results suggest that the 68% positioning errors of PDR and the proposed method are 8.77 m and 2.84 m, and the 95% positioning errors are 12.84 m and 3.98 m, respectively.These statistical data indicate that the proposed method has superior positioning accuracy compared to PDR.Outdoor field testing has verified the positioning effectiveness of the proposed method.

Figure 20 .
Figure 20.The CDF of the positioning results.
Magnetic field signals are stable and ubiquitous, and DTMB signals have strong penetration and stable transmission over a large range.To improve the positioning performance of PDR, this paper proposes a robust PDR integrating magnetic field signals and DTMB signals.Inspired by that magnetic field signals can be mapped into images, this paper innovative converts magnetic field sequence signals into images, and first uses SNN to learn the magnetic field landmark images of the environment, and then uses a learning SNN model for magnetic field landmark detection.Meanwhile, the DTMB signals are collected by the self-developed signal receiver, and the ranging estimation based on the DTMB signals is achieved.