Visual Exposure of Rock Outcrops in the Context of a Forest Disease Outbreak Simulation Based on a Canopy Height Model and Spectral Information Acquired by an Unmanned Aerial Vehicle

This research was focused on the study of visual exposure evolution in the locality of the Drátenická skála nature monument (in the Czech Republic) and the surrounding forest complex in terms of history and through modelling for further possible stand development. The local forests underwent conversion from a natural fir-beech composition to an intensive spruce monoculture with few insect pests or windbreak events to an actual bark beetle infestation. Historic maps, landscape paintings, photographs, and orthophotos served as the basic materials for the illustration of the past situation. Further development was modelled using canopy height models and spectral properties captured by unmanned aerial vehicles (UAVs). As an example, the possible situation of total mortality among coniferous spruce trees after a bark beetle outbreak was modelled. Other options and a practical use of such preprocessed data are, for example, a model for opening and transforming the stands around the rock as one of the ongoing outcrop management trends in the protected landscape area (PLA) of Žďárské vrchy.


Introduction
The authors would like to dedicate this article to the memory of Associate Professor Jan Lacina, who passed away on 21 March 2020. He was a Czech forester with the soul of poet and was a geobioecologist; lover of nature, birds, and landscape painting; and a tireless advocate of forest science. Many thanks to him for providing the inspiration for this article and helping with the search for paintings by artists from Žďárské vrchy-a landscape he liked a lot.
Visual exposure is a habitat characteristic that might not seem important at first. However, vegetation cover in localities with protruding rock outcrops could be the key factor in the life conditions for animal and plant species that occur specifically in this environment. This study targeted the protection of rock outcrops, which are considered to be specific microecosystems and are biodiversity hotspots facing the threats and impacts of ongoing climate change [1,2]. They are refuges for flora and fauna from the surrounding adversely human-affected landscape and have unique characteristics for life [3]. These areas remove unfavourable edaphic, nutrient, temperature, and other climatic conditions, and along with providing topographic heterogeneity, they also serve as suitable habitats for endemic species that require a specific environment [4].

Aim of the Study
In this study, we primarily aimed to describe the vegetation development around Drátenická skála, as it is one of the most attractive localities that has been influenced by humans or by natural processes in the past. We also investigated its possible future state, which could be decisive for the suitability for many organisms. Using geographical information systems (GIS) analyses, we modelled the possible visual exposure of the stands after a controlled intervention or a natural outbreak [5] and compared the situation over the last centuries through the use of:  landscape painting and historical maps,  historical and current photographs and orthophoto, and  unmanned aerial vehicle (UAV) remote-sensing data and digital terrain and surface models.

Landscape Painting
Painting is the oldest method of depicting reality. It began in prehistory, when our ancestors painted on cave walls and other artefacts. The ochre lines on the stone from the Blombos Cave (South Africa) are probably the oldest known human drawing. They have been dated to Middle Stone Age and determined to be 73,000 years old [6]. Other thousand-year-old paintings of animals on cave walls are known from France, Germany, and Spain [7]. We can speculate about their purpose, but there are several explanations. They could have served as messages about the presence of the depicted animals as the main source of a meal [8], as the oldest schemes or maps, or even for ritualistic purposes [9]. In connection to the evolution of the human brain, these drawings could also have been a form of personal expression. Over time, painting became a form of art with many directions and styles. Landscape painting seems to be a unique style that developed in the 17th century in the Netherlands and reached its prime during the 19th and 20th centuries. This then led to the development of several artistic styles, such as romanticism, realism, modernism, and impressionism, which were used by painters such as Jacob Ruisdael, El Greco, Carlos de Haes, Frederic Edwin Church, Pablo Picasso, Paul Cézane, Claude Monet, and Vincent van Gogh. In Czech lands, the history of landscape painting began roughly in the 19th century, as our artists were influenced by French Barbizon and other European schools. The most significant Czech landscape painters included Maxmilián Haushofer, Antonín Mánes, Julius Mařák, Otakar Nejedlý, František Kaván, Antonín Slavíček, Oldřich Blažíček, and Josef Jambor [10].

