Topological Abstraction and Planning for the Pursuit-Evasion Problem

Similar documents
Cell-Like Maps (Lecture 5)

Sheaf Theory Approach to Distributed Applications: Analysing Heterogeneous Data in Air Traffic Monitoring

Sheaf-Theoretic Stratification Learning

THE RISE, FALL AND RISE OF SIMPLICIAL COMPLEXES. Andrew Ranicki (Edinburgh) aar Warwick 2nd February, 2018

Homology cycle bases from acyclic matchings

Recent Advances and Trends in Applied Algebraic Topology

Motion Planning with Uncertainty

Observing Information: Applied Computational Topology.

What is a sheaf? Michael Robinson. SIMPLEX Program

Math 635: Algebraic Topology III, Spring 2016

Sheaf Theory: The Mathematics of Data Fusion

Robotic Motion Planning: Cell Decompositions (with some discussion on coverage and pursuer/evader)

Constraint Satisfaction Problems (CSPs)

ACLT: Algebra, Categories, Logic in Topology - Grothendieck's generalized topological spaces (toposes)

Lecture 0: Reivew of some basic material

66 III Complexes. R p (r) }.

Reach Sets and the Hamilton-Jacobi Equation

What is Computer Vision?

Topological space - Wikipedia, the free encyclopedia

Ian Mitchell. Department of Computer Science The University of British Columbia

Metrics on diagrams and persistent homology

Simplicial Global Optimization

Clustering algorithms and introduction to persistent homology

Topological Perspectives On Stratification Learning

Negations in Refinement Type Systems

Computing the Betti Numbers of Arrangements. Saugata Basu School of Mathematics & College of Computing Georgia Institute of Technology.

The Cyclic Cycle Complex of a Surface

VoI for adversarial information structures

3. Persistence. CBMS Lecture Series, Macalester College, June Vin de Silva Pomona College

Data Analysis, Persistent homology and Computational Morse-Novikov theory

Topological Classification of Data Sets without an Explicit Metric

Topological Inference via Meshing

Scan Scheduling Specification and Analysis

Western TDA Learning Seminar. June 7, 2018

CSC Discrete Math I, Spring Sets

Simplicial Complexes: Second Lecture

Hierarchical Multi-Objective Planning For Autonomous Vehicles

4.2 Simplicial Homology Groups

Lecture 18: Groupoids and spaces

Classifying Spaces and Spectral Sequences

The Borsuk-Ulam theorem- A Combinatorial Proof

Topological Data Analysis - I. Afra Zomorodian Department of Computer Science Dartmouth College

Braid groups and Curvature Talk 2: The Pieces

Chapter 11. Topological Spaces: General Properties

THE DOLD-KAN CORRESPONDENCE

Chapter 15 Introduction to Linear Programming

Constructive Homological Algebra V. Algebraic Topology background

Homework #6 (Constraint Satisfaction, Non-Deterministic Uncertainty and Adversarial Search) Out: 2/21/11 Due: 2/29/11 (at noon)

On the Topology of Finite Metric Spaces

2 Review of Set Theory

Robot Motion Planning

Capturing an Evader in a Polygonal Environment with Obstacles

15-780: MarkovDecisionProcesses

ECE276B: Planning & Learning in Robotics Lecture 5: Configuration Space

A Little Point Set Topology

Mapping and Pursuit-Evasion Strategies For a Simple Wall-Following Robot

