Research Article  Open Access
Precise Vehicle Position and Heading Estimation Using a Binary Road Marking Map
Abstract
Road markings are always present on roads to guide and control traffic. Therefore, they can be used at any time for vehicle localization. Moreover, they can be easily extracted by using light detection and ranging (LIDAR) intensity because they are brightly colored. We propose a vehicle localization method using a 2D road marking grid map. The grid map inserts the map information into the grid directly. Thus, an additional process (such as line detection) is not required and there is no problem due to false detection. We obtained road marking using a 3D LIDAR (Velodyne HDL32E) and binarized this information to store in the map. Thus, we could reduce the map size significantly. In the previous research, the road marking grid map was used only for position estimation. However, we propose a positionandheading estimation algorithm using the binary road marking grid map. Accordingly, we derive more precise position estimation results. Moreover, position reliability is an important factor for vehicle localization. Autonomous vehicles may cause accidents if they cannot maintain their lane momentarily. Therefore, we propose an algorithm for evaluating map matching results. Consequently, we can use only reliable matching results and increase position reliability. The experiment was conducted in Gangnam, Seoul, where GPS error occurs largely. In the experimental results, the lateral root mean square (RMS) error was 0.05 m and longitudinal RMS error was 0.08 m. Further, we obtained a position error of less than 50 cm in both lateral and longitudinal directions with a 99% confidence level.
1. Introduction
Precise vehicle localization is essential for driving safety in autonomous vehicles. Generally, a submeter level of localization performance is required for the autonomous vehicle to maintain its own lane, and it should be estimated with a high reliability of 95% [1] or more. Usually, Global Positioning System (GPS) is used for vehicle localization. However, GPS does not provide the position accuracy required for autonomous vehicles. Especially in urban areas, even if an expensive GPS (RTK)/INS is used, a position error of more than 1 m can occur. To solve this problem, many studies that improve position accuracy by combining various sensors (such as GPS, inertial measurement unit (IMU), vision, and light detection and ranging (LIDAR)) have been researched continuously.
Recently, many studies have been conducted to correct GPS error by generating maps based on LIDAR [2] because LIDARbased maps can be generated precisely and accurately. Li et al. [3] expressed the 3D environment on both sides of the road as a type of a depth map called road DNA. In addition, Im et al. [4] detected corners of vertical structures such as buildings and generated a cornerbased landmark map. However, the above methods have disadvantages in that they cannot be used when there are no vertical structures, such as a highway.
Therefore, we focused on the map using road markings. Road markings are always available because they exist everywhere the vehicle goes. Moreover, because of the use of bright colors, it is easy to detect them using the intensity of LIDAR, and there is a potential for robust localization performance because of the use of letters and symbols as well as lines such as lanes or stop lines. Therefore, road marking mapbased vehicle localization has been studied in this work. Levinson et al. [5, 6] and Wolcott and Eustice [7] performed vehicle localization on LIDAR and camera based on the intensity road marking grid map, respectively. Intensity road marking grid mapbased localization provides high accuracy. However, the intensity of LIDAR varies depending on the road surface, so the performance may deteriorate depending on the road surface condition of a given day. Further, the intensity grid map has the disadvantage of large map size. Schreiber et al. [8] used LIDAR to create a line map using road marking (such as lane, stop line, and curb). A stereo camera was used to estimate the vehicle position. A line map has the advantage of small map size. However, the estimation performance is relatively low and the road marking of letters and symbols is difficult to detect. Therefore, it cannot be used with such robust information.
In this paper, we proposed a method for the estimation of vehicle position and heading based on a binary grid map using road markings. The road marking grid map is a twodimensional (2D) map that stores the road marking information scanned by the LIDAR directly on the grid. Therefore, since the detection process is not necessary, complex letters and symbols can be stored and used for localization. However, the grid map has disadvantages that the size increases exponentially according to the grid size and map scale. Therefore, we generated the binary road marking information by binarizing the intensity of the LIDAR scan data and stored this information in the grid map. The binary road marking information has only 1bit information. Thus, we can reduce the map size significantly compared to the intensity road marking grid map. In addition, the problem of varying the intensity level according to the road surface condition is solved through intensity binarization. Most of the grid maps using road marking information are in the form of an occupancy grid map [9]. Therefore, each grid has a probability. However, the binary road marking grid map has only 1bit information. Thus, 0 indicates no road marking and 1 indicates road marking. Therefore, the map size becomes more compact. Moreover, since it contains only 1bit information instead of probability, map matching is possible by the simple correlation method.
Because road markings are on the road, the map matching performance deterioration in a situation such as traffic congestion is a disadvantage. Therefore, we used the accumulated scan data [10, 11] to generate a precise local map and match with the global map. Consequently, the problem of poor performance under traffic congestion was solved. In addition, the position and attitude for the precise local map were derived by accumulating incremental values of position and attitude derived from the iterative closest point (ICP) algorithm based on the geometric relationship between LIDAR scan data.
Vehicle localization using existing 2D gridbased road marking maps [12] only estimated the vehicle position. However, it is necessary to estimate not only the position but also the heading to ensure driving safety in autonomous vehicles. Therefore, we proposed a positionandheading estimation algorithm based on the binary road marking map. In the case of position estimation, fast and robust position estimation performance is derived using fast Fourier transform (FFT) based phase correlation [13, 14]. Further, the heading is derived by converting the grid map to the points and matching between the map points and the binary road marking points. The highly reliable position and heading matching results are derived using the algorithm for the evaluation of map matching result, and we used only these results for precise vehicle localization.
The contributions of this paper are as follows: (i)Proposed a method of the binary road marking grid map generation and vehicle position and heading estimation method using this map(ii)Derived reliable and accurate estimation results using the binary road marking grid map
We propose a new type of the road marking grid map using the binary information that is not a traditional road markingbased occupancy grid map. In addition, we propose the vehicle position and heading estimation algorithm using this map. Vehicle localization using the existing 2D gridbased road marking map only estimates the position. However, we propose an algorithm to estimate the heading in addition to the position. The RMS position error is approximately 0.10 m in both lateral and longitudinal directions in urban areas, and the estimation results are derived with a position error within 50 cm with 99% reliability. Accordingly, it is confirmed that very accurate and precise position and heading estimation results can be derived by using a binary road marking grid map.
The composition of this paper is as follows. In Section 2, we introduce a method for binarizing LIDAR intensity. Section 3 describes how to generate the binary road marking grid map, and Section 4 describes how to generate a precise local map by accumulating scan data. In Section 5, the map matching and vehicle localization algorithm are described. The experimental result is described in Section 6. Finally, the conclusion is given in Section 7.
2. Method of Road Marking Extraction Using Intensity Binarization
Road markings are the devices on paved roadways that provide guidance and information to drivers. They are made with bright paint so that the driver can distinguish them easily. Therefore, they can be easily extracted by using the intensity of LIDAR. Specifically, the threedimensional (3D) LIDAR [15] used in this paper provides 8bit intensity information. Thus, various values of intensity information can be obtained.
As shown in Figure 1, the 3D LIDAR produces 3D scan data by scanning. Figure 2 shows the intensity of the road scan data among the 3D scan data. A road marking on the road is distinguished from the asphalt.
In this paper, we binarized 8 bits of intensity to 1 bit. In the case of intensity, it is possible to express an image in gray scale, as shown in Figure 2, because it provides information with 8 bits. Thus, the intensity can be regarded as the brightness of the image, and binarization is possible by applying the threshold. For binarization, it is necessary to set the threshold value of the intensity. To solve this, Otsu’s thresholding method [16] is used. Otsu thresholding is a method to find the optimal threshold using the brightness distribution of all pixel values in the image. Using this method, we can find the threshold even if the intensity is changed in every scene.
However, even when scanning the same road marking, there is a difference in intensity among the vertical layers of the 3D LIDAR [17]. Figure 3 shows the intensity of the same road marking scanned by a 3D LIDAR for each vertical layer. The vertical layer means the number of vertically scanned 3D LIDAR channels. In this research, we used 32channel 3D LIDAR and 12 channels are used for scanning the road marking. As shown in this figure, intensity values are different depending on the vertical layer. Thus, if we binarize to the threshold of the entire layer, the road marking cannot be extracted properly. Therefore, we applied thresholding and binarization to each layer.
Figure 4 shows the result of binarizing intensity data of Figure 2. Lowintensity values, such as those for asphalt, are eliminated, and only the road marking is extracted. Thus, the road marking can be extracted from 3D LIDAR scan data.
3. Binary Road Marking Map Generation
3.1. Graph Optimization and Map Generation
A global map is essential for localization. In addition, a precise vehicle trajectory is required to generate a precise global map. In general, when using simultaneous localization and mapping (SLAM) [18], errors accumulate as the vehicle progresses. Therefore, when a global map is generated, the road marking does not match.
Figure 5 shows the binary road marking grid map generated by using the GPS (RTK)/INS (NovAtel RTK/SPAN system) trajectory. As shown in this figure, it can be confirmed that the road marking is not consistent. Therefore, we performed graph optimization using graph SLAM [19, 20]. Graph optimization is a method of setting constraints for each node and optimizing the trajectory using edge relationships between nodes. The incremental pose information outputted from the ICP scan matching algorithm was used as the edge measurement of the graph.
Figure 6 shows the result of graph optimization. As shown in this figure, the edge relation connected by black lines is set for all nodes and the optimization is performed using them. Nodes connected by black lines are the trajectories of different laps, and loop closure is performed by deriving edge relationships between these nodes. After that, the trajectory is optimized like a red dot. Figure 7 shows the binary road marking grid map generated by using the Graph SLAMoptimized trajectory. The binary road marking grid map was generated by accumulating the road markings extracted by the binarization process. Unlike the case in Figure 5, the road markings are consistent. In addition, the optimized trajectory is assumed to be the ground truth.
3.2. Advantage of Using a Binary Road Marking Map
The map used for vehicle localization using the difference in intensity is the intensity road marking grid map, and it is stored by using 8bit intensity values. Generally, the intensity road marking grid map is used in the form of a probability map [5, 6] that stores the mean and variance of the intensity. Further, due to the characteristic of the 2D grid map, when the grid size is smaller and the map scale is larger, the map size becomes larger. However, if the intensity is binarized and stored in a map, the size becomes very small as the information is compressed.
Table 1 gives a comparison of map size. The probability map stores the intensity mean and variance, and the binary map stores the binarized road marking information. As shown in Table 1, the map size of the binary map is approximately under 2% of the probability map. Unlike the probability map, in which the map size grows exponentially as the grid size becomes smaller, the map size of the binary map is very small and the amount of increase in size is not large as the resolution increases. In addition, the binary map consists of 0 s and 1 s. The characteristic of a road marking grid map is that there are many more 0 s in the map. Therefore, if we remove 0 and save it as a sparse matrix, the map size is further reduced and can be reduced to 0.4% of the probability map as shown in Table 1.

