A Robust Real Time Direction-of-Arrival Estimation Method for Sequential Movement Events of Vehicles

Parameters estimation of sequential movement events of vehicles is facing the challenges of noise interferences and the demands of portable implementation. In this paper, we propose a robust direction-of-arrival (DOA) estimation method for the sequential movement events of vehicles based on a small Micro-Electro-Mechanical System (MEMS) microphone array system. Inspired by the incoherent signal-subspace method (ISM), the method that is proposed in this work employs multiple sub-bands, which are selected from the wideband signals with high magnitude-squared coherence to track moving vehicles in the presence of wind noise. The field test results demonstrate that the proposed method has a better performance in emulating the DOA of a moving vehicle even in the case of severe wind interference than the narrowband multiple signal classification (MUSIC) method, the sub-band DOA estimation method, and the classical two-sided correlation transformation (TCT) method.


Introduction
Intelligent transportation and unmanned systems have made new demands to the parameter estimation of sequential movement events of vehicles, such as noise insensitive, scalable, and portable implementation, etc. Fortunately, microphone arrays bring new potential solutions to deal with these problems. Microphone arrays are widely used in civilian and military fields, such as noise reduction and signal enhancement [1][2][3][4], hearing aid [5,6], outdoor wind noise measurement [7,8], acoustic source localization and tracking [9][10][11][12], weapon classification and shooter localization [13,14], and so on. In the last few years, due to the development of the micro-electro-mechanical system (MEMS) technology, small aperture arrays have shown their potential superiority in terms of portability and easy deployment, especially in the field. Therefore, small aperture arrays have received more and more attention, both in industry and academia.
Recently, a small aperture MEMS microphone array system has been introduced in our earlier work [15] for direction finding for the sequential movement events of vehicles. Moreover, we have conducted a series of researches about small aperture microphone arrays and obtained some positive results, for example, acoustic target intrusion detection method [16,17], acoustic target classification method [18], and acoustic source localization method [19,20]. However, the robustness of these methods needs to be further improved in the field, especially under the influence of wind noise. To get a good bearing tracking performance, a robust direction-of-arrival (DOA) estimation method is crucial. In this paper, we will propose a wind noise robust DOA estimation method based on a small aperture microphone array system.
The mechanism and power spectrum of outdoor wind noise have been deeply studied in many works [21][22][23][24][25]. Morgan and Raspet have concluded, "The dominant source of wind noise in outdoor microphones is the pressure fluctuations that are caused by the velocity fluctuations of the incoming flow". References [21] and [22] have indicated that wind induced noise focuses on low frequency, especially below 1 kHz. Moreover, the acoustic detection and tracking of ground vehicles also use the signal of low frequency band. Therefore, it is a challenging problem to design a wind noise robust DOA estimation method for the sequential movement events of vehicles.
From [23], wind noise changes the characteristics of the low-frequency bands of the detected acoustic signal quickly and temporally. It is difficult to separate the target acoustic signal and the wind noise from a single channel acoustic signal, both in the time domain and the frequency domain. Many researches show that the wind noise is nearly incoherent [8,25]. Fortunately, we can distinguish them by computing the magnitude squared coherence (MSC) of two selected microphones.
Independent from apertures, the MSC of the wind noise (0~4 kHz) is close to zero, while the MSC of the corresponding target acoustic signal is close to 1 [26]. Due to this difference, the spatial coherence is used to select the frequency, which is less affected by the wind noise and to estimate the DOA of a moving vehicle in [19]. Moreover, the MSC is extended from wideband to sub-band, and meanwhile the sub-band magnitude squared coherence (SMSC) DOA estimation method is designed for the vehicle bearing tracking. However, we found that it is not enough to estimate DOA robustly only by using the sub-band with maximum SMSC value alone.
In this paper, we propose a multiple high SMSC sub-bands weighting strategy and integrate it into the moving vehicle DOA estimation method of our previous work [20]. As we know, most of the DOA estimation algorithms can be divided into two classes: Maximum Likelihood (ML) based and covariance matrix decomposition based. In recent years, most of the DOA estimation research works are concentrated on ML based algorithms such as [27][28][29][30]. Although the ML based methods can get good estimation accuracy when the signals are correlated, it needs a huge computation and long run time, which is not proper for a real-time embedded platform using in the wild environment. For another class of DOA estimation methods, we have compared their advantages and disadvantages in our previous work [19]. As a result, the MUSIC method is selected to estimate the DOA after sub-bands are weighted in our developed microphone array system.
To verify the proposed method, we have carried out field experiments to check out the actually achievable performance. The results show that the proposed method has better robustness to the uncorrelated wind noise than the narrow-band MUSIC method, sub-band method [20] and classical wideband TCT method [31].
The remainder of this paper is organized as follows: Section 2 describes sequential data modeling for moving vehicles bearing tracking. Section 3 proposes the SMSC and briefly introduces the famous MUSIC method and the incoherent signal-subspace method, which would be used in the fourth section and the fifth section. Section 4 mainly analyzes the SMSC characteristics of the acoustic signals of the moving vehicles and the wind noise. Section 5 describes our proposed weighted sub-band DOA estimation method. Section 6 first introduces our designed small aperture MEMS microphone array system and then shows the field experiment results to verify the good performance of our proposed method. Finally, Section 7 concludes this paper.

