Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
AULD: Large Scale Suspicious DNS Activities Detection via Unsupervised Learning in Advanced Persistent Threats
Next Article in Special Issue
Real-Time Efficient Relocation Algorithm Based on Depth Map for Small-Range Textureless 3D Scanning
Previous Article in Journal
Gas Void Fraction Measurement of Gas-Liquid Two-Phase CO2 Flow Using Laser Attenuation Technique
Previous Article in Special Issue
Customized Mobile LiDAR System for Manhole Cover Detection and Identification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Target for Accurate Terrestrial Laser Scanning and Unmanned Aerial Vehicle Point Cloud Registration

1
Faculty of Civil and Geodetic Engineering, University of Ljubljana, Jamova cesta 2, 1000 Ljubljana, Slovenia
2
Leica Geosystems AG, Heinrich-Wild-Strasse, CH-9435 Heerbrugg, Switzerland
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(14), 3179; https://doi.org/10.3390/s19143179
Submission received: 31 May 2019 / Revised: 15 July 2019 / Accepted: 17 July 2019 / Published: 19 July 2019

Abstract

:
The main goal of our research was to design and implement an innovative target that would be suitable for accurately registering point clouds produced from unmanned aerial vehicle (UAV) images and terrestrial laser scans. Our new target is composed of three perpendicular planes that combine the properties of plane and volume targets. The new target enables the precise determination of reference target points in aerial and terrestrial point clouds. Different types of commonly used plane and volume targets as well as the new target were placed in an established test area in order to evaluate their performance. The targets were scanned from multiple scanner stations and surveyed with an unmanned aerial vehicle DJI Phantom 4 PRO at three different altitudes (20, 40, and 75 m). The reference data were measured with a Leica Nova MS50 MultiStation. Several registrations were performed, each time with a different target. The quality of these registrations was assessed on the check points. The results showed that the new target yielded the best results in all cases, which confirmed our initial expectations. The proposed new target is innovative and not difficult to create and use.

1. Introduction

Three-dimensional point clouds can be achieved through either active or passive remote sensing systems [1]. Active sensors, such as laser scanners, use the laser light to directly provide a precise and dense point cloud, which is properly scaled. Passive sensors, such as digital frame cameras, use feature matching in overlapping imagery (i.e., the automated identification of conjugate features in the involved images) to derive 3D information. With known interior and exterior orientation of the camera, a 3D point cloud can be computed with the intersection of image rays that belong to the image coordinates of conjugate points [2]. Interior orientation parameters can be obtained with camera calibration. Exterior orientation parameters can be derived through either a direct or indirect georeferencing procedure. A thorough review of image-matching algorithms is given in Reference [3]. The computer vision community introduced an approach for creating point clouds from overlapping images, known as structure from motion (SfM). In SfM, features (key points) are obtained from images using feature detectors and are matched to their descriptors. Different feature detectors, descriptors, and their combinations were tested in References [4,5,6]. In the next step, a sparse bundle adjustment is performed in order to calculate the camera parameters and transform image coordinates into 3D points. The sparse point cloud is then intensified using dense image-matching algorithms [7,8]. Structure from motion simultaneously estimates the exterior orientation parameters of the involved images and derives the 3D coordinates of corresponding features in those images. In the absence of control information, SfM-based point clouds are usually referred to an arbitrary reference frame [1]. The transformation of the point cloud in the reference coordinate system can be achieved by the use of ground control points (GCPs). These must be included as constraints in bundle adjustment rather than by just applying a 3D similarity transformation from an arbitrary reference system after SfM [9]. The SfM method has a large computational cost associated with image matching, and many efforts have been undertaken to reduce that cost. For example, a skeletal camera network was embedded into an SfM to enable efficient 3D scene reconstruction from a large set of unmanned aerial vehicle (UAV) images [10].

1.1. Overview of Point Cloud Registration Methods

To register point clouds acquired from different scenes, there are two main phases: The initial (coarse) alignment and the fine alignment [11,12,13]. Coarse alignment is best achieved by matching conjugate 3D features that have been extracted from two datasets [1,14]. Conjugate features could be points, lines, or planar features. The pairs of matched features are used to calculate the rigid body transformation from the “source” to the “target” point cloud. In cases where the scales of the registered point clouds are not uniform, a 3D similarity transformation is more appropriate [15]. In the next step, the registration process is refined using the fine alignment method. The most widely used fine registration method is the iterative closest point (ICP) algorithm [16] and its variants. The iterative closest point optimizes the transformation between two point clouds by minimizing the square errors between the correspondences. In its basic form, the number of points is selected from the point set to be registered, and the points corresponding to these points in the reference point set are then identified. The transformation is obtained by minimizing the distances between these pairs of points. Because of the irregular nature of point clouds, point-to-point correspondence cannot always be assumed [1]. Different variants of the ICP have been introduced for such situations. Improved methods based on ICP were reviewed in Reference [12]. Due to their automatic registration procedure, the variants of ICP algorithms are commonly used for registering point clouds [17,18,19,20,21,22,23,24,25].
The iterative closest point assumes a good estimation of the initial location of the point clouds to be registered. Therefore, the coarse registration already needs to provide a fairly accurate alignment of the two datasets. The complexity of coarse registering two datasets lies in the robust identification of conjugate features that must be spread as much as possible across the area. Feature extraction and a matching step can be approached by manual point-based matching, target-based matching, or feature-based matching [14]. The manual method requires the user to manually select matching points in both point clouds. This approach is not always reliable, mainly because of difficulties in selecting appropriate matching points. The target-based method aims at inserting objects (targets) into the scene, which can be easily recognizable and matchable in different point clouds. The feature-based method aims at automatically extracting and reliably matching 3D features from different datasets. According to the kind of feature used for registration, the methods can be put into three categories: Point-, line-, and plane-based methods [26].

1.2. Target-Based Registration

Point-based registration could be using point targets, individual points within the data, or key points that have characteristic attributes [1]. Target-based registration uses artificial targets, which are set up within the field of view of the data acquisition systems and are designed to facilitate their automatic identification from the derived point clouds.
Artificial targets are often used in registration procedures or when georeferencing a point cloud [27]. When deciding upon the suitability of an individual target, we need to consider its stability, rigidity, and surface reflectivity [28,29]. For targets that are used in laser scanning, it is characteristic that, in addition to the geometry of points, point intensity information is also used to determine the position of the target in the point cloud. The strongest signal should be returned from the center of the reflective target [30]. Targets are made from different materials. Retroreflective materials can be easily recognized even if very long distances are considered. Retroreflective targets reflect the incident beam along the same direction [31]. The varying reflectivity of materials causes different intensity values [32]. Target reflection properties also cause an offset in distance measurements [33,34]. Experimental results, such as those presented in Reference [31], have revealed the presence of large systematic range errors when scanning retroreflective target surfaces. A precise target-based positioning can only be achieved after modeling the systematic errors based on carefully designed calibration tests. The color of the targets is also important, because different colors produce different reflected intensities. A weighted total least squares method for precise black and white target measurement was proposed in Reference [35]. Artificial targets differ in shape, and geometrically they can be divided into plane and volume targets. Plane targets are more diverse, and certain equipment manufacturers have supplied us with a few characteristic ones. In addition, targets that are neither retroreflective nor made from special materials can be used for point cloud registration if they have a definite template shape [36]. Volume targets are designed in simple geometric forms, which are easy to model, such as spheres, cylinders, and cones. The use of volume targets for registering point clouds was studied and written about in References [37,38,39,40,41]. Rabbani and van den Heuvel [42] used planes, cylinders, and spheres as registration primitives. Their methodology of finding possible conjugate features was complemented by geometric constraints. Standard deviations when registering point clouds using their method were lower than those from ICP. The size and shape of the target should correspond to the maximum distance from the scanner, thus minimizing systematic errors due to the increase in laser spot size and assuring a high-density sampling of the surface [43]. The most important characteristic of the target is that it enables the simple modeling of its form, which in turn gives us precise coordinates of the reference target point in the various point clouds that we wish to merge. A reference target point is a geometrically defined point that can be uniquely defined in a point cloud, and in the case of plane targets also on images (center of a circle or a sphere, apex of a cone, etc.).

1.3. Feature-Based Registration

