Analysis of Optimal Flight Parameters of Unmanned Aerial Vehicles (UAVs) for Detecting Potholes in Pavements

: Pavement maintenance seeks to provide optimal service conditions. Before maintenance, it is necessary to know the condition of the pavement by inspection, a crucial step in deciding on the repair to be carried out. In this sense, unmanned aerial vehicles (UAVs) seem to be an economic substitute compared to the ground laser scanner for pavement inspection tasks. This research seeks to develop a method to measure potholes using 3D models generated with photographs acquired by a UAV and process them using a software based on the Structure from Motion-MultiView Stereo (SfM–MVS) technique. The contribution of this document is the proposal of recommendations for the acquisition of photographs for the realization of the models. To develop these recommendations, an experiment was carried out to evaluate the accuracy in the reconstruction of 3D models using images obtained from the variation and combination of ﬂight planning parameters and data capture. Then, to validate these recommendations, a bumpy section of pavement was modeled using the SfM–MVS method. The results show that for heights of 10 and 15 m the use of this methodology is applicable for the measurement of the width and depth of potholes. Author Contributions: This paper represents the results of teamwork. E.R.-C., S.V.-Q. and F.M.-L.R. designed the research methodology. E.R.-C. and S.V.-Q. carried out the literature review, methods, experiments and results. All of the authors worked on the discussions and conclusions the manuscript. Finally, F.M.-L.R. and E.A. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.


Introduction
The maintenance of roads is a procedure that seeks to deliver optimal service conditions for the routes, generating the least possible impact on the regularity of vehicle flow during operation [1]. Prior to maintenance, it is necessary to know the physical and functional condition of the road, which requires extensive inspection work, focused on knowing its level of damage. It aims to provide measures to assess the magnitude of various types of pavement deterioration, delivering the overall condition of the pavement and providing an estimate of the amount of resources needed to preserve the infrastructure [2]. Inspection is crucial to decide the type and extent of repair a road should receive, and thus provide the desired level of service to maintain acceptable conditions for vehicle traffic [3].
Traditional manual road inspection is time consuming, labor intensive and subjective [2]. Some automated methods, which use equipment with cameras or laser technology, have also been implemented in pavement inspection, greatly improving the efficiency and objectivity of inspection [4]. However, multiple trips may be required to cover the entire width of the road with the inspection vehicle. In addition, the on-the-spot inspection process affects traffic flow, which is particularly unfavorable on high-traffic roads. Unmanned aerial vehicles (UAVs) have the advantage of high 2 of 33 flexibility, relatively low cost compared to inspection vehicles, easy maneuverability and less work on the ground [5]; therefore, they are promising in the task of assessing the condition of the pavement.
On the other hand, recent advances in both high-resolution optical cameras and image processing techniques have allowed three-dimensional (3D) models of UAV images to be derived with sufficient accuracy and high efficiency, making image-based 3D models a possible substitute or more economical supplement to 3D laser scanning. Currently, 3D models can easily be generated from images obtained with UAVs. There are several existing 3D modeling programs on the market which can do this automatically.
This article focuses on the measurement of potholes in asphalt pavements, namely depth, width and volume, using a 3D model derived from UAV images. This document aimed to use photographs obtained with a UAV platform and subsequently process them with the Structure from Motion-MultiView Stereo (SfM-MVS) method, in order to perform a better inspection with the primary goal of detecting potholes in pavements. The main contribution of this document is the proposal of recommendations for measuring the geometric characteristics of potholes based on 3D models.
In order to develop these recommendations, an experiment was conducted to evaluate the accuracy in the reconstruction of 3D models using images obtained from the variation and combination of different parameters of flight planning and data capture (digital photographs). Specifically, variations in the combinations of camera angles (vertical-oblique), image overlaps and flight height, relative to a study surface that simulates deteriorating pavement with bumps of different sizes, with the specific purpose of extrapolating the results to actual pavement. The accuracy was evaluated from the comparison of the actual values against the measurement made in each of the 3D models reconstructed from the data obtained by each flight, according to each of the defined mission plans. Bumpy pavement is modeled using the SfM-MVS technique with photographs acquired using the recommendations and via a UAV platform. A comparison of the geometric characteristics (depth, width and volume) of pavement potholes-measured in a 3D model versus those collected in field inspection-demonstrates that the use of the recommendations when extrapolated to real cases is applicable in the practice of pavement inspection.

Materials and Methods
For the development of this research, the methodology in Design Science Method (DSRM) was used as the base [6][7][8][9]. Tool components and activities were added to show the research done in more detail. In general, the research methodology is organized into five stages: (1) Identification of observed problems and motivations; (2) Goal definition of the potential solution; (3) Design and development; (4) Demonstration; and (5) Evaluation. Figure 1 specifies the activities and research tools that describe the scenario for each stage.
In the first stage, a literature review was carried out to identify emerging technological tools and problems used in pavement inspection and to define the advantages, disadvantages and potential of using a UAV platform in conjunction with the SfM-MVS method. The literature review was carried out through the web databases of Web of Science and Scopus. In the second stage, the research objective was defined. In the third stage, a solution was designed for the application of a project with UAV-SfM-MVS. This project seeks to materialize the previous scenario. This solution is expected to deliver a clear workflow for the acquisition and processing of images obtained through a UAV platform using the SfM-MVS technique. In the fourth stage, a study surface was generated, which represents different sizes of potholes on a pavement. This study surface sought to demonstrate the expected use of the solution created through its 3D modeling by multiple combinations of variables. Finally, the conclusions obtained were translated into practical recommendations for the solution.
In the fifth stage, the solution and recommendations obtained in the previous stage were evaluated through a comparison with a real case study. Finally, at this same stage, the fulfilment of the main research objective is questioned, in addition to future work related to the solution. Appl. Sci. 2020, 10,   Comparative analysis of the accuracy of the 3D modeling with the actual measurements of the study surface. (

5) Evaluation
Evaluation of the relative accuracy of a real case as a result of recommendations made through the analysis of a combination of variables applied to a study area.
Generation of recommendations to improve accuracy with the use of UAVs -SfM-MVS Generation of study surface representing different types of potholes in a pavement (

4) Demonstration
Use of a study surface for the creation of 3D models resulting from the combination of parameters that influence the accuracy of the modeling to later validate the results in a real application. Generation   In the first stage, a literature review was carried out to identify emerging technological tools and problems used in pavement inspection and to define the advantages, disadvantages and potential of using a UAV platform in conjunction with the SfM-MVS method. The literature review was carried out through the web databases of Web of Science and Scopus. In the second stage, the research objective was defined. In the third stage, a solution was designed for the application of a project with UAV-SfM-MVS. This project seeks to materialize the previous scenario. This solution is expected to deliver a clear workflow for the acquisition and processing of images obtained through a UAV platform using the SfM-MVS technique. In the fourth stage, a study surface was generated, which represents different sizes of potholes on a pavement. This study surface sought to demonstrate the expected use of the solution created through its 3D modeling by multiple combinations of variables. Finally, the conclusions obtained were translated into practical recommendations for the solution. In the fifth stage, the solution and recommendations obtained in the previous stage were evaluated

