Analysis and Precision of Light Detection and Ranging Sensors Integrated in Mobile Phones as a Framework for Registration of Ground Control Points for Unmanned Aerial Vehicles in the Scanning Technique for Building Information Modelling in Archaeological Sites

: The protection of heritage sites is one of the keys that our civilisation presents. That is why great efforts have been invested in order to protect and preserve movable and immovable property with a certain historical value, as is the case of archaeological sites scattered throughout the territory of southern Iberia (Spain) in the form of dolmens and negative structures dug into the ground, constituting a good sample of the megalithic culture in southern Spain. To study, manage and preserve these archaeological monuments, considered a set of cultural assets, various techniques and methodologies are required to facilitate the acquisition of three-dimensional geometric information. The Scan-to-BIM approach has become one of the most up-to-date work exponents to carry out these objectives. The appearance of LiDAR techniques, and recently their incorporation into smartphones through integrated sensors, is revolutionising the world of 3D scanning. However, the precision of these techniques is an issue that has yet to be addressed in the scientiﬁc community. That is why this research proposes a framework, through experimental measurement, comparison and knowledge of the limitations of this technology, to know the precision of the use of these smartphones, speciﬁcally the iPhone 13 Pro, as a measurement element to establish points of control with the aid of photogrammetry by unmanned aerial vehicles (UAVs) in archaeological sites. The results demonstrate a residual uncertainty of ± 5 mm in the capture of GCPs from the mobile phone’s LiDAR light detection and ranging sensor, and there was a deviation of the measurements in a range between 0 and 28 m of distance between the GCPs of (0.021, 0.069) m.


Introduction
The Industrial Revolution 4.0 is changing the way the industry and construction sectors are perceived, making the use of operational and information technology increasingly important. The paradigm of the fourth industrial revolution involves data management and the interconnection between machines, objects and processes [1]. In this context, digitisation covers both architectural and archaeological fields and is an integral process in the convergence towards industrial digitisation. In particular, the geometric characterisation of the complex forms present in archaeology and related fields implies a detailed and exhaustive description of the geometry of the objects that make up the configuration of the structures.
The idea of conservation of Cultural Heritage has to do with qualification systems where relationships of internal order (principles governing geometric shape) and relationships of shapes and volumes are described. Therefore, geometric characterisation implies the use of 3D digitisation techniques and instruments and interactive environments that bring us closer to a Historic Building Information Model (HBIM) [2] based on geomatic techniques that generate precise and detailed data on the shape, surface, texture and volume, among other geometric properties of the object.
In this sense, to acquire records in the complex framework of complex shapes that occur in archaeology, it is necessary to have certain technologies that have a high acquisition price, such as terrestrial laser scanning (TLS) technology.
However, portable devices such as mobile phones or tablets have emerged since 12 October 2020 that incorporate a sensor with a LiDAR (light detection and ranging) system in different versions of the iPhone, making them accessible to many people and professional fields. This Apple release in 2020 integrated a sensor designed for Augmented Reality (AR) without specifications on relative and absolute accuracies. Despite this, the versatility and applicability of the system challenged the geomatics and telemetry specialist community [3].
On the other hand, Structure from Motion (SfM) short-range photogrammetry has a highly developed space in the fields of architecture and archaeology [4][5][6] due to the enormous amount of data it provides. When this process is developed by unmanned aerial vehicles (UAVs), it acquires greater applications, especially in the registration of building roofs, where the surfaces of the zenithal planes are important since terrestrial laser scanning is not capable of covering them.
Fields such as the dynamic study of the earth's surface have been explored in various works, an example the one developed by Duró et al. [7], who tried to study the alterations in erosion on lake shores by means of UAVs, placing the ground control points (GCPs) in the alluvial plain. This type of study has been carried out in other work related to civil engineering [8] or safety at construction sites [9].
UAV-based photogrammetry typically uses ground control points to aid in measurements and establish precise scales in the model. If the points are made from a total station, georeferencing is indirect, referring to methods that assign world frame coordinates to 3D measurements collected in a relative local reference frame.
In the case of ground control points (GCPs) that have been taken with a global navigation satellite system (GNSS) receiver, the GCPs will be georeferenced in a universal navigation system. In the field of archaeology, the precision of unmanned aerial flights is an important factor, which involves flight plan indicators and georeferencing [10].
The advantages of unmanned aerial vehicles with respect to the terrestrial laser scanner are well known, due to the low cost of some flight equipment and the speed of data acquisition, although this aspect is relative if precision is currently considered for personal laser scanners [11] and the speed of their records and the taking of measurements for the GCPs.
Regarding the study of the precision of integrated mobile phones, one of the first studies was carried out in the field of geoscience [12] on small objects making comparisons with photogrammetry techniques. The choice of fourteen samples presented a singular morphology with sharp edges, and the field samples for its registration took measurements that did not exceed fifty centimetres in length.
In the case of three-dimensional models for small objects studied by some authors [13], it is stated that the limitation of the method is found at a length of 5 m. However, the case studies are varied, from the facade of a building to the sculpture and round bust of a piece of art.
One of the most interesting works on light detection and ranging (LiDAR) precision is the work of Costantino et al. [14], in which they perform a sampling of ten objects with an analysis of the quality of the cloud of points obtained in relation to the variation in the surface, loss of planarity and omnivariance. Another interesting work is the one carried out by Chase et al. [15], who studied the accuracy of the sensor integrated in the iPhone 12/13 Pro in an interior space with a length of 25 m and with different control points that took as reference the TLS and the total station (TS).
In addition, in the framework of architectural topography, Spreafico et al. [3] evaluated the LiDAR sensor for fast 3D mapping.
In summary, there are some very current studies focused on the precision of the aforementioned sensor, which range from the analysis of small objects to interior spaces of greater dimension, reaching its maximum expression at a length of 25 m. A basic question that the researchers expose in the work of Kršák et al. [16] is the low precision of 3D models from UAV photogrammetry. It is necessary to validate the sets of points obtained through algorithms to know in advance the level of precision that one is going to achieve in the records and if these are specially intended to achieve planimetry or to obtain results in the Scan-to-BIM process.
On the other hand, the possibility of using portable devices such as mobile phones or tablets in Pro versions that are relatively cheap and easy to transport presents a new challenge in the field of geomatics, replacing large-scale TS or TLS equipment. Therefore, an operator with a simple drone and a mobile phone could perform a suitable data capture record for Scan-to-BIM in archaeological sites.
For all these reasons, this research proposes a framework where, through measurement and comparison techniques and knowledge of the limitations of the experimental methods, the precision of using the iPhone 13 Pro as a measurement element to establish control points is determined. The aim is to create a framework that assists in photogrammetry through unmanned aerial vehicles (UAVs) in works related to archaeology and heritage architecture. For this purpose, a set of archaeological sites was investigated through different work procedures that will be seen below.

