An Optimized Probabilistic Roadmap Algorithm for Path Planning of Mobile Robots in Complex Environments with Narrow Channels

In this paper, we propose a new path planning algorithm based on the probabilistic roadmaps method (PRM), in order to effectively solve the autonomous path planning of mobile robots in complex environments with multiple narrow channels. The improved PRM algorithm mainly improves the density and distribution of sampling points in the narrow channel, through a combination of the learning process of the PRM algorithm and the APF algorithm. We also shortened the required time and path length by optimizing the query process. The first key technology to improve the PRM algorithm involves optimizing the number and distribution of free points and collision-free lines in the free workspace. To ensure full visibility of the narrow channel, we extend the obstacles through the diagonal distance of the mobile robot while ignoring the safety distance. Considering the safety distance during movement, we re-classify the all sampling points obtained by the quasi-random sampling principle into three categories: free points, obstacle points, and adjacent points. Next, we transform obstacle points into the free points of the narrow channel by combining the APF algorithm and the characteristics of the narrow channel, increasing the density of sampling points in the narrow space. Then, we include potential energy judgment into the construction process of collision-free lines shortening the required time and reduce collisions with obstacles. Optimizing the query process of the PRM algorithm is the second key technology. To reduce the required time in the query process, we adapt the bidirectional A* algorithm to query these local paths and obtain an effective path to the target point. We also combine the path pruning technology with the potential energy function to obtain a short path without collisions. Finally, the experimental results demonstrate that the new PRM path planning technology can improve the density of free points in narrow spaces and achieve an optimized, collision-free path in complex environments with multiple narrow channels.


Introduction
Automatic path planning is an indispensable technology for realizing the intelligence of mobile robots [1,2], which is mainly used to automatically plan a collision-free path for the mobile robot from its initial position to a final position in a complex environment. With the improvement of artificial intelligence in mobile robots, the specific requirements for the final path are gradually increasing, including smooth performance [3], shortest distance parameters [4], and minimum energy parameters [5]. However, finding an effective path without collision in the presence of obstacles has become the foundation of the later research, especially regarding complex workspaces with multiple narrow channels [6][7][8]. Therefore, the path planning through narrow channels has become a key focus of various researchers.
Many researchers have proposed various algorithms to solve path planning problems in complex environments. In the early days, the workspace was divided into grids to find an available path. A* and D* are typical algorithms, but the division into a grid can affect the appearance of obstacles and the direction of movement of the robot [9][10][11]. Meanwhile, inspired by nature, the Artificial Potential Function (APF) algorithm can effectively implement path planning for mobile platforms by constructing the distribution of potential energy fields, including the repulsive and attractive forces in a specified workspace [12,13]. With the development of swarm intelligence technology, researchers have continued to apply different swarm intelligence methods to solve path planning problems, including artificial bee colony [14], genetic algorithms [15], particle swarm optimization [16], and ant colony algorithms [17]. However, compared to other algorithms, these swarm intelligence algorithms may take longer to obtain efficient paths in complex environments. In addition, reinforcement learning methods have been gradually applied to the mobile robot path planning problem [18,19]. However, these algorithms need to obtain position information in real-time and adjust the corresponding state to achieve path planning. Therefore, the algorithm is affected by the performance of its own sensors and the distribution of obstacles. Recently, path planning methods based on random sampling have been widely used and studied, due to their high sampling efficiency and success rate. The most typical methods are rapidly-exploring random trees (RRT) and the probabilistic roadmaps method (PRM), which require use of the sampling points to describe the free space [20,21]. The RRT algorithm randomly selects a free point at each step and uses a separate query process to solve the path planning problem [22]. In contrast, the PRM algorithm constructs the free space and uses a multi-part query method to obtain the path [23]. Compared to the RRT algorithm, the PRM algorithm can quickly achieve a path with the shortest distance by establishing the free space in advance. Therefore, we take the PRM algorithm as the core algorithm and optimize the critical steps of the PRM algorithm by combining the APF algorithm.
The critical steps of the PRM algorithm are divided into two modules: the learning process and the query process [24,25]. The learning process is mainly to construct the free space and obstacle space using different point sets, which are determined by randomly sampling points in the workspace and using collision detection technology [26]. The collision-free line represents the local path, which is composed of the free point set. Therefore, the free space of the specified workspace consists of free points and collisionfree lines. The critical step in the path query process is using the local planner to plan the global path without collision in the presence of obstacles. The PRM algorithm can construct the free points and collision-free lines in advance for a specified environment, significantly reducing the required time for collision detection and optimizing the efficiency of path planning [27]. However, after years of in-depth research, the PRM algorithm still has a low success rate for solving path planning problems in complex environments with narrow channels, which is due to the inability of various sampling methods to describe narrow access.
To improve the performance of the PRM algorithm, researchers have proposed various methods to optimize the map construction technology for free space in a complex environment with narrow channels. The first optimization technique proposed the use of various random sampling principles to avoid the inability to achieve full-space distribution, including uniform random sampling technique [28], quasi-random sampling technique [29], and random walk sampling techniques [30]. These random sampling principles mainly increase the density of the sampling points in the narrow channel, according to the distribution of random points. However, it cannot be significantly increased, due to the low proportion of the narrow channel in the overall space. The second approach involves improving the distribution and sampling density in the narrow space, which is the main content of most articles related to optimizing the traditional PRM algorithm. For example, in [31], the Gaussian sampling technique has been proposed to transform the obstacle points into the free points of the obstacle boundary by adding random perturbations. The author of these papers [32,33] has also proposed the bridge sampling technology, which determines the collision-free midpoint position as the sampling point of the narrow passage by sampling between different obstacle points. In [34], Ye has determined the sampling points of narrow passage by using an artificial potential function between different obstacles. The second type of improvement method directly increases the density of free points by sampling near the narrow channel. However, the free points in the narrow channels are limited, due to the large number of obstacles and the significant increase in the surrounding environment. Meanwhile, the free points close to the obstacles increase the probability of collision with the mobile robot.
In order to solve the problem of path planning in complex environments with narrow channels, these above-mentioned studies have mainly focused on how to quickly increase the density of free points in the narrow channel and reduce the collision probability. In this paper, we proposes a new algorithm by effectively combing the APF algorithm and the typical PRM algorithm, in response to the discussed research. The primary process of the new algorithm is based on the PRM algorithm, while the APF algorithm is added to improve the critical operation. Therefore, the new PRM algorithm is not a simple combination of the two different path planning algorithms. The specific contributions of this paper are as follows: (1) To ensure the existence of narrow passages and optimize the number and distribution of free points in the free workspace, we extend the obstacle based on only the diagonal distance of the mobile robot, ignoring the safety distance in the movement. After that, we optimize the distribution of the sampling points using the quasi-random sampling principle in the workspace. Then, we re-classify the sampling points based on the potential energy function, in order to ensure the safety of the mobile robot during movement. (2) To improve the density of random points in the narrow channel, each obstacle point is first subjected to a random force and extended from eight directions. Then, the free points in the narrow channel are determined by combining the characteristics of narrow passages with the potential function. (3) To reduce the time required for path planning, we sample the bidirectional A*, instead of the typical local path planning planner in the learning process of the traditional PRM algorithm. (4) To shorten the distance of the path and ensure that distance from obstacles is maintained, the basic path obtained by the query process requires re-optimization by pruning techniques and potential functions, in order to improve the superiority of the final path.
The remainder of this paper is structured as follows: the Section 2 provides a the specific review of the PRM algorithm in the trajectory planning of mobile platforms in workspaces. In the Section 3, the specific combination of the APF algorithm and PRM algorithm is described in detail, including three kinds of classification of space points, the method of moving obstacle points to a narrow environment, and the optimization of local motion line segments. The Section 4 describes the optimization of the local planner used in the query process, and illustrates the specific steps for optimizing the final path in detail. Meanwhile, the principle for combining the path pruning technology with the artificial potential is described. The Section 5 details various experiments to prove that the new sampling method can effectively improve the density of free points in narrow channels. Furthermore, it is shown that the new PRM algorithm can effectively shorten the required time and reduce the length of the final path. The Section 6 summarizes the paper and describes specific future research directions.

