1. Introduction
Forest fires are very frequent crisis situations, especially in dry or arid landscapes [
1]. The prediction of forest fire occurrence depends on knowledge of the factors that affect the fires and on the technologies that facilitate the monitoring and modeling the spread of the fire. Ganteaume et al. [
2] analyzed the most common human and environmental factors driving forest fire ignition. The primary factors that directly cause forest fires are natural (lightning strikes, seismic and volcanic activity, etc.) or human (carelessness and activities such as arson, slash-and-burn agriculture, fire-fallow cultivation, machinery sparks, discarded glass bottles or cigarette butts, military activity, etc.). The factors that determine fire spread are as follows: forest type and structure (distances between trees,
DBH, canopy height, tree crown density, etc.); meteorological conditions (precipitations, temperature, wind speed, air humidity, cloudiness, soil moisture, etc.); topographic (morphological shapes of terrain, orientation of relief slopes, etc.); geological and pedological (underground structure, soil structure, and terrain surface color); and season and time of day, which determine the amount of available sunlight and temperature, etc. The technologies that facilitate the monitoring and modeling of a forest structure and the spread of forest fires include the sensor types for vegetation data collection and forest structure determination and technologies for monitoring and modeling the spread of fire. Blair, Rabine, and Hofton [
3] described the Laser Vegetation Imaging Sensor (LVIS), which operates at altitudes of up to 10 km aboveground and is capable of producing data for topographic mapping with dm accuracy and vertical height and structure measurements of vegetation. The LVIS instrument is also suitable for subcanopy ground elevation mapping. Lim et al. [
4] described many of the initial studies of the application of LiDAR for forestry focused on verifying through statistical analysis that LiDAR could be used to accurately measure forest attributes. The focus has been on canopy tree heights given the nature of this attribute as a predictor variable for other forest attributes, such canopy density. Hyyppä et al. [
5] analyzed existing algorithms and methods of airborne laser scanning that are used for extraction of the canopy height and individual tree information. Aschoff and Spiecker [
6] described an algorithm for detecting trees in a semiautomatic way. Gobakken and Næsset [
7] analyzed the effects of forest growth on laser-derived canopy metrics. Carson et al. 2004 [
8], Ahlberg et al. [
9], and Su et al. [
10] provided the overview of LiDAR applications in forestry. By comparing these methods based on laser scanning, it can be stated that, at present, an approximate forest structure for modeling the movement of rescue vehicles can be determined. At the metric density of the DSM (CHM), the error in tree positioning can reach values in decimeters, sometimes up to meters, depending on the structure and type of the canopy. At the decimeter density of the DSM points out, it is possible to calculate the tree position errors in centimeters to decimeters. When comparing the possibilities of using LiDAR and aerial optical images, both methods have advantages and disadvantages. Aerial images provide both spatial and image information, but they do not allow, unlike LiDAR, full automation to determine the forest structure.
However, the use of LiDAR and aerial optical images may be problematic in the area of fire because of the clouds or smoke generated by the fire. In these situations, radar methods can be used to measure forest parameters. The mapping of forest units by radar is described, for example, by Martoni et al. in [
11]. Kugler et al. [
12] compared the LiDAR and radar methods for determining the heights of the forest in three areas: boreal, temperate, and tropical. The correlations achieved confirm the possibility of combining the use of both forest mapping methods. Additionally, Cazcarra-Bes et al. [
13] described the possibilities of the horizontal and vertical forest structure mapping from radar using data obtained by synthetic aperture radar tomography. The use of radar methods is, however, limited in terms of the accuracy of the determination of the characteristics of individual trees. Furthermore, Landsat or Sentinel 2 global satellite data can be used to monitor forests before and during a fire. Sentinel 2 with a multispectral instrument (MSI) with 13 spectral channels in the visible/near-infrared (VNIR) and shortwave infrared spectral range (SWIR) and three bands for vegetation mapping can provide the crisis management with actual data in the shortest possible time, especially during a forest fire. The data accuracy (about 20 m) does not allow a more accurate mapping of the internal forest structure—see, e.g., Puletti et al. [
14].
Technologies for monitoring and modeling the spread of a fire are divided into stages: prediction, during the fire, and post-fire [
15]. Milz and Rymdteknik [
16] described the technologies of detection and the spread of the forest fires by using satellite-borne remote sensing techniques. However, the technologies of fire monitoring and distribution are limited by the availability of up-to-date data from satellites, planes, UAVs, or terrestrial observations. We also need to know the prediction of the spread of a fire to deploy rescue vehicles. Koo et al. [
17] described possible solutions using a physical model for the forest fire spread rate. This model successfully evaluated wind and slope effects of a fire on forest vegetation.
The above-mentioned factors and technologies are very important for the teams (fire brigades, military units, health services, and police) that are deployed to rescue people and reduce the damage during forest fires. Remote sensing support is very important for rescue units when they are moving across vegetation before and during a fire and, also, for the decision to deploy aircraft. We can use LiDAR and aerial image data to create a navigation analysis for rescue (fire brigade or military) vehicles—see also [
18,
19,
20]. These data can be supplemented by active fire scenes using infrared sensors or aerial or UAV images. Among the most effective data sources for Cross-Country Movement (CCM) navigation and optimal pathfinding across a forest are LiDAR data and the products of its analysis. A prerequisite for the success of this analysis is an up-to-date picture of vegetation data obtained by laser scanning. This precondition is especially crucial for forest stands, where data become quickly outdated due to vegetation growth.
The primary focus of this article was the LiDAR data update for the forest stand structure, a simulation of the creation of a forest structure with the subsequent creation of a model for navigating the movement of a rescue vehicle between trees as obstacles in the terrain. The reason for designing the method of detecting the current forest structure was that LiDAR data in the Czech Republic is gradually becoming obsolete as a result of tree growth. The following procedure was chosen: (1) Selection of the most common type of forest stands in the territory of the Czech Republic with the predominant spruce tree (Picea abies). Obtaining LiDAR data characterizing DSM with a density of 1 × 1 m. (2) Obtaining inventory data on the growth of spruce trees from MENDEL University, Brno. (3) Detection of DSM accuracy by geodetic and photogrammetric method. (4) Corrections of tree heights due to DSM density and tree growth. (5) Creation of forest structure by random distribution. (6) Selecting a simulated area where a fire could occur (older, drier forest). (7) Calculating the simulated shortest route for a particular vehicle (outside the area of the fire). The research results presented in this paper represent a new methodology of updating a digital surface model (DSM) or canopy height model (CHM) using the equations of tree growth and vegetation growth curves. DSM and Digital Elevation Model (DEM) data evaluated for forestry passability were scanned by LiDAR in 2013. CHM data and field measurements were used for determining the approximate forest structure (tree height,
DBH, and stem spacing between trees). The described methods were tested on a spruce forest stand composed only from one type of tree—Sitka spruce (Picea abies), situated approximately 300 m south of the village of Brno-Utechov (see
Figure 1), where the heights of trees in the Krtiny Training Forest Enterprise (TFE) area were detected and measured. This spruce forest was chosen because of the availability of a series of aerial photos and LiDAR data. Additionally, the forest is highly representative, as it contains the tree species most commonly found in many Central European countries. The current age of this forest is about 30 years.
The approximate forest structure with the calculated tree density (stem spacing) was used for modeling fire brigade vehicle maneuvers between the trees. DBH data were used in cases where it was easier for the vehicle to override the trees rather than maneuver between them. DBH data were also used for the calculation of the distances between trees. Due to the availability of a DSM with a density of 1 × 1 m, it was impossible to precisely determine the locations of individual trees, so a random simulated forest structure was created based on the number of trees per hectare.
The article describes the methodology of calculation of the rescue vehicle movement using simulated areas of a burning forest. In the case of a real fire, we can use the above-mentioned LiDAR data (DSM data) or current data from different sensors. The type of sensors and the accuracy of the data obtained will have a significant impact on the terrain analysis and search algorithm for optimal rescue vehicle routes. Only some scattered and low-resolution data of the fire can help. If the optimal route for a special vehicle (fire-resistant rescue vehicle, tank, etc.) that will move through a burning forest is to be calculated, we will need detailed vegetation and elevation data with a meter or decimeter resolution for reconnaissance of fallen trees, boulders, etc.
3. Results
The result of the creation of a forest structure from data obtained from the original forest (see
Figure 8) by generating random tree positions is shown in
Figure 9. The size of the analyzed forest area was 140 × 80 m (11,200 m
2). The length of the vehicle’s passage was 212 m. The direct path between the starting point and the target is shown by a black line. This path is generally unrealistic due to the tree stem obstacles (displayed as green points). All possible paths (blue closed Voronoi polygons) that match the tree distances and vehicle parameters have been computerized using Dijkstra algorithm and displayed in
Figure 8 using our own software tools. Unfinished Voronoi polygons (ending between trees) are nonbinding paths where the width of the vehicle does not allow passage between trees. We can choose any of these blue passable routes, but only one will be the shortest (fastest)—the red highlighted route. This route traverses around (between) the burning forest polygons. The simulation of the polygons displaying the fire areas was done completely at random by adding the points around the impassable zones (orange areas). These areas can be complemented e.g., by satellite images or aerial photos. If we wanted to avoid these risky places, we would have to create a security zone around the burning polygons—so-called buffers—using GIS tools.
There are also displayed the routes inside the areas of fires (blue lines inside the orange polygons)—see
Figure 9. These routes can be used later when the fires end, but they are primarily not included into the calculation of the shortest route. For some types of vehicles, such as tanks, we can also choose the route through the burning area and calculate the route segments inside the orange polygons. The influence of other elements of the terrain (slope gradient, soil properties, terrain surface roughness, forest paths, etc.) are not calculated.
The above-mentioned result of seeking an optimal forest path partially affected by a fire may be modified in case when the forest structure is regular triangular or rectangular. There are displayed all the possible routes and the shortest route—red line in the triangular forest structure in
Figure 10.
The methodology for finding the optimal rescue vehicle path was based on verifying the input data. Determining a forest structure from DSM data can be very unreliable, especially when LiDAR data are obsolete. Therefore, the use of growth curves of trees and derived vegetation parameters were used. These parameters were determined on a relatively small area. Verification of the vegetation parameters by photogrammetric and geodetic pathways lasted several months. The results of the presented model are valid for the designated coniferous forest. The general application of the optimal vehicle route determination will depend on the type of trees (coniferous, deciduous, and mixed). The author assumes that the presented model for finding the optimal route of a vehicle will be better utilized with the development of mapping methods aimed at determining the exact coordinates of the trees.
4. Discussion
The described methodology for determining the possibility of moving the fire brigade vehicles in forest vegetation can be used if tree position data or forest structure (generated from photogrammetric data or from LiDAR data) are available. In both cases, the same algorithm can be used to find the optimum forest path. In case we have more precise tree coordinates (from terrestrial or remote sensing sensors), the calculated route of the vehicle will be more reliable. Although the methods of directly determining the exact tree position by remote sensing data are constantly developing, the forest structure is often determined using DSM (CHM) methods. This is due to the financial cost of the high density of LiDAR data, as well as the personnel demand for data acquisition using photogrammetric methods. The quality of the photogrammetric evaluation depends on the scale of the images and the evaluator’s experience. The main problem is to target the marker at the tree’s top point, which is above the tree trunk. The accuracy of tree position evaluation is higher for coniferous trees than for leafy vegetation. The disadvantage of the photogrammetric method is the lower performance of manual evaluation compared to the possibility of automated evaluation of LiDAR data. LiDAR methods are faster than photogrammetric methods, and they allow a more efficient assessment of the forest structure and determination of the possibilities of vegetation passability without a manual evaluation. LiDAR methods can also be better combined with other remote sensing data sources (infrared, multispectral, radar, etc.). For example, an infrared spectrum can be used to map environmental and fire temperature characteristics, and at night, multispectral imagery can be used to classify species, and radar data is appropriate for mapping a burning forest covered with smoke or clouds. However, these methods have a disadvantage when scanning vegetation with a low density of DSM elevation points (smaller than 1 m × 1 m) and in the case of DTM data absence. On the other hand, the repeated photogrammetric evaluation of the representative forest stands and the data from DSM could bring about a new approach for forest growth analysis and could be a sufficient method for DSM updating in different growth conditions of forest stands.
The results of this experiment showed that this method is fully applicable for the DSM generated from LiDAR data. The method can be appropriately implemented as a relatively inexpensive updating tool for GIS technology between two laser scanning campaigns of a territory. This method can also be refined using the growth curves of individual types of trees. Forest growth characteristics are very important due to the age of LiDAR DSM data. The results of photogrammetric measurements from aerial images taken at consecutive time intervals and statistical calculations show that the growth curves of the trees are initially steeper, but vegetation growth later slows. It is also necessary to investigate the relationships between the natural environment factors and specific canopy growth. The above-mentioned DSM updating method could be used for many applications, e.g., in forestry, military, etc.—see, e.g., [
29,
30,
31,
32,
33]. It should be noted that the tree height correction values decrease with the increasing density of the LiDAR data. At the DSM density 1 m × 1 m, the average correction is approx. +6 m. At the DSM points density of 1 dm × 1 dm, it is possible to estimate the average height corrections of spruce trees in decimeters, depending on the age of the vegetation. Height corrections of the DSM can significantly affect the computationally generated forest structure and, hence, the vehicle motion models. The resulting model of forest crossing by a vehicle will depend, to a large extent, on the quality of the forest structure data. This study focused on a spruce forest—the predominant tree species in Central Europe. In general, it can be said that the species of vegetation may be variable in different forest groups. From this point of view, the study presented in this article can be considered as partially applicable. Using LiDAR/DSM data, the determination of the deciduous forest structure and positioning of the tree trunks will be more difficult, especially due to the crown surface diversity. From this point of view, it can be assumed that the model of finding the optimal vehicle path through the deciduous forest will be less reliable. The success of these models will largely depend on the resolution, coverage, and actuality of LiDAR data, as well as on the accuracy of the forest fire localization data. It should be noted that the use of this methodology in practice has a number of limitations resulting from data that cannot include all objects in the forest, such as lying trees, stones, low tree branches, etc. [
34].
Tree branches can be an important obstacle to the movement of rescue vehicles. It mainly concerns young forests or deciduous forests, where branches are thicker and located below the ground. In coniferous stands, the lower branches of older trees are dry and thin and do not represent a major obstacle for heavy wheeled or tracked vehicles. Below is
Table 1, containing measured data of tree branching; the lowest branches were about 1–2 cm thick. The measurement of tree canopy branching was performed only on trees for which resistance tensile forces were measured, not on all the trees in the area.
The problem is how to get the lower branch data. For this purpose, we plan to use LiDAR data with a resolution in cm [
35,
36,
37] and use the last but one reflection for this measurement. Additionally, terrestrial LiDAR could help to solve this problem—we tested it on a small area in March 2022—see
Figure 11 below.
The author assumes that, in the near future, it will be possible to solve the coordinates of trees, their DBH, and the characteristics of tree crown branches.
The spread of fire, depending on a number of factors, can be very variable, and actual data from burning areas will not always be available. Additionally, visibility can be significantly affected by smoke and the daytime. It should be noted that the calculation of the optimal vehicle route was based on the width of the vehicle. The reliability of route determination also depends on other vehicle parameters, such as vehicle length and height and minimum turning radius. This model did not even include a case where the vehicle would go back (e.g., in case of a spreading fire). For the more accurate calculations of a vehicle route, the impact of the side slope should be also considered. Notable was also the driver’s ability to overcome difficult terrain and to maneuver between trees in crisis situations.
5. Conclusions
The primary aim of this article was to introduce the theoretical aspects, methods, and results of modeling the possibilities of firefighting rescue vehicle mobility in forest areas during fires using remote sensing data. The main result of the presented research is the methodology of the forest structure creation from DSM data, updating due to the growing vegetation parameters, as well as the proposal of the methodology of finding the optimum path of the vehicle to cross the forest, which is considerably more difficult than navigation on the roads.
Although the presented methods are approximate and their applications depend on a number of other factors, the author of the article believes that the presented methods and research results will be applicable in relation to the severity of damage caused by fire. The author also expects further developments of the vehicle navigation methodology in forest regions and the calculations of other factors influencing the search for optimal rescue vehicle routes (low vegetation, lower branches of trees, inclination of slopes, soil influence, terrain surface, etc.). It will also be important to develop the theory and modeling of fire spread in forest areas and to link these models to rescue vehicle navigation.