Next Article in Journal
Robust Pedestrian Dead Reckoning Based on MEMS-IMU for Smartphones
Previous Article in Journal
Optimal Sub-Band Analysis Based on the Envelope Power Spectrum for Effective Fault Detection in Bearing under Variable, Low Speeds
Previous Article in Special Issue
Evaluating Muscle Activation Models for Elbow Motion Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heading Estimation of Robot Combine Harvesters during Turning Maneuveres

by
Md Mostafizar Rahman
1,2 and
Kazunobu Ishii
3,*
1
Graduate School of Agriculture, Hokkaido University, Kita-9, Nishi-9, Kita-ku, Sapporo 060-8589, Japan
2
Department of Food Engineering and Tea Technology, Shahjalal University of Science and Technology, Sylhet 3114, Bangladesh
3
Research Faculty of Agriculture, Hokkaido University, Kita-9, Nishi-9, Kita-ku, Sapporo 060-8589, Japan
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(5), 1390; https://doi.org/10.3390/s18051390
Submission received: 9 March 2018 / Revised: 17 April 2018 / Accepted: 17 April 2018 / Published: 1 May 2018
(This article belongs to the Special Issue Smart Sensors for Mechatronic and Robotic Systems)

Abstract

:
Absolute heading is an important parameter for a robot combine harvester or a robot tracked combine harvester, especially while it is turning, but due to the rapid turning of robot combine harvesters, its inertial measurement unit gives a gyro measurement bias that causes heading drift. Our research goal is to estimate the absolute heading of robot combine harvesters by compensating this gyro measurement bias during non-linear turning maneuvers. A sensor fusion method like the extended Kalman filter combined with the tracked combine harvester dynamic model and sensor measurements was used to estimate the absolute heading of a robot combine harvester. Circular, sinusoidal and concave shapes were used to evaluate the estimated heading produced by the sensor fusion method. The results indicate that the estimated heading is better than measured heading which was calculated from the integration of yaw rate gyro measurements, and the root mean square errors (RMSEs) for estimated headings are smaller than the measured headings. In practics, the target of this paper is thus the estimation of a heading or absolute heading that is bias compensated, and can be further used to calculate the exact crop periphery for automatic path planning of robot combine harvesters.

1. Introduction

Unlike wheeled vehicles, tracked vehicles are widely popular due to their non-linear contact characteristics between the tracks and the ground, which allows them to be operated under adverse field conditions for agricultural production, and allows turning at high speed with a small turning radius or higher steering command. Under autonomous conditions, the tracked vehicle has a global positioning system and inertial sensor for providing the vehicle state and direction, but when this tracked vehicle is turning, the inertial sensor reading has measurement uncertainties. Since the vehicle direction is more important for autonomous guidance and other navigation purposes, it is necessary to compensate this uncertainty of the inertial sensor measurements.
The main controlling feature of auto guidance is to steer the vehicle to follow a desired path automatically, which requires a proper guidance system able to detect the vehicle’s position and direction, create proper steering signals, and steer the vehicle according to the signals [1].There are different guidance sensing systems, including global positioning system, inertial sensors, geomagnetic direction sensors, machine vison and laser scanners which are used to find the control parameters of the autonomous vehicles such as heading and offset [2,3,4,5]. The functional characteristics of each sensor provides the desired information, which contains erroneous measurement readings due to noise, measurement errors, and time delays. In general, a single sensor is not able to provide enough information, whereas multi-sensor integration can provide more useful information, which is more helpful and informative than what can be observed using a single sensor. This information needs to be fused in a way that reduces sensor uncertainties and the additional task of interpretation must be performed [6]. There are different approaches used for sensor fusion to obtain the position and direction of vehicles, such as complementary filters [7] or Kalman filters with various architectures [8,9,10,11,12,13], particle filters or sequential Monte Carlo methods [14,15,16,17,18]. Since the dynamic motion of a vehicle is non-linear, a non-linear dynamic model and extended Kalman filter are commonly used for navigation purposes. For instances, Rohac [19] addressed a cost-effective solution by using a non-linear observer and extended Kalman filter with commercial grade inertial sensors, and studied the performances of different approaches to obtain navigation solutions with robustness to GNSS outages. Noguchi [20] developed a guidance system based on the real time global positioning system (RTK-GPS), geomagnetic sensor (GDS) and machine vision by using an extended Kalman filter method which provided the most appropriate vehicle heading in real time. Alatise and Hancke [21] used an extended Kalman filter with inertial sensors and camera to estimate the position and direction of a mobile robot. With a non-linear extended Kalman filter, a tracked vehicle model and sensor measurements are used to estimate the trajectory, orientation and other soil parameters for small-scale tracked vehicles [22].
An autonomous tracked combine harvester developed by Zhang [23] was used to cut wheat and paddy rice in real time by using an appropriate harvesting map which is created based on a real time global positioning system (RTK-GPS) and inertial measurement units (IMUs); before harvesting, the outside crop near to headland is cut twice or thrice by the tracked combine harvester, and during this time the sensor measurements are logged for making a navigation map. Generally, the tracked combine harvester makes turns at moderate to high speed with a small turning radius, which is very popular with farmers, and this turning position is represented by a circle marked in Figure 1. During this turn, the measured heading of a tracked combine harvester contains drift errors, which result from the IMU gyro measurement bias. To get an absolute heading, this bias needs to be compensated by using a tracked combine harvester model and sensor fusion method. In this research, the main contributions are the estimated absolute heading of tracked combine harvester during non-linear conditions (especially high speed turns with small turning radii) by using a tracked combine harvester dynamic model which was developed by us based on a real time global positioning system (RTK-GPS) and inertial measurement units (IMU). For more details readers should see [24]. In practice, this estimated heading can further be used to obtain the exact crop periphery for calculating the harvesting map of a robot combine harvester (for more details see [25]).