Linear and planar features have been repeatedly used as registration primitives for the coarse alignment of point clouds [1]. Researchers use those 3D features exclusively or in combinations to register point clouds. A methodology for registering photogrammetric datasets and light detection and ranging (LiDAR) data using 3D straight-line features instead of points was reported in Reference [15]. In Reference [44], line features were detected in both terrestrial and aerial laser scanner point clouds of built-up areas. The detected features were then used for registration in separate rotation and translation steps without any prior knowledge. The solutions were not stable, as the building interiors and exteriors could not be distinguished from line representation. In Reference [45], spatial curves were used as matching primitives to calculate the initial transformation parameters between the scanned point clouds of freeform surfaces (e.g., statues, cultural heritage artifacts). Al-Durgham and Habib [18] proposed an association matrix-based sample consensus approach for the registration of terrestrial laser scans using linear features. Their approach produced better results in fewer trials compared to the traditional random sample consensus (RANSAC) algorithm. A closed-form procedure for the coarse registration of 3D point clouds using automatically extracted linear features was presented in Reference [1]. The proposed method was capable of solving the registration problem in the presence of noisy measurements and providing reliable coarse alignment for a successful fine registration. The main limitation of the proposed approach was the reliance on manual identification of corresponding linear features in the point clouds. Bosché [14] automatically extracted planes to register terrestrial laser scanner (TLS) point clouds with industrial project 3D models. He calculated rigid transformation using two matches of nonparallel vertical planes and one match of horizontal planes. Three nonparallel planes may not always be present in construction scans, in which case the proposed approach fails. An alternative “2 planes + 1 point” registration approach was proposed in those cases. In Reference [46], terrestrial laser scanner (TLS) point clouds were registered with the method, which is based on virtual tie points that are generated by intersecting triples of scene planes. Planes in the point clouds were detected using RANSAC, which is often used to segment the planes in the point cloud. Other methods are also in use (e.g., a robust principal components analysis to segment linear and planar structures in point clouds was presented in Reference [47]). Xu et al. [48] reported an automatic marker-free method for fast and coarse registration between point clouds using the geometric constraints of planar patches, combining the voxel structure and a RANSAC-based strategy for selecting the corresponding plane pairs. The voxel structure suppressed the negative effects of outliers and unevenly distributed densities of points. In addition, using voxels instead of points to represent the data largely improved the efficiency of processing. Based on the voxel structure, a geometric constraint using small planar patches was designed to rapidly estimate the local coordinate frame of a scan in an urban scene. Compared to conventional ways of extracting planar surfaces with model fitting and region growing, the proposed planar constraint was generated in a more efficient way with a lower computational cost. In Reference [49], planar features were detected within unregistered scans. By registering planes instead of points, they reduced the amount of data that need to be processed by the hardware, and they successfully registered several thousand laser scans. Several researchers benefited by using combinations of extracted features to register point clouds. The least squares approach, which matches arbitrarily oriented 3D surface patches and curves, was presented in Reference [50]. In Reference [51], LiDAR point clouds were registered by employing point, line, and plane features. Depending on the scene geometry, each one of the features can be exclusively used or the features can be combined for calculating the transformation parameters. Due to combined measurements, the feature-based transformation model can afford a higher degree of flexibility and offer higher accuracy of point cloud registration tasks. Experimental results from this work showed that linear and planar primitives yield reliable estimates of the transformation parameters. Yao et al. [17] developed a RANSAC-based algorithm to register laser scans in a pairwise fashion, in which groups of linear and planar features are used as the registration primitives. They reported that the accuracy, robustness, and matching effectiveness of their registration method using planar features were better than with linear features. However, the algorithm was sensitive to the presence of repetitive patterns, and the registration of outdoor scenes usually failed. A non-iterative solution to estimate rigid transformation parameters using vector geometry from hybrid geometric features (including points, lines, and planes) was proposed in Reference [52]. The method does not require an initial alignment, and the estimated transformation parameters are at the same level of quality when using the classic least squares approach, but with an improved computational performance in real-field applications. A fully automatic registration procedure between mobile mapping and nadir airborne imagery to improve the localization quality of mobile mapping platforms in GNSS-denied areas was proposed in Reference [53]. They used common and clearly distinguishable elements, such as road markings, manholes, or curbstones, as conjugate features. The same authors expanded their procedures by including oblique imagery for the registration of both datasets [54]. They extracted façades in sparse point clouds created from terrestrial image correspondence and then used them as projection surfaces for both image datasets. The registration was conducted by a template matching approach. Besides geometry, semantic information about 3D features can contribute to the efficiency of point cloud registration. The registration of terrestrial point clouds, which uses semantic information to automate geometric feature matching, was proposed in Reference [55]. The semantic-aided feature matching concept leads to a more robust and efficient registration process. Yang et al. [26] proposed a marker-free and multiview registration method for large-scale urban scene point clouds based on semantic feature points. The main contribution of the proposed method was the detection of a small set of semantic feature points, matching them using both geometrical constraints (three-point scheme) and their semantic information (category and direction). The proposed method worked well even in cases where the overlap was only 15% between neighboring scans.
Besides providing the initial parameters for fine registration, coarse alignment is also important when there is not enough overlap between the registered point clouds. In fine registration, the main objective is to achieve the maximum overlap of two point clouds [12]. Even feature-based methods require overlap between point clouds: For example, at least three pairs of conjugate surfaces should be present in point clouds to be registered [26]. Problems with registration can occur when using multitemporal 3D point clouds obtained from data acquired at different times. The surface may change over time, and the point clouds might not be comparable. To overcome this problem, a multisensor and multidata integration approach performed on a common primary data level was presented in Reference [56]. Targets can be the appropriate primitives for the registration of datasets that are acquired at different times only if one can guarantee that the targets are set up at the same locations in all surveys [1]. An automatic technique to register multitemporal datasets from UAVs was presented in Reference [57]. The authors proposed a reference image block coregistration algorithm that uses reference anchor images to constrain the orientation of images acquired at different times. In order to produce good results for the registration, the anchor images must be well distributed and not include the changing area.

1.4. A New Target for TLS and UAV Image Point Cloud Registration

In this article, we focus on the registration of point clouds measured with TLS and produced from photogrammetric images taken from UAVs. Two different approaches can be adopted for the registration of 2D image-based and 3D laser point clouds [1]. The first approach is based on a 2D-to-3D alignment strategy, in which registration is achieved using identified features in the 2D images and 3D point cloud in question [58]. The second approach is based on a 3D-to-3D registration strategy. More specifically, a point cloud can be generated from 2D imagery through the automated identification of conjugate points in overlapping imagery. Then, the registration problem is solved through the alignment of image-based and laser scan point clouds [1].
Many examples that combine TLS and UAV image point clouds can be found in the scientific literature [11,13,19,21,23,24,59,60,61,62,63]. By reviewing that literature, we discovered that only a few studies have used special targets to register point clouds. To our knowledge, researchers have mostly relied on feature-based methods or fine registration algorithms or have registered point clouds merely on GCPs that were signalized by targets that were separated for both technologies. If common targets are used for TLS and UAV point clouds, their shape is planar.
Our aim was to design a target that would be appropriate for engineering tasks in smaller areas, where we cannot achieve high overlap between point clouds or there are not enough objects to successfully perform feature-based or fine-registration methods. The positional accuracy of TLS is usually much higher than that of SfM [23], and the resolution of the registered point clouds can vary. Therefore, the target should enable us to precisely measure the reference target point in TLS and UAV image datasets. In an innovative approach, we constructed a unique target that combines the properties of plane targets, which are usually used for the signalization of GCPs in photogrammetric images, and volume targets, which are often used to register laser scans. The orientation of planes whose normal vectors cover all main directions in space contributes to the successful registration of point clouds [55,64]. Hence, the “new target” (as we name it in this article) is composed of three perpendicular planes. Thus, we can describe it as a “three-plane target”. In an extensive experimental study, we compared and analyzed the quality of aerial and terrestrial point clouds that were registered using different commonly used plane and volume targets as well as the new target, hoping to establish whether the new target can yield better results.

2. Materials and Methods

2.1. New Target

In our study, we conceived and created a target that can be used for georeferencing and registering point clouds that are produced from UAV images and TLS point clouds. It comprises two parts: A normal UAV target (a black circle on a white background) and an additional upper part, which is used for TLS (Figure 1a). This target enables the determination of the reference target point in both datasets. The target is composed of 3 square aluminum panels, each measuring 40 cm × 40 cm × 0.2 cm, and a connecting element in the form of a discus with a diameter of 10 cm and a width of 0.8 cm. All panels are painted matte white, as this provides good reflection with its coarseness and shades.
In the event of photogrammetric surveying with a UAV, only the lower part of the target is used, i.e., the square panel, with a circle 30 cm in diameter marked on it with black matte vinyl foil. The reference target point is defined by the center of the circle.
When using TLS, the upper part of the target is attached. The rotating connecting element is attached to the UAV target, and two vertical, interconnected rectangular panels are attached to the connecting element (Figure 1b). The reference target point is defined by the cross-section of all three panels, which physically coincides with the center of the black circle. The alignment of the cross-section of the vertical panels with the center of the black circle on the lower part of the target is ensured by the connecting element, which has a rotating axis on the lower plane. This axis is attached to the hole in the center of the black circle in the horizontal panel. The upper plane of the connecting element has two rectangular slots that serve to stabilize and centralize the vertical panels. The vertical panels have circles with diameters of 15 cm cut into them, which enables us to measure the points on all four sections of the vertical panels while scanning and thus provides us with a better definition of the planes. The rotating connecting element enables us to turn the vertical panels while scanning so that the holes are turned toward the scanner. This improves the quality of the scanned points on the target, because the incidence angles on both vertical panels are lower and almost identical.
The dimensions of the target are primarily determined by the size of the black circle, which is often used for signalizing GCPs in photogrammetric UAV surveys. The circle diameter of 30 cm ensures that the target will be imaged with a sufficient number of pixels to perform automatic image measurements for the surveys at up to 100 m of altitude (depending on the camera’s focal distance). We limited the dimensions of the holes in the vertical plates to 15 cm in order to ensure the robust construction of the target. We scanned the target in different positions from the 10-m distance with a resolution of 1 mm and established that sufficient scan points would be produced at more distant sections of the vertical panels. A larger target would be more appropriate for the TLS point cloud registration, but the chosen dimensions of the target enable its easy transportation.

2.2. Determining the Coordinates of the Reference Target Points

When working with plane targets, the reference target point is the center of the circle. In spheres, the reference target point is the center of the sphere, whereas in cones, it is represented by the apex of the cone. In our experiment (Section 2.5), the coordinates of the reference target points were determined in all three datasets (referential measurements taken with a Leica Nova MS50 MultiStation (henceforth MS50), the UAV images or point clouds, and the TLS point cloud) and were later used to register TLS and UAV point clouds. Prior to determining the coordinates of the reference target points in the point clouds, we manually cut out the point cloud of the individual target.

2.2.1. Determining the Coordinates of the Reference Point for the New Target

