Next Article in Journal
Genomic Association Mapping of Apparent Amylose and Protein Concentration in Milled Rice
Previous Article in Journal
Morphological and Yield Parameters, Dry Matter Distribution, Nutrients Uptake, and Distribution in Strawberry (Fragaria × ananassa Duch.) cv. ‘Elsanta’ as Influenced by Spent Mushroom Substrates and Planting Seasons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Path Optimization Algorithm for Multiple Unmanned Tractors in Peach Orchard Management

1
National Engineering Research Center for Information Technology in Agriculture, Beijing 100097, China
2
Information Technology Research Center, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
3
School of Computer, Electronics and Information, Guangxi University, Nanning 530004, China
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(4), 856; https://doi.org/10.3390/agronomy12040856
Submission received: 16 February 2022 / Revised: 28 March 2022 / Accepted: 28 March 2022 / Published: 30 March 2022

Abstract

:
In order to improve the management efficiency of peach orchards, this paper considers the cooperative operation scheme of multiple unmanned tractors. According to the actual situation, this paper constructs the path planning model of multiple unmanned tractors in a standard peach orchard, designs the objective function to optimize the total turning time and total operating time according to the tractor driving parameters, and solves it by improving the differential evolution algorithm. Aiming at the premature convergence problem, the permutation matrix is introduced to represent the driving paths of multiple unmanned tractors. Then, the dynamic parameters are adopted to make the parameters change with the number of iterations, and the elite selection strategy is used to eliminate the redundant feasible solutions. An Adaptive Elite Differential Evolution (AEDE) algorithm suitable for multi-tractor path optimization is proposed. The results show that, compared with the traditional Differential Evolution algorithm (Differential Evolution, DE), the total turning time and total operating time in the rectangular peach orchard optimized by AEDE are reduced by 3.34% and 0.87%, respectively. Compared with the block operation, the total turning time and total operating time of the AEDE-optimized rectangular peach orchard operation path were reduced by 37.37% and 9.47%, respectively. Experiments show that AEDE, which optimizes the operating path of multi-tractors in standard peach orchards, is able to improve the efficiency and reduce the operating time.

1. Introduction

The cultivation area of peaches in China is 55,317.23 (km2), and the planting scale ranks first in the world [1]. Peach orchard management includes weeding, ditching, fertilizing, and soil cultivation, which requires a large amount of labor, contradicting the increasingly serious aging of the rural labor force. At present, China is vigorously promoting the integration of agricultural machinery and agronomy and standardized farmland transformation. The integrated path planning systems of agricultural unmanned tractors can perceive and act independently [2]. By connecting different agricultural tools to carry out different operations in standard peach orchards, the labor force can be reduced and the planting management efficiency can be improved. Optimizing the driving path is a key to improving unmanned production efficiency and operation quality [3].
Spekken et al. mentioned that during agricultural machinery operation, most of the time was used for non-productive operations, such as turning in non-working areas. Reasonable path planning can reduce the working time of agricultural machinery and improve operation efficiency [4]. At present, the B-Pattern model [5] proposed by Bochtis was adopted as the path planning model in most studies. Based on the B-Pattern model, the tractor operation scheduling is optimized by intelligent optimization methods, such as genetic algorithm [6,7,8], simulated annealing algorithm [9,10], ant colony algorithm [11,12], domain search [13], tabu search algorithm [14,15], particle swarm optimization algorithm [16], and Clark–Wright algorithm [17,18]. In addition, there is research on operation path optimization based on probability. Utama et al. introduced domain search and elite selection to improve the effect of distribution estimation algorithms in agricultural machinery operation path planning [19].
Fabre et al. [20] transformed the path planning problem of agricultural machinery into a Hamiltonian problem and proposed a path optimization algorithm based on a greedy strategy to minimize the repetitive operation area. Jensen et al. transformed the path planning problem of agricultural machinery into a Traveling Salesman Problem (TSP) problem and proposed an algorithm based on the state-space search to reduce the length of non-operational paths [21]. Muthukumaran et al. adopted the hybrid dragonfly cuckoo search algorithm to optimize the path coverage based on the TSP model [22], and the computational efficiency and effect were better than most intelligent optimization algorithms. Postmes transformed the operation coverage optimization of agricultural machinery robots into VRP (Vehicle Routing Problem) and tried to solve the problem with various approximation algorithms [23]. In order to reduce the complexity of agricultural machinery path planning, Oksanenand et al. used trapezoidal decomposition and an incremental method to simplify the path coverage problem [24].
At present, the studies on autonomous driving path planning in farmland mainly focus on food production, including global path planning and local path planning. The traditional B-Pattern planning model is mainly used in the field where the spacing of operation lines is equal (operation line refers to the route that the tractor starts from one end of the orchard and ends at the opposite end, such as transplanting and harvesting in the field). When the unmanned driving system is applied to the management of a peach or-chard (such as ditching, weeding, fertilization, soil tillage), the operation is usually one-sided, resulting in uneven spacing between operation lines. Therefore, the path autonomous decision-making in the peach orchards cannot be directly applied to the B-Pattern model [25].
The spacing of operation lines in the peach orchard is unequal, and peach trees are densely distributed. This paper studies a multi-machine cooperative operation path optimization method that is suitable for standard peach orchard management, which provides a reference for the continuous and stable operation of unmanned systems.

