Next Article in Journal
Effects of Nanocylinders on the Whispering Gallery Modes in a Microcylinder
Next Article in Special Issue
A Highly Reliable and Cost-Efficient Multi-Sensor System for Land Vehicle Positioning
Previous Article in Journal
A High Performance Torque Sensor for Milling Based on a Piezoresistive MEMS Strain Gauge
Previous Article in Special Issue
Vision-Based Steering Control, Speed Assistance and Localization for Inner-City Vehicles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vehicle Detection Based on Probability Hypothesis Density Filter

Robotics and Embedded Systems, Technische Universität München, 80333 München, Germany
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(4), 510; https://doi.org/10.3390/s16040510
Submission received: 7 January 2016 / Revised: 29 March 2016 / Accepted: 31 March 2016 / Published: 9 April 2016
(This article belongs to the Special Issue Sensors for Autonomous Road Vehicles)

Abstract

:
In the past decade, the developments of vehicle detection have been significantly improved. By utilizing cameras, vehicles can be detected in the Regions of Interest (ROI) in complex environments. However, vision techniques often suffer from false positives and limited field of view. In this paper, a LiDAR based vehicle detection approach is proposed by using the Probability Hypothesis Density (PHD) filter. The proposed approach consists of two phases: the hypothesis generation phase to detect potential objects and the hypothesis verification phase to classify objects. The performance of the proposed approach is evaluated in complex scenarios, compared with the state-of-the-art.

1. Introduction

Traffic accidents are a major cause of death worldwide. A study by the World Health Organization (WHO) reports that an estimated 1.2 million people die in traffic accidents every year, and up to 50 million people are injured [1]. Autonomous driving thus becomes significantly important in order to prevent accidents in traffic scenarios. However, it is still quite challenging for autonomous driving in all scenarios. For automotive manufacturers, the technology behind autonomous driving has been continually refined as a long-term goal, whereas the Advanced Driver Assistance System (ADAS) has been proposed as a short-term development to gradually improve road safety. Numerous ADAS functions have been developed to help drivers avoid accidents, improve driving efficiency, and reduce driver fatigue, in which vehicle detection plays an important role.
Most approaches rely on vision techniques to first detect Regions Of Interest (ROI) and then classify vehicles [2]. Khammari et al. use the Adaboost classification to detect vehicles [3]. Miller et al. and Paragios et al. have also utilized filtering techniques to detect vehicles [4,5]. Meanwhile, vehicle profile symmetry and the corresponding shadows are used in Reference [6]. However, vision techniques suffer from light intensities.
LiDAR is also widely used in vehicle detection. In contrast to vision sensors, LiDAR is robust against light intensities and offers a range of information [7,8,9,10,11]. Teichman et al. use the log odds estimators to recognize objects, where the performance is demonstrated in a large scale environment. Dominguez et al. demonstrate a data fusion platform for tracking vehicles [12]. Compared with vision sensors, LiDAR measurement often suffers from data association issues.
This paper extends our previous work to detect vehicles by using information from LiDAR, where objects are represented by the position and shape parameters (more details would be explained later). In Reference [13], we used the Difference of Normal (DoN) operator and the Random Hypersurface Model (RHM) to cluster the points cloud data and estimate the shape parameters [14,15]. To avoid the data association issue, the Probability Hypothesis Density (PHD) filter is proposed to detect vehicles based on Random Finite Set statistics (RFSs) [16]. In RFSs, several approaches are developed to avoid the data association issue, including the PHD filter, the Cardinalized PHD (CPHD) filter [17] and the Bernoulli filter [18]. The PHD filter propagates the probability hypothesis density function over the single target state space, whereas the CPHD filter also propagates the distribution of the target numbers (cardinality). By using the CPHD filter, the system requires more complex implementations and achieves more reliability in cardinality estimation. As the main goal of this paper is to estimate the states in a speed-critical environment, the PHD filter is thus considered. Unlike the PHD or CPHD filter, the multi-Bernoulli filter propagates the posterior target density. Although it has the same complexity as the PHD filter, the performance is better in highly nonlinear environments (it does not require the additional clustering step for state estimation) [18]. As the proposed RHM could be linearly implemented, the PHD filter is considered as the cheapest solution. The estimated states are then classified by the Support Vector Machine (SVM) to eliminate the non-vehicle objects.
The contributions are summarized as: first and foremost, the proposed solution achieves high performance in the presence of unknown data association environments. Furthermore, the shape parameter is first proposed to classify objects.
This paper is structured as follows: Section 2 describes the random hypersurface model, as well as the probability hypothesis density filter for hypothesis generation. Section 3 introduces the support vector machine for hypothesis verification. Section 4 demonstrates the performance of the experiments in urban environments. Finally, the paper is concluded in Section 5.