Background
The maintenance of roads, asphalt or concrete, is a procedure that seeks to deliver optimal service conditions to roads, generating the least possible impact on the regularity of the vehicle flow while doing so [1]. Prior to the maintenance work, it is necessary to know the physical and functional condition of the road, requiring extensive inspection work, focused on knowing its level of damage [5]. It provides measures to assess the magnitude of various types of pavement deterioration, delivering the overall condition of the pavement and providing an estimate of the amount of resources needed to preserve the infrastructure [2]. Inspection is the crucial element in deciding the type and extent of the repair that a road should receive, in order to provide the desired level of service and to maintain acceptable conditions for vehicle traffic [3]. The condition of pavement can be studied from two perspectives: structural failures or functional failures. The first refers to the collapse or breakage of one or more layers of pavement, preventing it from withstanding loads on the surface. The second type of failure causes a deterioration in the surface of the pavement, namely unevenness, cracks and formation of potholes, decreasing the service and standard levels of the pavement, causing inconvenience to passengers and negative impact on vehicles [10,11]. A deterioration is an irregularity on the surface of the road which affects the comfort and safety of the vehicle driver. As seen in [12], there are five categories of impairment: (1) cracks by fatigue, block cracks, longitudinal cracks, cross cracks and reflected cracks; (2) deteriorated patches and potholes (it should measure at least 15 cm in diameter); (3) pavement surface deformations, which include ripples and hollowing; (4) surface defects consisting of all deformations that affect the top layer of the pavement; and (5) berm descent, emergence of fines and water, and separation between the berm and pavement.

Traditional Methods for Inspecting Pavement Deterioration
Manual methods were for many years the traditional approach to carrying out pavement inspection [13]. The manual method is to propose a sectioning of the road in different inspection units in its length and width, achieving a detailed review. The revised units are usually at least 15 m long and can be up to four lanes wide, being inspected from a moving or walking vehicle, where the inspector visually identifies defects, and the affected area is estimated as a fraction of the road surface; this data is transferred to paper forms that are taken to the site [14]. This manual approach is costly, laborious, dangerous, subjective and tends to lead to inconsistencies in measurements [15,16].
With the motivation of an automated approach and to minimize these inconsistencies and disadvantages, technologies and procedures have emerged for data collection or processing [2]. The first clear benefit of an automated approach is the reduction of errors associated with transferring data from paper forms used in manual methods. The first way to collect data was through analog images. This consists of photographing the surface of the pavement, usually with 35-mm film. Generally, a downward-facing camera and, if possible, one or more forward-looking cameras or in other configurations, is installed in an inspection vehicle for data capture. Another technology for collecting data is through digital images which can be processed directly by computers. Like analog technology, these cameras are installed in an inspection vehicle, obtaining high-resolution images, higher than normal analog cameras [17]. Methods that use area scanning cameras refer to a photo that represents an area of pavement defined by thousands of pixels. One last technology is that of 3D lasers. This technology uses a laser radar, which uses an exploration laser and reflector to measure reflection times across the pavement surface, thus establishing a three-dimensional pavement surface [18].
On the other hand, data processing technologies can be divided into two groups. In the first group, the process is mainly manual and involves a trained evaluator sitting in a workstation where the images of the pavement are systematically reviewed, and the different deteriorations are identified and classified. In addition, the process is very demanding, as testers must be able to coordinate the simultaneous use of several monitors while tracking observed deteriorations and entering those observations into a rating software [19]. In the second group, fully automated in the context of crack analysis, involves the use of digital recognition software capable of recognizing and quantifying grayscale variations that are related to cracks or cracks on a pavement surface [20]. Various authors [21][22][23][24][25][26] have published methods based on computational algorithms for automatic detection of deterioration; however, these have mainly focused on crack detection, with the exception of [27] which detects two types of deterioration.

New Technologies for Inspecting Pavement Deterioration
Among the technologies available for detection of deterioration in pavement, the laser scanner is one of the most recent [28]. It is based on two principles: flight time and phase change. Flight time sensors estimate the distance between the target and the instrument center by measuring the elapsed Appl. Sci. 2020, 10, 4157 5 of 33 time between the emitted and reflected signal, while phase shift sensors are based on the measurement of angular displacement between the emitted and reflected signal. The output obtained from this technology is a dense point cloud. The point cloud, in conjunction with algorithms, can be used to automatically detect road potholes [29]. In addition, the point cloud can be used for general road asset management due to the large amount of information it contains [30,31]. Lasers have been proven to provide high accuracy but it is an expensive resource [32]. In response to the above, coupled with recent advances in both high-resolution digital cameras and image processing techniques, it is now possible to derive three-dimensional models from images obtained, for example, from a UAV, with sufficient accuracy and high efficiency [33,34]. Therefore, image-based three-dimensional models can be an economical substitute for three-dimensional laser scanners [35].
A UAV is an aerial vehicle that can fly without a human pilot [36]. UAVs have become extremely popular, especially in civil applications, due to their low cost and practicality [36][37][38][39]. A growing area of interest for UAV applications is the periodic inspection and evaluation of infrastructures, which is currently an expensive and time-consuming task [40]. Research has already been carried out that evaluates the UAV as a technology for data capture in pavement inspection applications, with good results due to many advantages it has [10,16]. The operating cost for road network inspection is potentially lower than that of using monitoring vehicles. UAVs can approach scenes that would otherwise be difficult to access or even unattainable under safe circumstances [41]. However, at a certain flight height, the spatial resolution of these images limits the ability to detect certain pavement deteriorations, such as individual cracks, because their widths may be less than 1 cm. In addition, especially due to legal limits and battery life, this technology is not usually an alternative as the length of the road to map increases [35].
For 3D reconstruction from image processing, one of the techniques currently used is Structure from Motion (SfM). For reconstruction, SfM relies on images taken from many points of view [42]. Commonly, processor reconstruction is referred to as SfM. However, it is only a step in the flow of image processing (process for 3D reconstruction). Precisely, this step begins after a key point filtration, where the Scale Invariant Feature Transform (SIFT) algorithm and its variations are the most common approaches [43]. Typically, the SfM process ends with two outputs, a scattered point cloud, and camera parameters and image parameters. A scattered point cloud is an intermediate step in the production of much denser point clouds. These dense point clouds are produced by a process called Multi-View Stereo (MVS) from known intrinsic and extrinsic camera parameters. Comparatively, in a dense point cloud versus a scatter, the points are at least two orders of magnitude denser than the sparse point cloud [44].
Data processing begins with feature detection, by identifying common points in several different photos. Because of the fact that there are geometric distortions, it is necessary to find feature points, which are invariant to changes in scale and orientation in order to make matches across a wider region [45]. A filtration process is then applied to the characteristic or key points to ensure that only correct matches remain [43]. After SfM is executed, what is obtained from this step is a 3D reconstruction of the scene, camera positions and orientations, and the camera's intrinsic parameters. This reconstruction lacks absolute distance between camera poses and between points because scale cannot be recovered from a monocular camera alone [46]. Scaling and georeferencing in the monocular camera case require a minimum of three control points; that is, points raised by topography that serve as the basis for the georeferencing. Another alternative is by using known camera positions obtained from GPS measurements and an inertial measurement unit (IMU) [47]. Finally, the MVS process is responsible for rebuilding a dense point cloud (increasing the number of point clouds obtained from the previous process) [44].

