A Method of Planning Disaster Emergency Rescue Paths in Road-Free Environment

The occurrence of natural disasters has caused great loss of lives and properties that affected many people. To improve the efficiency of the emergency rescue efforts, the problem of road planning for emergency workers, vehicles, and other mobile objects should be overcome. Due to the characteristics of the emergency planning for rescue paths, timely, effective, and rapid efforts should be taken into effect based on geographical data in the road-free environment. This paper analyzes trafficability, researches both access ability and traffic efficiency, and uses traffic efficiency concerning the road-free environment as an influential weight and timeliness to improve the A∗ algorithm, which realizes the rapid planning of the mobile objects' emergency rescue paths in the road-free environment, providing decision-making services for the rescue path planning of disaster emergencies in the environment of no road network. Finally, the experiments were carried out using available data for Shuozhou City, Shanxi Province. The results show that the research can provide a reference for the rescue path planning of disaster emergencies and has contributed to theoretical and applicable aspects of the research.


Introduction
Typhoons, earthquakes, floods, and other natural disasters could occur, causing heavy losses to the lives and properties of many people [1]. Available technical means have not yet been able to effectively prevent and hinder the occurrence of these natural disasters; only by taking high-efficient disposal decisions and emergency response can efforts effectively mitigate losses caused by natural disasters [2]. To improve the efficiency of the emergency rescue and ensure the safety of the people affected by the disaster and distribute the supplies, it is first necessary to solve the path planning problem of disaster relief personnel, vehicles, and other mobile objects [3].
Planning researches of conventional emergency rescue paths [4][5][6][7] are mostly based on the identified road network, rely mainly on experience, and conduct qualitative or semiquantitative analysis based on time, cost, risk, and other aspects, which do not emphasize the emergency rescue characteristics such as the weak economy and strong timeliness. Moreover, earthquakes, mudslides, and other destructive natural disasters may damage the original road network and cause traffic disruption [8], especially plateau mountain area with fewer roads than plain areas, which has fewer alternatives and is more sensitive to natural disasters [9]. While the road network is destroyed or does not exist, decision-makers are required to fully understand the present geographic environment by analyzing the topographical data. Assessing accurately and quickly the prevailing situation and providing reasonable advice for rescue path planning of major disaster emergencies is so crucial that appropriate action can be taken on time [10]. erefore, planning disaster emergency rescue path based on topographical data is particularly important in a road-free environment.
Nowadays, researchers have carried out a series of exploratory researches on the pan planning problem under the road-free environment. e content of the research mainly includes environmental trafficability analysis and path search algorithm. Environmental trafficability analysis is the first prerequisite for implementing path planning. Liu et al. used a digital elevation model (DEM) to generate various curves to approach the earth's surface through certain rules and to construct mathematical models of 3D maps to describe off-road environments [11]. Yuan et al. used a digital elevation model to simulate complex geomorphologic shapes and established a mathematical model to study the effect of slope on the trafficability and efficiency of vehicle off-road mobility [12]. Wu et al. used the area dominant method to rasterize the surface properties, established the roughness evaluation index of the surface properties, introduced the window movement method to calculate the slope of the terrain in advance, and completed the trafficability analysis by establishing a taboo table [13]. Han et al. established a slope and soil taboo table according to the theoretical analysis and simulation experiments to determine passing characteristics of the terrain-vehicle relationship [14]. However, there are still some shortcomings in the current research: (1) e current research mainly considers the influence of local factors on the traffic situation of relief objects in charge of disaster, and the planned path is difficult to meet the actual needs, so the trafficability analysis should comprehensively consider the combined influence of multiple factors. (2) e traffic capacity of various maneuvering objects differs under the same road-free environment. For example, wheeled vehicles cannot pass through the swamp terrain, while tracked vehicles can pass at low speed. erefore, the differences in the attribute information of maneuvering objects should be considered when studying the path planning of maneuvering objects, but the available path planning research does not distinguish this, which cannot give full play to the performances of maneuvering objects. In the literature, several search algorithms and methods proposed to conduct scheduled relief processes under various disaster scenarios have been investigated [15][16][17][18][19][20][21].
A * path search algorithm is based on the established environmental model, searching between the starting point and the target point to obtain an optimal path. For the path planning problem under a road-free environment, the main algorithms currently include genetic algorithm, ant colony algorithm, artificial potential energy algorithm, and A * algorithm. e A * algorithm has the advantages of flexibility and a high success rate among them, which is more suitable for global path planning under a road-free environment. Some experts and scholars have currently put forward some improvements to the A * algorithm. Geo et al. [22] put forward a path planning navigation strategy based on the two-way smoothing A * algorithm, which improves the real-time of the algorithm by optimizing the cost function of the A * algorithm and turning its overall search direction to a twoway street. Liu et al. [23,24] proposed the adaptive controller designs and adaptive fuzzy synchronization of uncertain fractional-order chaotic and nonlinear systems. Jan et al. [25] studied cyber securities by using fuzzy modeling techniques. Xian et al. [26] optimized the A * algorithm by introducing three new methods of "guideline," "key point list," and "bidirectional search" to solve the storage space occupation problem of the A * algorithm [26]. Ma et al. used the key point selection strategy for quadratic programming to effectively remove multiple inflection points and redundant nodes, to improve the planning efficiency [27]. However, with the increase of spatial dimension, the computational complexity of the A * algorithm increases exponentially, and it will eventually face dimension disaster, showing poor scalability, which makes it not perform longdistance emergency rescue path planning tasks in the roadfree environment. Based on the characteristics of strong timeliness and the existence of no road network in emergency rescue path planning, this paper carries out research in the context of the shortcomings of path planning in a roadfree environment.
Firstly, the traffic analysis of the environment without a road network is carried out. e process includes the use of terrain and geomorphologic data to complete environmental modeling. Secondly, environmental trafficability and traffic efficiency are evaluated by formulating pass evaluation rules after jointly considering the interaction between environment and maneuvering objects and their impact on traffic. Subsequently, the A * algorithm is improved with the traffic efficiency as the influence weight and timeliness as the goal, which effectively improves the computational efficiency and realizes the rapid path planning in the road-free environment, providing decision-making services for the path planning of disaster emergency rescue. On this basis, the emergency rescue path planning system in the road-free environment is developed and experimentally verified. e organization of the manuscript is presented as follows: Section 2 deals with the analysis of environmental trafficability based on the grid method to quantify basic terrain and geomorphologic data. A * algorithm and the improved A * algorithm have been presented and proposed, respectively, in Section 3. With the A * algorithm, the heuristic search function is advanced. Section 4 presents data collection and processing, the phases of data analysis, implementation of the proposed method, and results. Section 5 concludes the research.

