Big Earth Data for Cultural Heritage in the Copernicus Era

Digital data is stepping in its golden age characterized by an increasing growth of both classical and emerging big earth data along with trans- and multidisciplinary methodological approaches and services addressed to the study, preservation and sustainable exploitation of cultural heritage (CH). The availability of new digital technologies has opened new possibilities, unthinkable only a few years ago for cultural heritage. The currently available digital data, tools and services with particular reference to Copernicus initiatives make possible to characterize and understand the state of conservation of CH for preventive restoration and opened up a frontier of possibilities for the discovery of archaeological sites from above and also for supporting their excavation, monitoring and preservation. The different areas of intervention require the availability and integration of rigorous information from different sources for improving knowledge and interpretation, risk assessment and management in order to make more successful all the actions oriented to the preservation of cultural properties. One of the biggest challenges is to fully involve the citizen also from an emotional point of view connecting “pixels with people” and “bridging” remote sensing and social sensing

SAR Sentinel 1 imaging and detection of palaeo-landscape features in the mediterranean area

The use of satellite radar in landscape archaeology offers great potential for manifold applications, such as the detection of ancient landscape features and anthropogenic transformations. Compared to optical data, the use and interpretation of radar imaging for archaeological investigations is more complex, due to many reasons including that: (i) ancient landscape features and anthropogenic transformations provide subtle signals, which are (ii) often covered by noise; and, (iii) only detectable in specific soil characteristics, moisture content, vegetation phenomenology, and meteorological parameters. In this paper, we assessed the capability of SAR Sentinel 1 in the imaging and detection of palaeo-landscape features in the Mediterranean area of Tavoliere delle Puglie. For the purpose of our investigations, a significant test site (larger than 200 km2) was selected in the Foggia Province (South of Italy) as this area has been characterized for millennia by human frequentation starting from (at least) the Neolithic. The results from the Sentinel 1 (S-1) data were successfully compared with independent data sets, and the comparison clearly showed an excellent match between the S-1 based outputs and ancient anthropogenic transformations and landscape features

Landslide monitoring using multi-temporal SAR interferometry with advanced persistent scatterers identification methods and super high-spatial resolution terraSAR-X Images