Image Acquisition and Processing
Structure from Motion (SfM) is a technique for generating a three-dimensional reconstruction from 2D images. Unlike traditional photogrammetry, which Bonneval says [48] is "the technique that aims to accurately study and define the shape, dimensions and position in the space of any object using essentially measurements made on one or more photographs," SfM uses computer vision algorithms to identify matching characteristics in a set of overlapping digital images, calculating the location and orientation of the camera. Based on these calculations, overlapping images can be used to reconstruct a structure using a 3D point cloud of the object, surface, or photo scene. It is possible to refine and densify the 3D structure generated by SfM using Multi-View Stereo (MVS), thus completing the SfM-MVS workflow, requiring photographs acquired by means of an optical sensor. For the reconstruction process to be successful, a large number of well-exposed photographs of the object, scene or surface of interest must be obtained at the stage of image acquisition, i.e., with a balance between shutter opening, shutter speed and ISO to correctly capture light. In addition, with sufficient resolution and overlap so that the computer vision algorithm can work effectively. The following are the steps of the method used for image acquisition and processing processes, which are schematically detailed in Figure 2.  The technique of acquiring photographs from a UAV comes from a block of photographs formed from parallel flight lines, flown in a snake pattern at a stable altitude, with a constant overlap and a vertical camera angle (90°) [49]. However, the integration of oblique photographs can reduce the systematic deformation resulting from inaccurate calculation to determine the internal geometry of the camera in the modern SfM-MVS photogrammetry [50,51]. There is evidence that oblique photographs contribute to the integrity of the point cloud reconstruction [52].
Some recommendations for the acquisition of images are as follows: each point of the area of interest should appear in at least three images acquired at different positions; capture a static scene; The technique of acquiring photographs from a UAV comes from a block of photographs formed from parallel flight lines, flown in a snake pattern at a stable altitude, with a constant overlap and a vertical camera angle (90 • ) [49]. However, the integration of oblique photographs can reduce the systematic deformation resulting from inaccurate calculation to determine the internal geometry of the camera in the modern SfM-MVS photogrammetry [50,51]. There is evidence that oblique photographs contribute to the integrity of the point cloud reconstruction [52].
Some recommendations for the acquisition of images are as follows: each point of the area of interest should appear in at least three images acquired at different positions; capture a static scene; use constant light, so that the color of the features does not vary; avoid overexposed or underexposed images; avoid blurry images, slow shutter speed or camera movement; and avoid transparent, reflective or homogeneous surfaces, as these surfaces cause difficulties in further processing [53].
An operational challenge of reconstruction through photographs captured from a UAV is to determine the optimal flight parameters required to achieve a better rebuild quality without excessively increasing flight time or processing. There are key aspects that determine the quality of the reconstruction of the image, mainly the solution of the reconstruction and the location accuracy of the points that rebuild the model. To find an adequate ratio between quality and efficiency, the UAV operator has a set of flight parameters to be adjusted; this includes the height and flight speed, image overlap and camera angle, among other technical sensor parameters such as ISO, shutter speed and shutter opening. The overlaps between images should be high enough to correctly produce 3D models, usually requiring front overlap of more than 75% (adjusted by varying the number of images per second) and the side overlap of more than 60% (critical variable adjusted in UAV flight path planning) [54]. The influence of lateral overlap on reconstruction is not well established, compared or tested despite its influence on flight efficiency [55].
For flight planning, it is important to consider the scale of the photographs, which will limit the detection and consequently the accuracy of the size of the objects in 3D reconstruction. This requires predefining a spatial resolution on the ground or the ground sampling distance called the Ground Sample Distance (GSD) used to define the resolution of the digital image, as it is the representation of the pixel in the terrain [48,56]. The correlation between the flight, sensor and image parameters implies inevitable trade-offs between the different factors. That is, features based on the UAV platform, such as altitude and its resistance (flight time), but also geometry and image quality, and processing time. For example, it is clear that the higher the altitude, the more area is covered on the ground and the more area you can fly on one battery charge. At the same time, higher altitudes lead to fewer images per unit area that need to be processed. Unfortunately, altitude also directly affects achievable GSD and therefore the details that can be detected from images. It follows that lower altitudes lead to more images per unit area, resulting in other compensation related to data processing. The higher the spatial resolution (increase in the accuracy) of the sensor, i.e., the larger the images and the more images are processed per unit area and the higher the processing times [57].
The workflow (2) of Figure 2 is based on the SfM-MVS approach to 3D reconstruction models from a set of images in which intrinsic calibration parameters (or internal parameters, are focal length and lens distortion) and extrinsic (used to describe the transformation between the camera and the world) are unknown.
To apply the SfM-MVS technique, it is key to start with the feature detection process. A feature is a relevant part of an image that can refer to specific structures in the image itself, ranging from simple structures such as points or edges to more complex structures such as objects. These must be detected by robust methods in order to be able to locate the same characteristics in successive images regardless of image rotation, scale or changes in lighting [45]. Once key points (relevant characteristics) have been located in each image, their matches must be determined in the successive images of the captured set [43,58]. SfM aims to simultaneously reconstruct the 3D scene structure, positions and orientations of the camera, along with intrinsic calibration parameters. Initially, camera positions are estimated roughly, as photographs with different orientation are mapped and correlated along with common feature sets, assuming that the camera's internal geometry is constant throughout the process [59]. The initial positions of the camera are refined iteratively as more and more solutions are available, which generates the initial point cloud [60]. The point cloud reference system is relative because at this point only relative locations of the camera and scene geometry are obtained [59].
If you need to get a georeferenced SfM model, there are two options: provide camera positions or use ground control points (GPPs). GPPs are artificial targets or natural objects within the area to be mapped with a known position, which can then be identified and used to orient the final model in the computerized work environment. Alternatively, georeferencing and scaling can be performed from known camera positions derived from differential GPS (dGPS) measurements and an inertial measurement unit [47,61]. Even the latter tools are complementary to using control points (a visually identifiable point, which is known to be exactly positioned) on land, as direct georeferencing is used to provide approximate camera locations and then uses external GCP to obtain a more accurate solution [47,62].
Before you apply MVS to the point cloud, there is an optional additional step that may be required in projects with large image sets. As the number of images increases, the computational load of such a focus increases rapidly [63]. RAM requirements increase with the number of images used in rebuilding and set a practical limit on the number of images that can match simultaneously. The solution to this is the grouping of images, that is, splitting a large project into chunks. In [63], there is detailed a pre-processing step known as "grouping views" for MVS (CMVS), which is a method by which the image set is broken down into overlapping view groups to allow dense MVS reconstructions to run in separate groups. The point clouds obtained from the SfM stage is used as a reference for the MVS algorithm, so that it densifies this point cloud, and thus provides a complete reconstruction of the scene [44]. With the dense point clouds containing a wide metric formation on the scanned surfaces, as well as that relative to their color and reflectivity of the material, it is possible to obtain different end products for analysis, such as orthomosaic and DSM.

Experiments
The purpose of the experiment is to evaluate the accuracy in the reconstruction of 3D models using images obtained from the variation and combination of different parameters of flight planning and data capture (digital photographs); specifically, variations in the combinations of camera angles (vertical-oblique), image overlaps and flight height, relative to a study surface that simulates a deteriorating pavement with bumps of different size, with the specific purpose of extrapolating the results to actual pavement. Accuracy was assessed from the comparison of the actual values versus the measurement made in each of the 3D models reconstructed from the data obtained by each flight, according to each of the defined plans.
The study surface corresponds to a rectangular prism of plaster, which has been treated with surface paint for the purpose of simulating the characteristic color of an asphalt pavement, as seen in the Figure 3. Seven concave surfaces of different sizes were made on the study surface, the characteristics of which are exposed in Table 1, each responding to different sizes and severities, according to [12], to assess the capture sensitivity of the sensor or camera. There are three levels of severity to classify potholes: low severity is a pothole with a depth less than 30 mm, mean severity is between a depth of 30 to 50 mm and high severity must be greater than 50 mm.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 34 reconstructions to run in separate groups. The point clouds obtained from the SfM stage is used as a reference for the MVS algorithm, so that it densifies this point cloud, and thus provides a complete reconstruction of the scene [44]. With the dense point clouds containing a wide metric formation on the scanned surfaces, as well as that relative to their color and reflectivity of the material, it is possible to obtain different end products for analysis, such as orthomosaic and DSM.

