1. Introduction
All of the European Union Member States are obliged to monitor the conservation status of natural habitats and species, as well as to undertake active protection measures to prevent deterioration of the status of habitats (EEA) [
1]. A common (European-wide) formal methodological premise relating to the monitoring of the Natura 2000 network, is the methodology for monitoring habitats listed in Annex I of the Habitat Directive (EC) [
2] at the level of the biogeographic region [
3]. The currently used traditional methods of monitoring Natura 2000 habitats and identifying the main threats are based on expert assessment carried out directly on the ground. Due to the increasing potential of photogrammetric and remote sensing technologies, the availability of remote sensing sensors and airborne platforms, and the intensive development of data analysis methods (data mining, machine learning), it is reasonable to develop tools to optimize current methods of identification and assessment of the conservation status of Natura 2000 habitats.
One of the threats to non-forest Natura 2000 habitats, significantly affecting their conservation status and prospects for their preservation in Poland, is secondary succession. This is related to the gradual cessation of the use of agricultural land that has occurred since the early 1990s. Currently the process of secondary succession can be analyzed using airborne laser scanning (ALS) data [
4,
5], however, obtaining knowledge about the past dynamics of this phenomenon (e.g., from the 1960s onward) is only possible using archival materials—mainly aerial photographs, which are often the sole reliable source of information about the condition of land cover in the past.
Until now, the most common method of obtaining information on the extent of trees and shrubs based on archival aerial photographs has been photointerpretation [
6,
7,
8,
9,
10,
11]. Another method, albeit less frequently used in this kind of work, has been stereodigitalization [
12,
13,
14,
15,
16,
17]. However, both methods are highly time-consuming, and hence, the search for possibilities for automation of the process of obtaining the extent of trees and shrubs, based on these photographs. The development of digital processing algorithms and image analysis has brought a number of new possibilities in this field. Spectral classification is a frequently used solution [
18,
19], but because it usually requires several spectral bands, it is not a universal method and does not allow for the analysis of archival (digitized) analogue photographs (often gray scale images—panchromatic or infrared). In subsequent decades, color photographs in the optical (RGB) or color-infrared (CIR) range have become available, but they also contain only three spectral bands, which usually do not give satisfactory results in the pixel-based classification.
In the age of information technology and photogrammetry development, dense image matching (DIM) techniques could be used to automate this process. The most promising technique in recent years—semi-global matching (SGM) [
20,
21,
22,
23,
24]—has enabled obtaining xyz point clouds, characterized by high position accuracy, based on photogrammetric images [
25,
26]. These two algorithms are often compared [
27,
28]. An attempt at description of the state-of-art in image matching and the efficiency of DSM generation, with the use of different algorithms in commercial software packages, was carried out by EuroSDR [
29]. However, in the case of archival photographs, several conditions must be met in order to use these data. The previous study [
30] shows that the most important aspect is the scale of the photographs, which should not be less than 1:13,000 for analogue photos, or ground sampling distance (GSD), which should be less than 25 cm for digital images. In addition, the date of aerial photographs acquisition has a large impact on the quality of the final product—the trees and shrubs should be fully leafed when images are collected.
In the case of analysis going back many years, when the analysis is based on archival analogue photographs, the first stage must be the digitalization of analogue photos. Correct scanning plays a large role in the context of the usefulness of archival photos for further analysis. The next step is to define exterior orientation of the images in bundle adjustment (aerotriangulation). At this stage, it is necessary to know the coordinates of control points, i.e., the signaled points in the field visible in the images. In the case of archival images, there is often a lack of such information, which may impede proper aerotriangulation [
31]. The solution may be to measure characteristic points on an already-oriented model from later years, or on the ground, if it is possible to access an easily identifiable object that has not been moved over time. The very process of generating a dense cloud of points, which is then used to build a digital terrain model (DTM), can be carried out in various ways, depending on the algorithm used. There are also many types of commercial software dedicated for aerial images and for low-altitude images (e.g., Agisoft Photoscan, Pix4D, Trimble Inpho, Socet Set, LPS), as well as several open source tools (PRiME Stereo Match, MICMAC Tool, Visual SFM, StereoSGBM from the open source library OpenCV), that enable a creation of a dense cloud of points from archival photos. The algorithms implemented in them sometimes act differently, which has been the subject of many publications [
32,
33].
In practice, for the detection of the extent of trees and shrubs, the canopy height model (CHM), which is the difference between digital surface model (DSM) and digital terrain model, is the one most commonly used [
34]. In the case of entire areas, including also those that are not covered by vegetation, we can refer to a normalized digital surface model (nDSM). Then, all the areas which are above the assumed height threshold on CHM are cut off, with an assumption that the higher objects are vegetation. Of course, such an assumption may be adopted only for undeveloped areas. The most commonly used cut-off threshold is 2 m [
21,
35,
36,
37,
38].
The main advantage of point clouds generated by DIM is the exact geometric description of the top surface of the stand [
39]. In turn, their disadvantage is the lack of height measurements under the tree crowns level, which makes it difficult to correctly determine the height of trees. It is usually necessary to obtain DTM from another source. However, this is less important in assessing the extent of trees and shrubs, which is why, in areas characterized by flat terrain, point clouds filtering can be used to generate either DTM or ALS data, because they accurately represent the terrain in the forest areas where the variability of the terrain surface is not dynamic in time.
The main objective of our study was to assess the possibility of using digital height models, created using DIM techniques, to determine the extent of trees and shrubs, based on archival aerial photographs. As mentioned above, to determine the sites occupied by trees and shrubs, after the generation of CHM, a height threshold, above which the areas will be defined as being covered by medium and high vegetation, should be adopted. Because in the previous studies [
21,
35,
36,
37,
38], forest areas have been studied and the focus has not been on the earlier stages of the secondary succession process (trees and shrubs lower than 2 m), this study has attempted to analyze the impact of the CHM cut-off threshold on the effectiveness of determining the extent of the occurrence of smaller trees and shrubs. This article also analyzes the possibilities of using two types of software for point clouds generation, and consequently DSM and nDMS for vegetated areas (CHM): Agisoft Photoscan and Trimble Inpho.
In the next two sections, the study area, the methodology and the datasets used in the experiment are introduced. In the fourth section, the results of point clouds quality analysis (
Section 4.1) that originate from applying various approaches and different software commonly used in vegetation monitoring are compared, with respect to ALS data, to achieve the best results of surface modeling of vegetation in the succession area. The second part of this section (
Section 4.2) is the analysis of using a DSM obtained with a selected method for mapping the area of succession, considering data sources derived from DIM.
2. Study Area
The research was carried out in the Olsztyńsko-Mirowska Refuge, a NATURA 2000 protected area (PLH240015), which is located in the south of Poland, in the Silesian Voivodeship, near Czestochowa city (50°45” N; 19°17” E) (
Figure 1). The Olsztyńsko-Mirowska Refuge is an enclave of natural and semi-natural ecosystems among the highly urbanized areas of Silesia and Częstochowa industrial districts. This area is situated entirely within the Orle Gniazda Landscape Park. The refuge includes a group of limestone hills (moors) with numerous karst forms, such as caves, monadnocks, and karst hoppers. This area is characterized by large habitat diversity. The most important habitats are non-forest habitats associated with limestone rocks with numerous, rare and endangered, thermophilic species of plants and invertebrates (including the species listed in Annex II of the Habitats Directive—the Scarce Large Blue (
Phengaris teleius)). A number of species here reach the northern end of their range. In total, in the Olsztyńsko-Mirowska Refuge, there are 13 types of habitats identified in Annex I of the Habitats Directive (including 4 priority habitats) and 12 species of plants and animals listed in Annex II of the Habitats Directive [
39].
The threat to the habitats of this area is primarily the process of secondary succession, which is the result of the abandonment of pastoralism, with the consequential lack of grazing and mowing of meadows, which have gradually ceased since the 1990s. In order to protect valuable habitats in the Olsztyńsko-Mirowska Refuge, under the LIFE11NAT/PL/432 project “Protection of valuable natural non-forest habitats typical of the Orle Gniazda Landscape Park” (2012–2016), a number of active protection measures have been implemented, such as grubbing of shrubs and trees and controlled sheep grazing [
40].
The part of the Olsztyńsko-Mirowska Refuge selected for analysis covers a region which is varied in terms of terrain relief and land use, and the character of the secondary succession process, which has enabled a number of possible problems to be identified that can occur in the automation of the process of determining the extent of trees and shrubs succession. The study area encompasses an area of about 25 km2.
4. Results
This section presents the results of assessments for DSMs derived from DIM point clouds based on reference data from ALS and evaluation of selection of threshold height values in DSM to define a reliable range of medium and high vegetation.
4.1. Quality Assessment of the DSMs Derived from Point Clouds
In this section, the completeness and correctness of the DSMs derived from point clouds are presented.
The analysis of the quality of the products is divided into two parts:
the verification of the correctness of the height of individual trees, easily distinguishable in archival photographs in DIM-DSMs compared to the height recorded in the ALS data;
the verification of the correctness of the height of individual trees growing in high canopy closure in DIM-DSMs, compared to the height recorded in the ALS data.
In order to assess the usefulness of the obtained point clouds for the delineation of trees and shrubs, the following test was carried out. For two sets of points representing vegetation, the heights were determined based on the normalized digital surface model (nDSM) obtained on the basis of DIM-DSM from several variants and DTM from ALS. The example of this analysis for both datasets is shown in
Figure 3, where nDSMs and orthophotomaps are presented with some reference objects in three sets of photos with different acquisition times.
For each point (specimen), the ratio of the DIM-nDSM value in each of the examined periods to the reference value from the ALS-nDSM from laser scanning acquired in 2016, was expressed as a percentage (formula 1):
where:
v—variant of DIM-nDSM,
d—data acquisition date.
Due to the natural succession of vegetation, it is quite difficult to have reference heights from the past. It must, therefore, be assumed that with the analysis of increasingly older sets of data, the percentage of height will decrease. In this analysis, the proper result is the value indicating the existence of a specimen classified as, a) a single tree, or b) a tree in a high canopy closure.
The range of values of the reference height percentage that qualified as a positive detection of a specimen in the DIM-nDSM model is specified in
Table 3. The values were selected based on data observation, as well as consultation with a research team of botanists.
The analysis of the results revealed a significant variation in the obtained results. The result of tree detection based on the elevation model was not always visible, despite the tree outline being clearly noticeable on the orthophotomap (
Figure 3). The results of tree detection for each date and point cloud generation variant are included in
Table A1 and
Table A2. The visualisation of these results is presented in
Figure 4. In the case of individual trees, detected based on nDSM from point clouds generated in Trimble Inpho, in the vast majority of cases the best result was given by variant v3, but the differences were very slight. This allows us to conclude that it was increasing the density of the point cloud that contributed the most to obtaining the good result. Detection values, however, range from 18% to 99% for different datasets. The date of image acquisition has the biggest impact on such a result. This is best demonstrated by the images from 2012 and 2015 having similar GSD (24–25 cm), acquired with the same camera (UltraCamXP). The detection results are 22% and 99%, respectively, while the date of the flights is 25 March and 8 August. In March, the vast majority of deciduous trees do not usually have developed leaves and the trees’ stems and branches are hardly noticeable in the image. The lighting conditions and the contrast of the objects appearing in the image also have a significant influence on the correctness of the generated point cloud, as well as the GSD/photo scale. As regards the “high” variant, in Agisoft Photoscan, it returned the same or better result than the “medium” variant for all of the years, except 2003. It also typically had a worse detection rate, by a few to several dozen percentage points, in comparison to the v4 variant in Trimble Inpho which was characterized by similar point cloud density.
In the case of trees with high canopy closure, very good results were also obtained in Trimble Inpho. Approximately 55% to 100% of objects were qualified as properly detected, following the thresholds presented in
Table 3. Variant v3 provided a much better result than the others, only for the first date of image acquisition. However, the difference was almost 20%. The worst results were observed for images from 1982 and the best for 2015. The newest images allowed for a correct detection of single trees, however, v3 was the only variant that returned the average height equal to 100% of the height from ALS. Even though the v3 variant was not always the best in the detection of trees with high canopy closure, it always provided the highest average height for detected trees. For 2015, it was 100% accurate and it can be assumed that it would always be the closest to the truth, as bigger density allows for better treetop modeling.
In the case of Agisoft Photoscan software, the higher density in “high” approach did not always provide a better result. A big difference was noted in the case of high canopy closure in the images from 2003 and 2009, which can be clearly seen in
Figure 4. The algorithm was not able to determine characteristic points for the larger areas of forests, with too high a resolution of the photos and many areas looking similar. As a result, the point clouds sometimes had holes in the place of forests (no data). Reducing the resolution by 4 times (the quality: “medium” instead of “high” parameter) helped to diminish this problem. Thus, one can observe the advantage of the algorithm used in Trimble Inpho over that of Agisoft Photoscan, as in the first algorithm the problem occurred on a much smaller scale.
Analyzing
Figure 4, it can be seen that the detection of trees was strongly influenced by numerous factors related to image parameters or season of image acquisition. The result of trees detection in point cloud and elevation models was not easy to predict. Despite the right image parameters and the right season of image collection, fewer single trees were detected in the image from 2003 than in the 2015 image. Such a problem with the detection of trees was negligible in the case of trees growing in high canopy closure (i.e., forest, copse). Analyzing the average height of detected trees showed that the growing trend was much more visible for trees in high canopy closure. The trend was defined by determination coefficient of 0.92, while for single individual trees the determination coefficient was only 0.64 in the case of v3 variant. Although these two datasets cannot be compared, due to different conditions (light availability, neighboring plants, soil, etc.), in the case of both sets, certain trends can be noticed. The problem of using dense image matching was simply its failure to detect single trees and shrubs, which is of the first signs of vegetation succession.
4.2. Accuracy Assessment of the Extent of Trees and Shrubs for Different Threshold Values
In order to determine the impact of the cut-off threshold on the correctness of determining the extent of trees and shrubs, threshold values of 1.0, 1.25, 1.5, 1.75 and 2.0 m were tested. Visualizations of the results for each date for all analyzed area are presented in the figures included in the
Figure A1,
Figure A2,
Figure A3,
Figure A4 and
Figure A5. The obtained accuracy coefficients for determining the extent of trees and shrubs for each threshold value are summarized in
Table 4, and presented in
Figure 5. It contains a comparison of the following metrics: overall accuracy (OA), recall, precision, Cohen’s kappa coefficient and F1-score. In addition, the area size resulting from errors of omission (EO) and commission (EC) was analyzed (
Table 5). A comparison of the size of these areas with reference to the trees and shrubs reference masks (in %) is shown in
Figure 6.
As a result of the analysis of the calculated parameters of accuracy and the visual assessment of the extent of trees and shrubs obtained for each variant of the height cut-off threshold, there was no compliance of the formal assessment and visual assessment. The largest problem were the areas of rocky outcrops and inliers, which were incorrectly detected as trees and shrubs, which resulted in difficulties of indicating the optimal cut-off threshold value. Therefore, it was decided to reassess the accuracy of determining the extent of trees and shrubs after excluding these problematic areas. The comparison of results for both variants is presented in the form of graphs in
Figure 5. Their analysis shows that the obtained accuracy increases for the subsequent analyzed threshold values and this is true for each of the archival materials examined, with the exception of the photographs from 2009, which were acquired in the leafless period. The cut-off threshold providing the best results depends on the date (parameters and quality) of the analyzed archival materials. The greatest influence of the threshold values on the accuracy of the extent of trees and shrubs determination was observed in older photographs from 1971 (BW, 1:18,000) and 1996 (RGB, 1:26,000), where significant overestimation of trees and shrubs areas is visible. For photographs from 1971, the OA varies in the range of 0.850–0.903, kappa 0.512–0.663, F1-score 0.587–0.707, and for data from 1996, the OA varies in the range of 0.833–0.927, kappa 0.546–0.742, and F1-score from 0.648 to 0.748. The highest accuracy was obtained with 1.75 m and 2 m thresholds (
Figure 5). For both of these thresholds, the smallest area of erroneously detected trees and shrubs was also recorded (
Figure 6).
With regard to photographs from 2003 (RGB, 1:13,000) and 2009 (RGB, 24 cm), it can be seen that OA, kappa and F1-score are at a similar level for all the threshold values (respectively for 2003 OA: 0.914–0.924, kappa: 0.776–0.793, F1-score: 0.833–0.844, for 2009 OA: 0.898–0.902, kappa: 0.760–0.772, F1-score: 0.831–0.842). However, the smallest area of erroneously detected trees and shrubs was obtained only for the threshold of above 1.5 m for data from 2003, and below 1.25 m for data from 2009. Analyzing the results obtained on the basis of the latest data from 2015, it was found that OA, kappa, F1-score assume values at a similar level starting from the cut-off threshold of 1.25 m (OA: 0.936–0.942, kappa: 0.867–0.878, F1-score: 0.921–0.926). The smallest area size of incorrectly detected trees and shrubs was found at the 1.75 m threshold.
5. Discussion
In this section the results are discussed, with reference to investigations reported in the literature on the use of dense image matching in DSM for defining the extent of trees and shrubs.
5.1. Analysis of Results Obtained Using Agisoft Photoscan and Trimble Inpho
Based on the analysis of the results, it can be stated that generating point clouds based on archival images matching enables, in some cases, a detection of individual trees, and very often trees in forest stands with high canopy closure. Such detection is important because without it, trees and shrubs areas would be negligible in a DSM model. The results of the experiment aimed at the detection of single trees varied much, from a few percent to almost 70% for images taken in proper conditions and at appropriate times, while for trees in high canopy closure, the results of detection were better: from 40% to almost 100% of trees were detected from DSM. Such varied results are caused by different dates of photograph acquisition, the different image quality and algorithm of image matching. In relation to the obtained results, it can be stated that the photos with a scale of about 1:13,000 and digital images with a GSD equal to 25 cm acquired at the time of vegetation development (May—September) are the appropriate material for the analysis of vegetation succession.
When recommending parameters for the analysis of forest stands with high canopy closure, similar guidelines should be taken into account as those given when analyzing individual trees and their detection. However, a scale of less than 1:26,000 and a date beyond the development of vegetation do not necessarily exclude an analysis from use, as evidenced by the results obtained for the photographs from 1996 and 2012, respectively. The scale or GSD value consequently define the density of a generated point cloud. Denser point clouds provide a better chance of detection of a tree, as well as allowing for a more accurate, high measurement. Such a statement is also true for LiDAR-based DSM resolution [
45]. Interpolation in DSM generation can also have an impact on model quality [
28].
To generate clouds, it is much better to use professional photogrammetric software equipped with more powerful algorithms of dense image matching. In this research, Trimble Inpho with Semi Global Matching algorithm and elements of feature-based matching provided better results in the strategy of point cloud generation, which were more robust regarding problems with similarity and homogeneity of area. Due to the high resolution of point clouds generated in Agisoft Photoscan, areas with no data appeared. However, both approaches with Agisoft software provided slightly worse results than all the approaches with Trimble Inpho.
5.2. The Influence of Cut-Off Height Thresholds on Determining the Extent of Trees and Shrubs
As already mentioned, the difficulty in automating the procedure of generating the extent of trees and shrubs with the use of DIM technique, lies in the appropriate selection of the cut-off threshold, so as not to lose important information about the occurrence of earlier stages of succession. In various studies, the height cut-off threshold of 2 m [
21,
35,
36,
37,
38] or 3 m [
46,
47] has been used, but these studies have examined forest areas and have not focused on the earlier stages of the succession process. In addition, the majority of scientific works have concentrated on the processing of newer aerial photographs acquired by digital cameras, not by analogue cameras. Therefore, it was decided that this study would analyze the impact of various height cut-off thresholds on the correctness of determining the extent of forested and bushy areas for each archival material.
As a result of the analysis, it was found that for each archival material, slightly different results were obtained. In the case of the oldest data—from 1971, at a scale of 1:18,000—the highest accuracy (F1-score = 0.707, and kappa = 0.663) was achieved at the threshold of 2 m. However, 1.75 m is the threshold at which the total area of places mistakenly indicated as trees and shrubs, as well as omitted places, is the smallest (
Figure 6). At this threshold, the F1-score is 0.653 and Cohen’s kappa is 0.630. Unfortunately, the total EO + EC area for this threshold is close to 68% of the reference mask, and at the threshold of 2 m, it is 65% of the reference mask, with EC errors being the main factor. EO errors are ca. 5–10%. Such large EC errors are mainly due to the radiometric quality of photos from this period. In some places, overestimation of the trees and shrubs area resulted from the presence of small rocks, as well as from the presence of cereals on the arable fields (photos registered in July).
The greatest impact on the accuracy of determining the extent of shrubs and trees based on photographs from 1996 was their scale (1:26,000). The highest F1 value was obtained for the 1.75 m threshold, and Cohen’s kappa for the 2 m threshold. The total area size resulting from EO and EC errors was the smallest at the threshold of 2 m, but the difference between the area obtained with this value and the value of 1.75 m was small (35% vs. 36% in relation to the reference mask). The EC error for the threshold of 1.5 m was less than 20% of the reference mask, and the EO error was within 20–25% (
Figure 6). In
Figure 7, the effect of decreasing the threshold value on the increase of EC errors can be seen, especially in the area with rocky outcrops and inliers. In the case of these photographs, it should also be mentioned that due to the smaller scale, their photointerpretation (preparation of reference data) was difficult, so some of the EC errors might have resulted from the incompleteness of the reference data.
The next photos—from 2003, 2009 and 2015—were data of similar parameters of scale/GSD. Using the photographs from 2003, high values of F1 for the 1.25 m cut-off threshold (F1 = 0.840) were obtained. For higher height cut-off values, the F1-score increased slightly (up to 0.844 for the 1.75 m and 2 m thresholds). The same trends were observed for Cohen’s kappa coefficient values. The total area EO and EC for the threshold values from 1.25 m was 28% of the reference area (
Figure 6). EO and EC errors were at a similar level (several percent of the reference mask area), and none of them had a dominant influence on the accuracy assessment. This was due both to the scale of the photographs (1:14,000) and their radiometric quality, as well as the date of obtaining the photos (May). During this period, cereals in the fields, as well as natural herbs, had not yet grown high. Similar observations were made for data from 2015. The only difference was that the error values of EO and EC were two times lower (5–6% of the surface of the reference mask). The highest F1 value for the 1.75 m cut-off threshold (F1 = 0.926) was obtained for 2015. Nevertheless, for the threshold of 1.5 m, F1 was only 0.002 lower, and for 1 m, it was lower by 0.005. Analyzing other accuracy parameters, as well as visually assessing the obtained extent of trees and shrubs, one can assume that the 1.5 m threshold—in the case of 2015 data—allowed for obtaining the best results. Errors from revaluation resulted mainly from the late date of acquiring these images (August), and consequently, the fact that vegetation on arable fields and herbaceous plants had already grown high, so at the height threshold of 1 to 1.5 m, some patches of this vegetation could have been detected as shrubbery areas.
The data from 2009 were acquired in Spring, when deciduous trees and shrubs were not fully leafed. As a result, the highest value of F1-score (0.842) was obtained for the 1 m cut-off height threshold. Cohen’s kappa coefficient for this threshold was also the largest. This was the only period when the accuracy of determining the extent of trees and shrubs clearly decreased with the increase of the threshold values. In this period, larger EO errors (around 20% of the reference mask) than EC errors were observed. This is obviously due to the date of acquiring photos—herbaceous plants had not yet grown, so did not affect the overestimation of results, while the crowns of trees and shrubs had poor foliage.
Figure 6 clearly shows the influence of the quality of archival photos on the obtained results. The worst results (the biggest EO + EC errors) were obtained for 1971 (B/W, 1:18,000, low radiometric quality, with scratches on the film), and then for data from 1996 (RGB, 1:26,000). The best results were achieved for photos from 2015 obtained with a digital camera (RGB, GSD = 25 cm). Slightly lower accuracy was achieved with photographs from 2003 (B/W, 1:13,000) and 2009 (RGB, 1:14,000). The ratio between the EC and EO areas was different for both of these terms, which resulted from the phenological period in which the images were acquired.
Figure 7 shows, for a fragment of the analyzed area, the impact of the cut-off threshold values on determining the extent of trees and shrubs. A worse detection of low and relatively small trees and shrubs can be observed, which affects the error of omission (EO). For each of the dates, errors are evident that result from the inclusion of rock outliers and inliers in the designated extents. When CIR photos are examined, spectral information can be used to remove non-vegetative objects. However, one must bear in mind the fact that some of the rocks could be covered with vegetation and it will not always be possible to eliminate them. In such a situation, one should use the current ALS data from the leafless period (early spring), after a detailed (manually supported) point cloud classification, which would make it possible to obtain an accurate DTM.
5.3. Comparison of the Results of the Best Variants with Reference Data
As the final spatial extents of trees and shrubs for each period were selected, those extents were characterized by high accuracy (F1-score) and a small total area EO + EC. The resulting extents of trees and shrubs obtained for each dataset using DIM techniques for part of the study area is shown in
Figure 8. The comparison of the correctly indicated area with the area resulting from errors of omission (EO) and error of commission (EC), for each of the analyzed periods, is presented in the graph (
Figure 9). The graph shows results in the form of absolute (
Figure 9a) and relative values (
Figure 9b). The increase in the area covered by succession vegetation is clearly visible. In addition, the highest error of commission is observed for data from 1971—as much as 42%. In subsequent periods, the surface rate resulting from excess errors (EC) ranges from 5% to 14%. The differences in this error rate resulted from the scale of photos (small scale for 1971 and 1996) and the quality and dates of acquiring photographs (1971). The impact of the period of time is also visible in the erroneous underestimation of the area occupied by trees and shrubs, and this can be seen particularly for the data from 2009, where the underestimation error is 26% (data was obtained in early spring).
Analyzing
Figure 9, it can be concluded that the error of determining the wooded and shrubbed areas using DIM algorithms can be estimated at 15–30% depending on the date and technical parameters of archival photos. A lower rate of errors is observed in the case of photographs acquired by digital cameras, which can be seen when comparing the results obtained for 2015 images (registered using a digital camera) and 2003 photos (registered using an analogue camera) (
Figure 8).
6. Conclusions
In the presented research, an assessment of the possibility of using DIM algorithms to delimitate the extents of forested and shrubbed areas, based on archival aerial photographs with different technical parameters, was carried out. The research was carried out using this technology for monitoring the dynamics of the process of succession of trees and shrubs. In order to determine the effectiveness of DIM techniques, the algorithms implemented in Agisoft Photoscan and Trimble Inpho software were tested. Finally, the effect of nDSM height threshold on the correctness of determining the extent of trees and shrubs was investigated.
By analyzing the contemporary quality of aerial photographs captured with a digital camera, it is possible to detect an object using nDSM derived from DIM data, with a very high probability of reaching nearly 100%. The DIM algorithm, including feature-based approach eliminating errors related to homogeneity and similarity of the area, also allows for determining the correct height of trees, with a very high accuracy of only a few percent difference. These conclusions apply both to trees growing singly, and in stands with high canopy closure. The situation is slightly different as regards archival photos. The quality of detection and determining the height of trees is affected by the quality of photos and their parameters, which in the past were less favorable, as well as the period of obtaining data—unfavorable in a leaf-off situation and favorable for leaf-on. These circumstances significantly reduce the ability to detect individual trees. Detection of trees in high canopy closure was not influenced by these factors. Algorithms of image matching also allowed for a reliable determination of trees through all the years of image acquisition. In the case of detection of both trees growing singly and in high canopy closure, their heights showed an increasing trend, analyzing the subsequent acquisition dates. A larger dispersion of values was recorded for trees growing singly.
In the conducted research, a better dense image matching performance was observed, with parameters providing higher point cloud density and low smoothing of the model. Considering algorithms compared in the software used, dense image matching from Trimble Inpho and Agisoft Photoscan worked comparably for modern images when parameters were properly set. When archived data were used, better results were achieved for Trimble Inpho image matching algorithms. Worse image quality, lower image resolution and unfavorable period of data collection required an image matching strategy that involved steps based on both feature- and area-based approaches.
As a result of the analysis of the use of the generated elevation model, it was found that for archival photographs of good radiometric quality and a scale of not less than 1:13,000 or maximum GSD of 25 cm, the optimal cut-off threshold for determining the succession extent of trees and shrubs was 1.25–1.75 m. With this threshold, the error of determining the extent of trees and shrubs can be estimated at 10 to 30%, depending on the technical parameters of the images and the date of their acquisition. An important element in defining the cut-off threshold is the phenological development of vegetation during aerial images collection.
For spring images, the cut-off may be lower (even 1 m), because herbaceous vegetation is not yet developed and is of relatively small size. The spring term is particularly recommended for determining the range of succession of coniferous trees and shrubs. When choosing the threshold value, both the acquisition date of archival photos and the specificity of land cover of the analyzed area should be taken into account, and in particular, the possibility of occurrence of high herbaceous plants, which may hinder the process of obtaining the correct extent of trees and shrubs.
The results of the present study have served as recommendations in the development of methodology for analyzing the dynamics of trees and shrubs succession process using archival aerial photographs, implemented by the project leader, MGGP Aero company.