Open access peer-reviewed chapter

Path Tracking of a Wheeled Mobile Manipulator through Improved Localization and Calibration

Written By

Tao Song, Fengfeng (Jeff) Xi and Shuai Guo

Submitted: 04 February 2018 Reviewed: 18 June 2018 Published: 05 November 2018

DOI: 10.5772/intechopen.79598

From the Edited Volume

Applications of Mobile Robots

Edited by Efren Gorrostieta Hurtado

Chapter metrics overview

1,530 Chapter Downloads

View Full Metrics


This chapter focuses on path tracking of a wheeled mobile manipulator designed for manufacturing processes such as drilling, riveting, or line drawing, which demand high accuracy. This problem can be solved by combining two approaches: improved localization and improved calibration. In the first approach, a full-scale kinematic equation is derived for calibration of each individual wheel’s geometrical parameters, as opposed to traditionally treating them identical for all wheels. To avoid the singularity problem in computation, a predefined square path is used to quantify the errors used for calibration considering the movement in different directions. Both statistical method and interval analysis method are adopted and compared for estimation of the calibration parameters. In the second approach, a vision-based deviation rectification solution is presented to localize the system in the global frame through a number of artificial reflectors that are identified by an onboard laser scanner. An improved tracking and localization algorithm is developed to meet the high positional accuracy requirement, improve the system’s repeatability in the traditional trilateral algorithm, and solve the problem of pose loss in path following. The developed methods have been verified and implemented on the mobile manipulators developed by Shanghai University.


  • mobile manipulator
  • localization
  • tracking
  • path following

1. Introduction

Recently, mobile manipulators have been used in various industries including aerospace or indoor decoration engineering, which requires a large workspace [1, 2, 3, 4]. The said mobile manipulator consists of an industrial manipulator mounted on a mobile platform to perform various manufacturing tasks such as drilling/riveting in aerospace industry or baseline drawing in decoration engineering. Wheeled mobile platforms with Mecanum wheels that can easily navigate through crowded spaces due to their omnidirectionality with a zero turning radius are commonly used. Path tracking is one of the important issues for mobile manipulators, in particular for performing manufacturing tasks. This chapter addresses this issue from the aspect of localization and calibration.

Localization is a key functionality of mobile manipulator in order to track and determine its position around the environment [5]. Many methods are proposed to address this issue [6]. They can be divided into two categories: absolute localization and relative localization.

Absolute localization relies on detecting and recognizing different features in the environment to obtain the position and posture. The features can be normally divided into two types: artificial landmarks and natural landmarks. Compared with the natural landmarks, artificial landmarks have advantages of high recognition, which will lead to high accuracy. There is no cumulative error problem when a localization method based on artificial landmarks is used. The key challenge is to identify and extract the needed information from the raw data of landmarks. For the relative localization, dead reckoning and inertial navigation are commonly carried out to obtain the systems’ position. It does not have to perceive the external environment, but the drift error accumulates over time.

Researchers have proposed several solutions, such as fuzzy reflector-based localization [7] and color reflector-based self-localization [8]. The main drawbacks of these two methods are that the anti-interference ability is poor and the computation is huge. To solve these problems pertaining to the localization method based on artificial landmarks, Madsen and Andersen [9] proposed a method using three reflectors and the triangulation principle. Betke and Gurvits [10] proposed a multichannel localization method with the three-dimensional localization principle and the least squares method. Because of the unavoidable errors in position and angle measurement of reflectors, the use of only a trilateral or triangular method will not achieve high accuracy [11]. Nevertheless, there are still many challenges for mobile manipulator working in industrial environments such as aerospace manufacturing or decoration engineering, which requires high maneuverability and high accuracy at the same time. The stationary industrial manipulator has high repeated localization accuracy, which the mobile manipulator cannot achieve. This chapter focuses on the improvement of path tracking of a mobile manipulator through enhanced localization combined with calibration. Calibration is required for the system kinematic model established based on nominal geometry parameter to improve motion accuracy. Muir and Neuman [12] proposed a kinematic error model for Mecanum wheeled platform and applied actuated inverse and sensed forward solutions to the kinematic control. Wang and Chang [13] carried out error analysis in terms of distribution among Mecanum wheels. Shimada et al. [14] presented a position-corrective feedback control method with vision system on Mecanum wheeled mobile platform. Qian et al. [15] conducted a more detailed analysis on the installation angle of rollers. An improved calibration method is presented in this chapter to improve the tracking accuracy of a mobile manipulator.


