Three-Dimensional Modeling and Non-Invasive Diagnosis of a Huge and Complex Heritage Building: The Patriarchal Basilica of Santa Maria Assunta in Aquileia (Udine, Italy)

<jats:p>Three-dimensional modeling and non-invasive diagnosis are fundamental prerequisites for planning reliable assessments of the seismic vulnerability, renovation and conservation of heritage buildings. In the case of multi-layered, huge and complex ancient constructions, various problems can be encountered in the early design phases of interventions, among which there is often a lack of an adequate building documentation. Such issues include drawings that are without a metric scale, not detailed, not updated or not reflecting the real situation. In addition, the fragility of these constructions requires an accurate census of every sign of deterioration in order to prepare an ad hoc intervention for the site. As an example, in this paper, the results of a survey regarding the Patriarchal Basilica of Santa Maria Assunta (Aquileia, Italy) are reported. The basilica has a rich history of about two thousand years. in which each era has marked the actual architectural layout with its own culture and art. The result is an intricate association of complementary and/or interdependent elements that make the building very complex. Given the need to obtain accurate documentation, a realistic representation and a simulation of the criticalities of the structure, which previously did not exist or were not sufficiently accurate, a multi-methodological and multi-scale diagnosis was performed. In detail, the ground-penetrating radar (GPR) technique was applied to verify the presence of structures still buried under some internal surfaces, and a topographic survey, terrestrial laser scanning (TLS), and structure for motion (SfM) aerial and terrestrial photogrammetry were integrated for the detailed survey of the entire internal and external macro-structure. The resulting outcome provided the comprehensive information needed for preparing projects for the preservation, management and restoration of the basilica and the buildings connected to it.</jats:p>

section_content: Three-dimensional modeling and non-invasive diagnosis are the starting point for a highly accurate, detailed and reliable assessment of the status quo of a heritage building. Therefore, the work planning must be carefully prepared by evaluating the purpose of the survey, the degree of detail to be achieved, the correct tools to use for gathering the site data, the timing of data acquisition and processing, the multiple data integration strategy of different formats and the type of post-production output to be processed. In addition, documenting huge and complex buildings, with numerous architectonical details, and narrow or not very accessible interiors, is a challenge considering the correct planning of the survey, the processing of a large amount of dissimilar data and the delivery of cost-effective solutions that are easy to use and interpret. It follows that collaborative work among heterogeneous teams with different types of expertise is required. In recent years, For these reasons, an integrated use of the two methods is recommended during the analysis of complex cultural heritage considering the specific necessities of each diagnostic study. Regarding built cultural heritage, many authors report applications of indirect surveys to identify phenomena of degradation occurring over time such as cracking of walls and columns, presence of humidity, and irregularity of the structure in plan and height [6] [7] [8] [9] [10] [11] [12] [13] [14] [15]. 

Geophysical methods are very useful for assessing the presence of underground structures, addressing conservation and stability issues of architectural monuments, allowing inspection of soil foundations and mechanical properties of structural elements [16] [17] [18] [19], and exploring internal and superficial structures of precious and delicate targets such as statues, wall paintings and mosaics inside historical buildings [20] [21] [22]. Among the nondestructive techniques [23], ground-penetrating radar (GPR) is widely used thanks to the miniaturization of the instrumentation, the high investigation resolution and the minimal impact on the surfaces to be analyzed. In fact, many papers report successful applications of the method for different purposes [24] [25] [26] [27] [28] [29] [30] [31] [32] [33]. 

Following the trend of research in building diagnosis, in this paper the results of a survey of the Patriarchal Basilica of Santa Maria Assunta (Aquileia, Italy) (Figure 1 ) are reported. The work was carried out within the framework of a collaboration agreement between the Superintendence of Archaeology, Arts and Landscape of Friuli Venezia Giulia and the University of Molise, with the aim of planning targeted restoration interventions and assessing the seismic hazard and the vulnerability of the structure (the activities were financed within the framework of the Ministerial Decree of 19 February 2018). The intervention was realized in compliance with the "Guidelines for the evaluation and reduction of the seismic risk of the cultural heritage" issued by the Italian Ministry of Culture (Ministerial Decree of 14 January 2008). They provide for a series of actions summarized in the following phases:

-Cognitive investigation through which the current state of the building is defined with plano-altimetric and structural surveys and a damage and deformation state assessment of the structure. -Historical-critical analysis that guides the designers in the reconstruction of the current state of stress in the light of the changes and events that have affected the building over time. -Mechanical characterization of materials. -Definition of the levels of knowledge and the consequent factors of confidence of the structures. -Structural analysis using numerical models created based on the acquired data. -Critical evaluation of the results obtained. 

The aim of the research was to reach a more in-depth level of knowledge of the surveyed subject in the framework of the first step of the previous workflow. To this end, standard methodological approaches were adopted as follows:

-External topographic survey with total stations aimed at creating a topographic grid with a monograph of fiducial points. The results provided highly detailed documentation, which did not exist until now (a partial survey was carried out through a TLS survey in 2012 [34] ), of the complex architectural heritage, and highlighted several weaknesses of the buildings. 

-Definition of the levels of knowledge and the consequent factors of confidence o structures. -Structural analysis using numerical models created based on the acquired data. -Critical evaluation of the results obtained. 

The aim of the research was to reach a more in-depth level of knowledge o surveyed subject in the framework of the first step of the previous workflow. To this standard methodological approaches were adopted as follows:

-External topographic survey with total stations aimed at creating a topographic with a monograph of fiducial points. The results provided highly detailed documentation, which did not exist until (a partial survey was carried out through a TLS survey in 2012 [34] ), of the com architectural heritage, and highlighted several weaknesses of the buildings.
section_name: Introduction

section_name: The Patriarchal Basilica of Santa Maria Assunta in Aquileia (Udine, Italy)

section_content: The Patriarchal Basilica of Santa Maria Assunta is located in Aquileia (Udine, I The city is located in an alluvial morphological context, the Lower Friulian Plain, in w subsidence processes are in progress, as in the entire arch of the Po valley, which ext from Veneto to Emilia-Romagna. From a lithological point of view, loose and satu
section_name: Geological Background

section_content: The Patriarchal Basilica of Santa Maria Assunta is located in Aquileia (Udine, Italy). The city is located in an alluvial morphological context, the Lower Friulian Plain, in which subsidence processes are in progress, as in the entire arch of the Po valley, which extends from Veneto to Emilia-Romagna. From a lithological point of view, loose and saturated sandy-silty surface horizons are recognized, following a mostly cohesive horizon, which gives place to the underlying thickened sandy beam [35]. The site under analysis falls within the seismic micro-zoning mesh indicated as low hazard [36]. Despite this, Aquileia is surrounded by significant seismogenic structures. In fact, the Carnic-Friuli region represents one of the most seismically active sectors in all of northern Italy [37]. 

It is probable that Aquileia was affected by the propagation of distant epicenters with high magnitudes. Considering the historical sources (Catalog of Strong Earthquakes in Italy in the most recent version [38], Catalog of earthquakes and tsunamis from 11th to 15th century [39] and the Parametric Catalog of Italian Earthquakes [40] ), regarding the seismic effects in Aquileia and its basilica, the information concerns only the earthquakes of 988 and 1348.
section_name: Geological Background

section_content: The city was founded in 181 BC as a Roman colony and it became one of the largest cities of the Roman Empire under the Emperor Diocletian, having a river port and splendid public buildings. Starting from 313 AD, following the Edict of Milan, which decreed an end to religious persecutions, the first cult building in Aquileia dedicated to the Virgin and to the saints Ermacora and Fortunato was built according to the will of Bishop Theodore [41]. It was built in a peripheral area of the city on remains of Roman residences [42] and consisted of two parallel rooms (oriented from west to east, both about 37 × 20 m), decorated with floor mosaics and connected by a transversal room covered by cocciopesto (brick with lime mortar), oriented from south to north (Figure 2a ). The gaps in the decorations still visible today suggest the position of the pillars, of which there must have been six in the main halls and four in the central hall. The entrance to the structure was from the east. Various small service rooms, the changing room and the baptistery were located north of the transverse hall. Today it is possible to admire inside the Patriarchal Basilica of Santa Maria Assunta the splendid floor mosaics from this phase [43, 44] (Figures 3 and 4 ), which are the largest early Christian mosaics in the Western world (760 m 2 ). 

The northern mosaic (Figure 3a ), partially destroyed during the construction of the bell tower, mainly contains floral motifs and representations of animals (Figure 3b-d ). Perhaps the best-known depiction is the fight between the rooster and the tortoise that are, respectively, an allegory of Christ and a symbol of the dark powers of evil (Figure 3a ). 

The southern mosaic is divided into four parts, starting from the actual entrance (Figure 4 ). In the first, various portraits of donors, nodes with crossed ellipses and animals appear. In the second area, there are male and female portraits, depictions of the seasons, enclosed in clip-on medallions, and animals. In the third part, in the central panel, the allegorical scene of winged Victory with crown and palm is represented. Finally, the fourth span, which concludes the cycle of representations, consists of a single mosaic carpet, which represents a sea full of fish, with the story of Jonah, a Jewish prophet, sent by God to preach in Mesopotamia [45]. Among the most interesting images is Jonah praying with his hands raised to save the ship and crew from the storm, or in the jaws of the sea monster, or while being spat on by the monster and while resting under a pergola of vines of pumpkin. All around, sea waves, various fish, octopuses, mollusks and even ducks are represented. 

In the following centuries, the basilica was rebuilt four times, superimposing the new buildings on the remains of the previous ones. In the mid-4th century, the north hall was enlarged (Figure 2b ) and the four-sided portico in front of the facade was created (Figure 2c ). The basilica was divided into three naves by twenty-eight columns and was connected, through the baptistery, to the south hall of the Theodorian basilica. On the sides of the baptistery, two large symmetrical rectangular rooms were decorated with mosaic floors dated between the end of the 4th century and the first decades of the 5th century [46]. Of the northern one, only small fragments remain; the southern one, on the other hand, preserves a large surface, with figures of animals and elaborate geometric motifs. 

Subsequently, at the time of Bishop Chromatius (388-407), the southern hall was also enlarged and new buildings, including the current baptistery, were built (Figure 2d ) [47]. Attila destroyed the northern hall during the siege in 452 and it was not rebuilt. In the first half of the 9th century, the patriarch Maxentius began to renovate the ancient buildings by inserting a semicircular central apse, two side apses and the crypt [48]. 

After the earthquake of 988, the patriarch Poppone started a radical restoration of the complex in Romanesque style [49]. The facade and the altar were rebuilt and the wooden roof was added. The bell tower was built in the 11th century [50]. In the second half of the 12th century, the patriarch Ulrich of Treffen added the frescoes in the crypt [51, 52] (Figure 5 ). In 1348, after an earthquake, the patriarch Marquardo di Randeck replaced the round arches of the central nave with pointed arches, and in the 15th century the wooden ceiling in the shape of a ship was built, which can still be observed today.
section_name: Historical Background

section_name: Survey Methodologies

section_name: Digital Survey

section_content: In order to acquire the coordinates of points in space in a local reference system established in advance, a survey with a total station was carried out using the MTS 805R instrument. The activities were preceded by a careful territorial survey considering the geometry of the area and the location of the elements to be surveyed in the territory. A topographical network was therefore created with a framing of vertices organized in primary benchmarks set in such a way as to minimize the sequence of measurements. This was enabled by considering the largest possible areas using the polygonal scheme, setting up a chain of closed points around the building in which each one is visible from the previous and following points, and bypassing the visual obstacles. To avoid the propagation of errors in the measurements, for each station point, cross-checks were carried out with known points (building edges, road crossings, etc. ) which allowed redistribution and compensation of inaccuracies. Finally, secondary polygonal schemes were then connected to the main polygonal near the elements to be surveyed in detail.
section_name: Topographic Survey

section_content: For several years, the use of TLSs has become a widespread part of surveying activities in architectural and archaeological fields in order to obtain homogeneous and consistent numerical data on which to perform all section and projection operations. Laser scanning enables the automatic and direct measurement of 3D points. Specifically, laser beams are actively emitted and a sensor captures the backscattered signal to obtain a large 3D description of a surface using point clouds [53, 54]. The indiscriminate acquisition of point clouds by the instrument conditions the known critical act and the choice of what one wants to represent. The critical act has shifted from the data acquisition stage to the processing stage through several steps including preprocessing of data, registration of clouds, deletion of superfluous data, transformation of clouds into surfaces, determination of viewing planes and final representation. At the same time, searching for the complete scan of a monument is a difficult quest in view of the complexity of the monument under consideration. TLSs and registration algorithms developed in recent times (Iterative Close Point algorithms, ICPs) make it possible to extend the limits of surveying in a complex area from planimetric and elevation points of view [55]. Thus, it is now possible to
section_name: Terrestrial Laser Scanner (TLS)

section_content: For several years, the use of TLSs has become a widespread part of surveying activities in architectural and archaeological fields in order to obtain homogeneous and consistent numerical data on which to perform all section and projection operations. Laser scanning enables the automatic and direct measurement of 3D points. Specifically, laser beams are actively emitted and a sensor captures the backscattered signal to obtain a large 3D description of a surface using point clouds [53, 54]. The indiscriminate acquisition of point clouds by the instrument conditions the known critical act and the choice of what one wants to represent. The critical act has shifted from the data acquisition stage to the processing stage through several steps including preprocessing of data, registration of clouds, deletion of superfluous data, transformation of clouds into surfaces, determination of viewing planes and final representation. At the same time, searching for the complete scan of a monument is a difficult quest in view of the complexity of the monument under consideration. TLSs and registration algorithms developed in recent times (Iterative Close Point algorithms, ICPs) make it possible to extend the limits of surveying in a complex area from planimetric and elevation points of view [55]. Thus, it is now possible to imagine and carry out surveys of large structures, articulated on different planes and with a complex distribution of internal structures, while maintaining high quality standards in the overall representation of the monument. 

Current instrumental resolutions are capable of visualizing millimeter details even at a distance of several tens of meters. In the case of the Basilica of Aquileia, two phasedifference laser scanners, the FARO Focus 3D MS120 and FARO 3D X330, the latter with a longer range for hard-to-reach areas, were used outside and inside the basilica complex. In total, more than 250 scans were conducted with a resolution of one point every 6 mm at 10 m distance. 

To facilitate the data-recording algorithms, redundant scans were performed from different viewpoints within the basilica complex. The management of field acquisitions is now mainly based on basic topography and redundant overlaying of laser data aimed at maximum coverage of the areas to be surveyed and minimizing the shaded areas still present in the upper and inner parts of the monument [56]. 

The monument complex was segmented into several bodies in which the scanning groups were performed. The boundary areas between the different bodies were used to record the various parts of the complex and ensure a complete working of the entire basilica. 

To obviate the well-known problems of error propagation, which are inevitable on articulated routes of this kind, the integrated use of the total station allowed fixing some known points on the scans of the central nave of the basilica and the adjacent archaeological complex. This system permitted the definition of the relative position of the mother scans in the virtual space on which all others were docked through manual point recognition and the application of ICP and bundle adjustment algorithms.
section_name: Terrestrial Laser Scanner (TLS)

section_content: Part of the work was carried out using SfM-based digital terrestrial and aerial photogrammetry systems. These techniques are based on algorithms that extract the notable points from single photos, deduce the photographic parameters and cross the recognizable points on several photos, detecting the coordinates in the space of the points themselves [57, 58]. Given the complexity of the building, the choice of the two selected methods enabled preliminary photo modeling work of the exterior architectural part of the basilica, taking into account the shaded areas of the upper part also present in the TLS acquisitions. In order to obtain information related to the exterior high part and the roof, acquisitions from a drone were carried out to complete the information of the exterior. The surveying was carried out using a remotely piloted aircraft (Yuneec H520, camera with CMOS sensor (1 ), fixed focal length of 26 mm and resolution of 20 million pixels) with a manual flight control. The images (nadiral and oblique) were acquired through a circular flight around the monument by varying the flight altitude and the distance from the basilica [59]. The acquisitions of the exterior perimeter of the basilica were carried out in order to integrate the information of the elevations and to have enough information for the generation of the orthophotos. The external orientation of the photographic captures was performed by taking advantage of the laser scanner information of the exterior and, in some cases, the targets positioned and recorded by the total station. Despite having different resolutions, the data matching was always organized with an alignment error of less than 0. 025 cm. 

Ground photogrammetry was achieved using a Mirrorless Sony Alpha 5100, equipped with a CMOS sensor (23. 5 × 15. 6) having 24 million pixels and a focal length of 16-50 mm, positioned on a telescopic pole. Frames were acquired along horizontal bands, ensuring an average overlap of 80% both horizontally and vertically. 

In this way, two data sets were obtained consisting of 1. 015 frames for the drone survey and 5. 211 frames for the survey with the digital ground camera. 

Data processing of images was performed with Agisoft Metashape software (LLC. St. Petersburg, Russia) [60] using the structure from motion (SfM) technique [61, 62]. Considering the large size of the data set, the 3D model was processed by dividing it into 8 blocks (church tower, north, east, south, and west parts of the basilica, crypt, mosaics and baptistery). For each block, data processing followed precise steps:

-Aligning images: In this step, the position of the camera is found, the orientation (interior and exterior parameters) of each photo is determined and the model is built with a discrete point cloud. The texture algorithm allows applying photographic images to the model, thereby obtaining a photorealistic 3D model. Finally, it was possible to generate georeferenced orthophotos, Digital Elevation Models (DEMs) relating to the mosaics, the elevations and the vector representation of the sections (using Cloud Compare software (Paris, France)).
section_name: SfM Aerial and Terrestrial Survey

section_content: The data acquired by the different systems are heterogeneous data with completely different information, in terms of both resolution and accuracy. Data integration was enabled by different software tools in the various steps related to data processing. In particular, the most delicate step was related to the roto-translation between the different parts of the basilica within a single reference system. Considering the laser and total station data are more accurate depending on the organization of the fieldwork, some point clouds of the laser were oriented with respect to the total station data. On these, all others were recorded with both manual and automatic point recognition systems. From the laser point clouds, second-order points were recorded to integrate data derived from architectural and aerial photogrammetry of the exterior [63]. Differences in resolution between the two systems resulted in increased registration error between the two parts. Cloud Compare software was used to provide a further estimate of the metric quality of the clouds generated between photogrammetric and laser scanner data. By importing the two clouds in Cloud Compare, it is possible to reduce the deviation and compact the clouds using the "Cloud to Cloud" command. The algorithm allows having a statistical estimate (Gauss function), in which a histogram highlights the deviation intervals between the two clouds. The result is obtained by calculating the mean of the deviations and the standard deviation. 

Among the various algorithms that were used for the representation of architectural structures, special consideration was given to the X-ray visualization of JRC Reconstructor 3D software (Brescia, Italy). X-ray visualization was used to represent discontinuities of apparent contours and edges of structures (Figure 6 ). The use of this visualization system allowed overcoming the problem of mesh processing in architectural representation, taking into account that all apparent contours are visualized by the system. In order to distinguish those to be represented from those that should not be displayed, the correct positioning and distance of the cutting plane are crucial. After creating a section plane, it is necessary to associate a reference user coordinate system (UCS) and a projector that will be used to cut the model. Everything that falls inside the section line should be represented; what is outside, on the other hand, should be eliminated. Care must also be taken in setting up the distance of the projector to avoid representing information that is not relevant to the model. The result of this visualization enables the very accurate highlighting of details that mesh surfaces are unable to detect; the point cloud in most cases highlighted architectural moldings and other details that are often difficult to represent very well. 

The result is obtained by calculating the mean of the deviations and the standard deviation. 

Among the various algorithms that were used for the representation of architectural structures, special consideration was given to the X-ray visualization of JRC Reconstructor 3D software (Brescia, Italy). X-ray visualization was used to represent discontinuities of apparent contours and edges of structures (Figure 6 ). The use of this visualization system allowed overcoming the problem of mesh processing in architectural representation, taking into account that all apparent contours are visualized by the system. In order to distinguish those to be represented from those that should not be displayed, the correct positioning and distance of the cutting plane are crucial. After creating a section plane, it is necessary to associate a reference user coordinate system (UCS) and a projector that will be used to cut the model. Everything that falls inside the section line should be represented; what is outside, on the other hand, should be eliminated. Care must also be taken in setting up the distance of the projector to avoid representing information that is not relevant to the model. The result of this visualization enables the very accurate highlighting of details that mesh surfaces are unable to detect; the point cloud in most cases highlighted architectural moldings and other details that are often difficult to represent very well.
section_name: Integrated Data Management

section_content: GPR is an active electromagnetic method. During data acquisition, the transmitterreceiver antenna is pulled along the surface and signals are sent into the ground. Depending on the features of the material, the radar pulse can be reflected, attenuated, and diffused or totally dissipated. In the case of reflections, the received signal is related to the transmitted signal as a function of the travel velocity [4]. 

The GPR survey involved the floor mosaic of the southern hall of the Theodorian Church inside the modern basilica, the floors of the crypt and of the lateral chapels, and some portion of walls. The instrumentation used was the IDS RIS-K2 georadar equipped with a TRMF (200-600 MHz) multi-frequency monostatic antenna for the ground survey and a 3 GHz antenna for high-resolution surveys of walls. Radar reflections on each line were recorded as 16-bit data in a time window of 90 ns and 40 ns (respectively, for the 200-600 MHz and 3 GHz antennas), acquiring 512 samples per radar scan at 25 scans per mark (unit/marker, 1). During data acquisition, readings were performed continuously on lines arranged in regular grids. On the floors, the spacing between lines was set to 0. 5 m, while on the walls the spacing was 0. 1 m. Standard techniques were used during data
section_name: GPR

section_content: GPR is an active electromagnetic method. During data acquisition, the transmitterreceiver antenna is pulled along the surface and signals are sent into the ground. Depending on the features of the material, the radar pulse can be reflected, attenuated, and diffused or totally dissipated. In the case of reflections, the received signal is related to the transmitted signal as a function of the travel velocity [4]. 

The GPR survey involved the floor mosaic of the southern hall of the Theodorian Church inside the modern basilica, the floors of the crypt and of the lateral chapels, and some portion of walls. The instrumentation used was the IDS RIS-K2 georadar equipped with a TRMF (200-600 MHz) multi-frequency monostatic antenna for the ground survey and a 3 GHz antenna for high-resolution surveys of walls. Radar reflections on each line were recorded as 16-bit data in a time window of 90 ns and 40 ns (respectively, for the 200-600 MHz and 3 GHz antennas), acquiring 512 samples per radar scan at 25 scans per mark (unit/marker, 1). During data acquisition, readings were performed continuously on lines arranged in regular grids. On the floors, the spacing between lines was set to 0. 5 m, while on the walls the spacing was 0. 1 m. Standard techniques were used during data processing using the IdsGred [64] and GPR-SLICE 7. 0 software (Schwerzenbach, Swizzerland) [65], using the following techniques [66] : data editing and conversion, dc-drift removal via bandpass filtering, time-zero correction, background removal in the frequency domain and automatic gaining. Results are presented without converting time in depth as the complexity of layers made it difficult to estimate the velocity of propagation of waves, and thus apply the migration filter. 

Finally, radargram profiles taken over a survey grid representing 3D volume data sets were integrated to produce horizontal time slice maps in order to connect anomalies from closely spaced profiles. Time slice maps aid in the interpretation of archaeological data by effectively showing the size, shape, location and depth of buried targets. In each time slice map based on each individual grid, the relative normalization will display the own maximum and minimum grid values. Depending on the color table assumed, the strongest reflector and weakest reflector will be emphasized. In this work, a gray scale was used to avoid misinterpreting ghost features at the boundaries between two colors.
section_name: GPR

section_content: The 3D scan and the aerial and terrestrial photogrammetric survey allowed the creation of a 3D model of the interior and exterior of the basilica, the baptistery and the bell tower made up of approximately 386,733,987 million points (Figure 7 ). The enormous complexity of the basilica makes it impossible, from both a technical and a principal point of view, to represent an exhaustive set of every single architectural aspect of the monument. In this perspective, only some of the potential graphic representations that can be extracted from the analysis of the three-dimensional model obtained from the survey activities of the monument and its appurtenances are presented. In general, from a preliminary analysis, it was possible to define some constructive aspects, analyze the main components and define the main structural criticalities of the basilica. 

processing using the IdsGred [64] and GPR-SLICE 7. 0 software (Schwerzenbach, Swizzerland) [65], using the following techniques [66] : data editing and conversion, dc-drift removal via bandpass filtering, time-zero correction, background removal in the frequency domain and automatic gaining. Results are presented without converting time in depth as the complexity of layers made it difficult to estimate the velocity of propagation of waves, and thus apply the migration filter. 

Finally, radargram profiles taken over a survey grid representing 3D volume data sets were integrated to produce horizontal time slice maps in order to connect anomalies from closely spaced profiles. Time slice maps aid in the interpretation of archaeological data by effectively showing the size, shape, location and depth of buried targets. In each time slice map based on each individual grid, the relative normalization will display the own maximum and minimum grid values. Depending on the color table assumed, the strongest reflector and weakest reflector will be emphasized. In this work, a gray scale was used to avoid misinterpreting ghost features at the boundaries between two colors.
section_name: Results and Discussion

section_content: The 3D scan and the aerial and terrestrial photogrammetric survey allowed the creation of a 3D model of the interior and exterior of the basilica, the baptistery and the bell tower made up of approximately 386,733,987 million points (Figure 7 ). The enormous complexity of the basilica makes it impossible, from both a technical and a principal point of view, to represent an exhaustive set of every single architectural aspect of the monument. In this perspective, only some of the potential graphic representations that can be extracted from the analysis of the three-dimensional model obtained from the survey activities of the monument and its appurtenances are presented. In general, from a preliminary analysis, it was possible to define some constructive aspects, analyze the main components and define the main structural criticalities of the basilica. Figure 8 shows the horizontal section of the architectural complex extrapolated to the height shown in the lower left box of about one meter from the floor of the central nave. Since it is an X-ray view, it allows viewing everything below the section including the underground excavations of the northern hall of the Theodorian Church and the crypt, as well as the plan of the bell tower, the baptistery and the modern church. It represents a significant advance in the documentation previously available in terms of accuracy, detail and precision, and provides evidence of some considerable irregularities. In the figure, a grid has been superimposed on the plan in order to visually evaluate, for each cell, the presence of corners of the walls that are not entirely regular. 

Figure 8 shows the horizontal section of the architectural complex extrapolated to t height shown in the lower left box of about one meter from the floor of the central na Since it is an X-ray view, it allows viewing everything below the section including t underground excavations of the northern hall of the Theodorian Church and the crypt, well as the plan of the bell tower, the baptistery and the modern church. It represent significant advance in the documentation previously available in terms of accuracy, det and precision, and provides evidence of some considerable irregularities. In the figure grid has been superimposed on the plan in order to visually evaluate, for each cell, t presence of corners of the walls that are not entirely regular. The Basilica of Aquileia today is the result of an intricate interweaving of overlappi elements placed in different time intervals. A complete analysis of each historical co struction phase is not the objective of this work, which instead aims to highlight criti The Basilica of Aquileia today is the result of an intricate interweaving of overlapping elements placed in different time intervals. A complete analysis of each historical construction phase is not the objective of this work, which instead aims to highlight critical elements that could make the current structure vulnerable and endanger the conservation of the site and the safety of the many users who visit it daily. 

The plan of the modern church shows a structure with three naves with an unusual width of the lateral naves certainly due to the need to place walls and columns on old foundations [33]. In fact, the three naves have a greater width on the side of the entrance (left nave 7. 19 m, central nave 10. 80 m, right nave 7. 13 m) which tends to narrow towards the apse by about 0. 10 m (left nave 7. 17 m, central nave 10. 68 m, right nave 7. 09 m). Two rows of 10 columns support 11 arches and separate the 3 naves. A large central portal flanked by two small side doors guarantees the access. Two other openings open at the level of the seventh arch and are slightly asymmetrical. This area does not show particular irregularities in the geometry, while the easternmost part has evident differences. The two small apsidal niches just before the altar are not perfectly aligned and symmetrical with respect to the center of the nave. The body, composed of the altar, two side chapels, the transept, the apse and the sacristy (the room to the right of the altar), is rotated about -3 • towards the right side of the church. Furthermore, in this sector the non-perpendicularity of the walls can be observed, which is, for example, very evident in the Sant'Ambrogio Chapel (to the right of the ship). Externally, the absence of symmetry of the external buttresses is noted. The baptistery and the rooms connected to it are not perfectly aligned with the church but they too are slightly rotated towards the right side of the complex. The bell tower does not appear quadrangular and it is not parallel to the basilica, but it is rotated about 2 •. 

Figure 8 also indicates the location of a selection of sections and side elevations presented below. The bell tower has construction irregularities that make it visibly distorted (Figure 9 ). It is 69. 95 m high, 11. 50 m wide at the base, and tapers off in elevation where, at the level of the upper limit of the bell cell, it reaches a width of 10. 70 m. Although poles reinforced its foundations, in the Middle Ages there were already problems of statics which led to a necessary strengthening of the base with a large supplementary stepped plinth [44]. Today the tower is divided into five blocks, distinguishable by shaped horizontal bands, an octagonal lantern in simple masonry and a conical brick roof. The first basal block probably represents the only element dating back to the first construction phase, which took place in the period of the reign of Poppo, between 940 and 1060 [50]. It is stylistically different from the others; the surfaces are smooth and leveled and are decorated in the corners by pilasters 1. 15 m wide. The entrance, which was at ground level in the first phase of construction, is located on the southern side at a height of 5. 37 m (with respect to the walking surface) and is accessed by an external staircase together with the stepped reinforcement plinth [44, 50]. The deformations highlighted by the survey are of a constructive nature due to the numerous interventions that have been carried out over time. In detail, a slope of the main body from bottom to top in a north-south direction of 2 • appears (see sections AA'-CC' and EE'-GG') and in an east-west direction of 1 • (see sections BB'-DD' and FF'-HH)'. Furthermore, the bell cell, reconstructed in modern times with its arches, the octagonal lantern and the conical brick roof are not centered with respect to the center of gravity of the bell tower. They are shifted in a southwest direction, but straight with respect to an ideal line perpendicular to the horizontal plane, as if to balance the distortion of the shaft. The internal staircase is located in the northwest corner of the tower. It has a spiral structure and is lit by small slit windows. The openings are distributed irregularly, arranged in two crooked rows on elevations BB' and DD' and in a single row on elevations AA' and CC'. The riser of the steps is lower at the base of the bell tower and, vice versa, higher at the top. This generates more twist as the bells are reached. The shaft of the staircase follows the general inclination of the structure. 8 ), it is possible to notice that some walls are not perfectly perpendicula horizontal support plane. In particular, the entrance facade and the external wal right transept show an inclination of about 1° outwards. The column of the left t (vertical side elevation I-I, Figure 10 ) is inclined about 2° towards the apse. By comp Analyzing Figures 10 11 12 13 14 15 (vertical side elevations (I-N) in the positions indicated in Figure 8 ), it is possible to notice that some walls are not perfectly perpendicular to the horizontal support plane. In particular, the entrance facade and the external wall of the right transept show an inclination of about 1 • outwards. The column of the left transept (vertical side elevation I-I, Figure 10 ) is inclined about 2 • towards the apse. By comparison, the right transept column is straight (vertical side elevation K-K', Figure 12 ). The vaulted ceiling of the central nave and the ceiling of the left transept have a considerable depression in the central portion, of approximately 0. 5 m and 0. 15 m, respectively, at the maximum point (Figure 11 ). The most consistent fractures affect the central part of the left nave (Figure 10 ), the internal part of the main facade (Figure 15 ) and the walls of the transepts (Figure 14 ). The fractures develop vertically; no horizontal fractures are documented. The survey of these critical issues represents a current state of the art and can be used for monitoring of millimeter deviations over time. As anticipated, because of targeted inspections, according to the needs of future activities, similar representations can be obtained in the most suitable formats. All this information is contained in the acquired data set. 

Remote Sens. 2023, 15, x FOR PEER REVIEW 15 of 25 the right transept column is straight (vertical side elevation K-K,' Figure 12 ). The vaulted ceiling of the central nave and the ceiling of the left transept have a considerable depression in the central portion, of approximately 0. 5 m and 0. 15 m, respectively, at the maximum point (Figure 11 ). The most consistent fractures affect the central part of the left nave (Figure 10 ), the internal part of the main facade (Figure 15 ) and the walls of the transepts (Figure 14 ). The fractures develop vertically; no horizontal fractures are documented. The survey of these critical issues represents a current state of the art and can be used for monitoring of millimeter deviations over time. As anticipated, because of targeted inspections, according to the needs of future activities, similar representations can be obtained in the most suitable formats. All this information is contained in the acquired data set. the right transept column is straight (vertical side elevation K-K,' Figure 12 ). The vaulted ceiling of the central nave and the ceiling of the left transept have a considerable depression in the central portion, of approximately 0. 5 m and 0. 15 m, respectively, at the maximum point (Figure 11 ). The most consistent fractures affect the central part of the left nave (Figure 10 ), the internal part of the main facade (Figure 15 ) and the walls of the transepts (Figure 14 ). The fractures develop vertically; no horizontal fractures are documented. The survey of these critical issues represents a current state of the art and can be used for monitoring of millimeter deviations over time. As anticipated, because of targeted inspections, according to the needs of future activities, similar representations can be obtained in the most suitable formats. All this information is contained in the acquired data set. the right transept column is straight (vertical side elevation K-K,' Figure 12 ). The vaulted ceiling of the central nave and the ceiling of the left transept have a considerable depression in the central portion, of approximately 0. 5 m and 0. 15 m, respectively, at the maximum point (Figure 11 ). The most consistent fractures affect the central part of the left nave (Figure 10 ), the internal part of the main facade (Figure 15 ) and the walls of the transepts (Figure 14 ). The fractures develop vertically; no horizontal fractures are documented. The survey of these critical issues represents a current state of the art and can be used for monitoring of millimeter deviations over time. As anticipated, because of targeted inspections, according to the needs of future activities, similar representations can be obtained in the most suitable formats. All this information is contained in the acquired data set. Figures 16 and 17 show some drawings related to the crypt. In particular, a reflectance visualization (Figure 16b ) and a mirrored visualization (Figure 16c ) of the 3D model are reproduced, a survey detail that has not been previously realized. Figure 17c represents an orthophoto of the vaulted ceiling. The crypt, probably the oldest part of the basilica, is not perfectly centered with respect to the body of the basilica. The plan highlights the construction of a semicircular apse, enlarged towards the west until it assumes a U shape, included in an external quadrangular structure. The access to the room is via two openings on the sides of the presbytery on the western side where there is a small niche in the center. The vaulted ceiling, enriched by frescoes dating back to the 12th century illustrating the legend of the evangelist Marco, is supported by six columns. 

Remote Sens. 2023, 15, x FOR PEER REVIEW Figures 16 and 17 show some drawings related to the crypt. In particular, a reflectance visualization (Figure 16b ) and a mirrored visualization (Figure 16c ) of the 3D model are reproduced, a survey detail that has not been previously realized. Figure 17c represents an orthophoto of the vaulted ceiling. The crypt, probably the oldest part of the basilica, is not perfectly centered with respect to the body of the basilica. The plan highlights the construction of a semicircular apse, enlarged towards the west until it assumes a U shape, included in an external quadrangular structure. The access to the room is via two openings on the sides of the presbytery on the western side where there is a small niche in the center. The vaulted ceiling, enriched by frescoes dating back to the 12th century illustrating the legend of the evangelist Marco, is supported by six columns. some drawings related to the crypt. In particular, a reflectance visualization (Figure 16b ) and a mirrored visualization (Figure 16c ) of the 3D model are reproduced, a survey detail that has not been previously realized. Figure 17c represents an orthophoto of the vaulted ceiling. The crypt, probably the oldest part of the basilica, is not perfectly centered with respect to the body of the basilica. The plan highlights the construction of a semicircular apse, enlarged towards the west until it assumes a U shape, included in an external quadrangular structure. The access to the room is via two openings on the sides of the presbytery on the western side where there is a small niche in the center. The vaulted ceiling, enriched by frescoes dating back to the 12th century illustrating the legend of the evangelist Marco, is supported by six columns. The DEM of the mosaic surfaces of the crypt of the excavations, the main basilica and the Süd-Halle shows some irregularities, which suggest the presence of underground and very shallow archaeological structures (Figure 18 ). The DEM was generated based on the dense cloud. This option was preferred to the mesh model, as it provides results that are more accurate and allows for faster processing, since the mesh generation step can be skipped. A pixel size of 2. 25 mm/pix was chosen according to the average ground sampling resolution of the original images. 

In detail, in the crypt of the excavations (A in Figure 18 ), a variation in elevation of 71. 8 cm was registered. While the mosaic around the tower is placed at a lower level of the whole scale in the range between -61. 8 cm and -31. 5 cm, the western one has higher elevations (in the range between -31. 5 cm and 10 cm). In addition, we note a rectilinear convex anomaly (1 in Figure 18 ), parallel to the portion of the wall in the light on the right of the image (dashed black in Figure 18 ) and a concave rectangular anomaly (2 in Figure 18 ). The latter has dimensions compatible with a burial. The presence of buried structures probably prevented the manufacturers of the mosaics from laying a homogeneous mosaic on the ground with respect to the altitude, or conditioned the subsidence of the deposits over the years, opposing resistance to the lowering of the floor compared to areas without archaeological remains. 

In the Süd-Halle (B in Figure 18 ), the DEM shows overall variations of about 48 cm. It clearly highlights already known discrepancies due to the presence of eleven tombs that cut through the mosaic which materialize with depressions in the model and the absence of mosaic tiles (medium elevations). The highest altitudes are recorded in the points where the mosaic is preserved. 

The analysis of the DEM of the floor mosaic of the basilica (C in Figure 18 ) also shows the distribution of irregularities attributable to older structures present below the current surface. In detail, the surface of the mosaic is not perfectly flat and a difference in height of about 0. 67 m was found between the highest points (towards the entrance to the basilica, n. 7 in Figure 18 ) and the lowest points (near the altar). In the DEM, four buried walls clearly emerge (numbers 4-6, indicated with black arrows in Figure 18 ). These deformations are arranged in a regular manner forming an orthogonal grid in synchrony with the most ancient archaeological structures found on the site, allowing us to imagine the delineation of very large chambers. At the point indicated with the n. 7, even if the dividing walls are not clearly highlighted, the area with the greatest elevation has regular outlines. 

The GPR results confirm the presence of buried structures, even if they are not perfectly comparable to the surface DEM as they relate to different subsoil portions (Figure 19 ). The GPR map shows an articulated distribution of amplitude maxima of the electromagnetic signal, which seem to define regular rooms, with trends perfectly in line with the current building. In the crypt of the frescoes, on the altar in the right portion, in the two chapels and in the easternmost potion of the central and southern nave, the maxima of amplitude are clear and not very thick. Conversely, at the entrance to the basilica in the investigated western part, the anomalous zones describe highly disturbed regions. From the comparison with the DEM (Figure 20 ), these areas fall within the zone indicated with n. 7. The four elements mentioned in DEM shallow analysis (nn. [3] [4] [5] [6] are partially connected to the GPR anomalies in depth, supporting the hypothesis of the presence of walls. In addition, the wall section n. 3 seems to join up with the point where an abrupt variation in the wall stratigraphy is noted on the wall of the left aisle and, at the specular point on the right aisle, the GPR results show a probable similar variation in the wall texture (Figure 21 ).
section_name: Results and Discussion

section_content: In this paper, the issues related to the status quo of a huge and complex heritage building, such as the Patriarchal Basilica of Santa Maria Assunta in Aquileia, were investigated in an organic and exhaustive way, validating standard methodological workflows. Surveying was not an easy task. The data acquisition, integrated processing and interpretation of the results required a large amount of work and high precision at every stage. The final product was the acquisition of detailed documentation, which was previously non-existent or approximate, and is useful as a starting point for planning reliable assessments of the seismic vulnerability, renovation and conservation of heritage buildings. 

The integration of TLS and aerial and terrestrial photogrammetry allowed a complete geometric survey of the construction, defining with high precision the plano-altimetric characteristics of all the masonry elements. In the absence of previous documentation, it was possible to identify the position and dimensions of beams, pillars, stairs and partitions, and roof slabs, determining their type and vertical section. Some construction anomalies were detected in the structure, as highlighted in the bell tower, and some weaknesses were identified in the ceilings of the naves and in the walls of the entrance and transept. The crack pattern was also evaluated by classifying the lesions from a geometric point of view (extension and width). In addition, detailed representations, which were not previously elaborated, were produced for some environments. The integrated analysis of GPR and DEM results led to hypothesizing about the presence of buried structures belonging to the oldest construction phases under the investigated floors. 

For the sake of brevity, only a part of the possible graphical representations was presented here. Although these observations may be regrettable considering the age of the structure and its many construction phases superimposed over time, the data obtained are new and useful for continuing the analysis and restoration processes. Following the Italian legislation on seismic vulnerability, this work will be useful for designers and engineers to identify the constituent elements of the resistant body, placing particular attention on the construction techniques, the construction details and the connections between the elements. Furthermore, starting from a visual survey to the programming of further destructive and non-destructive tests, it will represent the basis for the exact identification of the materials by evaluating both the level of degradation and the mechanical properties. The well-known construction path of the artefact and the changes it has undergone over time, together with the results obtained, will help to identify the possible areas of discontinuity of the material and guide the structural analysis using numerical models.
section_name: Conclusions

section_name: Acknowledgments:

section_content: We are very thankful to Simonetta Bonomi, Head of the Superintendence of Archeology, Arts and Landscape of Friuli Venezia Giulia for involving us in this important project.
section_name: Acknowledgments:

section_content: Funding: This research was funded by SUPERINTENDENCE OF ARCHEOLOGY, ARTS AND LANDSCAPE OF FRIULI VENEZIA GIULIA (Funding number: CUP F37E18000180001, CIG Z3C3151094 ) and the APC was funded by UNIVERSITY OF MOLISE (Paolo Mauriello, Project Responsible for University of Molise ).

section_content: Data Availability Statement: Data can be requested from the authors.

section_name: Conflicts of Interest:

section_content: Author Contributions: Conceptualization, A. A., M. C., R. G., V. G. and P. M. ; methodology, A. A., M. C., R. G., V. G. and P. M. ; investigation, A. A., M. C., R. G., V. G. and P. M. ; data curation, A. A., M. C., R. G., V. G. and P. M. ; writing-original draft preparation, A. A., M. C., R. G., V. G. and P. M. ; writing-review and editing, A. A., M. C., R. G., V. G. and P. M. ; project administration, P. M. ; funding acquisition, P. M. All authors have read and agreed to the published version of the manuscript.

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: