Investigating the Impact of Triangle and Quadrangle Mesh Representations on AGV Path Planning for Various Indoor Environments: With or Without Inﬂation

: In a factory with different kinds of spatial atmosphere (warehouses, corridors, small or large workshops with varying sizes of obstacles and distribution patterns), the robot’s generated paths for navigation tasks mainly depend on the representation of that environment. Hence, ﬁnding the best representation for each particular environment is necessary to forge a compromise between length, safety, and complexity of path planning. This paper aims to scrutinize the impact of environment model representation on the performance of an automated guided vehicle (AGV). To do so, a multi-objective cost function, considering the length of the path, its complexity, and minimum distance to obstacles, is defined for a perfect circular robot. Unlike other similar studies, three types of representation, namely quadrangle, irregular triangle, and varying-size irregular triangle, are then utilized to model the environment while applying an inflation layer to the discretized view. Finally, a navigation scenario is tested for different cell decomposition methods and an inflation layer size. The obtained results indicate that a nearly constant coarse size triangular mesh is a good candidate for a fixed-size robot in a non-changing environment. Moreover, the varying size of the triangular mesh and grid cell representations are better choices for factories with changing plans and multi-robot sizes due to the effect of the inflation layer. Based on the definition of a metric, which is a criterion for quantifying the performance of path planning on a representation type, constant or variable size triangle shapes are the only and best candidate for discretization in about 59% of industrial environments. In other cases, both cell types, the square and the triangle, can together be the best representation.


