Research Article  Open Access
PCABased Line Detection from Range Data for Mapping and LocalizationAiding of UAVs
Abstract
This paper presents an original technique for robust detection of line features from range data, which is also the core element of an algorithm conceived for mapping 2D environments. A new approach is also discussed to improve the accuracy of position and attitude estimates of the localization by feeding back angular information extracted from the detected edges in the updating map. The innovative aspects of the line detection algorithm regard the proposed hierarchical clusterization method for segmentation. Instead, line fitting is carried out by exploiting the Principal Component Analysis, unlike traditional techniques relying on least squares linear regression. Numerical simulations are purposely conceived to compare these approaches for line fitting. Results demonstrate the applicability of the proposed technique as it provides comparable performance in terms of computational load and accuracy compared to the least squares method. Also, performance of the overall line detection architecture, as well as of the solutions proposed for linebased mapping and localizationaiding, is evaluated exploiting real range data acquired in indoor environments using an UTM30LXEW 2D LIDAR. This paper lies in the framework of autonomous navigation of unmanned vehicles moving in complex 2D areas, for example, being unexplored, full of obstacles, GPSchallenging, or denied.
1. Introduction
Among the technical challenges which drive the research activities carried out in the field of Unmanned Aerial Vehicles (UAVs), a major issue is to improve their level of autonomy. As discussed in [1], this is equivalent to strengthen their capabilities of autonomous guidance, navigation, and control, taking into account factors like environment and mission complexity, situational awareness, and realtime implementation.
For instance, in the case of largescale outdoor scenario, autonomous and safe navigation is ensured by the classical sensor fusion architectures integrating an Inertial Navigation System (INS) with GPS, typically indicated as GPSINS [2, 3], which has been extensively exploited by researchers considering fixedwing [4, 5], helicopter [6, 7], and multirotor [8] UAVs. Nevertheless, a wide range of both military and civil applications, for example, urban and indoor surveillance, infrastructure monitoring, and exploration, require a microUAV (MAV) to be able to navigate in more complex environments, such as unknown areas, full of static, and/or mobile obstacles in which the GPS signal may be completely absent (GPSdenied, e.g., indoor) or unreliable due to multipath, absorption, and jamming phenomena (GPSchallenging, e.g., urban or natural canyons). Hence, alternative solutions have to be found, which involve the use of exteroceptive sensors that can be active, such as RADAR [9], LIDAR [10], ultrasonic rangefinders [11], and ultrawideband (UWB) positioning system [12], passive, for example, monocular [13] and stereovision [14] cameras operating in the visible band of the electromagnetic spectrum, or hybrid, like RGBdepth cameras [15], as they simultaneously acquire passive RGB images and depth images of the same scene in an active way [16].
Indeed, the use of these sensors has provided a boost to the capability of navigating autonomously in complex environments by performing Simultaneous Localization and Mapping (SLAM) [17, 18], that is, the realtime process by which a mobile marine/ground/aerial robot, moving within an unknown area, computes its own trajectory (localization) while simultaneously building a map of that environment (mapping). However, the implementation of SLAM algorithms on board MAVs still presents significant technical challenges, mainly due to the limited resources available on board for sensing and computation, as well as their complex and fast 3D dynamics [19]. On this basis, the adoption of featurebased approaches for mapping is extremely promising as it can lead to a sparse but compact representation of the observed environment, with higher speed of execution and reduced amount of data storage with respect to alternative occupancy grid techniques [20]. Among the different technological solutions suitable for SLAM (such as visionbased [21], LIDARbased [22, 23], and RGBdepthbased [24]), LIDAR are more robust than passive vision systems to operate in indoor environment, especially in case of low or nolight conditions. Also, they directly provide range data. However, LIDARbased SLAM presents significant technical challenges, mainly related to the lower resolution and limited texture of the available measurements with respect to images provided by passive vision sensors. In order to deal with these issues, strong research efforts have been recently addressed to the possibility of augmenting LIDARbased SLAM exploiting intensity information [25] or detection of robust features [26]. Instead, a different investigated solution consists in integrating passive vision and active LIDAR sensors [27, 28] in order to exploit complementary advantages and compensate their drawbacks.
In this framework, the aim of this paper is to present a technique for fast and robust detection of line features from range data provided by a 2D LIDAR. Indeed, these sensors, due to their limited size and weight, are particularly suited for installation on board MAVs operating in indoor areas, for example, officelike environments, but also outdoor due to their limited sensitivity to ambient light variations. Line extraction algorithms are typically characterized by two steps, namely, segmentation and line fitting. The former identifies the number of potential lines and associates the measured points to each candidate, while the latter allows extracting the line parameters from the aggregates of points. A review of the techniques for line detection from 2D LIDAR data, mostly used in the literature, is presented in [29]. Six different approaches to the segmentation task, namely, the splitandmerge algorithm [30, 31], the lineregression algorithm [32], the incremental algorithm [33] (also known as linetracking [34]), RANSAC [35], the Hough Transform [36], and the ExpectationMaximization (EM) algorithm [37], are compared in [29] in terms of speed, complexity, correctness, and precision by means of experimental tests on real scan data collected in a largescale office environment. In [29], each one of these techniques is further accelerated by implementing clusterization, which allows coarsely dividing a raw laser scan into contiguous groups (clusters). With regard to line fitting, [29] states that stateoftheart approaches rely on the least squares (LS) method.
The proposed line detection algorithm provides solutions for both segmentation and line fitting. The issue of developing and testing robust techniques for segmentation has critical importance due to the challenges of dealing with clutter in 2D environments. In this respect, an original twolevel clusterization algorithm is presented which includes a merging check to avoid the possibility that an edge in the scene is split in two or more parts only due to the measurements noise. This solution for segmentation can also be adapted to represent both nonlinear and piecewise linear structures. An algorithm based on the Principal Component Analysis (PCA) [38] is proposed for line fitting. Specifically, it relies on the capability of the PCA to estimate the principal directions of the 3D datasets acquired by a 2D LIDAR. PCA has already been used to detect lines from row and column edges extracted in 2D images [39]. Regarding applications with range data, PCA is used in [40] to identify feature points belonging to planar structures and in [41] to recognize echoes produced by polelike objects and to extract the direction from roadshaped point clouds in 3D urban environments. In both [40, 41], data from 3D LIDAR systems, for example, the Multiplatform Mobile Laser Scanning [42], are processed. Such a system is not suitable for application on board MAVs. PCAbased methods have also been widely exploited in the open literature for pose estimation applications in different scenarios [43, 44]. Although PCA is not a new concept for line fitting, a detailed analysis of its performance when used with LIDAR range data is not available in the literature. On this basis, performance of the proposed PCAbased line detection algorithm is assessed by means of both numerical simulations and experimental tests. A specific analysis is carried out to compare the PCA to the LS approach for line fitting in terms of computational load and accuracy level. This is important to demonstrate the applicability of the proposed approach for implementation on board MAVs, since LS is a standard and extremely fast method to perform linear regression given a set of range 3D data.
An additional aim of this paper is to demonstrate the possibility of the PCAbased line detection algorithm to be used within a 2D SLAM process. In this respect, it is worth outlining that many UAV applications (e.g., indoor inspection and monitoring) involve the necessity of flying in 2D environments. In this context, the capability of estimating the horizontal motion and simultaneously a 2D map is relevant, also considering that the height above ground can be measured by using additional sensors suitable for the installation on board MAVs (e.g., ultrasonic rangefinders) or also by redirecting downward a narrow portion of the scan emitted by the 2D LIDAR using mirrors [19]. Hence, a linebased mapping algorithm is here presented. Specific attention is addressed to the procedures foreseen for handling the features to be stored in the updating map, as well as for improving the accuracy of the position and attitude solution provided by a localization algorithm (localizationaiding). A SLAM algorithm using a probabilistic PCA algorithm for line detection from 2D LIDAR data is presented in [45]. However, results obtained in [45] applying offline the algorithm on real data put into evidence errors in the estimated maps due to the complexity of correctly merging clusters from subsequent point clouds. This justifies the need of further investigations to demonstrate reliability and effectiveness in performing this task. Another example of localizationaiding exploiting line detection is discussed in [46]. Specifically, the relative displacement and inclination between corresponding lines in consecutive frames are integrated within an Extended Kalman Filter. Instead, in this paper, correspondences between new lines extracted from the actual frame and lines already stored in the updating map are determined using a more robust approach based on multiple angular and position metrics, and they are adopted to directly correct the vehicle’s position and the heading angle. However, it is worth outlining that the development and performance analysis of a complete SLAM framework based on the proposed algorithms for mapping and localizationaiding is out of the scope of this paper.
This paper is organized as follows. Section 2 describes in detail the proposed line detection algorithm operating on range data, as well as the linebased mapping technique within which the algorithm is implemented, including the approach proposed for localizationaiding. Section 3 presents the numerical simulations which are carried out to compare PCA and LS methods for line fitting in terms of computational load and accuracy level. The experimental setup conceived to collect real data and the results obtained applying the PCAbased line detection algorithm as well as an example of full implementation of localization and mapping are described in Section 4. Finally, Section 5 contains the conclusion and indications concerning possible future works.
2. PCABased Line Detection
The PCAbased algorithm proposed for line detection is conceived to operate receiving as input only range measurements provided by a 2D LIDAR. Typically, the data collected by these sensors are given in output according to a representation in polar coordinates, that is, the scan angle (), which is the angular position in the LIDAR FieldofView (FOV) of each emitted ray with respect to a reference direction (usually it is the boresight one), and the range (), which is measured by exploiting the TimeofFlight (TOF) principle. However, the measured 2D point clouds will be expressed in the following also using the representation in rectangular coordinates ( and ) which can be obtained by applying
Once a 2D point cloud is acquired by a 2D LIDAR, the PCAbased line detection algorithm starts with the segmentation step, which is carried out by performing a clusterization process characterized by two hierarchical levels. At the higher level, the scan is coarsely subdivided into main clusters which represent clearly separated portions of the observed scene. This is done by looking for peaks of the interpoint distance (i.e., the distance between consecutive points in the scan) as a function of . These peaks, indicated as breakpoints, represent the beginning of any individual cluster which is found moving along the scanning direction, and they are recognized as those locations where the interpoint distance is larger than a specific threshold (DTh). The correct value of DTh can be easily selected, considering the specifications of the available sensor in terms of angular resolution and operating range, in order to avoid separating points belonging to same main cluster. For instance, the UTM30LXEW produced by Hokuyo is characterized by an angular resolution of 0.25° and a maximum range of 30 m. Hence, DTh is set to 15 cm considering that the maximum interpoint distance at a range of 30 m is 13.1 cm. Figure 1 shows an example of implementation of this step of the algorithm over a point cloud acquired within a 2D corridor environment, using the UTM30LXEW LIDAR.
Four separate main clusters are identified in Figure 1. Before moving to the next phase, it is worth mentioning that if a cluster is too small (typically less than 5 points) it is neglected from the line detection process. Once the main clusters are separated from each other, a 2ndlevel (lowlevel) clusterization is required to deeply search for subclusters potentially corresponding to different linear structures in the observed scene. Specifically, this step looks for local peaks of the two derivatives with respect to of the functions and by setting, as new breakpoints, those locations (i) at which either (2a) or (2b) is satisfied:
In the equations above, where is the number of points in the considered main cluster, the threshold assigned to identify the local peaks of and is the sum of their mean and standard deviation computed over the cluster. This original approach is particularly convenient to better describe quite complex indoor environments where, for instance, the walls are characterized by frequent corners due to the presence of doors, windows, or other obstacles. By looking at Figure 2, it is clear how the number of breakpoints increases after the 2ndlevel clusterization is applied with respect to Figure 1.
Of course, due to the measurement noise, some of the new breakpoints do not correspond to the beginning of a different linear structure. This issue is solved, after line fitting, by implementing an additional step necessary to merge adjacent lines corresponding to the same linear structure. Before entering this detail, it is necessary to focus on the line fitting step, which relies on the PCA, that is, a technique able to determine the principal directions of a multidimensional dataset by analyzing the eigenvectors and eigenvalues of its covariance matrix (). In this case, given a generic subcluster (composed of measurements), can be estimated by applying
Hence, a line is extracted from the subcluster if the ratio, , between the two eigenvalues of (which is a measure of the elongation of the pointset) is larger than a specific threshold (ETh). Indeed, if the pointset does not have a linear structure, the value of will not be far from 1. So, the choice of ETh mainly affects the capability to declare a subcluster as line depending on its length, and it must be done considering the range measurement accuracy of the available LIDAR. Of course, the larger ETh is, the more elongated the selected pointset must be to be recognized as line. For instance, if ETh is set to 100, the probability to identify a 75 cm linear structure at a distance of 10 m (assuming a std of 4 cm in the direction perpendicular to the linear one) is 1.35% (evaluated over 10^{5} simulations). If ETh is increased to 380 this probability goes to zero. Specific numerical simulations (not discussed for the sake of brevity) have demonstrated that values of ETh around a few tens are adequate to robustly extract even very small linear structures (down to 10 cm) from the observed scene.
If the condition on is satisfied, the parameters suitable to univocally identify the line can be computed, for example, on the basis of the classical Cartesian representation shown in where and are the coordinates of a generic point of the line, is the angular coefficient, and is the known term. As the direction of the line () is given by the eigenvector (whose components are and ) corresponding to the maximum eigenvalue of , is derived usingwhile is obtained by assuming that the centroid of the subcluster (whose coordinates are and ) belongs to the line, as shown in
Of course, the two ends of the line must be assigned in order to obtain a finite line segment (which models the observed linear structure). This is done by projecting the first and last points of the subcluster in the line direction. Finally, the previously mentioned merging step allows joining two adjacent subclusters if the corresponding lines satisfy the conditions highlighted inwhere is the orientation of the direction of the segment which links their centroids, is an angular threshold, and the subscripts “1” and “2” are used to distinguish between the two lines. Also, the symbols and are used to indicate the intersection and union operators from set theory, respectively. Small values of (e.g., 0.05°) allow the method to recognize even slight variation in the orientation of the observed linear structure. As a result, Figure 3 is obtained by applying this line fitting approach (setting ETh to 50 and including the merging step) to the clusterized dataset shown in Figure 2.
2.1. LineBased Mapping and LocalizationAiding
The algorithm proposed for linebased mapping is conceived to build a 2D map of the observed environment in an incremental way, that is, by progressively adding and updating lines as they are detected through subsequent acquisitions of range data. Of course, in accordance with the concept of SLAM, an adequate solution in terms of position and attitude must be provided by a localization algorithm. Indeed, the first step consists in applying this pose solution to convert the acquired point cloud from the sensor reference frame (SRF) to the inertial frame (e.g., EastNorthUp or NorthEastDown with fixed origin) in which the trajectory of the moving platform and the map of the observed scenario are determined. Afterwards, the PCAbased line detection algorithm is implemented and the extracted line segments are used to update the map. Specifically, each of this segment is a potential candidate to contribute to the map generation process, but a decisionmaking step is required to establish whether or not it corresponds to a linear structure already contained (at least partially) in the map. To this aim, each of the new detected lines is compared to all the existing ones in the map using two metrics, namely, and . The former represents the difference between the two compared lines in terms of orientation () which can be easily obtained from (5). The latter, instead, is a distance metric derived by summing two contributions as shown in is the difference between the two compared lines in terms of the distance from the map origin () which is computed using
Instead, is the distance between the centroids of the two compared lines, whose coordinates are and , respectively. can be derived using
Given a line detected in the current dataset, if it is not possible to find an already stored segment for which both metrics are kept below a minimum threshold assigned taking localization errors into account (e.g., 15° and 1 m, resp.), it is classified as a new element of the map. Conversely, the detected line is matched to the line in the map for which is the minimum. Indeed, among the two metrics, is certainly the most important to ensure correct line matches, as in many indoor scenarios (e.g., corridors); it is likely to find more than one linear structure with the same orientation ( is similar) but located at different positions. Once all the possible matches are found, the corresponding elements of the map must be updated. This is done by projecting the two ends of the new line on the direction of the old one, so that the two updated ends are the ones for which the length of the line is the maximum.
The matches between new detected lines and old elements of the map can be used for localizationaiding, that is, to improve the accuracy of the position and attitude estimates provided by a localization algorithm. Here, a strategy for localizationaiding is proposed which exploits angular information extracted from the previously mentioned line matches. Specifically, if is the total number of line matches at a generic time (), the heading solution at that time can be corrected as shown inwhere the subscripts “OLD” and “NEW” indicate the original (provided by the localization algorithm) and updated estimates of the heading angle (), respectively, while the correcting factor () is the weighted average of the angular offsets between the matched lines, as shown in
The values of the weights () in (12) are determined using the distance metric as shown in (13), thus giving more relevance to those matches which are closer in terms of location occupied in the map
At this point, also the position can be updated usingwhere is the previous time instant and the subscripts “OLD” and “NEW” are again used to indicate the original (provided by the localization algorithm) and updated estimates of the vehicle’s position ().
3. Numerical Simulation Results
This section presents description and results of the numerical simulations carried out to assess absolute performance of the PCA when used for line fitting and to get a comparison with the traditional LS method [47], for which a fast version is implemented [48]. The main aspects which are deemed relevant to this analysis are the attained accuracy level and the computational load.
Since the goal of these simulations is to assess line fitting performance of the PCA method over sets of points distributed according to a linear pattern, a realistic reproduction of the operation of a range sensor (e.g., 2D LIDAR) is not carried out. In this respect, the simulations simply rely on the generation of a set of points, randomly distributed along a linear direction identified by the angular coefficient (), so that the point coordinates ( and ) can be determined usingwhere is the point separation and is the th extraction from a normal distribution with zero mean and standard deviation equal to . An example of simulated dataset is shown in Figure 4.
The effect on the results of considering different settings for and is analyzed. Indeed, is a measure of the angular resolution of the range data (although its variation can also reproduce the effect of modifying the distance of the sensor from the scene given a fixed resolution), while simulates the sensor noise in the measured distance. On the other hand, is set to 0 as line fitting performance is independent of the orientation of the linear dataset. As a consequence of this choice, the length of the simulated linear dataset (), which also has an impact on line fitting performance, is determined by the selected values of .
Once the setting parameters of the simulated range data are chosen, both PCA and LS methods are applied obtaining the corresponding lines, each one identified by an angular coefficient ( and , resp.) and a constant term ( and , resp.). Finally, the line fitting accuracy is evaluated as the standard deviation of the distance of the assigned points from the estimated line ( and , resp.), according to
It is worth mentioning that, for any selected configuration of the setting parameters, the results are averaged on 10000 simulations thus having the possibility of carrying out meaningful statistical analyses. Firstly, different values of are considered, ranging from 5 mm to 10 cm, thus being consistent with typical 2D LIDAR specification in terms of angular resolution (e.g., 0.1°–0.25°) taking also their typical operational range (up to a few tens of meters) into account. Instead, is kept fixed to 4 cm, which is also in line with typical range noise of a 2D LIDAR. The analyzed configurations and the corresponding simulation results are summarized in Table 1 where the computational time ( and , resp.) is expressed in milliseconds (ms).

