Efﬁciency Analysis and Improvement of an Intelligent Transportation System for the Application in Greenhouse

: In view of the future lack of human resources due to the aging of the population, the automatic, Intelligent Mechatronic Systems (IMSs) and Intelligent Transportation Systems (ITSs) have broad application prospects. However, complex application scenarios and limited open design resources make designing highly efﬁcient ITS systems still a challenging task. In this paper, the optimal load factor solving solution is established. By converting the three user requirements including working distance, time and load into load-related factors, the optimal result can be obtained among system complexity, efﬁciency and system energy consumption. A specialized visual navigation and motion control system has been proposed to simplify the path planning, navigation and motion control processes and to be accurately calculated in advance, thereby further improving the efﬁciency of the ITS system. The validity of the efﬁciency calculation formula and navigation control method proposed in this paper is veriﬁed. Under optimal conditions, the actual working mileage is expected to be 99.7%, and the energy consumption is 83.5% of the expected value, which provides sufﬁcient redundancy for the system. In addition, the individual ITS reaches the rated operating efﬁciency of 95.86%; in other words, one ITS has twice the ability of a single worker. This proves the accuracy and efﬁciency of the designed ITS system.


Introduction
With the aging of the global population and the deepening of urbanization, the contradiction between the loss of agricultural population, rising labor costs and the demand for agricultural production and supply has become an important issue faced by the sustainable development of agricultural production [1]. The development of high-efficiency, high-quality, low-cost intelligent mechatronic systems that can replace human operations is an important means of dealing with the current situation from an engineering perspective [2]. Agricultural labor intensity, high production costs, only harvesting costs can be accounted for more than 25% of total production costs [3,4]. The labor intensity in operations such as pest control is also very high, and directly threatens the health of workers [5]. Automated Guided Vehicle (AGV) as an important part of ITS can effectively reduce the risk of workers in hazardous operations [6], help increase production efficiency and deal with the lack of labor due to the aging of the population [7]. A variety of agricultural robots, including grazing [8], farmingfarming [9], fruit [10] collection, and sorting [11] robots, have been applied.

•
We have provided an optimal load factor solving solution. By converting the three user requirements including working distance, time and load into load-related factors, the optimal result can be obtained among system complexity, efficiency and system energy consumption. The specifications of the main components of the ITS such as drive structure, power and battery components are constrained to a limited range according to the load factor and the work scenario. The specification of constraints helps other researchers determine the design of ITS and the choice of key components.

•
We propose a special weighted connected graph structure to support map modeling and navigation applications. Based on this, a specialized visual navigation and motion control system has been proposed. With a fixed control action and boundary-defined Proportion Integral Differential (PID) control algorithm, the repeat path accuracy of the error < 0.02 m can be obtained without relying on third-party positioning. This makes prior path planning, accurate mileage and energy consumption prediction possible. It can help researchers and designers further improve the efficiency of the ITS system. • We present an efficient example of an ITS system and its detailed design based on the proposed method with optimal calculation results. Test results include working paths and mileage, speed and energy consumption are described in detail under different test conditions, including conditions beyond the rated load. The performance of the system, which is verified in the test environment, is described in detail.
The working environment, requirements, and design goals are described and three known input parameters in detail in Section 3. In Section 4, the overall design framework of the ITS is explored according to the working environment, and then the requirements are converted into three key constraints of load-volume ratio, working distance and work efficiency. In Section 5, a specialized visual navigation and motion control system has been proposed to simplify the path planning, navigation and motion control processes. The ITS prototype designed according to this constraint is elaborated, which is one of the main contributions of this paper. In Section 6, an improved ITS system is implemented as a case, and its detailed design parameters are calculated based on instance greenhouse environment. In addition, then, the validity of the proposed constraints and the efficiency of the ITS prototype are verified by actual verification. Section 8 concludes this paper.