Typical PRM Algorithm
The typical PRM algorithm is a well-known global path planning algorithm that can ensure that a robot finds a collision-free path from its specified initial position to a corresponding final position. In the typical PRM algorithm, the description of the free space is directly described by free points and collision-free line segments, which can effectively reduce the occupation of the workspace for path planning. Meanwhile, the robot can input the offline environment map for pre-processing before the actual work, in order to reduce the required time and improve the map construction efficiency. Therefore, the PRM algorithm has been successfully applied to the actual path planning problem of complex environments. The PRM algorithm is a multi-query random sampling algorithm which includes two fundamental processes: the learning process and the query process. These aspects must be analyzed and described, in order to improve the performance of the typical PRM algorithm.

The Learning Process
The primary task of the learning process is to construct a free space in the workspace through sampling techniques. The actual workspace of the mobile robot is denoted as C s , which generally includes: the obstacle space, denoted by C obs , and the free space, denoted by C f ree . In the typical PRM algorithm, the free space is constructed by the learning process and can be expressed as C f ree = [nodes, edges], where the nodes represents the sampling points in the free space without collision, which are also called free points, and the edges represents the connection lines in the space without collision with obstacles.
To realize the essential functions in the learning process, the specific steps for obtaining free points are shown in Figure 1a. First, we need to initialize the critical parameters for obtaining the free point set, including: setting the free point sets and the obstacle set in the free space as empty, setting the maximum number of free points as N max , and adding the initial and final positions into the free point set. Next, we choose a sampling point in the environment using the random sampling principle and set the index of the random point as 1 (i = 1). Then, the random point is submitted to the collision detection function. If it does not collide, the point is added to the free point set; otherwise, the point is put into the obstacle point set. Finally, we judge whether the sampling task is finished. When the value is smaller than the number of the free point set, we continue the sampling task, setting the index as (i = i + 1); otherwise, the sampling task is complete. While describing the free space using the free point set, the number and the distribution state of the sampling points may seriously affect the final free space of the workspace. Due to the slow development of sampling principles, researchers have mainly obtained sampling points using random sampling principle, ensuring that the distribution of sampling points covers the whole workspace by increasing the number of random points. After completing the description of the free space using the free point set, we directly use collision-free line segments between free points to represent effective local paths. These collision-free lines can ensure that the robot can perform motion in free space.
The specific steps of building the collision-free line segments are shown in Figure 1b. First, we select the first free point P i (i = 1) and another free point P j (j = 1) in the free space point set, in order to construct a line segment. Next, the line segment is submitted to collision detection to build the collision-free line. When a collision occurs along the line, we directly discard the line segment and continue to another sampling point P j (j = j + 1); otherwise, we add the two endpoints of the path to the edge set. Then, we judge whether the point (P j ) is the end point of the free point set. If the value of j is less than N max , we assign the random sampling index to j = j + 1 and continue to construct a new line segment, based on another free point P j ; otherwise, the point indicates the endpoint of the free point set. Finally, we judge whether the point (P i ) is the end point of the free point set. If the value of i is less than N max , we assign the random sampling index to i = i + 1, and continue to construct a new line segment based on another free point P i ; otherwise, this point indicates the endpoint of the free point set, and the free line segment sampling is complete.
When the number of the free points in free space is determined, the efficiency of collision detection and the total number of line segments may seriously affect the construction time of free line segments. To shorten the required time for collision detection, researchers haved used the binary collision detection methods, instead of the traditional incremental collision detection method. Meanwhile, researchers have also used local line segment connections, instead of overall line segment connections, in order to reduce the number of overall free line segments and the impact on path construction. Through optimization in these two aspects, the required time for constructing the free line segments can be significantly shortened, compared to the typical learning process.