It can be clearly stated that the computational time provided by the PCA (largely less than 1 ms) is comparable to the result of the LS approach. Also, there are no significant differences in the accuracy level. Indeed, the difference between the root mean squares (rms) of the estimation error provided by the two methods is always negligible if compared to the absolute performance attained by the PCA. This result is not affected by variation of density of the simulated dataset, which is represented by . As expected, the longer the dataset is (for a given value of ), the higher the associated value of becomes, and consequently the closer PCA and LS accuracy levels get. Indeed, the difference reduces of two orders of magnitude, that is, from 10^{−1}% to 10^{−3}%, as goes from 1 m to 10 m. Additional numerical simulations are carried out to show how the results presented above can be affected by a variation in the level of measurement noise. Specifically, is set to 5 cm and is set to 20, so that is 1 m, thus reproducing a quite sparse dataset characterized by a linear spatial distribution. Instead, five different values of are considered ranging from 5 mm (bestcase) to 20 cm (worstcase). Simulation results are collected in Table 2.

Before discussing the absolute accuracy level attained by the PCA and the comparison with the LS performance, it is worth outlining that three of the selected values of , that is, 5 mm, 1 cm, and 5 cm, are consistent with the nominal levels of performance that scanning LIDAR, exploiting the TOF principle for distance measurement, are able to ensure. Indeed, highperformance systems may be characterized by range accuracies from around 1 cm down to even submillimeter level, while, for the majority of scanning LIDAR, never exceeds the interval going from 3 cm to 5 cm. Conversely, the remaining values to which is set in the simulations (i.e., 10 cm and 20 cm) are representative of worstcase conditions which may only occur if the sensor receives toomuch energy due to, for instance, direct sun illumination in outdoor environment. On this basis, results in Table 2 show that the PCA is able to carry out line fitting with the same accuracy level as the LS method, even operating on sparse data characterized by large nonnominal values of the range measurement noise.
Finally, it is interesting to carry out further simulations to evaluate the impact on line fitting accuracy caused by the presence of outliers in the considered dataset. Since we are just dealing with the problem of line fitting, outliers must be modeled as measured points which have been erroneously attributed to a given cluster. Specifically, outliers are generated by randomly extracting (according to a uniform probability distribution) a percentage of points () in the dataset for which the displacement from the reference direction is equal to . Again, it is important to stress that even if these simulated data are not representative of real laser scan data, they can be considered adequate for the analysis of line fitting performance as the segmentation necessary to separate clusters belonging to different linear structures is not addressed. For this analysis, is set to 2.5 cm and is set to 200, so that is 5 m, while is 4 cm. Simulation results are collected in Table 3, where four different values, namely, 1%, 5%, 10%, and 20%, are considered for .