2. Materials and Methods

2.1. System Platform and Sensors

A model AG1100 tracked combine harvester (Yanmar Co., Ltd, Japan) was used as a prime research platform, which was equipped with a Real Time Kinematic Global Positioning System (RTK-GPS) and Inertial Measurement Unit (IMU) as shown in Figure 2. This tracked combine harvester is converted into a robot combine harvester which is fully controlled by a CAN bus, and the speed of the tracked combine harvester ranges from −2.0~2.0 m/s, approximately. This harvester can be used to harvest cereal crops like paddy rice and wheat. The VN-100 IMU (VectorNav Technologies, Dallas, Texas, USA) is a miniature, high performance Inertial Measurement Unit and Attitude Heading Reference System (AHRS) which is used to measure the heading and angular rate of the tracked combine harvester. Its dynamic accuracy for the yaw rate is 1.0° RMS, heading’s dynamic accuracy: 2.0° RMS and the gyro has an alignment error: ±0.05°, noise density: <0.0035°/s/√Hz, in-run bias stability: <10°/h. This IMU, consisting of three accelerometers, three magnetometers and three rate gyros, provides nine axes, which are sampled at a rate of 200 Hz through a USB serial port. On the other hand, the update rate of the RTK-GPS is up to 20 Hz, although for logging of the RTK-GPS output, the update rate is maintained at 5 Hz or 10 Hz and the Baud rate was 115,200, where the IMU also matches that update rate. The RTK correction signal was calculated from a Virtual Reference Station (VRS) via an Internet connected to the on-board computer that logs the data from the RTK-GPS receiver through an RS232C serial port. These sensor measurements were logged into the control PC of the tracked combine harvester as a GGA header for RTK-GPS and VNIMU header for IMU. The GGA header from the RTK-GPS sensor gives the position of the tracked combine harvester, whereas the VNIMU header provides bias-uncompensated IMU measurements such as magnetic, acceleration and angular rate measurements, temperature and pressure. In this research, only position and bias uncompensated angular rate data were considered to estimate the absolute heading of robot combine harvester under non-linear conditions.

2.2. Tracked Combine Harvester Model