Analysis of Environmental Trafficability
Environmental trafficability is the basis of path planning under a road-free environment, and it evaluates the trafficability and traffic efficiency of ground maneuvering objects by analyzing the traffic impact of terrain geomorphological data and the performance properties of mobile objects [28]. In this paper, the grid method is used to quantify the basic terrain and geomorphologic data to complete the environmental modeling. Traffic impact is considered based on a relationship between traffic and multiple influencing factors and formulating pass evaluation rules. Afterward, the environmental traffic and traffic efficiency of the maneuvering object are finally calculated.

Modeling of Environment.
Environment modeling is the chief prerequisite to realizing path planning under a road-free environment, which means building the corresponding mathematical model through processing and analyzing environmental information for path search and optimization [29]. e visual graph method, grid method, and so on are generally used to construct environmental models [30]. e grid method is widely used in environmental modeling among them since its simple modeling process, which is conducive to calculation and execution, its strong adaptability to the changeable environment, and its ability to represent different obstacles are the advantages concurrently. e main idea of the grid method is to decompose the environmental information discretely and represent it by grid pixels of a certain scale.
To conduct grid modeling, the following steps are followed: converting the vector data of the environment is planned to raster data, then digitally encoding and merging various types of terrain and geomorphologic information corresponding to the grid pixels are conducted, respectively, and finally, the attribute information of the processed raster data is stored as binary files. is storage method greatly improves the efficiency of reading and using raster information, guaranteeing the rapid generation of path planning schemes. e specific process of environmental modeling is shown in Figure 1.