2. Path Optimization Model of Unmanned Tractor

The distribution of peach trees in the standard peach orchard is shown in Figure 1. The peach orchard is divided into Upper headland (UH), working area, and Lower headland (LH). The working area is the area where the tractors manage the peach trees, and the headland is the area where the tractor turns. Suppose there are n rows of peach trees in the peach orchard—vectors xk and yk represent the x-coordinate and y-coordinates of each peach tree in row k, respectively. Additionally, xmk represents the x-coordinate of the peach tree m in row k (counted from bottom to top, consistent with the positive direction of Y-axis), where 1 ≤ kn, 1 ≤ mnk, nk represents the number of peach trees in row k. The coordinates of peach trees at both ends of each row are represented by n × 2-dimensional matrices U and L, and Uk,1 and Uk,2 respectively represent the horizontal and vertical coordinates of peach trees at the junction of headland UH and working area in row k, i.e., x n k k and y n k k . Similarly, Lk,1 and Lk,2 respectively represent the horizontal and vertical coordinates of the peach trees in row k at the junction of the headland LH and working area, i.e., x 1 k and x 1 k . The boundary of the headland and the working area in the standard peach orchard can be determined according to the matrices U and L.
In order to ensure that peach trees have enough space for growth and management, the spacing of peach trees is much larger than the width of tractors. As shown in Figure 1, the tractor needs to be driven at least twice to manage a row of peach trees, and probably does not work in sequence on both sides of the peach trees. Therefore, when the tractor works in the standard peach orchard, the driving path is not the equidistant operation line in the traditional B-Pattern model. If the peach orchard has n rows of peach trees, the tractor needs to make at least 2n to complete the management of all peach trees. The dashed line in Figure 1 represents the track of the tractor when it is working in the standard peach orchard. Although it is also composed of several parallel lines, they are not equidistant from each other.
Unmanned tractors need to make turns in the field before working on the next row of peach trees. There are three main types of turns, including the Π-Turn, Ω-turn, and T-turn. As the T-turn requires the tractor to reverse and stop, it will not only consume time, but also aggravate the tractor wear. This paper considers only the Π-Turn and Ω-turn, as shown in Figure 2a,b. In Figure 2, w and d represent the horizontal and vertical distance of two adjacent rows of peach trees, respectively, and r represents the turning radius of the tractor.
Based on the above path planning model, this paper designs a path optimization scheme for multiple unmanned tractors. The design of a multi-tractor driving route can avoid incomplete operation or repeated operation. The simplest scheme is to divide the peach orchard into several sub areas, and then assign different tractors to operate in the divided areas. The schematic diagram of block operation is shown in Figure 3.
Block operation is the operation mode of an S-type operation for multiple tractors, which requires multiple turns and has low operation efficiency. Reasonably arranging the operation sequence of the tractor and reducing the Ω-turn is conducive to improving the overall operation efficiency.

Optimization Model of Minimum Operating Time for Unmanned Tractor