Case Study and 3D Processing
The megalithic ensembles of Los Gabrieles (Valverde del Camino), Spain, are one of the most relevant spaces in terms of quality and number of megalithic monuments. They are located in the province of Huelva (Spain) (Figure 1a). Huelva's megalithism has developed as a cultural phenomenon exclusive to the death of the Copper Age [17]. It was after the discoveries by Obermaier in 1924 of the Dolmen del Soto that this type of construction became popular, especially due to the appearance of grave goods and ceramic artefacts.
As a result of the study of funerary architecture, researchers and archaeologists have agreed that two types of culture were established in Huelva: On the one hand, there is the group of shepherds of the Neolithic tradition who inhabited the Andévalo and the Sierra regions and, on the other hand, there are the tombs with false domes that relate to the culture of the Millares [17].
The territorial context in which these deposits are found is very similar to the rest of those found in the Sierra de Huelva and Seville regions, where mountain orography and forest spaces predominate and where the georeferenced data in a Geographic Information System have been processed using the QGIS software [18], which allowed us to define priorities in the cartographic information (Figure 1b).
Massive Data Capture Systems (MDCS) like aerial photogrammetry processes are ideal for capturing geometry from nadiral planes, which ground-based laser scans cannot reach [19]. Therefore, from a technical point of view, it improves the capture of the process that we call Scan-to-BIM, that is, the obtaining of construction information models of historical buildings or, in this case, archaeological sites. Massive Data Capture Systems (MDCS) like aerial photogrammetry processes are ideal for capturing geometry from nadiral planes, which ground-based laser scans cannot reach [19]. Therefore, from a technical point of view, it improves the capture of the process that we call Scan-to-BIM, that is, the obtaining of construction information models of historical buildings or, in this case, archaeological sites.
For the specific aerial photogrammetry work at these aforementioned sites, we used two small dolmens with an open morphological structure and small dimensions located in the Los Gabrieles complex in the municipality of Huelva. Thus, we could observe the process of photographic shots through two intelligent flights; the image in Figure 2 is a shot at a height of 12 m, and in Figure 3, a second intelligent flight mode based on zenithal and oblique exposures that determine a minimum height of 24 m can be seen. The 3D photogrammetric models that are extracted from the UAVs make it possible to collect a large amount of data in a short time, but like any system, it has its limitations, despite covering an overhead space that does not capture the terrestrial geometry. To obtain an estimate of the accuracies achieved, one of the most common methods from a scientific point of view is to make a comparison between the SfM model and the TLS. For the specific aerial photogrammetry work at these aforementioned sites, we used two small dolmens with an open morphological structure and small dimensions located in the Los Gabrieles complex in the municipality of Huelva. Thus, we could observe the process of photographic shots through two intelligent flights; the image in Figure 2 is a shot at a height of 12 m, and in Figure 3, a second intelligent flight mode based on zenithal and oblique exposures that determine a minimum height of 24 m can be seen. The 3D photogrammetric models that are extracted from the UAVs make it possible to collect a large amount of data in a short time, but like any system, it has its limitations, despite covering an overhead space that does not capture the terrestrial geometry. To obtain an estimate of the accuracies achieved, one of the most common methods from a scientific point of view is to make a comparison between the SfM model and the TLS.   The aerial photogrammetry point cloud was processed using v. 1.8.5 of the Agisoft Metashape [20] software. In the case of the study of the dolmens, seven control points were taken using wooden targets of 37.5 × 25.3 cm strategically placed (i.e., uniformly distributed) on the ground. For the measurements of the control points, measurements were made with a LEICA Flexline TS02 total station. To cover the entire surface, the GCPs were distributed using diagonal geometry and outside the surface of interest typical of the dolmens. The coordinate system was established through axes parallel and perpendicular to some established axes. The local XYZ coordinates set were 100, 100 and 10 m, respectively. Subsequently, the coordinates of the elements in each space were recorded to achieve a uniform set of points.
Accordingly, after importing the images to a Chunk, Agisoft Metashape performed an automatic calibration process using an exchangeable image file (EXIF) [21], detecting each of the images. The processing method used is the one described by Westoby et al. [22], which is carried out from the alignment of the photographs to the postprocessing with mesh generation. The processing parameters of said software were set to "high accuracy" for dense cloud extraction. Once the dense cloud extraction results for the two processes were obtained-one dense cloud with images at a height of 12 m, obtaining a set of 7,046,235 points, and another processing at a height of 24 m with oblique images, obtaining a dense point cloud of 2,281,026 points-the processing was also repeated for Dolmen number 3, which is very close to number 4 (they are separated by approximately ninety metres). The aerial photogrammetry point cloud was processed using v. 1.8.5 of the Agisoft Metashape [20] software. In the case of the study of the dolmens, seven control points were taken using wooden targets of 37.5 × 25.3 cm strategically placed (i.e., uniformly distributed) on the ground. For the measurements of the control points, measurements were made with a LEICA Flexline TS02 total station. To cover the entire surface, the GCPs were distributed using diagonal geometry and outside the surface of interest typical of the dolmens. The coordinate system was established through axes parallel and perpendicular to some established axes. The local XYZ coordinates set were 100, 100 and 10 m, respectively. Subsequently, the coordinates of the elements in each space were recorded to achieve a uniform set of points.
Accordingly, after importing the images to a Chunk, Agisoft Metashape performed an automatic calibration process using an exchangeable image file (EXIF) [21], detecting each of the images. The processing method used is the one described by Westoby et al. [22], which is carried out from the alignment of the photographs to the postprocessing with mesh generation. The processing parameters of said software were set to "high accuracy" for dense cloud extraction. Once the dense cloud extraction results for the two processes were obtained-one dense cloud with images at a height of 12 m, obtaining a set of 7,046,235 points, and another processing at a height of 24 m with oblique images, obtaining a dense point cloud of 2,281,026 points-the processing was also repeated for Dolmen number 3, which is very close to number 4 (they are separated by approximately ninety metres).

Experimental Process Equipment
One of the unmanned aerial vehicles with enormous potential for field work in inspection and photogrammetric records is the DJI Mavic 2 Enterprise Dual drone, whose takeoff weight does not reach 910 g and has a GPS+GLONASS GNSS system with an operating frequency of 2.400-2.4835/5.725-5.850 GHz (Figure 4). It features two cameras contained in the same gimbal: (i) an infrared M2ED camera with an uncooled Vox Microbolometer and an HFOV lens of 57 • with f/1.1 aperture and a 12-micron pixel size; (ii) an M2ED visual camera with a 1/2. 3 CMOS sensor with 12 M effective pixels and an FVO lens of approximately 85 • that is equivalent to a 35 mm format with f/2.8 aperture. Both sensors have output in .jpeg format and MP4 video format. takeoff weight does not reach 910 g and has a GPS+GLONASS GNSS system with an operating frequency of 2.400-2.4835/5.725-5.850 GHz (Figure 4). It features two cameras contained in the same gimbal: (i) an infrared M2ED camera with an uncooled Vox Microbolometer and an HFOV lens of 57° with f/1.1 aperture and a 12-micron pixel size; (ii) an M2ED visual camera with a 1/2. 3″ CMOS sensor with 12 M effective pixels and an FVO lens of approximately 85° that is equivalent to a 35 mm format with f/2.8 aperture. Both sensors have output in .jpeg format and MP4 video format. The operator, through the Smart Controller, can see the two spectra that present different views in each of the sensors that are in the gimbal of the drone. The DJI has four intelligent flight plan modes, including the "oblique plan", but the minimum flight height in this plan is 25 m, which implies a GSD as ground sampling distance of 0.76 cm/px for images acquired with the UAV. The other configured flight plan is a flight mode in which a minimum height of 12 m can be achieved by reducing the GSD to 0.43 cm/px.
In the experimental work, the TLS was used as an auxiliary record. With the appearance of new geomatic technology equipment such as Personal Laser Scanning (PLS), which simplifies data capture and the manageability of records in architectural spaces, it was decided to carry out a first survey through BLK360, which uses digitalisation technology in waveform (WFD) with a maximum scanning speed of 360,000 points/second.
The result obtained is a point cloud of 1,350,000 points, confined to the space of the archaeological site. The use of two massive data acquisition techniques allows us to guarantee that the entire surface of the dolmens is covered. For this reason, another part of the data capture techniques based on Structure from Motion (SfM) images is well known by the scientific and academic community [23] and allows for the construction of 3D models widely used in the field of archaeological architecture [4]. The photogrammetry workflow The operator, through the Smart Controller, can see the two spectra that present different views in each of the sensors that are in the gimbal of the drone. The DJI has four intelligent flight plan modes, including the "oblique plan", but the minimum flight height in this plan is 25 m, which implies a GSD as ground sampling distance of 0.76 cm/px for images acquired with the UAV. The other configured flight plan is a flight mode in which a minimum height of 12 m can be achieved by reducing the GSD to 0.43 cm/px.
In the experimental work, the TLS was used as an auxiliary record. With the appearance of new geomatic technology equipment such as Personal Laser Scanning (PLS), which simplifies data capture and the manageability of records in architectural spaces, it was decided to carry out a first survey through BLK360, which uses digitalisation technology in waveform (WFD) with a maximum scanning speed of 360,000 points/second.
The result obtained is a point cloud of 1,350,000 points, confined to the space of the archaeological site. The use of two massive data acquisition techniques allows us to guarantee that the entire surface of the dolmens is covered. For this reason, another part of the data capture techniques based on Structure from Motion (SfM) images is well known by the scientific and academic community [23] and allows for the construction of 3D models widely used in the field of archaeological architecture [4]. The photogrammetry workflow must manage a data plan based on measurement control points so that they are incorporated into Metashape's Agisoft software and, arriving at the intermediate postprocessing process of the images, include the measurements that scale the final model.
Control points are widely used for formal architectural control, in which an analytical model can be made after geometric analysis. Dos Santos et al. [24] revealed the complexity of the dimensional control of the control points in the georeferencing of a record with a terrestrial laser scanner through a scenario of interior and exterior spaces. And in this complexity, several authors [25] carried out a study that examined the georeferencing and geometric quality of a SfM record through multiple combinations of ground control points (GCPs). The conclusions established the proportionality of the control points in a uniformly distributed distribution. In our case study, the control points were taken by the LEICA Flexline TS02 total station with a precision of 2 mm [26]. The full description of the photogrammetry processing methodology used is described in previous works [27,28], leaving, in this case, the workflow well defined. Regarding the process carried out for the orientation of the base station, it was fixed by means of the local XYZ coordinates established at a single point of the project, which were taken as 100, 100 and 10 m, respectively.
Regarding the iPhone 13 Pro, its technical characteristics are the same as those of its predecessor, the 12 Pro version, as described by Monsalve et al. [29], who are a reliable source since Monsalve et al. quote their predecessors Luetzenburg et al. [12]. The LiDAR light detection and ranging sensor of this device operates at 8XX nm wavelengths, using vertical-cavity surface-emitting lasers to emit the laser, and the direct time of flight of the pulses is measured with a Single-Photon Avalanche Diode (SPAD). The iPhone 13 Pro captures a mesh through a mobile app such as "Polycam", which requires a paid subscription to export the data. Another well-known and, in this case, free app is "3d Scanner App". There are other applications, such as Scaniverse [30], but the most used are those specified above that export the data into various point cloud files such as .e57, .ply or .pts, among other files.
In the case of the "3d Scanner" application, it presents six capture modes in which LiDAR Advanced stands out as a prominent mode where the resolution of the point cloud can be set in a range of 5 to 50 millimetres. When measuring with Polycam, this software captures a record through a triangular mesh on the surface of the objects. In this case study, since the objective is to measure the different distances between targets, used as GCPs, the field on view is based on following the trajectory marked by the points. It must be remembered that the iPhone 13 Pro with a LiDAR sensor is equipped with a barometer, a three-axis gyroscope and an accelerometer, which the data capture software uses to determine the three spatial coordinates in the local reference system.