2. System modeling

The wheeled mobile manipulator, as shown in Figure 1, is built with a manipulator onto a wheeled mobile platform with four Mecanum wheels. This system aims to carry out fuselage drilling/riveting tasks at assembly stations in an adaptive and flexible way in the aerospace industry.

Figure 1.

The wheeled mobile manipulator for fuselage drilling/riveting.

This system needs to localize in real time during machining process. The global frame {XG,YG,ZG} is attached to the ground in the environment. The platform frame {XP,YP,ZP} is attached to the center of the mobile platform. The tool frame XMYMZM is attached to the tool of the manipulator. Two laser range finders are equipped and their frames XL1YL1ZL1 and XL2YL2ZL2 are attached to the left-front and right-rear of the mobile platform, respectively. The vision frame XVYVZV is attached to the industrial camera.

The position and posture of the system in the global frame can be defined as


where xG and yG are the positions in two directions, respectively, and θG is the azimuth angle, as shown in Figure 2.

Figure 2.

The position and posture of system.


3. Accuracy analysis of a wheeled mobile manipulator

3.1. Problem formulation

Figure 3 shows the kinematic model of the mobile platform with Mecanum wheels. According to kinematic analysis [16], the motion equation can be obtained as


Figure 3.

The kinematic model of a Mecanum wheeled mobile platform.

where V=VXVYω0T, D0=14R0R0R0R0R0R0R0R0R0l0+LR0l0+LR0l0+LR0l0+L, W=ω1ω2ω3ω4T.

VXVYω0T is defined as the linear and angular speeds of the platform; L is the half distance from the front axle to the rear axle as shown in Figure 3; l0 is the transverse distance from the wheel centers to the platform center line; R0 is the radius of the Mecanum wheel; andω1, ω2, ω3, and ω4 are angular velocities of the four wheels, respectively.

The Mecanum wheel and its roller are shown in Figure 4. The roller is fitted on the edge of the wheel at a certain angle (generally 45°) with its central axis. Each roller can rotate freely around the central axis. The wheel relies on the friction between the roller and the ground to move. The material of roller’s outer rim is usually rubber, which will deform under pressure between ground and the wheel.

Figure 4.

Mecanum wheel and roller.

Figure 5 shows the distribution of roller deformation. F and T are the force and driving torque on the roller, respectively. The radiuses of wheels reduce differently under different pressures. The deformation zone and the position of the wheel center change with the action of driving torque and the shifting [17]. Furthermore, with such deformation, R0 and l0 will change.

Figure 5.

Roller deformation.

As shown in Figure 6, based on the kinematic analysis and consideration on deformation of the roller, Eq. (2) can be revised as follows: [18].


where D1=14R1R2R3R4R1R2R3R4R1l1+LR2l2+LR3l3+LR4l4+L; here D1 is defined by R1, R2, R3, and R4, which are individual radiuses of the four Mecanum wheels, respectively. In order to improve the system motion accuracy, the relative errors Ri between Ri and R0, and relative errors li between li and l0 (i=1,2,3,4) should be obtained to revise matrix D1 (Figure 6).

Figure 6.

Motion model of the mobile platform.

3.2. Error modeling

By multiplying time t, Eq. (3) becomes a displacement equation as


where X=Vt is XYθT. θ=Wt is θ1θ2θ3θ4T. Individual geometric parameters including l1, l2, l3, l4, R1, R2, R3, and R4 are variables in Eq. (4). The relative errors can be determined from


leading to


where ∆X=∆X∆Y∆θT, ∆B=R1R2R3R4l1l2l3l4T, and


With the least squares method, the geometric errors can be solved as


Ri is generally defined as the tolerance resulting from manufacturing or assembly. With the deformation of the roller and li limited to hh and jj (h=3mm, j=5mm) respectively, li and Ri can be defined as:


where li and Ri are in l0jl0+j and R0hR0+h, respectively, and r is a random number between 0 and 1. The angular speeds of all wheels are set to 20rad/s, and the time t is set to 1 s.

The rank of matrix G in Eq. (7) is generally 3, so G is a full rank matrix and ∆B can be obtained. The following steps can be carried out. First, the displacement error measurement is analyzed. In order to obtain ∆B, it is needed to obtain the displacement error matrix ∆X first.

The system is moved in the YG direction and stopped every 2 s to obtain its position and posture XkGYkGθkG. The principle of measurement is shown in Figure 7 and X can be calculated as


Figure 7.

The principle of measuring displacement.

According to Eq. (2), XT, YT, and θT are theoretical values which can be obtained as follows: XT=Vxt, YT=Vyt, θT=0. Finally, through the experiment, the displacement errors can be obtained.


Substituting Eq. (11) in to Eq. (7), ∆B can be determined.

The Monte Carlo analysis is applied to deal with 50 samples. While X has been given in Eq. (11), the corresponding B should also satisfy its own tolerance, i.e., 3Ri3 and 5Li5.

The results are given in Figure 8. These data are averaged to form a new result: B=0.5411.2370.6051.0550.4580.6830.0480.683T. A process capability index is then used to estimate the value of B. Process capability refers to the ability of a process to meet the required quality. It is a measure of the minimum fluctuation of the internal consistency of the process itself in the most stable state. When the process is stable, the quality characteristic value of the product is in the range μ3σμ+3σ, where μ is the ensemble average of the quality characteristic value of the product, and σ is the standard deviation of the quality characteristic value of the product.

Figure 8.

The values obtained with the Monte Carlo method.

There are two kinds of situations for which the process capability index has to be solved. First, μ=M as shown in Figure 9: USL is the maximum error of quality and Lsl is the minimum error of quality; M=USL+LSL/2; μ is the average value of process machining; and Cp indicates the process capability index, and Cp=USLLSL/6σ.

Figure 9.


Second, μM as shown in Figure 10: Cpk indicates the process capability index, and Cpk=USLLSL/6σMμ/3σ. Only process capability index greater than 1 is valid. As 3∆Ri3 and 5∆Li5, M=0 and μM, the following results can be obtained: Cpk∆R1=1.58, Cpk∆l1=0.57, Cpk∆R2=1.06, Cpk∆l2=0.58, Cpk∆R3=1.54, Cpk∆l3=0.63, Cpk∆R4=1.17, andCpk∆l4=0.5.

Figure 10.


The values of R1, R2, R3, and R4 are closer to the real values than the values of l1, l2, l3, and l4. Now, taking B1a=0.5411.2370.6051.0550000T, B1b=00000.4580.6830.0480.683T in Eq. (5), it can be seen that X1=0.4119,40.0003T, X1b=000.0003T. According to the change rate relative to Eq. (10), the influence of R1,R2, R3, and R4 is much bigger than that of l1, l2, l3, and l4. Although the values of l1, l2, l3, and l4 are not accurate enough for the real values, B1 is valid.

The interval analysis is also carried out. In Eq. (6), the value ofG is uncertain because of the change of ∆Ri and ∆li in 33 and 55, respectively. Now, the increment of ∆Ri and ∆li is set to 1 mm; there are 7 values of ∆Ri, and 11 values of ∆li, so there are 74×114 groups of combinations of matrix G. Take all the combinations of matrix G to Eq. (6) to obtain ∆B and exclude the cases for which ∆Ri and ∆li are not in 33 and 55, respectively. Thus, 277 groups of ∆B can be obtained. As shown in Figure 11, the average values of R1, R2, R3, R4, l1, l2, l3, and l4, lead to a new B2=0.2331.4670.9320.8240.20.2130.0680.039T, which is close to B1.

Figure 11.

The values obtained in interval analysis.