Table 2 shows an example of saving a binary road marking grid map. Only the grid position where the road marking is stored can be saved because the map format is 1 bit. Therefore, it is not necessary to save the map as a 2D full matrix format and it is possible to store a 2D road marking grid map efficiently.

When intensity is binarized, loss of information occurs. Therefore, we might think that in using the binary map for map matching, localization performance is poor compared with the intensity map. On the contrary, matching performance may deteriorate because intensity has more information. For example, the value of intensity changes according to the road surface condition [21]. Thus, when the road surface gets wet, the intensity value is larger than that for a dry road surface.
Figure 8 shows the intensity map of dry and wet road surfaces. Further, Figure 9 compares the intensity count and threshold of Figure 8. Because the total intensity amount of the two cases is different, the count is normalized. As shown in Figure 9, the intensity threshold is larger on the wet road surface. Therefore, this difference in intensity causes poor localization performance. On the other hand, if the intensity is binarized, the road marking can be extracted regardless of the road surface condition. Consequently, more robust localization performance can be obtained.
(a)
(b)
Figure 10 shows the binary map of dry and wet roads. It can be seen in this figure that there is almost no difference in the extraction results between the two maps. Thus, using a binary road marking grid map, we can use road marking information regardless of road surface condition. In addition, when using the binary road marking grid map, it is obvious whether it is road marking or not. Therefore, more clear matching results can be obtained in map matching.
(a)
(b)
4. Precise Local Map Generation
In this section, we describe how to generate a precise local map [22]. It is generated by using the result obtained by binarizing the road marking in Section 2. Generally, the interval between the vertical layers of 3D LIDAR is larger than horizontal scan interval. Therefore, the local map generated by only one scan data has insufficient longitudinal direction information. In addition, the road marking is located under the vehicle. Thus, the problem is that only little road marking information can be extracted in traffic congestion situations. To solve this problem, we generated a local map by accumulating extracted road marking data for a certain period and defined it as a precise local map.
Figure 11 shows the flow chart of precise local map generation. First, data is accumulated by combining the previous data with the current position, pose increment, and LIDAR scan data. Next, a precise local map update is conducted using the accumulated data. The accumulated data are used to generate the precise local map at the next time. In addition, the sliding window method [10] was applied to the data accumulated.
In the case of local map generation using data accumulation, position and pose errors are also accumulated. In particular, in the case of commercial GPS/DR (dead reckoning), considerable errors may accumulate even in a short time. Therefore, we used the position and pose increment derived from the ICP scan matching. The ICP scan matching is matched using a precise 3D LIDAR point cloud. Thus, we can obtain more precise position and heading increment values.
Figure 12 shows an example of position and pose increment derived using ICP. The precise local map is generated by accumulating the scan data at the position and pose increment obtained from each node. In this case, the scan data is the road marking information extracted from the binarization process of Section 2.
As shown in Figure 13, it can be seen that as the data accumulates, the road marking information generated in the local map becomes sufficient, especially longitudinal direction. Consequently, sufficient road marking information can be obtained for matching even in traffic congestion situations.
5. Map Matching and Vehicle Localization Algorithm
In this section, we describe the estimation of position and heading using the map matching algorithm. The map matching result is used as the final measurement value of the navigation filter only for good results by executing the evaluation of the matching result algorithm. The navigation filter uses the Kalman filter (KF). The specific vehicle localization algorithm is shown in Figure 14.
5.1. Map Matching Algorithm
In this part, we describe the map matching algorithm. is the map matching result between the global map and the local map. The global map is taken from predicted vehicle position using the time update, and the local map is generated at the current vehicle position. Therefore, we can obtain error correction values between the predicted vehicle position and current vehicle position from map matching. Further, the predicted vehicle position was derived by using the increment value of position and heading from ICP scan matching algorithm.
5.1.1. Position
FFTbased phase correlation [13, 14] is used to estimate the vehicle position. Correlation is a method mainly used for matching two data sets. In particular, phase correlation is a technique used mainly for image registration. Because phase information is separated and correlated in the frequency domain, phase correlation has robustness to noise. Further, unlike the serialsearchbased correlation, which sets the search area and performs individual correlation in the area, the FFTbased correlation can match the large image at high speed because of correlated images simultaneously. The frequency domain transformation of the global and local maps for phase correlation is as follows: where and are the global map and local map, respectively. and are transformations in the frequency domain using FFT. The global map is taken from the predicted vehicle position, and the local map is generated at the current vehicle position. Therefore, matching result is represented as the error correction value from the predicted vehicle position. Thus, if and are the error correction values, can be expressed as follows:
For phase correlation, we need to know the phase difference between two images. Thus, the normalized cross power spectrum is calculated, and the equation is as follows:
Inverse Fourier transform (IFT) can be used to obtain the final phase correlation from equation (3).
The matching result through phase correlation is shown in Figure 15.
As shown in Figure 15, sharp peaks are extracted using phase correlation. However, if the road marking information of the local map is insufficient, the peaks do not appear or multiple peaks may occur. To solve this problem, we applied a Gaussian filter. By setting the region of interest (ROI) using a Gaussian filter, noise and multiple peaks can be removed. The ROI is set by considering the vehicle speed based on the current estimated position using . The Gaussian filter equation is as follows:
Figure 16 shows the result of applying a Gaussian filter. Noise and multiple peaks are removed through the Gaussian filter, and a sharp peak is extracted. The result of the map matching is as follows:
Since the map used for map matching is a grid map, the estimated value oscillates if the maximum value of is taken as the matching result because it is as uncertain as the grid size. Therefore, by using the peak and value around the peak as in equation (6), the uncertainty is reduced and the oscillation is prevented.
5.1.2. Heading
The ICP algorithm [23] was used to estimate the vehicle heading. The binary road marking grid map is the stored road marking information in the grid. Therefore, for point matching, the global map and local map are converted to points based on the grid center.
The road marking in the grid is converted to a point, as shown in Figure 17. The local map used in position estimation is a precise local map using accumulated data. However, in the case of heading, accurate error estimation cannot be performed because the error is mixed by accumulating data. Therefore, when the heading is estimated, the local map is generated using the current scan data and matched with the global map.
The ICP algorithm is executed in two stages: (1)Find the corresponding point between the two data sets(2)Compute a transformation matrix that minimizes the distance between corresponding points
Repeat the above two steps until it is below the threshold to derive the final transformation matrix. In addition, is set to prevent fault matching by outliers. Algorithm 1 shows the ICP algorithm.