The Query Process
The primary purpose of the query process in the PRM algorithm is to use local path planning to find an effective path by querying the collision-free paths constructed in the learning process, thus obtaining a practical path from the initial position to the final position. As the collision-free local path is constructed through the learning process, the robot does not need to build local collision-free paths in the entire free space. The performance of the local path planner affects the time required for the query process in the PRM algorithm. At present, the principal local path planners in the query process include the A* and D* algorithms. Compared with the D* algorithm, the A* algorithm adds a heuristic function to determine the optimal path, and has higher efficiency for path planning in static workspaces. Therefore, the A* algorithm has become the mainstream local path planner in the query process of the typical PRM algorithm. In completing the query process, we need to evaluate the basic performance of each free point using the typical A* algorithm, which is expressed as: where f (x) represents the performance value of the specified free point; g(x) is the operating cost function, which represents the Euclidean distance (actual cost) from the initial position to the specified position; and h(x) is the Heuristic function representing the distance from the current position to the final position. We can effectively determine the collision-free trajectory using the A* algorithm by adding the initial and target points to the map. The specific steps of the querying process in the typical PRM algorithm are shown in the Figure 2. First, we create the opened-loop list and closed-loop list to save the relevant parameters. The closed-loop list is used to store the free point that will not be considered again, but the opened-loop list is used to store the whole free point considered to find the shortest path. We also assigned the initial parameters to the opened-loop list, including the index of the initial point, the performance parameters, and the number of closed points. Next, we find the overall nearest free points that can form a local path with the initial query point, and estimate the relevant performance parameters of these nearest free points. Then, we select the free point (S 1 ) with the smallest sum of cost functions as the the nearest point with the best performance, and delete the previous query point from the open-loop list. Finally, we check whether the nearest point (S 1 ) has reached the target point; if it does not reach the target point, we add the previous point into the closed-loop list, set the query point as the nearest node (S 1 ), and continue to query the next free point with the best performance; otherwise, we have found the global path and the querying process is complete. Path planning for a mobile robot in a complex environment can be realized through the above learning and query processes. However, finding a collision-free motion path in a narrow passage may be impossible, due to the randomness of the free points in the learning process. However, when the number of sampling points is increased, the sampling time and the number of sampling points may inevitably be wasted. To ensure the mobile robot's safety in the actual movement, we hope that there is a certain distance between the robot and the obstacle to avoid collision. In this paper, we propose an optimized PRM path planning algorithm by effectively combining the artificial potential energy function, in order to solve the above problems.

The Optimized Learning Process
This section mainly details the optimization of the learning process to obtain free points in complex environments, especially inside narrow passages. Then, we improve the collision-free line segment in free space, in order to make it more suitable for the motion performance of the mobile robot. The optimization content of the learning process mainly includes re-defining the initial sampling point, optimizing the obstacle point, and improving the collision-free line segment. The specific steps of sample point acquisition and collision-free path construction in the learning process are shown in Figure 3.