Experiments
The purpose of the experiment is to evaluate the accuracy in the reconstruction of 3D models using images obtained from the variation and combination of different parameters of flight planning and data capture (digital photographs); specifically, variations in the combinations of camera angles (vertical-oblique), image overlaps and flight height, relative to a study surface that simulates a deteriorating pavement with bumps of different size, with the specific purpose of extrapolating the results to actual pavement. Accuracy was assessed from the comparison of the actual values versus the measurement made in each of the 3D models reconstructed from the data obtained by each flight, according to each of the defined plans.
The study surface corresponds to a rectangular prism of plaster, which has been treated with surface paint for the purpose of simulating the characteristic color of an asphalt pavement, as seen in the Figure 3. Seven concave surfaces of different sizes were made on the study surface, the characteristics of which are exposed in Table 1, each responding to different sizes and severities, according to [12], to assess the capture sensitivity of the sensor or camera. There are three levels of severity to classify potholes: low severity is a pothole with a depth less than 30 mm, mean severity is between a depth of 30 to 50 mm and high severity must be greater than 50 mm.    Data capture from the study surface was performed in a place free of objects that may cause shadow. In addition, these were made on sunny days, starting the capture of images at 12:52 P.M. until 14:52 P.M. (hours with less shadow), without the presence of clouds that could alter the color tones in the images.
The combination of analyzed variables includes the incorporation of oblique photographs into the set of photographs taken perpendicular to the terrain (vertical), overlap between photographs and flight height. The values of the variables and the justification of that value are detailed below: • Incorporation of oblique photographs into a set of vertical photographs; due to the lack of research on the effect of incorporating vertical photographs with oblique images and their influence on the accuracy of the 3D model in cases of pavement inspection, this research was carried out in order to evaluate the accuracy of the results obtained through this configuration, which combined the cases of 90 Vertical photos with any combination have twice as many images compared to 90 • , as it only has the corresponding set of photos. Considering that the area of the study surface is smaller than the area that will capture the UAV in a single image, the area to be captured for 3D reconstruction will be given by images containing a fraction or the entire study surface for each of the variable combinations set out above.
For the development of the experiment, the Phantom 4 Pro UAV was used as a digital imaging tool, the flight characteristics of which are shown in the following Table 2. Considering that the area of the study surface is smaller than the area that will capture the UAV in a single image, the area to be captured for 3D reconstruction will be given by images containing a fraction or the entire study surface for each of the variable combinations set out above.
For the development of the experiment, the Phantom 4 Pro UAV was used as a digital imaging tool, the flight characteristics of which are shown in the following Table 2.

Results and Discussion
In this section, the results are displayed in four parts: (1) different flight plans to assess the influence of the angle on the variation of the error; (2) different flight plans to assess the influence of overlap and height on the variation of the error; (3) the selection of optimal parameters; and (4) the case study to validate the recommendations and results. The analysis of the results was made by the relative error presented by each pothole according to its actual characteristics.

Angle Variation
In Table 3, the series are presented (named C1-C36 and associated with Figures 5-7) associated with different combinations of height and overlap. Each series includes the five angle combinations.  Figures 5-7 show, respectively, the results of the experiment carried out. The horizontal axis corresponds to the variation or combination of camera angle for photo acquisition, while the vertical axis is the average error for potholes for depth ( Figure 5), width (Figure 6), and volume ( Figure 7). The sub-figures correspond to the different severity levels of each of these parameters: (a) low, (b) medium, and (c) high. On the other hand, the curves correspond to the different combinations in Table 3, starting with C1 (lighter color, lower height of photograph) to C36 (stronger color, the higher height of photograph). Thus, for each set of graphs, when the curves are located in the lower areas of each graph, it will represent that this combination Cn has a lower percentage of error in the measurement with the method used, with respect to its real value. On the other hand, those located in the upper part will have errors close to or equal to 100%; that is, in the latter case, it has not been possible to completely digitally reconstruct the pothole and the parameter of interest has not been measured. Figure 5 displays the depth error values for each of the series, each including the variation and combination of camera angles for data acquisition. In Figure 5a, it is observed that the incorporation of an oblique angle in the acquisition of photographs produces lower error levels compared to models produced only with vertical photographs; more precisely, the reduction of error is 20%. Then, in Figure 5b,c the same behavior governs, i.e., the inclusion of an oblique angle decreases the error; however, for these severity levels (medium and high), the reduction of the error is much greater, on the order of 30%. Figure 6 displays the width error values obtained from measuring the models of each series with different combinations and angle variations. In Figure 6, it is observed that for this case the inclusion of angle has a minimum impact (less than 2%) for cases with the least error, which are the relevant ones, since those manage to reliably measure the geometric characteristics of potholes. For higher height cases the reduction of the error when incorporating oblique angles is up to 60%, with errors around 50-60%. Figure 7 displays the volume error values, obtained from the measurement of the models of each series with combinations and angle variations. According to the figure, the horizontal axis is the variation or combination of camera angle for photo acquisition, while the vertical axis in Figure 7a-c is the average error of measuring pothole volume by severity level, low in Figure 7a, middle in Figure 7b and high in Figure 7c. In Figure 7, for the best case, the reduction from considering only the vertical angle (90 • ) is 39%.