The same method is used to solve the process capability index of ∆B2. One obtains: Cpk∆R1=3.34, Cpk∆l1=0.68, Cpk∆R2=1.22, Cpk∆l2=0.59, Cpk∆R3=2.5, Cpk∆l3=0.63, CpkR4=1.73, andCpkl4=0.6. The values of R1, R2, R3, and R4 in ∆B2 are closer to the real ones than the values of l1, l2, l3, and l4 in ∆B2. As the influence of R1, R2, R3, andR4 is bigger than that of l1, l2, l3, and l4, ∆B2 is valid.

The result is verified as well. ∆B1 is used to revise the parameters of D1 in Eq. (2) as


By setting two sets of four wheel speeds W1=10π/2110π/2110π/2110π/21T and W2=4π/74π/74π/74π/7T, the displacement errors can be computed as S1=D1D0W2t and S2=D2D0W2t. The results of the two correction methods are almost identical in theory.

The experiment with four movements is shown in Figure 12(a), while the actual movement is shown in Figure 12(b). It can be found that S1 and S2 are close to the measured displacement errors, so both ∆B1 and ∆B2 are satisfied to revise the matrix D1.

Figure 12.

(a) Schematic diagram of measuring and (b) the actual measuring.


4. Localization

4.1. System configuration

The localization component of mobile manipulator includes two laser range finders and a number of reflectors that are of cylindrical shape placed in the environment. Each reflector is covered by a reflective tape with high reflectivity (Figure 13).

Figure 13.

Localization system for the mobile manipulator.

4.2. Dynamic tracking

To achieve accurate localization, the system should have the ability to perceive external information through extracting the reflector features. In this research, as shown in Figure 13, there are n reflectors, and the feature of each reflector Bii=12n is extracted from the raw data. The feature extraction algorithm consists of three steps: (i) filtering and clustering, (ii) identification, and (iii) feature extraction.

The first step is filtering and clustering. The raw data obtained by each finder are a set of discrete data sequences γλii=12n. γ is the distance from the target point to the finder. ∅ is the polar angle. λi is the intensity value of the ith data point. In the process of data analysis, the outlier points that are contaminated by noise will be filtered out.

The density of the collected data points is proportional to the distance from the target point to the laser range finder. To improve the efficiency of the feature extraction process, an adaptive clustering method is adopted as given by Eq. (12). Unless the distance between two data points is less than the threshold δ, these data points are clustered for one reflector.


where γi1 is the distance value of the i1th data point, is the angle resolution of the laser range finder, β is an auxiliary constant parameter, and σγ is the measurement error. The values of parameters σγ and β are given as 0.01 m and 10°, respectively.

The second step is identification. After clustering, the data set can be correlated to each observable reflector. Each reflector data set is then used to calculate the position of the reflector in the laser range finder frame [19]. Let λδ be the reflected intensity threshold and DDδD+Dδ be the diameter range of a reflector, where D is a nominal reflector diameter and Dδ is the tolerance. The values of λδ and Dδ are selected based on the actual situation. Considering that Wc represents a set after clustering, i.e., Wc=γλii=mn, for each reflector, the measured diameter Dc is calculated as Dc=γn2+γm22γnγmcosnm, where γm and γn are the beginning and end data of a reflector set, respectively. When the set Wc satisfies the following two conditions


then the set Wc is established for all reflectors.

The third step is feature extraction. The feature extraction algorithm of a reflector extracts its central position γL,BβL,B in the laser range finder frame XLYLZL, where γL,B and βL,B are the distance and azimuth of each reflector, respectively. In the process of extracting the feature of the circular reflector, only a small portion of the entire circle was scanned with noise, so the accuracy of the fitting result would be low if a general least square circle fitting method is used. Since the radius of the reflector is known, the known radius is used as a constraint to improve the accuracy.

First, the value of βL,B is obtained from nm consecutive data points of the reflector set


As shown in Figure 14, Mi represents the ith data and OL2B¯ is a line between the reflector center B and the laser range finder center OL2; the projected angle of the line OL2B¯ in the laser range finder frame is LβB, while the angle between the line OL2B¯ and the line OL2Mi¯ is θi. The distance from Mi to B is approximately the radius of the reflector, i.e., BMi¯=D/2.


Figure 14.

Extraction of the center of the cylindrical reflector.