2. Hypothesis Generation

In the hypothesis generation phase, LiDAR measurements are filtered based on the Random Hypersurface Model (RHM). To further extend RHM in the presence of unknown data association scenarios, the Gaussian Mixture Probability Hypothesis Density (GMPHD) filter is proposed. Notice that objects are tracked and estimated in the 2D Cartesian coordinate system, whereas the depth information is unnecessarily required.
The result of the generation phase is then utilized for the verification phase to eliminate non-vehicles.

2.1. Random Hypersurface Model (RHM)

As illustrated in Figure 1, in a 2D Cartesian coordinate system, a point is considered as a scaled point with the factor within the range [ 0 , 1 ] drawn on the surface. Thus, the RHM is defined as:
S ( b ¯ k ) denotes the surface which consists of both the shape parameter b ¯ k and the center c k . Each point that lies on the surface is represented as a scaled boundary when s is drawn from [ 0 , 1 ] :
c k + s · S ( b ¯ k ) c k
In Figure 2, r ( ϕ ) denotes the distance function calculated from the center to boundary on angle φ in the polar coordinate system.
Assuming r ( ϕ ) consists of the shape parameter b ¯ k and the center c k , the surface is represented as
S ( b ¯ k ) = { s · r ( b ¯ k , ϕ ) · e ( ϕ ) + c k | ϕ [ 0 , 2 π ] , s [ 0 , 1 ] }
where e ( ϕ ) : = cos ϕ sin ϕ and r ( b ¯ k , ϕ ) denote the unit vector and the radial function in the form of the Fourier series, respectively. Due to its periodic proprieties, the Fourier series expansion of degree N F becomes
r ( b ¯ k , ϕ ) = a k 0 + j = 1 N F { a k j cos ( j ϕ ) + b k j sin ( j ϕ ) }
where b ¯ k is given by
b ¯ k = [ a k 0 , a k 1 , b k 1 , , a k N F , b k N F ] T
If φ is fixed, Equation (3) is represented as
r ( b ¯ k , ϕ ) = R ( ϕ ) · b ¯ k
where
R ( ϕ ) = [ 1 , cos ( ϕ ) , sin ( ϕ ) , , cos ( N F ϕ ) , sin ( N F ϕ ) ]
Notice that a low number of Fourier coefficients encode rough information of the surface, whereas a larger number of coefficients give more details.

Bayes Filter

The RHM represents the shape information by using the Fourier coefficients, and the Bayes filter is utilized to calculate the corresponding parameters.
  • Process model
Assuming the state x k denotes the Fourier descriptors b ¯ k and does not drift against time, the process model is described as
x k = A k x k 1 + w k
where A k and w k denote the identity matrix and the process noise, respectively.
  • Measurement model
As illustrated in Figure 1, a single measurement y k from LiDAR is originated from a surface boundary point z k with scaled factor s,
y k = f ( z k , s ) + v k , f ( z k , s ) S ( b ¯ k )
where v k denotes the measurement noise.
Using Equation (2), Equation (8) becomes
y k = f ( z k , s ) + v k = s · r ( b ¯ k , ϕ ) · e ( ϕ ) + c k + v k : = h ( x k , v k )
which maps the relationship between the state b ¯ k and the measurement y k .
Based on Equation (5), the measurement model is represented as
y k = s · R ( ϕ ) · b ¯ k · e ( ϕ ) + c k + v k
with algebraic manipulations on Equation (10), we get
| | y k c k | | 2 = s 2 · | | R ( ϕ ) · b ¯ k | | 2 + 2 s R ( ϕ ) b ¯ k e ( ϕ ) T v k + | | v k | | 2
The measurement model is thus acquired as:
0 = s 2 · | | R ( ϕ ) · b ¯ k | | 2 + 2 s R ( ϕ ) b ¯ k e ( ϕ ) T v k + | | v k | | 2 | | y k c k | | 2
where a pseudo measurement 0 is used to model the relationship between the state, the scaled factor, the measurement and its noise. The state x k convergences to the true Fourier descriptor b ¯ k by updating with a large number of measurements.
Notice that the proposed measurement model is implemented in the 2D Cartesian coordinate system and is only effective on the backside of the target. In addition, the depth information from LiDAR measurement is unrequested.
Figure 3 exhibits the estimated result with a cross target, in which the state x and measurement y are represented by 13 Fourier coefficients and 2D Cartesian coordinates, respectively.