Optiming the Workspace Map and the Initial Point Classification
To complete the path planning task, the traditional algorithm first simplifies the environment into a typical two-dimensional map. To ensure the safety of the mobile robot, the traditional method expands the obstacle space according to the diagonal length of the mobile platform and the safe distance during movement. We used the typical PRM algorithm to obtain the free workspace, which is mainly composed of free points and collision-free line segments, where the free points serve as the basis for forming the collision-free line segments. However, when facing a complex environment with various narrow channels, the above obstacle expansion methods may lead to path planning failure, as the maximum width of the narrow channel may become smaller with the expansion of obstacles. Therefore, the difficulty of obtaining free points in narrow channels based on random sampling methods is increased. For narrow passages that are only suitable for mobile robots, the above obstacle expansion method leads to connections between obstacles. Therefore, some narrow channels cannot be described in the map, resulting in the path planning failure in the complex environment with narrow channels. Therefore, the obstacle expansion method must be optimized to display all narrow channels suitable for traversal by the mobile robot.
Considering the path planning problem in a complex environment with narrow channels, we first re-define the obstacle expansion method in the work environment to retain all narrow channels suitable for the mobile robot. Then, we replace the traditional random sampling technology with typical standard sampling techniques, in order to enhance the distribution of the sampling points. Finally, we decompose the sampling points into three categories, thus reducing the number of the free point, which can ensure the safety performance of the final acquisition path.
(1) Re-define the workspace map In traditional path planning algorithms, researchers generally deal with obstacle expansion according to the safety distance and the diagonal length of the mobile robot when the environment for the mobile robot does not include various narrow channels. For complex environments with various narrow channels, some narrow channels suitable for mobile robot may disappear in the map, especially those that just meet the movement conditions of the mobile robot. This is caused by the coincidence of the data of obstacles during map zoom-out or the obstacle expansion process. In order to obtain all narrow channels for the mobile robot, we must perform the obstacle expansion using the diagonal parameters of the mobile platform. With this method, some narrow channels can assigned as suitable for the mobile robot, while the rest are deemed as not suitable for its motion. We can retain and display all narrow passages where the mobile robot will not collide with obstacles, even if the channel is just large enough to place the mobile platform in stationary state. Therefore, we first directly expand the obstacles based on the diagonal length of the mobile robot, ignoring the movement safety distance.
Through the above method, we can simply obtain the complex working environment as a two-dimensional map with narrow channels. To ensure the safety of the mobile robot located at the sampling point, we then re-classify the sampling points based on the safe distance limit during movement. Meanwhile, we add safe distance constraints to construct the collision-free line segments and optimize the querying process.
(2) Enhance the distribution of the sampling points After obtaining the typical two-dimensional map, it is necessary to conduct multi-point sampling in the environment through use of random sampling principle, for which it is necessary to determine the type of sampling points. As a typical spatial sampling technique, the random sampling method has been widely used for the PRM map construction step. This method can realize effective sampling of multi-dimensional environments. However, due to the randomness of the sampling points, it is necessary to set a higher sampling number to ensure full-space sampling, which can waste time and reduce the efficiency of the PRM algorithm. Meanwhile, the random sampling method defined by Monte Carlo method cannot be applied to the sampling of complex environments with narrow channels. A low number of sampling points cannot effectively distribute multiple free points in the channel, while a high number of sampling points can increase the density of free points, resulting in concentrated free points in free space. Therefore, we must enhance the sampling point distribution to adequately realize sampling of the whole working environment.
To improve the uniformity of the sampling points over the whole space, we adopt the quasi-random sampling principle, instead of the typical random sampling principle, for the free space. The quasi-random sampling principle defined by the Monte Carlo method is currently mainly used to solve global optimization problems, which considers the discrete degree of the random point distribution. This new method has been shown to possess the same efficiency as the typical random sampling method. The sampling points based on the quasi-random sampling principle are also random, but the dispersion degree is higher in the spatial point set. Halton sets are some of the most typical point sets obtained by quasi-random sampling. Once the number of random sampling points is set, this method can ensure the uniform distribution of all sampling points throughout the entire workspace, which means that all sampling points are uniformly distributed inside the whole working environment. In the specified workspace, the sampling point set can be obtained by the Halton sets defined by the quasi-random sampling principle, as has been shown [35]. In brief, we choose d distinct primes P 1 , P 2 , · · · , P d , and the set's i th point is given by: where P i represents the sampling points in the sampling space and r p (i) is obtained by writing the digits of the p-ary notation for i in reverse order, which can be expressed explicitly as: where i is a real number written in p-ary notation. We set the number of sampling points to 100 and carried out random sampling in the specified workspace using the two methods. The sampling results are shown in Figure 4. By comparing the different distribution results of random points, we found that: (1) Some random points obtained by the random sampling principle in the free space have a shorter distance than the ordinary distance between sampling points based on the quasi-random sampling principle, as shown in Figures 4 and 5. Therefore, randomly sampled points based on the random sampling principle may be concentrated together, as can be seen from Figure 4. (2) Due to the randomness of the random sampling principle, some free spaces without any sampling points occur, as shown in Figure 4, while areas surrounded by sampling points do not have a large space, as shown in Figure 5. Therefore, the quasirandom sampling principle can effectively reduce the concentration of random points, compared with the typical random sampling method. The optimized method evenly distributes all the sampling points in the overall space, in order to improve the efficiency of the learning process.  (3) Re-classify the sampling points Through the new obstacle expansion technology, we obtain a variety of narrow channels that can hold the mobile robot at rest. To meet the actual movement requirements, we need to take the safe distance of the robot's movement into account in the path planning process. To avoid collisions between the mobile robot and the obstacles during movement, we optimize the free points obtained by the quasi-random sampling principle. Considering the safety distance during actual movement, we decompose the traditional sampling points into three categories: free points, adjacent points, and obstacle points. The obstacle point set is the same as the typically defined one, meaning that the sampling points fall into the set of random points in the obstacle space. However, the free point set not only needs to fall into the free space, but also requires the shortest distance between these points and the nearest obstacle to be greater than a specified distance. Then, the adjacent point set is the rest of the random points in the free space, which represents random points in free space whose shortest distance to the obstacle is less than the specified distance. Comparing the traditional classification method in Figure 6 and the optimized classification method in Figure 7, the adjacent points cannot guarantee safe distance limits between the mobile robot and obstacles. Therefore, the adjacent points need to be discarded, and cannot be used to construct the collision-free line segments.  The three different point sets proposed in this paper are not only related to the distribution in free space, but also to the closest distance to obstacles. According to the distance constraint between the obstacle and the sampling point, we use the artificial potential energy function to accurately classify the sampling points, which requires the analysis and description of the repulsive force generated by obstacles in the workspace. The artificial potential field is typically used to construct an effective map and find the shortest path, according to the natural laws of attraction and repulsion. However, the artificial potential field in this paper mainly optimizes the critical steps in the PRM algorithm. In the sampling process, we use only the repulsive function alone to construct the repulsive force field, in order to classify all sampling points. For any specified workspace, we can construct a potential energy map based on the basic information of the map and the potential function.
Let a random sampling point in the free space be represented as P o , and the i th obstacle point be expressed as P osi . Then, the distance between the two points is: where (x o , y o ) is the position of the random point in the two-dimensional environment, (x osi , y osi ) represents the position of the i th obstacle, and d represents the nearest straightline distance from the random point to the obstacle. In the potential function, the repulsive force will acts only on the nearby points around the obstacle. Therefore, the repulsive force on the points far away from the obstacle's position is equal to zero. The repulsive force at any point in the space can be specifically expressed as: where k ri represents the repulsion scale factor, d ri represents the minimum distance between the point and the i th obstacle, and d set represents the influence distance of each obstacle. For the workspace with obstacles shown in Figure 8, a corresponding potential map containing only contained repulsive forces can be established through the above potential energy function, as shown in Figure 9. Once the points in the free space are far away from the obstacle (greater than d min ), the value of the repulsive force is zero, based on the potential map in this figure. In a specified field, the repulsive force increases with decreasing distance from the obstacle. When the point falls on the boundary or the sounding of the obstacle, the repulsive force is at its largest. Therefore, using the value range of the potential energy function can effectively enhance the classification of all sampling points.  After creating a three-dimensional potential map using the artificial potential function, we can use the quasi-random sampling principle to obtain random sampling points from the workspace. These sampling points can then be decomposed into three different sets, according to the specified potential force, which can be expressed as: where p f i , p n i , and p o i represent free points, adjacent points, and obstacle points, respectively; S f , S n , and S o represent the set of free points, the set of adjacent points, and the set of obstacle points, respectively; F() represents the repulsive force based on the artificial potential energy function; and F max represents the maximum repulsive force generated by the specified distance. For any specified workspace, we can obtain three different point sets by combining the quasi-random sampling principle and the optimized classification technique proposed in this paper. The solution process for the the optimized classification technique is illustrated in Figure 3a. First, we use the artificial potential function to solve the repulsive force and then construct the three-dimensional potential energy map based on the obtained repulsive force. Next, we determine the three different point sets and the total number of overall sampling point sets. Then, we obtain one sampling point, based on the quasi-random sampling principle, and estimate the repulsive force at the sampling point. Finally, the corresponding point is classified into the specific point set, based on the potential function, and we determine the total number of sampling points. If it is less than the initial setting number, we continue to acquire new sampling points; otherwise, the sampling process is complete.
For the workspace shown in Figure 8, the above method was used to construct free space and obstacle space through the three different sampling point sets, and the specific results are shown in Figure 10. In this figure, blue points represent free space points, pink points represent adjacent space points, and red points represent obstacle points. Among these three types of sampling points, free points can be used to construct free space, obstacle points can become sampling points in narrow passages by optimizing these obstacle points (as proposed in the next subsection), and adjacent points are closer to obstacles only can be abandoned. We add these three types of random points into the two-dimensional workspace to directly observe the distribution of these sampling points. Figure 11 specifically shows the classification results of these random points in a two-dimensional working environment.