Limitations of the Methods Used
Stringent and reproducible measurements based on SfM must proceed in a flow order (established by researchers such as James et al. [31]) that determines the parameters to be taken into account, and any evaluation of the quality of the data must include a comparison with the coordinates of the verification points, the surfaces or the measures of length. A split test to determine systematic errors is also admissible. A split test aims to produce two sets of data with different shapes in the survey design. For example, it seems interesting that measurements carried out between the difference the GCPs and those obtained by a total station are analyzed in terms of the differential between the three coordinate axes (X, Y, Z), with which a data set could also be taken into account when comparing absolute distances between point units.
In the measurements carried out, there is an uncertainty added to the method that we could call residual uncertainty. This uncertainty stems from the absence of points recorded by the UAV on the crosshairs of the target and, likewise, from the absence of points on the crosshairs of the data captured with the iPhone 13 Pro scan and the TLS. In most cases, the best estimate is the arithmetic mean of the dimensional control of the control points in the georeferencing of a record with a terrestrial laser scanner through a scenario of interior and exterior spaces. And in this complexity, several authors [25] carried out a study that examined the georeferencing and geometric quality of a SfM record through multiple combinations of ground control points (GCPs). The conclusions established the proportionality of the control points in a uniformly distributed distribution. In our case study, the control points were taken by the LEICA Flexline TS02 total station with a precision of 2 mm [26]. The full description of the photogrammetry processing methodology used is described in previous works [27,28], leaving, in this case, the workflow well defined. Regarding the process carried out for the orientation of the base station, it was fixed by means of the local XYZ coordinates established at a single point of the project, which were taken as 100, 100 and 10 m, respectively.
Regarding the iPhone 13 Pro, its technical characteristics are the same as those of its predecessor, the 12 Pro version, as described by Monsalve et al. [29], who are a reliable source since Monsalve et al. quote their predecessors Luetzenburg et al. [12]. The LiDAR light detection and ranging sensor of this device operates at 8XX nm wavelengths, using vertical-cavity surface-emitting lasers to emit the laser, and the direct time of flight of the pulses is measured with a Single-Photon Avalanche Diode (SPAD). The iPhone 13 Pro captures a mesh through a mobile app such as "Polycam", which requires a paid subscription to export the data. Another well-known and, in this case, free app is "3d Scanner App". There are other applications, such as Scaniverse [30], but the most used are those specified above that export the data into various point cloud files such as .e57, .ply or .pts, among other files.
In the case of the "3d Scanner" application, it presents six capture modes in which LiDAR Advanced stands out as a prominent mode where the resolution of the point cloud can be set in a range of 5 to 50 millimetres. When measuring with Polycam, this software captures a record through a triangular mesh on the surface of the objects. In this case study, since the objective is to measure the different distances between targets, used as GCPs, the field on view is based on following the trajectory marked by the points. It must be remembered that the iPhone 13 Pro with a LiDAR sensor is equipped with a barometer, a three-axis gyroscope and an accelerometer, which the data capture software uses to determine the three spatial coordinates in the local reference system.

Limitations of the Methods Used
Stringent and reproducible measurements based on SfM must proceed in a flow order (established by researchers such as James et al. [31]) that determines the parameters to be taken into account, and any evaluation of the quality of the data must include a comparison with the coordinates of the verification points, the surfaces or the measures of length. A split test to determine systematic errors is also admissible. A split test aims to produce two sets of data with different shapes in the survey design. For example, it seems interesting that measurements carried out between the difference the GCPs and those obtained by a total station are analyzed in terms of the differential between the three coordinate axes (X, Y, Z), with which a data set could also be taken into account when comparing absolute distances between point units.
In the measurements carried out, there is an uncertainty added to the method that we could call residual uncertainty. This uncertainty stems from the absence of points recorded by the UAV on the crosshairs of the target and, likewise, from the absence of points on the crosshairs of the data captured with the iPhone 13 Pro scan and the TLS. In most cases, the best estimate is the arithmetic mean Ǭ of the n observations, according to Equation (1). of the n observations, according to Equation (1).
Control points are widely used for formal architectural control, in which an analytical model can be made after geometric analysis. Dos Santos et al. [24] revealed the complexity of the dimensional control of the control points in the georeferencing of a record with a terrestrial laser scanner through a scenario of interior and exterior spaces. And in this complexity, several authors [25] carried out a study that examined the georeferencing and geometric quality of a SfM record through multiple combinations of ground control points (GCPs). The conclusions established the proportionality of the control points in a uniformly distributed distribution. In our case study, the control points were taken by the LEICA Flexline TS02 total station with a precision of 2 mm [26]. The full description of the photogrammetry processing methodology used is described in previous works [27,28], leaving, in this case, the workflow well defined. Regarding the process carried out for the orientation of the base station, it was fixed by means of the local XYZ coordinates established at a single point of the project, which were taken as 100, 100 and 10 m, respectively.
Regarding the iPhone 13 Pro, its technical characteristics are the same as those of its predecessor, the 12 Pro version, as described by Monsalve et al. [29], who are a reliable source since Monsalve et al. quote their predecessors Luetzenburg et al. [12]. The LiDAR light detection and ranging sensor of this device operates at 8XX nm wavelengths, using vertical-cavity surface-emitting lasers to emit the laser, and the direct time of flight of the pulses is measured with a Single-Photon Avalanche Diode (SPAD). The iPhone 13 Pro captures a mesh through a mobile app such as "Polycam", which requires a paid subscription to export the data. Another well-known and, in this case, free app is "3d Scanner App". There are other applications, such as Scaniverse [30], but the most used are those specified above that export the data into various point cloud files such as .e57, .ply or .pts, among other files.
In the case of the "3d Scanner" application, it presents six capture modes in which LiDAR Advanced stands out as a prominent mode where the resolution of the point cloud can be set in a range of 5 to 50 millimetres. When measuring with Polycam, this software captures a record through a triangular mesh on the surface of the objects. In this case study, since the objective is to measure the different distances between targets, used as GCPs, the field on view is based on following the trajectory marked by the points. It must be remembered that the iPhone 13 Pro with a LiDAR sensor is equipped with a barometer, a three-axis gyroscope and an accelerometer, which the data capture software uses to determine the three spatial coordinates in the local reference system.