The reference coordinates of the new target were defined with a classic polar measurement with MS50 and the use of a miniprism in both faces.
Image matching is used to measure the coordinates of the reference target points on images in an aerial photogrammetric block: This can be performed, for instance, with least-squares matching [65]. We used Agisoft PhotoScan software, v. 1.4.1 [66], to measure the image coordinates of the reference target points. The spatial coordinates of the reference target point were calculated from the known parameters of the interior and exterior orientation of the camera, obtained in the SfM process.
In the TLS point cloud, the reference target point is defined by the cross-section of the three planes, one horizontal and two vertical. We created an automatic algorithm for calculating the coordinates of the reference target point. The algorithm is presented in Figure 2 and was practically realized in programs we wrote in the MATLAB environment.
The input data are represented by a cut-out point cloud of the target and its vicinity (Figure 3). With the primary segmentation performed using RANSAC [67], we located the points within the input point cloud that belonged to the three planes (Figure 4). Following primary segmentation, we eliminated the systematic errors that appeared on the horizontal plane (description in the following paragraph) and performed secondary RANSAC segmentation. In order to find three different planes, RANSAC was performed repeatedly. After the first plane was segmented, the points that belonged to that plane were removed from the point cloud. The process was repeated until the third plane and its points were located. The parameters used in the process of primary and secondary RANSAC segmentation are presented in Figure 2. Parameters are denoted as follows: p is the probability that inliers will be randomly chosen in at least one of the iterations; u is the expected percentage of inliers in a point cloud; m is the minimum number of points that uniquely define the model; t is the threshold, with points within it being considered inliers; and N is the number of iterations.
Figure 4 shows that only points on the white part of the lower plane are segmented in primary RANSAC segmentation. The target’s lower panel is scanned at a large incidence angle. Several authors have studied the influence of the incidence angle and the reflectivity of a material on reflectorless TLS distance measurements. The noise level increases with higher incidence angles and range, which in turn influences the precision of measuring the distance [68]. The offset of the distance measurement is also connected to the target reflection properties [32,34] and the algorithms used to calculate the distance, which may vary depending on the type of scanner [69]. Some authors have established that materials with lower reflectivity produce more noise [70,71]. In our case, we noticed elongated distances on the black circle of the horizontal panel. The effect is shown in Figure 5 (in order to create the figure, one of the targets was rotated to a vertical position and colored by height). In order to reduce the influence of the described errors on the parameters of the horizontal plane, we excluded the points that lay within a 15.5-cm radius from the approximate position of the reference target point (the radius of the black circle is 15 cm) of the horizontal plane point cloud obtained through primary segmentation.
This was followed by secondary segmentation of every plane segment obtained by primary segmentation (separately), this time using a lower tolerance (2 mm) for the point distances from the plane in RANSAC. Figure 6 shows that this eliminated a few points around the intersection of the vertical planes, which were most likely burdened by the multipath.
The plane parameters are defined by minimizing the orthogonal distances of points in the plane. Due to the thickness of the panels and the manner in which the target was created, the cross-section of the vertical planes does not lie on the reference target point. In order to calculate the exact position of the reference target point, we needed to move the vertical planes by half of the thickness of the panel (1 mm) in the direction of the normal vector of the plane, away from the scanner. As a last step, the coordinates of the cross-section of the three planes were calculated.
Within the TLS data, the coordinates of the reference point of the new target were defined in the point cloud of each scanner station individually. The reason behind this can be found in the way the target is used and the procedure for defining the reference target point. As the upper part of the target is always turned toward the scanner station, the registered point cloud from all scanner stations shows a group of points on the vertical planes in various positions above the horizontal plane, which does not permit the use of the developed algorithm when determining the reference target point. The final coordinates of the reference point for the new target were calculated as the mean values, obtained from georeferenced point clouds from individual scanner stations. The precision of defining the final coordinates of the reference point of the new target is also burdened by the quality of registration of the TLS point clouds, which is why this was estimated with the use of the standard deviation [72] of the coordinates of the reference target point from n scanner stations.

2.2.2. Determining the Coordinates of the Reference Target Points on the Remaining Used Targets

In this section, we refer to the targets that we used in our experimental study. The details of those targets are given in Section 2.5.
The coordinates of the reference point for the retroreflective target were defined through referential measurements with an MS50 and in the same way as for the new target in the UAV point cloud (Section 2.2.1). We determined the coordinates of the reference point for the retroreflective target in the TLS point cloud using the procedure described in Reference [73]. This procedure stipulates that the coordinates need to be determined by an evaluation of the precision of the coordinates by error propagation law.
The coordinates of the reference target points of all volume targets, both spheres and the cone, were defined with the least squares method, which is described in Reference [74]. The parameters of the geometric forms were obtained by solving the mathematical model. In the stochastic part of the adjustment, the covariance matrix was used to calculate the precision of all parameters. In order to calculate the coordinates and the precision of the coordinates of the reference target points, we wrote new scripts in the MATLAB environment.

2.3. Point Cloud Registration

The point cloud registration was performed with a Helmert transformation of the point cloud from the source coordinate system into the target coordinate system of a different point cloud. Equation (1) of the transformation is defined by 7 parameters:
  • 3 translations in the directions of the coordinate axes: Δ x , Δ y , and Δ z ;
  • 3 rotations around the coordinate axes: ω x , ω y , and ω z ;
  • a change in the scale: m .
Equation (1) is
[ x t r a n s y t r a n s z t r a n s ] = m · R · [ x s o u r c e y s o u r c e z s o u r c e ] + [ Δ x Δ y Δ z ] ,
where ( x s o u r c e ,   y s o u r c e ,   z s o u r c e ) are the coordinates in the source coordinate system, and ( x t r a n s ,   y t r a n s ,   z t r a n s ) are the transformed coordinates in the target coordinate system. The rotations around the coordinate axes form the rotation matrix R. The procedure to solve Equation (1) is described in Reference [43] (pp. 14–18).
We calculated the parameters of the transformation on the basis of the coordinates of the reference target points of the individual tie points (the parameters were calculated separately for each tie point type). We used the specified transformation parameters to transform the check points—signalized with the retroreflective targets—into the target coordinate system.
The following transformations from the “source” to the “target” point cloud were performed:
  • Transformation of the UAV point cloud into the TLS point cloud;
  • Transformation of the UAV point cloud into referential data, measured with MS50;
  • Transformation of the TLS point cloud into referential data, measured with MS50.
The calculation of the transformation parameters and the data transformation were performed with the scripts we wrote in the MATLAB environment.

2.4. Quality of Data Registration

The quality of the registration was evaluated on the basis of the differences between the coordinates of the retroreflective check points in the target coordinate system and the transformed coordinates of the check points. As a measure of accuracy, we used the root mean square error (RMSE), which is often used as a measure for spatial data accuracy [75].
The coordinate differences ( Δ x , Δ y , Δ z ) on the check points were calculated with Equation (2):
[ Δ x Δ y Δ z ] = [ x t a r g e t y t a r g e t z t a r g e t ] [ x t r a n s y t r a n s z t r a n s ] ,
where ( x t a r g e t ,   y t a r g e t ,   z t a r g e t ) are the coordinates of the check point in the target coordinate system, and ( x t r a n s ,   y t r a n s ,   z t r a n s ) are the coordinates of the check point from the second dataset transformed into the target coordinate system.
Equations (3) and (4) and the coordinate differences were used to calculate the 2D and 3D positional differences Δ:
Δ 2 D = Δ x 2 + Δ y 2 ,
Δ 3 D = Δ x 2 + Δ y 2 + Δ z 2 .
The RMSE was calculated by coordinate axes x, y, z and their 2D and 3D positions with Equations (5)–(7):
R M S E x = i = 1 n Δ x i 2 n ,   R M S E y = i = 1 n Δ y i 2 n ,   R M S E z = i = 1 n Δ z i 2 n
R M S E 2 D = i = 1 n Δ 2 D i 2 n = R M S E x 2 + R M S E y 2 ,
R M S E 3 D = i = 1 n Δ 3 D i 2 n = R M S E x 2 + R M S E y 2 + R M S E z 2 ,
where n represents the number of check points.

2.5. Experimental Setup

In order to analyze the appropriateness of artificial targets for the registration of TLS and UAV point clouds, we created a test area (Section 2.5.1) and scanned it with TLS and photographed it with a UAV. We placed artificial targets of various shapes and sizes throughout the test area. The first types of targets (Figure 7a–c) were used for georeferencing TLS and UAV point clouds:
  • Riegl cylindrical retroreflector with a diameter of 10 cm for registering TLS point clouds from different scanner stations;
  • Leica HDS 6′’ Tilt & Turn target used as a GCP for georeferencing TLS point clouds;
  • Photogrammetric GCP (a black circle with a diameter of 27 cm on a white background) for georeferencing UAV aerial images.
The second types of targets were volume targets (Figure 8), which were used as tie points for registering TLS and UAV point clouds:
  • Small sphere, diameter of 15 cm;
  • Large sphere, diameter of 20 cm;
  • Right circular cone, diameter of 40 cm;
  • New target.
The spheres and cones were painted matte white, and then various colored stickers were attached to the bodies in a random pattern. This was performed so that the featured points could be better identified with the SfM algorithms.
We evaluated the quality of registering the point clouds on check points made from retroreflective targets, which were created especially for this study (Figure 7d). The target is in the form of a black plane, 40 cm × 40 cm × 0.5 cm, in the center of which a circle made from retroreflective foil, measuring 30 cm in diameter, is stuck. The retroreflective target enables precise measurements in TLS and UAV point clouds.
The precise measurements of all tie point targets, GCPs, and check points were performed with an MS50, which enables laser scanning and also classic geodetic measurements.

2.5.1. Field Measurements and Initial Data Processing

Within the test area, which measures approximately 55 m × 55 m, we stabilized GCPs, check points, and tie point targets in five groups (Figure 9 shows groups G1–G5, marked by ellipses). We analyzed their appropriateness for registration of the point clouds in our research. Each group included the following targets:
  • Tie points (small and large sphere, cone, and new target);
  • Photogrammetric GCP;
  • Cylindrical retroreflector for registering TLS point clouds;
  • Retroreflective check point target.
The groups were organized throughout the test area so that one group of targets was in every corner of the test area and one was in the center of the area (Figure 9). Apart from the targets that were stabilized within the groups, an additional 4 photogrammetric GCPs for georeferencing the UAV point cloud and 4 Tilt & Turn targets for georeferencing the TLS point cloud were stabilized outside of the groups.
Target stabilization within a group is shown in Figure 10. The retroreflective check point targets were raised by approximately 10 cm on the outer side of the area. This reduced the incidence angle of the TLS measurement rays and thus increased the quality and number of measured points on the target without deforming the shape of the circle on the aerial images.

2.5.2. Referential Measurements

We used the MS50 measurements to define the local coordinate system, which we treated as the reference coordinate system. The coordinates of the centers of the plane targets (new targets, retroreflective check point targets, photogrammetric GCPs, and Tilt & Turn targets) were defined with the classic polar measurement method on both faces by signalizing the centers of the targets with a miniprism or, in the case of the Tilt & Turn target, with direct sighting of the center of the target. When we were determining the coordinates, we took into account the meteorological parameters (temperature and air pressure), which we measured with a Meteo Station HM30, and entered the data into the MS50. The volume targets (small spheres, large spheres, and cones) were scanned at a resolution of 2.5 mm × 2.5 mm. The MS50 enables measurements on the prism (single measurement) with an angle and distance accuracy of 1” and 1 mm + 15 ppm, respectively. The declared accuracy and precision of scanning in the 1000-Hz mode is 1 mm at 50 m.

2.5.3. Terrestrial Laser Scanning and Processing of the TLS Point Cloud