Optimizing the Obstacle Points by Artificial Potential Field
In the typical PRM algorithm, the obstacle points are generally discarded directly, and the number of sampling points in the narrow channel is increased by increasing the number of total sampling points. In the optimized PRM algorithm based on Gaussian sampling, the obstacle points are moved to the surroundings of the obstacle by adding disturbance, inevitably increasing the number of free points in the narrow channel. Another optimized PRM algorithm, based on bridge sampling, directly adds the midpoint of the two obstacle points collected by different obstacles to the narrow channel. In the sampling processing of the narrow channel, the above methods must perform sampling multiple times,while the number of sampling points in the narrow channel cannot be directly increased. In this paper, the obstacle points are directly optimized by combining the characteristics of obstacle points and narrow passages, thereby effectively increasing the sampling number of free points in the narrow channels.
To optimize the obstacle points, we need to directly move the obstacle points into the free space to become typical free points. Meanwhile, we can effectively transform the obstacle points into adjacent points through the effective combination of the force of the random point in space and the specified minimum repulsion force. First, we set a random force on the obstacle point to determine the initial movement direction, as shown for point 1 in Figure 12. Second, the obstacle point is required to move, gradually increasing its distance according to the initial force direction. Finally, the optimized obstacle point can be obtained by comparing the minimum repulsive force for determining adjacent points and the repulsive force of the obstacle point moving along the force direction. Therefore, the repulsive force and the movement distance of the obstacle point can be expressed as: where − → F represents the direction of the random force on the obstacle point, F 0 represents the initial repulsive force of the obstacle point, F S represents the repulsive force after moving some distance, and d represents the moving distance of the obstacle point.
The obstacle point can be guaranteed to become a free point after moving to free space under a random force; however, the free point may be transformed into a narrow passage or a free space far away from the obstacles. To increase the density of sampling points in the narrow channel, we must actually guarantee that the free point falls in the narrow channel. As a narrow passage is composed of at least two obstacles, the obstacle point will inevitably collide twice during movement in the narrow passage. Similarly, two numerical transformations of the repulsive force must occur, compared to the specified repulsive force during movement along the direction of the virtual repulsive force. Based on the characteristics of narrow passages, the corresponding positions determined by the numerical transformations must be in the narrow passage; otherwise, they are outside the narrow passage. The specific sampling results are shown in Figure 12. The specific constraints for adding these free points in a narrow passage can be expressed as follows: where f () represents the repulsive force at any point, P 1 represents the first adjacent point, P 2 represents the second adjacent point, d represents the moving distance, and F set indicates the maximum repulsive force. An obstacle point optimized by direct movement cannot be guaranteed to transform into an adjacent point, due to the uncertainty of the random force; for example, see Figure 12.
To improve the maximum possible points falling into the narrow passage, we continue to add other motion directions, different from the basic direction of the random force. Then, the obstacle point performs the motion transformation in these directions at the same time. Finally, we select the direction that first occurred in these two numerical transformations as the best motion direction, and regard the transformation points as two results after optimizing the obstacle point. The specific solution results are shown in Figure 13.  According to the basic principle detailed above, the critical technology of obstacle point optimization lies in the determination of the random force direction and the change in repulsive force of the obstacle point moving along the eight directions. After acquiring the obstacle point set, we can count the total number of obstacle points and optimize each point as free points in narrow channels. The specific solution process for optimizing each obstacle point is shown in Algorithm 1. In the process, we first assign an random force to the obstacle point and re-assign a random point whose force is zero. Next, we solve the actual angle of the random force and solve in all eight directions, according to the initial angle of the random force. Finally, we expand the obstacle point in eight directions gradually, and count the state transformation between the repulsive force of the expansion point and the specified repulsive force. When two state transformations occur first, we take the corresponding two positions as the optimization result of the obstacle point. In the potential map based on the above workspace, we used the above obstacle point optimization technology to move all obstacle points into the free point set; the specific results are shown in Figure 14. According to the distribution of the sampling points, we can see that (1) most obstacle points were effectively moved from the obstacle space to the narrow channel, effectively increasing the density of the sampling points in the narrow channel, and (2) some obstacle points may be misjudged when the obstacle points fall in the surroundings of the workspace, because the boundary of the workspace is set as an obstacle. However, the probability of sampling in the surrounding boundary by the quasi-random sampling method is low. Therefore, the above phenomenon occurred only for a few obstacle points, thus having a low impact on the optimization of all obstacle points. The optimized results of all obstacle points were mapped into the actual workspace to directly observe the free point distribution results, as shown in Figure 15. According to the distribution of sampling points in the narrow channel before and after the optimization of the obstacle point set, it can be seen that the obstacle point optimization technology proposed in this paper can significantly improve the sampling density of free points in the narrow passage.

Improving Collision-Free Line Segments Using the Repulsive Force
After adopting the quasi-random sampling method and obstacle point optimization processing, the density of sampling points in the narrow passage was significantly increased, and the distribution of sampling points in the free space did not appear to be concentrated. Next, we performed connection processing, based on the above-mentioned free point set. This process aims to construct the local path through line segments that do not collide with the obstacle space. Some researchers have directly connected one free point with all other points in the free space and performed collision detection to obtain line segments. Furthermore, the two free points and line segments that do not collide may be used together to build a free space. In this process, the number of line segments is large after connecting all points, and the local path needs to be segmented during collision detection. Therefore, the number of paths and the length of the total paths will inevitably lead to a high time being required to construct the free space. To date, researchers have reduced the number of overall segments by designing a specified range for one free point for connection processing and improving the collision detection performance by adopting the binary collision detection technique, instead of the incremental collision detection technique.
In the actual construction of line segments, the connection state in the large free space is less complex than the connection state in the narrow passage. It is easy to collide with the surrounding environment when connecting the free points in the free space to free points in the narrow passage. At the same time, the mobile robot does not want to move at these positions near to the obstacle, even if it does not collide. Considering the above-mentioned actual motion requirements, we propose a new type of connection state, considering the three essential types of sampling points and potential energy functions, in order to avoid collision and to stay away from the obstacle space.
The sampling points obtained by the optimized sampling technology include free point sets, adjacent point sets, and obstacle point sets. In the actual line segment construction process, we need to ensure that all sampling points in the free space do not collide and that the distance between a given free point and the nearest free points is greater than a specified distance. Therefore, we can only select the free points as sampling points during the connection process, such that the adjacent and obstacle point sets need to be discarded. Meanwhile, we connect the specified free point with the other free points in the deigned specified space to construct line segments, in order to optimize the time required for connection processing. According to the repulsive function, the other free points in the deigned specified space can be specifically described as: where P o represents a free point in the free space, p oi represents the other free points in the deigned specified space, d oi represents the distance between the two free points, and d set represents the radius limit for the specified space. The positions of two endpoints meet the collision-free conditions when constructing a line segment. Moreover, the line segment constructed by the two points needs to ensure that no collision with obstacles occurs, and the distance between the line segment and the nearest obstacle should be no less than the specified distance. According to the repulsion function, the distance between the line segment and the nearest obstacle needs to be limited, in order to completely avoid the possibility of collision, which can be expressed as: where edge represents a free collision-free line segment in space, maxF represents the maximum repulsive force of the line segment, F() represents the repulsive force at any point, and F set represents the repulsion based on the shortest distance to the obstacle. After obtaining the sampling points of the free space, we used the typical connection technique in the PRM algorithm and the optimized connection technique proposed in this paper to construct line segments in free space, as shown in Figure 3b. The distribution results of the collision-free line segments are shown in Figures 16 and 17. Due to the randomness of movement when optimizing the obstacle points, we directly selected the free space to observe the line segments, and found that: (1) The optimized connecting technique can effectively avoid collision with obstacles by limiting the distance between the line segment and the nearest obstacles.
(2) Some free points may not be able to construct a collision-free line segment, due to the limitation of the maximum repulsive force, which can reduce the number of total line segments and optimize the required time for collision detection. (3) Even if the force used for optimizing the obstacle point is random, the maximum repulsive force can effectively ensure that the free line segment is kept outside the specified range from the obstacle, especially in the narrow passage.