Figure 3 shows the free body diagram of dynamic model for the tracked combine harvester moving on a general plane, turning to the left or counter-clockwise [24,26,27]. Its acceleration is in the positive x c ,   y c and φ directions. The external thrusts and resistive forces acting on the tracked combine harvester are   F R , F L and R R , R L , respectively. The value f y indicates the lateral friction force due to the effect of lateral soil shear. Figure 3 is represented in the global reference frame XYZ, which indicates the tracked combine harvester turns around an instantaneous center of rotation (ICR). The angle β is called sideslip angle and is determined from the velocity V c and the longitudinal axis x c of the tracked combine harvester. It is assumed that the normal pressure distribution along the track is non-uniform, and the coefficient of lateral resistance μ is not constant. The instantaneous center of rotation must shift forwards of the tracked combine harvester centroid by the amount of D, as shown in Figure 3. This longitudinal shifting D depends on the tracked combine harvester lateral acceleration [27]. D is required to develop a net lateral force that accelerates the tracked combine harvester towards the instantaneous center of rotation, and also minimizes the resistive yawing moment [28]. The centrifugal force F c acting on the tracked combine harvester is given in Figure 3. We also described this tracked combine harvester model based on sensor measurements so that it can be further used for autonomous navigation purposes [24].
For a tracked combine harvester of mass m and a moment of inertia about the center of mass I , the equations of motion can be written in the body reference frame by using Equations (1)–(3), respectively:
m x ¨ c = F R + F L R R R L F c sin β
m y ¨ c = F c cos β μ m g
I φ ¨ = [ ( F R R R ) ( F L R L ) ] B 2 M r
Using a kinematic model of the tracked combine harvester, the velocity V c is calculated based on each track’s velocity including the slip of tracks from the tracked combine harvester expressed by Equation (4). The RTK-GPS positions used to calculate the actual turning radius R by using the least square method [24] is given by Equation (6), which was derived from the general circle Equation (5), where ( a , b ) and C are the center circle and coefficient:
V c = V R ( 1 S r ) + V L   ( 1 S l ) 2
( a x ) 2 + ( y b ) 2 = R 2
R = ( a 2 + b 2 C )
Now, Equations (1) and (2) are integrated to transfer the longitudinal and lateral direction velocities of the tracked combine harvester to the local coordinates. In order to run this tracked combine harvester in real time, the velocities in the longitudinal and lateral direction in the harvester coordinates are expressed in a global reference frame by the Equation (7):
[ X ˙ Y ˙ ] = [ cos φ sin φ sin φ cos φ ] [ x ˙ c y ˙ c ]
By using the Equation (7), the tracked combine harvester model equations can be rewritten as Equation (8):
[ X ˙ c Y ˙ c φ ˙ ] = [ ( 1 m ( F R + F L R R R L F c sin β ) cos φ 1 m ( F c cos β μ m g ) sin φ ) Δ t ( 1 m ( F R + F L R R R L F c sin β ) sin φ + 1 m ( F c cos β μ m g ) cos φ ) Δ t ( [ ( F R R R ) ( F L R L ) ] B 2   M r ) Δ t I ]
The discrete-time process model is developed from the continuous-time process model, where the state space equation is obtained by integrating the continuous-time equations over the interval from t k to t k + 1 . The discrete-time model is an approximation of the continuous-time model. Now, the discretization of the Equation (8) can be written by the following Equation (9):
[ X ˙ k + 1 X k + 1 Y ˙ k + 1 Y k + 1 φ ˙ k + 1 φ k + 1 ] = [ X ˙ k + ( 1 m ( F R + F L R R R L F c sin β ) cos φ 1 m ( F c cos β μ m g ) sin φ ) Δ t X k + X ˙ k Δ t + ( 1 m ( F R + F L R R R L F c sin β ) cos φ 1 m ( F c cos β μ m g ) sin φ ) Δ t 2 2 Y ˙ k + ( 1 m ( F R + F L R R R L F c sin β ) sin φ + 1 m ( F c cos β μ m g ) cos φ ) Δ t Y k + Y ˙ k Δ t + ( 1 m ( F R + F L R R R L F c sin β ) sin φ + 1 m ( F c cos β μ m g ) cos φ ) Δ t 2 2 φ ˙ k + ( [ ( F R R R ) ( F L R L ) ] B 2   M r ) Δ t I φ k + φ ˙ k Δ t + ( [ ( F R R R ) ( F L R L ) ] B 2   M r ) Δ t 2 2 I ]

2.3. RTK-GPS and IMU Fusion Algorithm

