Incremental Sampling-Based Algorithms and Stochastic Optimal Control on Random Graphs

This project develops incremental sampling-based motion planning algorithms and stochastic optimal control methods on random graphs to address computational bottlenecks in high-dimensional autonomous systems

Motion planning and trajectory generation problems are at the core of all autonomous and robotic systems. Although – at least in theory – such problems can be solved using optimal control or dynamic programming, the computational complexity for realizing these solutions is still prohibitive for many real-life problems, especially for high-dimensional systems. Recently, randomized algorithms have been proposed, which promise to remedy the “curse of dimensionality” associated with more traditional, deterministic methods.

In recent years, sampling-based motion planning algorithms (SBMPs) have become popular due to their ability to handle higher dimensions and kino-dynamic constraints. Incremental sampling based algorithms, in particular, Rapidly-exploring Random Trees (RRT), RRT* avoid apriori discretization of the search space and build a connectivity graph online by generating random samples from the search space. However, these randomized methods are hindered by three computational bottlenecks, namely:

  • The need for an intelligent sampling strategy to avoid redundant exploration
  • Avoiding an expensive collision checking module
  • Efficient ways of implementing the so-called “local steering” problem connecting any two vertices in the connectivity graph.

In this research, we propose ideas from classical heuristics, dynamic programming and machine learning to address these issues.

The starting point is the RRT# algorithm, which rewires the graph globally using value iteration, so that every vertex is optimally connected to the root at the end of each iteration rather than only locally repaired. The planner was rewritten in C/C++ and contributed to the Open Motion Planning Library, where it can be compared directly against the other planners in the literature. We then went further and optimized the locations of the samples as well as the wiring between them. The resulting Deformable RRT (DRRT) treats each new sample as a local perturbation: the branch from the root to the new sample is optimized before the change is propagated through the tree, which avoids a global gradient descent at every iteration. To measure how these algorithms degrade with dimension we introduced the hypercube diagonal experiment, and on it DRRT converges orders of magnitude faster than comparable sampling-based planners for any problem of dimension greater than three.

Three stages of an RRT# search in a planar environment with rectangular obstacles: as the iterations proceed the tree becomes progressively sparser away from the solution while the path from the start to the goal is refined RRT# on a planar problem. Global rewiring by value iteration keeps every vertex optimally connected, so the search concentrates on the region that can still improve the solution.

Because these planners work in configuration space rather than the plane, they apply directly to manipulators and other high-dimensional systems. The videos below show RRT# planning for anthropomorphic robots with six and twelve degrees of freedom, and DRRT planning for a manipulator.

Anthropomorphic robot motion planning with six degrees of freedom.

The same problem with twelve degrees of freedom.

Manipulator planning with the Deformable RRT algorithm.

Six rendered views of a dual-arm manipulator over a table of cups, arranged as three start configurations on the top row and the three corresponding goal configurations below Start and goal configurations for a dual-arm planning problem, where the two arms must reach their targets without colliding with each other or with the objects on the table.

The second thread attacks the exploration problem: where to draw the next sample. A good strategy has to find some initial solution quickly and then concentrate on the region that can still improve it. The non-parametric informed sampler does this by keeping a subset of well-dispersed kernel vertices and generating new samples in their vicinity, scoring each candidate by a combination of its heuristic value, the local density of the graph and a count-based exploration term, and steering away from directions where collisions have already been recorded. For planning on general cost-maps, where the cost varies continuously rather than being a binary obstacle test, the relevant region sampler uses the cost-to-come values already stored at the graph vertices, together with an admissible heuristic, to bound the cost of any solution passing through a given point, and samples only where that bound beats the current solution. Later, using ideas from reachability analysis, we defined a time-informed set containing every trajectory better than the incumbent, proved that exploring outside it is redundant, and built a library of reachable sets that makes sampling inside it cheap.

Schematic of the non-parametric informed sampler: kernel vertices along the current path inside a dashed L2 informed set, with candidate sampling directions at one vertex and previously recorded obstacle directions excluded Non-parametric informed sampling. New samples are drawn near well-dispersed kernel vertices, inside the informed set, and away from directions in which collisions have already been seen.

Two planning problems on continuous cost-maps rendered as greyscale images, a terrain map on the left and a smooth potential field on the right, each with the search tree in green and the resulting path in magenta Relevant region sampling on a terrain cost-map (left) and a potential cost-map (right). The tree spreads through the low-cost regions instead of the whole space.