The term “operation line” refers to the route that the tractor starts from one end of the orchard and ends at the opposite end. E represents the set of all operation lines, V is the set of end points of all operation lines, and the driving path is represented by a directed graph G = {V,E}. Graph G is mainly composed of vertices and connecting edges of parallel lines, so it is connected. E v i , v j represents the edge from the vertex vi to vj of the operation line, where vi, vjV. R represents the set of tractor paths, and Rk represents the driving path of the k-th tractor. Rk is composed of a vertex set <R1k, R2k, …, R m k k >, where mk represents the number of vertices of path Rk.
To improve the management efficiency of the standard peach orchard, it is considered to formulate a cooperative scheme for multiple tractors to ensure the shortest working time. As shown in Equation (1), this paper establishes an objective function to minimize the total turning time and working time of multiple tractors, which ignores the time of tractors entering and leaving the orchard.
t a l l = min t t u r n , t w o r k z t t u r n + ( 1 z ) t w o r k M
where,
tturn—total turning time of multiple unmanned tractors
twork—total operating time of unmanned tractors
M—the number of unmanned tractors
tall—weighted average of total turning time and total operating time
z—the weight variable represents farmers’ attention to the total turning time and the total operating time
tturn can be calculated according to Equation (2):
t t u r n = k = 1 M i = 1 | R k | 1 c R i k ,   R i + 1 k D
where,
c R i k ,   R i + 1 k D (D ∈ {UH,LH})—turning time for unmanned tractor to travel from Rik to Ri+1k,
c R i k ,   R i + 1 k UH —the time cost of changing an unmanned tractor from Rik to Ri+1k in headland UH.
c R i k ,   R i + 1 k LH —the time cost of changing an unmanned tractor from Rik to Ri+1k in headland LH.
twork can be calculated according to Equation (3):
t w o r k = max k ( i = 1 | R k | c R i k , R i k W + i = 1 | R k | 1 c R i k ,   R i + 1 k D )
where, c R i k , R i k W —operating time of tractor in the i-th operation line, which is calculated by the operation line length and tractor speed.
Assuming that only one tractor is allowed to work in each operation line, each operation line is traversed once, and constraints (4) and (5) can be obtained. Under the constraint of this model, the number of operation lines is 2n, and the number of turning paths is 2n-M.
m = 1 M R m = 0
E m = 1 M R m
When an unmanned tractor enters the next operation line from one operation line, Ci,j describes the time cost of the tractor from line i to line j, and its definition is shown in Equations (6) and (7).
C i , j = ( c i . j UH c i . j LH )   ( i j )
C i , j = ( c i . j W c i . j W )   ( i = j )
When i = j, Ci,j represents the work time in the i-th operation line. When ij, Ci,j means turning from the i-th operation line to the j-th operation line. At this time, there are two possibilities: ① the tractor changes the operation line at the headland UH; ② the tractor converts the operation line in the headland LH. Considering that peach trees may be asymmetrically distributed on both sides of the headland, the time cost of changing the operation line at headland UH may be different from that at headland LH; that is, ci,jUH may not be equal to ci,jLH. Based on the working time of the tractor in the operation line and the turning time of the tractor in the headland, the matrix C′ representing the time cost of the tractor changing the operation line is constructed in this paper:
C = ( c 1 , 1 W c 1 , 1 W c 1 , 2 UH c 1 , 2 LH c 1 , | E | UH c 1 , | E | LH c 2 , 1 UH c 2 , 1 LH c 2 , 2 W c 2 , 2 W c 2 , | E | UH c 2 , | E | LH c | E | , 1 UH c | E | , 1 LH c | E | , 2 UH c | E | , 2 LH c | E | , | E | W c | E | , | E | W )  
The calculation of ci,jUH and ci,jLH is related to the turning model used in path planning. In order to reduce turning time and improve operation efficiency, this paper adopts the Π-Turn model and Ω-Turn model. Assuming that ci,jD, D ∈ {UH,LH} represents the time cost of changing from line i to line j at the headland, and the time cost of the Π-Turn model and Ω-Turn model can be calculated by Equations (9) and (10), respectively.
Π i , j D = π r v t + ( d i , j D ) 2 + ( w i , j D 2 r ) 2 v s
Ω i , j D = ( 3 π a cos ( 1 ( d i , j D ) 2 + ( w i , j D 2 r ) 2 8 r 2 ) ) r v t
where,
di,jD—vertical coordinate difference between the i-th operation line and the j-th operation line at the endpoint near the headland D.
wi,jD—abscissa coordinate difference between the i-th operation line and the j-th operation line at the endpoint near the headland D.
vt—the turning speed of the tractor at the end of the headland.
vs—the straight speed of the tractor at the end of the headland.
di,jD can be calculated according to di,jD = |XiDXjD|. XiD and XjD represent the x-coordinate of the j-th operation line. After calculating all Uk,1 ± W(1 ≤ kn), XiU and XjU are sorted, and similarly, XiL is obtained, where W is the working width of the tractor.
According to W and r, ci,jD can be obtained:
c i . j D = { Π i , j D w i , j D 2 r Ω i , j D w i , j D < 2 r
Based on the number |E| of operation lines, the number M of unmanned tractors, and the cost matrix C′, the objective function is obtained as shown in Equation (1), where the set of traveled operation lines by each tractor is used as the variable.
A reasonable path optimization scheme for multi-unmanned tractors can reduce the total operating time and improve work efficiency [26]. However, this problem belongs to the TSP [27,28], and it is difficult to calculate the optimal solution in a short time. Differential evolution algorithm is an intelligent optimization algorithm based on population evolution proposed by Storn [29,30]. It has the characteristics of fewer operators and strong search ability. The solution can be obtained through mutation, crossover, and selection to update the population [31,32,33]. This paper designs a differential evolution algorithm to obtain multi-tractor paths.

3. Multi-Tractor Path Optimization Algorithm

3.1. Differential Evolution Algorithm

Differential evolution algorithm needs to balance the global search and local search in the iterative process [34]. A global search is to search different spaces that have not appeared. A local search means the ability to improve the solution quality according to the information of the current solution. Differential evolution algorithms can obtain effective solutions through mutation and crossover operations. However, when differential evolution algorithm pays attention to local search, it is easy to fall into local optimization, resulting in premature convergence. On the other hand, if the differential evolution algorithm focuses on global search, the convergence will be very slow [35].
The mutation is also called mutation operator, which calculates the position of the new individual according to the individual position of the last iteration and the optimal position of all individuals. Differential evolution algorithm mainly has six mutation operators, and Equation (12) gives a common mutation operator [36].
κ i t + 1 = κ i t + F ( κ b e s t t κ i t ) + F ( κ r 1 t κ r 2 t )
where,
t—iteration number.
i—individual ID in algorithm population.
κit—the position of the i-th individual in the t-th iteration, |E|-dimensional column vector.
κbestt—the individual position with the best fitness in the t-th iteration, |E|-dimensional column vector.
F—a scaling factor that controls the scale of the difference.
The crossover of differential evolution algorithm is shown in Equation (13), where there is a random variable subject to the uniform distribution. When the random variable is less than the crossover probability, accept the mutation operation, otherwise, keep the original value.
κ i , j t + 1 = { κ i , j t + 1 r a n d < C R κ i , j t otherwise
where,
CR—crossover probability.
rand—random number obeying uniform distribution [0,1).
The selection can be represented by Equation (14). If the fitness of the new solution is less than that of the original solution, the new solution is accepted, otherwise the original solution is retained.
κ i t + 1 = { κ i t + 1 f ( z i t + 1 ) < f ( z i t ) κ i t otherwise
where, f(κ)—objective function.
Although the differential evolution algorithm can obtain the minimum driving cost solution of the unmanned tractors in a limited time, it has the problem of premature convergence and insufficient exploration of solution space. In addition, the mutation operator generates a sequence of real numbers, not an integer sequence in a given range. Multi-tractor cooperation requires a finite number of integers to indicate the number of unmanned tractors and the current scheduling state. The traditional differential evolution algorithm cannot be directly used to solve the multi-tractor path optimization problem. When solving this kind of problem, they usually optimize the continuous variables first, and then indirectly obtain the discrete operation line travelling scheme through the round rather than directly optimizing the discrete operation line travelling scheme. Therefore, this paper provides a better balance between global search and local exploration by improving the mutation operator and crossover operator to optimize the differential evolution algorithm so as to obtain an efficient multi-tractor cooperative path optimization scheme.

3.2. Improved Differential Evolution Algorithm

3.2.1. Mutation Operator

A permutation matrix is not only an invertible matrix but also an orthogonal matrix. It can greatly reduce the convergence time by transposing when finding the inverse. In the cooperative path optimization problem for multi-tractors, the tractor sequence κi is one of the permutations composed of 1~|E|, and κi is a column vector. According to the definition of the permutation matrix, the column vector xi obtained by left multiplying a permutation matrix is still one of the permutations composed of 1~|E|; that is, it is still a legal operation sequence after the transformation of the permutation matrix. Therefore, this paper introduces a permutation matrix to optimize the mutation operator.
Let the position κ of the i-th individual in the t iteration be:
κ i t = P i t α
where,
Pit—the permutation matrix of the individual i in the t-th iteration, with the size of |E|×|E|.
α—|E|-dimensional column vector, the expression is [1 2 … |E|]T.
Substitute Equation (15) into mutation operator (Equation (12)). Since α is a non-zero vector, the simplified operator (16) is obtained:
P i t + 1 = P i t + F ( P b e s t t P i t ) + F ( P r 1 t P r 2 t )
However, when the permutation matrices are added or subtracted from each other, it is difficult to ensure the orthogonality of the matrix, so it can not be calculated directly by Equation (16).
Two permutation matrices are still permutation matrices after point multiplication, and the inverse cost is very low. Point multiplication can be used instead of plus; dot inverse matrix can be used instead of subtraction; let F = 1; Equation (16) can be rewritten as:
P i t + 1 = P i t ( P b e s t t ( P i t ) 1 ) ( P r 1 t ( P r 2 t ) 1 )
However, in Equation (17), the time complexity of calculating the inverse matrix is at least the cubic, and the calculation cost is large, which is not conducive to the fast acquisition of planning paths. The inverse operation of the orthogonal matrix is equivalent to the transpose operation, and the transpose operation is simpler. Therefore, this paper rewrites Equation (17) into Equation (18) with clearer mathematical meaning and less calculation:
P i t + 1 = P i t ( P r 1 t ( P r 2 t ) T ) ( P b e s t t ( P i t ) T )
This paper takes Equation (18) as the mutation operator of the improved differential evolution algorithm so as to optimize the operation path of the unmanned tractors.

3.2.2. Dynamic Parameters

In the traditional differential evolution algorithm, the crossover probability CR and scaling factor F are fixed and do not change with iteration. CR affects the probability of mutation. The smaller its value, the smaller the probability of mutation. Currently, the local search ability of the algorithm is strong. On the contrary, the greater the mutation probability, the stronger the global search ability of the algorithm. If F is fixed, the global search and local search cannot change adaptively, resulting in premature convergence and low search accuracy. As shown in Equation (19), the dynamic parameter scheme is adopted in this paper to make CR changes with the number of iterations. n is the number of operation lines, t is the current number of iterations, and T is the maximum number of iterations.
C R = 1 n + ( 1 1 n ) × ( 1 t T )
When t increases from 1 to T, the crossover probability CR gradually decreases to 1/n. The global search of the algorithm changes from strong to weak, and the local search changes from weak to strong. This can avoid premature convergence and improve search accuracy.
Similarly, too small a scaling factor F will also affect the search accuracy and convergence speed. This paper designs an adaptive function as shown in Equation (20). When t increases from 1 to T, the scaling factor F gradually increases to 1. The global search of the mutation operator changes from strong to weak, and the local search changes from weak to strong, which promotes the convergence of the algorithm.
F = 0.2 + 0.8 t T

3.2.3. Elite Selection

Elite selection can improve the convergence speed and local search performance. It is assumed that κn−2t, κn−1t, κnt are the first three individuals with the worst fitness in the t-th iteration and κ1t−1, κ2t−1, κ3t−1 are the first three individuals with the best fitness in the (t − 1)-th iteration. κn−2t, κn−1t, κnt can be replaced by κ1t−1, κ2t−1, κ3t−1 to realize the elite selection.
This paper introduces the permutation matrix, adopts dynamic parameters and elite selection strategy to improve the differential evolution algorithm, and proposes an Adaptive Elite Differential Evolution (AEDE) Algorithm. The algorithm flow chart is shown in Figure 4.
The multi-tractor cooperative path optimization method studied in this paper firstly generates several feasible driving paths randomly according to the input orchard information, tractor parameters, and the number of tractors. Then, the crossover probability CR and scaling factor F of the current iteration are calculated, and the mutation operator and crossover operator designed in this paper are used to generate a new multi-tractor cooperative operation path optimization. Then, the solution with the worst fitness is eliminated through the elite selection, and the three schemes with the best fitness in the last iteration are introduced. The scheme with better fitness is selected through the selection operation to replace the current scheme until the termination condition is reached. This scheme with the best current fitness is output as the optimal multi machine cooperative path planning scheme solved by the algorithm.

4. Experiment and Analysis

This paper took the John Deere 6B-1204 tractor commonly used in planting as the research object, and the total number of tractors was six. The working width of each tractor was W = 3.5 m, the minimum turning radius was rmin = 5.5 m, the working speed was vw = 1 m·s−1 in the working area, the straight speed was vs = 5 m·s−1 in the non-working area, and the turning speed was vt = 2 m·s−1. The peach orchard was rectangular, and the row spacing of peach trees was 4 × 5 m. As the peach tree distribution will affect the algorithm performance, this paper tests the path planning effect of the proposed model and the improved algorithm with different peach tree distributions. Figure 5 shows the distribution of peach trees in three rectangular standard peach orchards. Figure 5a shows the peach orchards with the rectangular distribution of peach trees. There were 20 rows of peach trees in the peach orchard, with 15 trees in each row. Figure 5b shows peach orchards with the approximately trapezoidal distribution of the peach trees, and Figure 5c shows peach orchards with an irregular distribution of peach trees. The experiment assumes that the tractor works in a peach orchard without obstacles and the driving process is continuous.
The parameters of the original differential evolution algorithm DE are as follows: population size N = 100, crossover probability CR = 0.1, scaling coefficient F = 0.5, and the maximum number of iterations is 500. Equation (1) is used as the fitness function. The experimental environment is Windows 10 and Python 3.7.

4.1. Analysis of Experimental Results

The evaluation indexes of the model adopt the total turning time decrease rate (TTTDR) and total operating time decrease rate (TOTDR). The higher the value, the better the optimization effect of the algorithm. TTTDRn represents the reduction rate of total turning time between AEDE optimized and control group when n tractors work together. Similarly, TOTDRn represents the reduction rate of operation time after AEDE optimization and comparison algorithms when n unmanned tractors work together. The calculation formulas of TTTDRn and TOTDRn are shown in Equations (21) and (22).
T T T D R n = T T T O n T T T I n T T T O n × 100 %
T O T D R n = T O T O n T O T I n T O T O n × 100 %
In addition, we also use the following indexes to evaluate the algorithm performance. TTTOn is the total turning time of the n combined after control group optimization. TTTIn is the total turning time of the n combined after AEDE optimization. TOTOn is the total operating time of the n combined after control group optimization. TOTIn is the total operating time of the n combined after AEDE optimization.

4.1.1. Performance Comparison of Elite Selection Algorithm

Using traditional DE as the control group, we compared it with the AEDE proposed in this paper in terms of total turning time and total operating time, respectively. Both DE and AEDE are intelligent evolutionary algorithms; different peach tree distributions do not affect the solution process. Therefore, this comparison is based on the rectangular distribution of peach trees. The results are shown in Table 1. According to Table 1, the total turning time of AEDE method proposed in this paper was reduced by 3.34%, and the total operating time was reduced by 0.87%. It shows that after improving DE, the total turning time and total operating time are optimized to a certain extent.

4.1.2. Comparison with Block Operation

Block is a common multi-tractor scheduling scheme, which divides orchards into different areas and is operated by a single tractor. Using block operation as the control group, this paper compares it with AEDE in terms of total turning time and total operating time, respectively. The results are shown in Table 2.
From Table 2, the total turning time of AEDE after optimization is reduced by more than 30%, and the total operating time is reduced by about 10%. This is because when working in block mode, the tractor needs to make frequent Ω-turns on the headland, which increases a lot of non-working time compared with the Π-Turn. AEDE avoids a lot of turning by optimizing the driving path of multiple tractors so as to reduce the total turning time and total operating time. In addition, AEDE can get good results whether the shape of the peach orchard is regular or irregular, which shows that AEDE is not sensitive to peach tree distribution and can be used to optimize peach orchard with complex peach tree distribution.

4.1.3. Comparison of Effective Operation Ability

Effective operation and average operation are important indexes to evaluate the cooperative operation efficiency of multiple tractors [37]. In a standard peach orchard, the effective operation refers to the ratio of the peach orchard area to the total operating time, and the average operation refers to the ratio of the effective operating capacity to the number of tractors. Figure 6 shows the changes of effective operating capacity and average effective operating capacity after block and AEDE optimization when the number of tractors working in rectangular orchards increases from two to six.
In Figure 7, the effective operating capacity shows an increasing trend because increasing the number of tractors reduces the total operating time in the orchard. AEDE reduces the total turning time by optimizing the turning path of the tractor, which in turn reduces the total operating time. When the number of tractor operations is two, the effective operation capacity and average effective operation capacity calculated by block operation are 0.0124 km2·h−1 and 0.0062 km2·h−1, respectively. The effective operation capacity and average effective operation capacity calculated by AEDE are 0.0138 km2·h−1 and 0.0069 km·h−1, respectively, which are 11.5% and 11.5% higher than block, respectively. When the number of tractors is six, the effective operation capacity and average effective operation capacity calculated by AEDE are increased by 8.6% and 8.6%, respectively, compared with block.

5. Conclusions

In order to improve the operating efficiency of tractors in standard peach orchards, this paper studies a path optimization algorithm based on differential evolution for multiple tractors.
(1) In the rectangular peach orchard, the total turning time and total operating time in the process of management operation are reduced by 3.34% and 0.87%, respectively. Compared with the non-optimized differential evolution algorithm, it shows that the performance of the algorithm studied in this paper is better.
(2) Compared with block operation, the total turning time and total operating time in the process of management operation obtained by this algorithm are reduced by 37.37% and 9.47%, and the effect is better. It not only improves the operation efficiency of the tractors in the standard peach orchard, but also reduces the operating loss of the tractors.
However, when the tractors meet in the turning area, they will automatically avoid obstacles by adjusting the path. As the adjusted path is shorter than the length of the turning path and working path, the mutual interference of tractors in the turning area is ignored in this paper. The time consumption of tractors entering and leaving the peach orchard is also not considered in the model. In further work, we plan to use the time window to construct the turning constraint so as to avoid the situation that tractors meet in the turning area. At the same time, a more accurate path optimization method for multi-tractor cooperative work is realized by setting fixed starting points and ending points of tractors.

Author Contributions

Literature search, figures, tables, and data collection, X.H.; conceptualization, methodology, and investigation H.W. and Y.L.; validation, H.W., Y.L. and X.H.; formal analysis, X.H.; writing—original draft preparation, all authors; funding acquisition, H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant number 2020YFD1100602, the China Agriculture Research System of MOF and MARA, grant number CARS-23-D07 and the Beijing Municipal Science and Technology Planning Project, grant number Z191100004019007.

Data Availability Statement

Not applicable.

Acknowledgments

The author would like to thank anonymous commentators and editors for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gao, G.; Xiao, K.; Jia, Y. A spraying path planning algorithm based on colour-depth fusion segmentation in peach orchards. Comput. Electron. Agric. 2020, 173, 105412. [Google Scholar] [CrossRef]
  2. Seyyedhasani, H. Using the Vehicle Routing Problem (VRP) to Provide Logistics Solutions in Agriculture. Ph.D. Thesis, University of Kentucky, Lexington, KY, USA, 2017. [Google Scholar]
  3. Moorehead, S.J.; Wellington, C.K.; Gilmore, B.J.; Vallespi, C. Automating orchards: A system of autonomous tractors for orchard maintenance. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) Workshop on Agricultural Robots, Algarve, Portugal, 7–12 October 2012. [Google Scholar]
  4. Spekken, M.; de Bruin, S. Optimized routing on agricultural fields by minimizing maneuvering and servicing time. Precis. Agric. 2013, 14, 224–244. [Google Scholar] [CrossRef]
  5. Bochtis, D.D.; Vougioukas, S.G. Minimising the non-working distance travelled by machines operating in a headland field pattern. Biosyst. Eng. 2008, 101, 1–12. [Google Scholar] [CrossRef]
  6. Nazarahari, M.; Khanmirza, E.; Doostie, S. Multi-objective multi-robot path planning in continuous environment using an enhanced genetic algorithm. Expert Syst. Appl. 2019, 115, 106–120. [Google Scholar] [CrossRef]
  7. Sch, A.; Fle, T.R.; Mohamed, S.; Uchiyama, N.; Sawodny, O. Coverage path planning for mobile robots using genetic algorithm with energy optimization. In Proceedings of the 2016 International Electronics Symposium (IES), Denpasar, Indonesia, 29–30 September 2016. [Google Scholar]
  8. Ryerson, A.F.; Zhang, Q. Vehicle path planning for complete field coverage using genetic algorithms. Agric. Eng. Int. CIGR J. 2007, IX, 1–11. [Google Scholar]
  9. Conesa-Muñoz, J.; Pajares, G.; Ribeiro, A. Mix-opt: A new route operator for optimal coverage path planning for a fleet in an agricultural environment. Expert Syst. Appl. 2016, 54, 364–378. [Google Scholar] [CrossRef]
  10. Conesa-Mu, N.; Oz, J.; Bengochea-Guevara, J.M.; Andujar, D.; Ribeiro, A. Efficient distribution of a fleet of heterogeneous vehicles in agriculture: A practical approach to multi-path planning. In Proceedings of the 2015 IEEE International Conference on Autonomous Robot Systems and Competitions, Vila Real, Portugal, 8–10 April 2015. [Google Scholar]
  11. Bakhtiari, A.A.; Navid, H.; Mehri, J.; Bochtis, D.D. Optimal route planning of agricultural field operations using ant colony optimization. Agric. Eng. Int. CIGR J. 2012, 13, 1–10. [Google Scholar]
  12. Moghadam, K.E.; Vahdanjoo, M.; Jensen, A.L.; Bochtis, D.D. An Arable Field for Benchmarking of Metaheuristic Algorithms for Capacitated Coverage Path Planning Problems. Agronomy 2020, 10, 1454. [Google Scholar] [CrossRef]
  13. Khajepour, A.; Sheikhmohammady, M.; Nikbakhsh, E. Field path planning using capacitated arc routing problem. Comput. Electron. Agric. 2020, 173, 105401. [Google Scholar] [CrossRef]
  14. Utamima, A.; Reiners, T.; Ansaripoor, A.H. Optimisation of agricultural routing planning in field logistics with Evolutionary Hybrid Neighbourhood Search. Biosyst. Eng. 2019, 184, 166–180. [Google Scholar] [CrossRef]
  15. Zhu, X.; Ding, Y.; Cai, X.; Wang, H.; Zhang, X. Optimal schedule for agricultural machinery using an improved Immune-Tabu Search Algorithm. In Proceedings of the 2017 36th Chinese Control Conference (CCC), Dalian, China, 26–28 July 2017. [Google Scholar]
  16. Wang, B.; Li, S.; Guo, J.; Chen, Q. Car-like mobile robot path planning in rough terrain using multi-objective particle swarm optimization algorithm. Neurocomputing 2018, 282, 42–51. [Google Scholar] [CrossRef]
  17. Seyyedhasani, H.; Dvorak, J.S. Using the Vehicle Routing Problem to Reduce Field Completion Times with Multiple Machines. Comput. Electron. Agric. 2016, 134, 142–150. [Google Scholar] [CrossRef] [Green Version]
  18. Seyyedhasani, H.; Dvorak, J.S.; Roemmele, E. Routing algorithm selection for field coverage planning based on field shape and fleet size. Comput. Electron. Agric. 2019, 156, 523–529. [Google Scholar] [CrossRef]
  19. Utamima, A.; Reiners, T.; Ansaripoor, A.H. Evolutionary estimation of distribution algorithm for agricultural routing planning in field logistics. Procedia Comput. Sci. 2019, 161, 560–567. [Google Scholar] [CrossRef]
  20. Fabre, S.; Soures, P.; Taix, M.; Cordesses, L. Farmwork path planning for field coverage with minimum overlapping. In Proceedings of the ETFA 2001 8th International Conference on Emerging Technologies and Factory Automation, Antibes-Juan les Pins, France, 15–18 October 2001. [Google Scholar]
  21. Jensen, M.F.; Bochtis, D.S.O.; Rensen, C.G. Coverage planning for capacitated field operations, part II: Optimisation. Biosyst. Eng. 2015, 139, 149–164. [Google Scholar] [CrossRef]
  22. Muthukumaran, S.; Ganesan, M.; Dhanasekar, J.; Loganathan, G.B. Path Planning Optimization for Agricultural Spraying Robots Using Hybrid Dragonfly–Cuckoo Search Algorithm. Alinteri J. Agric. Sci. 2021, 36, 2564–7814. [Google Scholar] [CrossRef]
  23. Postmes, J.; Botteghi, N.; Sirmacek, B.; Kamilaris, A. A System for Efficient Path Planning and Target Assignment for Robotic Swarms in Agriculture. In Proceedings of the IoT ‘21, 11th International Conference on the Internet of Things, St.Gallen, Switzerland, 8–12 November 2021; pp. 158–164. [Google Scholar]
  24. Oksanen, T.; Visala, A. Coverage path planning algorithms for agricultural field machines. J. Field Robot. 2009, 26, 651–668. [Google Scholar] [CrossRef]
  25. Bochtis, D.; Griepentrog, H.W.; Vougioukas, S.; Busato, P.; Berruto, R.; Zhou, K. Route planning for orchard operations. Comput. Electron. Agric. 2015, 113, 51–60. [Google Scholar] [CrossRef]
  26. Jensen, M.A.F.; Bochtis, D.; Sørensen, C.G.; Blas, M.R.; Lykkegaard, K.L. In-field and inter-field path planning for agricultural transport units. Comput. Ind. Eng. 2012, 63, 1054–1061. [Google Scholar] [CrossRef]
  27. Utamima, A.; Reiners, T.; Ansaripoor, A.; Seyyedhasani, H. Contemporary Approaches and Strategies for Applied Logistics, 10th ed.; IGI Global; University of Otago: Dunedin, New Zealand, 2018; pp. 261–283. [Google Scholar]
  28. Cheikhrouhou, O.; Khoufi, I. A comprehensive survey on the Multiple Traveling Salesman Problem: Applications, approaches and taxonomy. Comput. Sci. Rev. 2021, 40, 100369. [Google Scholar] [CrossRef]
  29. Storn, R.; Price, K. Differential evolution—A simple and efficient heuristic for global optimization over continuous spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  30. Meng, Z.; Liu, H.; Wang, H.; Fu, H. Optimal path planning for agricultural machinery. Trans. Chin. Soc. Agric. Mach. 2012, 43, 147–152. [Google Scholar]
  31. Mi, M.; Huifeng, X.; Ming, Z.; Yu, G. An improved differential evolution algorithm for TSP problem. In Proceedings of the 2010 International Conference on Intelligent Computation Technology and Automation, Changsha, China, 11–12 May 2010. [Google Scholar]
  32. Wang, X.; Xu, G. Hybrid differential evolution algorithm for traveling salesman problem. Procedia Eng. 2011, 15, 2716–2720. [Google Scholar] [CrossRef] [Green Version]
  33. Mei, M.; Xue, H.F.; Zhong, M.; Gu, Y. An improved differential evolution algorithm for TSP. Inf. Technol. 2011, 15, 2716–2720. [Google Scholar]
  34. Brest, J.; Boškovič, B.; Žumer, V. An improved self-adaptive differential evolution algorithm in single objective constrained real-parameter optimization. In Proceedings of the IEEE Congress on Evolutionary Computation, Barcelona, Spain, 18–23 July 2010; IEEE: Piscataway, NJ, USA, 2010. [Google Scholar]
  35. Zhang, L. Multiple Mutation Strategy of Adaptive Differential Evolution Algorithm Research. Int. Core J. Eng. 2021, 7, 1–9. [Google Scholar]
  36. Qin, A.K.; Huang, V.L.; Suganthan, P.N. Differential evolution Algorithm with strategy adaptation for global numerical optimization. IEEE Trans. Evol. Comput. 2009, 13, 398–417. [Google Scholar] [CrossRef]
  37. Seyyedhasani, H.; Dvorak, J.S. Reducing field work time using fleet routing optimization. Biosyst. Eng. 2018, 169, 1–10. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of multiple tractors working in standard peach orchard.
Figure 1. Schematic diagram of multiple tractors working in standard peach orchard.
Agronomy 12 00856 g001
Figure 2. Schematic diagram of the turning model for a tractor (a) Π-Turn; (b) Ω-Turn.
Figure 2. Schematic diagram of the turning model for a tractor (a) Π-Turn; (b) Ω-Turn.
Agronomy 12 00856 g002
Figure 3. Schematic diagram of partition operation.
Figure 3. Schematic diagram of partition operation.
Agronomy 12 00856 g003
Figure 4. Flowchart of improved differential evolution algorithm.
Figure 4. Flowchart of improved differential evolution algorithm.
Agronomy 12 00856 g004
Figure 5. Schematic diagram of the distribution of peach trees (a) Rectangular distribution; (b) Trapezoidal distribution; (c) Irregular distribution.
Figure 5. Schematic diagram of the distribution of peach trees (a) Rectangular distribution; (b) Trapezoidal distribution; (c) Irregular distribution.
Agronomy 12 00856 g005
Figure 6. Convergence curve of differential evolution algorithm (a) Comparison of convergence under rectangular distribution of peach trees; (b) Comparison of convergence under trapezoidal distribution of peach trees; (c) Comparison of convergence under irregular distribution of peach trees.
Figure 6. Convergence curve of differential evolution algorithm (a) Comparison of convergence under rectangular distribution of peach trees; (b) Comparison of convergence under trapezoidal distribution of peach trees; (c) Comparison of convergence under irregular distribution of peach trees.
Agronomy 12 00856 g006
Figure 7. Comparison of effective work capacity and effective work capacity per vehicle.
Figure 7. Comparison of effective work capacity and effective work capacity per vehicle.
Agronomy 12 00856 g007
Table 1. Performance comparision of algorithm.
Table 1. Performance comparision of algorithm.
Number of
Agricultural
Machinery
Total Turning TimeTotal Operating Time
TTTOn(s)TTTIn(s)TTTDRn(%)TOTOn(s)TOTIn(s)TOTDRn(%)
1558.29545.952.212838.292825.950.44
2507.00490.803.191417.891388.252.09
3541.10495.978.34954.71945.111.01
4509.67490.973.67703.96696.701.03
5458.38468.632.19556.68547.680.41
6467.13465.030.45487.44486.240.24
Mean506.93492.893.341159.831148.320.87
Table 2. Comparison of total turning time and total operating time.
Table 2. Comparison of total turning time and total operating time.
Peach
Orchard
Number of TractorsTotal Turning TimeTotal Operating Time
TTTOn/sTTTIn/sTTTDRn/%TOTOn/sTOTIn/sTOTDRn/%
Rectangle1839.3536.836.03119.32816.89.7
2816.9523.635.91548.41414.58.7
3794.5492.538.01077.2956.911.2
4772.0489.736.6763.0698.48.5
5749.6495.533.9605.9563.96.9
6732.5468.436.1528.3488.37.6
Mean784.1501.136.11273.71156.58.8
Trapezoid1839.3560.533.22799.32520.510.0
2816.9524.435.81484.41257.115.3
3794.5527.733.61053.2838.520.4
4772.0530.931.2755.0633.516.1
5749.6490.834.5605.9507.816.2
6732.1506.030.9528.3434.217.8
Mean784.1523.433.21204.41031.916.0
Irregular1837.8572.031.72477.82212.010.7
2815.4522.236.01227.71087.411.4
3793.1520.834.3852.8751.511.9
4770.5510.933.7602.6546.79.3
5748.8501.733.0477.9433.49.3
6731.2463.936.6416.1368.711.4
Mean782.8515.234.21009.2899.910.7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Han, X.; Lai, Y.; Wu, H. A Path Optimization Algorithm for Multiple Unmanned Tractors in Peach Orchard Management. Agronomy 2022, 12, 856. https://doi.org/10.3390/agronomy12040856

AMA Style

Han X, Lai Y, Wu H. A Path Optimization Algorithm for Multiple Unmanned Tractors in Peach Orchard Management. Agronomy. 2022; 12(4):856. https://doi.org/10.3390/agronomy12040856

Chicago/Turabian Style

Han, Xiao, Yanliang Lai, and Huarui Wu. 2022. "A Path Optimization Algorithm for Multiple Unmanned Tractors in Peach Orchard Management" Agronomy 12, no. 4: 856. https://doi.org/10.3390/agronomy12040856

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop