Jingtao Tang, Zining Mao, Lufan Yang, Hang Ma
Simon Fraser University
A time-optimal deterministic spatiotemporal planner for multi-robot motion planning.
We address the Multi-Robot Motion Planning (MRMP) problem of computing collision-free trajectories for multiple robots in shared continuous environments. While existing frameworks effectively decompose MRMP into single-robot subproblems, spatiotemporal motion planning with dynamic obstacles remains challenging, particularly in cluttered or narrow-corridor settings. We propose Space-Time Graphs of Convex Sets (ST-GCS), a novel planner that systematically covers the collision-free space-time domain with convex sets instead of relying on random sampling. By extending Graphs of Convex Sets (GCS) into the time dimension, ST-GCS formulates time-optimal trajectories in a unified convex optimization that naturally accommodates velocity bounds and flexible arrival times. We also propose Exact Convex Decomposition (ECD) to “reserve” trajectories as spatiotemporal obstacles, maintaining a collision-free space-time graph of convex sets for subsequent planning. Integrated into two prioritized-planning frameworks, ST-GCS consistently achieves higher success rates and better solution quality than state-of-the-art sampling-based planners - often at orders-of-magnitude faster runtimes - underscoring its benefits for MRMP in challenging settings.
Demonstration of the proposed PBS+TGCS method for MRMP: (a)(b) The solution visualized in (a) spatial coordinates and (b) space-time coordinates with the starts (squares), goals (stars), and trajectories (dot lines) of the robots highlighted in colors; (c) Low-level planner solving a TGCS program to generate a trajectory (blue solid line) through space-time collision-free convex sets (colored polyhedrons), which avoids trajectories of a high-priority robot (red solid line); (d) Magnified view of the congested region in (c), located at the center of the map.
We present the key idea of ST-GCS by demonstrating how to augment a space decomposition (i.e., a graph of convex sets) with an unbounded time dimension for time-optimal trajectories. By enforcing constraints on time flow and velocity bounds within a unified convex optimization, we obtain piecewise-linear, time-optimal trajectories without specialized time discretization.
■ GCS[1,2] solve a generalized shortest-path problem on a graph whose vertices are convex sets, determining which sets form the path (π) and the states (x, y) within each set, which jointly optimizes a chosen objective function. An example GCS (left) with its convex sets (right) in 2D space:
The following shows the Mixed-Integer Convex Program (MICP) using GCS for motion planning in time-invariant state spaces:
■ ST-GCS augments GCS with an additional time dimension. An example ST-GCS (left) with its convex sets (right) in 3D space-time:
ST-GCS adds constraints Eqn. (8) to prevent time reversal and Eqn. (9) to bound the velocity. The arrival time of the goal position is unconstrained (Eqn. (10)) as in GCS (Eqn. (5)). The following shows the MICP using ST-GCS for time-optimal spatiotemporal motion planning:
ECD "reserves'' a given piecewise linear trajectory on a graph G of space-time convex sets, thereby producing an updated graph G' whose convex sets are collision-free with respect to the reserved trajectory. Consequently, once some robots' trajectories are planned, they can be treated as dynamic obstacles for the remaining robots by applying ECD and then running ST-GCS on G'. In general, ECD takes as input a sequence of vertex-states tuples {(vi, xi, yi)}, where each trajectory segment xi -> yi lies entirely in a single convex set of the given graph G. ECD iteratively processes on each tuple (vi, xi, yi) to decompose a convex set into multiple convex sets that are space-time collision-free with the trajectory segment xi -> yi. An example of processing a tuple is shown below:
We now introduce two prioritized planning frameworks for MRMP that use ST-GCS. In both frameworks, robots plan one at a time using ST-GCS as a low-level trajectory planner to avoid collisions with higher-priority robots’ trajectories, which are incorporated as dynamic obstacles through ECD.
PBS Search Tree. Each node branches to two child nodes with two new priority orderings when finding a conflict.
■ Random-Prioritized Planning (RP) [3] randomly explores priority orders and plans each robot sequentially according to each priority order.
■ Priority-Based Search (PBS) [4] systematically explores priority orders to resolve collisions by searching a priority tree, where each node N contains a unique priority set ≺N of ordered pairs of robots and a set of n trajectories that respect the prioritized planning scheme specified by ≺N.
■ PBS+ST-GCS only reserves conflicting trajectories deemed high-priority via ECD, controlling the ST-GCS graph growth more effectively than RP+ST-GCS, which relies on total priority orders with sequential ECD, thus expanding the graph faster as the number of robots grows.
The following shows several examples of performance comparison between SP+T-PRM [5], SP+ST-RRT* [6] and our RP+ST-GCS and PBS+ST-GCS on 2D mobile robots (SP represents sequential planning at high level). Black regions are the (static/dynamic) obstacles, and colored regions are the space collision-free convex sets. The Sum-of-Costs (SoC; the sum of time costs of all trajectories), makespan (MS; the maximum time costs of all trajectories), and planning runtime (RT) are reported as metrics. Please check our paper for more detailed results.
Map empty (4 robot switching positions;
SP+T-PRM failed after 150 secs.)
Map empty (2 robots switching positions
w/ 2 dynamic obstacles; SP+T-PRM failed after 150 secs.)
SP+ST-RRT* (SoC=3.60; MS=0.90; RT=150s)
PBS+ST-GCS (SoC=3.60; MS=0.90; RT=0.16s)
SP+ST-RRT* (SoC=5.22; MS=3.01; RT=150s)
PBS+ST-GCS (SoC=4.69; MS=2.40; RT=8.95s)
Map simple (single robot w/ 4 dynamic obstacles;
SP+T-PRM failed after 150 secs.)
Map simple (4 robots switching positions w/ 8 dynamic
obstacles; SP+T-PRM failed after 150 secs.)
ST-RRT*
(SoC=MS=5.46; RT=150s)
ST-GCS
(SoC=MS=5.38; RT=0.51s)
SP+ST-RRT* (SoC=37.3; MS=15.4; RT=150s)
PBS+ST-GCS (SoC=22.6; MS=6.71; RT=106s)
Map simple (8 robots switching positions;
SP+T-PRM failed after 150 secs.)
Map complex (4 robot + 4 dynamic obstacles;
SP+T-PRM & SP+ST-RRT* failed after 150 secs.)
SP+ST-RRT* (SoC=29.7; MS=6.88; RT=150s)
PBS+ST-GCS (SoC=26.3; MS=4.89; RT=58.3s)
RP+ST-GCS (SoC=35.6; MS=13.6; RT=51s)
PBS+ST-GCS (SoC=28.2; MS=9.5; RT=17.5s)
The following shows 4 examples of PBS+ST-GCS solving MRMP with large-sized robot teams, where all other methods failed after 150 seconds runtime.
map simple (12 robots)
SoC=26.7; MS=5.22; RT=29s
map simple (16 robots)
SoC=43.1; MS=6.22; RT=68s
map complex (12 robots)
SoC=44.5; MS=10.7; RT=32s
map complex (20 robots)
SoC=79.8; MS=7.83; RT=79s
Formation and rearrangement of "IROS" characters via our PBS+ST-GCS:
[1] Marcucci, Tobia, et al. "Shortest paths in graphs of convex sets." SIAM Journal on Optimization 34.1 (2024): 507-532.
[2] Marcucci, Tobia, et al. "Motion planning around obstacles with convex optimization." Science robotics 8.84 (2023): eadf7843.
[3] Erdmann, Michael, and Tomas Lozano-Perez. "On multiple moving objects." Algorithmica 2 (1987): 477-521.
[4] Ma, Hang, et al. "Searching with consistent prioritization for multi-agent path finding." Proceedings of the AAAI conference on artificial intelligence. Vol. 33. No. 01. 2019.
[5] Hüppi, Matthias, et al. "T-prm: Temporal probabilistic roadmap for path planning in dynamic environments." 2022 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). IEEE, 2022.
[6] Grothe, Francesco, et al. "St-rrt*: Asymptotically-optimal bidirectional motion planning through space-time." 2022 International Conference on Robotics and Automation (ICRA). IEEE, 2022.