The Optimized Query Process
The query process of the PRM algorithm mainly obtains a successful path by querying the collision-free line segments. To reduce the required time for the query process and improve the effectiveness of the query process, we adopts the bidirectional A* algorithm, instead of the typical A* algorithm. Meanwhile, we combine a path pruning technique and a potential function to shorten the length of the final path and reduce the possibility of collision between the final path and the obstacle space. The specific principles for optimizing the query process are detailed in the following.

Optimizing the Query Process
After constructing the free space through the learning process, we need to solve the corresponding three performance estimated values of each free point to find the optimal path. Therefore, the complexity of solving the query process is high and the required time for solving the path planning is long. To optimize the performance of the PRM algorithm, we not only must ensure the superior performance of the original A* algorithm used in the query process, but also need to improve the speed of the local path planner. Therefore, we use the bidirectional A* algorithm to establish a local path planner both at the initial position and the final point position. The path performance evaluation functions of the two local planners can be specifically expressed as: (11) where f Init () represents the evaluation value of the A* algorithm from the initial position, f Last () represents the evaluation value of the A* algorithm from the final position, g Init () represents the Euler distance from the free point to the initial position, g Last represents the Euler distance from the free node to the final position, h Init () represents the Manhattan distance from the free point to the final position, and h Last () represents the Manhattan distance from the node to the initial position. The typical A* algorithm generally establishes a local path planner from the initial position. Once the final position is reached, the path planning is successful. In contrast, the bidirectional A* algorithm utilizes two local path planners, from both the initial and final positions. Once a common location point exists in the two closed lists established for the bidirectional A* algorithm, the path has been successfully queried. The specific process of the bidirectional A* algorithm is shown in Figure 18. First, we created four different lists to store the different free points in the query process, including Opened_List_Last, Opened_List_Init, Closed_List_Init, and Closed_List_Last. Closed_List_Init is an empty list used to store the path nodes queried from the initial position, Closed_List_Last is an empty list used to store the path nodes queried from the final point position, Opened_List_Init assigns the initial position and corresponding cost parameters, and Opened_List_Last assigns the final position and corresponding cost parameters. Next, we query the nearest free point from the initial point and determine the next point with best performance for the initial point. In this case, we obtain the neighborhood of the initial point and determine the nearest node (PS) with the best performance. After obtaining the best point (PS), we remove the point from Opened_List_Init and add the previous point into Closed_List_Init. Then, we query the nearest free point from the goal point and determine the next point with the best performance for the goal point. In this case, we query the neighborhood of the goal point and determine the nearest node (PF) with the best performance. After obtaining the best point (PF), we remove the point from Opened_List_Last and add the previous point into Closed_List_Last. Finally, we determine whether the closed list Closed_List_Init and the closed list Closed_List_Last contain the same point. If there is a matching point, the path planning is successful; otherwise, we assign the best parameter points to the initial and final points, in order to query the next two points with the best performance.