The third bottleneck is local steering, the two-point boundary value problem that must be solved every time a node is expanded. It has a closed-form solution only for simple linear systems, and solving it numerically becomes the dominant cost as the dynamics grow more complex. We therefore learned it instead: optimal trajectories computed offline become training data for a goal-conditioned state-feedback neural network, which predicts the connecting control online and serves as the steering routine inside a kinodynamic RRT*. Separately, the policy-iteration variant of RRT# was parallelized on an NVIDIA GPU. Profiling showed that invoking the CUDA kernel every iteration is dominated by host-to-device transfer overhead, so the kernel is instead called at strategic intervals, solving a larger dynamic programming problem less often — which runs significantly faster than the serial planner while producing solutions of the same quality.

Two views of a three-dimensional voxel reconstruction of an outdoor scene in false colour, each showing a planned trajectory in dark red threading between the obstacles Kinodynamic planning through a three-dimensional reconstructed environment, where the connections between samples are produced by the learned steering controller.

The algorithms were tested both in simulation and on hardware. A fleet of ten 1/10-scale RACECAR ground robots was assembled, each carrying its own sensing and an NVIDIA Jetson processor, with a ROS navigation stack used to run the planners developed here. The closed-loop variant, CL-RRT#, was also implemented on a Mars rover testbed at NASA’s Jet Propulsion Laboratory as part of their next-generation autonomous navigation work.

A 1/10-scale RACECAR ground robot: a four-wheeled chassis carrying a lidar, a depth camera and an on-board computer

Non-parametric informed sampling concentrating the search after an initial solution is found.

Autonomous navigation on the RACECAR platform in simulation.

CL-RRT# running within the next-generation AutoNav pipeline on a rover testbed at NASA JPL.

Sponsors

This project is sponsored by NSF.

Selected Publications and Presentations

1.

Arslan, O., and Tsiotras, P., “Incremental Sampling-based Motion Planners using Policy Iteration Methods,” 55th IEEE Conference on Decision and Control, Las Vegas, NV, Dec. 12-14, 2016, pp. 5004-5009.

2.

Arslan, O., Berntorp, K., and Tsiotras, P.,“Sampling-based Algorithms for Optimal Motion Planning Using Closed-loop Prediction,” International Conference on Robotics and Automation, Singapore, May. 12-14, 2017, pp. 4991-4996.

3.

Hauer, F., and Tsiotras, P., “Deformable Rapidly-Exploring Random Trees,” Robotics: Science and Systems, Cambridge, MA, July 12-16, 2017.

4.

Joshi, S. and Tsiotras, P., “Non-Parametric Informed Exploration for Sampling-Based Motion Planning,” IEEE International Conference on Robotics and Automation, Montreal, Canada, May 20-24, 2019, pp. 5915-5921.

5.

Joshi, S., and Tsiotras, P., “Relevant Region Exploration on General Cost-Maps for Sampling-Based Motion Planning,” IEEE/RSJ International Conference on Intelligent Robots and Systems, Las Vegas, NV, Oct. 25-29, 2020, pp. 6689-6695.

6.

Lawson, R. C, Wills, L., and Tsiotras, P., “GPU Parallelization of Policy Iteration RRT#,” IEEE/RSJ International Conference on Intelligent Robots and Systems, Las Vegas, NV, Oct. 25-29, 2020.

7.

Joshi, S., Hutchinson, S., and Tsiotras, P., “TIE: Time-Informed Exploration for Robot Motion Planning,” IEEE International Conference on Robotics and Automation, Xi'an, China, May 30-June 5, 2021.

8.

Joshi, S., Hutchinson, S., and Tsiotras, P., “TIE: Time-Informed Exploration for Robot Motion Planning,” IEEE Robotics and Automation Letters, Vol. 6, No. 2, pp. 3585-3591, 2021.

9.

Zheng, D., and Tsiotras, P., “Near-Optimal Finite-Time Feedback Controller Synthesis Using Supervised and Unsupervised Learning,” AIAA SciTech Forum, Guidance, Navigation, and Control, Nashville, TN, Jan. 11-15, 2021.

10.

Zheng, D., and Tsiotras, P., “Sampling-based Kinodynamic Motion Planning Using a Neural Network Controller,” AIAA SciTech Forum, Guidance, Navigation, and Control, Nashville, TN, Jan. 11-15, 2021.