Furthermore, as mentioned before, the position of the finder in the platform frame is xPL2yPL2. The position of reflector center B γPBβPB in the platform frame can be expressed as


The optimal triangulation localization algorithm based on angle measurement is carried out. A number of experiments were carried out on the reflector feature extraction algorithm before proposing a localization algorithm for the system. A single reflector was placed in the measuring range of the finder. The finder performed repeated scanning when it was placed at difference places. Two different feature extraction algorithms were tested to calculate the position of the reflector, and the results are shown in Figure 15. Algorithm A is the feature extraction algorithm proposed in this chapter, while algorithm B is the least square circle fitting method without a radius constraint. Apparently, the former one yields a better result.

Figure 15.

Effect diagram using different feature extraction algorithms.

After extracting the reflector center, the positions of all the reflectors G=γLBβBLii=12n can be obtained and used to calculate the distance measurement range and the angler measurement range Rβ of the reflector, as given below:


where γLBmaxG, γLBminG, βBLmaxG, βBLminG.

It can be seen from Figure 16 that the angle measurement accuracy is better than the distance measurement accuracy. Therefore, based on these results, this chapter proposes an optimal trilateral localization algorithm based on angle measurement. The idea is to use the azimuth angle βb of the reflector in the platform frame to find the optima γb by the cosine theorem. The global pose of the mobile manipulator is then calculated based on the triangular method. The algorithm details are given as follows.

Figure 16.

The position of single reflector in the laser range finder frame.

First, it is assumed that the finder can cover at least three reflectors at each position. After feature extraction, the positions of three reflectors B1, B2, and B3 are calculated as γPB1βB1P; γPR,B2βPR,B2; and γPR,B3βPR,B3, respectively. In the global frame, these positions can also be obtained as B1Gx1y1; B2x2y2G; and B3x3y3G. Since they are measured from the same finder, three circles must intersect at the same point OMPG, and the value of OMPGxMPGyMPG, represents the position of the mobile platform, as shown in Figure 17.

Figure 17.

Schematic diagram of the localization algorithm.

According to the cosine theorem, the relations between the above variables can be expressed by the following equations:


In Eq. (21), the known parameters βPB1, βPB2, and βPB3 are used to solve γPB1a, γPB2a, and γPB3a. Since the above equations are nonlinear, the steepest descent method is adopted. The three circle equations can be expressed as


The position OMPGxMPGyMPG of the system can be obtained in the global frame by solving the above equations. Since the actual position of the reflector in the global environment deviates from the theoretical position, these circles may not intersect at the same point. In order to minimize the difference between the calculated position of the system and the actual position, the least squares estimation principle is applied. Assuming that the coordinate of the reflector is BiGxiyii=12n, n is the number of reflectors detected by laser range finder. The position value OMPGxMPGyMPG of the system is calculated as




The posture of the system also includes an azimuth angle θG in the global frame. First, θiG is obtained from the ith reflector as


The azimuth angle θG of the system is the averaged value from all the reflectors, as in


The dynamic tracking algorithm is then carried out. Localization of the system is based on landmarks. The system needs to meet the following two conditions to complete real-time localization:

  1. The system requires real-time localization of the reflector in the environment.

  2. The localization system correctly matches the real-time observed reflectors.

During the matching process, the reflectors observed in real time are made to correspond to all reflectors in the last moment in the environment one by one to extract the effective reflectors [20]. The localization can be achieved.

During localization, owing to the fact that some of the reflectors are obscured by obstacles or confused with a highly reflective metal surface object, the loss of position of the system is observed. To address the above problems, a dynamic tracking algorithm is proposed as shown in Figure 18.

Figure 18.

Schematic diagram of the dynamic tracking algorithm.

After placing n reflectors in the global environment, the system actually observes q reflectors at the tth moment, and the coordinate value of the ith reflector in the platform frame is Mt,iγt,Biβt,Bi0iq. The sampling time of the laser range finder is 0.1 s. Therefore, the theoretical position value Mt1,j1γt1,Bjβt1,Bj of the jth reflector in the platform frame can be deduced by the position Ot1,Rxt1yt1 of the system at the t1th moment and the position OGjxjyj of the jth reflector in the global environment.