Related Works
Agriculture is a labor-intensive industry. With the degeneration of the demographic dividend and the rapid growth of the population, the pressure on agricultural production is increasing. As the most important and labor-intensive operation in agricultural production [5], the pesticides used in pest control have high toxicity and corrosive characteristics [19], which greatly threaten the health of production personnel. Strict and comprehensive protection measures, as well as a lower level of participation of practitioners, further increase the cost of such operations [20]. The United States Environmental Protection Agency (US-EPA) has finalized stronger standards "Agricultural Worker Protection Standard (WPS)" for people who apply for Restricted Use Pesticides (RUPs) [20,21]. These revisions to the Certification of Pesticide Applicators rule will reduce the likelihood of harm from the misapplication of toxic pesticides.
Agricultural automation and intelligent development are one of the effective means of solving these problems [18].The agricultural robot originated from agricultural machinery is a new type of multifunctional agricultural machinery. It is the result of adding intelligent attributes to agricultural machinery [22]. The earliest agricultural robots date back to the autopilot tractors of the 1970s [18]. Up to now, various agricultural robots have emerged, such as shearing robots, milking robots [8], transplanting robots [23][24][25], grafting robots [23,24], harvesting robots [10,26], and weeding robots [27]. According to the use and operation characteristics, the existing agricultural robots can be roughly divided into livestock management robots, field farming management robots, fruit and vegetable harvesting robots, seedling breeding robots, agricultural product sorting robots, etc.
Due to the complexity of agricultural operations, the complexity of work tasks, and the current level of software and hardware technology development, there are many types of agricultural robots and prototypes, but there are only a handful of cases that can effectively improve the efficiency of production. The more successful cases of agricultural robots are either similar to industrial lines or similar to large agricultural machines such as milking robots, agricultural sorting robots, self-driving tractors, weeding robots, etc. It must be said that, compared with the industrial robot represented by the dexterous mechanical arm [28], the highly intelligent transformers in the traditional concept of human beings [29], there is still a big gap in the development level of agricultural robots. Among the factors such as the agricultural work scene, the complexity of the work tasks and the current development level of hardware and software technology, the complexity of the work scene is the key factor [2,28]. The traditional agricultural production conditions with manual work as the main body do not match the requirements of the robot for the structured working environment, which is an important factor that restricts the agricultural robots from entering the laboratory and entering the farmland. Relative to the industrial structural environment, the special conditions of the sun in the agricultural production environment [12], the disordered crops and the soft and fragile working objects [30] are the objective problems faced by the robots in the agricultural environment [31].
Compared to the open farmland environment, the Dutch greenhouse has the potential to produce under a variety of weather conditions. These closed plant production systems, such as artificially illuminated and highly insulated plant factories, have offered perspectives for urban food production [13]. In terms of economic efficiency, the annual Net Financial Return (NFR) for the Dutch greenhouse is very high [32]. Even so, there are fewer agricultural robot systems available in the greenhouse. Compared with the complex agricultural environment, the industrial application environment is "compendious": the environmental closure and the highly structured operation. This is an important realization condition for the rapid development of industrial robots. Warehousing robots receive special attention among them. There are also a series of intelligent robot systems for greenhouse development, some of which are practical applications, such as Amazon KAVA, Richo M2, River System, Meituan Segway, Postmates, etc. [16,[33][34][35][36][37][38][39][40][41][42][43][44][45]. Comparison of major ITS systems in recent years is as shown in Table 1. As shown in Table 1, there are many robot systems that are used in different fields. However, the useful reference for design is very limited. In addition to different application environments and conditions, larger factors may be a limitation of commercial intellectual property. In the above robot system, only the Segway system [45] supports limited open source, and other products or prototype systems do not provide the necessary technical information. Although these systems provide some design references, the role of efficient ITS system design in the Dutch greenhouse application scenario is limited. Therefore, the author's team hopes to analyze the various factors and key conditions in the design one by one based on the project research experience and the basic conditions that need to be given. A normalized conditional formula is given giving detailed design and performance testing, providing a visual reference for other researchers.

An Application Case: ITS in a Greenhouse Environment
The application environment comes from an ITS design project of our team. Six Dutch greenhouses with an area of 6000 m 2 are deployed in the research base of the cooperative company. Figure 1 shows the schematic diagram of one of the standard greenhouses. As shown in Figure 1b, we can divide the entire space into three parts: 1 operating area, 2 working area, and 3 crop growing area. The operating area is mainly used for installation, temporary storage of equipment, and equipment turnover in the working area. Workers are also working in the area. The working area and the crop area are deployed at intervals, and there are n + 1 crop areas separated by n work areas. The length of each work area is l n (referred to as height H), the width is d p and the spacing between work areas is d. The width of the entire working interval is W. Taking pest control as an example, the role of pesticide spray control is limited because a large number of insect pests have light-shielding properties and gather behind the foliage. In addition, the role of other control measures is also limited, such as larval boards and other means. This requires workers to use equipment to spray pesticides from the bottom up to remove pests. Corrosive, highly toxic, volatile pesticides are bound to pose a threat to the health of workers, even when protective measures exist. In addition, carrying out the work in such a large range makes the workers have a very high working intensity. This is determined by the characteristics of the pesticide such as volatility and precipitation. The spray work must be completed within the specified time after the pesticide configuration is completed. Otherwise, its medicinal properties will be greatly weakened or even have side effects. Therefore, the use of highly efficient robots to perform similar tasks is not only economical, but also effective in protecting the health of workers.

Definition of Requirements and Input Conditions
It is challenging to design an ITS that meets user requirements and is optimal especially when there is a limited reference. This requires us to start the analysis from known and user-expected conditions. Firstly, we need to convert the user's requirements into key metrics based on the requirements analysis, as shown in Figure 2. The key metrics are then converted into sub-conditions that affect the ITS design, such as operating speed, load, and output torque, and the boundaries of these conditions are determined. Finally, the optimal values of key conditions are obtained under multiple conditions and their boundary constraints, and the detailed design scheme of ITS is determined. (1) Work efficiency. If the time required for a standard job is T norm and the actual completion time is T practical , the productivity e p can be defined as: This means that the more time is spent, the lower the efficiency. T norm is derived from actual needs and historical work experience, such as T norm = 5400 s in this scenario. This is determined by the nature of the liquid and the operating specifications. In addition, in the above conditions, the efficiency of a single worker is e p = 50% because one greenhouse requires two workers to work together. Thus, we can evaluate the efficiency of the designed ITS through e p .
(2) Work mileage. As shown in Figure 1b, the spraying operation only occurs in the working area, and we refer to such areas as "effective working areas". Obviously, for a given site, the "effective working distance" for each work area is known as l 1 , l 2 , · · · , l n . Therefore, the total effective working distance L can be described as the second condition: If the length of each workspace is the same and the length is H, then L = n × H. The moving distance in the actual job will far exceed this value, whether it is manual or ITS because it needs to work-interval within the interval, switching the working range, and returning the replenishment. In addition, the appropriate ITS system design and path planning system will help reduce these "extra" 'travels'-these will directly affect the actual time spent on the operation T practical , which will affect the work efficiency e p .
(3) Unit workload. Since the total volume of liquid medicine vol tank and the effective working distance L are known, the unit distance dose SV p can be calculated: Of course, this value is more appropriately calculated by professional agronomists based on different pests and crops. SV p will directly interfere with two key issues: (1) determine the specifications of the pumb deployed on the ITS; (2) affect the specification and continuous working time of vol ITS -this will further affect the mission planning of the ITS and the implementation process, and ultimately affect the efficiency of the work.

Multiple Conditional Constraints Reasoning
In this section, we will explore in detail and contraction the main conditions affecting efficient MIS design based on known input conditions and expectations. The first problem to be identified is to choose the appropriate overall design for the MIS. An integrated design that differs from existing split operations is more advantageous. As shown in Figure 1b, the hose line connects the separate spray assembly to the pesticide tank assembly in the existing split mode of operation. The main advantage of splite scheme is that the existing syrup supply system can be used directly. In addition, the complexity of the MIS system will be reduced because it does not need to carry a heavy tank and pump system itself. An integrated scheme that differs from existing split operations is more advantageous. However, the solution makes the pipeline system extremely complicated in practical applications. There are two reasons for this: firstly, the ITS needs to enter the working area so that the length of the pipeline will exceed l H . Secondly, when the ITS switches the working area, the tank also needs to move synchronously to the corresponding operating area, which will increase the complexity of the synchronous control system and the communication system. In comparison, the presticide tank, the sprinkler system, and the ITS system body are brought together in the integrated design. Its main advantage is its compact structure, no external equipment support, and good independence. The key questions here can be expressed using vol ITS , as shown in Figure 3. It can be seen that, when the fixed V moto is indirectly fixed to the running speed of the AGV, the larger volume is beneficial to reduce the working mileage and working time, which will directly improve the final working efficiency e p , reducing the distance of reentry due to lesser replenishment count. In addition, the side effect is that it directly increases the system's energy consumption and requires greater output torque-500 L is a big challenge for ITS design. In other words, a viable, efficient ITS system relies on the determination of the optimal vol ITS .

Load Constraint and Factor
In order to determine the appropriate vol ITS , we will search for favorable clues from known conditions: containers are a huge challenge, but not an incomplete design. However, in the feasible design, we will still find the favorable constraints: First, the width of the working area l d and the interval spacing will limit the width and length of the MIS and tank, respectively; then, the height and centroid of the tank will follow the volume-increasing and increasing, but the stability of the system will decrease. Finally, greater capacity means greater torque demand and energy consumption pressure.
(1) Load factor. SV p reflects an indirect relationship between tank and distance, that is, the working distance that the agent in tank can support should be an integer multiple of l H . It can effectively reduce the number of invalid round trips due to insufficient liquid pesticide. That additional replenishment will generate additional working mileage, which will directly affect the efficiency of ITS. It can be calculated that the amount of medication in a single work area is vol tank /n = SV p × l H . In addition, the volume of vol veh should be an integral multiple of it, as a key constraint of an efficient design. We defined m as load factor, and: where m ∈ [1, n] is an integer. In addition, the number of supplemental options Pst can be calculated when the pump position is fixed: In this way, we simplify the complex capacity problem to the relationship between the two scalars m and n.

Work Planning and Mileage Calculating
Through Equation (5), we can determine the number of times ITS returns to the pesticide tank, and we can also calculate when it is necessary to return the supply. Therefore, we can divide work planning into two parts: spray and replenishment, as shown in Figure 4. Correspondingly, the calculation of the mileage is also divided into two parts. (1) Mileage of spray operation.
Each sub spray sequence such as A 1 → B 1 → C 1 consists of three parts: 1 entering work area A 1 ; 2 performing spray B 1 ; 3 returning C 1 ready to enter the next work area. Equation (4) reflects the critical relationship between the operating distances L and vol ITS . In addition, this allows us to ensure that no pesticide shortages occur during the spray step B 1 . The spray operation distance L spary.i of each work area can be expressed as: The total spray operation distance L spary.total can be calculated because the working conditions l i = l H are fixed. In addition, we also incorporate the return step E into the process: (2) Mileage of refueling operation.
We avoided interrupting the spray operation in the middle of the working area according to Equation (5). At this point, the pesticide exhaustion at node n i , and the "return-replenishment-reset" process needs to be performed. By defining the serial number of the replenishment as j, we can calculate the line number of the dosing that occurred in i j via Equation (5): That is, the refueling must occur after one spray is completed. At this point, the ITS will return to the tank to perform the dosing, and its moving distance L spary.i j : The total distance of this replenishment operation L supply.total can be calculated: As a key condition, total working mileage L , regarded as an enhancement of the input condition L, can be expressed as:

Resistance with Load Capacity
Changes in vol ITS will directly affect the system's resistance F res , and also require output torque F out > F res to drive ITS. Consider the weight WT ITS of the ITS itself and its load WT load , the driving resistance F res can be calculated: where the e f is the empirical coefficient [49], mainly considering the deceleration of motor, transmission part of the loss. In addition, to consider when the work surface water stains will have a reduced friction coefficient, one needs to leave a margin for torque. For example, the efficiency of the primary reducer is about 85∼96%. Considering the power loss and redundancy caused by the variable friction coefficient, the effective power is about 70% of the theoretical, that is, e f ≈ 1.42. Select the maximum friction coefficient of the asphalt pavement µ s [50] in order to simplify the calculation. WT load ≈ vol ITS when the specific pesticide density is not considered. It can calculate the resistance data under the load and tyre size, according to Equation (12). This is to calculate the motor output torque and even motor drive system selection provided by the basis. Whether a reducer is needed to amplify the output torque depends on the actual size of F res . This will also limit the selection of the motor.

Speed Condition under Efficiency Constraints
Let us define the average speed of ITS to beV; then, the actual working time T practical can be calculated by Equation (11): Thus, if you want to ensure e ITS ≥ e p , ask for T practical ≥ T total , and we can get: It is more accurate to substitute the total fixed replenishment time T supply = vol tank /Q and the ITS state to determine the loss time T loss into Equation (14). A more accurate average speed can be expressed asV :V Then, according to the Equations (12) and (14), the output speed, torque and power of the motor can be pushed back to provide the basis for the selection and further detailed design. The output speed V PGR of the planetary gear reducer (PGR) and the output speed V moto of the motor can be calculated by Equation (16) with different loads according to the rotational speed formula. In addition, we further shrink the selection conditions of the drive system:

Structural Stability Conditions
A larger tank volume will increase the center of gravity of ITS and affect the stability of the system structure. In actual operation, it is also necessary to consider the ability of ITS to pass through obstacles. In addition, it also needs to consider the system stability under variable load conditions [51]. Calculated in a circular or equivalent circular tank with a radius of r t , tank height is h t and vehicle height is h v . This allows us to calculate the aspect ratio of vehicle diameter and total height, as shown in Figure 5. Excessive height and proportion will affect the stability of the system. Consider the capacity to run the system tilt angle θ sloped . Therefore, the height of the vehicle barycenter can be estimated. At this point, the vertical center of mass must be located within the vehicle support point. If the radius of the AGV is r v , and the safety height can be calculated taking into account the mass factor h mass :

Normalization and Optimal Solution of Multiple Conditional Constraints
Although the above conditions can converge the design conditions of the ITS to a large extent, the single is still insufficient to give a specific most preferred result. These conditions are all related to the capacity coefficient m and are divided into two categories: (1) proportional, as vol veh , L, F res ; and, (2) inversely, as Pst,V. Operating speedV and output torque F res are the two factors that have the greatest impact on ITS. The speed of workV can be replaced by the working mileage L when the upper limit of the total working time is fixed. In addition, they are all related to the volume of the pesticide tank-further related to the factor m. By mapping them to the same scalar space by normalization, although their units are different, the optimal solution can be calculated. Therefore, taking the factor m as a reference variable, we then employ the Matlab function "mapminmax(parameters, 0, 1)" to classify the main condition parameters into scalars and compare them.
Through the analysis, we can pass the factors in Sections 3 and 4 through a normalized heuristic method. ∑ Pst j=1 L spary.i j in Equation (11) can be simplified using the summation formula of the arithmetic progression S n = n(a 1 +a n ) 2 , where a 1 , a 2 , · · · , a n is an arithmetic progression. In addition, we can create an optimal evaluation function f (m) related to m because parameters other than m can be obtained through specific examples: When f (m) → 0, the optimal value of the result and m can then be calculated.

Efficient Visual Navigation and Motion Control System
Before giving the test site and other conditions, we need to briefly introduce the designed ITS navigation system because this will have a direct impact on actual operation and testing. Considering the operating environment in the greenhouse, the design principles of the navigation system are: simple, reliable, accurate and inexpensive. We have proposed and established a special navigation-control scheme which includes a map system, visual navigation system and navigation control system.

Specialized Undirected Weighted Connected Graph and Fixed Action Control
We propose a special map representation model called: specialized undirected weighted connected graph. As shown in Figure 6a, simplify the expression of the map by adding some restrictions, which is different from the general navigation map expression.

•
The relative positions between the nodes are limited to four types of up, down, left, and right. For example, for a node N i , the node adjacent to it is only allowed to appear in its 0 • , 90 • , 180 • and 270 • in four positions. This also limits the direction of ITS movement to four, which can significantly reduce the complexity of graph and control algorithms..

•
The association between the nodes is given a weight of w j according to the distance between the nodes. For example, the weight between nodes n 0 and n 1 is w 1 = d, and the weight between nodes n 1 and n 1 is w 1 = H. This weight value helps to control the distance heuristics in the algorithm.
This design can significantly reduce the cost of map modeling and path planning algorithms. Since the starting point and the target node are given, only one optimal/shortest path is generated. For example, the shortest path between n 0 and n 1 can be expressed as n 0 → n 1 → n 1 . Nodes can be expressed only by one edge, and there is no need to equidistantly set a large number of auxiliary nodes (a small number of auxiliary nodes that do not participate in the calculation help the system's own state check, avoiding the possibility of falling into the algorithm when moving long distances. In the problem of defects).
Based on this specialized graph navigation map, a "fixed action control" scheme was proposed as shown in Figure 6b. There are only four sports programs here: 1 Turn Left 90 • ; 2 Turn Right 90 • ; 3 Turn Around 180 • and 4 Forward. These four types of actions are further divided into two types of events to match the event processing mechanisms: forward and turn. Thus, whenever the vision system recognizes a node identifier, the algorithm will prioritize whether to ignore the point to continue straight or perform a "special turn". Take the n 2 → n 3 → n 1 → n 3 process as an example: First, AGV will perform a forward action until the center of the ITS p overlaps with N 3 ; then, the "left90" action will be triggered, i.e., the AGV will turn left and its azimuth is 90 • to point to n 1 . Then, go straight to the n 1 node and trigger the "trun180" action, i.e., ITS rotates 180 • in place and points to n 3 . Finally, go straight to the n 3 node to complete the task.
This control method is obviously not optimal, such as moving the farthest distance and consuming more time when cornering. However, this simple control method does not require complex algorithm support.

Vision-Based Navigation System
Each node Node or n i uses a "Quick Response Code (QR Code)" to mark it [17]. Thus, the ITS knows its current relative position, thereby constructing a complete path navigation system and supporting the concrete navigation algorithm. In addition, in a "QR Code Annotated Map" (QR-CAM) visual navigation system based on a specialized graph [52,53], as shown in Figure 7a, guide line connections are used between nodes to provide effective assistance in visual guidance and control. (1) We can determine the actual distance between the center point of the ITS, P, and the QR Code, using the "Pixel-Distance" formula established by camera calibration. ITS will trigger the "turn" event and perform the turn immediately after the l node distance. In addition, when the event is triggered, the cumulative state of the system is reset-the accumulated error is also set to zero, which increases the reliability of the navigation and control algorithms.
(2) Set 5 evaluation points S = {s 1 , s 2 , s o , s 3 , s 4 } on both sides of the image center. The distance between the evaluation points is ∆d. Thus, when the center of the ITS is near the guide line, its deflection angle can be calculated as θ 1 = arctan(∆d/l O ) and θ 2 = arctan(2 × ∆d/l O ). This allows a quick evaluation of the positional relationship between the ITS and the guide line, thus providing a p-value reference for the PID control (if the settings are appropriate, in an ideal environment, only proportional adjustment can be used to complete the effective control of the ITS).
(3) The specific control process is divided into three parts, as shown in Figure 7b. In addition, this process can be called "Setting-Find-Task Execution". At first, the staff sets the task target and work path through the terminal equipment such as the wired or wireless touch panel on the vehicle [54]. In addition, then, the vehicle will then complete the movement between the nodes in units of sub-paths and eventually reach the target node. For example, the n 0 to n 1 path can be split into n o → n 1 "forward", "turn right 90 • ", n 1 → n 1 "forward" as shown in Figure 6c. That is, all tasks can be planned and calculated in advance, just like Programmable Logic Controller (PLC), and each task can be completed in steps-this can greatly reduce the computational burden caused by dynamic programming and real-time control. The motion control only has two simple control processes: 2 "turn" and 3 "go forward". This ensures that the vehicle can move precisely along the guide line, inspired by θ 1 and theta 2 , based on the PID control algorithm. The deviation can be controlled at 0.5% or ±2 m, even if there are machining and assembly errors.
(4) Repeat the "move-search" process according to the planned route until all tasks have been completed or reached a "termination point". (Note: The user can manually set the termination node or when the ITS needs to return to the service station when it detects a failure [55].)

Qualification Calculation Based on Instance
The example parameters of the greenhouse environment can be given according to the greenhouse structure shown in Figure 1 and the user requirements and input conditions in Section 3. Table 2 lists instances of the main parameters. The above example parameters are brought into the formula one by one. Bring the parameters listed in Table 2 into the formula one by one, derived according to the theory. In addition, finally, determine the optimized ITS selection design.
1. First of all, the quality factor m ∈ [1,22] can be determined according to the working area count n = 22. 2. Then, the total working distance can be calculated according to the factor m that L ∈ [2376, 3300]. 3. In addition, the total effective working time limit can be calculated: When the three factors of workload, distance, and time are known, the resistance of the ITS system, F res , and the output torque of the motor can be calculated. As shown in Figure 8 by a two-coordinate graph ('plotyy'), here is a comparison of the resistance and the required output torque of three different radius drive wheels with different loads. It can be seen that as the load factor m increases: (1) The average output speed of the motor are 55, 72 and 88 rpm with different D wheel or D. Respectively, minor changes when m > 3. (2) However, the output torque is increasing rapidly, Fres ∈ (13.7∼21.9) N · m when m = 10. When m = 22, it grows to 28.4∼45.5 N · m, which is over the limit of ordinary DC motors and reducers.
Because the general motor is characterized by low torque and low torque, it is difficult to meet the drive of large loads. Therefore, it is necessary to use PGR to increase the output torque. The range of conditions for output torque can be determined by Figure 8 and the motor selection code. "Limitation conditions (3)" F out = 10 N · m when m > 10 to offset the resistance according to Equation (12). In addition, the unit N of F res or F out has been converted to N · m according to the power wheel radius.
The output torque is related to the driving wheel diameter, the motor torque and the deceleration ratio of the reducer. The ITS can only move when F res ≤ F out /D wheel (especially static state). In addition, the V PGR ∈ [50∼140] rpm and V moto ∈ [300∼1400] rpm based on drive wheel diameter 5-8 inchs.
According to the analysis of Section 4.6, the key conditions of F res and L are combined with factor m, and the results are as shown in Figure 9. Substituting Equation (18) into Equation (4) according to the given condition, we can see that m = 5.
Larger volume capacity provides more redundancy and shorter distance, such as nozzle loss during liquid spray. It should be noted that the excessive load can lead to a significant increase in the probability of ITS system failure and difficulty in diagnosing [56]. Compared to m = 5, when m = 6, the center shifts by only 0.08 m and the height of tank is about 0.65 m, which is still within the design tolerance. In addition, the working distance is reduced by 2.2%.
Based on the obtained optimal coefficient m , the structural stability of the system can be verified. The working area d p = 0.6 m and the bottom radius of the container is 0.4 ≤ r t ≤ 0.5 m. In addition, r v ∈ [0.5, 0.6] m. Bring these parameters to Equation (17) to calculate h mass = 1.45-2.25 m. In addition, considering that the liquid in the container may be shaken , the safety height should be set at least 80% of the limit value at design time. Therefore, the safety height h mass = h mass × 80% is 1.45-1.8 m.
Because at this time the bottom area of tank is about 0.196 m 2 , it can thus be limited to the car cabinet volume within 200 L, that is, m < 9. Although this result is larger than the optimal value (m = 5) we obtained in the normalization condition calculation, since the height of the vehicle is easy to estimate, the calculation condition can be contracted earlier by this method.

An ITS Instance for Greenhouse Spraying Application
We designed and implemented a complete ITS system based on the above analysis and given the greenhouse application scenario [55]. Figure 10a shows the schematic diagram of a full-featured prototype. Figure 10b shows a 3D view of the ITS. More information on subsequent versions can be found in 2017-the 15th China International Agricultural Trade Fair (http://www.xinhuanet. com/politics/2017-09/22/c_129710586.htm) and the 24th China Yangling Agricultural High-tech Achievements Expo (http://www.chinanews.com/cj/shipin/cns-d/2017/11-05/news739705.shtml). We can see the composition of the system for pesticide spraying applications by Figure 10a. It consists of nine parts, including a full-featured version of the entire ITS base system, and the final version is the result of this prototype's improvement and optimization. It can be divided into two parts from the perspective of use: (1) Pesticide spraying subsection. 5 Two 15 lwater tanks form the medicine box and the main load; 2 a self-priming booster pump (E-CHEN EC-RV-03L) is used to draw the liquid from the chamber and provide sufficient pressure to the 9 nozzle. 1 A solenoid valve system monitors the system pressure while controlling each nozzle. In addition, a variety of filtration systems are available for different liquids to prevent clogging of the above components.
(2) ITS and control system. In the prototype system, 8 the power system and ECU parts of the ITS are integrated into the chassis. The main processer 7 (NVIDIA TK1) and the 3 first-view visual sensor (wireless image transmission system based on GoPro Hero 5) are external. Although the 7 image sensor (Logitech C920, CH: Switzerland) for visual navigation is installed in the ITS chassis, its data are still processed by 3 [52].
Among them, 3 is used to identify information of each object in the environment, such as the plant or pest as 4 and is also used for remote control of the ITS by the user from the first perspective. 7 is mainly responsible for the identification of the node QR-Code and the guide line, and provides upper decision information for path planning and navigation. Eventually, this decision information will be sent to the ECU in the ITS for specific implementation. More details of the actual operation can be seen in the "video abstract".
The ITS (version "BigPan-III") shown in Figure 10b is an upgraded version of the prototype. All the necessary components, including the power supply, drive system and controller, are integrated together, and the high degrees of independence and autonomy are the advantages of this transitional version. Nine major modules are divided into two categories: (1) Control system. Compared to the prototype, the external intelligent processing device (6) is integrated into the ITS (10) after it has been redesigned and miniaturized. The image sensor (7) is replaced with a smaller camera module (45 (L)*15W (W)*15(H) mm) so that it can be mounted to the (15) position. In addition, the front and rear dual cameras allow ITS to achieve the same level of navigation in an environment that cannot be turned. The ECU (16) is updated to accommodate high temperature and high humidity environments. It integrates four currents, four voltages, three-axis accelerometers and three-axis angle sensors. This information will be combined with the rotary encoder (18) to perform attitude fusion calculations inside the Electronic Control Unit (ECU) (core chip is Freescale MC9S12X, Austin, TX, USA) to obtain the current speed of the ITS (relative) information such as position and power consumption. This information will all be uploaded to the intelligent controller (10) for decision, and the generated issued instructions will be sent to the (11) actuator to drive the actual movement of the ITS.
(2) Power and drive system. Two motors (13) (5PC120GU-24 rated power 120 W, 1800 rpm) are installed in a symmetrical manner, with the 5-inch drive wheel to form a differential drive system. The output section is equipped with a gear unit (12) (5GU10K supports 180 RPM output speed). This provides an output torque of >8.5 N·m at an output speed of 120 rpm (0.75 m/s. A battery (14) (48 V 12 Ah) allows the system to operate for at least one hour at maximum power consumption, while a second backup battery doubles the system's endurance to meet demanding mission's requirements.).
The specifications of an experimental ITS are shown in Table 3. Therefore, the total power consumption of the system is approximately P MAX = 2 × P moto + P ECU + P pump = 410 W. The total energy demand approximately 615 W · h. In practice, not all of the time is the maximum power-for example, the ECU and control computer average power consumption of about 15 W. While the pump system works less than half its working time, it can consider its average power consumption as 60 W. According to full load operation, the total power consumption of the system is aboutP total = 315 W and the total energy consumption is 472.5 W · h. Therefore, the optional lithium battery specification is 36 V· 16 Ah or 48 V · 12Ah. Except for the battery, the electricity consumption is lower than 36V, which still meets the safety requirements (GB/t 3805-2008) [57]. Considering the effect of actual loss, battery life and temperature on the battery, the battery with larger capacity should be selected in the actual selection. We can give the overall architecture of the ITS system when the two key components of the motor and battery are determined.

Test Environment Settings
We have selected a complete cabin as shown in Figure 1b for system testing. The main parameters are shown in Table 2. It is then converted to a digitized map of the XY markers to facilitate the description of the test results, as shown in Figure 11. In order to more clearly show the movement of the AGV under different loads, i.e., different load factor m. We define replenishment − work as the task group WS. Thus, when m = 6, WS 1 can be described according to Section 5: where a pesticide filling job Fill() needs to be performed before each working group starts or end in n 0 . For a 70 L/min replenishment pump, it takes about 430 s to complete the filling of vol tank = 500 L tank. In addition to m = 22, the AGV will return to the medicine box several times for replenishment. Since the total dose is the same, the total dosing time is still 430 s, but this time will be allocated to each replenishment process. Thus, the working distance of each group is L WS i , i ∈ [1, Pst] and the total distance L can be calculated: The motion control method proposed in this paper divides the motion of ITS into two types: linear and fixed turn. Therefore, through Section 6, we know that the spraying operation and the straight travel speedV 2 ≈ 0.67 m/s. In addition, within the rated load range, the acceleration both are fixed at 0.2 m/s 2 .

DAQ System
As shown in Figure 12 and Table 4, we built a DAQ system [55] to obtain the main operating parameters of the ITS runtime in order to accurately evaluate its work efficiency. An ADI AD7606 DAQ model (8-channel DAS with 16-bit) is integrated into the ECU for the acquisition and conversion of numerical data. Four of them are used for voltage detection, and the other four are used with Allegro ACS758LCB-050B linear hall current sensors (Worcester, MA, USA) for current detection. Two AVAGO HEDS-9140 encoders (San Jose, CA, USA) monitor the operating speed of the two drive wheels. An MPU-6050 triaxial accelerometer sensor is used to obtain ITS acceleration information for Inertial Measurement Unit (IMU) assisted positioning and measurement.

System Test and Result Comparison
The test mainly includes five items: working path, total mileage, speed, power and total power consumption. Based on load factor m = 1 → 8, eight sets of effective experiments were carried out here. It should be noted that the rated load of the given ITS instance is m = 6 , i.e., 138 Kg. The system's power performance will gradually decrease when the rated load is exceeded, such as the system load at m = 9. At the limit, the acceleration is only 0.03 m/s, and the maximum speed is only 0.1 m/s. This is much lower than the efficiency users need. Eight sets of valid experiments were carried out, and the corresponding load factor was set to m = 1, 2, 3, 4, 5, 6, 7, 8.
(1) The results of the work path are as shown in Figure 13. The task group represented by the different color trajectories shows the impact of different load factor m on the task flow. Obviously, a larger m makes the task less fragmented.
(2) The results of the center speed of the ITS are as in Figure 14. It can be seen that, due to the existence of the "fixed action control" described in Section 6.1, the speed curve of ITS is very similar to that of PLC, showing a high degree of regularity.
(3) The results of the work mileage comparison are as shown in Figure 15: It can be seen that, when the load factor m changes, it will cause changes in the WS process, the working path, and the work history, and exhibit similar changes. In addition, the change in total working distance is consistent with that shown in Figure 9, where the total working distance is the shortest when m = 6.
A summary of the above main test results are shown in Table 5.
(a) m = 1.          We can draw the following conclusions: • The trend of the above results is consistent with that shown in Figure 9, indicating that the optimal Equation (18) expression is valid for the design of the ITS system. • When m = 5, the main indicators of the system are optimal, and the m = 6 sub-optimal result is only 3% different. This provides some redundancy for system design.

•
The average error between the actual working mileage and the ideal calculated value is only 2%, which proves the accuracy of the navigation and control system designed in this paper. When m = 6, the actual working distance is 99.7% of the ideal value. Due to the presence of sensor errors and the response speed of the PID controller, the AGV is unlikely to reach exactly above the target node. It is possible for the AGV to stop during the previous sample-execution cycle or to stop during the next cycle. Thus, the AGV has a Deltad < 0.1 m error from the target node, and the accumulation of the error will cause a positive and negative deviation between the total mileage of the AGV and the ideal value.

•
In the above test of m, although the predicted value is exceeded by about 5%, the total system energy consumption is lower than the battery power supply capacity. Existing power supply systems are feasible, especially if a backup battery is enabled.

•
In terms of efficiency, e ITS = 95.86% when m = 6 is equivalent to two workers due to the maximum output speed limit, in other cases in which system efficiency is lower than e p , but still higher than the efficiency of a single worker e p . This proves the high efficiency of the ITS system designed in this paper.
The test results demonstrate the effectiveness of the design method presented in this paper. The proposed ITS system demonstrates features such as efficiency, accuracy, and energy control for the application in greenhouse.

Conclusions
In this work, a high-efficiency ITS system for greenhouse applications based on our proposed ITS efficiency formula was designed and validated. Designers can significantly reduce the range of critical components such as drive systems, energy supply systems, ITS mechanical frames, etc., by shrinking complex requirements to three main conditions. This provides guidance for the rapid design of high efficiency ITS systems required for composite applications. First, we present a real greenhouse application scenario. Then, according to the scenario, the designing process of the efficiency formula is elaborated. Finally, the efficiency of the designed ITS system is proved by experiments. In the future, the authors will further study the efficient ITS design issues in complex application environments (not limited to greenhouse applications); in addition, multi-ITS collaboration is an important research direction to further improve work efficiency.