2.2. Probability Hypothesis Density (PHD) Filter

As illustrated in Figure 3, the RHM estimates the shape information by using the Bayes filter. The key challenge for practice implementation is data association. A traditional estimator operates in measurement-to-target known scenarios, where all assignments are confirmed. However, LiDAR provides a large number of measurements without any association information. In previous work, the Difference of Normal (DoN) operator is utilized as a preprocessing procedure to cluster measurements. The miss detection is quite high since the clustering process may also eliminate objects. Hence, the PHD filter is proposed to track objects in presence of unknown data association scenarios.

2.2.1. Overview

The PHD filter is represented with the set-valued state and observation for multiple-object Bayesian filtering. All targets and observations are collected and represented in the set space, whereas the data association problem in traditional filtering domain is avoided. Figure 4 is a basic introduction of the RFS statistic. Compared to the single target filtering, the PHD filter relies on the random finite set statistics to process data in the set space level.

2.2.2. Mathematic Background

Considering the survived targets S k | k 1 , the spontaneous targets σ k and the spawned targets B k | k 1 , the set-valued state is described as:
X k = [ ζ X k 1 S k | k 1 ( ζ ) ] [ ζ X k 1 B k | k 1 ( ζ ) ] σ k
In addition, the set observation Z k consists of the reflections from both the targets θ k ( x ) and the clutters κ k .
Z k = [ x X k θ k ( x ) ] κ k
Similar to Bayesian estimator, the PHD filter is also divided into prediction and update processes. Notice that D and f k | k 1 ( x | ζ ) denote the posterior density and the transition function, respectively. ζ is the previous state.
Thus, the prediction is represented as:
D k | k 1 ( x ) = [ [ P S ( ζ ) f k | k 1 ( x | ζ ) + β ( x | ζ ) ] D k 1 ( ζ ) d ζ ] + γ k
The intensity function is updated based on the measurement set Z k :
D k ( x ) = ( 1 P D ) D k | k 1 ( x ) + z Z k P D g k ( z i | x ) D k | k 1 ( x ) κ ( z i ) + P D g k ( z i | ζ ) D k | k 1 ( ζ ) d ζ
where g k ( z i | x ) and P D denote the likelihood function and the detection probability, respectively.
Notice that the predict function in Equation (15) are affected by targets, which enter the scene ( γ k ), and survive from the previous time step P S and the spawn targets ( β ( x | ζ ) ).
The update function in Equation (16) corrects the prediction by using innovations from the observations. Notice that the clutter rate is also considered in the update function.
N ( k ) = Ψ D k | k ( x ) d x
Equation (17) exhibits that the integration of the intensity function represents the number of targets. Meanwhile, the intensity is not a probability density and thus unnecessarily sums up to 1 [19].
The PHD recursions have multiple integrals with no closed form representation. Thus, the most common approach is to use Gaussian Mixture (GM)-PHD approximations [20]:
f k | k 1 ( x | ζ ) = N ( x ; F k 1 ζ , Q k 1 )
g k ( z | x ) = N ( z ; H k x , R k )
where z and x denote the current measurement and state, whereas ζ denotes the previous state. A Gaussian distribution is represented as N ( · ; m , P ) with mean m and covariance P. F k 1 and Q k 1 denote the transition matrix and process covariance, respectively. H k and R k denote the observation matrix and observation noise covariance. Notice that the detection and survival probabilities are constant values:
P S , k ( x ) = P S , P D , k ( x ) = P D
Birth targets γ k are modeled as:
γ k ( x ) = i = 1 J γ , k ω γ , k ( i ) N ( x ; m γ , k ( i ) , P γ , k ( i ) )
where ω γ , k ( i ) , P γ , k ( i ) , m γ , k ( i ) and J γ , k denote the weight, covariance, mean and amount of the Gaussians.
Assuming the posterior intensity at time k 1 is a Gaussian mixture:
D k 1 ( x ) = i = 1 J k 1 ω k 1 ( i ) N ( x ; m k 1 ( i ) , P k 1 ( i ) )
Equation (15) to time k is also a Gaussian mixture
D k | k 1 ( x ) = P S i = 1 J k 1 ω k 1 ( i ) N ( x ; m S , k | k 1 ( i ) , P S , k | k 1 ( i ) ) + γ k ( x )
m S , k | k 1 ( i ) = F k 1 m k 1 ( i ) , P S , k | k 1 ( i ) = Q k 1 + F k 1 P k 1 ( i ) F k 1 T
and the Equation (16) at time k is calculated as
D k ( x ) = ( 1 P D ) D k | k 1 ( x ) + z Z k D D , k ( x ; z )
where
D D , k ( x ; z ) = j = 1 J k | k 1 ω k ( j ) ( z ) N ( x ; m k | k ( j ) ( z ) , P k | k ( j ) )
ω k j ( z ) = P D w k | k 1 ( j ) q k ( j ) ( z ) κ k ( z ) + P D l = 1 J k | k 1 w k | k 1 ( l ) q k ( l ) ( z )
q k ( j ) ( z ) = N ( z ; H k m k | k 1 ( j ) , H k P k | k 1 ( j ) H k T + R k )
m k ( j ) ( z ) = m k | k 1 ( j ) + K k ( j ) ( z H k m k | k 1 ( j ) )
P k ( j ) = [ I K k ( j ) H k ] P k | k 1 ( j )
K k ( j ) = P k | k 1 ( j ) H k T [ H k P k | k 1 ( j ) H k T + R k ] 1
The GMPHD filter addresses the data association challenge in contrast to the standard Bayes filter. The process model Equation (18) and measurement model Equation (19) are equal to the RHM Bayes filter in Equations (7) and (8), whereas the implementation process is different.