Angle Variation
In Table 3, the series are presented (named C1-C36 and associated with Figures 5-7 associated with different combinations of height and overlap. Each series includes the five angle combinations.            The results indicate that the inclusion of angles other than 90 • [-] have a positive impact on accuracy (error) for measuring the depth, width, and volume of potholes. However, this impact is low compared to the time-to-time effort to acquire this extra number of images, and the processing time is almost double. Likewise, there is no specific oblique angle or clear trend to ensure a reduction of the error. Therefore, working with vertical angles in this method remains the most efficient option for the acquisition of images with application in measurement of potholes on pavements, and in the same way remains consistent as a simple and fast method against more accurate technologies. Therefore, from now on, combinations with 90 • angle will be discussed without any other variation or combination. Figure 8 shows the values and behavior of the depth measurement error when varying height and overlap in vertical photo acquisition flights (90-degree camera angle) for different severity of potholes (low severity in Figure 8a, middle severity in Figure 8b and high severity in Figure 8c).  According to Figure 8a, it is observed that the error trend of all overlaps increases as the measurement height grows, which makes sense as the resolution of the photograph (GSD) decreases. In addition, the error variation decreases as the measurement height increases, indicating that errors converge due to their null interpretation, since the 100% error indicates that the variable to be measured in the model is not identified. On the other hand, there is no optimal overlap that applies to all of the heights. For 2 m, the best result was obtained for an overlap of 75%, with an error of 8.8%. For 5 m, the lowest error was 21.3% for a 90% overlap. For 8 m, the 85% overlap achieved the best results, resulting in an error of 37.1%. At 10 m, the minimum error was 57.7%, which is associated with an overlap of 80%, this error is 1.3% lower compared to an overlap of 75%, even more the error is reduced by more than 15% compared to overlaps of 85 and 90%. While, at 15 m, the 90% overlap has the lowest error (89.4%), an error that is 0.5% less compared to an overlap of 85%; moreover, the error is reduced by more than 5% compared to overlaps of 75 and 80%. However, the error is high (+89%). At 20 m, the least error (93.5%) was obtained with an overlap of 90%, the other overlaps had an error between 93.8-98.8%. At 25 m, there is a 100% error for the overlap of 80%, indicating that no depth could be measured for any of the potholes analyzed, while for overlaps of 75 and 90% the error is 87.5% for both and, for the overlap of 85%, an error of 98.8% was obtained. At 30 m, the overlaps of 80, 85 and 90% obtained a 100% error, which as stated above indicates that no depth was measured. Additionally, with the overlap of 75%, a 98.8% error was obtained, indicating that it was measured because the reconstruction of the 3D model created a slope in a pothole. At 40 m, all overlaps get an error of 100%, that is, no elevation was appreciated in the model; for this reason, the interpretation of the potholes in the model is like that of a circular horizontal dark spot. It can be observed that from 15 m that the error is high (about 90%) for the various overlaps, indicating that the analysis for higher heights is no longer representative for this type of measurement. However, it noted that the simulated potholes analyzed in the experiment would not present the same conditions in practice. In other words, since the experiment considers the value of width as twice the depth, due to its semi-spherical shape, it makes it a more conservative analysis, since the width of a pothole is much greater than its depth.

Variation in Overlap and Height
From Figure 8b, an optimal overlap is clearly observed for the heights of 2  According to Figure 8a, it is observed that the error trend of all overlaps increases as the measurement height grows, which makes sense as the resolution of the photograph (GSD) decreases. In addition, the error variation decreases as the measurement height increases, indicating that errors converge due to their null interpretation, since the 100% error indicates that the variable to be measured in the model is not identified. On the other hand, there is no optimal overlap that applies to all of the heights. For 2 m, the best result was obtained for an overlap of 75%, with an error of 8.8%. For 5 m, the lowest error was 21.3% for a 90% overlap. For 8 m, the 85% overlap achieved the best results, resulting in an error of 37.1%. At 10 m, the minimum error was 57.7%, which is associated with an overlap of 80%, this error is 1.3% lower compared to an overlap of 75%, even more the error is reduced by more than 15% compared to overlaps of 85 and 90%. While, at 15 m, the 90% overlap has the lowest error (89.4%), an error that is 0.5% less compared to an overlap of 85%; moreover, the error is reduced by more than 5% compared to overlaps of 75 and 80%. However, the error is high (+89%). At 20 m, the least error (93.5%) was obtained with an overlap of 90%, the other overlaps had an error between 93.8-98.8%. At 25 m, there is a 100% error for the overlap of 80%, indicating that no depth could be measured for any of the potholes analyzed, while for overlaps of 75 and 90% the error is 87.5% for both and, for the overlap of 85%, an error of 98.8% was obtained. At 30 m, the overlaps of 80, 85 and 90% obtained a 100% error, which as stated above indicates that no depth was measured. Additionally, with the overlap of 75%, a 98.8% error was obtained, indicating that it was measured because the reconstruction of the 3D model created a slope in a pothole. At 40 m, all overlaps get an error of 100%, that is, no elevation was appreciated in the model; for this reason, the interpretation of the potholes in the model is like that of a circular horizontal dark spot. It can be observed that from 15 m that the error is high (about 90%) for the various overlaps, indicating that the analysis for higher heights is no longer representative for this type of measurement. However, it noted that the simulated potholes analyzed in the experiment would not present the same conditions in practice. In other words, since the experiment considers the value of width as twice the depth, due to its semi-spherical shape, it makes it a more conservative analysis, since the width of a pothole is much greater than its depth.
From Figure 8b, an optimal overlap is clearly observed for the heights of 2, 5, 8, 15, 20, 25 and 30 m, while, for the heights of 10 and 40 m, its optimum is not apparent at first sight. At 2 m, the least error was obtained (4%), four with an overlap of 85%; the other overlaps (75, 85 and 90%) have an error difference that exceeds 0.6% compared to the overlap of 75%. At 5 m, the optimal overlap is 85%, because with this overlap and height, the error is the lowest, specifically 2.3%; this error is 2.3% lower compared to an overlap of 90%. At 15 m, the least error was obtained (30.2%) with an overlap of 80%, the other overlaps (75, 85 and 90%) have an error difference that exceeds 22% compared to the overlap of 80%. At 8 m and 80% overlap, a 6% error was obtained, which is the lowest-5.9% lower compared to the 75% overlap. At 10 m, the optimal overlap is 75%, because with this overlap and height the error is the lowest, specifically 2.6%; this error is 0.6% lower compared to an overlap of 80%. Simultaneously, the error is reduced by more than 39% compared to the overlaps of 85 and 90%. At 15 m, the least error was obtained (30.2%) with an overlap of 80%, the other overlaps (75, 85 and 90%) has an error difference that exceeds 22% compared to the overlap of 80%. At 20 and 80% overlap, a 57% error was obtained, which is the lowest and is 4% lower compared to the 75% overlap; even better, the error is reduced by more than 15.8% compared to 85 and 90% overlaps. At 25 m the least error was obtained (70.2%) with an overlap of 75%, the other overlaps (80, 85 and 90%) have an error difference that exceeds 14.7%, compared to the overlap of 75%. At 30 m, the least error was obtained (79.6%) with an overlap of 75%, the other overlaps (80, 85 and 90%) have an error difference greater than 16.1% compared to the overlap of 75%. At 40 m, the error of all overlaps exceeds 95%, so the analysis is no longer representative. As mentioned in the previous analysis, the potholes analyzed would not have the same conditions as in practice (width >> depth), so the analysis remains conservative. Additionally, when looking at the same figure, the error of average severity potholes is reduced compared to low severity bumps, because their dimensions are larger.
From Figure 8c, an optimal overlap is observed for the heights of 8, 15, 20, 25 and 40 m, while, for the heights of 2, 5, 10 and 30 m its peak is not apparent at first sight. At 2 m, the least error was obtained (2.1%) with an overlap of 75%, the other overlaps (80, 85 and 90%) have an error difference greater than 1.5% compared to the overlap of 85%. For 5 m, the least error was obtained (0.7%) with an overlap of 85%, which has 0.7% less error compared to a 90% overlap. At 8 m, the least error (6.2%) was obtained with an overlap of 80%, which has 3.7% less error compared to an overlap of 90%. Simultaneously, the error is reduced by more than 5.9% compared to overlaps of 85 and 90%. At 10 m, the least error (1.4%) is for overlaps of 75 and 80%, the other overlaps (85 and 90%) have an error difference greater than 10.8% compared to overlaps of 75 and 80%. At 15 m, the least error was obtained (12.7%) with an overlap of 80%, the other overlaps (75, 85 and 90%) have an error difference greater than 13.1% compared to the overlap of 80%. At 20 m, the least error was obtained (24.1%) with an overlap of 80%, which has between 9.4% to 12.4% less error compared to an overlap of 75 and 85%. Moreover, the error is reduced more than 36.6% compared to the 90% overlap. At 25 m, the least error was obtained (64.8%) with an overlap of 75%, the other overlaps (80, 85 and 90%) have an error difference greater than 8.4% compared to the overlap of 75%. At 30 m, the least error was obtained (82.5%) with an overlap of 90%, which has 2.2% less error compared to an overlap of 75%, and at the same time, the error is reduced by more than 11.8% compared to the overlap of 85% and 17.5% compared to the overlap of 80%. At 40 m, the least error was obtained (92.2%) with an overlap of 90%, while for other overlaps (75, 80 and 85%) 100% errors were obtained, indicating that these measurements are no longer representative for this analysis. The potholes analyzed for high severity meet the conditions to be characterized as potholes, in the sense of their depth and width dimensions. However, it is important to be reminded that in practice the width is much larger than the depth, so the analysis remains conservative. Additionally, it is observed that the high severity bump error is reduced compared to the average severity potholes, because their dimensions are larger. Figure 9 shows the values and behavior of the width measurement error when varying height and overlap in vertical photo acquisition flights (90-degree camera angle) for different severity of potholes (low severity in Figure 9a, middle severity in Figure 9b, high severity in Figure 9c).  According to Figure 9a and compared to the figure above, it is observed that the width measurement has a lower error than the bump depth measurement. At 2 m, the least error was obtained (3.6%) with an 80% overlap, which has 0.2% less error compared to the 85% overlap. At 5 m, the minimum error value (0.9%) was achieved with an overlap of 80%, which has between 1.4 and 2.3% less error compared to the other overlaps (75, 85 and 90%). For 8 m and an overlap of 80%, an error of 1.5% was found, which is the lowest for this height, comparing this error with respect to overlaps of 75 and 85% there is a difference between 2.3 and 3.5% less error respectively, even more, the error is reduced by 5.6% compared to an overlap of 90%. For 10 m, the least error was obtained (2.7%) with an overlap of 80%, which has 2.6% less error compared to an overlap of 75%, moreover, the error is reduced by more than 5.9% compared to overlaps of 85 and 90%. At 15 m, the least error (6.7%) it was obtained with an overlap of 80%, which has 5.7% less error compared to Z an overlap of 75%. Simultaneously, the error is reduced by more than 13.6% compared to overlaps of 85 and 90%. At 20 m, the variability of the error decreases compared to that obtained at 15 m. Additionally, for this height, the minimal error is obtained (16.2%) with an overlap of 80%, the other overlaps (75, 85 and 90%) have an error difference greater than 13.8% compared to the overlap of 80%. At 25 m, the error variation for the different overlaps decreases more compared to a height of 20 m. On the other hand, the least error (39.0%) was achieved with an overlap of 75%, which has 1.1% less error compared to an overlap of 80% and 4.4% less than the 90% overlap. Moreover, the error is reduced by 8.9% compared to an overlap of 85%. At 30 m, the least error was obtained (52.3%) with an overlap of 75%, for other overlaps (80, 85 and 90%) 100% errors were obtained, because the model did not rebuild potholes. At 40 m, all overlaps got 100% errors, which as mentioned earlier, is due to problems of modeling small bump sizes.
From Figure 9b it is observed that the measurement of widths has less error than the measurement of bump depth. In addition, the error obtained is reduced compared to the analysis for low severity potholes. It is observed that up to a 30-m measuring height, the error does not exceed 20% for the different overlaps and at 40 m, error values greater than 45% are reached. At 2 m, the least error was obtained (1%) with an overlap of 80%, which has 0.6% less error compared to the overlap of 85%, while for an overlap of 75 and 90% the error is reduced by 1.4 to 1.1%. At 5 m, the minimum error value (0.5%) was obtained with an overlap of 75%, which has between 0.1 and 0.6% According to Figure 9a and compared to the figure above, it is observed that the width measurement has a lower error than the bump depth measurement. At 2 m, the least error was obtained (3.6%) with an 80% overlap, which has 0.2% less error compared to the 85% overlap. At 5 m, the minimum error value (0.9%) was achieved with an overlap of 80%, which has between 1.4 and 2.3% less error compared to the other overlaps (75, 85 and 90%). For 8 m and an overlap of 80%, an error of 1.5% was found, which is the lowest for this height, comparing this error with respect to overlaps of 75 and 85% there is a difference between 2.3 and 3.5% less error respectively, even more, the error is reduced by 5.6% compared to an overlap of 90%. For 10 m, the least error was obtained (2.7%) with an overlap of 80%, which has 2.6% less error compared to an overlap of 75%, moreover, the error is reduced by more than 5.9% compared to overlaps of 85 and 90%. At 15 m, the least error (6.7%) it was obtained with an overlap of 80%, which has 5.7% less error compared to Z an overlap of 75%. Simultaneously, the error is reduced by more than 13.6% compared to overlaps of 85 and 90%. At 20 m, the variability of the error decreases compared to that obtained at 15 m. Additionally, for this height, the minimal error is obtained (16.2%) with an overlap of 80%, the other overlaps (75, 85 and 90%) have an error difference greater than 13.8% compared to the overlap of 80%. At 25 m, the error variation for the different overlaps decreases more compared to a height of 20 m. On the other hand, the least error (39.0%) was achieved with an overlap of 75%, which has 1.1% less error compared to an overlap of 80% and 4.4% less than the 90% overlap. Moreover, the error is reduced by 8.9% compared to an overlap of 85%. At 30 m, the least error was obtained (52.3%) with an overlap of 75%, for other overlaps (80, 85 and 90%) 100% errors were obtained, because the model did not rebuild potholes. At 40 m, all overlaps got 100% errors, which as mentioned earlier, is due to problems of modeling small bump sizes.
From Figure 9b, it is observed that the measurement of widths has less error than the measurement of bump depth. In addition, the error obtained is reduced compared to the analysis for low severity potholes. It is observed that up to a 30-m measuring height, the error does not exceed 20% for the different overlaps and at 40 m, error values greater than 45% are reached. At 2 m, the least error was obtained (1%) with an overlap of 80%, which has 0.6% less error compared to the overlap of 85%, while for an overlap of 75 and 90% the error is reduced by 1.4 to 1.1%. At 5 m, the minimum error value (0.5%) was obtained with an overlap of 75%, which has between 0.1 and 0.6% less error compared to the other overlaps (80, 85 and 90%). For 8 m and an overlap of 90%, an error of 1.1% was made. At 10 m, the least error was obtained (0.6%) with an overlap of 75 and 85%, the other overlaps (80 and 90%) have an error difference of 0.5% compared to overlaps of 75 and 85%. For 15 m, the least error (2.8%) was obtained with an overlap of 75%, which has 0.2% less error compared to overlaps of 80 and 90%. Moreover, the error is reduced by more than 5.9% compared to the overlap of 85%. At 20 m, the least error (3.5%) is obtained with an overlap of 80%, which is 1.2 and 1.3% less error compared to 90 and 75% overlaps respectively. The error is further reduced, up to 5.7% compared to the 85% overlap. For 25 m, the least error (6.9%) was obtained with an overlap of 80%, which has 1.9 and 2% less error compared to overlaps of 75 and 90%, respectively, while compared to the overlap of 85% the error is reduced by 7.1%. For 30 m, the least error was obtained (13.6%) with an overlap of 90%, which is between 2.3 and 3.1% less error compared to the other overlaps (75, 80 and 90%). At 40 m, the least error was obtained (47.0%) with an overlap of 85%, which has between 9 and 14% less error compared to the other overlaps (75, 80 and 90%).
From Figure 9c, it is observed that the measurement of widths has less error than the measurement of bump depth. The error obtained is reduced compared to the analysis for medium severity potholes. Additionally, note that for up to a 30-m measuring height, the error does not exceed 10.5% for the different overlaps and at 40 m, the overlap of 75% is the only one that exceeds 20% error, while the other overlaps do not exceed 14% error. At 2 m, the least error was obtained (0%) with an overlap of 80%, while for other overlaps the error was less than 1.6%. At 5 m, the minimum error value (0.5%) was obtained with an overlap of 80%, which has between 0.5 and 1% less error compared to the other overlaps (75, 85 and 90%). For 8 m and an overlap of 90%, an error of 0.3% was made, which is the lowest for this height, comparing this error with respect to overlaps of 75 and 80% there is a difference between 1.9 and 0.9% less error respectively, more so, the error is reduced by 2.7% compared to an overlap of 85%. At 10 m, the least error was obtained (0.3%) with an overlap of 80 and 85%, the other overlaps (70 and 90%) have an error difference of 0.4% compared to overlaps of 75 and 85%. At 15 m, the least error was obtained (0.8%) with an overlap of 80%, which has between 0.4 and 0.5% less error compared to overlaps of 90 and 75%, respectively, moreover, the error is reduced by 0.9% compared to the overlap of 85%. At 20 m, the least error is obtained with an overlap of 80%, which is 1.6% and is 0.3% lower compared to the overlap of 90%, moreover, the error is reduced by 1.5% to 2.4% compared to overlaps of 75 and 85%. At 25 m, the least error was obtained (2.8%) with an overlap of 80%, which is between 0.7 and 2.6% less error compared to the other overlaps (75, 85 and 90%). At 30 m, the least error was obtained (5.7%) with an overlap of 90%, which has 1.4% less error compared to an overlap of 80%. Moreover, the error is reduced by 4 to 4.3% compared to overlaps of 75 to 85%. At 40 m, the least error was obtained (11.4%) with an overlap of 90%, which has 1.8% less error compared to overlaps of 80 and 85%. Moreover, the error is reduced by 11.5% compared to an overlap of 75%. Figure 10 shows the values and behavior of the volume measurement error by varying the height and overlap in vertical photo acquisition flights (90-degree camera angle) for different severity in potholes (low severity in Figure 10a, middle severity in Figure 10b, high severity in Figure 10c).  In Figure 10a, a trend is observed of overlapping error increases as the measurement height grows. On the other hand, it is observed that the error at 10 m is greater than 60% for all overlaps, indicating that the experiment does not work well for this type of analysis. At 2 m, the minimum error was 4.1%, which is associated with an overlap of 80%, this error is 2% lower compared to an overlap of 85%, even more the error is reduced by more than 3.6% compared to overlaps of 75 and 90%. At 5 m, the 80% overlap has the lowest error (6.4%), an error that is 9.2% less compared to an overlap of 85%. Moreover, the error is reduced by more than 9.4% compared to 75 and 90% overlaps. At 8 m, the least error (36.6%) was obtained with a 90% overlap. For 10 m and with an overlap of 75% an error of 63.1% was obtained, which is the lowest for this height, since it has 5.5% less error compared to the overlap of 80%. Moreover, the error is reduced by 12.8 to 18.7% compared to the overlaps of 85 and 90%. At 15 m, the least error was obtained (81.9%) with an overlap of 90%, which is between 2.9 and 10.1% less error compared to the other overlaps (75, 80 and 85%). At 20 m, the error is 95% for an 80% overlap, which is high for further analysis of this type of measurement. So, from 20 m, an analysis of volume for low severity potholes loses value in practice, since the volume of the pothole fails to describe a semi-sphere, but rather tries to reconstruct a kind of cone, from the position in which the maximum depth was obtained. The behavior is similar to that observed in the pothole depth analyses, as the measured volume depends directly on how good the measurement has been with respect to depth and area characteristics.
From Figure 10b, a trend is observed of overlapping error increases as the capture height grows. On the other hand, it is observed that the error at 10 m is greater than 30% for overlaps of 75 and 80% and more than 50% for overlaps of 85 and 90%, indicating that the experiment does not work properly for this type of analysis. However, the error is reduced from what was observed in the previous analysis. At 2 m, the minimum error was 1.9%, which is associated with an overlap of 85%, this error is 1.8% lower compared to an overlap of 90%, even more so the error is reduced by more than 2.3% compared to overlaps of 75 and 80%. At 5 m, the 80% overlap has the lowest error (7.3%), an error that is 7.5% less compared to a 90% overlap. Moreover, the error is reduced by more than 7.8% compared to 75 and 85% overlaps. At 8 m, the least error (22.5%) was obtained with a 90% overlap. At 10 m, the least error was obtained (35.8%) with an overlap of 80%, which has 1% less error compared to the overlap of 75%, while for an overlap of 85 and 90% the error is reduced by 18  In Figure 10a, a trend is observed of overlapping error increases as the measurement height grows. On the other hand, it is observed that the error at 10 m is greater than 60% for all overlaps, indicating that the experiment does not work well for this type of analysis. At 2 m, the minimum error was 4.1%, which is associated with an overlap of 80%, this error is 2% lower compared to an overlap of 85%, even more the error is reduced by more than 3.6% compared to overlaps of 75 and 90%. At 5 m, the 80% overlap has the lowest error (6.4%), an error that is 9.2% less compared to an overlap of 85%. Moreover, the error is reduced by more than 9.4% compared to 75 and 90% overlaps. At 8 m, the least error (36.6%) was obtained with a 90% overlap. For 10 m and with an overlap of 75% an error of 63.1% was obtained, which is the lowest for this height, since it has 5.5% less error compared to the overlap of 80%. Moreover, the error is reduced by 12.8 to 18.7% compared to the overlaps of 85 and 90%. At 15 m, the least error was obtained (81.9%) with an overlap of 90%, which is between 2.9 and 10.1% less error compared to the other overlaps (75, 80 and 85%). At 20 m, the error is 95% for an 80% overlap, which is high for further analysis of this type of measurement. So, from 20 m, an analysis of volume for low severity potholes loses value in practice, since the volume of the pothole fails to describe a semi-sphere, but rather tries to reconstruct a kind of cone, from the position in which the maximum depth was obtained. The behavior is similar to that observed in the pothole depth analyses, as the measured volume depends directly on how good the measurement has been with respect to depth and area characteristics.
From Figure 10b, a trend is observed of overlapping error increases as the capture height grows. On the other hand, it is observed that the error at 10 m is greater than 30% for overlaps of 75 and 80% and more than 50% for overlaps of 85 and 90%, indicating that the experiment does not work properly for this type of analysis. However, the error is reduced from what was observed in the previous analysis. At 2 m, the minimum error was 1.9%, which is associated with an overlap of 85%, this error is 1.8% lower compared to an overlap of 90%, even more so the error is reduced by more than 2.3% compared to overlaps of 75 and 80%. At 5 m, the 80% overlap has the lowest error (7.3%), an error that is 7.5% less compared to a 90% overlap. Moreover, the error is reduced by more than 7.8% compared to 75 and 85% overlaps. At 8 m, the least error (22.5%) was obtained with a 90% overlap. At 10 m, the least error was obtained (35.8%) with an overlap of 80%, which has 1% less error compared to the overlap of 75%, while for an overlap of 85 and 90% the error is reduced by 18.6 to 20.1%. At 15 m, the minimum error value (39.1%) was obtained with an overlap of 80%, which has between 20.0 and 25.6% less error compared to the other overlaps (75, 85 and 90%). For 20 m and an overlap of 75%, an error of 64.6% was produced, which is the lowest for this height, comparing this error with respect to overlaps of 80 and 85% there is a difference between 0.8 and 2.8% less error respectively, more than that, the error is reduced by 11.6% compared to an overlap of 90%. At 25 m, the least error was obtained (71.3%) for an 85% overlap, which is between 21.6 and 26.4% less error compared to other overlaps (75, 80 and 90%). At 30 m, the error obtained was 96.5% for an overlap of 75%, which is the lowest and at the same time high considering its magnitude, so it loses practical value. At 40 m, the error is higher than that obtained at a height of 30 m. From the above it is observed that from 30 m it loses practical value.
From Figure 10c, a trend is observed of overlapping error increasing as the capture height grows. However, at 15 m with an overlap of 80%, the volume obtained was 1.8% lower than that obtained at 10 m with the same overlap. This is probably explained by problems in rebuilding with the software. On the other hand, it is observed that the error is similar to that made in the previous analysis, since at 10 m the error is greater than 30% for overlaps of 75 and 80% and more than 40% for overlaps of 85 and 90%, indicating that the experiment does not work properly for this type of analysis. At 2 m, the minimum error was 0.9%, which is associated with a 90% overlap, this error is 1.2% lower compared to an overlap of 85%, even more so the error is reduced by more than 1.9% compared to overlaps of 75 and 80%. At 5 m, the 80% overlap has the lowest error (5.5%), an error that is 7.5% less compared to a 90% overlap. Moreover, the error is reduced by more than 8.8% compared to 75 and 85% overlaps. At 8 m, the least error (24.6%) was obtained with an overlap of 85%. At 10 m, the least error was obtained (35.7%) with an overlap of 80%, which has 0.9% less error compared to the overlap of 75%. Moreover, the error is reduced by 5.6 to 7.8% compared to overlaps of 85 and 90%. At 15 m, and with an overlap of 80% an error of 33.9% was obtained, which is the lowest for that height, and which is 11.9% lower compared to the overlap of 75%. Moreover, the error is reduced by 25.2 to 27.7% compared to overlaps of 85 and 90%. For 20 m, the lowest error is 41.2%, which is associated with an 80% overlap, and has 10% less error compared to an overlap of 75%. Moreover, the error is reduced by 23.8 to 34.6% compared to overlaps of 85 and 90%. At 25 m, the least error was obtained (75.9%) with an 85% overlap, which is between 9.7 and 18.8% less error compared to the other overlaps (75, 80 and 90%). At 30 m, the minimum error (90.2%) was held for an overlap of 90%, which has no practical value. At 40 m, the error is high with respect to the value obtained at a height of 30 m. For the overlaps of 75, 80 and 85%, the volume could not be measured because no depth was recorded.

Selecting Optimal Parameters
From the selection by geometric characteristic type, i.e., width, depth, and bump volume, the overall combination that allows for the most reliable results must be selected, since the characteristics of potholes in a real application will be determined by a single model. In this way, optimal recommendations in depth, width and volume measurements were analyzed, always seeking to minimize the error of these types of measurement. Table 4 was based on the discussion presented in Section 3.2, the recommended overlap was chosen by severity level, as the use of these recommendations varies depending on the minimum level of severity to be measured. For heights of 2, 5 and 8 m, the experiments were excessively laborious. For these heights, it is difficult to ensure proper positioning and correct overlap between images, as the handling of the UAV at such a low height becomes imprecise due to the sensitivity of the control. In addition, the GPS does not deliver a correct position of the aircraft, which in some cases led to the necessary repetition of the shooting. Although capturing photographs at low altitude allows for more accurate 3D models, in practical terms and in order to optimize the flight and inspection (balance in time, quality, area observed, etc.), it is not recommended to perform the flights at the indicated heights. Therefore, it is not recommended to plan flights with heights less than 10 m. At 40 m, the error levels Table 5 shows longitudinal (B) and transverse distances (A) between photographs captured with the UAV to comply with overlaps; longitudinal (p) and transverse (q) for each of the heights.
These distances are used to calculate the flight lines needed to fully cover the pavement unit to be inspected. To know the number of flight lines for a given road width, Equation (1) should be used. An additional line is added to ensure full width coverage. Additionally, to have a reference of the number of photographs to be taken per flight line, it is recommended to use Equation (2). Four additional photos are added, two for the start and two for the end of the flight line, in order to have the ends properly overlapped. Naturally, the corners will have a lower number of photo coverage compared to the center, as can be seen in an example in Figure 11. That is necessary to guarantee a minimum of photographs in the corners of the pavement to have all the areas covered on the map.

Case Study
A study is carried out a real case study to validate whether recommendations can be extrapolated in the measurement of bumps of deteriorated asphalt pavements. The actual pavement section has 6 potholes with different dimensions (indicated in the Figure 12. Flexible pavement pothole for recommendation validation. Figure 12 shows the potholes by indicating by a horizontal arrow the largest width and a vertical arrow the smallest width.

Case Study
A study is carried out a real case study to validate whether recommendations can be extrapolated in the measurement of bumps of deteriorated asphalt pavements. The actual pavement section has 6 potholes with different dimensions (indicated in the Table 6. Flexible pavement pothole for recommendation validation.  Figure 12 shows the potholes by indicating by a horizontal arrow the largest width and a vertical arrow the smallest width.
The first analysis was carried out at a height of 10 [m] applying the recommendations detailed above. The results are shown in Table 7, in addition to Figure 13, where the reconstruction of a pothole in the case study can be observed.
A study is carried out a real case study to validate whether recommendations can be extrapolated in the measurement of bumps of deteriorated asphalt pavements. The actual pavement section has 6 potholes with different dimensions (indicated in the Figure 12. Flexible pavement pothole for recommendation validation. Figure 12 shows the potholes by indicating by a horizontal arrow the largest width and a vertical arrow the smallest width.  The first analysis was carried out at a height of 10 [m] applying the recommendations detailed above. The results are shown in Table 7, in addition to Figure 13, where the reconstruction of a pothole in the case study can be observed. Figure 13. Volumetric representation of a pothole in 3D reconstruction, based on Context Capture (Bentley) software. Table 7. Dimensions of potholes measured in model at 10 m and error.  (e) (f) Figure 12. Flexible pavement pothole for recommendation validation. The first analysis was carried out at a height of 10 [m] applying the recommendations detailed above. The results are shown in Table 7, in addition to Figure 13, where the reconstruction of a pothole in the case study can be observed.   The second analysis went to one of 15 [m], the results are shown in Table 8. It is important to note that potholes were marked for case study planning in order to clearly identify the measured points in the field, in order to minimize the measurement error of the evaluator.

Width Error Depth Error Volume Error
For the analysis at 10 m, the measurements of the height and width were obtained with high accuracy, compared to the other measurement variables, since with the help of the marks, these characteristics could be better determined. However, for depth measurement, the error was greater, with the worst case being 9.9% higher than the value of the recommendation; this can be explained by the irregular shape of the bumpy surface, i.e., the surface is not completely horizontal. However, the measurement error could be obtained with high accuracy (error less than 1 cm). On the other hand, the volume analysis presented a similar error with respect to the value of the recommendation. However, it remains a high error to consider it as a potential variable to be measured with this method.
For the 15 m height analysis, the width measurement was greater than the values obtained at 10 m, although for the pothole (c) a difference of 0.4% was presented. This is probably explained by the variability of measurement in terrain and the sensitivity of the study variable, which for a millimeter error in the measurement results in the error percentage presenting an atypical or unexpected behavior. In the measurement of depth and volume, the error behaved better than expected (according to the results of the study surface) which can be explained by the fact that these potholes have a large surface with respect to their depth; this certainly helps to rebuild the volume and obviously the depth.
As the above results demonstrate, the recommendations allow for a good measure of the geometric characteristics of pavement potholes. The level of error is about one centimeter, which is a promising methodology applied to the regular practice of engineering.

Conclusions
The literature review process allowed us to gather a history of the use of UAVs as a tool for data acquisition. This tool is positioned as an economical substitute compared to more sophisticated emerging technologies, such as the terrestrial laser scanner, used for pavement inspection. It also served as the basis for developing a workflow for photo acquisition and further processing using the SfM-MVS technique. It also understands the importance of variables that influence flight planning.
This research focused on developing a practical method for measuring geometric characteristics of potholes, using 3D models generated with photographs acquired with a UAV and processed through software based on the SfM-MVS technique. To develop this methodology, an experiment was conducted to evaluate the accuracy in the reconstruction of 3D models using images obtained from the variation and combination of different parameters of flight planning and data capture (digital photographs). Specifically, variations in the combinations of camera angles (vertical-oblique), image overlaps and flight height.
The experiment concludes that the incorporation of oblique photographs into a set of vertical photographs does not guarantee a significant reduction in error. Therefore, the research focused on error assessment for different combinations of height and overlap with vertical angle of photo acquisition. Within the geometric characteristics of the pavement studied, the width is the characteristic that has the lowest level of error, followed by depth and last by volume.
After the evaluation of the error, both for depth, width and volume, practical recommendations were generated according to the level of severity of potholes. It was obtained that the methodology is applicable for heights of 10 to 15 m, since at higher heights the error levels do not allow to represent well the characteristics sought, while at lower heights the process becomes extremely laborious when having to perform manual flights; additionally, the GPS becomes sensitive and imprecise at these heights. These recommendations were validated in a real case study, demonstrating that they are extrapolatable to the regular practice of engineering, specifically, in the inspection of asphalt pavements.
The future work of this research should focus on improving the accuracy of the models generated by incorporating ground control points (GCPs), mainly to apply it on road pavements that have milestones or benchmarks that are accessible prior to the acquisition of photographs. The GCPs improve accuracy because they are the first mooring points that improve relative accuracy levels while increasing the absolute position of the model. Moreover, another line of research could include the automatic extraction of characteristics of pavement models generated by SfM-MVS in conjunction with photographs acquired based on the recommendations of this research.