Optimizing Path Pruning Technology
We can obtain a reachable path from the initial point to the target point by adopting the above optimization query method. However, there may be redundant nodes between adjacent line segments of the initial path. By observing the optimized path after removing some redundant nodes, the corrected path can still ensure the effective motion of the mobile robot. The above phenomenon occurs as any free points may select the nearest free points to connect directly, rather than using the free points outside the specified space during the connection process. Therefore, optimizing the initial path obtained by the query process is necessary.
The path pruning technology is a typical and effective method for optimizing the initial path obtained by the query process. Typical path pruning technology mainly prunes the redundant nodes of the path based on the basic principle that the sum of the two sides of a triangle must be greater than the third side. Therefore, the path pruning technique can form a new optimized path closest to the shortest path. The initial path obtained by the query process is mainly described in terms of free points. Therefore, the path pruning technique can reduce the required number of free points and optimize the memory required to describe the optimal path. In practice, we must consider whether the connection between the first node and the second node collides in the process of processing any three connected node paths. We also need to ensure that the maximum potential energy function of the optimized node connection is less than the specified value, in order to avoid possible collisions. To ensure the excellent performance of the optimal path, we optimized the pruning technology by combining the triangle principle and the potential function. The basic principle for the optimization of the three connection nodes can be specifically expressed as: (12) where F() represents the repulsive force function, F max represents the maximum repulsive force of the line segment, A, B, C represents the three adjacent free points, and F set represents the specified repulsive force.
Assuming that the path obtained by the path query method is composed of N points P 0 , P 1 , P 2 , ..., P N , the path pruning technology used in this paper can be described explicitly, as in Algorithm 2. First, we initialize the data parameters, mark the initial point as Start = P 0 , and set the pruning path as Path_s = [P 0 ]. Second, we connect the free node (Next(Path(i = 2, :)) to the start point with a line segment, and perform collision detection on the line segment. When the line segment collides with an obstacle, we reset the Start point as the location of the previous point and re-construct the path. Then, we solve the maximum potential force of the overall line segment and use Formula (12) to determine whether the conditions are met. When the formula is satisfied, the free point can be optimized. In this case, we reset the Start point as the location of the previous point and re-construct the path. Otherwise, we reset the parameter i as i = i + 1 and mark the point (Next) by the new free point. Finally, we finish path pruning when the value of i is equal to the total number of points in the initial path. if Collision_LineStart,Next,map== false then 6: if max(Rep_Force_Line(Next, Fins)) > F set then According to the workspace in Figure 8, we used the typical path pruning technology and the optimized path pruning technology described above to re-optimize the initial path and obtain the optimized path; the specific results are shown in Figures 19 and 20. By comparing the different paths before and after optimization, it was found that: (1) The number of free points in the optimized path was significantly reduced, especially in the free space with a great area, where redundant free points were deleted, thereby effectively shortening the actual length of the path. (2) The optimized path still did not collide with the obstacles and satisfied the specified distance limit. The above results effectively illustrate that the optimized path pruning technology proposed in this paper can significantly enhance the actual performance of the initial paths obtained by the local path planner.

Simulations
In this section, we mainly verify the performance of the proposed optimized PRM path algorithm. We design various environments for path planning tasks, including single-channel path planning problems, multi-channel trajectory planning problems, and conventional environmental trajectory planning problems. The simulations designed in this section were run on a personal computer with a 3.40 GHz Intel-Core (i5-7500) CPU and 8 GB memory. The specific working map and experimental results are detailed in the following.

Single-Channel Path Planning Problems
In the first simulation, we chose two workspaces with single channels as simulation environments, as shown in Figure 21(a1-b1). The two workspaces have the same size (500 × 500), where the black space indicates the existence of obstacles, and the white space indicates the environment without obstacles. There is only one narrow passage in these workspaces, but that in Figure 21(a1) is a straight narrow channel, while the one in Figure 21(b1) is a winding narrow channel. To verify the effectiveness of the PRM algorithm proposed in this paper, we chose three typical PRM algorithms to conduct path planning for these workspaces. These three methods included the PRM algorithm based on random sampling, the PRM algorithm based on Gaussian sampling, and the PRM algorithm based on bridge sampling.
First, we compared the characteristics of the free points and obstacle points in these workspaces, in order to verify the superiority of the sampling principle in this paper. Therefore, these sampling points were obtained separately by the quasi-random sampling method and the random algorithm method. The specific results are shown in Figure 21(a2,b2) and (a3,b3). By comparing the distribution results of different sampling points in Figure 21(a2,b2) and (a3,b3), the random sampling principle led to the concentration of free points in free space, such that the obstacle space could not be entirely described by the sampling points. When increasing the number of sampling points, the concentration of sampling points in free space must be more obvious, but the number of sampling points in the narrow channel may be significantly increased. However, the method proposed in this paper randomly and uniformly sampled the obstacle space and free space, thus improving the distribution characteristics of sampled points.
To demonstrate the performance of the optimized sampling algorithm, we chose the distribution of the sampling points in the narrow channel and path planning results in two workspaces for analysis and comparison with the other PRM algorithms considered in this paper. First, we selected the different PRM algorithms to solve the path planning problem, and each method was performed 20 times. Then, we used statistical methods to obtain various parameters for the four PRM algorithms in the different workspaces, including the mean number of all sampling points (Ns), the mean number of free points (Nf), the mean number of the free point in the narrow channel (Nc), the mean time required for completing the sampling point (Ts), the mean proportion of free points in the narrow passage to the total free points (Pf), the mean time required for constructing free space (TL), and the success rate of path planning (Suc). These results are provided in Table 1. Finally, we recorded the query process time, as given in Table 2, where TQ1 indicates the time taken to obtain the path by a typical query method, and TQ2 indicates the time taken to obtain the path using the optimized query method.   (1) By setting the same number of sampling points, the optimized sampling technology proposed in this paper can increase the number of free spaces to find the collision-free path. The PRM algorithm based on random sampling principle directly sets the total number of free points and the obstacle points. The other two PRM algorithms mainly directly set the number of free points. (2) By increasing the number of sampling points, the number of sampling points in the narrow channel can be rapidly increased. This is because the sampling optimization technique in this paper can turn obstacle points into free points in narrow passages. (3) By comparing the number of free points in the overall workspace to that in the narrow channel, the quasi-random sampling principle and obstacle point optimization technique proposed in this paper can significantly increase the number of sampling points in the narrow channel. (4) By comparing the success rate of finding the path, it was found that the method proposed in this paper can enhance the success rate of path planning, even when using a lower number of sampling points than the other PRM algorithms. (5) By comparing the required time for the query process shown in Table 2, it can be seen that the method proposed in this paper took less time to find an effective path, effectively increasing the efficiency of the PRM algorithm.
We chose the above four typical PRM algorithms to intuitively describe the superior performance of the algorithm proposed in this paper. The typical PRM algorithms included the PRM based on random sampling, the PRM based on Gaussian sampling, and the PRM based on bridge sampling. To compare the performance of these PRM algorithms, we set the number of sampling points as 200, and obtained the actual distribution of the sampling point and the path planning results. By comparing the sampling results, it was found that: (1) Comparing the free space construction results in Figures 22(a1-d1) and 23(a1-d1), the optimized sampling technology in this paper directly increased the sampling number in the narrow passages. Meanwhile, the collision-free local path in free space ensured that the shortest distance to obstacles was within the specified range. (2) By comparing the path query results in Figures 22(a3-d3) and 23(a3-d3), the optimized query algorithm proposed in this paper could effectively obtain the initial path, showing no collision with the obstacle space. The distance from the nearest obstacle satisfied the specified distance constraint, ensuring the basic performance of the optimized path. (3) By comparing the results before and after optimizing the initial path in Figures 22(a3-d3) and 23(a3-d3), the path optimization technology used in this paper could effectively delete redundant points in the initial path, meet the distance limit for the obstacles, shorten the length of the path, and optimize the actual.

Multi-Channel Path Planning Problems
After completing the single-channel path planning problem, we designed a more complex path planning problem involving two multi-channel environments, in order to demonstrate the effectiveness of the proposed PRM algorithm. In these experiments, the size of the environments was still 500 × 500, where the first map was composed of two winding channels, and the second workspace included three narrow channels, as shown in Figure 24. To prove the effectiveness of the optimized PRM algorithm proposed in this paper, we chose one typical PRM algorithm and two optimized PRM algorithms to solve the path planning problem for a mobile robot in complex environments with multiple narrow channels.   We set the number of samples to be 100, 200, or 300, and used the three PRM algorithms to perform sampling in different workspace. Different results related to the learning process and query process are summarized below, and can be seen in Tables 3 and 4; the meanings of the variable in the two tables are the same as in the first experiment.  According to the comparison between the number of sampling points in the narrow channel and the actual results of different path planning problems in Tables 3 and 4, it was found that: (1) By increasing the number of sampling points, the number of free points in the narrow channel increased significantly in these PRM algorithms. However, the free points in the narrow channel could be increased directly by optimizing the obstacle points after obtaining the overall sampling points. Therefore, the optimized PRM algorithm proposed in this paper increased the number more than the other algorithms. (2) As the number of sampling points increased, the required time for the learning and query processes increased in these PRM algorithms. Nevertheless, the optimized PRM algorithm proposed in this paper could successfully complete path planning with a lower number of sampling points. (3) By comparing the success rate with similar time for constructing the free space, the solution success rate of the proposed algorithm was higher than that of the other algorithms within a similar time. For example, comparing Task4 with 100 points based on the new algorithm and 300 points based on the other algorithms in Table 3, the success rate of the optimized PRM algorithm reached 100%, while the success rate of the other algorithms was less than 80%. (4) By comparing the required time for the learning process, the sampling time was higher than for the other PRM algorithms. However, the other PRM algorithms could not effectively solve the path planning problem in a complex environment. For example, by setting the same initial number sampling points as 300 in the various PRM algorithms for solving Task3 and Task4, the corresponding solution success rate of the optimized PRM algorithm reached 100%, while the other PRM algorithms potentially could not solve the problem. (5) By comparing the required time for the query processes, as shown in Table 4, the required time for the optimized query process was shorter than that of other PRM methods. The optimized query process proposed in this paper can effectively improve the required time for finding the initial path in the narrow channel.
To directly observe the results of the three PRM algorithms, we set the number of sampling points to 300 and performed path planning in the complex workspaces with multi-narrow channels, as shown in Figure 25 and 26. Comparing and analyzing the specific results, we found that: (1) The density of free points obtained by Gaussian sampling was still near the obstacles, but the PRM algorithm based on bridge sampling reduced the density of free points in the narrow channels, due to the existence of wide channels in the environment. Similarly, the optimized PRM algorithm in this paper reduced the number of free points in the narrow channels due to the existence of wide channels. Compared with other methods, the number of free points obtained by the optimized PRM algorithm in the narrow channels was still the largest, as shown in Figure 25(a1) and 26(a1). (2) By comparing the distribution of the free points and collision-free line segments in Figure 25(a2-d2) and 26(a2-d2), the optimized PRM algorithm proposed in this paper effectively avoided collision with obstacles, thus improving the performance of the local path. (3) By comparing the paths obtained by the query process in Figure 25(a3-d3) and 26(a3-d3), there was a great difference between the final path obtained by the optimized PRM algorithm and the final paths obtained by other PRM algorithms, due to the deletion of different redundant points. However, the optimized PRM algorithm effectively ensured the a certain distance from the nearest obstacle boundary was maintained. Therefore, the path optimization technology proposed in this paper can effectively improve the performance of the initial path obtained in the query process.

Path Planning of Normal Workspace
After completing path planning in the single-and multi-channel scenarios, we verified the effectiveness of the optimized PRM algorithm in normal environments. Therefore, we selected three conventional workspaces to solve the path planning problem in a specific environment, as shown in Figure 27. These three workspaces were composed of multiple obstacles. The first map included multiple significant obstacles, with a working environment more complex than the other workspaces; meanwhile, the second and third maps included multiple grid obstacles. The size of the three workspaces was 500 × 500, the initial position was (10,10), and the final point position was (490,490). We use the optimized PRM algorithm proposed in this paper to solve the path planning problem, effectively demonstrating the effectiveness of the PRM algorithm proposed in this paper in normal environments. The results of the optimized PRM algorithm in three normal workspaces: (a1-a3) Sampling points in the three environments; (b1-b3) free space construction results in the three environments; and (c1-c3) query process results in the three environments. Figure 27 represents the critical results for solving the path planning of the three regular workspaces. In these figures, the blue dots represent free points, while the pink dots are adjacent points. Otherwise, the green line segment is the collision-free line segment in free space, the blue path is the initial path successfully obtained by the query process, and the red path represents the final optimized path. By analyzing the results of path planning for the three workspaces, we found that: (1) Comparing the distribution and the density of sampling points in Figures 27(b1,c1), the number of obstacles increased significantly, but the volume of obstacles decreased. Fewer free points could be obtained in the narrow channel by obstacle point opti-mization technology, leading to significantly reduced sampling probability in the workspace. However, the quasi-random sampling method can guarantee enough sampling points to construct the free space of the normal workspace effectively. (2) Comparing the actual distribution of the sampling points in Figure 27(a1), the external shape of the obstacle had changed. However, the obstacle point optimization technology could still transform the obstacle points into free points, effectively increasing the number of free points near the boundary of narrow channels. It also effectively increased the number of sampling points in free space. Therefore, the PRM algorithm proposed in this paper can still effectively deal with the path planning problem in normal environments. (3) Comparing the initial path and the optimized path in Figure 27(a3,c3), the optimized query process proposed in this paper obtained an initial path far from the obstacle space. Meanwhile, the optimized path pruning technology effectively removed the redundant points of the initial path without changing the distance to the obstacle, thus shortening the length of the path.

Conclusions
In this paper, we proposed a novel PRM which can effectively solve the problem of path planning for mobile robots in a complex environments with narrow channels. First, we expanded the obstacle using the diagonal distance of the robot to ensure the existence of all narrow channels in the working map. In addition" we obtained the sampling point based on the quasi-random sampling principle and reclassified the sampling point type, which can optimize the distribution of random points and ensures the safety of the free points. Second, the obstacle points were transformed from the obstacle space to the narrow passage through use of the potential energy function and the characteristics of the narrow channel, thus improving the density of the sampling points in the narrow passage. Meanwhile, the connection line in the free space was optimized to improve the local motion path obtained by the PRM algorithm in the narrow passage. Third, the bidirectional A* algorithm was chosen, instead of the typical local path planner, in order to reduce the time required by the query process. Furthermore, the potential energy function and path pruning were effectively combined to optimize the initial path obtained by the query process. Finally, we conducted many experiments, in order to demonstrate that the PRM algorithm proposed in this paper can reduce the collision probability of the robot while optimizing the length of the overall path, thus effectively solving the path planning problem for workspaces with narrow environments.
The new PRM algorithm proposed in this paper mainly realizes the path planning problem for mobile robots in global static state. However, the surrounding information of the complex environment typically changes with time in actual workspaces. Therefore, solving the path planning problem for dynamic environments by effectively combining the PRM algorithm and the motion of obstacles is considered a valuable future research direction.

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