As shown in Figure 18, the points of the local map are corrected through the ICP algorithm.
Equation (7) is used to obtain the heading correction value using the transformation matrix. Since the global map is based on the predicted position of the vehicle and the local map is based on the current vehicle position, is the correction value of the heading.
There are two reasons why heading estimation is required for vehicle localization. First, if the position estimation is accurate but the heading error is large, the direction of the autonomous vehicle will be deviated from the lane continuously. Next, the precise local map generated accumulated increment of position and heading. We obtain the increment of position and heading through ICP scan matching. However, to match with the global map, we need to add the heading value from GPS/DR. Therefore, if the GPS/DR heading error is large, the road marking information of the generated precise local map does not match the global map.
Figure 19 shows the heading error of GPS/DR and estimated heading result. Figure 20 shows the result of overlapping of the global and local maps at the true reference position. The local map is represented by black, and the global map is represented by gray. In Figure 20, the local map is generated by using GPS/DR heading angle and estimated heading result, respectively. If the heading angle is same with true reference, global and local maps must be matched each other. As shown in this figure, (a) is not consistent with each other, while (b) is almost consistent. Consequently, if the heading estimation is not conducted, map matching error increases. Therefore, heading estimation is required for vehicle localization.
(a)
(b)
5.2. Evaluation of Map Matching Result
Vehicle localization through map matching can lead to large errors if the matching results are wrong. In particular, in the case of autonomous vehicles, a wrong localization result may cause a big accident. Therefore, when the map matching result is wrong, it should be detected and the wrong matching result must not be used. In this paper, we conducted an evaluation of position and heading matching results. Accordingly, only good results were used for localization.
5.2.1. Evaluation of Position Matching Result
In this paper, we used two parameters to evaluate of the position matching result. One is the protection level (PL) of the matching result, and the other is the peaktoside lobe ratio (PSR).
PL is a parameter [24] used to determine the reliability of aircraft GPS. When the PL exceeds a certain threshold, a false alarm indicates that the localization result is wrong. We applied that to the vehicle [25, 26]. The equation for the PL of a vehicle is as follows: where is the probability of missed detection and it is calculated as the Rayleigh inverse cumulative distribution function. The term can be estimated as the maximum eigenvalue of the covariance matrix.
PSR [27] is a parameter that determines the sharpness of the correlation peak. A sharp peak means that the matching result is reliable. Therefore, PSR can be a basis for evaluating of the matching result.
Figure 21 shows an example of the peak and side lobe. The peak means the peak of correlation, and side lobe means noise around the peak.
As shown in equation (9), the PSR can be obtained by using the peak and side lobe and the sharpness of the peak can be determined by the ratio of the peak to the side lobe.
The above two parameters are used to evaluate the position matching result. We set the threshold of the two parameters to evaluate the position matching result.
Figure 22 shows an example of the threshold setting for PL and PSR, where the blue points denote the matching result. In the case of PL, the lower the value and the higher the PSR, the higher the reliability of the matching result will be. Thus, out of the four ranges shown in Figure 22, (a) gives highly reliable results. Therefore, we used map matching results in the range (a) as a measurement in the navigation filter and did not use the results of the remaining range.
5.2.2. Evaluation of Heading Matching Result
The heading matching result is also evaluated by two parameters. One is the number of points, and the other is the distance ratio of the corresponding points. In the case of heading matching, since the ICP algorithm is used, the number of points plays an important role in convergence and derivation of the transformation matrix. In addition, the distance of the corresponding points can also be a determining factor for the matching performance because the ICP algorithm is used to find the transformation matrix with the closest distance of corresponding points. The distance between the corresponding points was derived using the nearest neighbor search (NNS) algorithm. Next, we extracted corresponding points whose distances are within a certain threshold. Finally, the ratio between the extracted points and the total number of points was used as a parameter for evaluating the heading matching result. We set the threshold of the two parameters and used the result in the range as a measurement value in the navigation filter.
5.3. Navigation Filter Configuration
In this paper, the KF [28] is used as a navigation filter for vehicle localization. The KF is the most common filter used for localization. To use this filter, it is necessary to define the state vector and input vector.
In equation (10), is the state vector that represents the position and heading of the vehicle in the global frame and is the input vector that refers to the increment of the position and heading derived from the ICP scan matching algorithm. The state variables are predicted by the time update. The time update is as follows:
In equation (11), we need to initialize the state vector and covariance matrix before the time update. The initial state vector is obtained from the GPS/DR position and heading angle. However, since the GPS/DR position and heading error are large, map matching is performed once to correct position and heading. Because the position and heading error are large, they cause a large residual in the measurement update. After the initialization, the time update is conducted by the input vector. The increments from ICP scan matching are based on the vehicle frame. Therefore, we need to transform to the world frame through the transformation matrix . After the time update, we derive and . is the predicted state vector, and is the predicted covariance matrix.
In equation (12), we defined observation matrix and measurement residual. The measurement residual is derived from map matching.
Equation (13) is the measurement update. After the measurement update, we derive and . Further, and are used for the next state estimation.
6. Experimental Result
In this section, we analyze the performance of precise vehicle localization using a binary road marking grid map through experiments.
Figure 23 shows the experimental area. The experiment was carried out in Gangnam, Seoul, which is surrounded by buildings, and was carried out in a traffic congestion situation. In addition, data collected from two laps were completed in the experimental area and the length of the experimental area was approximately 5 km. Figure 24 is types of road markings in the experimental area. As shown in this figure, there are many kinds of road markings on the road. Therefore, it is possible to obtain robust localization performance by using them. We used a commercial GPS/DR (MicroInfinite CruzCore DS6200) and 32channel 3D LIDAR (Velodyne HDL32E). Generally, the localization method using 3D LIDARbased road marking grid maps has intensity information in the grid [6, 7]. Therefore, we compared vehicle localization results using the binary road marking map and the intensity road marking map. Through this, we can demonstrate the advantages of using the binary road marking map mentioned in Section 3, as well as the comparison of localization performance using intensity road marking map and binary road marking map.
Table 3 shows the ratio of the matching result reliability of the intensity map and binary map. It is calculated as the ratio between the entire map matching result and the reliable matching result derived from the evaluation of the position matching result in Section 5. As shown in this table, map matching using a binary map is more reliable than using an intensity map and better matching performance can be expected.