Limitations of the Methods Used
Stringent and reproducible measurements based on SfM must proceed in a flow order (established by researchers such as James et al. [31]) that determines the parameters to be taken into account, and any evaluation of the quality of the data must include a comparison with the coordinates of the verification points, the surfaces or the measures of length. A split test to determine systematic errors is also admissible. A split test aims to produce two sets of data with different shapes in the survey design. For example, it seems interesting that measurements carried out between the difference the GCPs and those obtained by a total station are analyzed in terms of the differential between the three coordinate axes (X, Y, Z), with which a data set could also be taken into account when comparing absolute distances between point units.
In the measurements carried out, there is an uncertainty added to the method that we could call residual uncertainty. This uncertainty stems from the absence of points recorded by the UAV on the crosshairs of the target and, likewise, from the absence of points on the crosshairs of the data captured with the iPhone 13 Pro scan and the TLS. In most cases, the best estimate is the arithmetic mean Ǭ of the n observations, according to Equation (1).
One of the parameters to take into account in the study of precision is the quality of the point cloud. This factor will be established by means of a verification procedure of the density of the cloud of points, measured in cm 2 , and could be determined from the data set obtained from one of the targets, which will be the reference of the measurement model. Costantino et al. [14] determined three parameters that influence the quality of the point cloud: (i) the variation in the surface, (ii) the planarity and (iii) the omnivariance. These three parameters determine aspects that can influence the construction of a 3D model, that is, the transformation of the cloud of points to a digital mesh, but for the precision of the data sets that we are going to observe, the most important thing is to know the omnivariance, that is, the degree of inhomogeneity of the cloud of points in the three dimensions (Equation (2)), or what we could specify as the density of the cloud of points.
O λ = (λ 1 ·λ 2 ·λ 3 ) 13 (2) In our experiment, we analyse each of the results that show the density of points in a fixed element that, as already mentioned, has as its object a rectangular wooden target of 948.75 cm 2 strategically placed on the surface of the ground as GCPs. The procedure consists of taking a series of 12 samples from the points closest to the filial cross of the target on the cloud of points and doing an analysis and average of the results obtained in each of the samples involved in the analysis. The sets of points of the aerial photogrammetry series at 24 m will be taken first, followed by 12 m, whose denominations of the set of points obtained are Dsfm_24 and Dsfm_12. In the same way, we work with the set of TLS points, and this procedure will also be carried out in the cloud of points from the scan using the sensor integrated in the mobile phone. From this analysis, the random error can be naturally inferred when taking the measurements of each of the samples.
To determine the density of the 3D point cloud, CloudCompare [32] software was used, in which elements are segmented to allow one to observe the factors already exposed by Rizali et al. [33], such as roughness, reflectivity strength, RGB colours and texture.
After the analysis, the following results were obtained (illustrated in Figure 5): The Dsfm_24 obtained a value of 0.3573 pt/cm 2 , the Dsfm-12 obtained a value of 1.2057 pt/cm 2 , the mobile phone sensor recorded 4.0737 pt /cm 2 and, finally, the TLS obtained a value of 7.9841 pt/cm 2 , obtained through surveys carried out in Dolmen no. 3. The values are reflected in red on the graph in Figure 5. Thus, from the values of the 12 samples that determine the distance between the points, the box plots of each of the selected samples were obtained. Of all the surveys, the one with the greatest dispersion was the one from the TLS, an the one that performed best was the record of the iPhone 13 Pro, although its uncertain expressed in the range of metric units was very similar to the values of the TLS. Clearl the density values obtained from aerial photogrammetry between 24 and 12 m are ve different, with the 24 m survey yielding greater uncertainty.
Point density is one of the parameters that determines the quality of the geometry a 3D mesh. In the modelling process, the number of points and their distribution can d termine the morphological parameters of the objects and their final result. The resultin Of all the surveys, the one with the greatest dispersion was the one from the TLS, and the one that performed best was the record of the iPhone 13 Pro, although its uncertainty expressed in the range of metric units was very similar to the values of the TLS. Clearly, Point density is one of the parameters that determines the quality of the geometry of a 3D mesh. In the modelling process, the number of points and their distribution can determine the morphological parameters of the objects and their final result. The resulting triangular mesh for 3D modelling varies according to the point density parameters, the modelling algorithms [34,35] and the shape of the object surface [23]. The dot density is a parameter that is measured according to the surface area. Therefore, the dot density of the dot array surface was evaluated per square centimetre, and the pattern taken was a rectangular bull's eye. The best way to check the spatial resolution of the point cloud, that is, the 3D Euclidean distance (see Equation (3)) between the closest points, is to sample the set of points. The distance between the points and the box plots of each of the selected samples was obtained as follows: where d E is the Euclidean distance between points in space and x, y and z are the Cartesian coordinates of those points.

Experimental Methodology and Results
In the first phase, an alignment is carried out with the Hand-Held Mobile Laser Scanner (HMLS), a portable mobile laser scanner that obtains a dense point cloud exported in a .ply file in the same coordinate system, and then a network of control points is measured with the total station (TS) using LEICA Flexline TS02 equipment. Therefore, throughout the work, different tests were used to obtain a set of data related to the comparison between the measurements of the total station (TS) and the set of points obtained through the iPhone 13 Scan Pro.
In the second phase of the experiment, a comparison is carried out to find out the geometric precision between the Hand-Held Mobile Laser Scanner and the SfM UAV using ground control points taken with the LEICA Flexline TS02 equipment. All of this estimation processing has been carried out using the open-source software CloudCompare. This analysis will determine the behaviour of the HMLS in the morphological structures of archaeological sites that involve organic forms and complex geometries.
Finally, in a case study different from the previous ones, we proceed to evaluate the precision of the HMSL that serves as a reference for capturing the control points on the ground for photogrammetry using a UAV. In this case, the point cloud record obtained from the TLS is taken as a reference. A comparison between the two of them will determine new data on the performance of the systems in the area of the archaeological site being studied. The 3D model of the terrestrial laser scanner is used as a reference to determine the coincidence points and their dispersion in the system used.