Introduction
Automated guided vehicles (AGVs) have been used for several decades [1,2]. They mostly rely on predefined routes to move from one point to another. However, the flexibility and agility introduced by industry 4.0 [3] are pushing these technologies towards a complete autonomous navigation capability [4]. The required autonomy forces the AGVs to know the environment prior to moving intelligently and safely. Therefore, one of the most important tasks for these vehicles is to model the navigation environment as accurately as possible [5]. In indoor navigation on a planar floor of a factory, a 2D model of the environment, including the specifications of the obstacle-free area, is required to perform path planning. Abstracting, encoding, and saving this type of information is called 'mapping' in mobile robots [6]. A mobile robot uses this map repeatedly to find a path between the current position and the goal position. The essential information for constructing the environment model is usually obtained from raw sensor data with the help of data structures. There are different types of data structures for saving and reading necessary spatial information in navigation. After saving the information, path planning uses this Robotics 2022, 11, 50 2 of 30 information via queries (by means of the data structure) to perform a task. Therefore, the data structure types, spatial models, and the resulting map all play decisive roles in the performance of the global planners.
One of the most well-known mapping methods, which has gained special attention due to its simplicity and robustness, is the occupancy grid technique [7]. This method is composed mainly of uniform grids with simple square shapes, which are made cell by cell, independently. Other shapes such as triangles (with application in polygonal mesh [8], reconstruction [9], etc.) can be used for the discretization of the environment. Squares and triangles are the most famous cell types in a planar application. The analogy between the occupancy grid map and a triangular mesh is the creation of a discretized environment model using the cell concept. Their difference, though, lies on the shape and number of cells, as well as their connection rules (topology) in forming a complete and accurate free space model. By utilizing such a model, a global path planner, e.g., cell decomposition method [10], can be used to find an efficient path between two points.
The efficiency of the path planning algorithm is highly dependent on the environment model, which in turn is impacted by the representation type (cell's shape, number, and topology). For instance, according to the number of cells, it may take a shorter or a longer time to find a path. Moreover, the realized path may have different lengths and distances to the obstacles considering the cell shape and size in the environment. An efficient path should consider all of these essential parameters (such as length, distance, and computational time) [11]. Another aspect worth noting is that, when the conditions of an environment change (new distribution of obstacles), a new representation (discretized model) may lead to a more efficient path. However, investigating the influence of environment model representation on generating a suitable path has escaped the attention of many researchers. Works such as [12,13] have tried to define the environment's difficulty based on parameters such as obstacle shape and distribution or distance to the nearest obstacles. Nevertheless, they did not consider the effect of a navigation scenario (position of start and goal points) and map representation, and they did not propose any solution to consider this difficulty. Recently, valuable work such as that of [14] tried to relate this difficulty to the robot's traversal time in navigation benchmarking. However, this work focused on the overall performance of the navigation platform and does not consider the global planner effect. In a similar task, [14] only considers the effect of cell shape and size on computational efficiency and does not consider their effect on the path length and safety. In addition, ref. [13] defines some metrics and compares different path planning methods only on grid-type cells. Other works such as [15], which uses an adaptive mesh for riskaware path planning, do not mention important parameters such as distance to the nearest obstacle and the effect of robot size. When the environment is cluttered, to our knowledge, there is no straightforward method to select the most appropriate navigation environment approach. For consideration of robot size in path planning, there are various methods. The first strategy applies offset to the boundary of the obstacles, then uses this scaled view as an input for discretization [16]. This strategy is unsuitable when it is necessary to modify the map according to changes in the environment or when the real model of the obstacles is necessary for manipulation. The second strategy is used in the game industry (famous as navigation mesh [17]) and checks each triangle for the input and output edge's length as a possibility of robot passage [18,19]. The third strategy is map generation based on the actual boundary, then applying the inflation layer on the discretized view (the offset or inflation layer is the robot's radius plus any optional value for safety, such as the method implemented in the ROS (Robot operating System) [20]). In addition, it is usual to build an inflated map on top of the metric map in robotic applications [21]. Nevertheless, no specific study has been carried out to select the best environment representation based on the effect of this inflation. Therefore, the main contributions of this paper are:

•
Showing that the type of environment representation has a direct impact on path planning, and quantifying this by a metric.
• Presenting the results showing the pros and cons of each representation when considering the inflation layer By considering the points discussed, this paper explores the effect of the environment representation on the performance of an AGV. In this regard, a multi-objective cost function is developed, and three kinds of cell decomposition (quadrangle, irregular triangle, and varying-size irregular triangle) are used for environment modeling. Furthermore, the robot size is considered through an inflation layer imposed on the model. To the best of the authors' knowledge, this is one of the first attempts to explore the impact of different cell decomposition techniques and robot size on the performance of an AGV. A specific platform is used to generate the three representations with different applicable sizes [22]. Then, a performance comparison is performed for each representation using a fixed scenario obtained by a global planner [13,14].
The rest of the paper is organized as follows. Section 2 depicts the problem formulation and assumptions. Section 3 introduces some common concepts concerning cells and neighbors regarding all the representations. In Section 4, the methodology for all representations is explained by a mathematical formulation. Section 5 discusses all the results, corresponding to the simulation tests on various environments and representations. Section 6 summarizes the outcomes and puts forward future directions.

Problem Formulation and Assumption
An AGV (Automated guided vehicle), denoted by A, navigates in a factory as workspace W s . The workspace contains obstacles, B i , with i ∈ {1, 2, 3, . . . , N o } and N o , the number of obstacles. The obstacle boundaries (Γ(B i )) are generated by a LIDAR (Light Detection and Ranging) placed on the AGV [23]. These can be polygonal or curvilinear curves. The W s is divided into W B , the union of all obstacles, and W free , the free space. S and G are the start point and the goal point of A, respectively, in the W free . According to the distribution of B i s and the size of A, various paths (τ) can be traversed to get to G from S. From an economic point of view, τ must be realized quickly, be as short as possible, and keep a safe distance to B i s. If τ does not satisfy one of these properties, AGV cannot be used reliably for a long time. In this regard, an optimal path (τ opt ) is the one that has the best of the properties mentioned. It is typically determined by a global planner with the help of encoding W free to an environment map (M reps ). The subscript "reps" denotes the encoding types or representations, which have a significant impact on finding τ opt .
This study assumes a perfect circular robot without any sensor noise and environment modeling uncertainty, in order to be able to map spatial information into various representation formats. The index for comparing different representations is a multi-objective cost function (f) which takes the length of τ opt , the complexity of finding the τ opt (correlated with running time), and minimum distance to obstacles into account [24]. Furthermore, no post-processing on the generated path, such as path smoothing [25], has been performed, since this work studies the behavior of various representations on a global planner (as in the deliberative approach concept, in which information is complete before starting the navigation) [26]). The next section discusses common concepts regarding the environment representation types used in mapping and navigation.

Environment Modelling for Navigation
The simplest required information to perform navigation in an environment is the W free . In this regard, one of the best-known methods is cell decomposition, in which W free is converted into digitized cells or a graph of cells. Various shapes can be used for the discretization of the environment. In the planar application, squares and triangles are the most famous cell types. Another practical cell type is the irregular triangular mesh (ITM), used extensively in computer graphics and gaming [17,18] as well as outdoor robotics applications [27]. Using these cells for navigation needs the definition of the neighbor types as either common-edge and common-node neighbors. In a quadrangle (QUAD) cell in a general condition (far from the obstacle), the number of common-edge neighbors and common-node neighbors is four and four. However, in an ITM, these numbers are generally three for common-edge neighbors or a variable number, depending on the degree of nodes.
As is clear, a QUAD cell with the edges parallel to the global frame of an environment cannot be flexible enough to show the boundary of obstacles in the environment except for shapes parallel with the global frame. This limitation causes some inaccuracy in the definition of the area around obstacles and presents limitations to accessing this area. Solving this problem requires a fine square, which, on the other hand, increases the computational burden. An ITM cell, with the flexibility of an edge's direction, can accurately show the boundaries of obstacles. Nevertheless, their irregular shape requires the saving of topological and geometrical information, which makes the map larger. Therefore, each cell shape has its own pros and cons. These features can affect the performance of the algorithm that uses the encoded information of these cells.
The cell decomposition method takes a point as the representative of a cell in a navigation task in order to model the W free . The cell's center can be a simple selection for square cells, in order to calculate properties such as the nearest distance to an obstacle, or as a reference for the position of cells and for moving between cells. After obtaining a representative for each cell, a path can be defined as the lines connecting the start and goal points through the center of cells, if they exist. For a fair comparison between QUAD and ITM, the center of the triangle is selected as the representative of each cell in this study (other methods like that of [28] consider the middle of the edges, with lower computational time but with a longer path). Optimal searching in this graph is performed according to criteria such as the cost and the heuristic of each cell with respect to start and goal points. Since the building block of this graph is a cell, the selection of its size is an effective parameter for the construction and searching of the graph. Small cells give more information than large cells after digitizing, but increase the computational effort of the agent. In addition to cell size, the rules of neighboring, which are different for each cell shape, can change the graph. Hence, finding a good path presents the challenge of finding contradictory parameters. The selection of a suitable cell shape and size can help find a reasonable solution. In the next section, the methodology for studying the effects of these various representations will be described.

Methodology
The performed analysis in this study is divided into two parts. First, the effect of representation is considered regardless of the robot size. Second, the robot size is considered using the concept of inflation layer. This concept generates the map based on the real boundary. Then, it applies an inflation layer to the discretized view. In fact, the offset or inflation layer is the radius of the robot plus any optional value for safety [20]. This strategy has already been implemented in the ROS (Robot Operating System) [20]. In what follows, the decomposition methods are described.

Pure Decomposition
The geometry and the physical size of the robot are not considered in this section. Furthermore, the discretized cell is general and does not depend on the specific type of representation. For navigation on a planar surface by cell decomposition method, it is necessary to encode and abstract the environment information into a map representation, M reps .
where c i denotes the i th cell and implies that the union of all discretized cells (finite number of N) can approximately estimate the free space W free . Each path between S and G points in the map is defined as τ, which is a set of cells: τ = c s , c 1 , c 2 , . . . , c g , g ≤ N Robotics 2022, 11, 50

of 30
A set of cells as a path is divided into several segments S ij between successive cells c i and c j in τ. This leads to: τ = S s1 , S 12 , S 23 , . . . , S ij , . . . , S ig , g ≤ N, i = g − 1 (3) where S ij is generally a curve connecting two cells (c i and c j ) together with all the possible points inside the cells. This segment will be a building unit of movement for AGV according to its controller. So, each movement between c i and c j is a line segment between their centers. However, linear movement between the centers of two neighbors does not always guarantee safety near the obstacles. Therefore, it is necessary to check the feasibility of movement for some cells near obstacles. To do so, Mov c i , c j is defined as: If Mov c i , c j is true, it shows that the movement between the centers of c i and c j is feasible. Otherwise, there is the possibility of collision between agents and obstacles. In the geometrical view, if the combination of neighbor's vertices is convex, then the linear segment will be inside the union of cells, and the movement is safe. Nevertheless, if the combination of c i and c j 's vertices is concave, then the movement is not permitted. This is because it increases the possibility of collision with an obstacle. To include this extra condition into τ: τ = c s , c 1 , c 2 , . . . , c g , ∀ c i and c j =⇒ Mov c i , c j = True (5) Checking the movement feasibility is simple for QUAD cells as the common-edge neighbor produces a convex shape. However, in ITM, it is possible to have a concave shape even for common-edge neighbors near obstacles. Hence, the movement feasibility needs to be checked in this mesh. This extra checking enhances the computational effort of this mesh type [29].
Among all possible τ, finding τ opt is of great interest. According to the discretized nature of the path, it is possible to convert cells and movement between them into vertices and edges of a graph. Then a graph search algorithm can be used to find τ opt between the start and goal cells. A*, as a popular graph search method [30], is employed in this study. This method uses Euclidean distance for calculating the cost and heuristic criterion between start and goal [31]. The cost of cell Among all possible τ, finding τ is of great interest. According to the discretized nature of the path, it is possible to convert cells and movement between them into vertices and edges of a graph. Then a graph search algorithm can be used to find τ between the start and goal cells. A*, as a popular graph search method [30], is employed in this study. This method uses Euclidean distance for calculating the cost and heuristic criterion between start and goal [31]. The cost of cell ℊ(c ) is given by: the nearest obstacle is lower than a specific value. With this definition, the resulting map considering this inflation will be:

given by:
Among all possible τ, finding τ is of great interest. According to the discretized nature of the path, it is possible to convert cells and movement between them into vertices and edges of a graph. Then a graph search algorithm can be used to find τ between the start and goal cells. A*, as a popular graph search method [30], is employed in this study. This method uses Euclidean distance for calculating the cost and heuristic criterion between start and goal [31]. The cost of cell ℊ(c ) is given by: ℊ ( ) = ∥ cent( c ), cent(par(c )) ∥ (6 the nearest obstacle is lower than a specific value. With this definition, the resulting map considering this inflation will be: M = ∀c (i = 1,2, … N); ∑c ≈ W ; cent(c ) ⊄ inflation layer} (16 where cent(c i ) stands for the center of the cell, Symbol a, b is the Euclidean distance between point a and point b, which are the center of two successive cells, and par(c i ) is the parent of c i . This cost starts from S and, in each sequence, accumulates the Euclidean distance between the current cell and its parent in a graph. The heuristic criterion of the The result of A* is a τ with the minimum cost between two cells containing start and goal points. This τ is shown by τ opt : τ opt = c s , c 1 , c 2 , . . . , c g , ∀ c i and c j , where L τ denotes the length of a path. As stated earlier, in this study, the utilized multi-objective cost function is composed of three metrics. These are the complexity of A* algorithm for finding τ opt , length of τ opt , and minimum distance between AGV and obstacles when it follows τ opt . Herein, the complexity (C reps ) is defined by counting the number of visited (In A*, visited cells go to a list named open list, so C reps = size (Open list)) nodes during A*.

C reps = Number of visited cells in A *
The length of the optimal path (L τ opt, reps ) is calculated as follows: The minimum distance to the nearest obstacle is calculated with the help of the Euclidean distance grid, as explained in [32,33].
where DT(p) is a function for mapping each point of free space p into a cell c RG,k in the reference grid (RG), which is a Euclidean distance transform on a 200 × 100 grid. c RG,k contains the value of distance to the nearest obstacle as d k . Therefore, D min,reps , as the minimum distance between all agent's distances to the nearest obstacles, is given by: In triangular representation, as shown in Figure 1, there are some free points near obstacles that are not placed inside any cells in the RG. In this situation, the Euclidean distance between the point and other nodes and edges, which are tagged as the obstacle, is calculated. Then, the minimum value is returned as the minimum distance. This value may be smaller and larger than the previous distance value of the path and must be compared to the current minimum distance.
where L denotes the length of a path. As stated earlier, in this study, the utilized multi-objective cost function is composed of three metrics. These are the complexity of A* algorithm for finding τ , length of τ , and minimum distance between AGV and obstacles when it follows τ .
Herein, the complexity (C ) is defined by counting the number of visited (In A*, visited cells go to a list named open list, so C = size (Open list)) nodes during A*.

C = Number of visited cells in A *
The length of the optimal path (L , ) is calculated as follows: The minimum distance to the nearest obstacle is calculated with the help of the Euclidean distance grid, as explained in [32,33].
where DT(p) is a function for mapping each point of free space p into a cell c , in the reference grid (RG), which is a Euclidean distance transform on a 200×100 grid. c , contains the value of distance to the nearest obstacle as d . Therefore, D , , as the minimum distance between all agent's distances to the nearest obstacles, is given by: In triangular representation, as shown in Figure 1, there are some free points near obstacles that are not placed inside any cells in the RG. In this situation, the Euclidean distance between the point and other nodes and edges, which are tagged as the obstacle, is calculated. Then, the minimum value is returned as the minimum distance. This value may be smaller and larger than the previous distance value of the path and must be compared to the current minimum distance. After calculating the three normalized metrics, a linear weighted sum approach is utilized to define the cost function for each representation. After calculating the three normalized metrics, a linear weighted sum approach is utilized to define the cost function for each representation. where LN reps is the normalized optimal length, CN reps is the normalized complexity for finding optimal length, DN reps is the minimum distance of agent to the nearest obstacle when following the optimal path, L τ opt,reps is the optimal length obtained by (10), C reps is the calculated complexity from (9), D min,reps is the minimum distance of the optimal path from (12), and L max , C max , and D min are the maximum lengths between start and goal, maximum complexity to find this maximum length path, and the minimum possible distance between agent and obstacle, respectively. Finally, the multi-objective cost function (f reps ) is formulated as: where ω 1 , ω 2 and ω 3 are the weights. For simplicity, equal importance is considered for all weights: ω 1 = ω 2 = ω 3 = 0.333.

Decomposition of Inflation Layer
In this sub-section, the robot size is integrated into the inflation layer. The inflation layer is defined as follows: where r AGV is the radius of AGV, and SD is an optional value for a safe distance between the AGV and an obstacle. According to (15), the inflation layer is formed by all the points in W free whose distance to the nearest obstacle is lower than a specific value. With this definition, the resulting map considering this inflation will be: Among all possible τ, finding τ is of great interest. According to the discretized nature of the path, it is possible to convert cells and movement between them into vertices and edges of a graph. Then a graph search algorithm can be used to find τ between the start and goal cells. A*, as a popular graph search method [30], is employed in this study. This method uses Euclidean distance for calculating the cost and heuristic criterion between start and goal [31]. The cost of cell ℊ(c ) is given by: the nearest obstacle is lower than a specific value. With this definition, the resulting map considering this inflation will be:

Results and Discussion
This section presents the results regarding the AGV performance in different conditions. After this, the considered case studies and the utilized platform are first explained. Subsequently, the results concerning environment representation and the inflation layer are discussed.

Platform for Experiment
The intended comparative study is expected to clarify the best representation (M reps ) for each obstacle's distribution (W s ) according to the above-explained multi-objective cost function (f). In this regard, a fixed scenario (a path planning between S and G) is utilized for various M reps and W s . The following types of representation (M reps ) are considered in this study: • QUAD mesh (GRID Nx × Ny): four cell sizes (GRID 100 × 50, GRID 120 × 60, GRID 160 × 80, and GRID 200 × 100) are considered for this mesh where Nx and Ny are the number of cells in horizontal and vertical directions. • ITM: four cell sizes (ITM-700, ITM-500, ITM-300, ITM-100) are also considered for this mesh. ITM-100 refers to a triangle with an average of 100 mm edges. For some environments where this was not applicable, larger triangles and smaller ones were used. • Varying-size ITM (VITM): Concerning this mesh, ITM-100-700 is utilized. This size refers to a triangular mesh with an average edge size of 700 mm, which has been refined to size 100 mm in some places. When confronting narrow passages, VITM-50-700 or VITM-30-700 have also been utilized.
Concerning various factory environments, the following W s are taken into account (see Figure A1): In this study, a research platform developed by our team (ERICCA, Equipe de Recherche en Integration Cao-CAlcul) is utilized to perform the mesh generation for different environments. This platform is based on the Unified Topological Model [22] and automatic meshing and remeshing tools, as explained in [34][35][36][37]. The computer for using this platform and simulation was Intel Core i7-4790 RAM 32GiB. Since the main task is to study the effect of the discretized cell types and sizes on the generated path by the global planner, the required map preprocessing is performed by the operator. For instance, an STP file (the STP is a standard format for exchange of modelling data) is used to define the boundary edges of obstacles and the exterior frame of the environment. This preprocess is similar to changing a point cloud, as sensor data, into a continuous curve of the obstacle's boundary. Using the STP file as an input, various representations (QUAD, ITM, and VITM) with different sizes can be produced by the mentioned platform. A rectangle area of 10,000 × 5000 (mm) is utilized as the environment for all the performed simulations. Various types of obstacle distribution (single and multi-obstacles) are carried out regularly and randomly in this environment. A fixed scenario, navigation between points (300, 300) as the start and (9600, 4400) as the goal, is utilized for all the tests. The reason for selecting this path is that AGV needs to turn around the obstacles to pass the optimal length. After the generation of each map, the UTM platform makes it possible to perform global path planning and find the length, distance, and complexity in order to calculate the cost function for different representations. The running time of the utilized search algorithm is calculated by means of the number of visited cells in the graph search. This strategy is selected in order to be independent of coding details while comparing different cases.

Result-Environment and Representations
Consider an environment with a single medium obstacle as an example, which has been discretized by the three representations (QUAD, ITM, and VITM) with different sizes. Figure 2 represents the optimal path obtained by A* solely for one size of each representation. The scenario is to find a path between the start (300, 300) and goal (9600, 4400) points without considering inflation. Table 1 reports the length, minimum distance to obstacles, complexity, and cost of each case. From this table, ITM-700 has reached the lowest cost followed by VITM-100-700 while the GRID 200 × 100 has led to the highest cost. Therefore, ITM-700 and VITM-100-700 perform better than GRID 200 × 100 in terms of finding an optimal path in an environment with medium-size obstacles in parallel with the global frame. an optimal path in an environment with medium-size obstacles in parallel with the global frame.   Figure A2. The considered environments are the different distributions variously mentioned in Section 5.1. Moreover, different cell sizes of QUAD, ITM, and VITM, as explained earlier, are considered. The complete result is presented in Table A1, and for simplicity, only the best representation of each environment (minimum f values among others) is shown in Table 2. According to this table, a coarse ITM (like ITM-700 or ITM-500) seems to be a good candidate for almost all environments except the hybrid environment where a VITM (VITM-100-700) shows better performance. It is worth reminding the reader that the detailed results, as in Table 1, are presented in the Appendix for all   Figure A2. The considered environments are the different distributions variously mentioned in Section 5.1. Moreover, different cell sizes of QUAD, ITM, and VITM, as explained earlier, are considered. The complete result is presented in Table A1, and for simplicity, only the best representation of each environment (minimum f values among others) is shown in Table 2. According to this table, a coarse ITM (like ITM-700 or ITM-500) seems to be a good candidate for almost all environments except the hybrid environment where a VITM (VITM-100-700) shows better performance. It is worth reminding the reader that the detailed results, as in Table 1, are presented in the Appendix for all environments. The discussed results have been obtained by considering the pure decomposition effect without an inflation layer. In the next sub-section, the results concerning the effect of robot size or inflation layer are explored.

Result-Inflation Effect on the Representations
In practice, it is necessary to consider the robot size in path planning. This is normally done by imposing an inflation layer around the obstacles. In this regard, the tests performed in the previous sub-section are repeated herein, considering the inflation layer.
Firstly, the test regarding the environment with medium-size obstacles is repeated by considering a robot size diameter of 354 mm (turtle bot 2 [38]). Figure 3 shows the obtained optimal path for this new case. The red layer around the obstacle and exterior boundary is not accessible by the robot. Therefore, the movements between the neighbors are affected by this layer, and the results will be changed. The new values considering this limitation are given in Table 3 for the three different representations. From this table, VITM-100-700 reaches the best performance ( f reps : 0.009) compared to the other two cases. This shows that considering the inflation layer influences the choice of the best representation, due to the change in complexity.
Secondly, all the previously considered environments of different types and sizes are repeated, considering the inflation layer. The obtained results are shown in Table 4. In this new analysis, as opposed to the case in which inflation was ignored, it is seen that many ITM representations have not been able to find the path, and they have been replaced by the VITM and QUAD representations. Furthermore, it has been necessary for the VITM to use finer triangles in the narrow passages that have not been blocked completely. Decreasing the size of triangles (increasing their number) can increase the chance of finding more triangles whose centers are inside the free space, and movement between them and their children is possible. Some of the generated paths are presented in Figure 4. According to Table 4, coarse ITM shows the best performance for a single obstacle or hybrid environment, with enough wide space (for Map_1, Map_3, and Map_14). Regarding the corridor and regular multi-obstacles, the coarse grid is the best method. In irregular multi-obstacles and a curvilinear passage, a VITM indicates the best performance. Robotics 2022, 11, x FOR PEER REVIEW 11 of 25 Secondly, all the previously considered environments of different types and sizes are repeated, considering the inflation layer. The obtained results are shown in Table 4. In this new analysis, as opposed to the case in which inflation was ignored, it is seen that many ITM representations have not been able to find the path, and they have been replaced by the VITM and QUAD representations. Furthermore, it has been necessary for the VITM to use finer triangles in the narrow passages that have not been blocked completely. Decreasing the size of triangles (increasing their number) can increase the chance of finding more triangles whose centers are inside the free space, and movement between them and their children is possible. Some of the generated paths are presented in Figure 4. According to Table 4, coarse ITM shows the best performance for a single obstacle or hybrid environment, with enough wide space (for Map_1, Map_3, and Map_14). Regarding the corridor and regular multi-obstacles, the coarse grid is the best method. In irregular multi-obstacles and a curvilinear passage, a VITM indicates the best performance.  To find the effect of inflation on each type of representation, Table 5 shows the average cost function (f values) for all environments in two conditions, without inflation and with an inflation layer. Their ratio shows that ITM is more sensitive to the inflation layer, and misses its benefit. QUAD representation shows better behavior when using inflation. To compare the complexity of finding a path for each type of representation, we have calculated the average complexity (number of visited cells when graph searching) of each type of representation in Table 6. If any size could not give a path, a smaller size has been used to produce one. Overall, VITM shows the lowest complexity as compared to other types.

Results-Discussion
Finding the best representation for an AGV inside a factory according to different environment encounters during its navigation is an important part of the use of autonomous agents. Different representations give different path. A good path has a short length with a maximum distance to obstacles that can be calculated quickly. Without knowing this path, a mobile robot can be an unfeasible solution for a factory, along with increasing energy consumption due to a less than optimal path, lowering efficiency by inappropriate planning time and increased cost due to the high risk of an inadequate distance to obstacles (high probability of collision). By considering all these important parameters of path planning, this work defined a single criterion and checked it for different types and sizes of representation for the various environments inside a factory.
The obtained results show that for various types of environments, such as single or multi obstacles and corridor-like, some level of coarse ITM is the predominant strategy for achieving a compromise solution for length, safety, and computation time; this verifies the work in [14], which shows that refinement is not necessary for some motion planning. However, what has not been considered is that, when considering the inflation layer around obstacles, VITM and GRID representations are better. It is clear that, by increasing the resolution or decreasing the cell size, the complexity increases, and in most cases, the nearest distance to the obstacles decreases, with some average improvement in length. According to this result, when the robot is of a fixed size in a static environment, and it is possible to apply first the robot size to the map, it is better to use the ITM method with coarse size like [16] in which inflation has been applied as offset during the first mapping. However, when there are the following situations, it is better to use VITM or GRID representations:

•
In multi-robot sizes conditions, when is not possible to apply the offsetting method. • Some changes in the environment make it necessary to apply online inflation to the map, and not by offsetting the obstacles in the primary map.

•
There are narrow passages for a robot which makes ITM representation useless.
This benchmark also demonstrated that varying the size of the triangular mesh with a different strategy for refinement could be a competitive method when ITM cannot find the answer. Of course, this needs to be more refined when using the elements in an area with narrow passages. The other benefit of VITM is the more straightforward implementation relative to the quadtree method for neighbor finding, such as regular triangulation.
Robotic navigation in a real factory with multi AGVs with non-circular shapes is a more complicated task. This paper surveyed the effect of environment and map representation types on path planning considering a single circular AGV. Although a circle can cover any shape as a peripheral circle, it cannot obtain a path in some narrow passages. According to the robot's footprint, this limitation leads to consideration of a more complex definition of the inflation layer. In addition, to show the clear behavior of various cell shapes in response to A* search graph, we did not consider any path smoothing, which causes a non-optimal path from the point of energy consumption.

Conclusions and Future Work
In using an AGV as a feasible solution in a factory in response to the fourth industrial revolution, it is necessary to consider a flexible multi-objective path planner to give better answers according to the changes in the environment. Hence, this paper investigates the impact of two important factors, environment model representation and robot size, on the performance of an AGV. In this regard, a multi-objective cost function, considering length, complexity, and minimum distance to obstacles, is formulated for a perfect circular robot with a fixed path. Subsequently, three types of representation, namely QUAD, ITM, and VITM, with different cell sizes are employed to model the environment. Moreover, an inflation layer is imposed on the discretized view to check its influence on the performance. A* is utilized to find the optimal path for the different cell representations and sizes with and without the inflation layer. The obtained results show that a coarse ITM is a good candidate for a relatively unchanging environment with fixed robot size and wide space. However, when it is possible to repeatedly change the factory's plan and use different 'robot sizes, VITM and GRID representations give better performance by applying an inflation layer, especially in narrow spaces. In comparing triangle and square representations, in about 59% of industrial environments, constant or variable size triangle shapes are the only best for the discretization of that specific environment.
Future work can pursue the following paths, which will bring important practical application in the field of AGV environment modeling and navigation. Consideration of dynamic environments will enrich the optimization function with more complex criteria using a solution such as right-hand traffic, and applying some simultaneous constraints on the navigation of a narrow passage. Furthermore, the results of this paper can be used as an efficient map server that selects the best representation according to the environment. This can bring a smarter AGV platform to industry.