During the moderate to high speed turns of the tracked combine harvester, the yaw rate measurement from the IMU sensor will give a bias error that creates a heading drift error, which requires compensation for estimating the absolute heading. For this reason, the model given in Equation (9) was used to estimate the heading of the tracked combine harvester by using a non-linear extended Kalman filter. According to the given Equation (9), the state vector of tracked combine harvester was defined as the Equation (10):
X k + 1 = [ X k + 1 , Y k + 1 , X ˙ k + 1 , Y ˙ k + 1 ,   φ ˙ k + 1 ,   φ k + 1 , φ g k + 1 , b k + 1 ] T
where:
  • X k + 1 = Tracked combine harvester position in the east direction at time t k + 1 ;
  • Y k + 1 = Tracked combine harvester position in the north direction at time t k + 1 ;
  • X ˙ k + 1 = Tracked combine harvester velocity in the X-direction at time t k + 1 ;
  • Y ˙ k + 1 = Tracked combine harvester velocity in the Y-direction at time t k + 1 ;
  • φ ˙ k + 1 = Yaw rate at time t k + 1 ;
  • φ k + 1 = Heading of tracked combine harvester at time t k + 1 ;
  • φ g k + 1 = Heading from RTK-GPS at time t k + 1 ; and
  • b k + 1 = Yaw rate gyro measurement bias at time t k + 1 .
The system measurement vector is defined by the following Equation (11):
Z k + 1 = [ X g p s , Y g p s ,   φ ˙ i m u ,   φ g p s ] T
where:
  • X g p s = Measurement position in the east direction from RTK-GPS;
  • Y g p s = Measurement position in the north direction from RTK-GPS;
  • φ ˙ i m u = Yaw rate gyro measurement from the IMU; and
  • φ g p s = Measurement heading from the RTK-GPS.
Now, the extended Kalman filter equations in the following Equations (12) and (13) are considered for the prediction and correction steps; which is maintained in a cyclic manner. The initial state covariance matrix P k , process noise covariance matrix Q k = d i a g . [ 0.002909839 0.002909839 0 0 0.000781853 0.001350055 0.00508 0 ] and measurement noise covariance matrix R k = d i a g . [ 0.00000297878 0.00000297878 0.000000438461 0.010348 ] are used for this extended Kalman filter, which were obtained by the static and dynamic test of RTK-GPS and IMU sensors:
Prediction:
Predicted   state , X k + 1 = f ( X k , u k ) Predicted   covariance , P k + 1 = A k P k A k T + Q k Predicted   measurement , Z k + 1 = h ( X k + 1 )   }
Correction:
Estimate   state , X ^ k + 1 = X k + 1 + K v k Estimate   covariance ,   P ^ k + 1 = P k + 1 ( I K H k ) }
where:
  • Innovation, v k = Z k + 1   h ( X k + 1 )
  • Kalman gain, K = P k + 1 H k T S 1
  • Innovation covariance, S = H k P k + 1 H k T + R k
The Jacobian matrices of A k and H k in the prediction step are the following, in which A k was obtained from the partial derivatives of each state vector by using the Equation (9), and H k was obtained from the partial derivatives of each measurement vector.
A k = f x = [ f 1 x 1 f 1 x m f n x 1 f n x m ]   and ,   H k = h x = [ h 1 x 1 h x m h n x 1 h n x m ]

3. Results and Discussion

3.1. Trajectory of Tracked Combine Harvester

For the evaluation of the estimated heading for the tracked combine harvester by compensating the yaw rate gyro measurement bias under non-linear conditions, the tracked combine harvester was run in circular and sinusoidal shapes, and during this time, the sensor measurements were logged, and the running trajectory of the tracked combine harvester can be obtained by using RTK-GPS position. Figure 4 shows the circular trajectories of the tracked combine harvester from the different input steering commands that were at 10° and 15°, respectively. On the other hand, the sinusoidal trajectories of the tracked combine harvester were made by a series of input steering commands such as ± 20 ° and ± 40 °, respectively, as shown in Figure 5. Based on these circular and sinusoidal trajectory sensor measurements, the absolute heading of the tracked combine harvester was estimated by the extended Kalman filter combined with the tracked combine harvester model.

3.2. Estimated Heading of Circular Trajectory