Analysis of the Distances between TS Targets and iPhone 13 Pro
For this section, three case studies were used, the first of which was Dolmen no. 3, in which five ground control points were used ( Figure 6). The data set was exported as a file .ply and TS data were exported in .txt. The data obtained from the deviation were taken to a graph to establish the existence of a direct correlation between the deviation and the measurements obtained in the images of Figure 7, 9 and 11. In this first data output, a total of ten samples were taken with a maximum distance of 9.7367 m and a minimum distance between targets of 3.0432 m. The absolute differential values obtained ranged from 20.47 cm (points 3-6) to 1.69 cm (points 1-6). This analysis is performed to be able to predict the changes and local trends observed after the field work is carried out.
which five ground control points were used ( Figure 6). The data set was exported as a file .ply and TS data were exported in .txt. The data obtained from the deviation were taken to a graph to establish the existence of a direct correlation between the deviation and the measurements obtained in the images of Figures 7,9 and 11. In this first data output, a total of ten samples were taken with a maximum distance of 9.7367 m and a minimum distance between targets of 3.0432 m. The absolute differential values obtained ranged from 20.47 cm (points 3-6) to 1.69 cm (points 1-6). This analysis is performed to be able to predict the changes and local trends observed after the field work is carried out.   The following work was carried out in an aerial photogrammetry of the Palma del Condado train station (Figure 8), in which seven control points were taken. In total, this analysis was carried out using twelve samples with a maximum distance of 14.1996 m and a minimum distance between targets of 2.8557 m. The absolute differential values obtained ranged from 1.98 cm (points 5-6) to 0.471 cm (points 4-5).
The last work (Figure 9) was carried out on an almost horizontal plane with little unevenness to check how the LiDAR sensor of the mobile phone works in a parking lot of the Higher Technical School of Building Engineering. Also in this work, twelve samples were used with a maximum distance of 23.3692 m and a minimum distance of 2.3680 m. The absolute differential values obtained ranged from 23.1087 cm (points 1-8) to 2.3560 cm (points 7-8). Therefore, the results of the following studies will determine the correla- The following work was carried out in an aerial photogrammetry of the Palma del Condado train station (Figure 8), in which seven control points were taken. In total, this analysis was carried out using twelve samples with a maximum distance of 14.1996 m and a minimum distance between targets of 2.8557 m. The absolute differential values obtained ranged from 1.98 cm (points 5-6) to 0.471 cm (points 4-5).  The following work was carried out in an aerial photogrammetry of the Palma del Condado train station (Figure 8), in which seven control points were taken. In total, this analysis was carried out using twelve samples with a maximum distance of 14.1996 m and a minimum distance between targets of 2.8557 m. The absolute differential values ob- The last work (Figure 9) was carried out on an almost horizontal plane with little unevenness to check how the LiDAR sensor of the mobile phone works in a parking lot of the Higher Technical School of Building Engineering. Also in this work, twelve samples were used with a maximum distance of 23.3692 m and a minimum distance of 2.3680 m. The absolute differential values obtained ranged from 23.1087 cm (points 1-8) to 2.3560 cm (points 7-8). Therefore, the results of the following studies will determine the correlation between the accuracy of the measurements and the distance between targets.       From the processed data set, the range of the deviation was determined in a box plot to find out the dispersion of the values obtained Figure 12. The distance between previously selected characteristic points represented the differences between the three-point clouds in the experiment carried out at the School, the Station, and Dolmen no. 3. The standard deviation results were lower in the Station experiment since a value of 0.02125 m was obtained, compared to the two values obtained in the School and Dolmen no. 3 experiments, which were 0.06921 and 0.06343 m, respectively. The results are shown in Table 1, configured by the parameters of the total number of samples in each case study and their corresponding values. Therefore, after analysing the experiment results for the three sets of points, we can observe that the school experiment was the one that obtained the worst results. However, the maximum length value of 23.3692 m between targets 1 to 8 shows that the accuracy of the values decreases at greater distances in this case. The applications that are developed in photogrammetry for surveying purposes, land construction or Scan-to-BIM in the field of archaeology have shown that GCPs greatly increase the accuracy of a dense point cloud [36]. Three targets with a black and white photogrammetric target are enough to establish a scale model, although it is already known that there are other parameters that can influence the quality of the point cloud in UAV photogrammetry, such as the flight height, the orientation of the camera in the flight plan process and the quality of the image that the drone gimbal has [37,38]. The objectives of this part of the work are to perform a scan with the LiDAR light detection and ranging sensor of the HMLS mobile phone and, through an analysis process with the algorithms available from CloudCompare, to make a comparison.

Analysis of the Deviation between Photogrammetry Point Clouds Using an UAV and the iPhone 13 Pro
Once the SfM point cloud data were obtained using the UAV, it was first scaled with the measurements taken in the Agisoft Metashape software and, subsequently, aligned with the registration data obtained from the light detection and ranging sensor (mobile phone LiDAR) by selecting common point pairs between these two-point clouds in CloudCompare. The three pairs of points (R0-A0, R1-A1 and R2-A2) were located at the centre of the ground control points, which, as already specified, used a 10 mm thick coloured chipboard base in black and white, yielding a final alignment RMS value of 0.00580 m.
Subsequently, the Iterative Closest Point (ICP) algorithm was applied to try to automatically optimize the alignment in order to obtain better results, but this process did not yield positive results, so the coupling between both point clouds was maintained by searching for pairs of adjacent points in the two data sets, and then the transformation parameters between them were calculated [39].
To determine the deviation between clouds of points, the C2C tool, which calculates the distance from each point, is applied. In this case, the point cloud of the aerial photogrammetry was taken as a reference to compare with the closest point in the point cloud of the iPhone 13 Pro.
In point cloud-to-point cloud comparison procedures, it is likely that there are outliers that are not visible to the naked eye but may affect the results. The best way to detect them is to conduct a first comparison analysis and then eliminate those points that can generate noise and alter the final results. These work filters in the cloud of points establish a certain adjustment that allows for the evaluation of the results in their entirety [40].
In aerial photogrammetry over forests or rocky massifs where there is abundant vegetation, the absence of records in forests is very frequent, especially between the ground and the base of the tree crown [41]. The occlusion of the points that record the tree trunk in the process of recording naridal photogrammetry means that the SfM model, using a UAV, cannot reconstruct its position. Therefore, in the open dolmen records in the field, occlusions may appear, as shown in Figure 13. The advantage of using a mobile phone LiDAR light detection and ranging sensor is the possibility of getting closer to the details of the morphology of the environment of the archaeological site, which, in the case of the SfM using a UAV, if the flight plan is not processed correctly, can affect the final results. The traditional way of characterising an archaeological complex through TLS involves certain difficulties [42], especially in horizontal planes of the ground, such as eliminating nondetailed parts, inaccessibility and bias in the collection of overhead data. photogrammetry was taken as a reference to compare with the closest point in the point cloud of the iPhone 13 Pro.
In point cloud-to-point cloud comparison procedures, it is likely that there are outliers that are not visible to the naked eye but may affect the results. The best way to detect them is to conduct a first comparison analysis and then eliminate those points that can generate noise and alter the final results. These work filters in the cloud of points establish a certain adjustment that allows for the evaluation of the results in their entirety [40].
In aerial photogrammetry over forests or rocky massifs where there is abundant vegetation, the absence of records in forests is very frequent, especially between the ground and the base of the tree crown [41]. The occlusion of the points that record the tree trunk in the process of recording naridal photogrammetry means that the SfM model, using a UAV, cannot reconstruct its position. Therefore, in the open dolmen records in the field, occlusions may appear, as shown in Figure 13. The advantage of using a mobile phone LiDAR light detection and ranging sensor is the possibility of getting closer to the details of the morphology of the environment of the archaeological site, which, in the case of the SfM using a UAV, if the flight plan is not processed correctly, can affect the final results. The traditional way of characterising an archaeological complex through TLS involves certain difficulties [42], especially in horizontal planes of the ground, such as eliminating nondetailed parts, inaccessibility and bias in the collection of overhead data. Once the differential between both point clouds has been processed, it is possible to extract information from matching points, where the CloudCompare software [39] determines the C2C deviation on a colour scale that represents the deviation of the compared model.
The distance between the point sets was calculated through a comparison between SfMUAV and HMLS, Figure 14 analysis, and Figure 15 histogram of the previous figure. The parameters that yield the results of the comparison process they expressed in Table  2: are the root mean square error (RMS), the standard deviation, the minimum distance, the maximum distance, the mean distance and the maximum estimation error. The deviation between sets of similar points has two characteristics: (i) the high presence of points in values close to zero in relation to other values that are especially far from the value 0.04 Once the differential between both point clouds has been processed, it is possible to extract information from matching points, where the CloudCompare software [39] determines the C2C deviation on a colour scale that represents the deviation of the compared model.
The distance between the point sets was calculated through a comparison between SfM UAV and HMLS, Figure 14 analysis, and Figure 15 histogram of the previous figure. The parameters that yield the results of the comparison process they expressed in Table 2: are the root mean square error (RMS), the standard deviation, the minimum distance, the maximum distance, the mean distance and the maximum estimation error. The deviation between sets of similar points has two characteristics: (i) the high presence of points in values close to zero in relation to other values that are especially far from the value 0.04 m and the rest of the intervals; and (ii) the high standard deviation according to the formulation of the points along those intervals (Equation (4)).
where n is the sample size, x i are the points in the intervals and x is the average sample value.  Figure 14.
where n is the sample size, are the points in the intervals and x � is the average sample value.  Figure 14.