If the difference between the theoretical value Mt1,jγt1,Bjβt1,Bj and the observed value Mt,iγt,Biβt,Bi is less than η, then the ith observed reflector is an effective reflector and is considered matched with the reference reflector Bj as follows:


Therefore, a set A of effective reflectors at the tth moment can be obtained, and A=Mt,oγt,Boβt,BoMt,iγt,Boβt,Bii=01q. Taking the mobile manipulator’s moving speed into account, the value of η depends on the actual situation. Through the use of the optimal triangulation localization algorithm based on angle measurement, the pose Ot,Rxtyt of the mobile manipulator can be calculated at the tth moment.

4.3. Experimental verification

The experimental data are obtained by the LMS 100 laser range finder with a scanning range of 270° and an angular resolution of 0.25°. The experimental platform is shown in Figure 19. The outside of the reflector is wrapped by reflective tape. In the experimental environment, five reflectors are placed around the system. Their global coordinate values are (0,995); (0, 0); (0, 1774); (2905, −2449); and (3956, −2032), and the unit is mm, as shown in Figure 20.

Figure 19.

Experimental platform.

Figure 20.

Experimental environment.

The optimal triangulation method based on angle measurement is used for validation by the repeatability of the system. In the stationary state of the system, the environment is scanned by finders. Each result is indicated by a red dot in Figure 21. The repeatability obtained by the trilateral method is nearly 18 mm, while the repeatability of the optimal method is only 9 mm. It can be shown that the optimal method is better than the traditional method.

Figure 21.

Repeatability localization of the system at the same location.

The mobile manipulator moves in the direction of the arrow in Figure 19, and each time the system moves a certain distance, the localization system will perform an experiment, i.e., it will use the left rear finder to calculate the current position. An average of 30 samples is taken for each experiment.

Figure 22 shows the results of static localization accuracy. The maximum distance error is 18 mm and the maximum angle error is 2°, which satisfies the localization requirement of the system.

Figure 22.

System localization error.

The mobile manipulator moves in the designated route, and it needs to constantly calculate and record its own position in the moving process. As shown in Figure 23, the trajectory of the moving system based on the localization method is smoother.

Figure 23.

Tracking results.

This chapter demonstrates the feasibility of a tracking and localization algorithm for mobile manipulators. The following conclusions can be drawn from this study: (i) In the detection of a reflector in the laser range finder frame, the angle repeatability of the reflector is better than that of the distance repeatability based on the feature extraction algorithm; (ii) The repeatability localization accuracy using the optimal triangulation method based on the angle measurement is nearly 9 mm, which is better than that of the trilateral method; (iii) The localization error of the system is 18 mm, which satisfies the localization requirement of system. Improvements in the location method based on reflectors, such as optimizing the layout of reflectors and the map of reflectors selection strategy for localization, are still needed.


5. Summary

In this chapter, through analyzing the roller deformation of the Mecanum wheel, the changed parameters of the motion equation of mobile system are found. The relative variation of the parameters in the motion equation of the Mecanum motion platform is solved by Monte Carlo analysis and interval analysis. Using the relative variation of the parameters to revise the motion equation, the displacement errors in different spaces in theory are solved for and compared with the measured displacement errors. From the comparison, both the methods are found to satisfy the system’s requirement. Then, the feasibility of a tracking and locating algorithm for mobile manipulator is demonstrated. The following conclusions can be drawn from this study: (i) In the detection of a reflector in the laser range finder frame, the angle repeatability of the reflector is better than that of the distance repeatability based on the feature extraction algorithm; (ii) The repeatability localization accuracy using the optimal triangulation method based on the angle measurement is nearly 9 mm, which is better than that of the trilateral method; (iii) The localization error of the system is 18 mm, which satisfies the localization requirement of system. Improvements in the localization method based on reflectors, such as optimizing the layout of reflectors and the map of reflectors selection strategy for localization, are still needed.

The method in this chapter is also used in the research of MoMaCo (Mobile manipulator in Construction), which can draw baseline for architectural decoration engineering as shown in Figure 24. The application result also verified the effectiveness of the method.