Figure 6 indicates the measured and estimated heading of circular trajectories (left figures) as well as the difference of these headings from the reference GPS heading (right figures) which was obtained by the linear regression of the noisy trend GPS heading. The heading from GPS measurements is indicated by a blue line, where the measured and estimated headings are marked by the black and red lines. In general, a full circle rotation is counted to 360°, and when the tracked combine harvester was run in the field with a constant steering command, therefore, a full trajectory of circle will be 360°. That why, the left Figure 6a,b indicate that the output headings were bounded by 360°. According to these figures, the red line tried to follow the blue one, which means the estimated heading based on the extended Kalman filter was consistent with the GPS heading. In addition, the heading difference analysis is shown in right Figure 6a,b, which indicates the difference between the measured and estimated headings from the reference GPS heading. Since, the GPS heading provides the exact direction of the vehicle that is computed from the exact position of the tracked combine harvester, the linear regression of the GPS heading was used as a reference GPS heading for the evaluation of estimated headings. The result indicated out that the estimated heading (red line) matches the reference GPS heading for all circular trajectories as shown in Figure 6a,b (right figures) rather than the measured heading.
From the error analysis of these headings, the RMS errors between the measured or estimated heading and the reference GPS heading were calculated. The RMS errors for the measured heading are 6.74 for 10° steering angle and 4.18 for 15° steering angle, whereas the RMS errors for estimated heading are 1.93 for 10° steering angle and 1.58 for 15° steering angle, respectively. The RMS errors for the estimated heading are lower than the measured heading for all circular trajectories, indicating that the heading drift errors caused from the yaw rate gyro measurement bias are compensated.

3.3. Estimated Heading of Sinusoidal Trajectory

Figure 7 shows the measured and estimated heading (left figures) and heading difference between the measured or estimated heading and the non-linear regression of the GPS heading (right figures). The purple color line in the left figures indicates the non-linear regression of the GPS heading called the reference GPS heading which was obtained by using the Gauss-Newton algorithm. The Gauss-Newton algorithm is a simple method for solving any non-linear regression problem, and the equation for non-linear regression analysis was considered as Y i = A s i n ( X i T + θ ) + a X i + b , where A , T ,   θ , a and b are the amplitude, period, phase shift angle, phase shift (horizontally) and phase shift (vertically), respectively. On the other hand, the measured and estimated headings were represented by the black and red color lines. According to Figure 7a,b (left side), the estimated heading is better than the measured heading and it tried to match with the reference GPS heading. The heading difference for the estimated heading and measured heading was obtained from the reference GPS heading as shown in the right Figure 7a,b. From the error analysis of measured and estimated heading regarding the reference GPS heading, the RMS errors for the measured heading are 2.39 for ±20° steering angle and 3.96 for ±40° steering angle, whereas the RMS errors for estimated heading are 2.29 for ±20° steering angle and 3.48 for ±40° steering angle, respectively. The RMS errors for the measured and estimated heading do not have a big difference because the reference GPS heading is sometimes matched with the noisy trend GPS heading, and sometimes not matched. If the reference GPS heading matches with the exact GPS heading, the RMS error of measured and estimated heading will change, and the RMS error of the estimated heading will be lower than that of the measured heading.

3.4. Estimated Heading of Concave Polygon Field

The heading estimation method based on the tracked combine harvester dynamic model and the extended Kalman filter can also be evaluated for any polygon-shaped field, because this polygon field with crops will be harvested by the robot combine harvester, which is our main research target. Here, only a concave polygon field as shown in Figure 8 was used for this estimation. The circle marked on the concave polygon field (top figure) indicates the turning position of the tracked combine harvester. In these circle positions, the tracked combine harvester cuts the surrounding crop to calculate the exact crop periphery for obtaining the working path of the robot combine harvester by forward and backward movement, and sometimes turns at high speed with low turning radius are considered, which creates a yaw rate gyro measurement bias. This bias creates the measured heading drift error that is important to compensate. Figure 8 (bottom figure) shows a heading of a concave polygon field where the black and red lines are marked for the measured and estimated headings, which was analyzed with respect to the reference GPS heading (blue color line). This reference GPS heading was obtained from the moving average of the noisy trend GPS heading. The estimated heading is better matched with the reference GPS heading than the measured heading in Figure 8 (bottom figure). In addition, the heading difference between the measured heading or estimated heading and the reference GPS heading was calculated, as shown in Figure 9. The RMS error was calculated from the measured or estimated heading and the reference GPS heading, and from the error analysis, the RMS errors for measured and estimating heading were 14.93° and 5.17°, respectively. Since, the RMS error of estimated heading is lower than the measured heading, the estimated heading shows more consistency with the reference GPS heading. The results also describe that the measured heading drift error was compensated by the extended Kalman filter and the tracked combine harvester dynamic model.