Conversion Grid Vector Information.
Vector conversion grid refers to the conversion of vector geometry information to grid pixels, providing data support for the implementation of the path planning algorithm in the roadfree environment, realized by calling open-source Geospatial Data Abstraction Library (GDAL).

Encoding of Environmental Attribute.
Environmental attribute encoding refers to encoding the environmental attribute information corresponding to raster data to the corresponding integer number. Types of attribute field data of grid data are divided into floating-point and string. Floating-point fields can be directly assigned to pixels after conversion, such as water depth, vegetation spacing, and slope data, while fields of string type, such as vegetation, geology, and disasters, need to be encoded as appointed integer numbers and then assigned to pixels, for example, the tree in the vegetation type, encoded as a number 1 by the attribute. When data is read, the pixel value is 1, so the vegetation type can be extracted as the tree.

Storage of Information
Overlay. Information overlay storage refers to combining and storing the integer data encoded by environmental attributes for representing various types of influencing factors of traffic conditions, for example, data 10000000, single-bit stored disaster data, tenbit stored geological data, and one hundred-bit stored vegetation type data. For slope data, it occupies three digits and retains one decimal place. For example, 20.5 degrees is stored as 205, of which 005 represents 0.5 degrees. ese three digits occupy thousands, ten thousand, and one hundred thousand of integer data. By analogy, all traffic condition factor data are merged into long integer data and stored in a binary file.

e Rules of Passing Evaluation.
e passing evaluation rule implies that the path planning scheme can be used and practical and thus plays a very important role in path planning under the road-free environment.
is paper comprehensively considers the pass-through relationship between the environment model and maneuvering object model and the environment model and kinematic constraint and thus puts forward the pass evaluation index system, which has an impact on environmental accessibility. Its influencing factors are presented in Table 1. While the environmental model refers to the perception of the environment, including the surface cover and the geometry of the surface, the maneuvering object model refers to its basic physical properties, and the maneuvering kinematic constraint of an object describes its permissible state of motion.
Evaluation factors in the evaluation index system are not independent. Some factors interact with each other. For example, the width of the car body determines the minimum spacing that can pass through the trees, and the chassis height affects its wading capacity so that it determines the maximum average depth of water that can pass through and the climbing ability determines the maximum slope that can pass through the terrain. In addition, geologic and disaster types can also affect the speed of vehicles.
In this paper, the pass evaluation rules have been formulated with combinations of the effects mentioned above, which are presented in Table 2. In the table, the influencing weight W is the weight of all the influencing factors, and the greater the weight is, the greater the impact it has on the passing efficiency of the ground maneuvering object. e velocity coefficient fv is used to measure the extent of each influencing factor's affection to the travel speed of the maneuvering object, which can be adjusted flexibly according to the actual situation, with a value in the range of [0-1]. e velocity coefficient of 1 indicates that the maneuvering object can pass at the maximum velocity, and a velocity coefficient of 0 means that it is impassable.
In the rules of general assessment setup in this paper, the main factors affecting the passing efficiency of maneuvering objects are vegetation, water system, geology, disaster, and the DEM. Based on the principle of analytic hierarchy process, the influence degree of each factor on the traffic efficiency is obtained by utilizing data, expert scoring, and investigation in the form of the judgmental matrix employing a pairwise relationship between the five factors, which is presented in (1) [31][32][33].
Computational Intelligence and Neuroscience CR � CI RI .
After calculation, the corresponding eigenvectors of matrix A are obtained, as shown in formula (2), and then by formulas (3) and (4), λ max � 5.292 and CI � 0.073 are obtained and n in the formula refers to the number of influencing factors. e RI value is the consistency test value of the matrix, and the value of the fifth-order matrix is 1.12. After substituting it into equation (5), CR � 0.065 < 0.10 is obtained, indicating that the judgmental matrix A satisfies the consistency test, so the influence weights of vegetation, water system, geology, disaster, and the DEM have been found to be 0.087, 0.110, 0.275, 0.357, and 0.171, respectively. e wheeled transport vehicle as an emergency rescue material is utilized for the actual situation.
is paper evaluates the velocity coefficient fv of the maneuvering object employing the impact of factors such as the collection of data, scores obtained by experts, and analysis of data [31][32][33]. Table 2 presents all values.