Remote Sensing
The phenomena of modelling the Earth´s terrain shapes began with the development of geographical information systems (GIS) and remote sensing (RS). One of these terrain-modelling technologies is light detection and ranging (LiDAR), which is built on the principle of a laser beam reflection from the Earth's surface. Laser scanners are active systems that emit light pulses and record their energy. These readings are carried out by a specially constructed airplane. Other components of this system are the length measurement unit, the control and recording unit, the IMU (inertial measurement unit)/global navigation satellite systems (GNSS) unit, and the imaging unit [11]. By calculating the measured values of the distance and time between the emission and reception of the laser beam and determining the reflection place, the X, Y, and Z coordinates of the points are exactly determined, and a point cloud is generated. Using several interpolating methods, the models of surface and terrain can be calculated [12]. Currently, in the Czech Republic, we are able to work with the first-generation digital surface model (DSM 1G), which has an accuracy of 0.7 m for vegetation and 0.4 m for buildings [13,14], and the fifth-generation digital terrain model (DTM 5G), which has an accuracy of approximately 0.18 m in open terrain and 0.3 m in terrain covered by vegetation [13,15].
Currently, there is a growing interest in the use of unmanned aerial systems and the structurefrom-motion (SfM) algorithm in forestry. This technology has been most studied for use in forest inventories [16]. Typically, the created canopy height models and spectral information are used to estimate forest variables [17]. A further application field is the monitoring and assessment of forest health conditions. For this purpose, a combination of 3D models with 2D spectral properties are often used and obtained using multispectral sensors and are also used in the previous cases [18,19]. Complex studies dealing with the UAV-based estimation of forest variables such as the dominant height, Lorey's mean height, stem density, basal area, and timber volume have yielded the work of Puliti et al. [20] and Tuominen et al. [21]. In temperate forests, Mikita et al. [22] determined the biophysical forest parameters using a combination of UAV and close-range photogrammetry. Methods for UAV-based forest inventory can be divided into area-based approaches (ABA) [23] and individual tree crown (ITC) approaches [24]. ITC enables the detection of trees with an accuracy of 25% to 90% [25,26], to classify them by tree species with an accuracy of up to 95% [27] and to measure tree heights with an RMSE ranging from 0.5-2.84 m [25,28].
Currently, the role of remote sensing in forestry is growing, especially for the purpose of assessing the forest health status. Remote sensing enables the rapid and objective evaluation of a wide area of forests. Using multi-or hyperspectral sensors, it is possible to determine the content of chlorophyll, water, and leaf pigments [29]. The acquisition of discrete spectral data using UAVs is possible, because many lightweight sensors have been developed in recent years. However, compared to RGB cameras, multispectral cameras are characterised by lower resolutions [16]. Unlike agriculture, where spectral information is already routinely used, there are only a few studies in forestry dealing with the use of spectral data acquired by UAVs. Näsi et al. [30] used a combination of the ITC approach and hyperspectral data as a cost-effective alternative to field-based monitoring of bark beetle infestation. According to Minařík and Langhammer [31], the red edge and NIR part of the electromagnetic spectrum are the most useful for stress monitoring. The normalized difference vegetation index, especially in time series and in conjunction with the random forest classifier, may be one of the best predictors for mapping physiological stress [32]. The main reason for using UAVs is above all the high flexibility, which is especially advantageous today, when forests are changing rapidly due to climate change and related diseases. The minor reason for creating the UAV/SfMbased canopy height model was the high point density compared to the general LiDAR product, DTM 5G. It should be noted that the minimum LiDAR point density required for the creation of canopy height models is approximately 5 pts per square metre, and the average density of available LiDAR data for this locality is 1 point per 10 square metres [33].
In regard to specific rock outcrop research, there are many studies aimed at rock outcrop modelling using different remote sensing technologies or terrestrial laser scanning; in some cases, the studies are also connected with structural geological surveys [34][35][36][37][38]. With increasing technologies, outcrop modelling is not limited only to the Earth but is also performed within the solid-state planets or moons of the solar system [39]. We can also mention the previous research of the authors' team, which compared different 3D models of boulders [40]. In the locality of Žďárské vrchy, few studies dealing with outcrop management have been performed; one of the few from recent years was a study of Malinská skála and identification of its outcrops [5].