2.3. RHM–GM–PHD Filter

In Section 2.2, the GMPHD filter is introduced for dealing with unknown data association issues. Notice that the standard GMPHD filter operates in condition that one reflection is received for each target per frame, called “Point Target (PT) tracking”. In the vehicle detection scenario, a large number of measurements would be collected from the surface of a single object, called “Extended Target (ET) tracking”. Therefore, the GMPHD filter should be redesigned for dealing with extended targets solely relying on LiDAR measurement.
The process model and measurement model are similar to the RHM Bayes filter in Equations (7) and (8). The prediction equation of the ET–GM–PHD filter are also the same as the standard GMPHD filter. The measurement update formulas for the ET–GM–PHD filter is introduced as:
D k | k ( x ) = L Z k ( x ) D k | k 1 ( x | Z )
and the pseudo-likelihood function L Z k ( x ) is defined as
L Z k ( x ) = 1 ( 1 e γ ( x ) ) P D ( x ) + e γ ( x ) P D ( x ) p Z k w p W p γ ( x ) | W | d W z W ϕ z ( x ) λ k c k ( z )
where λ k c k ( z ) is the mean number of clutter measurements, c k ( z ) is the spatial distribution of the clutter, notation p Z k denotes that p partitions the measurement set Z k into non-empty cells W, notation W p denotes that W is a cell in the partition p, w p and d W denote the non-negative coefficients for each partition and cell, and ϕ ( x ) denotes the same likelihood function for a single measurement in Equation (12).
Here,
w p = W p d W p Z k W p d W
and
d W = δ | W | , 1 + D k | k 1 [ e γ γ | W | P D z W ϕ z λ k c k ( z ) ]
where δ i , j is the Kronecker delta function and | W | is the number of measurements in cell W. More details of the implementation process could be found in [21,22].
Hence, the ET–GM–PHD filter tracks the potential objects by solely relying on LiDAR without any association or cluster process. The estimated states represent the potential objects and would be filtered again by the support vector machine to eliminate non-vehicle objects.

3. Hypothesis Verification

To eliminate the outliers, the support vector machine is utilized to classify the vehicle and non-vehicle Fourier coefficients.

3.1. Support Vector Machine (SVM)