Trafficability and Passing
Efficiency. Trafficability and passing efficiency are the results of environmental trafficability analysis, which can be calculated using the established passing evaluation rules after environmental modeling and can be used for path searching in a road-free environment.

Trafficability.
In the passing domain model, the passable value is used to represent the trafficability of the maneuvering object, the 0 value is inaccessible, the non-0 value is accessible, and the values between 0 and 1 are used to distinguish different maneuvering objects, for example, vegetation.
e values 1 and 2 of the vegetation part in Figure 2 represent different vegetation types, respectively. Finally, the stack analysis based on the array of traffic domains of each influencing factor is carried out to obtain the binary traffic map for path planning, and the traffic value is either 0 or 1. When the traffic conditions corresponding to each influencing factor are accessible, the traffic value is equal to 1, indicating that the environmental traffic conditions are assessed as good; when the traffic condition corresponding to the influence factor is impassable, the pixel is an impassable area, and the traffic value takes 0.

Passing Efficiency.
Passing efficiency can be measured based on the time used by the maneuvering object through the raster pixel. e calculation of the passing time needs to take into account the influence of environmental factors on the speed of the maneuvering object's movement and establish the formula for calculating the passing speed. According to the size of the pixel and the speed of the moving pixel, the calculation formula of the transit time of adjacent pixels is constructed to support the path transit time calculation. e maximum travel speed of the maneuvering object is V Max , the velocity coefficient is denoted by fv (see Table 2), and the passing speed of a single influencing factor denoted by Vx is defined by W vegetation , W water , W geology , W DEM , and W disaster represent the weight of vegetation, water system, geology, the DEM, and disaster factors in all the factors affected, respectively. When equation (7) is written, we abbreviate the words with the first letter of each word, namely, V, W, G, DEM, and D, respectively. e maneuvering object in each cell passing velocity coefficient fv t can be defined by (based on the actual cell influence factors) By substituting equations (5) into (6), the velocity V of the maneuvering object in each pixel can be expressed by Figure 3 depicts the maneuvering object that can move in eight directions in the current cell. Let R be the pixel length and width, V be the current pixel velocity, and V 1 be the forward pixel velocity. e passing time for adjacent cells is calculated as follows: 2.Time of oblique distance : Computational Intelligence and Neuroscience 5

The Improved A * algorithm
3.1. e Conventional A * Algorithm. A * algorithm that uses a heuristic function for path search is suitable for path planning with known global environment information. A * algorithm uses cost function for path planning whose cost function is defined by In equation (11), F(n) represents the cost assessment from the starting grid point to the target point via any node n, G(n) is the actual distance assessment from the starting point s to the grid n, and H(n) represents the European distance assessment from the grid n to the endpoint e. In equation (12), H(n) is defined as follows: where e x and e y in the formula are the x coordinates and y coordinates of the endpoint e and the n x and n y are the x coordinates and y coordinates of the current point n.
e simplified version of the conventional A * search algorithm is shown in Figure 4. e initial cell is used as a search cell, then searching and calculating the F(n) value of the adjacent child node is conducted, and all the F(n) values of the searched nodes in the list are recorded. e list is arranged by the F(n) value and determines the priority of the nodes. e smaller the F(n) value is, the higher the priority is searched. Repeating this process is necessary until the path planning is completed, so the final path is generated.
In Figure 4, the upper-left corner number of each grid pixel represents G(n), the upper-right corner number represents H(n), the lower-right corner number represents F(n), the lower-left corner number X represents the X search pixel selected according to F(n) value. e arrow designates a movement from the subpixel to the parent pixel.