Analysis of the Deviation between the Cloud of Points Obtained from SfM UAV Using GCP by HMLS and TLS
It is evident that the direct processing that is carried out can bring knowledge of the results closer, taking into account that the main objective of the research work is to know the limitations of the experimental methods in relation to the precision of the use of the iPhone 13 Pro as a measurement element to establish control points on the ground with the help of photogrammetry using unmanned aerial vehicles (UAV) in archaeological sites.
For this case, a point cloud was processed using UAV image capture with ground control points, in which the "Scala Bars" parameters were entered with the dimensions obtained by the LiDAR light detection and ranging sensor of the mobile phone. For both analyses of the sets of points, a partial part of the area surrounding the archaeological site was selected, segmenting both previously aligned sets.
According to studies carried out in other works, remote sensing-based image capture (SfM) using remote sensing is reaching great popularity in the characterisation of landscapes with centimetre-level precision. For example, in landslide studies [43] (where it was shown that the use of multi-rotor UAVs improves the point cloud by 16% compared to the use of fixed-wing UAVs) or in the heritage preservation field [44] (in a study on a Monastery, the quality parameters of the point cloud were evaluated using what the author calls the Remotely Piloted Aircraft System (RPAS) and the TLS). However, in the field of archaeology, there are studies [45] that only approach this investigation using SfM without UAV. Hence, it can be deduced from the previous approaches that the deviation analysis between SfM UAV by means of GCP by HMLS and TLS is innovative. Figure 16 shows the deviation between aerial photogrammetry using intelligent flight at 24 m and TLS developed on four BLK360 stands using LiDAR technology. The results show a standard deviation of 0.066 m and a maximum distance of 0.83 m.
the help of photogrammetry using unmanned aerial vehicles (UAV) in archaeological sites.
For this case, a point cloud was processed using UAV image capture with ground control points, in which the "Scala Bars" parameters were entered with the dimensions obtained by the LiDAR light detection and ranging sensor of the mobile phone. For both analyses of the sets of points, a partial part of the area surrounding the archaeological site was selected, segmenting both previously aligned sets.
According to studies carried out in other works, remote sensing-based image capture (SfM) using remote sensing is reaching great popularity in the characterisation of landscapes with centimetre-level precision. For example, in landslide studies [43] (where it was shown that the use of multi-rotor UAVs improves the point cloud by 16% compared to the use of fixed-wing UAVs) or in the heritage preservation field [44] (in a study on a Monastery, the quality parameters of the point cloud were evaluated using what the author calls the Remotely Piloted Aircraft System (RPAS) and the TLS). However, in the field of archaeology, there are studies [45] that only approach this investigation using SfM without UAV. Hence, it can be deduced from the previous approaches that the deviation analysis between SfMUAV by means of GCP by HMLS and TLS is innovative. Figure 16 shows the deviation between aerial photogrammetry using intelligent flight at 24 m and TLS developed on four BLK360 stands using LiDAR technology. The results show a standard deviation of 0.066 m and a maximum distance of 0.83 m.
Certainly, from the point of view of geometric understanding, the highest values are found in the registration points of the treetops, but these data are not representative. The mostly green area that corresponds to values of 0.31 m refers to occlusion points on the floor of the Dolmen recorded by the TLS. Figure 17 shows the histogram of the study, where 90% of all points have a coincidence of less than 0.032 m. The coincidence between Figures 16 and 17 is consistent since most of the points are in the consolidated values in the blue colour range. The quantitative data provided by the analysis software is represented in Table 3. Figure 16. Analysis of the difference between SfM UAV and GCP by HMLS and TLS. Result of applying the C2C algorithm. Distance distribution colour map. In the colour scale, red represents the maximum distance, expressed in metres, and blue represents the minimum distance.
Certainly, from the point of view of geometric understanding, the highest values are found in the registration points of the treetops, but these data are not representative. The mostly green area that corresponds to values of 0.31 m refers to occlusion points on the floor of the Dolmen recorded by the TLS. Figure 17 shows the histogram of the study, where 90% of all points have a coincidence of less than 0.032 m. The coincidence between Figures 16 and 17 is consistent since most of the points are in the consolidated values in the blue colour range. The quantitative data provided by the analysis software is represented in Table 3.
Analysing the morphological structure of archaeological sites and the geomorphology of the environment requires accurate topographic information, but when deposits are inaccessible, operators must resort to technologies that revolutionise the current digitalisation market. The presence of trees and shrubs should not be considered an obstacle to characterising these structures. In the processing and analysis stages, it was observed that the SfM UAV point cloud at 12 m behaved differently from the SfM UAV point cloud at 24 m. For this reason, it was decided to open a new analysis section.  Figure 16. Analysis of the difference between SfMUAV and GCP by HMLS and TLS. Result of applying the C2C algorithm. Distance distribution colour map. In the colour scale, red represents the maximum distance, expressed in metres, and blue represents the minimum distance. Figure 17. Histogram of Figure 16. Histogram units: metres (X axis) and number of points (Y axis). Table 3. Table of deviation between the SfMUAV and TLS sets of points from the survey in Figure 16.

Deviation (σ) (m)
Min Analysing the morphological structure of archaeological sites and the geomorphology of the environment requires accurate topographic information, but when deposits are inaccessible, operators must resort to technologies that revolutionise the current digitalisation market. The presence of trees and shrubs should not be considered an obstacle to characterising these structures. In the processing and analysis stages, it was observed that the SfMUAV point cloud at 12 m behaved differently from the SfMUAV point cloud at 24 m. For this reason, it was decided to open a new analysis section.

Analysis of the Deviation between Photogrammetry Point Clouds Using a UAV in Intelligent Flight at a Height of 12 m and 24 m
Very recent studies emphasise the study of processing through image capture records using UAVs at different distances on the Z axis. These studies have focused on forestry and agriculture to analyse the morphology of trees, essentially volume analysis. In other words, most studies [46] fall within the framework of forestry analysis and geology studies [25], in which different flight scenarios are exposed with different angles of the UAV camera in complex landscapes with high tree density and structured layered morphology, while other works fall within the framework of knowledge about the accuracy and reconstruction of extensive topographical environments [47]. Thus, it seems interesting that, in the context of archaeological sites, the difficulties that we can find in processing field work in different intelligent flight plans are exposed. For this, a study including the processing within the context of the dolmen was carried out, taking the set of SfMUAV points at 24 and 12 m, as can be determined in Figure 18.
The data revealed that 71% of the points are coincident at values lower than 0.122 m. This is an important deviation, considering that many points are in the green zone with deviations of 40 cm. This analysis can be seen in the histogram of Figure 19.   Figure 16.

Analysis of the Deviation between Photogrammetry Point Clouds Using a UAV in Intelligent Flight at a Height of 12 m and 24 m
Very recent studies emphasise the study of processing through image capture records using UAVs at different distances on the Z axis. These studies have focused on forestry and agriculture to analyse the morphology of trees, essentially volume analysis. In other words, most studies [46] fall within the framework of forestry analysis and geology studies [25], in which different flight scenarios are exposed with different angles of the UAV camera in complex landscapes with high tree density and structured layered morphology, while other works fall within the framework of knowledge about the accuracy and reconstruction of extensive topographical environments [47]. Thus, it seems interesting that, in the context of archaeological sites, the difficulties that we can find in processing field work in different intelligent flight plans are exposed. For this, a study including the processing within the context of the dolmen was carried out, taking the set of SfM UAV points at 24 and 12 m, as can be determined in Figure 18.  The data revealed that 71% of the points are coincident at values lower than 0.122 m. This is an important deviation, considering that many points are in the green zone with deviations of 40 cm. This analysis can be seen in the histogram of Figure 19. Figure 18. Analysis of the difference between SfMUAV at 12 and 24 m. Result of applying the C2C algorithm. Distance distribution colour map. In the colour scale, red represents the maximum distance, expressed in metres, and blue represents the minimum distance. The presence of deviations around forty centimetres that represent the majority of the treetops, shrubs and, especially, the upper part of the orthostats, that is, the stone blocks that make up the course of the wall of the dolmen, can be observed. From this line, it can be inferred that there is a presence of systematic vertical deviations. Such errors, which in typical UAV usage data are expressed as a vertical "dome" of the surface, are the result of a combination of nearly parallel imaging directions and an imprecise correction for the distortion of the radial lens [48].
The error in the systematic vertical deviations can be corrected with the additional capture of oblique images. By separating the context from the site and focusing on the dolmen itself, the evidence for deviations is more significant, and a greater appreciation of the results can be obtained. In this second data output, a segmentation of a rectangular prism is carried out on Dolmen no. 3, in which the set of points of the 24 m SfMUAV is taken as a reference.
In Figure 20, it can be seen that the systematic vertical deviations referred to in the studies by James and Robson [48] appear in red. The presence of deviations around forty centimetres that represent the majority of the treetops, shrubs and, especially, the upper part of the orthostats, that is, the stone blocks that make up the course of the wall of the dolmen, can be observed. From this line, it can be inferred that there is a presence of systematic vertical deviations. Such errors, which in typical UAV usage data are expressed as a vertical "dome" of the surface, are the result of a combination of nearly parallel imaging directions and an imprecise correction for the distortion of the radial lens [48].
The error in the systematic vertical deviations can be corrected with the additional capture of oblique images. By separating the context from the site and focusing on the dolmen itself, the evidence for deviations is more significant, and a greater appreciation of the results can be obtained. In this second data output, a segmentation of a rectangular prism is carried out on Dolmen no. 3, in which the set of points of the 24 m SfM UAV is taken as a reference.
In Figure 20, it can be seen that the systematic vertical deviations referred to in the studies by James and Robson [48] appear in red.  The resulting histogram ( Figure 21) shows the number of points per distance. This is counted in absolute values and in metric units. The distribution parameters corresponding to the maximum, minimum, standard deviation and mean values are not really significant here, but the most relevant piece of data is that the area of points in red in the range of colours has a range between 0.15 and 0.25 m. Figure 20. Analysis of the difference between SfMUAV at 12 and 24 m. Result of applying the C2C algorithm. Distance distribution colour map. In the colour scale, red represents the maximum distance, expressed in metres.
The resulting histogram ( Figure 21) shows the number of points per distance. This is counted in absolute values and in metric units. The distribution parameters corresponding to the maximum, minimum, standard deviation and mean values are not really significant here, but the most relevant piece of data is that the area of points in red in the range of colours has a range between 0.15 and 0.25 m.