In Figures 25 and 26, the blue line represents the position error of the GPS/DR, the red line represents the position error of the localization result using the intensity map, and the black line represents the position error of the localization result using the binary map. As shown in Table 4, the lateral and longitudinal RMS position errors using the intensity map are greater than those using the binary map because the map matching using the intensity map has an ambiguity greater than the map matching using the binary map. Intensity data provides 8 bits of information. Therefore, the difference in the intensity level between the global map and local map affects the localization performance. In particular, when we use 3D LIDAR, the intensity level depends on the vertical layers, even though the same road marking is scanned. On the other hand, the binary map is made of binarized intensity. Therefore, it has information only in the form of 0 s and 1 s. Therefore, the above ambiguity problem does not occur.

In Figures 27 and 28, the blue line represents the CDF of the position error result using the intensity map and the red line represents the CDF of the position error result using the binary map. Table 5 gives the position error with a 99% confidence level. It also shows that the error with the intensity map is larger than that with the binary map. As shown in Table 5, the lateral error is under 0.5 m with a 99% confidence level using the binary map. It means that the lanekeeping condition is satisfied for autonomous driving. On the other hand, when using the intensity map, the lateral error is more than 0.5 m with a 99% confidence level. In addition, the maximum error is 0.70 m for the intensity map and 0.26 m for the binary map. Therefore, it is possible to maintain the lane when we use a binary map for vehicle localization, while lane departure may occur when we use an intensity map. Also, the longitudinal error is under 0.5 m with a 99% confidence level using the binary map. That is a very good result considering that the longitudinal error condition for safe autonomous driving is generally 1 m. On the other hand, the longitudinal error with a 99% confidence level exceeds 1 m when using the intensity map. In particular, the maximum error is 2.08 m for the intensity map and 0.55 m for the binary map. Vehicle localization using intensity map can create a dangerous situation due to a large instantaneous error. On the other hand, when we use a binary map for vehicle localization, the vehicle can be driven more safely because of the small error.