The laser scanning was performed with a Riegl terrestrial laser scanner VZ-400, which was controlled with the RiSCAN PRO software. The scanner has a declared scanning precision and accuracy of 3 mm and 5 mm, respectively. We performed the scanning with 4 stations (Figure 9) positioned throughout the test area. We created an overview scan on every scanner station and a detailed scan of the tie points and retroreflective check points with a resolution of 1 mm × 1 mm on the target. The software has its own way of scanning TLS tie targets (cylindrical retroreflector and flat Tilt & Turn target) and a built-in algorithm, with which it calculates the coordinates of the tie points in the coordinate system of an individual scanner station. Based on the coordinates of the cylindrical retroreflectors, we registered all four TLS scans with RiSCAN PRO software. The registered TLS point cloud was georeferenced into the reference coordinate system with the coordinates of the Tilt & Turn targets, which were measured with the MS50. The registration precision of the scanner stations is shown with the standard deviation of the registered coordinates of the tie points (Table 1). The TLS point cloud georeferencing precision was 3.6 mm.

2.5.4. Photogrammetric Surveying and the Creation of UAV Point Clouds

A UAV DJI Phantom 4 PRO with an FC6310 camera with a focal length of 8.8 mm was used for photogrammetric surveying. In the UAV point cloud, the number of points on an individual target depended on the size of the target and the altitude of the flight, which at the given focal length and size of the image sensor determined the ground sample distance (GSD) of the aerial images. In order to check the suitability of the individual artificial targets at various GSDs, we performed three surveys with the UAV at various altitudes, which gave us the following approximate average GSDs of aerial images: 0.5, 1.0, and 2.0 cm. Only nadir images were taken in all three surveys.
We used the obtained images to create three photogrammetric image blocks (UAV20, UAV40, and UAV75), which were in turn used to create UAV point clouds. The images were recorded with an 80% longitudinal and 70% transversal image overlap. The photogrammetric block UAV20 is comprised of images recorded in one flight direction, at an altitude of 20 m above ground level. The blocks UAV40 (Figure 11) and UAV70 comprise images recorded from two perpendicular flight directions, which differ in their altitudes by between 6 and 8 m. This improved the inner geometry of the block of aerial photographs in bundle block adjustment, which in turn improved the vertical accuracy of the model [76]. The recording characteristics of the individual block are presented in Table 2.
The individual blocks were photogrammetrically processed with PhotoScan, which uses SfM to create dense point clouds (Figure 12). The external orientation of the images and the georeferenced dense point clouds were calculated with the coordinates of the photogrammetric GCPs measured by the MS50. The quality of the georeferencing was assessed with the check points, i.e., retroreflective targets and new targets. The accuracy of the georeferencing UAV point clouds is shown in Table 3, in which the RMSE is given in cm. The results show that the R M S E 2 D on the check points is within 1 cm for all missions, whereas the R M S E 3 D is within 2 cm, which we believe to be a good result. However, what we can also see is that the RMSE does not necessarily increase with a higher flight altitude. The reason for this is that the RMSE value depends on a number of input factors (terrain morphology, camera calibration, number and distribution of GCPs, image block composition, etc.) and not only on the flight altitude. This issue was thoroughly investigated and explained in Reference [77].

3. Results

Section 3.1 presents the precision of defining the coordinates of the reference target points of the artificial targets, which we used as tie points for registering TLS and UAV point clouds or as check points for evaluating the performed registration. Section 3.2 describes the quality of registration of the UAV and TLS point clouds, taking into account the types of tie points used.

3.1. Precision of Coordinates of Reference Target Points

Table 4 and Table 5 provide 2D and 3D precisions, respectively, of the coordinates of the reference target points for the referential measurements with MS50 for the TLS point cloud and three UAV point clouds (UAV20, UAV40, UAV75). The precision of the coordinates in the tables is defined by the standard deviation, provided in cm. The exceptions are the accuracies of the coordinates of the reference points for the new and retroreflective targets in the UAV point clouds. These are provided in cm with RMSE, calculated from the differences between the coordinates of the referential measurements with MS50 and the coordinates of these points, estimated by aerial triangulation. In Table 4 and Table 5, we marked tie points with the following abbreviations: “sSph” for small sphere, “LSph” for large sphere, and “Cone” and “newT” for the new target. The retroreflective check point target is marked by the abbreviation “rRef”.
As expected, Table 4 and Table 5 show that the reference target points were defined with the greatest precision for the dataset of the referential measurements with MS50. Similar precision for the coordinates of the reference points for volume targets was achieved in the TLS dataset. Somewhat poorer, in the range of just under 5 mm, were the reference points of the new and retroreflective targets defined in the TLS dataset. When compared to the precision of the coordinates of the remaining reference target points, the coordinates of the centers of the spheres were much more precisely defined in all datasets. The reason behind this lay in the use of RANSAC for the segmentation of the point cloud of the sphere points, where, by selecting the tolerance value, we eliminated the blunders. The estimated precision of the calculated coordinates of the centers of the spheres was thus overrated. The centers of the large spheres were defined with the greatest precision in all UAV point clouds, regardless of the spatial resolution. The importance of the spatial resolution of the images is shown in the poor definition of the reference target points in aerial photogrammetric datasets. At a GSD of 1 cm and 2 cm (UAV40 and UAV75), the worst defined were the apexes of the cones, whereas at a GSD of 0.5 cm (UAV20), the plane targets were the worst defined. Roughly speaking, we could say that the volume targets were more precisely defined than the plane targets when the highest spatial resolution was used, whereas at a lower spatial resolution the situation was reversed (if we excluded the sphere targets, whose estimated precision was overrated). In the UAV75 point cloud produced from images with a spatial resolution of 2 cm, the small sphere with a diameter of 15 cm failed to provide an appropriate number of quality points to determine the center of the target. Three of the five targets revealed merely 4 or 7 points.

3.2. Quality of Registration

The quality of registering TLS and UAV point clouds was estimated on the basis of the RMSE of check points, signaled with a retroreflective target. As a comparison, Section 3.2.1 presents the RMSE prior to the performed registrations with the analyzed tie points.

3.2.1. Discrepancies of TLS and UAV Point Clouds on the Check Points Prior to Registration

The UAV point clouds were georeferenced into the reference coordinate system on the basis of photogrammetric GCPs. The TLS point cloud was georeferenced with the Tilt & Turn targets. Both types of georeferencing points were measured with the MS50. Table 6 shows the RMSE calculated from the coordinate differences of the retroreflective check points, defined in UAV and TLS point clouds. Figure 13 and Figure 14 show RMSEs calculated prior to the registration of TLS and UAV point clouds with the marking “Georef”.

3.2.2. Discrepancies of the TLS and UAV Point Clouds on the Check Points after the Registration

All UAV point clouds (source coordinate system) were transformed into the TLS point cloud (target coordinate system) as described in Section 2.3. In order to define the transformation parameters, a different type of tie point was used each time. We performed 12 transformations, the qualities of which are represented by R M S E 2 D and R M S E 3 D in Table 7 and Table 8, respectively. The graphic representations of R M S E 2 D and R M S E 3 D are shown in Figure 13 and Figure 14, respectively. In Table 7, Table 8, Table 9 and Table 10, tie points are marked with the following abbreviations: “sSph” for small sphere, “LSph” for large sphere, and “Cone” and “newT” for the new target. The retroreflective check point target is marked by the abbreviation “rRef”.
Figure 14 shows that the best 3D accuracy in all UAV point clouds was achieved with the use of the new target. When using aerial datasets with spatial resolutions between 0.5 cm and 2 cm, the 3D accuracy of the new target was within 1 cm. The remaining volume targets all provided a 1-cm spatial accuracy ( R M S E 2 D ) when using aerial images with a GSD of at least 1 cm (UAV20 and UAV40; see Figure 13). The R M S E 3 D shows that when volume targets were used, the registration by height tended to be problematic. The worst performance was delivered by the small sphere, which was expected due to its size. Poorer quality 3D registration of the TLS and UAV point clouds was also delivered when the cone was used as a tie point.
When the obtained accuracies are compared to the accuracies from previous georeferencing attempts (Table 6), it is clear that the registrations with all types of analyzed targets at all GSDs of aerial images provided a worse or comparable accuracy, with the exception of the new target, which provided higher accuracy with all GSDs of aerial images.
The suitability of the new target for registering TLS and UAV point clouds is graphically presented in Figure 15, which shows the overlap of the TLS and UAV40 point clouds once georeferencing was performed (prior to registration, see Figure 15a) and after the transformation of the UAV40 point cloud into the TLS point cloud with the use of the new target (Figure 15b). Figure 15 shows that georeferencing with the use of GCPs, measured in the same coordinate system, did not ensure appropriate alignment of the point clouds obtained with different types of technology, even when the GCPs were defined with great precision. However, we successfully aligned TLS and UAV point clouds with the use of the new target as a tie point.

3.2.3. Registration of TLS and UAV Point Clouds on Referential Measurements

The coordinates of the retroreflective check point targets and the new targets were better defined in the dataset of the referential measurements with MS50 than they were in the TLS dataset (Table 4 and Table 5). With the transformation of the UAV and TLS point clouds into referential measurements (target coordinate system), performed with MS50, we checked whether better results of registration could be obtained as a result of the higher precision of the coordinates of the reference target points in the referential measurements data. We performed 16 registrations. In order to define the transformation parameters, a different type of tie point was used each time. The quality of the registrations is represented by R M S E 2 D and R M S E 3 D in Table 9 and Table 10, respectively. The graphic representations of R M S E 2 D and R M S E 3 D are shown in Figure 16 and Figure 17, respectively. We calculated the RMSE from the coordinate differences of the retroreflective check point targets (measured with MS50) and the transformed coordinates of the retroreflective targets in TLS and UAV point clouds. In Figure 16 and Figure 17, RMSEs calculated prior to the registration of TLS and UAV point clouds (Table 6) are shown with the marking “Georef”.
Once again, the new target yielded the best results in all datasets when registering UAV point clouds with referential measurements. With the use of the new target, the R M S E 3 D of the check points was below 1 cm, while the R M S E 2 D ranged between 0.24 and 0.36 cm. The greatest discrepancies on the check points could be found in the results obtained with the small sphere. R M S E 3 D , when registering with the large sphere, ranged between 1.73 cm and 4.92 cm, and when the cone was used, between 2.03 and 6.68 cm. Within the UAV point clouds for all volume targets (spheres and cones), defining the z coordinate was the most problematic, which resulted in a much higher R M S E 3 D compared to R M S E 2 D . When used with UAV point clouds, the new target was the only target, which in all cases yielded better results compared to the results of prior georeferencing. The new target yielded the most accurate results ( R M S E 3 D = 0.49 cm) when registering the TLS point cloud with referential measurements. Similar accuracy was achieved with the small and large spheres.