Total No. of Questions : 18] [Total No. of Pages : 02. M.Sc. DEGREE EXAMINATION, DEC First Year COMPUTER SCIENCE.

Computational Topology in Neuroscience

Final Exam Practice Fall Semester, 2012

Directed algebraic topology and applications

Topology - I. Michael Shulman WOMP 2004

Sheaf. A sheaf is a functor from from the category of simplicial complexes. Example: σ F(σ), a vector space ρ σ L : F(ρ) F(σ), a linear map

Solving Pursuit-Evasion Problems with Graph-Clear: an Overview

4. Simplicial Complexes and Simplicial Homology

Random Simplicial Complexes

HOMOLOGICAL SENSOR NETWORKS

Relational Database: The Relational Data Model; Operations on Database Relations

Random Simplicial Complexes

A New Combinatorial Design of Coded Distributed Computing

A deterministic action is a partial function from states to states. It is partial because not every action can be carried out in every state

AUTOMORPHISMS OF THE GRAPH OF FREE SPLITTINGS

Homological Algebra and Data

1.2 Numerical Solutions of Flow Problems

Reach Sets and the Hamilton-Jacobi Equation

Homology and Persistent Homology Bootcamp Notes 2017

Pacific Journal of Mathematics

Advanced Robotics Path Planning & Navigation

Numerical schemes for Hamilton-Jacobi equations, control problems and games

STANLEY S SIMPLICIAL POSET CONJECTURE, AFTER M. MASUDA

Using Hybrid-System Verification Tools in the Design of Simplex-Based Systems. Scott D. Stoller

A GRAPH FROM THE VIEWPOINT OF ALGEBRAIC TOPOLOGY

TOPOLOGICAL DATA ANALYSIS

Combinatorial Algebraic Topology and applications to Distributed Computing

Topological Inference via Meshing

Discrete Mathematics Lecture 4. Harper Langston New York University

What is high-dimensional combinatorics?

Statistical topological data analysis using. Persistence landscapes

Detection and Mitigation of Cyber-Attacks using Game Theory

An Evolution of Mathematical Tools

Simplicial Objects and Homotopy Groups

Crossing Families. Abstract

DISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane CMU /858B Fall 2017

4 Basis, Subbasis, Subspace

Combinatorics I (Lecture 36)

Simplicial Complexes of Networks and Their Statistical Properties

Topological Invariance under Line Graph Transformations

Topology and Topological Spaces

arxiv: v1 [math.at] 8 Jan 2015

Introduction to Computational Manifolds and Applications

Lecture 3: Motion Planning 2

Transcription:

Topological Abstraction and Planning for the Pursuit-Evasion Problem Alberto Speranzon*, Siddharth Srivastava (UTRC) Robert Ghrist, Vidit Nanda (UPenn) IMA 1/1/2015 *Now at Honeywell Aerospace Advanced Technology This material is based upon work supported by Defense Advanced Research Project Agency under contract No. HR0011-15-C-0043 Sponsored by the Defense Advanced Research Project Agency Defense Sciences Office (DSO) Program: DARPA Seedling Program Purchase request No.: HR0011513129, Program code: DMSSD Issued by DARPA/CMO under contract No. HR0011-15-C-0043 The views and conclusions contained in this presentation are those of the authors and should not be interpreted as representing the official policies, either expressly or implied, of the Defense Advanced Research Projects Agency or the U.S. Government. 1

Complex Planning Problems: Hierarchical Control Consider a complex planning problem: Mission level problem what the whole system should do Motion level problem what the single sub-systems do Control level problem generate kinematically/dynamically feasible motions It is generally convenient to approach the design considering a hierarchical control structure We consider pursuit-evasion problem (P/E) as a case study HIGH-LEVEL Refinement Planning over Mission Planner Sheaves Abstraction * e 2 * * p 2 p 3 LOW-LEVEL Planning over Motion Planner the Nerve Cplx * p 1 * e 1 * * Refinement Controller Abstraction 2

Complex Planning Problems: Hierarchical Control Consider a complex planning problem: Mission level problem what the whole system should do Motion level problem what the single sub-systems do Control level problem generate kinematically/dynamically feasible motions It is generally convenient to approach the design considering a hierarchical control structure We consider pursuit-evasion problem (P/E) as a case study HIGH-LEVEL Refinement Planning over Mission Planner Sheaves Abstraction * e 2 * * p 2 p 3 LOW-LEVEL Planning over Motion Planner the Nerve Cplx * p 1 * e 1 * * Refinement Controller Abstraction 3

Assumptions of P/E Problem in this Talk Multiple evaders in the domain Unknown number of evaders in the domain Different sensor capabilities of the pursuers Sensor footprint Multiple pursuers to detect one evader Limited range of pursuers sensors p 1 p e 2 2 p 3 e 1 Communication constraints: not all-to-all Deterministic models: position of the pursuers and sensing Known environment but not restricted to 2D, with obstacles 4

Contents Sequential decision making Pursuit/Evasion and topological abstractions Decision making over sheaves Generalization of the P/E with sheaves Nerve complex and downward refinement Conclusions 5

Sequential Decision Making We consider the special case: deterministic, fully observable, i.e. an MDP Initial state Given transition system T = (S, A, T, S goal, S init ): S: Set of states, A: set of action symbols T: S x A S, deterministic transition function Possible actions State Desired set of states S goal Compute: strategy for reaching a desired state Forms of solutions: Finite state machine Policy: State Action Goal state 6

P/E Abstractions: State/Action Different level of abstractions have been considered in the problem of P/E: Continuous versions: generally posed as a differential game requiring the solution of Hamilton-Jacobi-Bellman-Isaac equation. Solutions available only for simple models of P s and E s and environments Discrete versions: generation of a discrete representation, typically as a graph grid or visibility graph. Planning is then computed as path along such graph. Ad-hoc abstraction using a grid Several assumptions: Configuration space is the same as the map space Sensor footprint matches the graph/grid resolution Downward refineability holds Result in post-facto engineering for fixing the problem with on-the-fly replanning 7

P/E New Perspective on Abstraction We consider a topological abstraction of the problem Pursuer space P Evader space E Ad-hoc abstraction using a grid Topological Abstraction Pursuer space P 2 connected components with 3 holes Evader Space E 4 connected components 8

P/E New Perspective on Abstraction We consider a topological abstraction of the problem Pursuer space P Evader space E Ad-hoc abstraction using a grid Topological Abstraction Pursuer space P 2 connected components with 3 holes Evader Space E 4 connected components The topological abstraction is minimal as it encodes the critical information about the problem 4 connected areas where evaders can be hiding. However, how can we make it useful to generate actionable plans? We need a way to encode more (metric) information about the overall problem Sheaf as a way to capture information necessary for planning 9

Sheaf, Restriction Maps and Global Sections A sheaf is a mathematical object that stores locally-defined data over a space and allows to infer about global properties of the space from such data EXAMPLE: 10

Sheaf, Restriction Maps and Global Sections A sheaf is a mathematical object that stores locally-defined data over a space and allows to infer about global properties of the space from such data EXAMPLE: 11

Sheaf, Restriction Maps and Global Sections A sheaf is a mathematical object that stores locally-defined data over a space and allows to infer about global properties of the space from such data EXAMPLE: Stalks Base Space Restrictions maps, are maps between two cells σ and τ where σ is a face of τ EXAMPLE OF RESTRICTIONS MAPS: 1 2 A 1 12 : linear map between σ = {1} and τ = 12 A 2 12 : linear map between σ = {2} and τ = {12} 12

Simple (Familiar) Example Consider discrete points (0-simplexes) at certain time instances Assign vector spaces to each 0-simplex And to each 1-simplex as well I A n 1 R k A n R k I R k I R k R k A n+1 n-1 n n+1 In this particular example, the sheaf co-homology represents the evolution of a linear time varying (or invariant if all the restriction maps are constant) system: u n+1 = A n u n 13

Sheaf Theory: Rich Set of Tools There are a number of compelling tools suited to reason about abstractions based on sheaves GLOBAL SECTION FUNCTOR, H 0 ( ) Collates all solutions to the constraints imposed by the sheaf stalks and restriction maps; gives an algebraic form to the solution set. In P/E, H 0 ( ) classifies connectivity of evasion paths. COHOMOLOGY, H ( ) H ( ) characterizes constraint satisfaction as a function of both the domain topology and the algebra of the constraints. In P/E, this gives necessary and/or sufficient criteria for evasion. PUSHFORWARDS/PULLBACKS These operations transform sheaves from one base space to another In P/E, these allow for dynamic updates to domain/parameters. HOM AND TENSOR PRODUCTS HOM classifies relationships between sheaves; convolves sheaf data. In P/E, HOM gives a collation of constraints; acts as sensor fusion. 14

Sheaf Abstraction of P/E Problems: Evader Sheaf 15

Sheaf Abstraction of P/E Problems: Evader Sheaf Action: (up, right) (Merging) 16

Sheaf Abstraction of P/E Problems: Evader Sheaf Action: (up, right) (Merging) Action: (up, right) (Splitting + Death) 17

Sheaf Abstraction of P/E Problems: Evader Sheaf Action: (up, right) (Merging) Action: (up, right) (Splitting + Death) Action: (down, up) (Born) 18

Sheaf Abstraction of P/E Problems: Evader Sheaf Action: (up, right) (Merging) Action: (up, right) (Splitting + Death) Action: (down, up) Compress evasion region to connected components Use H 0 ( ; R) to have sheaf take values in real vector spaces Use (persistent) co/homology to track evader components Note that at this level of abstraction we are considering an abstract state where each state encodes the number of connected components of E Here we depict one sheaf for one set of actions, but we have multiple sheaves for each possible set of actions (Born) 19

Sheaf Abstraction of P/E Problems: Evader Sheaf Action: (up, right) Action: (up, right) Action: (down, up) In this simple setting, the sheaf captures the evolution of the connected components of the evader space Restriction maps in this context (can be made) are maps between vector spaces 20

Sheaf Abstraction of P/E Problems: Evader Sheaf A 1 A 2 A m Sheaf 1 Sheaf 2 Sheaf k For each path in this search space we can associate a sheaf corresponding to a sequence of actions π = {A i1, A i2, A i3,, A i n } Evader capture in this context corresponds to the fact that dim H 0 ; R = 0 or equivalently that the (evader) sheaf has no global section 21

Simpler Example: Global Section in P/E up down right Note that in this case there is no continuous path from left to right This indicates that there are no global sections in the sheaf and equivalently that the sequence of actions corresponds to a capture Thus we have that in the P/E problem, the evader is captured if there exists a set of actions such that the corresponding sheaf has no global sections 22

Simpler Example: Global Section in P/E up down right Note that in this case there is no continuous path from left to right This indicates that there are no global sections in the sheaf and equivalently that the sequence of actions corresponds to a capture Thus we have that in the P/E problem, the evader is captured if there exists a set of actions such that the corresponding sheaf has no global sections Sheaves provide a very general framework to approach this problem as we can encode arbitrary data structures on the sheaves thus capturing richer settings 23

So What Have We Done? AGENT MODEL ABSTRACTION AGENT MODEL ABSTRACTION OVER A DOMAIN NERVE CONSTRUCTION EVADER SHEAF AND PLANNING REDUCE TRANSITION SYSTEM THROUGH STATE ABSTRACTION TRANSITION SYSTEM 24 45

Extensions Using Sheaf Theoretic Methods Domain Size & Topology Sensing & evasion sheaves Evader Capture Criteria E 1 = P 1 P 2 + P 1 P 3 + P 4 E 2 = P 2 P 3 + P 2 P 4 E 3 = P 1 P 2 P 3 + P 1 P 3 P 4 Augment sheaf stalks; build escape and capture complex Communication Network P 4 P 1 Factor into base space as communication complex P 3 P 2 Temporal Dynamics VS Sheaves of semigroups; duality 25

Capture Constraints: Propositional Logic E 1 = P 1 P 2 + P 1 P 3 + P 4 E 2 = P 2 P 3 + P 2 P 4 E 3 = P 1 P 2 P 3 + P 1 P 3 P 4 Each evader is assigned its own capture criterion in the form of a Boolean expression involving (some of) the pursuers In this example, the first evader can be caught only if one of the three following conditions it satisfied it must be seen simultaneously by pursuers P 1 and P 2, or seen simultaneously by pursuers P 1 and P 3, or seen by pursuer P 4 alone On the other hand, the third evader needs to be simultaneously seen by pursuers P 1 and P 3 along with one of P 2 or P 4. In this manner, we can encode vastly asymmetric models of pursuer and evader capabilities. MULTIPLICATION AND : ADDITION OR 26

Evasion and Capture Sheaf Idea: Encode Boolean capture criteria among pursuers/evaders as a complex. Use this complex as the base of an Evasion Sheaf, to enforce AND operations. Define a secondary Capture Sheaf, to enforce OR operations. Given: Sensed (coverage) regions C i, t in R n as a function of time t (action sequence) and pursuer P i positions; Boolean statements for capture criteria among pursuers/evaders. Base space: product of action sequence (R) with Escape Complex DISCUSSION FOLLOWS THE SIMPLEST TIME-INDEPENDENT CASE NOVEL CONSTRUCTION The Escape Complex is the subcomplex of the full simplex on the pursuer set {P i } defined as the maximal subcomplex of the simplex not intersecting the (open) cells defined by the Boolean pursuer relations. 27

Escape Complex The Escape Complex is the subcomplex of the full simplex on the pursuer set {Pi} that consists of the largest subcomplex of the simplex not intersecting the (open) cells defined by the Boolean pursuer relations. EXAMPLE: If the capture criterion is E = P 1 P 2 + P 1 P 3 = P 1 (P 2 + P 3 ) (capture requires participation of P 1 and either P 2 or P 3 ) P 1 E = P 1 + (P 2 P 3 ) P 1 ESCAPE COMPLEX FULL SIMPLEX P 3 P 2 P 3 P 2 28

Capture Criteria: AND EXAMPLE: If the capture criterion is E = P 1 P 2 + P 1 P 3 = P 1 (P 2 + P 3 ) Let the domain be a simple interval with four nodes (evader regions) At a fixed time let the pursuer coverage regions be: P 1 P 2 P 3 Build the evasion sheaf as follows: R 1 Stalks = vector space whose dimension over vertices records size of uncovered region P 1 R 2 R 4 R 2 Sheaf maps = inclusions identifying common regions P 3 P 2 29

Capture Criteria: AND Cohomology over the escape complex tracks the AND operations in escape regions In this simple example, one easily computes that dim H 0 = 2. One global section comes from what lies in the complement of P 1 s coverage P 1 c P 1 R 1 R 2 R 4 R 2 P 2 c P 3 c P 3 P 2 The other global section is supported over the P 2 -P 3 edge and is generated by the intersection of complements of P 2 and P 3 coverage regions This yields the AND operations in the factored Boolean expression E = P 1 P 2 + P 1 P 3 = P 1 (P 2 + P 3 ) 30

Capture Criteria: OR Express the OR as a cohomology over another complex, the Capture Complex In this simple example, we need to intersect capture by P 1 with capture by (P 2 OR P 3 ) Consider the simplicial complex generated by these two conditions in the AND. P 1 c P 1 R 3 R 4 P 2 c P 3 c P 2 + P 3 R 3 Repeat the construction: take complements of the regions implicit in H 0 of the evasion sheaf; use these as stalks of the Boolean capture sheaf, then compute H 0. In this example, dim H 0 (CAP; H 0 (ESC) ) = 2 the two locations of capture. In general, this yields the OR operations in the Boolean expression E = P 1 P 2 + P 1 P 3 = P 1 (P 2 + P 3 ) 31

Communication Constraints We can accommodate complex constraints on how the pursuers are permitted to communicate with one another P 4 P 1 These constraints are represented by a simplicial complex whose vertices index the pursuers. An edge implies communication that is pairwise higher collaboration comes from higher simplices. P 3 P 2 All planning and coordination become subordinate to the structure of this communication complex: in the maximal case, one has full communication and a full simplex; in the worst case, there are only vertices and no communication. We assume a model where communication leads to symmetric sharing of information about coverage; this corresponds to taking an OR operation with coverage. 32

Communication Complex Idea: Encode communication constraints among pursuers as a simplicial complex Given: Sensed (coverage) region C i in R n as a function of pursuer P i ; Communication constraints as complex on pursuer set. Base space: product of time axis (action sequence) R with Communication Complex Stalks: begin with R n C i, i.e. the non-sensed regions where evaders can hide; convert these to vector spaces whose dimensions over vertices give size of uncovered region P 4 P 1 P 3 P 2 P 3 and P 4 can communicate P1, P2 and P3 can communicate with each other 33

Communication Complex: Example EXAMPLE: Consider a collection of three pursuers viewing a circular domain. COMMUNICATION COMPLEX R 4 FULL SIMPLEX P 1 P 1 P 3 R 6 P 1 R 6 R 2 R 3 P 3 P 2 P 2 P 3 R 6 P 2 TOTAL COMMUNICATION COVERAGE REGIONS COMMUNICATION SHEAF (stalk dimension ~ evasion set) In this example, the cohomology H 0 = 0, since the lower P 2 -P 3 edge has trivial kernel. In other words, when everyone communicates, the evasion region vanishes no escape. 34

Communication Complex: Example EXAMPLE: Consider a collection of three pursuers viewing a circular domain. COMMUNICATION COMPLEX R 4 P 1 P 1 P 3 R 6 P 1 P 2 comms broken R 3 R 3 P 3 P 2 P 2 P 3 P 2 TOTAL COMMUNICATION COVERAGE REGIONS COMMUNICATION SHEAF (stalk dimension ~ evasion set) In this example, the cohomology H 0 = R 4 when P 2 cannot communicate, meaning that an evader can hide from P 2 (in three locations) or from P 1 and P 3 (in one). 35

Temporal Constraints [Ghrist & Krishnan 15] All of our data are allowed to vary with time. The directionality of time makes it difficult to classify evasion paths as topological constructs. sensed region, time-varying true evasion path (global section of evader sheaf) Action sequence axis false evasion path (violates causality) Action sequence axis 36

Complex Planning Problems: Hierarchical Control Consider a complex planning problem: Mission level problem what the whole system should do Motion level problem what the single sub-systems do Control level problem generate kinematically/dynamically feasible motions It is generally convenient to approach the design considering a hierarchical control structure We consider pursuit-evasion problem (P/E) as a case study HIGH-LEVEL Refinement Planning over Mission Planner Sheaves Abstraction * e 2 * * p 2 p 3 LOW-LEVEL Planning over Motion Planner the Nerve Cplx * p 1 * e 1 * * Refinement Controller Abstraction 37

Nerve Complex Consider a collection U = U α of (compact) subsets of a topological space The nerve N U of U is built as follows: The k-simplices of N(U) correspond to nonempty intersection of k + 1 distinct element of N(U) EXAMPLE: 3 U 1 U 2 U 3 U 6 1 5 6 U 4 U 5 2 4 38

Abstract Model Agent Motion Sensing Agent Abstract Model Given: Sensing model of agents Motion model of agents We can build low-level model that capture various details: θ = π θ = π θ y x 39

Abstract Model Action Definition for 2D Environments Time is discretized Actions are determined from a minimal cover of the 1-step sensing envelop a 8 a 1 a 2 a 7 a 3 a 2, a 3 a 4, a 5 a 1 a 6 a 6 a 5 a 4 A = {a 1, a 2,, a M } 40

State Abstraction via Nerve s 0 a 1 a 2 a 3 a M s 1 s 2 s 3 s N a 1 a 2 a M Generate the transitions system T from action A and state S for a given environment Nerve complex is used to determine the number of connected components corresponding to various configurations 41

State Abstraction via Nerve s 0 a 1 a 2 a 3 a M s 1 s 2 s 3 s N a 1 a 2 a M Generate the transition system T from action A and state S for a specific environment For example: in the figure any of the three positions of the green pursuer or the blue pursuer will lead to the same number of connected components 42

State Abstraction via Nerve s 0 a 1 a 2 a 3 a M s 1 s 2 s 3 s N a 1 a 2 a M Generate the transition system T from action A and state S for a specific environment A sub-tree/sub-branch of the full T can be collapsed into a single state and action based on the number of connected components obtained using the nerve construction: A 1 = {a i1, a i2, a i3, } A 1 43

Advantages of Nerve Construction Observations: Transition system T needs to be constructed only once Transition system T can be very large however we never use T to compute a plan Transition system T is much smaller and suitable for planning Nerve construction generalizes the grid construction that is thus encompassed with this model Downward refinability easier to achieve within this more general model 44

Summary AGENT MODEL ABSTRACTION AGENT MODEL ABSTRACTION OVER A DOMAIN NERVE CONSTRUCTION EVADER SHEAF AND PLANNING REDUCE TRANSITION SYSTEM THROUGH STATE ABSTRACTION TRANSITION SYSTEM 45

Conclusions ``Semantic abstraction of the pursuit-evasion problem using connected components of the evader space Connected states to actions through using a sheaf Topological properties (global sections) of the sheaf lead to evader capture conditions Nerve based construction to abstract low-level control laws Left To Do (Everything! ) Global section is a too coarse notion: how can we select branches that are most promising? How do we capture uncertainty? Computational aspects. 46