Results in Table 3 show that the PCAbased line fitting technique has comparable performance to LS and ensures cmlevel accuracy even setting to 20%, which in the analyzed case corresponds to 40 outliers over the total number of 200 points composing the dataset.
4. Tests on Experimental Data
The scope of this section is to present the results obtained applying the proposed algorithms and techniques on real range data. To this aim, it is necessary to describe the experimental setup conceived for data collection as well as to clarify in which scenarios the experimental tests are carried out. As regards the setup, the selected 2D LIDAR is the UTM30LXEW produced by Hokuyo, for which the main specifications and a full characterization can be found in [49] and [50], respectively. This sensor is installed on a portable platform together with one autopilot, that is, the Pixhawk produced by 3D robotics, which is used to obtain reference information regarding the attitude of the platform, and one embedded board, that is, the Nitrogen6X produced by BoundaryDevices. This latter component is used to register data from both the Pixhawk (using USB connection) and the LIDAR (using Ethernet connection) by exploiting the corresponding nodes of the Robot Operating System (ROS) [51], that is, the mavros and the urg_node, respectively. This ensures simultaneous recording of range data from the LIDAR and inertial data from the Pixhawk, together with their timestamps, within the same bagfile. The data acquisition is commanded by a laptop exploiting a WiFi connection. A conceptual representation of the experimental setup is shown in Figure 5.
This setup has been first used to perform static acquisition of range data within indoor cluttered environments considering both limitedsize areas, like corridors or offices, and wider scenarios such as halls. Figure 6 shows a portion of the scaledplanimetry of the ground floor of Building 83 (which is part of the School of Engineering) inside the Cranfield University Campus. In this planimetry, red polygons are used to highlight the areas where LIDAR data have been collected.
(a)
(b)
(c)
(d)
(e)
Each horizontal sequence of images in Figure 7 shows partial (clusterization) and final (line fitting) outputs obtained implementing the PCAbased line detection algorithm over point clouds collected within the previously defined indoor scenarios. These results demonstrate the robustness of the proposed line detection approach. Indeed, it has been able to extract a large number of linear structures characterized by largely variable lengths, that is, ranging from a few meters down to a few tens of centimeters. The standard deviation of the distance between the points belonging to each cluster and the corresponding detected line is of cmlevel. Also, the algorithm appears to be able to deal with the presence of clutter in the observed scene. This is particularly clear by looking at Figures 7(b), 7(c), and 7(e).
4.1. Localization and Mapping Experiments
Dynamic acquisitions have also been realized with the experimental setup transported by a human operator, thus partially reproducing the 3D dynamics which characterizes the flight of MAVs. The 2D maze shown in Figure 8 is the test area prepared for the dynamic experiment.
The major goal of these tests is the preliminary assessment of the algorithms presented in Section 2.1 for mapping and localizationaiding. This is done by implementing the entire localization and mapping process, in which the LIDAR/Inertial Odometry algorithm presented in [52] is in charge of the localization step. The structure of this localization and mapping process is summarized in Figure 9.
Once LIDAR and inertial data are collected exploiting the experimental setup described in the previous subsection, the localization and mapping algorithm are run offline in MATLAB environment. An example of application of this localization and mapping approach is shown in Figure 10.
By looking at Figure 10, it is clear that the linebased mapping technique allows generating a sparse but accurate representation of the observed scene. Indeed, all the edges of the real map (black lines) are adequately represented by the detected lines, in terms of length, location, and orientation. In this respect, the localizationaiding step has a fundamental role since it allows bounding potential drifts of the error in the estimated trajectory which could cause the edges of the estimated map to be characterized by position and/or angular offsets with respect to the real corresponding linear structure.
In order to better highlight the attained advantages, Figure 11(a) shows the correction provided by the localizationaiding step, according to (12), to the time behavior of the heading angle with respect to the solution given by the Pixhawk (). As a consequence, (14) can be applied thus making the estimated trajectory not to drift eastward as it occurs if the localizationaiding step is not implemented (as in [52]). This is highlighted in Figure 11(b).
(a)
(b)
5. Conclusions
This paper presented a new technique for line detection from range data (2D point clouds) provided by a 2D LIDAR, which is of interest to Unmanned Aerial Vehicles which need to carry out autonomous navigation applications, such as localization and mapping, in 2D cluttered environments.
As regards the two steps composing the algorithm, segmentation consisted of a hierarchicallevel clusterization approach suitable for quickly subdividing the measured point cloud into distinct aggregates of points representing potential candidates to be extracted as lines. Line fitting relied on the Principal Component Analysis (PCA) instead of exploiting the classical least squares (LS) linear regression as most of the stateoftheart approaches.
Numerical simulations were realized to assess line fitting performance of the PCA in comparison with a fast LS solution. The comparison was done in terms of accuracy and computational load, and the simulations considered the variability of dataset density and length, and robustness against noise and outliers. Although the algorithms were implemented offline in MATLAB, results allowed stating that the PCA is compatible with realtime implementation as it kept the runtime below 1 ms and close to the LS performance, even considering large clusters of 2000 points. Also, PCA ensured the same accuracy level of LS in spite of the variability of the characteristics of the simulated dataset.
Performance of the entire line detection process was evaluated over real LIDAR data acquired by using a purposely realized experimental setup. Runs over static acquisitions performed in different indoor environments, such as corridors, offices, and halls, demonstrated the effectiveness of the proposed method in extracting most of the linear structures in the scene with cmlevel accuracy and its robustness when dealing with the presence of high level of clutter.
In order to demonstrate the applicability of the proposed line detection algorithm for autonomous navigation in 2D environments, it was integrated within an innovative linebased mapping technique including also a new method for localizationaiding based on the extraction of angular information from the updating map. Dynamic acquisitions within a purposely prepared 2D test area were carried out to obtain a first performance assessment of the linebased mapping algorithm and to show the advantages provided by the introduction of the localizationaiding method. Indeed, the results showed the capability of this approach to generate a sparse but cmlevel accurate representation of the observed scene and to correct the heading and position solution provided, respectively, by an inertial unit and by a localization algorithm. Future works will be aimed at testing the proposed methods installing the experimental setup on board of a MAV, in order to fully assess performance in terms of localization and mapping capabilities.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was carried out in the frame of Program STARLinea 2, financially supported by UniNa and Compagnia di San Paolo.
References
 F. Kendoul, “Survey of advances in guidance, navigation, and control of unmanned rotorcraft systems,” Journal of Field Robotics, vol. 29, no. 2, pp. 315–378, 2012. View at: Publisher Site  Google Scholar
 R. L. Greenspan, “GPS and inertial integration,” in Global Positioning System: Theory and Applications, Volume 2, pp. 187–220, AIAA, Washington, DC, USA, 1996. View at: Google Scholar
 J. Farrell, Aided Navigation: GPS with High Rate Sensors, McGrawHill, New York, NY, USA, 2008.
 D. A. GrejnerBrzezinska, R. Da, and C. Toth, “GPS error modeling and OTF ambiguity resolution for highaccuracy GPS/INS integrated system,” Journal of Geodesy, vol. 72, no. 11, pp. 626–638, 1998. View at: Publisher Site  Google Scholar
 K. JongHyuk, S. Sukkarieh, and S. Wishart, “Realtime navigation, guidance, and control of a UAV using lowcost sensors,” in Field and Service Robotics, pp. 299–309, Springer, 2003. View at: Google Scholar
 J. N. Gross, Y. Gu, M. B. Rhudy, S. Gururajan, and M. R. Napolitano, “Flighttest evaluation of sensor fusion algorithms for attitude estimation,” IEEE Transactions on Aerospace and Electronic Systems, vol. 48, no. 3, pp. 2128–2139, 2012. View at: Publisher Site  Google Scholar
 H. J. Kim and D. H. Shim, “A flight control system for aerial robots: algorithms and experiments,” Control Engineering Practice, vol. 11, no. 12, pp. 1389–1400, 2003. View at: Publisher Site  Google Scholar
 J. Wendel, O. Meister, C. Schlaile, and G. F. Trommer, “An integrated GPS/MEMSIMU navigation system for an autonomous helicopter,” Aerospace Science and Technology, vol. 10, no. 6, pp. 527–533, 2006. View at: Publisher Site  Google Scholar
 A. F. Scannapieco, A. Renga, and A. Moccia, “Preliminary study of a millimeterwave FMCW InSAR for UAS indoor navigation,” Sensors, vol. 15, no. 2, pp. 2309–2335, 2015. View at: Publisher Site  Google Scholar
 D. Droeschel, J. Stuckler, and S. Behnke, “Local multiresolution representation for 6D motion estimation and mapping with a continuously rotating 3D laser scanner,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '14), pp. 5221–5226, Hong Kong, China, June 2014. View at: Publisher Site  Google Scholar
 N. Gageik, T. Müller, and S. Montenegro, “Obstacle detection and collision avoidance using ultrasonic distance sensors for an autonomous quadrocopter,” in Proceedings of the 1st Microdrones International Research Workshop (UAVWeek '12), Siegen, Germany, November 2012. View at: Google Scholar
 A. Alarifi, A. AlSalman, M. Alsaleh et al., “Ultra wideband indoor positioning technologies: analysis and recent advances,” Sensors, vol. 16, no. 5, article 707, 2016. View at: Publisher Site  Google Scholar
 J. Engel, J. Sturm, and D. Cremers, “Camerabased navigation of a lowcost quadrocopter,” in Proceedings of the 25th IEEE/RSJ International Conference on Robotics and Intelligent Systems (IROS '12), pp. 2815–2821, Vilamoura, Portugal, October 2012. View at: Publisher Site  Google Scholar
 M. Achtelik, A. Bachrach, R. He, S. Prentice, and N. Roy, “Stereo vision and laser odometry for autonomous helicopters in GPSdenied indoor environments,” in Unmanned Systems Technology XI, vol. 7332 of Proceedings of SPIE, Orlando, Fla, USA, April 2009. View at: Google Scholar
 A. S. Huang, A. Bachrach, P. Henry et al., “Visual odometry and mapping for autonomous flight using an RGBD camera,” in Proceedings of the 2012 International Symposium on Robotics Research, Taipei, Taiwan, August 2012. View at: Google Scholar
 S. Song, S. P. Lichtenberg, and J. Xiao, “SUN RGBD: a RGBD scene understanding benchmark suite,” in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition (CVPR '15), pp. 567–576, IEEE, Boston, Mass, USA, June 2015. View at: Publisher Site  Google Scholar
 H. DurrantWhyte and T. Bailey, “Simultaneous localization and mapping: part I,” IEEE Robotics and Automation Magazine, vol. 13, no. 2, pp. 99–110, 2006. View at: Publisher Site  Google Scholar
 T. Bailey and H. DurrantWhyte, “Simultaneous localization and mapping (SLAM): Part II,” IEEE Robotics and Automation Magazine, vol. 13, no. 3, pp. 108–117, 2006. View at: Publisher Site  Google Scholar
 A. Bachrach, R. He, and N. Roy, “Autonomous flight in unknown indoor environments,” International Journal of Micro Air Vehicles, vol. 1, no. 4, pp. 217–228, 2009. View at: Publisher Site  Google Scholar
 Z. Lu, Z. Hu, and K. Uchimura, “SLAM estimation in dynamic outdoor environments: a review,” in Proceedings of the 2009 Intelligent Conference on Robotics and Applications, pp. 255–267, Singapore, December 2009. View at: Google Scholar
 R. Munguía, S. Urzua, Y. Bolea, and A. Grau, “Visionbased SLAM system for unmanned aerial vehicles,” Sensors, vol. 16, no. 3, article no. 372, 2016. View at: Publisher Site  Google Scholar
 S. Kohlbrecher, O. Von Stryk, J. Meyer, and U. Klingauf, “A flexible and scalable SLAM system with full 3D motion estimation,” in Proceedings of the 9th IEEE International Symposium on Safety, Security, and Rescue Robotics (SSRR '11), pp. 155–160, IEEE, Kyoto, Japan, November 2011. View at: Publisher Site  Google Scholar
 J. Tang, Y. Chen, X. Niu et al., “LiDAR scan matching aided inertial navigation system in GNSSdenied environments,” Sensors, vol. 15, no. 7, pp. 16710–16728, 2015. View at: Publisher Site  Google Scholar
 K. Di, Q. Zhao, W. Wan, Y. Wang, and Y. Gao, “RGBD SLAM based on extended bundle adjustment with 2D and 3D information,” Sensors, vol. 16, no. 8, p. 1285, 2016. View at: Publisher Site  Google Scholar
 R. A. Hewitt and J. A. Marshall, “Towards intensityaugmented SLAM with LiDAR and ToF sensors,” in Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS '15), pp. 1956–1961, IEEE, Hamburg, Germany, October 2015. 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, article no. 1268, 2016. View at: Publisher Site  Google Scholar
 J. Zhang and S. Singh, “Visuallidar odometry and mapping: lowdrift, robust, and fast,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '15), pp. 2174–2181, IEEE, Seattle, Wash, USA, May 2015. View at: Publisher Site  Google Scholar
 E. López, R. Barea, A. Gómez et al., “Indoor SLAM for micro aerial vehicles using visual and laser sensor fusion,” in Robot 2015: Second Iberian Robotics Conference, pp. 531–542, Springer, 2016. View at: Publisher Site  Google Scholar
 V. Nguyen, A. Martinelli, N. Tomatis, and R. Siegwart, “A comparison of line extraction algorithms using 2D laser rangefinder for indoor mobile robotics,” in Proceedings of the IEEE IRS/RSJ International Conference on Intelligent Robots and Systems (IROS '05), pp. 1768–1773, Edmonton, Canada, August 2005. View at: Publisher Site  Google Scholar
 L. Zhang and B. K. Ghosh, “Line segment based map building and localization using 2D laser rangefinder,” in Proceedings of the 2000 IEEE International Conference on Robotics and Automation, pp. 2538–2543, San Francisco, Calif, USA, April 2000. View at: Google Scholar
 C. Berger, “Toward rich geometric map for SLAM: online detection of planes in 2D LIDAR,” Journal of Automation Mobile Robotics and Intelligent Systems, vol. 7, pp. 36–41, 2013. View at: Google Scholar
 K. O. Arras and R. Y. Siegwart, “Feature extraction and scene interpretation for mapbased navigation and map building,” in Proceedings of the Symposium on Intelligent Systems and Advanced Manufacturing, pp. 42–53, Pittsburgh, Pa, USA, October 1997. View at: Publisher Site  Google Scholar
 J. Vandorpe, H. V. Brussel, and H. Xu, “Exact dynamic map building for a mobile robot using geometrical primitives produced by a 2D range finder,” in Proceedings of the 1996 IEEE International Conference on Robotics and Automation, pp. 901–908, Minneapolis, Minn, USA, April 1996. View at: Google Scholar
 A. Siadat, A. Kaske, S. Klausmann, M. Dufaut, and R. Husson, “An optimized segmentation method for a 2D laserscanner applied to mobile robot navigation,” in Proceedings of the 3rd IFAC Symposium on Intelligent Components and Instruments for Control Applications, Annecy, France, 1997. View at: Google Scholar
 M. A. Fischler and R. C. Bolles, “Random sample consensus: a paradigm for model fitting with applications to image analysis and automated cartography,” Communications of the ACM, vol. 24, no. 6, pp. 381–395, 1981. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Li and E. B. Olson, “Extracting generalpurpose features from LIDAR data,” in Proceedings of the 2010 IEEE International Conference on Robotics and Automation, pp. 1388–1393, Anchorage, Alaska, USA, May 2010. View at: Google Scholar
 S. T. Pfister, S. I. Roumeliotis, and J. W. Burdick, “Weighted line fitting algorithms for mobile robot map building and efficient data representation,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 1304–1311, Taipei, China, September 2003. View at: Google Scholar
 S. Wold, K. Esbensen, and P. Geladi, “Principal component analysis,” Chemometrics and Intelligent Laboratory Systems, vol. 2, no. 1–3, pp. 37–52, 1987. View at: Publisher Site  Google Scholar
 Y.S. Lee, H.S. Koo, and C.S. Jeong, “A straight line detection using principal component analysis,” Pattern Recognition Letters, vol. 27, no. 14, pp. 1744–1754, 2006. View at: Publisher Site  Google Scholar
 C. Ulas and H. Temeltas, “A fast and robust featurebased scanmatching method in 3D SLAM and the effect of sampling strategies,” International Journal of Advanced Robotic Systems, vol. 10, no. 11, p. 396, 2012. View at: Google Scholar
 M. Lehtomäki, A. Jaakkola, J. Hyyppä, A. Kukko, and H. Kaartinen, “Detection of vertical polelike objects in a road environment using vehiclebased laser scanning data,” Remote Sensing, vol. 2, no. 3, pp. 641–664, 2010. View at: Publisher Site  Google Scholar
 A. Kukko, H. Kaartinen, J. Hyyppä, and Y. Chen, “Multiplatform mobile laser scanning: usability and performance,” Sensors, vol. 12, no. 9, pp. 11712–11733, 2012. View at: Publisher Site  Google Scholar
 R. Opromolla, G. Fasano, G. Rufino, and M. Grassi, “Large space debris pose acquisition in closeproximity operations,” in Proceedings of the 2nd IEEE International Workshop on Metrology for Aerospace (MetroAeroSpace '15), pp. 491–496, IEEE, Benevento, Italy, June 2015. View at: Publisher Site  Google Scholar
 Z. Cao, Y. Sheikh, and N. K. Banerjee, “Realtime scalable 6DOF pose estimation for textureless objects,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '16), pp. 2441–2448, IEEE, Stockholm, Sweden, May 2016. View at: Publisher Site  Google Scholar
 E. Brunskill and N. Roy, “SLAM using incremental probabilistic PCA and dimensionality reduction,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 342–347, Barcelona, Spain, April 2005. View at: Publisher Site  Google Scholar
 Y. Gao, S. Liu, M. M. Atia, and A. Noureldin, “INS/GPS/LiDAR integrated navigation system for urban and indoor environments using hybrid scan matching algorithm,” Sensors, vol. 15, no. 9, pp. 23286–23302, 2015. View at: Publisher Site  Google Scholar
 J. F. Kenney and E. S. Keeping, “Linear regression and correlation,” in Mathematics of Statistics, chapter 15, part 1, pp. 252–285, Van Nostrand, Princeton, NJ, USA, 3rd edition, 1962. View at: Google Scholar
 http://mathworld.wolfram.com/LeastSquaresFitting.html.
 Hokuyo, https://www.hokuyoaut.jp/02sensor/07scanner/download/products/utm30lxew.
 P. Demski, M. Mikulski, and R. Koteras, “Characterization of hokuyo utm30lx laser range finder for an autonomous mobile robot,” in Advanced Technologies for Intelligent Systems of National Border Security, pp. 143–153, Springer, Berlin, Germany, 2013. View at: Google Scholar
 M. Quigley, K. Conley, B. P. Gerkey et al., “ROS: an opensource Robot Operating System,” in Proceedings of the ICRA Workshop on Open Source Software, 2009. View at: Google Scholar
 R. Opromolla, G. Fasano, G. Rufino, M. Grassi, and A. Savvaris, “LIDARinertial integration for UAV localization and mapping in complex environments,” in Proceedings of the International Conference on Unmanned Aircraft Systems (ICUAS '16), pp. 649–656, Arlington, Va, USA, June 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Roberto Opromolla 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.