3.2.
e Improved A * Algorithm. Based on the comprehensive consideration of the efficiency of maneuvering object passage and the efficiency of algorithm searching, an improved A * algorithm is proposed. e main improvements can be summarized as follows: (1) e traffic cost of different maneuvering objects is affected by the surface type and the traffic capacity of the maneuvering object itself; this paper improves the original heuristic function, and the formula is as follows: G q (n) denotes the sum of the costs from the starting point s of the maneuvering object to the grid n, and the generation value of each pixel is obtained by the passing velocity coefficient fv t in equation (7). e formula is defined by (2) e improved A * algorithm takes into account the elevation difference, and the relatively flat terrain is preferred to move forward. H t (n) is obtained by the Euclidean distance and elevation difference, and the calculation is shown in equation (15), where H(n) is the Euclidean distance from the current point n to the endpoint e , which plays a pointing role. e calculation is shown in equation (12). H d is the absolute value of the difference between the parent node elevation H f and the current node elevation H C , which is used to evaluate the terrain fluctuation, and the calculation is shown in equation (16).
(3) Improving the A * search algorithm leads the model to be simplified as shown in Figure 5. e initial cell is used as a search cell, and the F(n) value of the adjacent child node is searched and calculated. en, the cell with the smallest F(n) value is selected as the current search cell. Repeating this process is needed  until the path planning is completed and the final path is generated.
In Figure 5, the left upper corner number of each grid pixel represents G q (n), the right upper corner number represents H t (n), the right lower corner number represents F sh (n), the left lower corner number X represents the X search pixel, and the arrow points from the subpixel to the parent pixel are found.
In this simplified model, G q (n) and H t (n) are represented by distance. When compared with the simplified model of the conventional A * search algorithm, the improved A * algorithm is seen to be not optimal in the planning path distance effect, but the computational efficiency is improved. Particularly, with the increase of the data scale, the advantage of its computational efficiency will be more obvious.
e detailed steps of the improved A * algorithm are presented in Table 3.

Experiments and Analysis
To verify the validity of the trafficability analysis and algorithm, the rescue path planning system of an environmental emergency for a road-free environmental emergency is developed, and experimental analysis is carried out.

Geographical Environment Data.
is experiment constructs the nonroad network environment through the topography data of Shuozhou City, Shanxi Province. Shuozhou City is fully characterized as a loess-covered mountain terrain plateau with complex and diverse natural conditions and obvious transitional nature, which is suitable for simulating complex nonroad network environments. e data used include the raster DEM data and vector data of vegetation, water system, geology, disaster, and residential area depicted in Figure 6.

Data of Maneuvering Object.
e maneuvering object is called a wheeled emergency relief vehicle, and the performance data of the maneuvering object are presented in Table 4.

Results and Analysis.
e environment is modeled using the provided geographic environment data. e roadfree environment model has a grid resolution of 2 meters and about 200 million grid cells. Subsequently, the experimental data are processed by the established pass evaluation rules to generate the passing map of wheeled emergency rescue  Figure 7 depicts the result of the trafficability analysis of the wheeled emergency relief vehicle. Figure 7(a) denotes the binary pass diagram of the maneuvering object, where the black area is the accessible area and the white area is the inaccessible area. When compared with Figures 6(b), 6(d), 6(e), and 6(h), the areas of the residents' land, water system, geological disaster area, and slope value not meeting the prevailing conditions which are determined as impassable by the rules of passing evaluation can be seen. Figure 7(b) depicts the passing efficiency map, which reflects the passing efficiency of the maneuvering object through the color depth of the passable area. e closer to black the color is, the higher the passing efficiency of the maneuvering object is. On the contrary, the closer to white the color is, the lower the maneuvering object's passing efficiency is. When compared to Figure 6, the mountain and hilly areas are more difficult to pass, so the color is close to white, while the more accessible in flatlands are close to black.

Trafficability Results and Analysis.
Covering existing road data on the road map is depicted in Figure 6. e red area is the road data of the area. Besides, Figure 8 depicts most of the existing roads that are in the passable area, which suggests that the maneuvering object passage map generated by this method has reliability.

