Integrated Geomatic Techniques for Georeferencing and Reconstructing the Position of Underground Archaeological Sites: The Case Study of the Augustus Sundial (Rome)

: A large part of the archaeological remains still to be discovered and excavated are not in remote and depopulated areas of the earth but are often beneath urban centres that have buried them with centuries of debris and later constructions. Excavating in these contexts is much more complex than digging in rural or sparsely inhabited areas instrument itself. The aim of this work is to study and thus define the combination and integration of existing geomatic techniques for this specific field of application.


Introduction
Some of the most important archaeological remains in the world are found underground in urban centres with a high population density where the succession of demolitions and reconstructions in subsequent periods has caused the progressive burial of archaeological remains. For example, in the centre of Rome, a specific rule of the Papal State imposed for centuries determined that new buildings had to be built on Roman foundations where possible, and over the centuries this has further complicated the present investigation and recovery of Roman structures. This stratification of building on buildings of different eras has, however, not only occurred in Rome or Italy but also in many other European and non-European countries, such as Egypt, Peru and Mexico, to name but a few; therefore, specific procedures are of great interest. These procedures could allow for "surgical" interventions as circumscribed and limited as possible in urban areas. This field of study, which starts from the analysis and comparison of existing technologies and methodologies, could lead to the definition of methods and/or specifications of any dedicated devices.
This paper illustrates the methodologies by which, with an integrated use of Geomatics techniques, it is possible to reconstruct the orientation and the absolute position in a global datum of poorly lit and difficult to access hypogeum spaces. The techniques are illustrated and discussed in general, analysing and evaluating the obtainable accuracy. Reference is then made to a case of application consisting of the survey of a portion of the so-called "Augustus sundial" in the subsoil of Rome.
There are two specific purposes of the sundial relief: 1. To support with great geometric accuracy the evaluation of the feasibility of continuing the excavations, considering that a mere 5 mm inaccuracy in the positions of the visible points (that presently stretches for around seven metres) can lead to an uncertainty of about 0.4 metres in the most distant part still to be discovered, which would be tens of metres from the part already excavated, so lower accuracy could be inadequate; 2. To allow the reconstruction of the observable geometric characteristics of the sundial in order to evaluate the possibility of a rigorous mathematical and astronomical reconstruction of the sundial's functioning during a subsequent phase of research.
More generally, the aim of this research is to evaluate various geomatic techniques for this specific field of application, among which the comparison between laser scanning techniques (more consolidated) and photogrammetric techniques based on the more recent and less consolidated "structure from motion" (SFM) approach that is of particular interest [1].
The survey of an archaeological artifact such as the Sundial of Augustus, with its historical, archaeological and scientific importance unique in the world and in history, requires results that are at least as accurate as those that were used in its realisation. The task is obviously much more complex today, since the sundial itself is no longer on the walking surface of an unbuilt area, as it was at the time of its realisation, but about seven metres below the current ground level in a densely urbanised area [2,3].
The problem of georeferencing archaeological sites in an absolute reference system, for the purpose of reconstructing important parameters such as the orientation of roads or construction axes with respect to the north, is normally dealt with through GNSS measures for the framing phase, then integrated for the detail survey by other geomatics techniques (e.g., total stations, digital photogrammetry or TLS-Terrestrial Laser Scanning) [4][5][6]. Other methodologies could be used, such as handheld mobile laser [7,8], but the achievable accuracy cannot be adequate to our aim.
For the GNSS geodetic network, the quality of which is decisive because it is the first phase of the survey, relative positioning measurements are carried out with respect to external framing networks [9] and/or local known points, performed with post-processing or real time techniques, depending on the degree of accuracy to be achieved.
The detailed survey can be performed using several techniques to obtain precise and reliable results. In this work two different approaches were carried out and compared; moreover, the final results can be integrated [10][11][12]. The TLS has the capacity to scan and map the environment in a significantly shorter time, especially if compared to the traditional archaeological survey techniques [13,14]. On the other hand, photogrammetry is widely used for modelling and documenting large and complex archaeological sites [15,16], even when located in critical situations such as in underground [17] and underwater environments [18].
The connection to external reference networks and known vertices is relatively easy when the archaeological site is outdoors, with a good visibility of the satellites and the possibility of carrying out GNSS measurements directly on the site. The operations are considerably more complex when the site is made up of underground compartments, in which GNSS measurements cannot be performed directly. In these cases, it is necessary to provide a two-level framing network: a GNSS geodetic network of absolute georeferencing outside, and a secondary network measured with topographic instruments (e.g., total stations) to connect the underground spaces to the GNSS external framing network [19].
If the underground-external connection is long and complex, due to the conformation of the underground spaces to be surveyed and the other rooms and corridors that give access to them, the geodetic network structured on two levels will be necessarily very complex. Such a complex network, often comprising a large number of vertices, requires a very careful design and execution in order to limit the propagation of measurement errors, which can lead to insufficient accuracy for the survey purposes [20].
The article will first describe the historical, cultural and scientific importance of the Augustus Sundial, after which the two survey campaigns and the specific methodologies used will be described and the estimated error parameters will be presented; the surveys will then be compared with each other and the final deviations described.