4. Conclusions

This research deals with the tracked combine harvester dynamic model, sensor measurements and the extended Kalman filter, which are used to estimate the absolute heading of a tracked combine harvester under non-linear condition by compensating the IMU yaw rate gyro measurement bias which influences the heading drift errors. Different sets of experiments such as circular, sinusoidal and concave polygon trajectories were executed to evaluate this estimation method. The results for different trajectories revealed that the IMU yaw rate gyro measurement bias was compensated based on the extended Kalman filter and the tracked combine harvester dynamic model, and the absolute heading was simultaneously determined, which is better than the measured heading. The RMS errors of the estimated heading are lower than those of the measured heading for all trajectories. Therefore, this estimation method can be used to estimate the absolute heading when the tracked combine harvester will makes turn at high speed and with high order steering commands in order to cut wheat and paddy rice near headlands to calculate the exact crop periphery for the development of a harvesting map, for more details see [25].

Author Contributions

M.M.R. conducted the research idea, experiments, analysis the experimental results, writing, reviewing and editing the manuscript; K.I. provided the initial research concepts, supervision and reviewing the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Q.; Reid John, F.; Noguchi, N. Agricultural vehicle navigation using multiple guidance sensors. In Proceedings of the International Conference on Field and Service Robotics, Pittsburgh, PA, USA, 29–31 August 1999. [Google Scholar]
  2. Benson, E.R.; Reid, J.F.; Zhang, Q. Machine vision based steering system for agricultural combines. In Proceedings of the ASAE Annual Meeting, Sacramento, CA, USA, 31 July–2 August 2001. [Google Scholar]
  3. Choi, J. Development of Guidance System Using Local Sensors for Agricultural Vehicles. Ph.D. Thesis, Graduate School of Agriculture, Hokkaido University, Sapporo, Japan, 25 March 2014. Available online: http://hdl.handle.net/2115/56098 (accessed on 25 March 2014).
  4. Ishii, K.; Terao, H.; Noguchi, N. Studies on Self-learning Autonomous Vehicles (Part-1). J. Jpn. Soc. Agric. Mach. 1994, 56, 53–60. [Google Scholar]
  5. Kise, M.; Noguchi, N.; Ishii, K.; Terao, H. Field mobile robot navigated by RTK-GPS and FOG (Part-2). J. Jpn. Soc. Agric. Mach. 2001, 63, 80–85. [Google Scholar]
  6. Adla, R.; Yousef, B.; Al-Holou, N. Multi sensor Data Fusion, Methods and Problems. In Proceedings of the International Conference on Parallel and Distributed Processing Techniques and Applications, Las Vegas, NV, USA, 22–23 July 2013; Volume II. [Google Scholar]
  7. Reinstein, M.; Kubelka, V. Complementary filtering approach to orientation estimation using inertial sensors only. In Proceedings of the IEEE International Conference on Robotics and Automation (ICRA), St. Paul, MN, USA, 14–18 May 2012. [Google Scholar]
  8. Alam, M.; Moreno, G.; Sipos, M.; Rohac, J. INS/GNSS localization using 15 state extended Kalman filter. In Proceedings of the International Conference in Aerospace for Young Scientists, Beijing, China, 12–13 November 2016; pp. 425–435. [Google Scholar]
  9. Azam, S.E.; Chatzi, E.; Papadimitriou1, C.; Smyth, A. Experimental validation of the Kalman-type filters for online and real-time state and input estimation. J. Vib. Control 2017, 23, 2494–2519. [Google Scholar] [CrossRef]
  10. Rezaeian, A.; Zarringhalam, R.; Fallah, S.; Melek, W.; Khajepour, A.; Chen, S.-K.; Litkouhi, B. Cascaded Dual Extended Kalman Filter for Combined Vehicle State Estimation and Parameter Identification. In Proceedings of the SAE World Congress and Exhibition, Detroit, MI, USA, 16–18 April 2013. [Google Scholar]
  11. Simanek, J.; Reinstein, M.; Kubelka, V. Evaluation of the EKF-based estimation architectures for data fusion in mobile robots. IEEE ASME Trans. Mechatron. 2015, 20, 985–990. [Google Scholar] [CrossRef]
  12. Simon, D. Kalman filtering with state constraints: A survey of linear and nonlinear algorithms. IET Control Theory Appl. 2010, 4, 1303–1318. [Google Scholar] [CrossRef]
  13. Zihajehzadeh, S.; Loh, D.; Lee, T.J.; Hoskinson, R.; Park, E.J. A cascaded Kalman filter-based GPS/MEMS-IMU integration for sports applications. Measurement 2015, 73, 200–210. [Google Scholar] [CrossRef]
  14. Gustafsson, F.; Gunnarsson, F.; Bergman, N.; Forssell, U.; Jansson, J.; Karlsson, R.; Nordlund, P.-J. Particle filters for positioning, navigation and tracking. IEEE Trans. Signal Process. 2002, 50, 425–437. [Google Scholar] [CrossRef]
  15. Martino, L.; Read, J.; Elvira, V.; Louzada, F. Cooperative parallel particle filters for online model selection and applications to urban mobility. Digit. Signal Process. 2017, 60, 172–185. [Google Scholar] [CrossRef]
  16. Míguez, J. Analysis of selection methods for cost-reference particle filtering with applications to maneuvering target tracking and dynamic optimization. Digit. Signal Process. 2007, 17, 787–807. [Google Scholar] [CrossRef]
  17. Sotak, M.; Sopata, M.; Kmec, F. Navigation systems using Monte Carlo method. In Proceedings of the 6th International ESA Conference on Guidance, Navigation and Control Systems, European Space Agency, Loutraki, Greece, 17–20 October 2005. [Google Scholar]
  18. Urteaga, I.; Bugallo, M.F.; Djuri’c, P.M. Sequential Monte Carlo Methods under model uncertainty. In Proceedings of the IEEE Statistical Signal Processing Workshop (SSP), Palma de Mallorca, Spain, 26–29 June 2016. [Google Scholar]
  19. Rohac, J.; Hansen, J.M.; Alam, M.; Sipos, M.; Johansen, T.A.; Thor, I.; Fossen, T.I. Validation of nonlinear integrated navigation solutions. Annu. Rev. Control 2017, 43, 91–106. [Google Scholar] [CrossRef]
  20. Noguchi, N.; Reid, J.F.; Benson, E.; Will, J.; Stombaugh, T. Vehicle Automation System Based on Multi-Sensor Integration; Paper 983111; ASAE (American Society of Agricultural Engineers): St. Joseph, MI, USA, 1998. [Google Scholar]
  21. Alatise, M.B.; Hancke, G.P. Pose estimation of a mobile robot based on fusion of IMU data and vision data using an extended Kalman filter. Sensors 2017, 17, 2164. [Google Scholar] [CrossRef] [PubMed]
  22. Dar, T.M.; Longoria, R.G. Slip estimation for small scale robotic tracked vehicles. In Proceedings of the American Control Conferencem, Marriott Waterfront, Baltimore, MD, USA, 30 June–2 July 2010. [Google Scholar]
  23. Zhang, Z. Development of a Robot Combine Harvester Based on GNSS. Ph.D. Thesis, Graduate School of Agriculture, Hokkaido University, Sapporo, Japan, 25 March 2014. Available online: http://hdl.handle.net/2115/56100 (accessed on 25 March 2014).
  24. Rahman, M.M.; Ishii, K.; Noguchi, N. Study on Tracked Combine Harvester Dynamic Model for Automated Navigation Purposes. Adv. Robot. Autom. 2017, 6, 178. [Google Scholar] [CrossRef]
  25. Rahman, M.M.; Ishii, K.; Noguchi, N. Automatic Path Planning of Robot Combine Harvester based on Optimum Harvesting Area of Convex and Concave Polygon Field. IFAC J. Syst. Control. under review.
  26. Kitano, M.; Jyozaki, H.A. Theoretical analysis of steerability of tracked vehicles. J. Terramech. 1976, 13, 241–258. [Google Scholar] [CrossRef]
  27. Wong, J.Y. Theory of Ground Vehicles; A Wiley-Interscience Publication: New York, NY, USA, 2008. [Google Scholar]
  28. Le, A.T.; Rye David, C.; Durrnat-Whyte Hugh, F. Estimation of track-soil interactions for autonomous tracked vehicles. In Proceedings of the IEEE International Conference on Robotics and Automation, Albuquerque, NM, USA, 20–25 April 1997; pp. 1388–1393. [Google Scholar]