Figure 24.

MoMaCo (mobile manipulator in construction).


  1. 1. Guo S, Tao S, Xi F(J), Mohamed RP. Tip-over stability analysis for a wheeled mobile manipulator. Journal of Dynamic Systems, Measurement, and Control. 2017;139:054501
  2. 2. Tao S, Xi F(J), Guo S, Lin Y. Optimization of a mobile platform for a wheeled manipulator. Journal of Mechanisms and Robotics. 2016;8:061007
  3. 3. Tao S, Xi F(J), Guo S, Tu X, Li X. Slip analysis for a wheeled mobile manipulator. Journal of Dynamic Systems, Measurement, and Control. 2017;140(2):021005-021005-12
  4. 4. Guo S, Jin Y, Bao S, et al. Accuracy analysis of omnidirectional mobile manipulator with Mecanum wheels. Advances in Manufacturing. 2016;4:363
  5. 5. Guo S, Fang T-T, Tao S, Xi F-F, Wei B-G. Tracking and localization for omnidirectional mobile industrial robot using reflectors. Advances in Manufacturing. 2018;6(1):118-125
  6. 6. Larsson U, Forsberg J, Wernersson A. Mobile robot localization: Integrating measurements from a time-of-flight laser. IEEE Transactions on Industrial Electronics. 1996;43(3):422-431
  7. 7. Buschka P, Saffiotti A, Wasik Z. Fuzzy landmark-based localization for a legged robot. In: IEEE/RSJ International Conference on Intelligent Robots and Systems, vol. 2; 2000. pp. 1205-1210
  8. 8. Jang G, Lee S, Kweon I. Color landmark based self-localization for indoor mobile robots. In: IEEE International Conference on Robotics & Automation; 2002. pp. 1037-1042
  9. 9. Madsen CB, Andersen CS. Optimal reflector selection for triangulation of robot position. Robotics and Autonomous Systems. 1998;23(4):277-292
  10. 10. Betke M, Gurvits L. Mobile localization using landmarks. IEEE Transactions on Robotics and Automation. 1997;13(2):251-263
  11. 11. Kai OA, Tomatis N, Jensen BT, et al. Multi-laser range finder on-the-fly localization: Precision and reliability for applications. Robotics and Autonomous Systems. 2001;34(2–3):131-143
  12. 12. Muir PF, Neuman CP. Kinematic modeling for feedback control of an omnidirectional wheeled mobile robot. In: IEEE International Conference on Robotics and Automation; 1987. pp. 1772-1778
  13. 13. Wang Y, Chang D. Preferably of characteristics and layout structure of Mecanum comprehensive system. Journal of Mechanical Engineering. 2009;45(5):307-310
  14. 14. Shimada A, Yajima S, Viboonchaicheep P, et al. Mecanum-wheel vehicle systems based on position corrective control. In: IEEE Annual Conference on Industrial Electronics Society; 2005
  15. 15. Qian J, Wang X, Xia G. Velocity rectification for mecanum wheeled omni-direction robot. Machine Tool and Manufacturing Technology. 2010;11:42-45.
  16. 16. Huang L et al. Design and analysis of a four-wheel omnidirectional mobile robot. In: 2nd International Conference of Autonomous Robots and Agents; 2004
  17. 17. Song JB, Byun KS. Design and control of a four-wheeled omnidirectional mobile robot with steerable omnidirectional wheels. Journal of Robotic Systems. 2004;21(4):193-208
  18. 18. Asama H et al. Development of an omni-directional mobile robot with 3 DOF decoupling drive mechanism. In: IEEE International Conference on Robotics and Automation
  19. 19. Liu YM. Three-dimensional optical positioning method research [dissertation]. Shenyang University of Technology; 2003
  20. 20. Jia W, Zhou W, Kaiser J. Efficient algorithm for mobile multicast using any cast group. IEE Proceedings: Communications. 2001;148(1):14-18

Written By

Tao Song, Fengfeng (Jeff) Xi and Shuai Guo

Submitted: 04 February 2018 Reviewed: 18 June 2018 Published: 05 November 2018