Protected Landscape Area Žďárské vrchy
The highland of Žďárské vrchy is situated in the Central Czech Republic, Vysočina and Pardubice regions (Figure 1) on the "roof of Europe", where the main divide of the European Elbe-Danube rivers runs through. The area is well-known for its specific landscape characteristics, which comprise typical geomorphological forms, rock formations, deep forests, scattered vegetation, stone pastures, flowery meadows, river systems, and preserved natural habitats completed by elements of folk architecture and harmonious coexistence of humans. All these phenomena were and are endangered by human expansion and misuse of the entrusted nature wealth, which started with Middle Ages colonisation-the first written indication of a settlement is from 1420 [41,42]. The protected landscape area (PLA) was proclaimed in 1970 by an announcement of the Ministry of Culture of the Czechoslovak Socialistic Republic and covers over 70,000 ha [43]. The subject of protection was defined by this announcement, which says [44]: "The mission of the locality is the protection of the local landscape, its looks and typical features to create a harmonious environment. Within these typical landscape features, the specific surface formations include the watercourses and water areas, design and use of forest and agricultural land funds, vegetation cover and wildlife, housing estate structure, urban composition of housing estate, local folk architecture and considerable or dominant forms of production." The following purposes are considered to be the most important and seem to be the most related to our topic: preservation and manging the natural forest ecosystems, increasing the ecological stability, and conserving typical abiotic natural features and landscape characteristics.

Geology and Geomorphology
The PLA consists of several geological units of the Czech Massif formed during the Varisk folding in the Palaeozoic (Devonian and Carboniferous). In the Central and Northwestern parts, we can find the largest one, Svratka Crystalinicum, with migmatites, double-daize orthogneisses, and mica schists. The most significant features are outcrops of coarse grain augen orthogneisses with amphibolites and scarns. Minor features were formed in Strazec Mouldanubicum with paragneisses, amphibolites, crystalline limestone, and serpentinites, and then, Polička Crystalinicum with gneisses, Iron Mountains pluton with granites and porphyries, Hlinsko zone with phyllites, Ransko massif with eruptive rocks and iron ores, and finally, the area of the Czech Cretaceous Formation with lime sandstones, claystones, and marlites. During the next geological period-the Pleistocene-the outcrops were influenced by cryogenic weathering, creating karstic and stone sediments and fallow fields. Alluvial river terraces and moors were formed during the Holocene [45][46][47].
The local landscape is well-known for its characteristic large ridges and widespread river valleys between them. This macro-relief has cryogenic microstructures, e.g., rock outcrops, tors, and frost cliffs surrounded by fallow fields and cryoplanation terraces. All these geomorphological phenomena are a significant part of the local geoheritage based on the nature of geodiversity. The importance of geoheritage has been rising recently as we have increasingly aware of the need for geosite (or geomorphosite) protection [48,49].

Flora and Fauna
The area was covered by dense primeval forest until medieval colonisation in the 13th century. Fragments of that natural forest have been preserved until the present day and are usually located around ridges with outcrops; all of them are the objects of conservation [42]. Several botanical surveys have been carried out, and they mention, for example, the following significant botanical species from many plant communities (e.g., [50]). The typical vegetation compositions of these stands are acidophilic spruce beechwoods (Luzulo-Fagion), with beech (Fagus sylvatica L.), spruce (Picea abies (L.) H. Karst), and fir (Abies alba Mill.) as the dominant wood species [51]. In some places, residual moors and waterline spruce woods are extant with other tree species-pines (Pinus sylvestris L. and P. rotundata L.). Representatives of other less dispersed stands and communities are preserved in several protected landscape sites, as is the herbal cover of rare plant species (Leucojum vernum L., Aconitum variegatum L., Vaccinium uliginosum L., Oxycoccus palustris Pers., Eryophorum vaginatum L., and Sphagnum sp.). The majority of these species were replaced by Norway spruce production forest. Additionally, on nonforest protected sites, such as standing and running water ecosystems, reeds, moor meadows, or heaths [34], it is possible to view a wide range of rare vascular plant species, e.g., Caricion sp., Carex sp., Sphagnum sp., Trichoporum alpinum L., Drosera rotundifolia L., Dactylorhiza majalis (Rchb.) P. F. Hunt et Summerh., and Nymphaea candida J. Presl [43].