Sequential Data Modeling for Bearing Tracking
Bearing tracking of sequential movement events of vehicles by a microphone array can be approximated as a sequential data modeling. Figure 1 shows the data observation scene for multiple targets passing a Microphone Array (MA) node. Assuming that the velocity of a vehicle (denoted as v) is uniform in the sampling instant time, where t c represents the moment that a target is passing the closest point of approach (CPA) and d represents the distance between the MA node and the lane center. Therefore, the DOA of a target passing the MA node can be described as Equation (1).
where t 0 , t 1 , . . . , t k represents multiple discrete observation time in the whole process of bearing tracking, and θ(t) is the observed azimuth at time t. Note that θ(t) is relative to the horizontal direction from left to right. Without any interference, an ideal DOA curve of three vehicles successively passing a MA node from left to right is shown in Figure 2.  As can be seen from Figure 2, the DOA ranges from 0 to 180 • in this situation. When a target is far at the left side, the DOA is nearly 180 • and it is gradually decreased to 0 when the target is travelling to the right side of the MA node. The ideal DOA curve is very smooth when there is no interference. However, in the real world, there always exists some interference, including instrument noise, wind disturbance, which greatly declines the DOA estimation precision. Therefore, it is especially important to improve the robustness of the DOA estimation method in the wild environment.

Sub-Band Magnitude-Squared Coherence (SMSC)
We first list some notations that will be commonly used in this paper.

•
The matrices are indicated by uppercase bold letters and vectors by lowercase bold letters.

•
The superscript * denotes the conjugate of a complex number.

•
The superscript H denotes the conjugate transpose of a matrix.

•
The superscript T denotes the transpose of a matrix.

•
The italic E[·] denotes the statistical expectation. • f denotes the frequency of Hertz and ω is the corresponding radial frequency, ω = 2πf.
The magnitude-squared coherence (MSC) between two signals x(t) and y(t) is defined as [26]: where p xy (f ) is the cross-power spectral density of x(t) and y(t). p xx (f ) and p yy (f ) are the auto-power spectral density, respectively. γ xy (f ) is the MSC value at frequency f for x(t) and y(t). The MSC is a frequently used signal processing technique that returns real values that are between 0 and 1 to indicate how relative between two time domain signals x(t) and y(t). If x and y are completely uncorrelated, γ xy (f ) = 0; if x is strictly correlative with y, then γ xy (f ) = 1. The sub-band MSC (SMSC) is used to evaluate the sub-band correlation. Suppose the wideband frequency of x(t) and y(t) can be divided into J uniformly distributed sub-bands, and then the ith sub-band MSC is defined as Equation (3) [20].
where f i is the central frequency of the ith sub-band. To calculate auto-power and cross-power spectral densities for x(t) and y(t), the Welch's averaged modified periodogram method [32] is typically used. The time-domain signals x(t) and y(t) are first divided into N time windows of the same length (i.e., with the same number of samples). Then, the auto-power and cross-power spectral densities in each window can be calculated by discrete Fourier transform (DFT). The final MSC and SMSC can be obtained by averaging all of the corresponding values of MSC and SMSC in N time windows.