Landslides are one of the most common and dangerous threats in the world that generate considerable damage and economic losses. An efficient landslide monitoring tool is the Differential Synthetic Aperture Radar Interferometry (DInSAR) or Persistent Scatter Interferometry (PSI). However, landslides are usually located in mountainous areas and the area of interest can be partially or even heavily vegetated. The inherent temporal decorrelation that dramatically reduces the number of Persistent Scatters (PSs) of the scene limits in practice the application of this technique. Thus, it is crucial to be able to detect as much PSs as possible that can be usually embedded in decorrelated areas. High resolution imagery combined with efficient pixel selection methods can make possible the application of DInSAR techniques in landslide monitoring. In this paper, different strategies to identify PS Candidates (PSCs) have been employed together with 32 super high-spatial resolution (SHR) TerraSAR-X (TSX) images, staring-spotlight mode, to monitor the Canillo landslide (Andorra). The results show that advanced PSI strategies (i.e., the temporal sub-look coherence (TSC) and temporal phase coherence (TPC) methods) are able to obtain much more valid PSs than the classical amplitude dispersion (DA) method. In addition, the TPC method presents the best performance among all three full-resolution strategies employed. The SHR TSX data allows for obtaining much higher densities of PSs compared with a lower-spatial resolution SAR data set (Sentinel-1A in this study). Thanks to the huge amount of valid PSs obtained by the TPC method with SHR TSX images, the complexity of the structure of the Canillo landslide has been highlighted and three different slide units have been identified. The results of this study indicate that the TPC approach together with SHR SAR images can be a powerful tool to characterize displacement rates and extension of complex landslides in challenging areasPeer ReviewedPostprint (published version

Interpretation of multispectral satellite data as a tool for detecting archaeological artifacts (Navkur Plain and Karamleis Plain, Iraq)

Contemporary studies of geographical space, including archaeological research, incorporate multiple spatial digital data. Such data provide an opportunity to extend research to large areas, and to objectify studies on the basis of quantitative data thus obtained and gaining access to the hard-to-reach study area. Examples of such data are satellite images at various spatial resolutions and in a wide spectrum of electromagnetic radiation (visible, infrared, and microwave). The authors made an attempt to use satellite images to analyze the areas of probable location of the Battle of Gaugamela (the Navkur Plain and the Karamleis Plain in Iraq). The photointerpretation was performed, enhanced by the multivariate processing of the multispectral image. The aim of the work was indicating the most likely places where the camp and the battle were located based on the visual interpretation of an array of satellite data. The adopted methodology of precise allocation of interpretative values to remote sensing materials for every detected artifact provided an opportunity to accumulate an extensive amount of information. It also provided the basis for a synthetic analysis regarding the methods of image processing on the one hand and the dates of recording on the other. It turned out that the season in which the photos are recorded is very important—although the best data for analysis turned out to be the autumn data (38% of all recognized artifacts), the use of data from three seasons increased the total number of indicated artifacts by as much as about 50% (the so-called unique detections). In addition, advanced image processing (such as principal component analysis and decorrelation stretch) turned out to be important, as it increased the number of areal artifacts by 31% compared to the interpretation of only photos in natural (true) color composite and false color composite (with near-infrared). The conducted analyses have confirmed the usefulness of high-resolution satellite data for archaeological applications, and the detected and described anomalies visible in satellite images are excellent material for selecting sites for detailed field research

Advanced pixel selection and optimization algorithms for Persistent Scatterer Interferometry (PSI)

Tesi amb diferents seccions retallades per dret de l'editorPremi Extraordinari de Doctorat, promoció 2018-2019. Àmbit de les TICGround deformation measurements can provide valuable information for minimization of associated loss and damage caused by natural and environmental hazards. As a kind of remote sensing technique, Persistent Scatterer Interferometry (PSI) SAR is able to measure ground deformation with high spatial resolution, efficiently. Moreover, the ground deformation monitoring accuracy of PSI techniques can reach up to millimeter level. However, low coherence could hinderthe exploitation of SAR data, and high-accuracy deformation monitoring can only be achieved by PSI for high quality pixels. Therefore, pixel optimization and identification of coherent pixels are crucial for PSI techniques. In this thesis, advanced pixel selection and optimization algorithms have been investigated. Firstly, a full-resolution pixel selection method based on the Temporal Phase Coherence (TPC) has been proposed. This method first estimates noise phase term of each pixel at interferogram level. Then, for each pixel, its noise phase terms of all interferograms are used to assess this pixel’s temporal phase quality (i.e., TPC). In the next, based on the relationship between TPC and phase Standard Deviation (STD), a threshold can be posed on TPC to identify high phase quality pixels. This pixel selection method can work with both Deterministic Scatterers (PSs) and Distributed Scatterers (DSs). To valid the effectiveness of the developed method, it has been used to monitor the Canillo (Andorra) landslide. The results show that the TPC method can obtained highest density of valid pixels among the employed three approaches in this challenging area with X-band SAR data. Second, to balance the polarimetric DInSAR phase optimization effect and the computation cost, a new PolPSI algorithm is developed. This proposed PolPSI algorithm is based on the Coherency Matrix Decomposition result to determine the optimal scattering mechanism of each pixel, thus it is named as CMD-PolPSI. CMDPolPSI need not to search for solution within the full space of solution, it is therefore much computationally faster than the classical Equal Scattering Mechanism (ESM) method, but with lower optimization performance. On the other hand, its optimization performance outperforms the less computational costly BEST method. Third, an adaptive algorithm SMF-POLOPT has been proposed to adaptive filtering and optimizing PolSAR pixels for PolPSI applications. This proposed algorithm is based on PolSAR classification results to firstly identify Polarimetric Homogeneous Pixels (PHPs) for each pixel, and at the same time classify PS and DS pixels. After that, DS pixels are filtered by their associated PHPs, and then optimized based on the coherence stability phase quality metric; PS pixels are unfiltered and directly optimized based on the DA phase quality metric. SMF-POLOPT can simultaneously reduce speckle noise and retain structures’ details. Meanwhile, SMF-POLOPT is able to obtain much higher density of valid pixels for deformation monitoring than the ESM method. To conclude, one pixel selection method has been developed and tested, two PolPSI algorithms have been proposed in this thesis. This work make contributions to the research of “Advanced Pixel Selection and Optimization Algorithms for Persistent Scatterer InterferometryLes mesures de deformació del sòl poden proporcionar informació valuosa per minimitzar les pèrdues i els danys associats causats pels riscos naturals i ambientals. Com a tècnica de teledetecció, la interferometria de dispersors persistents (Persistent Scatter Interferometry, PSI) SAR és capaç de mesurar de forma eficient la deformació del terreny amb una alta resolució espacial. A més, la precisió de monitorització de la deformació del sòl de les tècniques PSI pot arribar a arribar a nivells del mil·límetre. No obstant això, una baixa coherència pot dificultar l’explotació de dades SAR i el control de deformació d’alta precisió només es pot aconseguir mitjançant PSI per a píxels d’alta qualitat. Per tant, l’optimització de píxels i la identificació de píxels coherents són crucials en les tècniques PSI. En aquesta tesi s¿han investigat algorismes avançats de selecció i optimització de píxels. En primer lloc, s'ha proposat un mètode de selecció de píxels de resolució completa basat en la coherència temporal de fase (Temporal Phase Coherence, TPC). Aquest mètode estima per primera vegada el terme de fase de soroll de cada píxel a nivell d’interferograma. A continuació, per a cada píxel, s'utilitzen els termes de la fase de soroll de tots els interferogrames per avaluar la qualitat de fase temporal d'aquest píxel (és a dir, TPC). A la següent, basant-se en la relació entre el TPC i la desviació estàndard de fase (STD), es pot plantejar un llindar de TPC per identificar píxels de qualitat de fase alta. Aquest mètode de selecció de píxels es capaç de detectar tant els dispersors deterministes (PS) com els distribuïts (DS). Per validar l’eficàcia del mètode desenvolupat, s’ha utilitzat per controlar l’esllavissada de Canillo (Andorra). Els resultats mostren que el mètode TPC pot obtenir la major densitat de píxels vàlids, comparat amb els mètodes clàssics de selecció, en aquesta àrea difícil amb dades de SAR de banda X. En segon lloc, per equilibrar l’efecte d’optimització de fase DInSAR polarimètrica i el cost de càlcul, es desenvolupa un nou algorisme de PolPSI. Aquest algorisme proposat de PolPSI es basa en el resultat de la descomposició de la matriu de coherència per determinar el mecanisme de dispersió òptim de cada píxel, de manera que es denomina CMD-PolPSI. CMDPolPSI no necessita buscar solucions dins de l’espai complet de la solució, per tant, és molt més eficient computacionalment que el mètode clàssic de mecanismes d’igualtat de dispersió (Equal Scattering Mechanism, ESM), però amb un efecte d’optimització no tant òptim. D'altra banda, el seu efecte d'optimització supera el mètode BEST, el que te un menor cost computacional. En tercer lloc, s'ha proposat un algoritme adaptatiu SMF-POLOPT per al filtratge adaptatiu i l'optimització de píxels PolSAR per a aplicacions PolPSI. Aquest algorisme proposat es basa en els resultats de classificació PolSAR per identificar primer els píxels homogenis polarimètrics (PHP) per a cada píxel i, alhora, classificar els píxels PS i DS. Després d'això, els píxels DS es filtren pels seus PHP associats i, a continuació, s'optimitzen en funció de la mètrica de qualitat de la fase d'estabilitat de coherència; els píxels classificats com PS no es filtren i s'optimitzen directament en funció de la mètrica de qualitat de la fase DA. SMF-POLOPT pot reduir simultàniament el soroll de la fase interferomètrica i conservar els detalls de les estructures. Mentrestant, SMF-POLOPT aconsegueix obtenir una densitat molt més alta de píxels vàlids per al seguiment de la deformació que el mètode ESM. Per concloure, en aquesta tesi s’ha desenvolupat i provat un mètode de selecció de píxels, i s’han proposat dos algoritmes PolPSI. Aquest treball contribueix a la recerca en "Advanced Pixel Selection and Optimization Algorithms for Persistent Scatterer Interferometry"Postprint (published version

Potentialities and limitations of research on VHRS data : Alexander the Great’s military camp at Gaugamela on the Navkur Plain in Kurdish Iraq as a test case

This paper presents a selected aspect of research conducted within the Gaugamela Project, which seeks to finally identify the location of one of the most important ancient battles: the Battle of Gaugamela (331 BCE). The aim of this study was to discover material remains of the Macedonian military camp on the Navkur Plain in Kurdish Iraq. For this purpose, three very high resolution satellite (VHRS) datasets from Pleiades and WorldView-2 were acquired and subjected to multi-variant image processing (development of different color composites, integration of multispectral and panchromatic images, use of principle component analysis transformation, use of vegetation indices). Documentation of photointerpretation was carried out through the vectorization of features/areas. Due to the character of the sought-after artifacts (remnants of a large enclosure), features were categorized into two types: linear features and areal features. As a result, 19 linear features and 2 areal features were found in the study area of the Mahad hills. However, only a few features fulfilled the expected geometric criteria (layout and size) and were subjected to field groundtruthing, which ended in negative results. It is concluded that no traces have been found that could be interpreted as remnants of an earthen enclosure capable of accommodating around 47,000 soldiers. Further research perspectives are also suggested

Integrating Remote Sensing and Geophysics for Exploring Early Nomadic Funerary Architecture in the “Siberian Valley of the Kings”

This article analyses the architecture of the Early Iron Age royal burial mound Tunnug 1 in the “Siberian Valley of the Kings” in Tuva Republic, Russia. This large monument is paramount for the archaeological exploration of the early Scythian period in the Eurasian steppes, but environmental parameters make research on site difficult and require the application of a diversity of methods. We thus integrate WorldView-2 and ALOS-2 remote sensing data, geoelectric resistivity and geomagnetic survey results, photogrammetry-based DEMs, and ortho-photographs, as well as excavation in order to explore different aspects of the funerary architecture of this early nomadic monument. We find that the large royal tomb comprises of a complex internal structure of radial features and chambers, and a rich periphery of funerary and ritual structures. Geomagnetometry proved to be the most effective approach for a detailed evaluation of the funerary architecture in our case. The parallel application of several surveying methods is advisable since dataset comparison is indispensable for providing context

Evaluation of Multi-frequency Synthetic Aperture Radar for Subsurface Archaeological Prospection in Arid Environments

The discovery of the subsurface paleochannels in the Saharan Desert with the 1981 Shuttle Imaging Radar (SIR-A) sensor was hugely significant in the field of synthetic aperture radar (SAR) remote sensing. Although previous studies had indicated the ability of microwaves to penetrate the earth’s surface in arid environments, this was the first applicable instance of subsurface imaging using a spaceborne sensor. And the discovery of the ‘radar rivers’ with associated archaeological evidence in this inhospitable environment proved the existence of an earlier less arid paleoclimate that supported past populations. Since the 1980’s SAR subsurface prospection in arid environments has progressed, albeit primarily in the fields of hydrology and geology, with archaeology being investigated to a lesser extent. Currently there is a lack of standardised methods for data acquisition and processing regarding subsurface imaging, difficulties in image interpretation and insufficient supporting quantitative verification. These barriers keep SAR technology from becoming as integral as other remote sensing techniques in archaeological practice The main objective of this thesis is to undertake a multi-frequency SAR analysis across different site types in arid landscapes to evaluate and enhance techniques for analysing SAR within the context of archaeological subsurface prospection. The analysis and associated fieldwork aim to address the gap in the literature regarding field verification of SAR image interpretation and contribute to the understanding of SAR microwave penetration in arid environments. The results presented in this thesis demonstrate successful subsurface imaging of subtle feature(s) at the site of ‘Uqdat al-Bakrah, Oman with X-band data. Because shorter wavelengths are often ignored due to their limited penetration depths as compared to the C-band or L-band data, the effectiveness of X-band sensors in archaeological prospection at this site is significant. In addition, the associated ground penetrating radar and excavation fieldwork undertaken at ‘Uqdat al-Bakrah confirm the image interpretation and support the quantitative information regarding microwave penetration

Remote Sensing and Geosciences for Archaeology

This book collects more than 20 papers, written by renowned experts and scientists from across the globe, that showcase the state-of-the-art and forefront research in archaeological remote sensing and the use of geoscientific techniques to investigate archaeological records and cultural heritage. Very high resolution satellite images from optical and radar space-borne sensors, airborne multi-spectral images, ground penetrating radar, terrestrial laser scanning, 3D modelling, Geographyc Information Systems (GIS) are among the techniques used in the archaeological studies published in this book. The reader can learn how to use these instruments and sensors, also in combination, to investigate cultural landscapes, discover new sites, reconstruct paleo-landscapes, augment the knowledge of monuments, and assess the condition of heritage at risk. Case studies scattered across Europe, Asia and America are presented: from the World UNESCO World Heritage Site of Lines and Geoglyphs of Nasca and Palpa to heritage under threat in the Middle East and North Africa, from coastal heritage in the intertidal flats of the German North Sea to Early and Neolithic settlements in Thessaly. Beginners will learn robust research methodologies and take inspiration; mature scholars will for sure derive inputs for new research and applications

ANCIENT MINING LANDSCAPES AND HABITATIVE SCENERIES IN THE URBAN AREA OF CENTOCELLE: GEOMATIC APPLICATIONS FOR THEIR IDENTIFICATION, MEASUREMENT, DOCUMENTATION AND MONITORING

This study, focused on the Archaeological Park of Centocelle, was carried out to test the functionality of different geomatic products for the identification and monitoring of complex archaeological evidences in urban areas. The method proposes a better knowledge of the archaeological context as a tool to favour a better protection, allowing the establishment of limits to urban enlargement in areas of respect. The test area is chosen because of the combined presence of hypogeal evidences related to Roman and pre-Roman exploitation of local litotypes and for the dense presence of archaeological vestiges at its surface, related to the inhabitation function of the zone in a period contemporaneous to the beginning of the quarrying activities. The methods used are the digital photogrammetry, 3D modelling, remote sensing interpretation and digital cartography. The protocol is then customized for the peculiarities of the area under study, considering both the underground structures and the ones at the surface. Archaeological features are identified by processing optical and SAR dataset to enhance the contrast of archaeological features from the background. Historical and recent DSM have been then compared to evaluate the evolutions of local topography. Concerning the study of the subterranean quarrying system in the area, a 3D model of one gallery was produced, with the aim to understand the type of ancient exploitation. A DTM of the toolmarks was then produced to understand the technological skills used for the exploitation of the local tuff and used as an indirect proof for chronological interpretation. A final trial of PSInSAR was addressed to test the method for monitoring the hypogeal levels. Several field prospections were executed, in order to first set the method properly and then validate the results