4. Discussion and Conclusions

This article analyzes the use of various artificial targets as tie points for registering TLS and UAV image point clouds. We used a small and large sphere, a cone, and a new “three-plane target” created especially for this study. The analysis was performed on data obtained by terrain measurements with a terrestrial laser scanner and UAV imagery. Surveys using UAV were performed at three different flight altitudes, so that we could define the use of various tie points in relation to the GSD of the images, from which the UAV point clouds were created.
We calculated the coordinates of the reference target points in both datasets for each tie point and estimated their accuracy. For our new target, we developed a new procedure for defining the coordinates of the reference target point in the TLS point cloud. The point clouds were registered with the use of seven-parameter transformation. The transformation parameters were calculated from the coordinates of the reference target points in TLS and UAV point clouds. We performed numerous registrations of TLS and UAV point clouds for every type of tie point individually. The quality of the performed registrations was estimated using the R M S E of the check points.
Artificial targets are useful as tie points for registering two point clouds when both point clouds enable precise definition of the coordinates of the reference target point. When point clouds from the same source are merged, we can usually ensure that the accuracy and precision of the coordinates of the tie points in both point clouds are similar. When registering point clouds from different sources, as was the case in our example, the accuracy and the precision of the coordinates of the tie points differ. The precision of the coordinates of all analyzed tie points in the TLS point cloud ranged between 1 mm and 4.5 mm, whereas in the UAV point cloud it ranged between 1 mm and 52 mm (Table 4 and Table 5). We can say that the coordinates of the tie points in UAV point clouds were roughly 10 times less precise. The worst precision was calculated for the z coordinate, which has also been ascertained by other authors [21,23,59]. Through a visual overview of the UAV point cloud, which we have not described yet in this article, we ascertained that the volume targets (spheres and cones) yielded a lot of noise, which distorted the geometry of the targets. Most likely, the reason for this lay in overexposed aerial images, as the white volume targets were brighter than their surroundings. Consequently, the volume targets in the created UAV point clouds were slightly flattened (Figure 18). A critical review and analysis of selected dense image-matching algorithms was presented by Remondino et al. [3]. The datasets they adopted for testing included terrestrial and aerial image blocks acquired with convergent and normal (parallel axes) images at different scales and resolution. The authors reported several reasons for problems in dense matching algorithms. These reasons were initial image quality (noise, low radiometric quality, shadows, etc.); poor image configuration; certain surface materials (shiny or textureless objects); or scene/object characteristics, such as the presence of shadows, sharp discontinuities, and small structures. This could result in noisy point clouds, smoothing effects, and/or difficulties in feature extraction. The fact that varying the illumination of images also affects image tie point matching was reported by Gerke, Nex, and Jende [6] in their representation of ISPRS benchmark results of multiplatform photogrammetry.
The precision of the coordinates of the reference target points in the UAV point cloud was dependent on the GSD of the images (Table 5) and their photographic quality, which directly influenced the geometrical quality of the created UAV point cloud. The image GSD could be improved with a lower flight altitude, which extends the time for terrain measurement and the processing of images with SfM [19]. The second possibility for improving the point cloud of the volume tie points is to use larger targets. By increasing the size of the targets, we would get more points on the target and the possibility of using the targets from greater flight altitudes. The downside of the larger-volume targets is that they are harder to transport and are less flexible to work with.
Using a sphere target measuring 20 cm in diameter, we could register the TLS point cloud and the UAV point cloud made from images with a GSD of 1 cm or less, with an accuracy of approximately 2.5 cm (Table 8). The small sphere, measuring 15 cm in diameter, yielded approximately three times worse accuracy. When a UAV point cloud is created from images with a GSD of approximately 1 cm, we suggest that a sphere target with a diameter of at least 20 cm be used for the registration of the TLS and UAV point clouds. The cone-shaped target yielded poorer results when registering point clouds UAV20 and UAV75 with the TLS point cloud ( R M S E 3 D ranging from 3.7 cm to 6 cm). However, it provided better results when registering the UAV40 point cloud (GSD of the images ≈1cm), where the accuracy was estimated with an R M S E 3 D of 1.5 cm.
The new target yielded the best 3D accuracy as a tie point for all UAV point clouds (GSD of the images ranging from 0.5 cm to 2 cm). The R M S E 3 D did not surpass 1.1 cm. The poorer quality of the UAV point clouds on the tie points did not influence the precision of the coordinates of the reference point for the new target. The reason for this lay in the way the coordinates are defined in the UAV dataset. They are calculated with the use of SfM from the image coordinates, which are obtained automatically with least squares matching. In our case, the image measurements enabled a more precise calculation of the coordinates of the reference target point than the modeling of the UAV point clouds did. With the latter, we calculated the reference target points for volume targets. In the TLS point cloud, the reference point of the new target is defined by the cross-section of the three planes, which coincides with the center of the black circle, which represents the reference target point in the UAV dataset. When compared to the volume targets, the precision of the coordinates of the reference target point in the TLS point cloud was slightly lower, but still within a few millimeters. The geometrical relation between the scanner and the target represents the main influence, as the incidence angle of the laser beam on the horizontal plane of the target can be very high when measuring from long distances. The achieved results confirm the good design of the target and its suitability for registering TLS and UAV point clouds whenever a high-quality registration is needed.
By analyzing the registration results, one can see that the RMSE of the registered point clouds was lower than the measurement precision of the new target. To assess the quality of the new target position, its 2D and 3D precision were calculated in all datasets. The new target was provided by the reference coordinates by measuring with the MS50 in both faces using a miniprism. The precision of the reference coordinates was assessed, taking into account the accuracy of the MS50 and using error propagation law. In this case, the new target coordinates were determined with an R M S E 3 D of 1.7 mm (Table 5). In the TLS point clouds, the new target positions were determined in each scan station individually. Due to the new target geometry, which allows for the rotation of vertical panels with holes toward the scanner, the coordinates of the new target were defined as the arithmetic mean of the quadruple-calculated coordinates. This was also related to the procedure for determining the precision of a new target in the TLS point cloud. We calculated it as a standard deviation of four-times-determined coordinates. The R M S E 3 D of the new target in the TLS point clouds was 4.3 mm (Table 5) and was related to the quality of the scan stations’ registration with Riegl cylindrical retroreflectors. The positions of new targets and their RMSEs in the UAV point cloud were determined using the automatic procedure in PhotoScan. The precisions were between 11.5 mm and 20.3 mm (Table 5) for three different flight altitudes. With respect to the measurement equipment, the precision of the new target positions in the various datasets was within the expected values. When we used the new target to register the UAV and TLS point cloud at the check points, we received an R M S E 3 D less than 7 mm in all registrations (Table 10). The reason for the R M S E 3 D being up to three times lower at check points than the R M S E 3 D of the new target positions in the UAV point clouds was the calculation mode. The R M S E 3 D was calculated from the differences in the positions of the check points in the target coordinate system after the transformation. The precision of the new target was therefore not included in the calculation of the quality of the registration.
One of the main shortcomings of the new target is that it cannot be used from further distances from the scanner. In that case, the incidence angles are too large, which affects the measurement accuracy of the horizontal plane of the new target. The incidence angle and the range affect the individual point signal-to-noise ratio, and a larger standard deviation of errors in the direction of the laser beam was observed for larger incidence angles in Reference [68]. Lichti [78] suggests using the a priori threshold of a maximum incidence angle of 65° for removing nonreliable measurements. Therefore, careful planning for TLS measurements is essential when using the new target. The targets can be set on the border of the scanning area and are preferably tilted toward the scan stations to improve the scanning geometry. Future work will be performed to model the range and incidence noise with the methods described in Reference [68]. Even if we can model the noise, we must prevent the horizontal plates of the targets from being perpendicular to the incident rays of the scanner, which is sometimes difficult to control in practical settings. In these cases, the calculation of the reference target point’s coordinates with a cross-section of three planes will not be possible because we will have no points for modeling the horizontal plane of the target.
A systematic error resulting from the new target design is the generation of mixed pixels [79] that appeared in the TLS point cloud, which are clearly visible in Figure 5. This data artifact was caused due to nonzero laser beam width. A pulsed laser scanner without full-waveform recording capability cannot discriminate between the returns from two surfaces separated by less than half the pulse length, so point measurements appeared in the range discontinuity region between the two surfaces [27]. Segmentation of the target planes using RANSAC removed all mixed pixels (Figure 4 and Figure 6). Nevertheless, the influence of the mixed pixels on the reference target point determination should be investigated in future work.
In our experiment, the new target and volume targets were manually cropped from the point clouds before determining the coordinates of the reference target points. This task could possibly be automatized for the new target, as in our case we first used GCPs to coarsely register the point clouds. As the image coordinates of the new targets were automatically measured with PhotoScan, and their 3D reference coordinates were calculated in SfM, we could use these coordinates to obtain conjugate point clouds of the new target in the TLS point cloud. However, due to the small dimensions of the used targets, we found it difficult to search them automatically in the point clouds. In the case of point clouds not being aligned, the automation of the registration process would be a demanding task, a challenge to be addressed and solved in the future.
By analyzing the available data, we ascertained that the image GSD had a relatively small influence on the accuracy of the registration of TLS and UAV point clouds, ranging between 0.5 and 2 cm. In contrast, the type of target that we used as a tie point influenced the quality of registration. The new target provided the highest accuracy in all the analyzed examples. It should be mentioned that better photographic quality of aerial images might yield higher-quality results when working with volume targets. Due to the various restrictions of projects, it is impossible to freely choose the time for terrain measurements in practice, which means that we do not always have optimal weather and light conditions. The new target proved to be the most optimal tie point for the registration of TLS and UAV point clouds in poorer conditions as well. We can also use it as the GCP for georeferencing the point cloud into the desired coordinate system.

Author Contributions