As exhibited in Figure 5, the SVM is proposed to obtain classifiers with good generalization [23]. The mathematical background is introduced as follows:
For x i R n with respect to the classification y i { 1 , + 1 } , i = 1 , , k , the hyperplane is defined as:
f ( x ) = w T x + b
to linearly separate each data. Notice that x , w and b denote the input vector, weight vector and the bias, respectively. Hence a maximum margin is found to separate positive class from negative class based on
f ( x ) = sgn ( w T x + b )
The calculation of the hyperplane is subjected to the following constraints:
Minimize 1 2 | | w | | 2
Subject to y i ( w i x i + b ) 0
and the classification performance relies on the optimization as:
Maximize i = 1 k α i 1 2 i , j = 0 k α i α j y i y j K ( x i , x j )
Subject to i = 1 k y i α i = 0 , α i 0 f o r i = 1 , , k
where K ( x i , x j ) and α i denote the Kernel function and Lagrange multiplier, respectively. Notice that if the data can not be separated linearly, the kernel function changes according to
K ( x i , x j ) = K ( x i , x j ) + 1 C δ i j
where C denotes the penalize parameter.

3.2. Implementation Detail

During the implementation, there are still issues in both the generation and verification phases.
  • ET–GM–PHD Implementation
To track objects in multiple frames, the alignment issue should also be considered. Thus the state to be tracked in the ET–GM–PHD filter is x k = [ v x , v y , x , y , ( b ¯ k ) T ] T , where [ v x , v y , x , y , ] T is the normal state of the object centroid in the 2D Cartesian coordinate system and ( b ¯ k ) is the shape parameter by using 13 Fourier descriptors. Although a higher number of Fourier descriptors returns more details of the surface, the computation performance and robustness is unsatisfied. In this paper, 13 Fourier descriptors are selected based on the experience during the experiment.
Each object follows the liner Gaussian dynamic Equation (18) with the following configurations:
F k = diag ( A k , I 13 ) , A k = I 2 0 2 I 2 I 2
Q k = diag ( C k , 0 . 03 I 13 ) , C k = I 2 0 2 2 0 2 2 1 4
where I n and 0 n denote the n × n identity and zero matrix, respectively. The measurement covariance is given using parameters diag [ 0 . 5 , , 0 . 5 ] .
Regarding the measurement model, it follows the non-linear Gaussian dynamic Equation (19) and is described as:
0 = s 2 · | | R ( ϕ ) · b ¯ k | | 2 + 2 s R ( ϕ ) b ¯ k e ( ϕ ) T v k + | | v k | | 2 | | y k m k | | 2
where s is a Gaussian function with mean 0.5 and variance 0.03. The birth intensity of the PHD filter is given according to the geometry information. To reduce the computation complexity, measurements which fall into the road regions are utilized to initialize the vehicles. The birth intensity is thus given by
γ k ( x ) = 0 . 1 N ( x ; m γ , k 1 , P γ , k 1 ) + 0 . 1 N ( x ; m γ , k 2 , P γ , k 2 ) + 0 . 1 N ( x ; m γ , k 3 , P γ , k 3 )
where m γ , k 1 = [ 0 , 0 , 0 , 1 , 2 , 0 , , 0 ] , m γ , k 2 = [ 0 , 0 , 5 , 1 , 2 , 0 , , 0 ] and m γ , k 3 = [ 0 , 0 , 5 , 1 , 2 , 0 , , 0 ] represents the vehicles in the middle (0,1), left (−5,1) and right (5,1), respectively. P γ , k 1 = P γ , k 2 = P γ , k 3 = diag [ 2 , 2 , 2 , 2 , 0 . 3 , , 0 . 3 ] .
  • SVM Implementation
The KITTI dataset is utilized for training the classifier, which provides a set of 5000 training frames with 1893 manually labeled objects (car, van, tram, misc, pedestrian, cyclist, trunk and so on). For each object, the original measurements (all objects are labeled with a 2D box, where measurement-to-target association is confirmed) are projected to the 2D Cartesian coordinate system around its original point. The 13 Fourier coefficients are then calculated by using the standard nonlinear Kalman filter [15]. Instead of multiple categories, objects are only considered as vehicles and non-vehicles (actually it is divided by cars and non-cars). Then, the SVM is trained by collecting Fourier coefficients from all calculated objects. In further evaluation, another 2000 frames of test data is utilized to guarantee the performance both quantitatively and qualitatively. The SVM is implemented in both the training and test phases without cross-validation.
  • Key Parameters and Open Issues