e Results of the Path Planning
Analysis. While Figure 9 shows a two-dimensional effect of the perspective planning, Figure 10 depicts the three-dimensional effect of it. rough visual analysis, the planning path avoids impassable areas and follows the valley. e improved A * algorithm is compared with the conventional A * algorithm, and the comparison data are shown in Table 5. e search time of the conventional A * algorithm increases exponentially and eventually faces the dimensional disaster; the phenomenon of computational stagnation or collapse occurs when the length of the search path grows. By comparing the search time data, the search efficiency of the improved A * algorithm is greatly improved, especially in long-distance path search. When comparing the path distance data, the path distance planned by the improved A * algorithm that is relatively longer is found, and the phenomenon is analyzed based on the algorithm design level. e improved A * algorithm with strong timeliness that takes into account the influences of both slope and the weight of passing efficiency as a goal to achieve the task will bypass the planning scheme when a high slope difficult to pass is encountered in the terrain. is can be supported by the data in Table 5. When comparing the total time data, it is found that the improved A * algorithm embodies the characteristics of strong timeliness. Besides, by comparing some path distances and estimated time data results, it is discovered that the estimated traffic time of A * is shorter, although the path distance is longer than the traditional A * algorithm. When most of the slope proportion and geological type proportion data are compared, it is found that the improved A * algorithm tends to choose a smoother and easy-to-pass terrain. In summary, the improved A * algorithm has a great theoretical significance and application value to the rapid planning of rescue paths for disaster emergencies in a road-free environment. e process for the improved A * algorithm Input: Environment digital model, starting point, and target point (1) Initialize the cell (2) If initialization is marked as 0, indicate that the grid has not been evaluated (3) Else initialization is marked as 1, indicating that it is represented as a point on the path or impassable (4) Create a container m_openList that holds the path node (5) Add the starting point to the m_openList as the current node (6) While (the current node is not the target point) (7) View the current node adjacent to passable child nodes (8) If (there are child nodes in m_openList) (9) ere is a cross in the path, so remove the cross section (10) Else if (child nodes exist passable) (11) Calculate the F sh (n) value of the adjacent passable subpoint; see formula (13) (12) Select the adjacent passable subpoint with the smallest F sh (n) value (13) Add to the m_openList and use it as the current point of the path (14) Else (15) Return the previous node and select the minimum F sh (n) value among the remaining child nodes (16) Falling back to the starting point means that no path has been found (17) Return failed to find a way (18) Return m_openList (19) Traverse the nodes on the path to get path information (20) Calculate the estimated passing time of the path by using formulas (9) and (10) Output: e path from the starting point to the target point

Conclusion
e rapid planning of the path of disaster relief personnel, vehicles, and other mobile objects in the road-free environment is of great significance to ensure the safety of the lives of the affected people and the supplies of materials. erefore, since planning those efforts in the current road-free environment looks more complicated than the areas having road network, optimizing the characteristics of emergency rescue path planning with maneuvering objects requires putting forward a better implementable method for disaster emergency rescue path planning in the road-free environment, which mainly contains the environmental trafficability analysis dealing with the issues of timeliness and traffic efficiency.
us, the improvement of the currently implemented A * algorithm is proposed in this paper, which is called the improved A * algorithm. e environmental trafficability analysis of the road-free network has laid a data foundation to enhance the A * algorithm. e improvement is constructed based on the issues that have not been visited yet, which are called the trafficability analysis and strong timeliness. To resolve those issues, the rapid planning of the emergency rescue paths that mobile objects follow in the road-free environment is conducted based on the trafficability of the area. us, the traffic efficiency and timeliness are improved by heuristic search optimization and grid method by focusing on the influential weights. On the other hand, the A * algorithm generally takes a longer time to find a rescue path than does the improved A * algorithm when long distance is a concern. Besides, the A * algorithm generally fails under the search for long distance or does not generate practical solutions.
By experimenting and analyzing the method using the data of Shuozhou City, Shanxi Province, the results show that the environmental trafficability data of the road-free network have certain rationality and reliability, while the search efficiency of the improved A * algorithm is greatly improved, the time-sensitive characteristics of the character and emergency rescue are greatly improved, and the influence of slope and weights of passing efficiency is taken into account by the improved A * algorithm, which makes the planned path more flexible and practical. e research has contributed a certain theoretical significance and application value, which can provide a reference for the rapid planning of disaster emergency rescue paths in the environment of no road network.
Although a set of preliminary feasible solutions have been formed in the planning and research work of disaster emergency rescue paths in the environment of no road network, further research is needed to enhance the proposed algorithm. e next round of research will be focusing on improving the rules of passing an evaluation, optimizing the algorithm of the path planning, and advancing the rapid planning ability of disaster emergency rescue paths in road-free environments.

Data Availability
e data used to support and prove the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.