T.U. and D.G. conceived the main research idea; T.U., Ž.R. and D.G. conceived and designed the experiment; T.U. and Ž.R. performed the onsite measurements; Ž.R. analysed the results and wrote the original draft; M.K.F. and D.G. supervised the project and provided scientific guidance throughout the project; T.U., M.K.F. and D.G. revised, expanded and improved the final manuscript.

Funding

The authors acknowledge the financial support of the Slovenian Research Agency (research core funding No. P2-0227 Geoinformation Infrastructure and Sustainable Spatial Development of Slovenia, and No. P2-0406, Earth Observation and Geoinformatics).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, F.; Habib, A. A Closed-Form Solution for Coarse Registration of Point Clouds Using Linear Features. J. Surv. Eng. 2016, 142, 04016006. [Google Scholar] [CrossRef]
  2. Kraus, K. Photogrammetry, Geometry from Images and Laser Scans, 2nd ed.; Walter de Gruyter: Berlin, Germany, 2007; ISBN 978-3-11-019007-6. [Google Scholar] [CrossRef]
  3. Remondino, F.; Spera, M.G.; Nocerino, E.; Menna, F.; Nex, F. State of the art in high density image matching. Photogramm. Rec. 2014, 29, 144–166. [Google Scholar] [CrossRef] [Green Version]
  4. Dahl, A.L.; Aanæs, H.; Pedersen, K.S. Finding the Best Feature Detector-Descriptor Combination. In Proceedings of the 2011 International Conference on 3D Imaging, Modeling, Processing, Visualization and Transmission, Hangzhou, China, 16–19 May 2011; pp. 318–325. [Google Scholar] [CrossRef]
  5. Urban, S.; Weinmann, M. Finding a good feature detector-descriptor combination for the 2D keypoint-based registration of TLS point clouds. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2015, II-3/W5, 121–128. [Google Scholar] [CrossRef]
  6. Gerke, M.; Nex, F.; Jende, P. Co-registration of terrestrial and UAV-based images—Experimental results. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, XL-3/W4, 11–18. [Google Scholar] [CrossRef]
  7. Micheletti, N.; Chandler, J.H.; Lane, S.N. Structure from motion (SFM) photogrammetry. In Geomorphological Techniques (Online Edition); Clarke, L.E., Nield, J.M., Eds.; British Society for Geomorphology: London, UK, 2015; Chapter 2, Section 2.2; ISBN 2047-0371. [Google Scholar]
  8. Piermattei, L.; Carturan, L.; De Blasi, F.; Tarolli, P.; Dalla Fontana, G.; Vettore, A.; Pfeifer, N. Suitability of ground-based SfM–MVS for monitoring glacial and periglacial processes. Earth Surf. Dyn. 2016, 4, 425–443. [Google Scholar] [CrossRef]
  9. Benassi, F.; Dall’Asta, E.; Diotri, F.; Forlani, G.; Morra di Cella, U.; Roncella, R.; Santise, M. Testing Accuracy and Repeatability of UAV Blocks Oriented with GNSS-Supported Aerial Triangulation. Remote Sens. 2017, 9, 172. [Google Scholar] [CrossRef]
  10. Xu, Z.; Wu, L.; Gerke, M.; Wang, R.; Yang, H. Skeletal camera network embedded structure-from-motion for 3D scene reconstruction from UAV images. ISPRS J. Photogramm. Remote Sens. 2016, 121, 113–127. [Google Scholar] [CrossRef]
  11. Persad, R.A.; Armenakis, C. Automatic co-registration of 3D multi-sensor point clouds. ISPRS J. Photogramm. Remote Sens. 2017, 130, 162–186. [Google Scholar] [CrossRef]
  12. Cheng, L.; Chen, S.; Liu, X.; Xu, H.; Wu, Y.; Li, M.; Chen, Y. Registration of Laser Scanning Point Clouds: A Review. Sensors 2018, 18, 1641. [Google Scholar] [CrossRef] [PubMed]
  13. Zang, Y.; Yang, B.; Li, J.; Guan, H. An Accurate TLS and UAV Image Point Clouds Registration Method for Deformation Detection of Chaotic Hillside Areas. Remote Sens. 2019, 11, 647. [Google Scholar] [CrossRef]
  14. Bosché, F. Plane-based registration of construction laser scans with 3D/4D building models. Adv. Eng. Inform. 2012, 26, 90–102. [Google Scholar] [CrossRef]
  15. Habib, A.; Ghanma, M.; Morgan, M.; Al-Ruzouq, R. Photogrammetric and Lidar Data Registration Using Linear Features. Photogramm. Eng. Remote Sens. 2005, 71, 699–707. [Google Scholar] [CrossRef]
  16. Besl, P.J.; McKay, N.D. A method for registration of 3-D shapes. IEEE Trans. Pattern Anal. Mach. Intell. 1992, 14, 239–256. [Google Scholar] [CrossRef]
  17. Yao, J.; Ruggeri, M.R.; Taddei, P.; Sequeira, V. Automatic scan registration using 3D linear and planar features. 3D Res. 2011, 1, 6. [Google Scholar] [CrossRef]
  18. Al-Durgham, K.; Habib, A. Association-Matrix-Based Sample Consensus Approach for Automated Registration of Terrestrial Laser Scans Using Linear Features. Photogramm. Eng. Remote Sens. 2014, 80, 1029–1039. [Google Scholar] [CrossRef]
  19. Xu, Z.; Wu, L.; Shen, Y.; Li, F.; Wang, Q.; Wang, R. Tridimensional Reconstruction Applied to Cultural Heritage with the Use of Camera-Equipped UAV and Terrestrial Laser Scanner. Remote Sens. 2014, 6, 10413–10434. [Google Scholar] [CrossRef] [Green Version]
  20. Zhang, R.; Li, H.; Liu, L.; Wu, M. A G-Super4PCS Registration Method for Photogrammetric and TLS Data in Geology. ISPRS Int. J. Geo-Inf. 2017, 6, 129. [Google Scholar] [CrossRef]
  21. Fugazza, D.; Scaioni, M.; Corti, M.; D’Agata, C.; Azzoni, R.S.; Cernuschi, M.; Smiraglia, C.; Diolaiuti, G.A. Combination of UAV and terrestrial photogrammetry to assess rapid glacier evolution and map glacier hazards. Nat. Hazards Earth Syst. Sci. 2018, 18, 1055–1071. [Google Scholar] [CrossRef] [Green Version]
  22. Liang, H.; Li, W.; Lai, S.; Zhu, L.; Jiang, W.; Zhang, Q. The integration of terrestrial laser scanning and terrestrial and unmanned aerial vehicle digital photogrammetry for the documentation of Chinese classical gardens—A case study of Huanxiu Shanzhuang, Suzhou, China. J. Cult. Herit. 2018, 33, 222–230. [Google Scholar] [CrossRef]
  23. Jo, H.Y.; Hong, S. Three-Dimensional Digital Documentation of Cultural Heritage Site Based on the Convergence of Terrestrial Laser Scanning and Unmanned Aerial Vehicle Photogrammetry. ISPRS Int. J. Geo-Inf. 2019, 8, 53. [Google Scholar] [CrossRef]
  24. Moon, D.; Chung, S.; Kwon, S.; Seo, J.; Shin, J. Comparison and utilization of point cloud generated from photogrammetry and laser scanning: 3D world model for smart heavy equipment planning. Autom. Constr. 2019, 98, 322–331. [Google Scholar] [CrossRef]
  25. Park, J.; Kim, P.; Cho, Y.K.; Kang, J. Framework for automated registration of UAV and UGV point clouds using local features in images. Autom. Constr. 2019, 98, 175–182. [Google Scholar] [CrossRef]
  26. Yang, B.; Dong, Z.; Liang, F.; Liu, Y. Automatic registration of large-scale urban scene point clouds based on semantic feature points. ISPRS J. Photogramm. Remote Sens. 2016, 113, 43–58. [Google Scholar] [CrossRef]
  27. Lichti, D.; Skaloud, J. Registration and calibration. In Airborne and Terrestrial Laser Scanning, 1st ed.; Vosselman, G., Mass, H.G., Eds.; Whittles Publishing: Dunbeath, UK, 2010; pp. 97–98, 111–113. [Google Scholar]
  28. Lichti, D.D.; Harvey, B.R. The effects of reflecting surface properties on time-of-flight laser scanner measurements. In Proceedings of the ISPRS Commision IV Symposium 2002 on Geospatial Theory, Processing and Applications, Ottawa, ON, Canada, 9–12 July 2002; Volume 34. [Google Scholar]
  29. Voegtle, T.; Wakaluk, S. Effects on the measurements of the terrestrial laser scanner HDS 6000 (Leica) caused by different object materials. ISPRS Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2009, XXXVIII-3/W8, 68–74. [Google Scholar]
  30. Lichti, D.D.; Stewart, M.P.; Tsakiri, M.; Snow, A.J. Calibration and testing of a terrestrial laser scanner. ISPRS Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2000, XXXIII-B5, 485–492. [Google Scholar]
  31. Pesci, A.; Teza, G. Terrestrial laser scanner and retro-reflective targets: An experiment for anomalous effects investigation. Int. J. Remote Sens. 2008, 29, 5749–5765. [Google Scholar] [CrossRef]
  32. Schmitz, B.; Holst, C.; Medic, T.; Lichti, D.D.; Kuhlmann, H. How to Efficiently Determine the Range Precision of 3D Terrestrial Laser Scanners. Sensors 2019, 19, 1466. [Google Scholar] [CrossRef]
  33. Pfeifer, N.; Dorninger, P.; Haring, A.; Fan, H. Investigating Terrestrial Laser Scanning Intensity Data: Quality and Functional Relations. In Proceedings of the 8th International Conference on Optical 3D Measurement Techniques, Zurich, Switzerland, 9–13 July 2007; pp. 328–337. [Google Scholar]
  34. Zámečníková, M.; Wieser, A.; Woschitz, H.; Ressl, C. Influence of surface reflectivity on reflectorless electronic distance measurement and terrestrial laser scanning. J. Appl. Geod. 2014, 8, 311–326. [Google Scholar] [CrossRef]
  35. Ge, X.; Wunderlich, T. Target identification in terrestrial laser scanning. Surv. Rev. 2015, 47, 129–140. [Google Scholar] [CrossRef]
  36. Akca, D. Full automatic registration of laser scanner point clouds. In ETH Zurich Research Collection; ETH Zurich: Zurich, Switzerland, 2003. [Google Scholar]
  37. Artese, G.; Achilli, V.; Salemi, G.; Trecroci, A. Automatic Orientation and Merging of Laser Scanner Acquisitions Through Volumetric Targets: Procedure Description and Test Results. ISPRS-Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2004, XXXV-B5, 210–215. [Google Scholar]
  38. Rabbani, T.; Dijkman, S.; Van den Heuvel, F.; Vosselman, G. An integrated approach for modelling and global registration of point clouds. ISPRS J. Photogramm. Remote Sens. 2007, 61, 355–370. [Google Scholar] [CrossRef]
  39. Franaszek, M.; Cheok, G.S.; Witzgall, C. Fast automatic registration of range images from 3D imaging systems using sphere targets. Autom. Constr. 2009, 18, 265–274. [Google Scholar] [CrossRef]
  40. Yun, D.; Kim, S.; Heo, H.; Ko, K.H. Automated registration of multi-view point clouds using sphere targets. Adv. Eng. Inform. 2015, 29, 930–939. [Google Scholar] [CrossRef]
  41. Huang, J.; Wang, Z.; Gao, J.; Huang, Y.; Towers, D.P. High-Precision Registration of Point Clouds Based on Sphere Feature Constraints. Sensors 2017, 17, 72. [Google Scholar] [CrossRef] [PubMed]
  42. Rabbani, T.; Van den Heuvel, E. Automatic point cloud registration using constrained search for corresponding objects. In Proceedings of the 7th Conference on Optical 3-D Measurement Techniques, Vienna, Austria, 3–5 October 2005; pp. 177–186. [Google Scholar]
  43. Vezočnik, R. Analysis of Terrestrial Laser Scanning Technology for Structural Deformation Monitoring. Ph.D. Thesis, University of Ljubljana, Ljubljana, Slovenia, 2011. [Google Scholar]
  44. Von Hansen, W.; Gross, H.; Thoennesen, U. Line-based registration of terrestrial and airborne LiDAR data. ISPRS Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, XXXVII-B3a, 161–166. [Google Scholar]
  45. Yang, B.; Zang, Y. Automated registration of dense terrestrial laser-scanning point clouds using curves. ISPRS J. Photogramm. Remote Sens. 2014, 95, 109–121. [Google Scholar] [CrossRef]
  46. Theiler, P.W.; Schindler, K. Automatic registration of terrestrial laser scanner point clouds using natural planar surfaces. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, 3, 173–178. [Google Scholar] [CrossRef]
  47. Maalek, R.; Lichti, D.D.; Ruwanpura, J.Y. Robust Segmentation of Planar and Linear Features of Terrestrial Laser Scanner Point Clouds Acquired from Construction Sites. Sensors 2018, 18, 819. [Google Scholar] [CrossRef]
  48. Xu, Y.; Boerner, R.; Yao, W.; Hoegner, L.; Stilla, U. Automated coarse registration of point clouds in 3D urban scenes using voxel based plane constraint. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2017, IV-2/W4, 185–191. [Google Scholar] [CrossRef]
  49. Wujanz, D.; Schaller, S.; Gielsdorf, F.; Gründig, L. Plane-based registration of several thousand laser scans on standard hardware. ISPRS-Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2018, XLII-2, 1207–1212. [Google Scholar] [CrossRef]
  50. Gruen, A.; Akca, D. Least squares 3D surface and curve matching. ISPRS J. Photogramm. Remote Sens. 2005, 59, 151–174. [Google Scholar] [CrossRef] [Green Version]
  51. Jaw, J.J.; Chuang, T.Y. Feature-based registration of terrestrial LiDAR point clouds. ISPRS Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, XXXVII-B3b, 303–308. [Google Scholar]
  52. Han, J.Y. A Noniterative Approach for the Quick Alignment of Multistation Unregistered LiDAR Point Clouds. IEEE Geosci. Remote Sens. Lett. 2010, 7, 727–730. [Google Scholar] [CrossRef]
  53. Jende, P.; Nex, F.; Gerke, M.; Vosselman, G. A fully automatic approach to register mobile mapping and airborne imagery to support the correction of platform trajectories in GNSS-denied urban areas. ISPRS J. Photogramm. Remote Sens. 2018, 141, 86–99. [Google Scholar] [CrossRef]
  54. Jende, P.; Nex, F.; Gerke, M.; Vosselman, G. A guided registration strategy employing virtual planes to overcome non-standard geometries—Using the example of mobile mapping and aerial oblique imagery. ISPRS-Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2018, XLII-2, 471–477. [Google Scholar] [CrossRef]
  55. Thapa, A.; Pu, S.; Gerke, M. Semantic feature based registration of terrestrial point clouds. ISPRS Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2009, XXXVIII-3/W8, 230–235. [Google Scholar]
  56. Zieher, T.; Toschi, I.; Remondino, F.; Rutzinger, M.; Kofler, C.; Mejia-Aguilar, A.; Schlögel, R. Sensor- and scene-guided integration of TLS and photogrammetric point clouds for landslide monitoring. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2018, XLII-2, 1243–1250. [Google Scholar] [CrossRef]
  57. Aicardi, I.; Nex, F.; Gerke, M.; Lingua, A.M. An Image-Based Approach for the Co-Registration of Multi-Temporal UAV Image Datasets. Remote Sens. 2016, 8, 779. [Google Scholar] [CrossRef]
  58. Huang, R.; Zheng, S.; Hu, K. Registration of Aerial Optical Images with LiDAR Data Using the Closest Point Principle and Collinearity Equations. Sensors 2018, 18, 1770. [Google Scholar] [CrossRef]
  59. Mancini, F.; Dubbini, M.; Gattelli, M.; Stecchi, F.; Fabbri, S.; Gabbianelli, G. Using Unmanned Aerial Vehicles (UAV) for High-Resolution Reconstruction of Topography: The Structure from Motion Approach on Coastal Environments. Remote Sens. 2013, 5, 6880–6898. [Google Scholar] [CrossRef] [Green Version]
  60. Tong, X.; Liu, X.; Chen, P.; Liu, S.; Luan, K.; Li, L.; Liu, S.; Liu, X.; Xie, H.; Jin, Y.; et al. Integration of UAV-Based Photogrammetry and Terrestrial Laser Scanning for the Three-Dimensional Mapping and Monitoring of Open-Pit Mine Areas. Remote Sens. 2015, 7, 6635–6662. [Google Scholar] [CrossRef] [Green Version]
  61. Seier, G.; Stangl, J.; Schöttl, S.; Sulzer, W.; Sass, O. UAV and TLS for monitoring a creek in an alpine environment, Styria, Austria. Int. J. Remote Sens. 2017, 38, 2903–2920. [Google Scholar] [CrossRef]
  62. Li, W.; Wang, C.; Zai, D.; Huang, P.; Liu, W.; Wen, C.; Li, J. A Volumetric Fusing Method for TLS and SFM Point Clouds. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 3349–3357. [Google Scholar] [CrossRef]
  63. Vasilakos, C.; Chatzistamatis, S.; Roussou, O.; Soulakellis, N. Comparison of Terrestrial Photogrammetry and Terrestrial Laser Scanning for Earthquake Response Management. In Intelligent Systems for Crisis Management; Altan, O., Chandra, M., Sunar, F., Tanzi, T.J., Eds.; Springer International Publishing: Basel, Switzerland, 2019; pp. 33–57. [Google Scholar]
  64. Dold, C.; Brenner, C. Registration of terrestrial laser scanning data using planar patches and image data. ISPRS Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2006, XXXVI-5, 78–83. [Google Scholar] [CrossRef]
  65. Gruen, A.W. Adaptive least squares correlation: A powerful image matching technique. S. Afr. J. Photogramm. Remote Sens. Cartogr. 1985, 14, 175–187. [Google Scholar]
  66. Agisoft PhotoScan. Available online: https://www.agisoft.com (accessed on 21 June 2018).
  67. Fischler, M.A.; Bolles, R.C. Random Sample Consensus: A Paradigm for Model Fitting with Applications to Image Analysis and Automated Cartography. Commun. ACM 1981, 24, 381–395. [Google Scholar] [CrossRef]
  68. Soudarissanane, S.; Lindenbergh, R.; Menenti, M.; Teunissen, P. Scanning geometry: Influencing factor on the quality of terrestrial laser scanning points. ISPRS J. Photogramm. Remote Sens. 2011, 66, 389–399. [Google Scholar] [CrossRef]
  69. Tan, K.; Zhang, W.; Shen, F.; Cheng, X. Investigation of TLS Intensity Data and Distance Measurement Errors from Target Specular Reflections. Remote Sens. 2018, 10, 1077. [Google Scholar] [CrossRef]
  70. Zámečníková, M.; Neuner, H.; Pegritz, S.; Sonnleitner, R. Investigation on the influence of the incidence angle on the reflectorless distance measurement of a terrestrial laser scanner. Iulia Dana Negula Violeta Poenaru 2015, 103, 208–218. [Google Scholar]
  71. Kregar, K. Optimizacija Postopkov Terestričnega Laserskega Skeniranja za Meritve Visoke Natančnosti—Optimization of Terrestrial Laser Scanning for High Precision Measurements. Ph.D. Thesis, University of Ljubljana, Ljubljana, Slovenia, 2016. [Google Scholar]
  72. Mikhail, E.M.; Ackermann, F.E. Observations and Least Squares, 1st ed.; Thomas, Y., Ed.; Crowell Company, Inc.: New York, NY, USA, 1976; ISBN 0-7002-2481-5. [Google Scholar]
  73. Kregar, K.; Grigillo, D.; Kogoj, D. High precision target center determination from a point cloud. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2013, II-5/W2, 139–144. [Google Scholar] [CrossRef] [Green Version]
  74. Urbančič, T.; Vrečko, A.; Kregar, K. The reliability of RANSAC method when estimating the parameters of geometric object. Geodetski Vestnik 2016, 60, 69–83. [Google Scholar] [CrossRef]
  75. Chai, T.; Draxler, R.R. Root mean square error (RMSE) or mean absolute error (MAE)? —Arguments against avoiding RMSE in the literature. Geosci. Model Dev. 2014, 7, 1247–1250. [Google Scholar] [CrossRef]
  76. Cramer, M.; Przybilla, H.J.; Zurhorst, A. UAV cameras: Overview and geometric calibration benchmark. ISPRS-Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2017, XLII-2/W6, 85–92. [Google Scholar] [CrossRef]
  77. Agüera-Vega, F.; Carvajal-Ramírez, F.; Martínez-Carricondo, P. Accuracy of Digital Surface Models and Orthophotos Derived from Unmanned Aerial Vehicle Photogrammetry. J. Surv. Eng. 2017, 143, 04016025. [Google Scholar] [CrossRef]
  78. Lichti, D.D. Error modelling, calibration and analysis of an AM–CW terrestrial laser scanner system. ISPRS J. Photogramm. Remote Sens. 2007, 61, 307–324. [Google Scholar] [CrossRef]
  79. Hebert, M.; Krotkov, E. 3D measurements from imaging laser radars: How good are they? Image Vis. Comput. 1992, 10, 170–178. [Google Scholar] [CrossRef]