During the ET–GM–PHD process, potential objects are collected in the set-valued state. The Fourier coefficients describe the shape information, and the position vector represents the location. The estimated objects are also shown by using bounding boxes, where the mean point uses the position and the width/height are calculated based on the Fourier coefficients (the Fourier coefficients represent the rough shapes in polar coordinate systems, in which the width and height of the boundary box are calculated by setting the φ equals to 0 and π 2 ). Since the PHD filter addresses the data association issue, the points cloud data is directly utilized to estimate the set-valued states. No further cluster process is required. Meanwhile, for each single measurement, the probabilities of both detection and survival are constant and no more than 1. In addition, the PHD filter may estimate close objects as single objects mainly due to the scale factor s. For the RHM model, measurements are considered as random draws from boundaries with different scaling factors in the range [ 0 , 1 ] . When objects are close to each other, it is quite challenging to distinguish them.
During the SVM process, the calculated coefficients from objects are utilized to eliminate the outliers. Since most cars have a similar width/length rate, the car-labeled objects are treated as vehicles and the rest are non-vehicles. Furthermore, the 13 Fourier coefficients have been found to be quite challenging for linearly separating vehicles and non-vehicles. To better train the SVM, the radial based kernel function is utilized as:
K ( x i , x j ) = exp ( | | x i x j | | 2 2 σ 2 )
where σ affects distributing complexity in the feature space.

4. Experiment Evaluation

To evaluate the approach quantitatively and qualitatively, the KIT dataset is utilized and compared with the state-of-the-art [24]. During the experiment, the proposed approach is implemented in Matlab with 2 Cores@3 GHz, and the average time is 5 s per frame.
Figure 6, Figure 7 and Figure 8 demonstrate the detection performance based on the proposed approach in one scenario. In Figure 6, a bicycle is fully observed in the middle of the road. On the left side, a parking car is observed with partial occlusion. On the right side, both car and pedestrian are fully observed. As illustrated in Figure 7, potential objects are detected based on the RHM–ET–PHD filter. It is observed that the proposed approach extracts the potential objects based on the geometry of the road, where the birth model plays an important role. The extracted objects are drawn by boundary boxes calculated by the set-valued state (the center point is based on the position, and the width/height are based on the Fourier coefficients). Afterwards, the SVM is utilized to eliminate the non-vehicle objects. Figure 8 shows the results of the verification phase. Due to the occluded issues, it is observed that the left vehicle is also eliminated.
Table 1 demonstrates the overall performance, in contrast to the state-of-the-art, in all scenarios from the KITTI dataset. Although there are also approaches using cameras, the evaluation focuses on the algorithms which only use LiDAR measurements.
As illustrated in Table 1, moderate, easy and hard denote the occlusion level of vehicles, with respect to partly occluded, fully visible and difficult to see. Among the references, the proposed approach achieves a high performance for the easy category and poor performance for both moderate and hard categories. In easy scenarios, all vehicles are fully observed and the corresponding reflections on surfaces are uniformly distributed. Hence, the proposed approach can track and estimate the states successfully, and the final classification process has high performance. In moderate and hard scenarios, the corresponding performance drops significantly in both the generation and verification processes. For PHD filter, although it detects the potential objects, the calculated Fourier coefficients are strongly influenced by the invisible measurements. For SVM classification, the training process mainly relies on the visible measurements to calculate Fourier coefficients.
Nevertheless, compared with the overall performance in easy scenarios, the proposed approach still improves almost 20 % 40 % .
As a summary, the contributions are concluded as follows: first and foremost, the proposed framework solely relies on LiDAR measurements for vehicle detection in the presence of unknown data association environments. Furthermore, the Fourier coefficient is first proposed for object classification and concluded with high performance for fully visible vehicles.

5. Conclusions

Vehicle detection is important for developing driver assistance systems. To address the data association problem that suffers from points cloud data, the Probability Hypothesis Density (PHD) filter is proposed in this paper. The proposed scheme utilizes contour information for classification. The evaluation results illustrate a high performance in contrast to the state-of-the-art techniques.
Future work focuses on the improvement of detecting occluded vehicles.

Acknowledgments

This work was supported by the German Research Foundation (DFG) and the Technische Universität München within the funding program Open Access Publishing.

Author Contributions