Drátenická Skála
Drátenická skála nature monument is situated in the Blatiny cadastral area, which has an altitude of 776 m and is the highest outcrop peak in the area. The rock belongs to Svratka Crystalinicum double-daize migmatites and orthogneisses. This 200-m rock ridge is considered to be a relict frost cliff with a few isolated tors. The main rock block consists of several tors (towers) with passes between them ( Figure 2). Sokolí věž is the highest (35 m), and the others include Orlí věž, Otoman, Zbojnická věž, and Indiánská věž. Krakonošova zahrádka is a ravine with rocks bedded alongside the foliation layers. Around the blocks, fallow rock fields and boulder streams spread. Additionally, many microrelief features can be found, such as abri, rock tunnels, rock windows, or honeycombs. Furthermore, some geological processes can be observed, such as well-developed joint systems, cryogenic weathering, and rock collapse. The locality is also strongly influenced by humans because of intensive rock climbing and visits by thousands of tourists during the year [54,55].

Baseline of the Study
Our study in the area of Drátenická skála was inspired by the approach of Jan Lacina, the Czech landscape ecologist, who conducted several studies of changes in the landscape through monitoring and assessment based on comparing the contemporary real visual situation of the localities with landscape paintings by Czech artists. His study was aimed, for example, at evaluating changes in the ecological conditions in the chosen localities in the Vysočina region [10]. Another study focused on the changes around contemporary water dam localities [59]. These ideas were the initial material for our concept and were complemented by computer modelling using GIS.
Another idea that was included in this study is linked with unique outcrop management plans in protected areas that consist of stand-thinning and rock exposure through the use of new visual axes [60,61]. This very slow and deliberate procedure should lead to previous biotope conditions before spruce afforestation, as well as to higher species diversity (e.g., Peregrine Falcon, Falco peregrinus L., nesting is expected), ecological stability, strengthening forest functions, and stabilising microclimatic functions important for rare species occurrence (e.g., [62][63][64]). This process is anchored in management plans of Malínská skála, Bílá skála, Černá skála, and Devět skal reservations [64][65][66][67]. However, these days, not only does the Czech Republic have to face the largest bark beetle infestation in recent history, the PLA Žďárské vrchy probably will not stay passed, and we must count on other endangering events. Thus, that intended management would be influenced by the outbreak to some extent, as the outcrops are not surrounded by the original stands of European beech (F. sylvatica), silver fir (A. alba), and sycamore maple (Acer pseudoplatanus L.) but, largely, by Norway spruce (P. abies). These stands are truly threatened by spruce bark beetle (Ips typhographus L.) and may go extinct in a few years.