The MUSIC Method
The MUSIC method is a well-known and classical signal parameters estimation method. Different from those methods that directly process the covariance matrix of received array data, such as minimum variance method [33,34], maximum entropy method [35], conventional beam forming method [36], it performs eigenvalue decomposition on the covariance matrix of the observed signal and divides this matrix into signal subspace and noise subspace. By making use of the orthogonality between these two subspaces, the MUSIC can estimate the incident direction of targets.
Consider P narrowband signals in the field far from the array are impinging on an M-element array, then the received array signals can be expressed as: . . . , s p (t)] T represents the incident signals of P targets. N = [n 1 (t), n 2 (t), . . . , n M (t)] T is the additive noise. A = [α(θ 1 ), α(θ 2 ), . . . , α(θ P )] is the array manifold matrix. Each column of A is also referred to as steering vector and represents the array response to each incident signal. The matrix element α ij depends on the ith array element, its position relative to the origin of the coordinate system and its response to the jth incident signal. Once the array structure is determined and the time delays between array elements, the array manifold matrix is almost fixed. Suppose the noise signals are Gaussian distributed and are irrelevant to the source signals, then the correlation matrix of received array data can be expressed as: Due to the mutual independence of noise and signal, the above covariance matrix can be divided into signal subspace and noise subspace, respectively. By performing eigenvalue decomposition of R, where U S ∑ S U H S represents the signal part and U N ∑ N U H N represents the noise part. Concretely, the covariance matrix can also be expressed as Ideally, the signal subspace composed of steering vectors is orthogonal to the noise subspace [37], which is expressed in the following equation: The MUSIC spectrum function is defined as: When Equation (7) is satisfied, the MUSIC spectrum P MUSIC at θ i is the largest. Theoretically, the MUSIC method is just using Equation (8) to distinguish different source incident signals by finding every peak of the MUSIC spectrum. Accordingly, each peak value relates to a DOA that can minimize the left part of Equation (7).
The above derivation is based on very ideal conditions and consumptions. In practical situations, the general method is to use some statistical techniques on the received signal. We can get many snapshots of the measured signal and use the statistical mean value to approximate the ideal value. Thus, the maximum likelihood estimation of the covariance matrix is given as: where L is the number of snapshots and X l represents the observed signal of the M-element array.
Similarly, by eigenvalue decomposition ofR, we can get the estimated noise subspaceÛ N . In practical situations, the steering vector α(θ) is not completely orthogonal toÛ N . However, we can search all of the possible θ to find the peak value of MUSIC spectrum. The final estimated DOA using MUSIC method can be expressed as

Incoherent Signal-Subspace Method (ISM)
The incoherent signal-subspace method (ISM) [38] is first applied to wideband DOA estimation by decomposing wideband signals into many narrowband signals. Then, narrowband DOA estimation methods can be used on each decomposed narrowband signal. The results from all of the sub-bands are combined to get the final wideband DOA estimation.
Suppose that the wideband signal is split into total J sub-bands in frequency domain. The MUSIC method is applied on each sub-band. Then, the estimated covariance matrix of the ith sub-band can be given by Equation (11).R Similar to MUSIC, by proceeding eigenvalue decomposition on every sub-band covariance matrix, each sub-band signal is divided into signal subspace and noise space. Besides, the sub-band MUSIC spatial spectrum can also be calculated according to Equation (8). Then, we can obtain the wideband MUSIC spatial spectrum by averaging all sub-bands spatial spectrums, which can be expressed as