All authors have contributed to the manuscript. Feihu Zhang performed the experiments and analyzed the data. Alois Knoll provided helpful feedback on the experiments design.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Global Status Report on Road Safety 2013. Available online: http://www.who.int/violence_injury_prevention/road_safety_status/2013/en/ (accessed on 1 October 2013).
  2. Leon, L.; Hirata, R. Vehicle detection using mixture of deformable parts models: Static and dynamic camera. In Proceedings of the 25th SIBGRAPI Conference on Graphics, Patterns and Images (SIBGRAPI), Ouro Preto, Brazil, 22–25 August 2012; pp. 237–244.
  3. Khammari, A.; Nashashibi, F.; Abramson, Y.; Laurgeau, C. Vehicle detection combining gradient analysis and AdaBoost classification. In Proceedings of the 2005 IEEE Intelligent Transportation Systems, Vienna, Austria, 13–15 September 2005; pp. 66–71.
  4. Sun, Z.; Bebis, G.; Miller, R. On-road vehicle detection using evolutionary Gabor filter optimization. IEEE Trans. Intell. Transp. Sys. 2005, 6, 125–137. [Google Scholar] [CrossRef]
  5. Paragios, N.; Deriche, R. Geodesic active contours and level sets for the detection and tracking of moving objects. IEEE Trans. Pattern Anal. Mach. Intell. 2000, 22, 266–280. [Google Scholar] [CrossRef]
  6. Sun, Z.; Bebis, G.; Miller, R. Monocular precrash vehicle detection: Features and classifiers. IEEE Trans. Image Process. 2006, 15, 2019–2034. [Google Scholar] [PubMed]
  7. Huang, L.; Barth, M. Tightly-coupled LIDAR and computer vision integration for vehicle detection. In Proceedings of the 2009 IEEE Intelligent Vehicles Symposium, Xi’an, China, 3–5 June 2009; pp. 604–609.
  8. Behley, J.; Kersting, K.; Schulz, D.; Steinhage, V.; Cremers, A. Learning to hash logistic regression for fast 3D scan point classification. In Proceedings of the 2010 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Taipei, Taiwan, 18–22 October 2010; pp. 5960–5965.
  9. Li, Y.; Ruichek, Y.; Cappelle, C. Extrinsic calibration between a stereoscopic system and a LIDAR with sensor noise models. In Proceedings of the 2012 IEEE Conference on Multisensor Fusion and Integration for Intelligent Systems (MFI), Hamburg, Germany, 13–15 September 2012; pp. 484–489.
  10. Xiong, X.; Munoz, D.; Bagnell, J.; Hebert, M. 3-D scene analysis via sequenced predictions over points and regions. In Proceedings of the 2011 IEEE International Conference on Robotics and Automation (ICRA), Shanghai, China, 9–13 May 2011; pp. 2609–2616.
  11. Li, Y.; Ruichek, Y.; Cappelle, C. Optimal Extrinsic Calibration Between a Stereoscopic System and a LIDAR. IEEE Trans. Instrum. Meas. 2013, 62, 2258–2269. [Google Scholar] [CrossRef]
  12. Dominguez, R.; Onieva, E.; Alonso, J.; Villagra, J.; Gonzalez, C. LIDAR based perception solution for autonomous vehicles. In Proceedings of the 11th International Conference on Intelligent Systems Design and Applications (ISDA), Cordoba, Spain, 22–24 November 2011; pp. 790–795.
  13. Zhang, F.; Clarke, D.; Knoll, A. LiDAR based vehicle detection in urban environment. In Proceedings of the 2014 International Conference on Multisensor Fusion and Information Integration for Intelligent Systems (MFI), Beijing, China, 28–29 September 2014; pp. 1–5.
  14. Ioannou, Y.; Taati, B.; Harrap, R.; Greenspan, M. Difference of normals as a multi-scale operator in unorganized point clouds. In Proceedings of the 2012 Second International Conference on 3D Imaging, Modeling, Processing, Visualization and Transmission (3DIMPVT), Zurich, Switzerland, 13–15 October 2012; pp. 501–508.
  15. Baum, M.; Hanebeck, U. Shape tracking of extended objects and group targets with star-convex RHMs. In Proceedings of the 14th International Conference on Information Fusion (FUSION), Chicago, IL, USA, 5–8 July 2011; pp. 1–8.
  16. Mahler, R. Multitarget Bayes filtering via first-order multitarget moments. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1152–1178. [Google Scholar] [CrossRef]
  17. Mahler, R. Approximate multisensor CPHD and PHD filters. In Proceedings of the 13th Conference on Information Fusion (FUSION), Edinburgh, UK, 26–29 July 2010; pp. 1–8.
  18. Vo, B.T.; Vo, B.N.; Hoseinnezhad, R.; Mahler, R. Robust multi-Bernoulli filtering. IEEE J. Sel. Top. Signal Process. 2013, 7, 399–409. [Google Scholar] [CrossRef]
  19. Mahler, R. Multitarget Bayes filtering via first-order multitarget moments. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1152–1178. [Google Scholar] [CrossRef]
  20. Vo, B.N.; Ma, W.K. The Gaussian Mixture Probability Hypothesis Density Filter. IEEE Trans. Signal Process. 2006, 54, 4091–4104. [Google Scholar] [CrossRef]
  21. Mahler, R. PHD filters for nonstandard targets, I: Extended targets. In Proceedings of the 12th International Conference on Information Fusion, FUSION ’09, Seattle, WA, USA, 6–9 July 2009; pp. 915–921.
  22. Han, Y.; Zhu, H.; Han, C. A Gaussian-mixture PHD filter based on random hypersurface model for multiple extended targets. In Proceedings of the 16th International Conference on Information Fusion (FUSION), Istanbul, Turkey, 9–12 July 2013; pp. 1752–1759.
  23. Cortes, C.; Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  24. Geiger, A.; Lenz, P.; Urtasun, R. Are we ready for Autonomous Driving? The KITTI Vision Benchmark Suite. In Proceedings of the 2012 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Providence, RI, USA, 16–21 June 2012; pp. 3354–3361.
  25. Wang, D.Z.; Posner, I. Voting for voting in online point cloud object detection. In Proceedings of the Robotics: Science and Systems, Rome, Italy, 13–17 July 2015.
  26. Plotkin, L. PyDriver: Entwicklung Eines Frameworks für Räumliche Detektion und Klassifikation von Objekten in Fahrzeugumgebung. Bachelor’s Thesis, Karlsruhe Institute of Technology, Karlsruhe, Germany, March 2015. [Google Scholar]
  27. Behley, J.; Steinhage, V.; Cremers, A. Laser-based segment classification using a mixture of bag-of-words. In Proceedings of the 2013 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Tokyo, Japan, 3–7 November 2013; pp. 4195–4200.
