In the plain of Kilkis prefecture (Greece), 44 km north of the city of Thessaloniki (), rises the stronghold of Lazaritsa Chorygi (41°0′52.45″ N 22°42′29.59″ E, ). It is visible from tens of kilometers away and is the result of a natural geological formation of bedrock. In prehistoric and historical times, the hill served as the citadel of the settlement that stretched along its slopes. It was of particular strategic importance from antiquity to Byzantine and modern times. The archaeological research findings date from the Late Bronze Age to the Middle Byzantine period. Fragments of pottery, mainly from the Classical period, have been found, while a bronze helmet from the 4th century BC with an engraved representation of lions on its front has been found after the desecration of a grave. An important find is a small statue of the goddess Aphrodite, dating from the late Hellenistic period. Also, a very important fact is that the settlement has not been linked to an ancient site. The current “ravaged” form of the hill’s slopes is the result of the 1985 operation to create land steps for its reforestation, which destroyed the antiquities [
1,
2].
. The location of Greece in Europe and the locations of the archaeological sites of the Lazaritsa Chorygi in Greece.
. The stronghold of Lazaritsa Chorygi (<b>left</b> figure approaching from southwest, <b>right</b> figure approaching from southeast).
The first objective of this paper is the creation of a Digital Surface Model (DSM) and Orthophotomosaic of the stronghold of Lazaritsa Chorygi, as this area of great archaeological/historical interest and special natural beauty has not yet been mapped.
During the First World War (1914–1918) and after the arrival of the armies of the East in 1915, the hill was fortified by French soldiers, who must have built military trenches near the hill and in the surrounding area. After the French moved further west, English soldiers [
1] settled on the hill. Close observation of the hill and the surrounding area using Google Earth Pro (GEP) images led to the discovery of military trenches located a short distance to the northwest, north, and northeast of the hill. This was carried out in the older GEP images, as in the most recent images, the structures are mostly covered by crops or soil.
Thus, another aim of this paper is to map these military trenches, which are almost absent from modern GEP images.
2.1. Unmanned Aircraft System and Google Erath Pro
The Unmanned Aircraft System (UAS) used in this work was the DJI Mini 3 Pro, which has a 48 MegaPixels (MP) RGB 1/1.3″ CMOS sensor [
3]. The manual image trigger allows the collection of images with a resolution of 48 MP. However, in the shots that were taken, the trigger was automatically activated (e.g., using timelapse), resulting in a 12 MP resolution for all images.
Google Earth Pro (GEP) is an open-source software with time-lapse images collected by satellite or airborne sensors, used by both the public and scientists in a variety of applications (mapping forest areas, urban areas,
etc.) [
4,
5,
6,
7,
8,
9] and policymakers,
etc. [
10,
11]. GEP images are only RGB composites (lacking universal Multispectral-MS images); therefore true DNs of the original images cannot be retrieved. Nevertheless, they allow for visual data collection and exploration, such as mapping, monitoring, and managing existing archaeological sites and landscapes, as well as the identification of marks of new archaeological remains beneath the soil [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25]. Survey areas can be anywhere in the world, and their extent is not subject to constraints, thus facilitating studies that conduct research at any scale. For the paper, GEP images from 2004 were utilized as military trenches are almost absent in modern images.
2.2. Flight Plans and Image Collection
A vertical collection of images would not allow a complete mapping of the hill, as the bedrock also creates surfaces that either develop vertically or are indented. Therefore, an additional collection of inclined and horizontal images were required ().
. The various flight heights for vertical, inclined and horizontal image collection (perspective from the west).
The three flights (vertical, inclined και horizontal capture, ) with the RGB sensor of the archaeological site of Lazaritsa Chorygi took place on 23/9/2024 from 11:30 a.m.–14:00 p.m. The wind speed was 1.5 m/s (1 Beaufort) and the air temperature was 22 °C.
In the first flight (), the image capture was vertical (~90°) with 70% side and front overlap images (). A total of 9 strips were implemented and the flight height was 40 m, steadily above the highest point (275 m) of the hill,
i.e., 315 m () above sea level (it was chosen that the UAS should not follow the terrain during image collection, as the altitude at a horizontal distance of 200 m varies up to 75 m). The distances of the objects from the sensor varied from 40 m to 115 m (due to the uneven relief), the total flight time was 10 min, and a total of 99 RGB images were collected.
. (<b>a</b>) Flight diagram with the positions of the capture points (Maps Data: Google Maps, image ©2024 Airbus, Maxar Technologies, Center of figure: 41°0′52.45″ N 22°42′29.59″ E). Start of image captures from the southeastern corner with a flight direction from south to north; (<b>b</b>) the layout of the vertical images and; the 3D perspectives of the collected images (<b>c</b>) from the southwest and (<b>d</b>) from the east.
.
UAS data collection.
Flight |
Image Collection |
Overlap Images |
Strips |
Flight Height above Sea Level (m) |
RGB Images |
First |
Vertical (~90ο) |
70% side and front |
9 |
315 |
99 |
Second |
With inclination (~45ο) |
70% side and front |
14 |
315 |
56 |
Third |
Horizontally (~0o) |
90% |
- |
260 (one flight), 275 (two flights) |
220 |
In the second flight () the images were captured at an inclination (~45°) with 70% side and front overlap images ( and ). A total of 14 strips (9 strips parallel to the X-axis, , and 5 strips parallel to the Y-axis, ) were implemented, and the flight height was again, consistently, 315 m () above sea level. The distances of the objects from the sensor ranged from 55 m to 145 m (due to the uneven relief), the total flight time was 16 min, and a total of 56 RGB images were collected.
. (<b>a</b>) Flight plan parallel to the Y-axis with the positions of the capture points (Maps Data: Google Maps, image ©2024 Airbus, Maxar Technologies, Center of figure: 41°0′52.45″ N 22°42′29.59″ E). Start of image captures from the northwestern corner with a flight direction from north to south; (<b>b</b>) the layout of the collected inclined images and; the 3D perspectives of the collected images; (<b>c</b>) from the southwest and (<b>d</b>) from the east.
Finally, in the third flight (), the image capture was performed horizontally (~0°) with ~90% overlap images. Here, three independent shorter flights were conducted. The hill does not have a circular plan view;therefore, an automatic circular flight centered on the center of the hill could not be performed so that the flight path would be uniform and the distance to the monument would be constant. In contrast, the hill has an ellipsoidal plan view, and therefore, manual flights should be performed (hence, the distances from the monument are not constant). In the first flight, the overall study object (hill slopes and bedrock) was captured, the flight height was 275 m () consistently above the sea surface, and the distances of the objects from the sensor ranged from 40 m to 50 m (). In the second flight, only the bedrock was captured, the flight height was 260 m () consistently above the sea surface, and the distances of the objects from the sensor ranged from 25 m to 35 m (). In the third flight, only the hilltop was captured. The flight height was 275 m () consistently above the sea surface and the distances of the objects from the sensor ranged from 5 m to 15 m (). The total flight time for all three flights was 18 min, and a total of 220 RGB images were collected.
. (<b>a</b>) Flight plan parallel to the X-axis with the positions of the capture points (Maps Data: Google Maps, image ©2024 Airbus, Maxar Technologies, Center of figure: 41°0′52.45″ N 22°42′29.59″ E). Start of image captures from the northeastern corner with a flight direction from east to west; (<b>b</b>) the layout of the collected inclined images and; the 3D perspectives of the collected images; (<b>c</b>) from the south and (<b>d</b>) from the southwest.
. First horizontal image collection flight; (<b>a</b>) the positions of the acquisition points (Map Data: Google Maps, image ©2024 Airbus, Maxar Technologies, Center of figure: 41°0′52.45″ N 22°42′29.59″ E) and; the 3D perspectives of the collected images; (<b>b</b>) from the south and (<b>c</b>) from the southeast.
. Second horizontal image collection flight; (<b>a</b>) the positions of the acquisition points (Map Data: Google Maps, image ©2024 Airbus, Maxar Technologies, Center of figure: 41°0′52.45″ N 22°42′29.59″ E) and; the 3D perspectives of the collected images; (<b>b</b>) from the southwest and (<b>c</b>) from the south.
. Third horizontal image collection flight; (<b>a</b>) the positions of the acquisition points (Map Data: Google Maps, image ©2024 Airbus, Maxar Technologies, Center of figure: 41°0′52.45″ N 22°42′29.59″ E) and; the 3D perspectives of the collected images; (<b>b</b>) from the south and (<b>c</b>) from the southeast.
From the above paragraphs, it is evident that image collection from different flight altitudes and various capture angles can collectively ensure comprehensive photographic documentation of the monument. The determination of the model’s spatial resolution and accuracy, as well as the completeness of its thematic information, remains pending.
2.3. Image Processing and Product Production
All images were imported into the Agisoft Metashape Professional© version 2.0.3 software (). The UAS has, among others, a Global Navigation Satellite System (GNSS) (compatible with GPS, Galileo, and BeiDou), which allows the recording of the coordinates (X, Y and Z) of the center of each image collected during the flight, with an accuracy of about 3 to 5 m. The coordinates were automatically converted from World Geodetic System 1984 (WGS84) to Greek Geodetic Reference System 1987 (GGRS87) in the Agisoft Metashape Professional© version 2.0.3 software, the aligning of images (align photos with high accuracy) was performed and at the same time a sparse point cloud model based on matching pixel groups between images was generated. Meanwhile, the Root Mean Square Error for X, Y and Z coordinates (RMSE
XYZ) for all the sensor locations was calculated [
26]. Then, the dense point cloud was created (build the dense cloud with high-quality and aggressive depth filtering). Next, the 3D mesh generation (build mesh) followed, where the point cloud was transformed into an actual 3D surface. The following step was conducted to build the texture (build texture),
i.e., the colored overlay of the generated 3D mesh. The last step was to generate a DSM and orthophotomosaic. In summary, RMSE
XYZ is 2.42 m, 67,846,306 points (Point Cloud), 6,924,069 faces (3D Model), a DSM with spatial resolution of 5.6 cm and an orthophotomosaic with a spatial resolution of 2.8 cm were created ().
. Workflow diagram for image processing using Agisoft Metashape Professional© version 2.0.3 software.
. (<b>a</b>) DSM and (<b>b</b>) orthophotomosaic (center of figures: 41°0′52.45″N 22°42′29.59″E).
2.4. Google Earth Pro
After collecting images using the UAS, the longitudinal Google Earth Pro© (GEP) version 7.1.8. 3036 images were studied for both the hill and the wider area. The visual observation of the northwest, north and northeast of the hill led to the identification in the oldest temporal images of the GEP (images for the years 2003 or 2004) of obvious remains of military trenches, which are now either completely absent or have been mostly covered by crops or soil (). There are no crop or soil marks in the contemporary images that would indicate the presence of covered structures. However, it is hypothesized that such structures may exist beneath the crops or soil.
. The location of military trenches in the GEP image; (<b>a</b>) of 2003 (Map Data: Google Earth, image Landsat/Copernicus, image ©2024, Maxar Technologies) and; (<b>b</b>) of 2023 (Map Data: Google Maps, image Landsat/Copernicus, image ©2024, Airbus). The almost total absence of structures is evident in the most recent image.
Initially, for the years 2003 or 2004, images of parts of the wider area were collected with a resolution of 8192 × 4771 pixels and at a suitable viewing scale that allows the collection of the spatial resolution available in GEP (from 0.5 to 1 m for the study area) of the images [
20]. A mosaic of all the images was then created and finally, the geometric correction (1st degree polynomial transformation) of the mosaic in GGRS87 was performed, utilizing the X and Y coordinates of the reference points (20 points were used in total) retrieved from the National Cadastre website [
27] (coordinate accuracy of about 1 m). The geometrically corrected mosaic has a final planimetric accuracy of 1.5 m. All processing operations described in this section were performed using Erdas Imagine 2015© version 15 software.
According to the above paragraphs of the image processing in Agisoft Metashape Professional© version 2.0.3, the hill slopes and the bedrock were mapped using 68 million points spaced approximately 5.5 cm apart. In addition, 7 million faces were created that fully capture the geometry of the objects with exceptional texture. At the same time, it is found that all vertical surfaces and recesses were captured (). Therefore, while the RMSE
XYZ is 2.42 m (which is acceptable at first), the bedrock and hill slopes were captured with high spatial resolution and completeness.
a shows the wider area to the northeast, north and northwest of the hill. b shows the relief of the study area (source: National Cadastre, Digital Terrain Model with spatial resolution 5 m), featuring the military trenches. They have a total length of 9 km.
. (<b>a</b>–<b>d</b>) 3D models of the hill and their indicative figures with (<b>e</b>) Points; (<b>f</b>) Mesh; (<b>g</b>) Textures.
. (<b>a</b>) The location of the military trenches. Southwest of the center of the figure is the location of Lazaritsa Chorygi (Map Data: Google Earth, image Landsat/Copernicus, image ©2003, Maxar Technologies); (<b>b</b>) The relief of the terrain in the study area (light tone of grey at high altitudes to dark tone of grey at lower altitudes), and the design (in yellow) of the military trenches. In the red frame, the areas of the military trenches are presented in the following figures. Center of figures: 41°1′19.80″ N 22°43′33.77″ E.
In this paper, Ground Control Points (GCPs) were not measured in the field or used during digital image processing to produce DSM and orthophotomosaic. This will be carried out in the next stage of the study (GCPs measurements, not new image acquisitions as follows), ultimately producing the above products with high spatial accuracy [
28]. This study will also include the use of Lidar in UAS [
29,
30,
31,
32,
33,
34,
35,
36] to optimally capture the military graves covered by the crops and determine their degree of preservation. Thus, it is important to determine both the completeness of the monument’s mapping images and whether the mapping was done with high spatial resolution (), as well as whether the results meet a reasonably acceptable level of spatial accuracy.
In addition to the fact that the bedrock and hill slopes were captured with high spatial resolution and completeness (mapped using 68 million points, spaced approximately 5.5 cm apart, and 7 million faces created to capture the geometry of the objects with exceptional texture fully), it was found that all vertical surfaces and recesses were recorded. As a result, while high resolution orthophotomosaics can be generated, it is also possible, using the 3D model, to identify areas in the bedrock that may be at risk of rockfalls in the future (and propose preventive structural reinforcements). Military trenches act as a barrier against invasions from the north (the direction of invasion is derived from the relief of the area, ). , and show, by way of illustration, the positions 1–3 of b. The types of military structure typology [
29,
30,
31,
32,
37,
38,
39,
40] observed are (e.g., b) fire trenches, communication trenches, shelters, front and support trenches, and strong points.
More specifically, in , we observe that the structures follow the relief of the terrain (light grey tones at higher altitudes transitioning to dark grey tones at lower altitudes, with a contour interval of 7.5 m). On the front line are the fire trenches. Parallel and at a short distance, the communication trenches are developed. At a second level and a higher altitude, the support line is observed with new fire trenches. The shelters are developed perpendicular to and between the above-mentioned fire trenches. A strong point, protected by fire trenches, is observed at one location in the southwest. In contrast, in , no strong point is observed, but successive fire trenches (following the terrain’s relief) are interconnected almost perpendicularly by communication trenches. Similar observations (as in ) are made in , with the difference that in this area, multiple successive fire trenches are identified (four parallel fire trenches are visible). In the remaining areas of that are not shown, similar spatial arrangements and construction typologies are observed, with minor variations, as seen in , and .
. (<b>a</b>) The military trenches of position 1 of b (Map Data: Google Earth, image Landsat/Copernicus, image ©2003, Maxar Technologies); (<b>b</b>) The relief of the terrain (light tone of grey at high altitudes to dark tone of grey at lower altitudes, the interval between contours is 7.5 m) and design (in yellow) of the military trenches. Center of figures: 41°1′27.04″ N 22°42′38.50″ E.
. (<b>a</b>) The military trenches of position 2 of b (Map Data: Google Earth, image Landsat/Copernicus, image ©2003, Maxar Technologies); (<b>b</b>) The relief of the terrain (light tone of grey at high altitudes to dark tone of grey at lower altitudes, the interval between contours is 7.5 m) and design (in yellow) of the military trenches. Center of figures: 41°1′27.26″ N 22°43′10.36″ E.
. (<b>a</b>) The military trenches of position 3 of b (Map Data: Google Earth, image Landsat/Copernicus, image ©2003, Maxar Technologies); (<b>b</b>) The relief of the terrain (light tone of grey at high altitudes to dark tone of grey at lower altitudes, the interval between contours is 7.5 m) and design (in yellow) of the military trenches. Center of figures: 41°1′21.40″ N 22°43′35.85″ E.
Due to the complex geometry of a monument, the process of collecting images for its mapping needs to be adapted. For the optimal mapping of the stronghold of Lazaritsa Chorygi (Greece) and its slopes, besides vertical images, inclined and horizontal images were collected. Thus, for a monument that had not been mapped until now, a rich data set and final products (DSM and orthophotomosaic) of high spatial resolution have become available. In the next phase, GCPs measurements will be carried out to produce DSM and orthophotomosaic of high spatial accuracy.
Unfortunately, in most countries where military structures were built during the Great War, there is no care for their preservation. These structures are being sacrificed when there is a need for housing development, construction, or agricultural exploitation. Thus, while we are aware of the problem today as we try to trace the landscape of the great battles of antiquity (e.g., the battle of Philippi in 42 BC, Greece, for the succession of the Roman Empire after the assassination of Julius Caesar [
41]), future generations, less than 100 years from now, will not be able to trace the traces of the military structures of the Great War. On the other hand, France immediately after the Great War designated a Zone Rouge (Red Zone) with restrictions on building, agriculture and forestry. As a result, the surface features of the Great War are, for the most part, extremely well preserved to this day [
29,
33]. In the wider area of the Lazaritsa Chorygi stronghold, Great War military structures of 9 km length were identified, which can be observed in the 2003 or 2004 GEP images, but worryingly are almost universally absent from modern images. In a next phase, Lidar in UAS will be utilized to optimally capture military trenches covered by crops or natural vegetation and determine their degree of preservation. At the same time, a multispectral sensor on UAS will be used to identify the marks of the covered military structures from the soil or crop [
28].
Thanks to Georgia Stratouli (Head of the Ephorate of Antiquities of Kilkis, Greece) and Nektario Poulakaki (Head of the Department of Prehistoric and Classical Antiquities and Museums of the Ephorate of Antiquities of Kilkis, Greece) for permission to collect data of the stronghold of Lazaritsa Chorygi.
Not applicable.
Not applicable.
No original images or raw data will be made available on the locations, as they concern archaeological/historical sites.
This research received no external funding.
The author declares that he has no known competing financial interests or personal relationships that could have appeared to influence the work re-ported in this paper.