Figure 1. A polygonal field with a circle representing the turning area of the tracked combine harvester.
Figure 1. A polygonal field with a circle representing the turning area of the tracked combine harvester.
Sensors 18 01390 g001
Figure 2. View of the robot combine harvester equipped with RTK-GPS and IMU sensors.
Figure 2. View of the robot combine harvester equipped with RTK-GPS and IMU sensors.
Sensors 18 01390 g002
Figure 3. Free body diagram of the tracked combine harvester dynamic model.
Figure 3. Free body diagram of the tracked combine harvester dynamic model.
Sensors 18 01390 g003
Figure 4. Circular trajectories of the tracked combine harvester at 10° and 15° input steering angles.
Figure 4. Circular trajectories of the tracked combine harvester at 10° and 15° input steering angles.
Sensors 18 01390 g004
Figure 5. Sinusoidal trajectories of the tracked combine harvester at ±20° and ±40° input steering angles.
Figure 5. Sinusoidal trajectories of the tracked combine harvester at ±20° and ±40° input steering angles.
Sensors 18 01390 g005
Figure 6. Measured and estimated headings (left figures) for circular trajectories as well as heading difference (right figures) (where, GPSH = GPS Heading, MH = Measured Heading, EH = Estimated Heading and GPSH (Reg) = Linear Regression of GPS Heading).
Figure 6. Measured and estimated headings (left figures) for circular trajectories as well as heading difference (right figures) (where, GPSH = GPS Heading, MH = Measured Heading, EH = Estimated Heading and GPSH (Reg) = Linear Regression of GPS Heading).
Sensors 18 01390 g006
Figure 7. Measured and estimated headings (left figures) for sinusoidal trajectories as well as heading difference (right figures) (where, GPSH = GPS Heading, MH = Measured Heading, EH = Estimated Heading and GPSH (Reg) = Non-linear Regression of GPS Heading).
Figure 7. Measured and estimated headings (left figures) for sinusoidal trajectories as well as heading difference (right figures) (where, GPSH = GPS Heading, MH = Measured Heading, EH = Estimated Heading and GPSH (Reg) = Non-linear Regression of GPS Heading).
Sensors 18 01390 g007
Figure 8. Outlook of a concave polygon field (top figure) and, measured and estimated heading (bottom figure) for the concave polygon field.
Figure 8. Outlook of a concave polygon field (top figure) and, measured and estimated heading (bottom figure) for the concave polygon field.
Sensors 18 01390 g008
Figure 9. Heading difference based on the reference heading for a concave polygon field.
Figure 9. Heading difference based on the reference heading for a concave polygon field.
Sensors 18 01390 g009

Share and Cite

MDPI and ACS Style

Rahman, M.M.; Ishii, K. Heading Estimation of Robot Combine Harvesters during Turning Maneuveres. Sensors 2018, 18, 1390. https://doi.org/10.3390/s18051390

AMA Style

Rahman MM, Ishii K. Heading Estimation of Robot Combine Harvesters during Turning Maneuveres. Sensors. 2018; 18(5):1390. https://doi.org/10.3390/s18051390

Chicago/Turabian Style

Rahman, Md Mostafizar, and Kazunobu Ishii. 2018. "Heading Estimation of Robot Combine Harvesters during Turning Maneuveres" Sensors 18, no. 5: 1390. https://doi.org/10.3390/s18051390

APA Style

Rahman, M. M., & Ishii, K. (2018). Heading Estimation of Robot Combine Harvesters during Turning Maneuveres. Sensors, 18(5), 1390. https://doi.org/10.3390/s18051390

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