Scanning for BIM in an Archaeological Context
The use of HBIM for archaeology arises in analogy and continuity with research and applications on BIM for historical buildings, a subject on which the scientific community has been debating for many years, with different regulatory repercussions from country to country. Nevertheless, there are several specificities in the use of BIM models for archaeology that call for a separate discussion.

Scanning for BIM in an Archaeological Context
The use of HBIM for archaeology arises in analogy and continuity with research and applications on BIM for historical buildings, a subject on which the scientific community has been debating for many years, with different regulatory repercussions from country to country. Nevertheless, there are several specificities in the use of BIM models for archaeology that call for a separate discussion.
As it is well known, BIM is, according to ISO 19650, a series of standards whose content is intended to support agents and operators that participate in any phase of the building cycle or infrastructure pathways through BIM methodology [49], which includes a part of process management. In the case of archaeology, this tool could help not only in the management of all the activities related to preservation, conservation and restoration but also in the "interpretation, presentation, access and public use of the material remains of the past" through virtual archaeology. In the international document "Principles of Seville" [50], aimed at providing practical applications of the principles of the "London Charter" [51], virtual archaeology is defined as "the scientific discipline that seeks to research and develop ways of using computer-based visualisations for the comprehensive management of archaeological heritage" [52], including "virtual restoration, anastylosys, reconstruction, and recreation" [53,54]. To this aim, "inventories, surveys, excavation work, documentation, and research" become fundamental knowledge moments to create two-and three-dimensional representations of the real objects, guaranteeing their "authenticity", in the meaning expressed in the Sevilla Charter, to be obtained through a rigorous and faithful process of representation of their real consistency, also using parametric software, which is fundamental to being true to the geometric and morphological complexity that characterises the stratified ancient heterogeneous structures [55]. All this makes the construction of an HBIM particularly complex [56]. However, research and experimentation in the field of three-dimensional and parametric modelling of these kinds of objects, especially through Scan-to-BIM and photogrammetric processes, have led to the identification of modelling workflows shared by the scientific community [57].
From this point of view, the use of HBIM models represents a valid aid in terms of integrated management, in a single three-dimensional work environment, of different types of data: geometric, stratigraphic, historical, material, environmental, geological, etc. The HBIM model could be considered a 3D database that can be queried in different ways [58] during the design, construction and valorisation phases. The possibility of linking parametric objects, more or less simplified according to different levels of geometric and informative detail, depending on the purpose, with relational databases to be implemented and updated continuously, guarantees that we have a tool for various users to exchange and share information, from archaeologists to conservation offices to scholars. Operators involved in a restoration project can consult the model to access related documents, such as degradation mapping, to use them directly on the archaeological site and plan the interventions to be carried out [59]. In the same way, the model can be consulted by experts interested in accessing iconographic material or archival sources inserted within the HBIM, allowing for new spatial correlations to be established and not only documentary ones. At the same time, there is the possibility of the models being consulted by tourists or simply curious people who want to live immersive experiences through new modes of digital narration, using, for example, augmented and virtual reality [60]. Finally, an HBIM could represent a valid support for public administrations, taking advantage of a valid and reliable verification and control system that, in the long term, provides a mapping of works on the territory and their state of conservation [61], thanks to the possibility of BIM/GIS integration.
It is for this reason that, in recent years, many researchers, at an international level, have dealt with different topics related to BIM for archaeology. The most frequent ones tackle the issues of parametric three-dimensional modelling [62,63], starting from complex data surveying [64], ontologies and sharable vocabularies [65], interoperability between different software [66] and stratigraphic analysis [67].
Currently, there is an evaluation method that combines the creation of a geometric model from a cloud of points with the set of points obtained from area photogrammetry. This method, known as C2M and implemented through CloudCompare, has been used by several researchers [42,68,69] in the context of Cultural Heritage. The direct geometric model is the result of the process of creating static 3D models with simple or complex surfaces [70]. In the field of archaeology, most of the time complex or very complex surfaces are involved, which means that field tests for their evaluation require more complete processes.