Figure 1. Random hypersurface model.
Figure 1. Random hypersurface model.
Sensors 16 00510 g001
Figure 2. Example of star-convex object.
Figure 2. Example of star-convex object.
Sensors 16 00510 g002
Figure 3. Bayesian estimation in a data association known environment.
Figure 3. Bayesian estimation in a data association known environment.
Sensors 16 00510 g003
Figure 4. Set-valued states and set-valued observations.
Figure 4. Set-valued states and set-valued observations.
Sensors 16 00510 g004
Figure 5. The concept of support vector machine (SVM)
Figure 5. The concept of support vector machine (SVM)
Sensors 16 00510 g005
Figure 6. Original data from image view.
Figure 6. Original data from image view.
Sensors 16 00510 g006
Figure 7. Result on hypothesis generation phase.
Figure 7. Result on hypothesis generation phase.
Sensors 16 00510 g007
Figure 8. Result on hypothesis verification phase.
Figure 8. Result on hypothesis verification phase.
Sensors 16 00510 g008
Table 1. Performance of the vehicle detection approach compared with the state-of-the-art.
Table 1. Performance of the vehicle detection approach compared with the state-of-the-art.
MethodOursVote3D [25]CSoR [26]mBoW [27]
Easy75%57%35%36%
Moderate15%48%26%24%
Hard3%43%23%18%
Average time5 s0.5 s3.5 s10 s

Share and Cite

MDPI and ACS Style

Zhang, F.; Knoll, A. Vehicle Detection Based on Probability Hypothesis Density Filter. Sensors 2016, 16, 510. https://doi.org/10.3390/s16040510

AMA Style

Zhang F, Knoll A. Vehicle Detection Based on Probability Hypothesis Density Filter. Sensors. 2016; 16(4):510. https://doi.org/10.3390/s16040510

Chicago/Turabian Style

Zhang, Feihu, and Alois Knoll. 2016. "Vehicle Detection Based on Probability Hypothesis Density Filter" Sensors 16, no. 4: 510. https://doi.org/10.3390/s16040510

APA Style

Zhang, F., & Knoll, A. (2016). Vehicle Detection Based on Probability Hypothesis Density Filter. Sensors, 16(4), 510. https://doi.org/10.3390/s16040510

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

Article Metrics

Back to TopTop