Analysis of Vehicle Acoustic Signal and Wind Noise
The measured vehicle acoustic signals that are acquired in the field are mainly composed of the target acoustic signal and ambient wind noise. References [21,25] have pointed out that the dominate source of wind noise in outdoor microphones is the pressure fluctuations caused by the velocity fluctuations of the incoming flow. This phenomenon can seriously affect the microphone and cannot be completely removed by a fine-designed wind shelter.
In order to analyze the main characteristics of vehicle acoustic signals and turbulent wind noise, several typical sample signals are collected from Chongming Island, Nanjing and other places of China mainland. To avoid the involving of other noise sources, the selected experimental places are flat, open, and far away from the residential areas. We construct a signal acquisition system with a laptop, DAQ signal acquisition card and a MEMS microphone array. The INMP504 (InvenSence, Sunnyvale, CA, USA), whose signal-to-noise ratio can reach up to 65 dBA, is used to build the microphone array. Moreover, the acquisition card NI9239 (National Instruments, Austen, TX, USA) can simultaneously perform a 4-channel 24-bit synchronous data collection, with no more than 50 kHz sampling rate. The actual sampling rate is set to be 8192 Hz because the acoustic energy of the vehicle mainly targets focuses below the frequency 2 kHz. The microphone array node is placed 10 m away from the road center as shown in Figure 1. The wind speed is recorded by an ultrasonic anemometer, which is deployed beside the array during the experiments.
The original acoustic signal and its power spectrum of a passing vehicle without wind are shown in Figure 3a,b, respectively. We divide the observed wideband signal, whose frequency ranges from 0 to 4096 Hz, into equally spaced 64 sub-bands. Figure 3c shows the SMSC of each sub-band acoustic signal from two different microphones. Correspondingly, the three kinds of signals, as shown in Figure 3, are presented in Figures 4 and 5, respectively, for the cases when there are only wind noise and when there are both vehicle signal and wind noise.   According to Figure 3, when the vehicle is approaching the microphone array, the signal amplitude is increasing and the spectrum range is widening gradually. The signal amplitude reaches the peak when the vehicle moves to the nearest position, and spectrum range also extends to the full band. The original signal amplitude decreases gradually while the vehicle is moving far away from the microphone array. At the same time, the spectrum range also becomes narrow slowly. When the vehicle is approaching but far away from the microphone array, there are only a few low-frequency sub-bands but with high SMSC, as depicted in Figure 3c. When the distance between the vehicle and the microphone array gets closer and closer, more and more sub-bands have higher SMSC value. Almost all sub-bands have relatively high SMSC when the vehicle arrives at the nearest point to the microphone array. When the vehicle is leaving far away from the microphone array, only low-frequency sub-bands have high SMSC again.
As shown in Figure 4, the power spectrum of wind noise in low frequency region changes rapidly at different times. However, the SMSC of wind noise is quite stable and very close to zero, especially in low frequency bands. Seeing from Figure 5, the vehicle acoustic signal is polluted by wind noise. In this case, the power spectrum is greatly damaged, especially in high frequency bands. However, the SMSC is not influenced that seriously. Although the SMSC of some sub-bands, especially the low frequency sub-bands, is nearly 0, the SMSC of most bands are still approaching to 1, which can easily be used to distinguish the source signal from wind noise.
From Figure 5, we can see that the wind noise can severely damage the original vehicle acoustic signal. When the selected frequency is greatly deteriorated, the DOA estimation result will deviate from the true value when using narrowband DOA estimation methods (such as narrowband MUSIC). Other classical methods, such as TCT, is relatively robust, the DOA estimation performance still degrades heavily when several non-continuous frequencies are damaged. The traditional sub-band method uses spatial spectrums of all the sub-bands. However, as shown in Figure 5c, some sub-bands are greatly damaged. These contaminated sub-bands can destroy the final DOA estimation result. Fortunately, there still exist some high SNR sub-bands. We can select these sub-bands with high SNR to estimate DOA which is similar to [9,39]. The main difference is that we use SMSC to evaluate the SNR. The larger the SMSC is, the higher the SNR is, and vice versa. Thus, we can search the SMSC of all the sub-bands and select a few top sub-bands with high SMSC. By weighting these sub-bands, we can output more accurate DOA estimation results. The next section will carefully present our proposed method.

The Weighted Sub-Band DOA Estimation Method
In Section 4, we analyze the spatially correlated characteristics of target acoustic signals under the wind noise. Inspired by the idea of wideband DOA method for ISM estimation, we also divide the target acoustic signal into several sub-bands in the frequency domain. The turbulent wind causes different sound pressure fluctuation on each sub-band. Therefore, the sub-band that is severely affected by the wind noise will have lower SNR and SMSC, while the sub-band is less affected by the wind noise will have rather higher SNR and SMSC. If the contaminated sub-band is used to estimate the DOA, the performance will be seriously deteriorated. In order to dismiss those low SNR sub-bands, we set a threshold TH for the SMSC. Thus, only those sub-bands that exceed the threshold are selected out. These sub-bands, which have relatively high SNR, are then weighted to calculate the spatial spectrums. The final estimated DOA of a moving vehicle can be found out by searching the spectrum peak on the selected sub-band spatial spectrums.
Suppose that there are K sub-bands selected for weighting from total J sub-bands, and then the weight of the selected ith sub-band ω i can be calculated by Equation (13): where γ xy (f (i) ) denotes the ith SMSC, which is calculated by Equation (3). When considering that the high SMSC usually has larger SNR and can give a high degree of confidence for DOA estimation, thus in our method, the sub-band with higher SMSC is assigned a larger weight. By weighting the sub-bands with high SMSC, it is inclined to get more accurate estimated DOA and to improve the robustness compared with other methods based on MUSIC. The final MUSIC spatial spectrum using the weighted sub-bands is then calculated by Equation (14): where α(f (i) ,θ) denotes the steering vector of the signal subspace of the ith sub-band. When all of the J sub-bands are selected out and the received signals have the same correlation on all sub-bands (all sub-bands have the same SMSC), then ω i = 1/J and the above function will have the same form of the typical ISM, as expressed in Equation (12). In other words, the proposed weighted sub-band DOA estimation method is a generalized method for sub-band DOA estimation and ISM is a special form of ours.
In particular, we adopt ISM to implement weighting process on multiple sub-bands in this paper. Similarly, the weighting process also can apply to the Coherent Signal-subspace Method (CSM). Since the MEMS microphone array aperture is only 68 mm, to avoid the spurious peaks during implementation of the MUSIC method, the spectrum range is limited to no more than 2.5 kHz. When compared with the previous work [19,20], the proposed method improves the robustness of the target DOA estimation using narrowband MUSIC in the wild environment.
The block diagram of the weighted sub-band DOA estimation method is shown in Figure 6. In accordance with Figure 6, our proposed method can be divided into the following steps: 1.
Calculate the frequency domain signals for all M arrays using FFT on the original observed time domain signals.