Figure 1. New target: (a) Target components; (b) upper part of the target, which is used for terrestrial laser scanning (TLS), and lower part of the target, which is used for unmanned aerial vehicle (UAV) surveys, are assembled with the connecting element.
Figure 1. New target: (a) Target components; (b) upper part of the target, which is used for terrestrial laser scanning (TLS), and lower part of the target, which is used for unmanned aerial vehicle (UAV) surveys, are assembled with the connecting element.
Sensors 19 03179 g001
Figure 2. Diagram of algorithm used to define coordinates of reference point for new target.
Figure 2. Diagram of algorithm used to define coordinates of reference point for new target.
Sensors 19 03179 g002
Figure 3. Input point cloud of the target.
Figure 3. Input point cloud of the target.
Sensors 19 03179 g003
Figure 4. Segmented point cloud of the new target, following primary segmentation.
Figure 4. Segmented point cloud of the new target, following primary segmentation.
Sensors 19 03179 g004
Figure 5. Elongated distances to the points on the black part of the horizontal plane.
Figure 5. Elongated distances to the points on the black part of the horizontal plane.
Sensors 19 03179 g005
Figure 6. Segmented point cloud of the new target, following secondary segmentation. Reference target point is marked by “×”.
Figure 6. Segmented point cloud of the new target, following secondary segmentation. Reference target point is marked by “×”.
Sensors 19 03179 g006
Figure 7. Targets for georeferencing TLS and UAV point clouds: (a) Cylindrical retroreflector for the registration of individual scanner stations’ point clouds; (b) Tilt & Turn target for georeferencing TLS point clouds; (c) photogrammetric ground control point (GCP); (d) retroreflective target that was used as a check point.
Figure 7. Targets for georeferencing TLS and UAV point clouds: (a) Cylindrical retroreflector for the registration of individual scanner stations’ point clouds; (b) Tilt & Turn target for georeferencing TLS point clouds; (c) photogrammetric ground control point (GCP); (d) retroreflective target that was used as a check point.
Sensors 19 03179 g007
Figure 8. Volume targets for registering TLS and UAV point clouds: (a) Small sphere; (b) large sphere; (c) cone; (d) new target.
Figure 8. Volume targets for registering TLS and UAV point clouds: (a) Small sphere; (b) large sphere; (c) cone; (d) new target.
Sensors 19 03179 g008
Figure 9. Distribution of targets, TLS scanner stations, and MS50 station throughout the test area.
Figure 9. Distribution of targets, TLS scanner stations, and MS50 station throughout the test area.
Sensors 19 03179 g009
Figure 10. Group of targets in the test area.
Figure 10. Group of targets in the test area.
Sensors 19 03179 g010
Figure 11. Photogrammetric survey: (a) Photogrammetric block UAV40 (oriented images and dense point cloud); (b) image sample (colored in red in (a)).
Figure 11. Photogrammetric survey: (a) Photogrammetric block UAV40 (oriented images and dense point cloud); (b) image sample (colored in red in (a)).
Sensors 19 03179 g011
Figure 12. Segment taken from a dense point cloud, created from a UAV40 photogrammetric block.
Figure 12. Segment taken from a dense point cloud, created from a UAV40 photogrammetric block.
Sensors 19 03179 g012
Figure 13. The R M S E 2 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point, given in cm. “Georef” represents the RMSE prior to the registration of the TLS and UAV point clouds.
Figure 13. The R M S E 2 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point, given in cm. “Georef” represents the RMSE prior to the registration of the TLS and UAV point clouds.
Sensors 19 03179 g013
Figure 14. The R M S E 3 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point, given in cm. “Georef” represents the RMSE prior to the registration of the TLS and UAV point clouds.
Figure 14. The R M S E 3 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point, given in cm. “Georef” represents the RMSE prior to the registration of the TLS and UAV point clouds.
Sensors 19 03179 g014
Figure 15. The TLS (blue color) and UAV40 (grayscale) point clouds: (a) Before registration; (b) after registration with the new target.
Figure 15. The TLS (blue color) and UAV40 (grayscale) point clouds: (a) Before registration; (b) after registration with the new target.
Sensors 19 03179 g015
Figure 16. The R M S E 2 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50.
Figure 16. The R M S E 2 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50.
Sensors 19 03179 g016
Figure 17. The R M S E 3 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50.
Figure 17. The R M S E 3 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50.
Sensors 19 03179 g017
Figure 18. (a) Cone in the TLS point cloud (left) and noise-distorted cone in the UAV point cloud (right); (b) large sphere in the TLS point cloud (left) and flattened large sphere in the UAV point cloud (right).
Figure 18. (a) Cone in the TLS point cloud (left) and noise-distorted cone in the UAV point cloud (right); (b) large sphere in the TLS point cloud (left) and flattened large sphere in the UAV point cloud (right).
Sensors 19 03179 g018
Table 1. Standard deviation of the registration on TLS scanner stations.
Table 1. Standard deviation of the registration on TLS scanner stations.
Scanner Stationσ3D (mm)
S11.0
S21.3
S31.7
S42.2
Table 2. Characteristics of the photogrammetric blocks. GSD: Ground sample distance.
Table 2. Characteristics of the photogrammetric blocks. GSD: Ground sample distance.
BlockFlight Altitude (m)Number of ImagesGSD (cm)
UAV20201120.49
UAV40Transversal 40, longitudinal 461751.09
UAV70Transversal 75, longitudinal 83502.06
Table 3. Root mean square error (RMSE) on GCPs and check points after the georeferencing of the photogrammetric image blocks. The RMSE is given in cm.
Table 3. Root mean square error (RMSE) on GCPs and check points after the georeferencing of the photogrammetric image blocks. The RMSE is given in cm.
Point Type UAV 20   R M S E 2 D UAV 20   R M S E 3 D UAV 40   R M S E 2 D UAV 40   R M S E 3 D UAV 75   R M S E 2 D UAV 75   R M S E 3 D
GCP0.581.030.872.010.391.38
check 0.541.140.851.930.431.77
Table 4. Two-dimensional precision of coordinates of reference target points (in cm). sSph: Small sphere; LSph: Large sphere; new T: New target; rRef: Retroflective check point target.
Table 4. Two-dimensional precision of coordinates of reference target points (in cm). sSph: Small sphere; LSph: Large sphere; new T: New target; rRef: Retroflective check point target.
Point TypeMS50TLSUAV20UAV40UAV75
sSph0.0040.0040.220.71/
LSph0.0030.0030.070.210.45
Cone0.020.010.150.390.70
newT0.170.290.500.840.39
rRef0.150.380.570.870.46
Table 5. Three-dimensional precision of coordinates of reference target points (in cm).
Table 5. Three-dimensional precision of coordinates of reference target points (in cm).
Point TypeMS50TLSUAV20UAV40UAV75
sSph0.0050.0050.240.76/
LSph0.0040.0030.080.260.55
Cone0.120.140.652.405.23
newT0.170.431.152.031.76
rRef0.150.451.141.821.79
Table 6. Root mean square error (RMSE) calculated from differences between coordinates of the check points in TLS and UAV point clouds, prior to the performed registration.
Table 6. Root mean square error (RMSE) calculated from differences between coordinates of the check points in TLS and UAV point clouds, prior to the performed registration.
RMSE (cm)UAV20UAV40UAV75
2D0.820.980.82
3D1.361.901.94
Table 7. The R M S E 2 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point (in cm).
Table 7. The R M S E 2 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point (in cm).
Tie PointUAV20UAV40UAV75
sSph0.920.64/
LSph0.960.921.66
Cone0.810.750.66
newt0.590.560.67
Table 8. The R M S E 3 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point (in cm).
Table 8. The R M S E 3 D on check points, following the registration of UAV and TLS point clouds, taking into account the used tie point (in cm).
Tie PointUAV20UAV40UAV75
sSph7.362.97/
LSph1.822.434.70
Cone6.031.503.74
newt1.030.681.00
Table 9. The R M S E 2 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50, provided in cm.
Table 9. The R M S E 2 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50, provided in cm.
Tie PointUAV20UAV40UAV75TLS
sSph0.620.99/0.42
LSph0.440.501.700.45
Cone0.890.700.740.39
newT0.270.240.360.41
Table 10. The R M S E 3 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50, provided in cm.
Table 10. The R M S E 3 D of the check points following the UAV and TLS point cloud registration into the referential measurements performed with MS50, provided in cm.
Tie PointUAV20UAV40UAV75TLS
sSph7.523.27/0.50
LSph1.732.524.920.53
Cone6.682.033.140.81
newT0.690.360.680.49

Share and Cite

MDPI and ACS Style

Urbančič, T.; Roškar, Ž.; Kosmatin Fras, M.; Grigillo, D. New Target for Accurate Terrestrial Laser Scanning and Unmanned Aerial Vehicle Point Cloud Registration. Sensors 2019, 19, 3179. https://doi.org/10.3390/s19143179

AMA Style

Urbančič T, Roškar Ž, Kosmatin Fras M, Grigillo D. New Target for Accurate Terrestrial Laser Scanning and Unmanned Aerial Vehicle Point Cloud Registration. Sensors. 2019; 19(14):3179. https://doi.org/10.3390/s19143179

Chicago/Turabian Style

Urbančič, Tilen, Žiga Roškar, Mojca Kosmatin Fras, and Dejan Grigillo. 2019. "New Target for Accurate Terrestrial Laser Scanning and Unmanned Aerial Vehicle Point Cloud Registration" Sensors 19, no. 14: 3179. https://doi.org/10.3390/s19143179

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

Article Metrics

Back to TopTop