Methodological Approach
With an effort to reach our aims described in Section 1.1. (Figure 3), the approach started by looking for the most appropriate outcrop locality within Žďárské vrchy; however, this choice depended on the required visual attractiveness, available landscape paintings, photographs, and other materials. After this exploration, Drátenická skála was selected as the most suitable outcrop area because of its high visual attractiveness and potential special management, as characterised above [68]. Due to the stand-thinning process and new visual axes that have not yet been designed, we focused on modelling the spruce bark beetle outbreak at the locality.
Drátenická skála was artistically rendered by several painters from the Vysočina region. The following authors and their work were included in the research (  The extent of the forest around the rock formation can also be estimated from the old maps. In the Czech lands, several historical maps are available, such as products from the 1st, 2nd, and 3rd military mappings, which served as the basic material for modern cartography in the 20th century ( Figure 5). We experienced considerable difficulty in locating relevant historical photography; nevertheless, the shot of haymaking in Milovy (1941) with Drátenická skála in the background was discovered (Figure 4). Additionally, historical airborne orthophotos offered quite different but also highly rewarding landscape views. Currently, this rock formation represents one of the most popular objects in the region for photographers and artists ( Figure 6).  Our remote-sensing analyses consisted of imaging the locality with a UAV. We used a fixedwing UAV: the senseFly eBee Plus, which was equipped with a Parrot Sequoia+ multispectral camera (both products manufacturer: senseFly Parrot Group, Rout de Genève 38, 1033 Cheseaux-sur-Laussane, Switzerland) [74]. This camera has a global shutter and four spectral bands with 1.2 Mpix resolution: green (550 nm ± 40 nm), red (660 nm ± 40 nm), red edge (735 nm ± 10 nm), near infrared (790 nm ± 40 nm), and a separate RGB camera w 16 Mpix resolution and a rolling shutter. The camera is also equipped with a sunshine sensor for automatic radiometric calibration. The separate inertial measurement unit and magnetometer were also included. The UAV flight mission was carried out at a mean height of 159.2 m above the terrain with 60% lateral and 80% longitudinal overlap in the resulting ground sample distance (GSD) of 15 cm/px for spectral bands. The distance between the photos within one flight line was 29 m, the flight line spacing was 77 m, and the single image coverage was 192×144 m. In total, 520 multispectral images were captured. A postprocessed kinematic (PPK) method was used to improve the image position accuracy. Differential corrections in the one-second sampling rate were created using a virtual reference station (VRS) based on the Czech network of GNSS permanent stations (CZEPOS). These differential corrections were imported into the eMotion software to improve GNSS accuracy. The resulting root mean square error (RMSE) of position coordinate error of the image taken by the UAV was 0.018 m for XY and 0.022 m for Z. The RMSE after block bundle adjustment was 0.005 m.
By processing the captured image data using a well-known workflow in Agisoft Metashape, the point clouds (Figure 7), detailed RGB orthophotos and orthomosaics for each spectral band were obtained. Using the SfM point cloud, the digital surface model (DSM) was calculated by the nearest neighbour interpolation method with a raster resolution of 0.1 m. The digital terrain model (DTM) was calculated by the same method at the same resolution by utilising the LiDAR point cloud data from the Czech State Administration of Land Surveying and Cadastre (ČÚZK). Subtracting these two raster models, the vegetation height was calculated, which is presented as the so-called canopy height model (CHM). With the inverse watershed segmentation method [75], individual tree crowns were identified automatically, as were the top and height of each tree. During this process, the neighbourhood circle radius for CHM smoothing by focal statistics was experimentally estimated. The best results were achieved using a 1.7-m radius. The pixels located lower than 3 m above the ground were removed from the model to eliminate low vegetation, which does not influence visual conditions of the locality. These raster data were converted to vector format to prepare them to align with orthophoto data. The CHM, treetop position, and height served to visualise the situation by the 3D model.
To determine coniferous and deciduous trees, two raster images: true orthophoto and normalized difference vegetation index (NDVI) of the forest stands around the rock were used. However, there are other indices that could be created from the Sequoia+ cameras that are useful for tree species determination, e.g., normalized difference red edge index (NDRE), enhanced vegetation index (EVI) or soil-adjusted vegetation index (SAVI), and canopy chlorophyll content index (CCCI).
The orthophoto was taken during the UAV flight and is a combination of red, green, and blue bands. For these analyses, iso cluster unsupervised image classification was conducted in ArcMap 10.7, because RGB orthophoto pixels were distributed into 20 colour classes. Then, compared with RGB orthophoto, these 20 classes were classified to represent the coniferous or deciduous trees based on the fact that each tree type has a specific spectral reflectance. Classes representing canopy shadows were erased. The Zonal Statistics tool served to determine whether the majority pixels within each individual tree crown area were coniferous or deciduous. Through this analysis, the crowns were classified as coniferous or deciduous.
The NDVI is one of the basic vegetation indices that shows how much chlorophyll is contained in a plant. The idea is based on the fact that chlorophyll strongly absorbs visible light (wavelengths 0.4-0.7 µm) during photosynthesis, whereas a specific plant cell arrangement strongly reflects the near-infrared light spectrum (0.7-1.1 µm). By this, it is possible to quantify the photosynthesis capacity of particular pixels expressed by NDVI values ranging between -1 and 1. Generally, we can say that a significantly higher intensity of reflected near-infrared light than absorbed red light represents dense vegetation (forests, draws, etc.), whereas the opposite case represents thin vegetation (grassland, withered, or damaged trees, etc.). We used the red and near-infrared bands obtained by a spectral UAV camera for NDVI calculation according to the following equation [76]: NDVI = (ρNIR-ρred)/(ρNIR+ρred). (1) Then, the NDVI raster was classified into coniferous and deciduous trees in the same way as the orthophoto.
To complete the visual conditions model around the outcrop, visibility analyses of the rock peak were also performed, as it was calculated for contemporary forest area and for the situation without coniferous stands.
To identify the visible rock borders, we counted that, in the rock area, there was minimal altitudinal difference between the DTM and DSM. Consequently, visible rock was extracted during the removal of low vegetation. To eliminate errors, manual control of the orthophoto was performed as the fastest and most suitable choice for gaining rock borders with the highest possible accuracy.

Historic and Artistic Material
In local galleries, three landscape paintings of Drátenická skála dated from 1946, 1951, and 2005; one photograph dated from 1941 ( Figure 4); and a historical orthophoto dated from 1953 and others from the turn of the 20th and 21st centuries ( Figure 5) were found. The situation in the 19th century was estimated from the second and third military mapping (1836-1852 and 1876-1878), and the maps of the first military mapping belonged to the 18th century (1764-1768). Last, the named maps were not rigorously positioned, and the rock formation was not specifically depicted either. It was not possible to determine the area of the forest around the outcrop. The two subsequent military mappings better determined the rock position, as well as the area of the surrounding forest. In particular, the second military mapping showed it as a relatively significant visual landscape feature ( Figure 5).
The stand area in the first half of the 20th century could be estimated with the help of the collected paintings and photographs. After a detailed survey and identification of a presumable author location, it was obvious that forests were formed by young planting of spruce, so the trees were not as high as they are currently, and the rock was more exposed. Additionally, in other historical photography, this trend was also obvious-for example, the outcrop of Devět skal ( Figure  4) was exposed completely, and moth infestation events were documented from the beginning of the 20th century, whether it was nun moth infestations (1920)(1921)(1922) or windbreaks in November 1930. Much more intensive spruce afforestation followed, and these monocultures contemporarily reached the top parts of many local rock formations, some of which are now completely covered by transformed forest and have artificial spatial and species composition. In the example of Drátenická skála, we observed a partial stand loss from the south, where there are shorter broadleaves. This situation is very close to the required result of the PLA management plan. When viewed from the other site, the rock is surrounded mainly by spruce monocultures or mixed stands with spruce dominance. This statement is also supported by both historical and contemporary orthophotos ( Figure 6).

GIS Modelling
By applying the GIS procedures described above, the DTM and DSM were created. The area of visible rock and the vegetation lower than 3 m were erased from the DSM, and the canopy height model (CHM) of the current stand extent was obtained (Figure 8).
By determining the individual tree crowns, assessing their heights, and categorising them as deciduous or coniferous, we obtained a detailed overview of the stand spatial arrangement and structure in the immediate surroundings of the rock formation.
For the tree classification, two types of spectral image-based data were used: suitably-coloured orthophotos and NDVI. However, we must expect some degree of inaccuracy, because some pixel classes represented both tree types at a similar rate, and shades also play an important role in iso cluster unsupervised classification (Figure 8).
For the following analysis of the spatial species composition, the area within a 250-m radius from the centre of the outcrops was evaluated in greater detail. Classification of both the orthophoto and NDVI achieved similar results at the stand level (Table 1 and Table 2), but significant differences at the tree level were discovered (Table 3 and Table 4). The UAV flight occurred in September, when the differences between the coniferous and deciduous leaves were evident (especially when we mainly distinguished spruce and beech with various green hues; the chlorophyll content in deciduous leaves also decreases during autumn). This was crucial to make the most accurate possible classification. It can be assumed that much worse results would be achieved during image capture in different seasons, especially in the case of orthophoto-based classification. Unfortunately, it was not possible to perform NDVI thresholding for further use. It is always necessary to reclassify the images, because the vegetation has a different chlorophyll content in different seasons. The NDVI value could also be influenced by the water content. Therefore, supervised classifications and machine-learning methods are more successful for more complex analyses, such as health assessments [32].
In terms of the basic wood composition, the stand comprised 42% deciduous trees and 58% coniferous trees. The average whole stand height was 22.2 m, with 16.5 m for deciduous trees and 25.8 m for coniferous trees. Individual deciduous crown areas were mainly larger than those of coniferous trees, and the lower tree trunk density in deciduous and mixed stands corresponded with this difference. Stand density was highest on the northern side, where spruce predominance was apparent, whereas, on the southern side, where deciduous and mixed stands prevailed, the dispersion and openness of the forest were obviously higher ( Figure 8). This situation could roughly correspond to the resulting state after the intended management plan measures. Comparing this classification with the forest stand map (for our purposes, this map was simplified into three categories according to the main stand type-coniferous, deciduous, or mixed), we can state that our NDVI-based classification corresponded to this map, because the coniferous stands were 85% coniferous trees, deciduous stands were 81% deciduous trees, and the proportions of the two tree types in mixed stands were approximately 55% vs. 45%, respectively. It follows that, by using the unsupervised classification, a quite accurate differentiation of coniferous and deciduous trees was reached ( Figure 9 and Table 2). The following step consisted of 3D modelling of the present stand appearance and then modelling the situation that could potentially occur as a result of a current bark beetle infestation if the nonoptimistic scenario was fulfilled and all the spruce trees died. It should be considered that not all coniferous species were spruce, as other species (larch, pine, and fir-approximately 4%, altogether) were present [77], and they were also erased from the model during this step. This scenario could be remodelled for the situation when the bark beetle, for example, did not affect the mixed stands. In Figure 10, the visualisation of this change is shown from one of the most attractive points of view in the comparison with the actual photography and landscape in the paintings.  Confusion matrices (Table 3 and Table 4) for classification of the NDVI and orthophotos were computed based on 200 stratified accuracy points. Visual investigation of the orthophotos was used as the ground truth classification. At the tree level, the NDVI-based classification was significantly more successful. However, both classifications had very balanced types of errors. Nevertheless, the error of omission within the classification of the deciduous trees was the highest (deciduous trees were more often mistakenly classified as coniferous) in both classifications.
With these basic data and models, the choice to design various possible situations was then quite simple to apply. In practice, the process could consist of the localisation of trees considered to be felled by global navigation satellite systems during the field survey. These localised trees were related to the canopy height model of the crowns, erased, and, in accordance with the process described above (Section 2.3.), the 3D model of this new planned visual tree composition was designed. By this approach, we will be able to choose the most appropriate measures leading to stand-thinning and conversion.
By using more detailed image classification, it would also be possible to determine various tree species and operate with better differentiation, which could enable us to focus on more specific situations, such as mortality in a particular species, as a consequence of other insect pests, from bacterial, viral, and fungal diseases or due to increased drought sensitivity. Practical use of this approach is suitable mainly for the care of similar specially protected areas that are valuable by their visual attractiveness or in the care of less extensive forest properties. For larger forest complexes, UAV imaging is still too uneconomical; however, it is possible to use freely available spectral satellite data with high (10-30 m) or very high resolutions (0.3-10 m), which are commonly used in current forest remote sensing. Additionally, a visibility model was calculated for both situations: with and without coniferous trees. We analysed from which places the rock peak (776 MASL) was visible when it was covered by contemporary forest stands and from which it would be visible with total conifer mortality. These models were compared within the 335,801-m 2 surrounding area polygon ( Figure 11 and Table 5). The peak was visible from 3.9% of the area (13,216.8 m 2 ), and in the case of a model without conifers that was identified by the NDVI-based classification, it was visible from 14.4% of the area (48,338.5 m 2 ).

Conclusions
In this study, the development of forest cover around Drátenická skála over the last few centuries was depicted and evaluated using historical maps, orthophotos, photography, and landscape paintings by local artists. The GIS analysis consisted of taking data by UAV scanning and then computing each tree location, crown area, and height. The tree type (deciduous or coniferous) determination was realised using orthophoto and NDVI image unsupervised classifications. Higher classification accuracy was demonstrated on the basis of the NDVI, especially at the tree level. These processes led to the creation of forest cover and visual rock formation models that corresponded well to the real situation that was documented by the forest stand map. Additionally, a model of the possible situation after a bark beetle calamity, which would kill all spruce trees, was created. Finally, the visibility analysis was added for both forest stand situations. Another modelling method was discussed in connection with the fact that the spruce monoculture growth is naturally considered to not be the best current method, and the current PLA management relies on the transformation and opening of the stands.
This study can contribute to effective management planning specifically in outcrop localities that are being rebuilt with the aim of protecting microclimatic life conditions for endangered species that are dependent on these specific habitats. This study can also be helpful during any outbreaks of tree species diseases that threatens to kill the stands. Vysočiny Polička for the permission to use the painting reproductions. Our thanks also belong to Lubomír Dajč for providing his aerial photographs.
Special thanks are owed to the organisers of the GIS Ostrava 2020-UAV in Smart City and Smart Region conference. This article was published in the special issue for this conference, and an abstract was published in the proceedings, which can be found at this link: http://gisak.vsb.cz/GIS_Ostrava/GIS_Ova_2020/proceedings/papers/gis20205e33e0d8b1b75.pdf.

Conflicts of Interest:
The authors declare no conflicts of interest.