2.
Select two channels of microphone signals and calculate their SMSC for all of the sub-bands in the frequency domain. Sort all of the SMSC in descending order, and select the best K sub-bands with K largest SMSC. K is the number of sub-bands, whose SMSC is greater than the threshold TH. TH is an empirical value.

3.
For the selected K sub-bands, calculate the weight for each sub-band.

4.
Perform eigenvalue decomposition to get the noise subspace for each sub-band. Estimate the number of sources according to the sources number estimation criteria, such as MDL criterion.

5.
Calculate the weighted MUSIC spatial spectrum from the K sub-bands and obtain the final DOA by spectrum peak searching.

Field Experiments
In order to verify the anti-wind noise performance of the proposed method, experimental studies have been performed on Chongming Island (the third biggest islands in China) and a suburban district around Nanjing. The field test scene is shown in Figure 7 when a vehicle is passing the MEMS microphone array system. Figure 8 shows the PC (Personal Computer) user interface, which can demonstrate the real-time DOA by UART in a LabVIEW 8.5 programming environment.  First, we will introduce our hardware design for the MEMS microphone array system. Then, the experimental results using this system is described.

Hardware Architecture of The MEMS Microphone Array System
Our MEMS microphone array system consists of the following four units: MEMS microphone array, preprocessing and sampling unit, processing or acquisition unit, and display and control unit. The block diagram of the system is shown in Figure 9. We choose a 4-element uniform circle array (UCA) with a diameter L = 68 mm. The MEMS microphone array converts the collected acoustic signals to electric signals and sends them to the preprocessing and sampling unit. This unit can complete various acoustic signal processing, including filtering and amplification, then through the analog digital converter (ADC) to complete the synchronous sampling conversion of multiple acoustic signals. There are two operating modes: processing or acquisition. When in mode 1, the target detection, classification, and DOA estimation is executed. When in mode 2, only four channel acoustic signals are sampled and sent out in real time. Accordingly, the display and control unit also has two operating modes. In mode 1, the target information is transmitted to a personal digital terminal device and the original data is transferred via USB cable to a portable laptop in mode 2.
As shown in Figure 10, our microphone array system includes a main board, four array element boards, and a radio frequency (RF) board. The main board and the array element board are connected through a Flexible Printed Circuit (FPC), and the RF board is connected with the main board through a pin connector. The main board adopts a DSP ADSP21479 (Analog Devices, Norwood, MA, USA) as the core processor, and adopts a four-channel, 16-bit ADC MAX11043 (Maxim Integrated Products, San Jose, CA, USA). Four array element boards are placed evenly on the edge of the process board. All of the array element boards have the same PCB (Printed Circuit Board) layout and routing design. This kind of design can minimize the orientation error that is caused by the interchannel inconsistency between microphone arrays. Besides, it is easy to replace any damaged MEMS microphone after a long-term working in the field. The microphone chooses INMP504 (InvenSence, Sunnyvale, CA, USA), of which the signal-to-noise ratio (SNR) can reach 65 dBA. The RF board is responsible for transmitting the detected target information in mode 1. Figure 10. Photograph of the micro-electro-mechanical system (MEMS) microphone array system, array aperture is 68 mm.