Figure 29 shows the heading error. The intensity map contains intensity information of all objects on the road. Thus, heading estimation cannot be done using the map matching algorithm described in Section 5 because we cannot extract only road marking information. Therefore, we compared the heading results between GPS/DR with those obtained by using the binary map. Table 6 shows the RMS heading error. As shown in Table 6, when we use the binary map, the heading performance can be improved. Moreover, the maximum heading error was 2.12° for GPS/DR and 1.45° for heading estimation using the binary map.

Thus far, we analyzed the localization result in the same road condition. However, common vehicle localization may have different road surface conditions. Therefore, we compared the matching performance for different road surface conditions, as shown in Figure 29. However, there is a problem that the ground truth cannot be set because the data are different from each other. Therefore, we analyzed the matching performance using the PSR. Because the PSR is a parameter that determines the sharpness of the correlation peak as described in Section 5, the higher the PSR, the higher the matching performance will be.
Figure 30 shows a comparison of PSR between different data, and Table 7 shows the average of PSR. The global map is a wet road surface, and the local map is a dry road surface. As shown in this figure, the PSR of the binary map is higher overall. Therefore, we can confirm that the matching performance of the binary map is better than the intensity map even in different road surface conditions.

Thus far, we analyzed the experimental result of the binary road marking grid mapbased precise vehicle position and heading estimation. The results show that the binary road marking map matching method can provide good position and heading estimation accuracy.
7. Conclusions
In this paper, we proposed the binary road marking grid mapbased precise vehicle localization method using 3D LIDAR. First, the intensity scan data obtained from 3D LIDAR was binarized using Otsu’s thresholding. After binarizing the intensity, we generated the binary road marking grid map and precise local map for map matching. The position and heading are estimated through the map matching and navigation filter. In the experimental result, the 2D lateral RMS error was approximately 0.05 m and longitudinal RMS error was approximately 0.08 m. In addition, a CDF analysis of the error showed that the lateral error was 0.21 m and longitudinal error was 0.36 m with a 99% confidence level. Thus, we confirmed that vehicle localization using the binary road marking grid map provides very precise and reliable estimation results.
The advantage of using the binary road marking grid map is as follows. First, a very small size map can be generated. A typical grid map increases the size exponentially when the grid size is smaller and map scale is bigger. However, the binary road marking grid map has only 1bit information. Therefore, map size can be reduced efficiently. Second, unlike the intensity road marking information where the level changes according to the road surface condition, binary road marking information provides reliable and accurate estimation results that can be used even if the intensity changes. Therefore, map matching using road marking is possible under any situation.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This research was supported by a grant (18TLRPC11326903) from the Transportation Logistics Research Program funded by the Ministry of Land, Infrastructure and Transport of Korean government.
References
 K. Ansari and Y. Feng, “Design of an integration platform for V2X wireless communications and positioning supporting CITS safety applications,” Journal of Global Positioning Systems, vol. 12, no. 1, pp. 38–52, 2013. View at: Publisher Site  Google Scholar
 J. Choi and M. Maurer, “Hybrid mapbased SLAM with RaoBlackwellized particle filters,” in 2014 IEEE Information Fusion (FUSION), 2014 17th International Conference, pp. 1–6, Salamanca, Spain, July 2014. View at: Google Scholar
 L. Li, M. Yang, C. Wang, and B. Wang, “Road DNA based localization for autonomous vehicles,” in 2016 IEEE Intelligent Vehicles Symposium (IV), pp. 883–888, Gothenburg, Sweden, June 2016. View at: Publisher Site  Google Scholar
 J. H. Im, S. H. Im, and G. I. Jee, “Vertical corner feature based precise vehicle localization using 3D LIDAR in urban area,” Sensors, vol. 16, no. 8, p. 1268, 2016. View at: Publisher Site  Google Scholar
 J. Levinson, M. Montemerlo, and S. Thrun, “Mapbased precision vehicle localization in urban environments,” Robotics: Science and Systems, vol. 4, p. 1, 2007. View at: Publisher Site  Google Scholar
 J. Levinson and S. Thrun, “Robust vehicle localization in urban environments using probabilistic maps,” in 2010 IEEE International Conference on Robotics and Automation (ICRA), pp. 4372–4378, Anchorage, USA, May 2010. View at: Publisher Site  Google Scholar
 R. W. Wolcott and R. M. Eustice, “Visual localization within LIDAR maps for automated urban driving,” in 2014 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2014), pp. 176–183, Chicago, USA, September 2014. View at: Publisher Site  Google Scholar
 M. Schreiber, C. Knöppel, and U. Franke, “LaneLoc: lane marking based localization using highly accurate maps,” in In 2013 IEEE Intelligent Vehicles Symposium (IV), pp. 449–454, Gold Coast, Australia, June 2013. View at: Publisher Site  Google Scholar
 M. Konrad, M. Szczot, and K. Dietmayer, “Road course estimation in occupancy grids,” in 2010 IEEE Intelligent Vehicles Symposium (IV), pp. 412–417, San Diego, USA, June 2010. View at: Publisher Site  Google Scholar
 Z. J. Chong, B. Qin, and T. Bandyopadhyay, “Synthetic 2D LIDAR for precise vehicle localization in 3D urban environment,” in 2013 IEEE International Conference on Robotics and Automation (ICRA), pp. 1554–1559, Karlsruhe, Germany, May 2013. View at: Publisher Site  Google Scholar
 D. Kim, T. Jung, and K. S. Yi, “Lane mapbased vehicle localization for robust lateral control of an automated vehicle,” Journal of Institute of Control, Robotics and Systems, vol. 21, no. 2, pp. 108–114, 2015. View at: Publisher Site  Google Scholar
 A. Y. Hata and D. F. Wolf, “Feature detection for vehicle localization in urban environments using a multilayer LIDAR,” IEEE Transactions on Intelligent Transportation Systems, vol. 17, no. 2, pp. 420–429, 2016. View at: Publisher Site  Google Scholar
 B. S. Reddy and B. N. Chatterji, “An FFTbased technique for translation, rotation, and scaleinvariant image registration,” IEEE Transactions on Image Processing, vol. 5, no. 8, pp. 1266–1271, 1996. View at: Publisher Site  Google Scholar
 A. Averbuch and Y. Keller, “FFT based image registration,” in IEEE International Conference on Acoustics Speech and Signal Processing, pp. 3608–3611, Orlando, FL, USA, May 2002. View at: Publisher Site  Google Scholar
 Velodyne LiDAR, I, HDL32E User’s Manual, 2012.
 N. Otsu, “A threshold selection method from graylevel histograms,” IEEE Transactions on Systems, Man, and Cybernetics, vol. 9, no. 1, pp. 62–66, 1979. View at: Publisher Site  Google Scholar
 A. Y. Hata and D. F. Wolf, “Road marking detection using LIDAR reflective intensity data and its application to vehicle localization,” in 2014 IEEE 17th International Conference on Intelligent Transportation Systems (ITSC), pp. 584–589, Qingdao, China, October 2014. View at: Publisher Site  Google Scholar
 C. C. Wang and C. Thorpe, “Simultaneous localization and mapping with detection and tracking of moving objects,” in 2014 IEEE International Conference on Robotics and Automation (ICRA), pp. 2918–2924, Washington Ds C, USA, May 2002. View at: Publisher Site  Google Scholar
 S. Thrun and M. Montemerlo, “The graph SLAM algorithm with applications to largescale mapping of urban structures,” The International Journal of Robotics Research, vol. 25, no. 56, pp. 403–429, 2006. View at: Publisher Site  Google Scholar
 R. Kümmerle, G. Grisetti, H. Strasdat, K. Konolige, and W. Burgard, “g2o: a general framework for graph optimization,” in 2011 IEEE International Conference on Robotics and Automation (ICRA), pp. 3607–3613, Shanghai, China, May 2011. View at: Google Scholar
 R. W. Wolcott and R. M. Eustice, “Robust LIDAR localization using multiresolution Gaussian mixture maps for autonomous driving,” The International Journal of Robotics Research, vol. 36, no. 3, pp. 292–319, 2017. View at: Publisher Site  Google Scholar
 K. W. Kim, B. H. Lee, J. H. Im, and G. I. Jee, “Intensity local map generation using data accumulation and precise vehicle localization based on intensity map,” Journal of Institute of Control, Robotics and Systems, vol. 22, no. 12, pp. 1046–1052, 2016. View at: Publisher Site  Google Scholar
 A. Segal, D. Haehnel, and S. Thrun, “GeneralizedICP,” in Robotics: Science and Systems, p. 435, Seattle, WA, USA, July 2009. View at: Publisher Site  Google Scholar
 J. Blanch, T. Walter, P. Enge et al., “Advanced RAIM user algorithm description: integrity support message processing, fault detection, exclusion, and protection level calculation,” in Proceedings of the 25th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2012), pp. 2828–2849, Nashville, USA, September 2012. View at: Google Scholar
 R. ToledoMoreo, D. Bétaille, and F. Peyret, “Lanelevel integrity provision for navigation and map matching with GNSS, dead reckoning, and enhanced maps,” IEEE Transactions on Intelligent Transportation Systems, vol. 11, no. 1, pp. 100–112, 2010. View at: Publisher Site  Google Scholar
 D. Margaria and E. Falletti, “The local integrity approach for urban contexts: definition and vehicular experimental assessment,” Sensors, vol. 16, no. 2, p. 154, 2016. View at: Publisher Site  Google Scholar
 E. O. Omidiora, S. O. Olabiyisi, J. A. Ojo, A. AbayomiAlli, F. Izilein, and P. I. Ezomo, “Mace correlation filter algorithm for face verification in surveillance scenario,” International Journal of Computational Science and Engineering, vol. 18, no. 1, pp. 6–15, 2013. View at: Google Scholar
 G. Bishop and G. Welch, “An introduction to the Kalman,” in Proceedings of SIGGRAPH, Course 8, pp. 27599–23175, Los Angeles, CA, USA, August 2001. View at: Google Scholar
Copyright
Copyright © 2019 KyuWon Kim et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.