Discussion
New data capture systems for three-dimensional representation are changing the way geomatics operators and topography experts work. Newer and easier-to-use tools appear, but they can call into question their precision compared to other, more experimental equipment.
The incorporation of LiDAR techniques, such as sensors integrated into mobile phones, is revolutionising the world of 3D scanning. One of the most positive aspects in terms of LiDAR smartphones is their working conditions, among which one can find the possibility of scanning different surfaces from a medium-range distance, including materials in the same scenario, and the processing of data in the "app" of the system lasting a few minutes.
But these record systems are not exempt from serious problems. For example, Costantino et al. [14], for the evaluation of an iPhone 12 Pro, showed that small-scale variations in objects of up to 2 cm of difference in different materials and the architectural space used with somewhat larger dimensions, such as the transept of a church, were detected with precisions of the order of 1 to 3 cm of variation. Some inertial navigation system drift and loss of planarity issues were also mentioned.
Luetzenburg et al. [12], who made measurements of a 130 m by 15 m wide cliff, compared the SfM-MVS with the iPhone 12 Pro point cloud, and the results obtained showed that for 80% of all recorded points, the maximum distance was 15 cm. Through the analysis of the graph of the percentage of precision with respect to the three coordinate axes, it should be emphasised in this work that Z (height) yields quite imprecise values.
The comparison between the integrated sensor of an iPad Pro and TLS was measured by Spreafico et al. [3], who determined that the precision was 2 cm and the accuracy was 4 cm.
One of the most interesting studies developed in the field of interior spaces is the work of Chase et al. [15], where the results in workspaces with a length between 1 and 13 m determined an average error of 3 cm in coordinate (X) and 7 mm in coordinate (Z) in fifteen different shots.
The main difference in the results between Luetzenburg et al. and Chase et al. [15] is the variability of the precision in coordinate (Z), and we believe that Chase's work occurs in conditions where the surface is slightly flat, so the iPhone sensor behaves extremely precisely in these measurements. On the other hand, when the morphology is abrupt and highly variable, precision changes are important.
One of the study samples that made him reevaluate the behaviour of the mobile phone sensor was the result that he showed in the work of Chase et al. [15]. Thus, we selected the results from Table 3 of his work and proceeded to graph them, observing the existence of a linear correlation between the equipment error and the distance. This process was carried out by validating the total station equipment as a reference on the iPhone. This graph is shown in Figure 22, in which it can be seen that there is no correlation between distance and relative precision, so it is determined that the differences may be related to other factors.
We can assume that the relationship between the distance between target points and their precision can be maintained in increasing order and thus predict the behaviour of the deviation between the values of length and precision of the LiDAR system from the linear regression line.  Table 3 of the work of Chase et al. [15].
On the basis of these studies, we carried out a section on the uncertainty associated with the measurement method, which has been called "residual uncertainty", by studying the density of points on the object to be measured. Interestingly, the method that registered the highest density of points was TLS, followed by the LiDAR light detection and ranging sensor of the mobile phone. Regarding the dispersion of the points in the measurements, the iPhone behaved quite stably; that is, in short measurements, the set of points is very distributive, reaching uncertainty values of a range of 0.005 m, that is, values of 5 millimetres.
Next, we generated an experimental methodology consisting of four sections. The first one had to do with the errors of the LiDAR light detection and ranging sensor of the  Table 3 of the work of Chase et al. [15].
On the basis of these studies, we carried out a section on the uncertainty associated with the measurement method, which has been called "residual uncertainty", by studying the density of points on the object to be measured. Interestingly, the method that registered the highest density of points was TLS, followed by the LiDAR light detection and ranging sensor of the mobile phone. Regarding the dispersion of the points in the measurements, the iPhone behaved quite stably; that is, in short measurements, the set of points is very distributive, reaching uncertainty values of a range of 0.005 m, that is, values of 5 millimetres.
Next, we generated an experimental methodology consisting of four sections. The first one had to do with the errors of the LiDAR light detection and ranging sensor of the mobile phone compared to the data obtained by the total station (TS) using the LEICA Flexline TS02 equipment. On the basis of this analysis, we proceeded to determine the existence of a linear correlation between the distance of the GCPs and the deviation that occurs with respect to the measurements with the TS. The statistical results of the comparison revealed that there is a slight relationship between the parameters analysed (R2) in Figures 7 and 8, with very low values of 0.04 and 0.09 m, respectively. The graph that presents some correlation is Figure 11, with values of 0.24.
One must think that the differences between the graphs of Chase et al. [15] and those elaborated in this investigation are justified in that the author takes distances as absolute values. From the subsequent analysis (Figure 12), there are three cases that comprise a standard deviation between 0.021 and 0.069 m. In the case of the station's field work study, the deviation value was quite satisfactory with only a 21 mm difference, but the rest of the studies presented highly significant variations.
Within the development of the field work, the case studies of the School and the Station were studied on slightly flat surfaces, which was not the case in the Dolmen case study, where the Z axis was more variable. But it is not ruled out that the analysis surface in a variable geomorphological environment increases the error value of the LiDAR light detection and ranging sensor of the mobile phone.
In the next study, an analysis of the difference between photogrammetry using UAV and HMLS was carried out. The results showed a standard deviation of 0.0889 m. However, this deviation was estimated to be due to non-alignment values since most of the scattered points were at the end of the CGPs, as shown in Figure 14. A significant value was the high presence of points in values close to 0 in relation to other values that were especially far from the value of 0.04 m and the rest of the intervals; that is to say, most of the points presented a differential of 4 centimetres. This would comprise an average interval of the results obtained in the previous investigation [15].
Regarding the analysis of the difference between SfM UAV and GCP by HMLS and TLS, there was a good degree of agreement. In other words, after meeting the requirements of the objectives set in this part of the study, the results were positive. In this case, the standard deviation parameter was lower than in the previous one, with a result of 0.0666 m. A significant value was the high presence of points at values close to 0, in which 90% of all points had a coincidence of less than 0.032 m. These values were lower than in previous studies. Therefore, when we process photogrammetry using a UAV at a height of 24 m with intelligent flight and nadiral and oblique images and use the iPhone 13 Pro to measure ground control points, we obtain an error against the TLS using the BLK360 32 mm. Thus, these are some results that can be assumed for Scan-to-BIM environments.
Also, the use of aerial photogrammetry at different levels and ways of proceeding is essential to knowing how an optimal 3D model can be obtained. Thus, an analysis of the difference between SfM UAV at 12 and 24 m with the intelligent flights provided by the drone used was determined with the next field work test. In this work, it was possible to determine the presence of systematic vertical deviations, in addition to the fact that these errors, in the typical data of the use of the UAV, are expressed as a vertical "dome" of the surface in the absence of oblique images. In this document, it was possible to quantify the deviation value of the orthostats, establishing a range between 0.15 and 0.25 m.
To provide comprehensive knowledge in the context of Scan-to-BIM applied to archaeological sites, all results have been graphed in Figure 23, where the representation in the range of millimetres of the slight accuracy (LOA) of the independent standard deviation is combined for each method used to acquire the rank of these measurements in USIBD [71]. The green points are sets of points randomly established to represent the record of a surface from a point cloud, where the actual surface is established as an imaginary line close to a value of zero. The second overlaid graph represents the mean deviation of four researchers and the results of the four experiments conducted in this work in a range expressed in metres. The equivalence is evident between the coordinates of the ordinate axes in red and black, although both are expressed in different units. It should be mentioned that the difference between the last survey and the rest of the tests is that each and every one of the tests was conducted with the LiDAR light detection and ranging sensor equipment of the mobile phone except for Survey_4, which determines the deviation between photogrammetry areas with heights of 12 and 24 m. . Average values with their errors in a graph with LOA range in the USIBD guide comparing the authors: Costantino et al. [14]; Luetzenburg et al. [12]; Spreafico et al. [3] and Chase et al. [15] with surveys carried out.

Conclusions
Disruptive techniques of massive data acquisition are presented as a new opportunity for work that completes the records of 3D digitisation and the implementation of Scan-to-BIM nowadays. Each method requires different data processing and is generally associated with specific software that processes the data until the point cloud is ready to take it to BIM. With this study, we have had the possibility of knowing the deviations shown by the use of the LiDAR light detection and ranging sensor on the iPhone 13 Pro.
Thus, in the data processing phase of this work, four comparison studies were carried out, along with a section in which an approximation to the Scan-to-BIM of archaeological sites was investigated. Through the study of aerial photogrammetry data obtained through two intelligent flight plans, the deviations generated by the measurements of a mobile phone could be compared, and its applicability to capture data from GCPs (ground control points) that serve to improve photogrammetry was assessed.
It must be taken into account that the analysis of the morphological structure of archaeological sites and their geomorphology is carried out in a complex context, such as on open sites where abrupt vegetation and groves abound. Most of these deposits are difficult to access, so having a drone and a mobile phone with these characteristics in one's backpack means that any expert can carry out 3D survey work with certain guarantees of Figure 23. Average values with their errors in a graph with LOA range in the USIBD guide comparing the authors: Costantino et al. [14]; Luetzenburg et al. [12]; Spreafico et al. [3] and Chase et al. [15] with surveys carried out.

Conclusions
Disruptive techniques of massive data acquisition are presented as a new opportunity for work that completes the records of 3D digitisation and the implementation of Scan-to-BIM nowadays. Each method requires different data processing and is generally associated with specific software that processes the data until the point cloud is ready to take it to BIM. With this study, we have had the possibility of knowing the deviations shown by the use of the LiDAR light detection and ranging sensor on the iPhone 13 Pro.
Thus, in the data processing phase of this work, four comparison studies were carried out, along with a section in which an approximation to the Scan-to-BIM of archaeological sites was investigated. Through the study of aerial photogrammetry data obtained through two intelligent flight plans, the deviations generated by the measurements of a mobile phone could be compared, and its applicability to capture data from GCPs (ground control points) that serve to improve photogrammetry was assessed.
It must be taken into account that the analysis of the morphological structure of archaeological sites and their geomorphology is carried out in a complex context, such as on open sites where abrupt vegetation and groves abound. Most of these deposits are difficult to access, so having a drone and a mobile phone with these characteristics in one's backpack means that any expert can carry out 3D survey work with certain guarantees of precision.
With the set of data obtained, data on the precision of this equipment can be inferred from these investigations, although the door is left open for other investigations that relate to the capture of records made with a LiDAR light detection and ranging sensor of a mobile phone, correlating the results with other parameters to clarify questions about how to reach precisions close to an LOA of 50 (range values between 0 and 1 millimetre).
It must be recognised that applying BIM to heritage is a form of digital registration that is open to different administrations and that, sometimes, the registration can be optimal in its performance. The path that leads us towards digitisation and the construction information model in historical buildings or archaeological sites consists of innumerable stages, one of the first being the generation of geospatial data. Knowing the requirements and specifications of the levels of precision allows professionals in the Architecture, Engineering, Construction and Ownership (AECO) industry to clearly articulate the objectives of the works. Therefore, precisely understanding the scope of the impact of component errors improves the results of the applicability of BIM.
In the future, through experimental work related to point cloud mapping in archaeological pieces and altarpieces, it will be possible to study how certain deformations produced by the LiDAR light detection and ranging sensor of the mobile phone influence linear structures, quantifying the deformation according to image capture.
Finally, it is also considered necessary to indicate that one of the most important objectives of this work has been fulfilled since a new scenario of possibilities for the applicability of new techniques in the management of the acquisition of massive records to produce increasingly accurate 3D models has opened up.