Uniqueness and Historic Introduction of the Site
This important archaeological element, linked to the Roman identity and its classical past, was built in the area of Campus Martius in 9 B.C., the area chosen for the construction of most of the urban plans of the imperial age [21]. Historically an unhealthy swamp, often subject to the Tiber flooding or used for military exercises and electoral rallies, Campus Martius quickly became the privileged area for new monuments celebrating the apotheosis of emperors and their families [22]. Augustus' mausoleum, the Ara Pacis and the solar meridian together with the obelisk, formed a monumental complex with rather a great symbolic value for those times (Strab. V 3,8).
The solar meridian should have proved the validity of the edict that Augustus himself promulgated to remedy the disastrous state of the Roman calendar. The Numan calendar was the first Roman calendar we have certain sources of. It was a lunisolar calendar that did not follow the Sun nor the Moon and it lost, on average, one day each year [23]. Moreover, due to the incorrect application of the intercalations, whose management was entrusted to the Pontifex Maximus, in 46 B.C. it anticipated the solar year by 90 days, a whole season. Caesar achieved the realignment by adding 67 days to the Numan intercalation. The annus confusionis, as 46 B.C. was renamed, was, therefore, composed of 15 months, for a total of 445 days (Suet. Iul. 40.1).
Caesar then opted for a solar calendar, the Annus Iulianus, composed of 365 days which were divided into 12 conventional months. To link the new calendar to the tropical year, whose duration had been estimated by Egyptian astronomers to be 365 days and 1/4, it was decided to intercalate one day every four years. With the doubling of February 24, the dies bis sextus-which gave the leap year its name-was obtained [24].
Despite the simplicity of this solution the Pontifex Maximus, which was already responsible for the ruin of the Numan calendar, it misinterpreted the Julian reform: the leap day was added every three years instead of four years. The error lasted from 45 to 9 B.C., the year of Augustus's solution: during this time there were 12 three-year intercalations, instead of nine four-year intercalations, which caused the loss of three days [24].
In Rome, the mistake was noticed, but they could only solve the problem after 36 years. For an ideal revision, they were in fact waiting for a year in which the "three-year intercalation", which was mistakenly used until then, would match with the "four-year intercalation", the solution provided by Caesar. The concomitance occurred three times: in 33 B.C. during the fratricidal war between Antonio and Octavian, when Lepidus, the third triumvirus accused of double-dealing, was exiled with the title of Pontifex Maximus. He still held this function in 21 B.C., when Augustus already controlled Roman politics. The perfect time to intervene in the calendar occurred in 9 B.C. because, with the death of Lepidus three years earlier, the title of Pontifex Maximus passed into the hands of Augustus, who now cumulated both political and religious powers. Legitimised by the sacred investiture, Augustus promulgated an edict that corrected the erroneous application of the Julian reform. He ordered the suspension of the leap year intercalations three times, a solution that guaranteed the recovery of the three lost days. The leap day was omitted in 5 and 1 B.C., then in 4 A.D.: with the intercalation of 8 A.D., the Julian reform finally came into effect [23]. Roman mathematicians had the task of devising a new type of instrument that would demonstrate to the world the reliability of the Augustan edict. In Naturalis historia, Pliny describes that on the solar meridian built in the Campus Martius, the shadow cast by an obelisk would indicate the succession of days, settling in correspondence of each single rule at noon each day. The shadow would eventually attest its return to the same point after a leap cycle. It would reach its maximum length at the winter solstice, when the Sun assumes, compared to the other half-days of the year, the least possible declination. On the contrary, it would gradually shorten, becoming the shortest at noon of the summer solstice (Plin. NH. XXXVI 72-73).
In addition to this, Pliny refers that the sundial unexpectedly stopped working just a few years after it was built, perhaps due to an earthquake or a flood of the Tiber. In any case it would seem that Pliny did not see the sundial himself, because he does not know the reason why it stopped working. Heslin suggests that the sundial floor, found covered with a waterproof layer of earthenware, began to serve as a water cistern [25].
As for the obelisk described by Pliny and used as a gnomon, it was brought to Rome in 10 B.C., together with a second obelisk, now visible in Piazza del Popolo (Rome). Moved from its original site, Heliopolis (Egypt), the obelisk had been carved from a red granite block almost 22 metres high, at the behest of Pharaoh Psamtic II, in the VI century B.C.. Once transported in Rome, its base was enriched by two identical inscriptions that highlight Augustus' devotion to the figure of the Sun god. On its top, a metal sphere was placed to direct the shadows on the correct rule of the solar meridian. After Pliny's testimony, it is believed that the obelisk fell in 1084, following the days of sacking that the Norman troops, led by Robert Guiscard, perpetrated in the streets of Rome.
Starting from the modern age, both the sundial and the obelisk are linked to the history of the religious complex of San Lorenzo in Lucina. From some notes left by Antonio Lelio we have knowledge that the obelisk was first glimpsed in 1512 in the house of a barber during the realisation of a sewer. In 1587, following the request made by Sixtus V for an archaeological investigation, Vacca and Lauro reported that the obelisk had been broken into several pieces and ruined by fire. In 1748, taking advantage of the demolition of some houses overlooking Largo dell'Impresa, Benedict XIV finally ordered the recovering of the obelisk. In order not to lose the memory of the place where it was found, a marble plaque was affixed above the portal of the new building. On that occasion, the architect Tommaso De Marchis made a detailed description of the monolith. The report is useful to reconstruct the obelisk height from the ground, without the metal gnomon, and its original position.
The obelisk was restored with the pieces of the Antonine column; then, it was raised and placed in the square of Montecitorio in 1792.
The history of the solar meridian finds is different. Many scholars of the 16th century (reinterpreted and spread a bad transcription of some notes written by a student of the humanist Pomponio Leto in 1484. This misunderstanding fuelled for centuries the theory of the existence of a great sundial that would have marked the daily hours and whose floor would have extended for most of the Campus Martius. Giovan Battista De Rossi, a 19th century archaeologist, found and published the original version of the Pomponians' notes and dismantled the theory of the big sundial. According to him, the student did not glimpse the floor of the solar meridian beneath the residence of the chaplains but far from the church and in the further south of the religious complex [26].
The part of the solar meridian we know today is due to the efforts of Edmund Buchner, president of the Deutschen Archaeologischen Instituts. In 1976, he suggested the possibility that the obelisk once acted as a gnomon for a vast horologium, despite the evidence offered by De Rossi [27]. Buchner undertook several archaeological investigations that led to positive results in 1979. A section of marble floor, with a sundial line and bronze letters, came to light at a depth of about six metres from the street level, at number 48 in Via di Campo Marzio. The pavement is composed of large square stones, it is 5.4 metres wide and it extends for a maximum length of 7.5 metres. Subsequent excavations, as well as the uncovered section of meridian itself, offer no clues about the possibility that it was part of a larger structure used for measuring the daily hours. The floor is, in fact, crossed throughout its length by a bronze line, placed in a south-north direction, interspersed with 28 rules [28]. As described by Pliny in Naturalis historia, each of these rules concretises a couple of days that the shadow would indicate on the way back and forth, between the summer solstice and the winter solstice.
The . It is assumed that the double distance between the two solstices was not only divided into 360 parts, each one corresponding to 1° ecliptic, but also into 12 sections of 30° or 30 notches each, representing the constellations of the Zodiac [29].
Two other inscriptions complete the section: the first, ΕΤΗΣΙΑΙ ΠΑΥΟΝΤΑΙ (ethesiai pauontai), divides the signs of Leo and Virgo and indicates the period of the year when the ethesies, which are periodic winds that blow towards the end of July, cease; the second inscription, ΘΕΡΟΥΣ ΑΡΧH (therous arke), identifies the point at the beginning of the summer, towards the end of Taurus, in late May. Auber links these inscriptions to the parapegmata, astro-meteorological calendars widespread in ancient Greece, which indicated different periods of the year through easily observable celestial and meteorological phenomena [30]. In one parapegmata in particular, the one made by the Greek astronomer Euctemon between the fifth and fourth century B.C., two sentences exactly coincide with the indications found on the sundial of Augustus: the end of the ethesii winds, with an advance of 12° of ecliptic and therefore 12 rules, and the beginning of the summer, with an advance of 2° of ecliptic.

Materials and Methods
In the present case, the survey of the underground cavity in which a section of the sundial is visible was carried out using a combination of satellite and terrestrial techniques due to the complexity of the site. The measures were repeated two times, in January 2017 and April 2018.
The methodological criteria for the execution of the measures and the processing of the resulting data, which can be assumed as a reference for a survey context as that described in the introduction, are presented in the following section.

GNSS Network
To provide coherent and comparable results with official cartography and astronomical data, the survey has to be referred to an absolute global datum. In Italy, the datum currently officially adopted for the mapping and surveys of public administrations (and consequently the de facto standard) is the ETRF2000 realisation of the ETRS89 datum, referring to the conventional date of January 1, 2008 (epoch 2008.0).
In order to insert the local GNSS network in that datum, it is necessary to connect the survey to several points whose coordinates are known in ETRF2000. In the present case, four permanent stations of the Lazio Region GNSS network were used (Figure 1a), the characteristics of which are summarised in Table 1 [31]. These points can be assumed to have an accuracy less than 1 cm in the ETRF2000 datum [32].  In the Piazza del Parlamento, the closest site to the Augustus sundial which, has a good satellite visibility, a local GNSS network (Figure 1b) was measured, consisting of three stations (1000, 2000 and 3000), in which three GPS+GLONASS dual frequency geodetic receivers (two Topcon GR-5 and one Topcon Legacy GGD with Legant-E antenna) have been placed. In this local network, the receivers were kept on acquisition for over 4 hours in static mode, with a sampling interval of 5 seconds and a 15° cut-off angle.
The vertices of the network (GNSS and total station, see Section 2.2) were materialised on the ground with markers. For each measurement session, a tripod was placed on each vertex; the tripods were not moved until the end of all measures. On each tripod was mounted a Leica/Wild precision tribrach (forced centring device), mounting on it alternatively a GNSS receiver, a total station or a precision target. Therefore, the horizontal and vertical centring accuracy in the network (GNSS and total station) can be estimated at about ±1 mm.
A post-processing calculation was then carried out on the acquired data with Topcon Tools software (version 8.2.3), computing the local network separately and together with the external network including the 4 stations of the Lazio GNSS network (ETRF2000), which were finally assumed as fixed. Figures 1 and 2 show the final configuration of the external segment and of the internal/underground segment of the network. From the post-processing calculation and the subsequent least squares adjustment, the three GNSS vertices in Piazza del Parlamento (Parliament's square) have been determined in the ETRF2000 datum (Figure 2a). For the horizontal coordinates, slightly over 1 millimetre were obtained; for the vertical coordinates, standard ellipses were obtained with semiaxes with a height of about 2 mm. Such uncertainty regions can be derived by the variance-covariance matrix of the adjusted coordinates of the local geodetic network [33].
In addition to the solution described above, two more have been developed: a second solution considering only the closest permanent station (ROUN), and a third considering all four Lazio stations but assuming only ROUN as fixed. The differences between the coordinates obtained from the three different solutions are very small (1-2 mm), so the first solution with the four permanent stations all constrained was taken as final.
The altimetric positions obtained for the vertices 1000, 2000 and 3000 are ellipsoidal heights in the ETRF2000 datum. To convert them to orthometric heights, the ITALGEO2005 national geoid model was utilised, interpolating the grid files provided by the Italian Military Geographical Institute. From this model, the mean geoid height in Piazza del Parlamento was 48.11 metres (Figure 3).

Local Geodetic Network
To connect the GNSS network measured in Piazza del Parlamento to the underground room where the remains of the sundial emerge, a three-dimensional local geodetic network was created (Figure 2a). It is a three-dimensional traverse that includes 10 vertices and has been measured twice, back and forth.
The vertices of this network are located as follows: • 1000, 2000 and 3000 in Piazza del Parlamento (coincident with three GNSS stations).
• 4000 in front of the entrance to the building including the sundial (Via di Campo Marzio). • 5000, 6000, 7000, 8000, 9000 and 10,000 inside the building (10,000 in the sundial room).
The sundial room is located about 7 metres underground with respect to external areas (Via di Campo Marzio, building entrance). To reach it from outside, the only possible path is through some narrow corridors and stairs going down. Thus, the network has a marked three-dimensionality, which is better visible in a perspective representation (Figure 3).
A network with such a particular conformation, characterised by short and sloped legs, involves special operative problems. In particular, a very high accuracy in centring is required, and the collimation targets must be of a type suitable for short distances. For this reason, Leica GMP101 mini prisms and Wild precision collimation targets GFZ11 were used.
The measurements on the 10-vertices network (traverse) were performed twice in both directions (back and forth). All total station measurements (horizontal/vertical angle, distance) were repeated three times in the straight position and three times in the reverse position of the instrument.
The GNSS vertices (1000, 2000 and 3000) were assumed as known and fixed on the coordinates resulting from the computation described at paragraph 2.1 (ETRF2000, UTM33N projection). The traverse adjustment was performed in two distinct modes: a-back and forth measures assumed as independent, processing a continuous traverse (stations were given different names in the back and forth path); b-all measures adjusted together (stations were given the same names in the back and forth path).
The coordinate differences between the two computations described above are summarised in Table 2.  The network adjustment has provided ellipse semiaxes of less than 1 cm: the detail referring to the solution b (assumed as final) is shown in Table 3. The unavoidable propagation of errors proceeding from the GNSS vertices to the underground spaces is visible; however, the uncertainty parameters are contained within 1 cm for the most distant points. From the stations 4000 to 10,000, a suitable number of targets (checkerboard opaque plates) were measured; such points are used for the orientation of the subsequent surveys (TLS and digital photogrammetry).
From station 10,000 in the underground room were directly measured the points of intersection between the longitudinal metal band of the sundial and its transversal notches (Figure 2b). These measures were executed on both sides of the band, in order to be able to compute the sundial axis coordinates by a mean of two points on its sides. After the computation sequence described above, all points (network vertices, targets, sundial axis) were determined in the ETRF2000 datum, UTM33N projection, and orthometric heights.

TLS Survey
A laser scanner survey was performed in the sundial room and in the parts of the building along the three-dimensional (3D) traverse path. The scans were performed with a FARO Focus 3D X130 phase-difference laser scanner [34]. The average density of the TLS measured points in the sundial compartment, due to the short instrument-object distances, is very high (mean distance between contiguous points of about 1.5 millimetre), giving a very detailed representation in the form of a highdensity point cloud.
In total, 12 scans (4 of which were repeated in the 2018 survey) were performed from different positions, proceeding from the outside of the building (near station 4000) down to the underground compartment, including the sundial remains.
The laser scanner data processing was performed with the Leica Cyclone version 9.1 software. The point clouds obtained from the individual scans were joined to one point cloud (Figures 4 and  5). The relative orientation between the point clouds was performed by means of an adequate number of targets (opaque plaques with a checkerboard pattern) and a set of six calibrated spheres of 145 mm in diameter. The target and sphere layout was organised so that between a scan and the adjacent one there were at least four targets and three spheres in common.
The absolute orientation of the assembled point cloud was obtained by means of the threedimensional coordinates of the targets, known from the previous phases of the survey (see Sections 2 and 3) in the ETRF2000 datum. Thus, the TLS survey results georeferenced in the ETRF2000 datum (UTM33N coordinates and orthometric heights).
The scans produced a three-dimensional model of the inner space including the visible part of the sundial (Figure 6). On this 3D model, the points of intersection between the longitudinal metal band of the sundial and the transversal notches have been measured, similarly to what had been done directly with the total station.   An evaluation of the accuracy of the laser scanner survey, and of the congruence between it and the topographic and GNSS measurements, can be obtained from the values of the absolute orientation residuals on the target coordinates, assumed as known. The results are presented in Table 4: the maximum residual on the target used for the orientation (marked with "On" in the table) is about 4 millimetres; for the control target (marked with "Off", not used for the orientation) the maximum value is 13 millimetres.  Table 5 shows the differences between the coordinates of the band axis points obtained from the direct survey (total station) and from the measures on the 3D TLS model. The differences are less than 5 mm, in good agreement with the accuracies of the methods described above.

Photogrammetric Survey
Beside the laser scans, a photogrammetric survey was carried out in the sundial room; this survey was performed by a NIKON D800E full frame calibrated camera equipped with a lens of 20 mm nominal focal length [35]. The photogrammetric project consists of 161 images taken, keeping the focus locked at a constant distance. The calibration procedure was carried out using a circular coded target calibration test-field, holding the same focus length of the survey [36]. All the images were oriented by means of the software Agisoft Metashape version 1.5.0, obtaining a maximum reprojection error of 1.1 pixel and an RMS (Root Mean Square) of 0.16 pixel. The scaling process was handled using some targets, measured by a total station, and by a high precision steel calibrated scale bar of 800 mm (± 0.05 mm). Specifically, four targets have been used as GCPs (Ground Control Points), with the aim to contain the deformations of the photogrammetric model; such deformations were estimated by means of two other targets used as CPs (Check Points) [37]. The maximum errors obtained on GCPs and CPs have been 3 mm and 8 mm respectively, while the error on the calibrated scale bar has been 0.5 mm; after that, the bundle block adjustment procedure was completed.
Once the photogrammetric model was correctly scaled, the image matching process was executed with the aim to obtain a dense cloud point: a different density can be attained from this procedure, according to the purpose of the survey. In this work, low quality was chosen for the building of dense cloud points: a total of 6.5 million points were obtained, corresponding to a density of two points per millimetre, more than sufficient for the survey purpose.
Ground control point and check point coordinates used for the photogrammetric survey were expressed in a local reference system. At a later time, the total station survey has been linked to ETRF 2000-RDN reference system described in section 2 (GNSS network).
The obtained 3D dense point cloud is shown in Figure 7.

Results
The comparison between the two 3D models obtained with laser scanner and photogrammetry was performed by CloudCompare software using seven pairs of homologous points [38]; such points have been selected on the walls and on the floor, where there was no noise in the laser point cloud caused by the presence of small pools of water that wrongly reflected laser signal.
The results of the comparison performed between the two 3D models after the registration phase do not reveal significant differences for the two point clouds; even with other registration options such as ICP with or without scale factor, the RMS on the relative distances is 0.007 metres (Figure 8).
In all cases, the only differences highlighted are mostly due to temporary items present during the surveys:

•
The calibrated scale bar used for photogrammetric survey.

•
The targets placed on the walls during laser scanner survey.

•
The tubes of the portable pump used for drying the floor from groundwater during laser scanner survey.

•
A bag hung on the eastern wall. Both 3D models have generated a digital orthophoto of the floor, where it is possible to determine the coordinates of the mark points relative to the gnomon shadow projected by the sun in the different days of the year. Such coordinates can be compared with those obtained by the total station ( Figure 9).

Discussion
As we explained in the introduction, the expected results of this work were twofold: the correct survey of the sundial site and the comparative evaluation of various techniques in the specific field of application. For this second task, particular attention has been given to the comparison between the most consolidated laser scanning techniques and those based on photogrammetry with the SFM approach more recent and less investigated.
The survey of the sundial allowed to obtain accuracies that will guarantee the reconstruction of the positions of the parts still to be excavated with accuracies from decimetric to submetric, which are sufficient to adequately design the feasibility of the excavation. In particular, it will be very interesting to verify whether some of the buildings in the area have based their foundations on the base of the sundial itself at the time of construction, as some alignments still observable today would suggest. Such overlapping, if verified, would make the excavation of some segments of the sundial very complex technically and probably impossible.
As far as the comparison between the geomatic methodologies used is concerned, it has been verified that the greatest accuracy is provided by a polygonal supported by temporary points surveyed immediately afterwards with double frequency GNSS receivers, because the materialisation of the ground points with topographic markers does not guarantee sufficient accuracy of the framing points. On the other hand, in the historical centre of a city of art, it is not possible to monumentalise control points with the classic forced centring pillars. The comparison between TLS and photogrammetry with the SFM approach shows that they are equivalent, both having pros (speed of execution and processing of clouds for TLS and easy survey and low-cost equipment for photogrammetry) and cons (high cost of instrumentation and logistic complication for TLS and longer processing phase for photogrammetry).

Conclusions and Further Developments
This work has shown how geomatic techniques, from the most consolidated to the most modern ones, can be a support for the planning of excavations in urban areas. The comparison between the various surveys lead to compare them reciprocally and also absolutely, considering that many of the surveys have been fully repeated using different settings.
During the first survey, the GNSS framing network, constrained by the presence of buildings, showed error ellipses on the single points of 1-2 millimetres. The total station survey, based on the same GNSS points, showed an accuracy with maximum deviations of one centimetre. Finally, the comparison of the laser scanning survey (always based on the same GNSS points) with the points measured with the Total Station showed millimetric differences. In the second survey, performed only by photogrammetry, the precision and accuracy were evaluated in 3 and 5 millimetres, respectively.
The comparison between the results of the first and second survey in terms of the point cloud is between 2.5 and 5 mm.
These results allow to estimate errors of a maximum of one centimetre for these types of reliefs, which, on the one hand, show a great achievable accuracy, while on the other hand, could be considered barely sufficient considering the original accuracy with which the sundial was built.
It must also be considered that even small errors on the measurement of the short distance (less than 7 m) currently observable can lead to much larger errors in the reconstruction of position and size of the entire sundial.
In the coming months, the results of this research will allow to define with greater accuracy the positions of any other portion of the sundial (if still existing) and to verify or rebut with greater certainty the various hypotheses on the functioning of the sundial itself. In fact, the position and the height of the obelisk will be obtained considering the distances among such marks together with celestial coordinates of the sun at the epoch of the sundial construction. This will be the task of a future paper still to be published.