Experiments Result
During the experiments, we compared the performance of our proposed method with narrowband MUSIC method [19], sub-band method [20] and classical TCT method. Each signal is split by sliding window, and the window length is 4096 with no overlap in 8192 Hz sampling frequency. The bandwidth of TCT, sub-band method, and weighted sub-band method are selected from (30~2500) Hz. Besides, to identify the useful sub-band of the signal, we check whether the spatial coherence is above the threshold. When considering the detecting distance of the vehicle, the threshold of weighted sub-band method is set to 0.5. Figure 11a shows the acoustic signal when the vehicle passes through the location of the microphone array system. Figure 11b shows the correlation of each sub-band in the process of vehicle, and the SMSC of the frame are approximately equal to 1, then the proposed weighted sub-band method is converted to the ISM method. Figure 11c shows the DOA estimation results of the four methods in the process of a vehicle passing under the condition that no wind disturbance. At this time, the four curves coincide perfectly. In order to test the performance of the proposed method in the presence of wind noise, we choose a small car as the test target, which emits less sounds and is more vulnerable to wind noise than a heavy or medium duty car. Observing the two original acoustic signals that are shown in Figures 11a and 12a, it is clear that the acoustic signal of the car is smaller. Figure 12b shows how the sub-band correlation of a small car changes during passing the UCA at level 4 wind power. An ultrasonic anemometer was placed near the microphone array to record the wind speed in real time. We record the moment when the target is passing the UCA and analyze the corresponding wind speed data. Figure 12c shows the variations of wind speed during the target passing the UCA. The performance to suppress wind noise of the aforesaid four DOA estimation methods is shown in Figure 13. As shown in Figure 12, the wind speed in the temporal interval 1 and the temporal interval 2 accelerates and drastically changes. As a result, the correlations between the sub-band in the low frequency parts are severely deteriorated, and the influence of each sub-band is quite different. The SMSC in this region varies randomly, from 0 to approximately 1, in comparison with Figure 11b where the target passing by, the SMSC of all the sub-bands in this temporal frequency interval are approximately equal to 1. Since the narrowband MUSIC method [19] chooses the fixed frequency within the bandwidth to calculate manifold matrix, the DOA estimation results will deviate significantly when the area near the fixed frequency is heavily contaminated by wind noise, which can be proved by the few large deviation in the Temporal interval 1 and Temporal interval 2 in Figure 13a.The sub-band method [20] selects the sub-band that has the largest SMSC for DOA estimation. However, the problem is that when the divided sub-band is too narrow, the result will have a large deviation, and when the divided sub-band is too wide, it cannot effectively suppress the turbulent wind noise interference in the discontinuous area of the sub-band. Thus, several significant transitions occur in Temporal interval 1 and Temporal interval 2 in Figure 13b. Figure 13c shows that when the target passing the CPA point, the curve of TCT method is the smoothest among all of the curves. Although there are some fluctuations in Temporal interval 1 and Temporal interval 2, it is obviously better than the narrowband MUSIC [19] and sub-band method [20]. Figure 13d shows the entire tracking process of the designed weighted sub-band method when the target passes through the microphone array. Even in the two intervals with severe wind noise disturbance, the curve is still very smooth and the only major fluctuation of the DOA occurs around the 16th seconds. From Figure 12b, we can find that the number of sub-bands of which the threshold is bigger than 0.5 in the bandwidth range is small, the limited number of sub-bands involved in the weighting process causes this problem. Figure 13e is actual DOA change curve of our experiment. Take Figure 13e as the base results, the distributions of the DOA errors for four methods are shown in Figure 13f. The absolute mean and standard variance of the error are also given in Table 1. As shown in these results, the weighted sub-band method provide the most robust and accurate (smallest standard variance) DOA estimates.

Conclusions
In this paper, we propose a robust DOA estimation method that is designed for the sequential movement events of vehicles. We carefully analyze the correlation characteristics of the vehicle signal, respectively, with and without wind disturbance. Therefore, by selecting the sub-bands with high SMSC, we obtain the sub-bands with high SNR. Referring to the idea of ISM, we calculate the MUSIC spectrum using these sub-bands to enhance the robustness of DOA estimation. In comparison with the narrowband MUSIC method, the sub-band method, and classical TCT method, the weighted sub-band method is the best one for the DOA estimation with the presence of wind noise. By selecting the proper number of sub-bands to participate in the weighted calculation, we can effectively control the amount of computation. Therefore, the proposed method can run on a resource-limited embedded processor, which has been verified in the designed small aperture microphone array system. Field experiments have also been conducted to further test the actual achievable performance for the proposed method. The results demonstrate that our proposed method has the best robustness when compared with other three methods in terms of the absolute mean error and standard variance.