3. Results
The implementation of the set goal is demonstrated through the study of the historical and cultural heritage object, the Citadel defensive complex located in Lviv, where the Stalag-328 prisoner-of-war concentration camp was situated during World War II.
Using the obtained interferograms generated from 36 satellite images acquired from the Sentinel-1 satellite over three years from 2020 to 2022, we were able to identify the extremes of vertical ground displacements in the research area. By discarding the indicators of minimal displacements, only the maximum displacements were considered. There were 24 points identified with such displacements. The range of vertical ground displacements at these points varies from 17 to 36 centimeters. Additional investigations led to the exclusion of 5 points located within the object's territory in the built-up area. It was hypothesized that the vertical displacements in the built-up area may have occurred due to specific repair works.
Furthermore, based on previous research [
35], the territory of prisoner shootings and fraternal burials was overlaid onto the modern topographic plan of 2007 at a scale of 1:5000. Therefore, the points of extreme vertical displacements were color-coded: green for those falling within the areas of shootings and burials, yellow for those falling on the paved infrastructure likely formed by repair works, and red for those falling on the undeveloped territory. Additionally, the current boundary of the historical and cultural heritage object, represented by a thick black line, is delineated on the plan (
Figure 6).
The analysis of vertical displacements can reveal potential threats to historical structures, such as precipitation, erosion, soil subsidence, or geodynamic processes. This allows appropriate measures to be taken to prevent potential damage or destruction. Information about vertical displacements can be used for planning restoration work and preserving historical structures. This enables effective restoration measures that take into account the potential impacts of geodynamic processes [
36,
37]. Maps of vertical displacements serve as documentation of the state of historical structures at a certain period in time. This is important for archiving and monitoring changes over time, as well as for further scientific research. Smoothed maps of maximum vertical displacements of the Earth's surface in the territory of the Citadel for the years 2020-2022 are presented in
Figure 7.
After determining the points of maximum surface displacements using radar interferometry, we proceeded with ground-based non-invasive investigations using ground-penetrating radar (GPR). Seven points were processed, that raised the most questions, and were located either on the boundary of the historical and cultural heritage object or beyond it. Anomalies of artifacts were deciphered on the radargrams at two points, one located on the current territory of the tennis court, and the other beyond the boundary of the current object.
Figure 8 illustrates the interval of object AOI1 at a depth of approximately 0.68 meters (t=13.75 ns). Attention should be drawn to the darker visualization parts (higher amplitude) on the left side in the center of the object, indicating a possible disturbance that may have occurred during the construction of the tennis courts.
Figure 9 illustrates the interval of the object at a depth of approximately 0.78 meters (t=15.63 ns). The termination of this object can be observed at AOI2. A rectangular object can be seen below the center of the image. The soil outside the 1944 fence is located on the right side.
Using the newly acquired data from radar and ground-penetrating radar surveys, the boundaries of the historical and cultural heritage site, as established by the Ministry of Culture of Ukraine to date, were overlaid onto the modern topographic map from 2007, at a scale of 1:1000. Additionally, adjustments were made to the boundaries in light of the newly discovered underground artifacts (
Figure 10).
For the validation of the conceptual model for refining the boundaries of historical and cultural heritage sites and establishing protection zones around these sites, a software module has been developed and programmed. This module processes the synthetic aperture radar (SAR) data of the surface part and the ground-penetrating radar (GPR) profiles of the underground part of the object [
38,
39]. The main challenge in developing this software module was to mathematically integrate the results obtained from SAR and GPR.
Mathematical reconciliation of radar and ground-penetrating radar (GPR) data can be achieved using various data processing methods, including filtering and signal processing techniques. One approach to reconciling the results involves utilizing filtering algorithms such as the Kalman filter. This method relies on computing the system state estimation based on processed information from all sources and modeling the system using stochastic differential equations. Another approach to reconciling the data involves using signal processing algorithms like correlation analysis. This method is based on comparing signals obtained from different sources and determining the degree of their correspondence. Additionally, other mathematical methods such as principal component analysis, classification methods, and pattern recognition can be employed to reconcile the data processing results.
Regardless of the approach, when reconciling the results of radar and GPR data processing, it is essential to consider the characteristics of each method and data source and ensure the quality and accuracy of data processing.
To combine the results of spaceborne radar interferometry from the Sentinel-1 satellite and GPR imaging at a specific point, a weighted approach can be utilized (1).
where:
R - the result of combining the two methods for a specific point on the surface;
α and β - weights that determine the importance of each data source in the combined result. These weights can be determined by experts or through optimization depending on the specific application;
Ci - the measurement result obtained from spaceborne radar interferometry at a specific point;
Cg - the result of ground-penetrating radar imaging at the same point;
This approach allows combining information from both sources to obtain a more complete and reliable representation of a specific point.
The formula for combining InSAR (Interferometric Synthetic Aperture Radar) and GPR (Ground Penetrating Radar) can be complex and depend on specific use cases. However, the general idea is to integrate information obtained from both sources to improve the accuracy of interpreting and analyzing geological structures or soil properties.
A possible general formula (2) could be:
Here, f is a function that may include various data processing operations, filtering, and other analysis methods to achieve a more accurate result.
For the creation of the software module, a technological scheme was proposed, as shown in
Figure 11.
To automate the described processes, a software module was developed in the VisualBasic programming language for refining the boundaries of historical and cultural heritage sites. Although VisualBasic is less suitable for programming cartographic elements compared to, for example, Python, its advantage lies in its object-oriented nature and full integration with the Windows operating system. By using additional modules available for VisualBasic, we were able to implement the assigned task.
For testing and validation of the created software module, graphical materials on the Lviv Citadel were used in raster and vector formats, namely: Citadel's Cadastral Map in GeoTIFF format, radarogram of radar scanning of individual areas of interest of the Citadel in GeoTIFF format, vector layer of digitized radarogram anomalies in DXF vector format, and the boundary of the historical and cultural heritage site in DXF vector format.
The program execution started with testing the first functional block "Data Loading." When the command "Load 3D Model" is executed, the GeoTIFF format CAD model of the Citadel is loaded.
Next, we proceed to the next step - loading the radiogram obtained from the ground-penetrating radar survey in GeoTIFF format. Then, commands to load vectorized anomalies based on the radiograms of ground-penetrating radar survey are executed in DXF vector exchange format.
In fact, the dataset loaded in the window in
Figure 12 is final for analysis, and we can proceed to the next functional block. The first command executed in the second block is "Boundary Analysis," and its execution is hidden, with the result being loaded into the computer's memory. That is, after executing this command, the picture inside the window remains unchanged. This command performs an overlay analysis of the loaded data, after which the module determines whether all additional elements fit within the existing object boundary. If so, the next commands are unnecessary. If not, the software rebuilds the boundary according to the newly identified elements.
The next step is executing the command for boundary reconstruction (
Figure 13), as a result of which the program visualizes the outcome of the previous step from its memory in the data field. After the reconstruction of the object's boundaries according to the new elements, the protective zone of the cultural heritage object is constructed (
Figure 14).
The result of the module execution is the reconstructed updated boundary of the object of historical and cultural heritage with an automatically generated protective zone (
Figure 14). This result can then be exported to the exchange vector format DXF and opened in any geographic information system or CAD system for further work.
It is worth noting that this module is advisable to apply to already registered objects of cultural heritage based on the centroids of polygons when there is a need for directly refining the boundaries of the object itself.
5. Conclusions
The refinement of boundaries for historical and cultural heritage immovable objects and the determination of their protective zones using remote sensing methods represent crucial steps in preserving and safeguarding these invaluable assets. Through the application of remote sensing techniques, such as satellite imagery and ground-penetrating radar, precise boundaries can be delineated, allowing for effective management and conservation efforts. This conceptual approach not only aids in the identification of potential threats to heritage sites but also facilitates the planning of archaeological excavations and scientific research initiatives aimed at further understanding the historical and cultural significance of these objects. By integrating remote sensing data into the geographical information system, comprehensive models can be developed to guide conservation strategies and mitigate risks associated with urbanization, tourism, and natural disasters. Overall, the utilization of remote sensing methods offers a valuable toolset for enhancing the protection and preservation of historical and cultural heritage immovable objects for future generations.
The conceptual and detailed models are proposed for the study of cultural heritage objects using remote and non-invasive methods, and a technological scheme for their implementation has been developed. Extremes of vertical displacements of the Earth's surface within the territory of the Citadel, an object of cultural heritage in Lviv, were determined based on interferograms created using the radar interferometry method from 36 Sentinel-1 satellite images spanning from 2020 to 2022. A total of 24 such points were identified.
The obtained extremes of vertical displacements were investigated using ground-penetrating radar imaging. Special attention was paid to two areas where anomalies were observed, one of which lies beyond the currently established boundary of the object.
A software module for refining the boundaries of cultural heritage objects and determining their protective zones has been developed and programmed. Additionally, a technological scheme for implementing the module development process has been provided.