**SYSTEM AND METHOD FOR CONTROLLING MOTION OF VEHICLE**

WANG, Yebin (Inc. 201 Broadway, Cambridg, Massachusetts ., 02139, US)

DAI, Jin (Inc. 201 Broadway, Cambridg, Massachusetts ., 02139, US)

*;*

**G05D1/02***;*

**A01D34/00***;*

**A47L11/40***;*

**B60W10/20***;*

**B60W30/06***;*

**B60W50/00***;*

**B62D15/02***A47L9/28*;

*B60W50/06*

US20160313735A1 | 2016-10-27 | |||

US8428820B2 | 2013-04-23 |

BERNHARD MULLER ET AL: "Continuous Curvature Trajectory Design and Feedforward Control for Parking a Car", IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY., vol. 15, no. 3, 1 May 2007 (2007-05-01), US, pages 541 - 553, XP055466566, ISSN: 1063-6536, DOI: 10.1109/TCST.2006.890289

STEPHAN SCHEUREN: "Prozessoptimierte Planung für kooperative mobile Roboter", 1 March 2014 (2014-03-01), pages Frontpg., 18, 72, XP002779995, Retrieved from the Internet

JIN DAI ET AL: "On existence conditions of a class of continuous curvature paths", 2017 36TH CHINESE CONTROL CONFERENCE (CCC), 1 July 2017 (2017-07-01), pages 6773 - 6780, XP055466558, ISBN: 978-988-1563-93-4, DOI: 10.23919/ChiCC.2017.8028426

[CLAIMS] [Claim 1] A system for controlling a motion of a vehicle from an initial state to a target state, comprising: a path planner to determine a discontinuous curvature path connecting the initial state with the target state by a sequential composition of driving patterns, wherein the discontinuous curvature path is collision-free within a tolerance envelope centered on the discontinuous curvature path; a path transformer to locate and replace at least one treatable primitive in the discontinuous curvature path with a corresponding continuous curvature segment to form a modified path remaining within the tolerance envelope, wherein each treatable primitive is a predetermined pattern of elementary paths; and a controller to control the motion of the vehicle according to the modified path. [Claim 2] The system of claim 1, wherein the path planner is a Reeds-Shepp (RS) path planner forming the discontinuous curvature path using a combination of a straight-line segment and a circular arc tangentially connected at discontinuous curvature junctions, further comprising: a converter to convert the sequential composition of driving patterns into a sequence of treatable primitives. [Claim 3] The system of claim 2, wherein the converter converts each RS driving pattern independently from each other., [Claim 4] The system of claim 2, wherein the converter inserts a new elementary path in the sequential composition of RS driving patterns or splits at least one existing elementary path in the sequential composition of RS driving patterns to modify the sequential composition and partition the modified sequential composition of RS driving patterns into a sequence of treatable primitives. [Claim 5] The system of claim 2, wherein the converter iteratively tests each tuple of elementary paths forming the discontinuous curvature path for being a treatable primitive and inserts a new straight-line segment staying inside the envelope of the discontinuous curvature path into a tuple not representing the treatable primitive. [Claim 6] The system of claim 5, wherein the converter inserts the new straight-line segment in the front of the tuple or after the second elementary path of the tuple. [Claim 7] The system of claim 1 , wherein the path planner imposes a constraint on the discontinuous curvature path requiring at least one of an initial and a final segment of the discontinuous curvature path being a straight-line segment, such that the discontinuous curvature path is a sequence of the treatable primitives. [Claim 8] The system of claim 1, wherein the treatable primitive includes a sequence of a straight-line segment connected to a first arc segment connected to a second arc segment, and wherein the path transformer replace the treatable primitive with a pair of clothoid turns. [Claim 9] The system of claim 1 , wherein the treatable primitive includes a sequence of a first arc segment connected to a second arc segment connected to a straight-line segment, and wherein the path transformer replace the treatable primitive with a pair of clothoid turns. [Claim 10] The system of claim 1 , wherein the path transformer replaces a sequence of an even number of circular arc segments connected to the straight-line segment with a sequence of the even number of clothoid turns. [Claim 11] The system of claim 10, further comprising: a converter to split a circular arc of the discontinuous curvature path to convert a sequence of an odd number of circular arc segments into the sequence of the even number of circular arc segments. [Claim 12] The system of claim 1 , wherein the path transformer determines parameters of the continuous curvature segment to satisfy a μ-tangency constraint and remain in the envelope of the discontinuous curvature path. [Claim 13] A method for controlling a motion of a vehicle from an initial state to a target state, wherein the method uses a processor coupled with stored instructions implementing the method, wherein the instructions, when executed by the processor carry out at least some steps of the method, comprising; determining a discontinuous curvature path connecting the initial state with the target state by a sequential composition of driving patterns, wherein the discontinuous curvature path is collision-free within a tolerance envelope centered on the discontinuous curvature path; locating and replacing at least one treatable primitive in the discontinuous curvature path with a corresponding continuous curvature segment to form a modified path remaining within the tolerance envelope, wherein each treatable primitive is a predetermined pattern of elementary paths; and controlling the motion of the vehicle according to the modified path. [Claim 14] The method of claim 13, wherein the path planner is a Reeds-Shepp (RS) path planner forming the discontinuous curvature path using a combination of a straight-line segment and a circular arc tangentially connected at discontinuous curvature junctions, further comprising: converting the sequential composition of driving patterns into a sequence of treatable primitives. [Claim 15] The method of claim 13, wherein the determining imposes a constraint on the discontinuous curvature path requiring at least one of an initial and a final segment of the discontinuous curvature path being a straight-line segment, such that the discontinuous curvature path is a sequence of the treatable primitives. [Claim 16] A non-transitory computer readable storage medium embodied thereon a program executable by a processor for performing a method, the method comprising: determining a discontinuous curvature path connecting the initial state with the target state by a sequential composition of driving patterns, wherein the — -discontinuous curvature patiVis collision-free within-a tolerance envelope centered on the discontinuous curvature path; locating and replacing at least one treatable primitive in the discontinuous curvature path with a corresponding continuous curvature segment to form a modified path remaining within the tolerance envelope, wherein each treatable primitive is a predetermined pattern of elementary paths; and controlling the motion of the vehicle according to the modified path. |

[Title of Invention]

SYSTEM AND METHOD FOR CONTROLLING MOTION OF VEHICLE

[Technical Field]

[0001]

This invention relates generally to path planning of motion of the vehicles, and more particularly to path planning systems and methods for automatically controlling a motion of a vehicle between different states.

[Background Art]

[0002]

Several control systems, employed by either autonomous vehicles or vehicles executing in autonomous-driving mode, predict future, safe motions, or paths, of the vehicle, both in order to avoid obstacles, such as other vehicles or pedestrians, but also to optimize some criteria associated to the operation of the vehicle. The vehicle is controlled along the path that is a part of the object trajectory or, more precisely, the projection of the trajectory onto the

two-dimensional plane of movement. To say it more figuratively, the path is the line connecting all geometrical coordinates of the moving vehicle.

[0003]

A fundamental objective of path planning is to construct a trajectory connecting two states, e.g., an initial state and a target state, of a moving object, such as autonomous or semi-autonomous vehicles executing in

autonomous-driving mode, and to control the motion of the vehicles along the trajectory. A "state" as used herein is meant to be a vector including the geometrical coordinates of the vehicle in the plane and an orientation of the vehicle at these geometrical coordinates. The target state can either be a fixed location, a moving location, a velocity vector, a region, or a combination thereof.

The surroundings, such as road edges, pedestrians, and other vehicles, are sensed

l by the sensors of the vehicle and/or are at least partially known by a priori given information.

[0004]

A path connecting an initial state and a target state is normally made up of several elementary paths (or EPT for short). In some situations, the path planning solely works with elementary paths of the same or different kinds. To that end, the path of a moving body is normally made up of a plurality of elementary paths that share the same junction points (abbreviated as JP). Thus, the junction points are special points of the path because each of them belongs to two adjoining EPTs.

[0005]

One special characteristic demanded of the paths is that states do not change erratically along the entire course or at the junction points. This is to say that the transition from one elementary path to an adjacent elementary path must be ensured to be continuous at the junction points of the elementary paths as well as both its first derivative (tangent) and its second derivative (curvature). The path is continuous at junction points when the transition from one elementary path to the adjacent elementary path can be performed by changing a steering and velocity of the vehicle in a linear or continuous manner.

[0006]

For example, pioneering work by Dubins and Reeds and Shepp

investigated patterns of a shortest path between two states of the vehicle. The shortest path, termed as Reeds- Shepp 's (RS) path, is a sequential composition of line segments and/or tangential circular arcs of a minimum/fixed turning radius. The RS path can be computed in a very efficient way, but the curvature along the RS path involving multiple segments is discontinuous. Such discontinuity of the curvature is undesirable in practice, since a vehicle has to stop and perform stationary steering, leading to unnecessary time delay and extra wearing of tires.

[0007] Some conventional methods prefer to base its path planning on clothoids because they are closely related to the rules of movement of steered motor vehicles. A clothoid is a special type of plane curves, which has its curvature, denoted by c, changing in a linear manner along the path. Clothoids have the added benefit of the length of the path being minimized when an initial state is connected with a terminal state.

[0008]

For example, the teaching of T. Fraichard and A. Scheuer, "From Reeds and Shepp's to continuous curvature paths," introduces a concept of a clothoid turn (abbreviated as CC turn) formed by two clothoids connected by a circular arc, and teaches -tangency conditions to ensure continuous curvature transitions among CC turns and straight lines. However, a continuous curvature path typically includes a sequence of multiple CC turns and lines, and it can be time consuming to find the parameters of the optimal sequence forming the continuous curvature path.

[0009]

Instead of computing continuous curvature paths directly, another category of methods construct a continuous curvature path by first computing an RS path.

For example, method described in U.S. 2016/03133735 first determines an RS path and tries to find a continuous curvature path which stays as close as possible to the RS path. This method computes the continuous curvature path by solving optimization problems with nonlinear constraints. However, the solution to this method requires iterative process, which can be computationally prohibitive for a number of practical applications. Also, the resultant optimization problem of that method is non-convex. A feasible solution to the non-convex problem is not guaranteed solved by the iterative process, and thus the construction of the continuous curvature path may fail. Hence, this method may not be suitable for a real-time control of the vehicles. Similarly, the method described in U.S. 8428820 improves the computation efficiency, but still entails computationally expensive iterative process to construct the continuous curvature path. Hence, this method also may not be suitable for a real time control of the vehicles.

[0010]

Accordingly, there is a need to construct a continuous curvature path for controlling a motion of a vehicle from an initial to a target state in a manner suitable for real-time control of the vehicles. This work performs continuous curvature path planning, with the hope to achieve similar computational efficiency and reliability as the RS path planning.

[Summary of Invention]

[0011]

It is an object of some embodiments to provide a path planning system and method to construct a continuous curvature path connecting an initial state to a target state of a vehicle. It is another object of some embodiments to provide such a system and a method that are suitable for real-time control of the vehicles. As used herein, the continuous curvature path is as a sequential composition of multiple elementary paths connected such that the transition from one elementary path to an adjacent elementary path is continuous at the junction points of the elementary paths as well as both its first derivative and second derivatives.

Examples of an elementary path include straight line, circular arc, splines, polynomials, trigonometric functions, clothoids, and clothoid turns. In contrast, a discontinuous curvature path is a path tin with at least two elementary paths are connected at discontinuous curvature junction.

[0012]

Some embodiments perform real-time continuous curvature (CC) path planning for vehicles and other robots considering obstacle avoidance. Some embodiments are motivated by the fact that Reeds-Shepp's (RS) based path planning remains unmatched in terms of computation efficiency and reliability when compared with various CC path planning results that consider obstacle avoidance. However, the RS path is discontinuous in its curvature. To that end, some embodiments disclose methods to achieve almost equivalent computational efficiency and reliability as RS path planning by post-processing RS paths. In some implementations, the post-processing resorts to geometric interpretation of -tangency conditions, and thus not only allows a majority of planning work done offline, but also enables development of convergence guaranteed

algorithms.

[0013]

Some embodiments post-processes discontinuous curvature paths to generate a CC path, which not only satisfies the CC property, but also stays inside a neighborhood of the discontinuous curvature path to maintain obstacle clearance so that collision detection for the CC path can be avoided. In order to ensure the CC property. Some embodiments exploit the geometric insights casted by -tangency conditions to process discontinuous curvature paths.

[0014]

Specifically, some embodiments are based on recognition that according to classes and patterns of RS paths, the discontinuous curvature paths have up to four types of discontinuous curvature junctions (DCJs) : SC, CS, CC, C \ C, with S , C , and I representing straight line segments, circular arcs, and change of velocity direction, respectively. L and R specify left and right turns with + or - denoting forward or backward motion, respectively. Therefore, four distinctive schemes are devised offline for post-processing all four types of DCJs. As long as one can resolve these four DCJs, then the entire RS path between two configurations q _{s } and q _{f } , including a combination of these DCJs, can be processed to produce a CC path between q _{s } and q _{f } .

[0015]

One embodiment further identifies and simplifies critical conditions to determine parameters of disclosed schemes. These critical conditions, derived on the basis of -tangency conditions, take forms of nonlinear equations. For the SC type of DCJs, the resultant nonlinear equation can be simplified and solved analytically, and thus done offline. For the CC and C | C types, these critical equations are further simplified offline, based on geometry, into a set of simplified nonlinear equations, which can be shown the existence of solutions, and solved analytically. Thus the proposed schemes guarantee the success of treating RS paths. The simplified nonlinear equations need to be solved online, and this invention teaches algorithms that can solved them efficiently.

[0016]

In one embodiment, the RS path planning should generate an RS path, between an initial configuration q _{0 } and a final configuration q _{f } , which has a ε -clearance (at least ε far away from obstacles in all directions). In another embodiment, the RS path planning generates an RS path, between q _{(s } ^{= } q _{0 } -c and f , where c is an offset vector along the backward motion of the robot at q _{0 } . As a result, the RS path between q _{0 } and q _{f } starts with a forward straight line motion with a certain length, and thus the post-processing always deals with RS paths beginning with S segments. The schemes, to treat SC,CC,C \ C DCJs, does not include backward motion, and the resultant CC path has simpler maneuvers.

[0017]

Accordingly, one embodiment discloses a system for controlling a motion of a vehicle from an initial state to a target state. The system includes a path planner to determine a discontinuous curvature path connecting the initial state with the target state by a sequential composition of driving patterns, wherein the discontinuous curvature path is collision-free within a tolerance envelope centered on the discontinuous curvature path; a path transformer to locate and replace at least one treatable primitive in the discontinuous curvature path with a corresponding continuous curvature segment to form a modified path remaining within the tolerance envelope, wherein each treatable primitive is a predetermined pattern of elementary paths; and a controller to control the motion of the vehicle according to the modified path.

[0018]

Another embodiment discloses a method for controlling a motion of a vehicle from an initial state to a target state, wherein the method uses a processor coupled with stored instructions implementing the method, wherein the

instructions, when executed by the processor carry out at least some steps of the method. The method includes determining a discontinuous curvature path connecting the initial state with the target state by a sequential composition of driving patterns, wherein the discontinuous curvature path is collision-free within a tolerance envelope centered on the discontiguous curvature path; locating and replacing at least one treatable primitive in the discontinuous curvature path with a corresponding continuous curvature segment to form a modified path remaining within the tolerance envelope, wherein each treatable primitive is a predetermined pattern of elementary paths; and controlling the motion of the vehicle according to the modified path.

[0019]

Another embodiment discloses a non-transitory computer readable storage medium embodied thereon a program executable by a processor for performing a method. The method includes determining a discontinuous curvature path connecting the initial state with the target state by a sequential composition of driving patterns, wherein the discontinuous curvature path is collision-free within a tolerance envelope centered on the discontinuous curvature path; locating and replacing at least one treatable primitive in the discontinuous curvature path with a corresponding continuous curvature segment to form a modified path remaining within the tolerance envelope, wherein each treatable primitive is a predetermined pattern of elementary paths; and controlling the motion of the vehicle according to the modified path. [Brief Description of Drawings]

[0020]

[Fig. 1A]

Figure 1 A shows a schematic of a vehicle equipped with a front-fixed steering wheel and fixed parallel rear wheels control by some embodiments.

[Fig. IB]

Figure IB shows an example of a continuous curvature path planning (CCPP) problem considered some embodiments.

[Fig. 1C]

Figure 1C summarizes classes and driving patterns of shortest length paths, which are also termed as Reed-Shepp (RS) paths.

[Fig. ID]

Figure ID shows an example of a clothoid turn used by some

embodiments.

[Fig. IE]

Figure IE shows a graph illustrating how curvature changes along a clothoid turn.

[Fig. IF]

Figure IF illustrates an example of / -tangency condition used by some embodiments.

[Fig. 1GA]

Figure 1GA illustrates - tangency conditions between a left forward clothoid turn and a right forward clothoid turn.

[Fig. 1GB]

Figure 1GB illustrates -tangency conditions between a left forward clothoid turn and a right backward clothoid turn.

[Fig. 2A]

Figure 2A shows a block diagram of a system including a motion planner for controlling a motion of a vehicle from an initial state to a target state according to some embodiments.

[Fig. 2B]

Figure 2B shows a system diagram of the system of Figure 2 A according to some embodiments.

[Fig. 3A]

Figure 3A shows a structure of a discontinuous curvature path determined by a path planner according to some embodiments.

[Fig. 3B]

Figure 3B shows a table with examples of different RS driving patterns determined by different embodiments of an RS path planner.

[Fig. 4A]

Figure 4A shows a table with examples of some predetermined pattern of elementary paths forming the treatable primitives.

[Fig. 4B]

Figure 4B shows an example of treating a discontinuous curvature path to form a sequence of treatable primitives according to different embodiments.

[Fig. 4C]

Figure 4C shows another example of treating a discontinuous curvature path to form a sequence of treatable primitives according to different

embodiments.

[Fig. 5A1]

Figure 5A1 shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5A2]

Figure 5 A2 shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5A3]

Figure 5A3 shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5A4]

Figure 5A4 shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5B]

Figure 5B shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5C]

Figure 5C shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5D]

Figure 5D shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5E]

Figure 5E shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5F]

Figure 5F shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 5G]

Figure 5G shows one of different examples of converting different RS patterns into a sequence of treatable primitives according to some embodiments. [Fig. 6]

Figure 6 shows a flow chart of a method for processing a sequence of driving patterns of a discontinues curvature path according to one embodiment. [Fig. 7A]

Figure 7A shows a schematic of discontinuous curvature junction (DCJ) treating schemes used by one embodiment.

[Fig. 7Bl] Figure 7B 1 shows one of schematics of treating discontinuous curvature junction SC to form a continuous curvature junction according to one embodiment.

[Fig. 7B2]

Figure 7B2 shows one of schematics of treating discontinuous curvature junction SC to form a continuous curvature junction according to one embodiment.

[Fig. 7B3]

Figure 7B3 shows one of schematics of treating discontinuous curvature junction SC to form a continuous curvature junction according to one embodiment.

[Fig. 7B4]

Figure 7B4 shows one of schematics of treating discontinuous curvature junction SC to form a continuous curvature junction according to one embodiment.

[Fig. 7Cl]

Figure 7C1 shows one of schematics of treating discontinuous curvature junction CC to form a continuous curvature junction according to one embodiment.

[Fig. 7C2]

Figure 7C2 shows one of schematics of treating discontinuous curvature junction CC to form a continuous curvature junction according to one embodiment.

[Fig. 7C3]

Figure 7C3 shows one of schematics of treating discontinuous curvature junction CC to form a continuous curvature junction according to one embodiment.

[Fig. 7C4]

Figure 7C4 shows one of schematics of treating discontinuous curvature junction CC to form a continuous curvature junction according to one embodiment.

[Fig. 7D1]

Figure 7D1 shows one of schematics of treating discontinuous curvature junction C|C to form a continuous curvature junction according to one embodiment.

[Fig. 7D2]

Figure 7D2 shows one of schematics of treating discontinuous curvature junction C|C to form a continuous curvature junction according to one embodiment.

[Fig. 7D3]

Figure 7D3 shows one of schematics of treating discontinuous curvature junction C|C to form a continuous curvature junction according to one embodiment.

[Fig. 7D4]

Figure 7D4 shows one of schematics of treating discontinuous curvature junction C|C to form a continuous curvature junction according to one embodiment.

[Fig. 8A]

Figure 8A shows a function diagram of an automated parking system according to one embodiment.

[Fig. 8B]

Figure 8B shows a general structure of the motion-planning system according to one embodiment.

[Fig. 9]

Figure 9 shows a schematic of a system according to one embodiment. [Fig. 10]

Figure 10 shows a block diagram of an automated parking system according to some embodiments. [Description of Embodiments]

[0021]

Figure 1 A shows a schematic of a vehicle equipped with a front-fixed steering wheel and fixed parallel rear wheels control by some embodiments. A point R is located at the mid of rear wheels of the vehicle. A pose of the vehicle is uniquely described by a triple (χ, γ, θ) where (x, y) represent the coordinates of R in an earth frame and Θ is an orientation angle of the vehicle with respect to the positive x -axis of the earth frame. The vehicle has a wheelbase b , and a steering angle ψ .

[0022]

A kinematic model of the vehicle is given by

where a curvature κ is an extra configuration parameter in addition to (χ, γ, θ) .

[0023]

Control inputs in the vehicle model (1) are u = (ν,σ) , where v is a driving velocity of the rear wheels and σ a steering rate. The relationship between φ ,

K and σ is established as follows:

tan < . φ

K = ,σ = κ = — .

b bcos Φ

[0024]

Assume that both forward and backward motions are allowed for the vehicle and the driving velocity is bounded, whereas the angle of the steering wheel is subject to mechanical constraints, i.e.,

\ ν \< ν _{η }^ φ \≤φ _{η }^ κ \≤κ _{η∞ } = ^ ^{Χ } ^ .

b (2)

[0025] The steering rate σ is also assumed to be bounded, i.e., \ ^{σ } \ ^{≤σ } _{ηαχ }· (3)

[0026]

Any path satisfies the vehicle model (1) has continuous curvature κ along the path, and thus is a continuous curvature (CC) path. According to some embodiments, a configuration or state of the vehicle is defined by q = (χ,γ,θ, κ) .

[0027]

Figure IB shows an example defining a continuous curvature path planning (CCPP) problem considered some embodiments. Given an initial state 101

<ls and the vehicle model (1) with constraints (2) (3), find a feasible path P _{cc } between q _{s } and q _{f } such that for a finite time T _{f } }

1. Boundary conditions P _{cc } (0) ^{= i }7 _{s } and ^{= i }7/ ;

2. Feasibility Vi e [0,7}] , P _{cc }(t) e R ^{2 } * S ^{l } x [-/c _{max }, K _{m } .

[0028]

Figure 1C summarizes classes and driving patterns of shortest length paths, which are also termed as Reed-Shepp (RS) paths. In the context of RS paths, C in Figure 1C stands for a circular arc and S stands for a straight line segment; while L and R specify left and right turns with + or - denoting forward or backward motion, respectively. Subscripts denote (absolute) angular values of certain circular arc and | represents a change of moving direction.

[0029]

Figure IB also includes one example RS path in a driving pattern

R + S + R + of class CSC - X connecting 101 and 102. The example RS path consists of three segments: a circular arc 105, a line 106 and a circular arc 107. The RS path is undesirable according to certain criteria because its curvature is discontinuous at points 101, 103, 104, and 102, which incur stationary steering, time delay, and extra tire wearing.

[0030]

Solving the exact CCPP is time-consuming. Instead, some embodiments perform the CCPP by seeking for sub-optimal CC paths on the basis of RS path patterns. Furthermore, some embodiments obtain CC paths that "mimic" the corresponding RS path patterns within a bounded error refereed herein as a tolerance envelope centered on the discontinuous curvature path. Towards this end, some embodiments define a ε neighbor of a given path P as follows

B(P, ε) = (J {(*', /) e R ^{2 } :| ( ^{χ }', /) - (x, y) \< ε)

(x,y)eP

where ε > 0 is the error bound and | X | denotes the standard 2-norm of vectors X.

[0031]

Accordingly, the sub-optimal CCPP can be viewed as the following

"continuous curvature conversion" problem.

Problem 1 Given an RS path Ρ^ connecting an initial configuration q _{0 } and a final configuration ^{a }f in C , find a continuous curvature path P _{cc } such that

(i) . Pcc(0) ^{= } 0 _{> }Pcc(l) ^{= i }7 _{/ } where ά (resp., q _{f }' ) is the initial (resp., final) configuration in the augmented configuration space C' with null curvature

K _{0 } = 0 (resp., *> = 0 );

(ii) . _{C } satisfies the model (l)-(3) for certain [ , σ] e U ;

(iii) . P _{C }c remains in the ^ -neighborhood B(P _{RS }, s) for a given ε > 0.

[0032]

Some embodiments are based on realization that Problem 1 leverages the design freedom v to post-process S paths into CC paths. The problem is well-posed, i.e., its solutions are guaranteed to exist, which can be loosely shown by the following argument. By treating the path length 5 as the independent variable, system ( 1 ) can be rewritten dx/ds = cos Θ, dyldx = sin Θ, d&ds = κ, didds = σΐν.

Let us define a velocity profile v _{cc } along P _{RS } , where v _{cc } takes a small value v around discontinuous curvature junctions (DCJs). This means that around the DCJs, didds is fairly large, and thus instantaneous curvature changes of P _{rs } can be approximated well. Each v _{cc } gives rise to a CC path P _{cc } , which can be obtained by integrating system (1). By choosing different v , one can obtain a family of CC paths P _{cc }(v) , which are parameterized by v and has continuous dependence on it. As v— » 0 , P _{cc } (v )→ P^ . Hence, given a neighborhood Β^Ρ^, ε) , one can always find v > 0 such that as long as 0 < v < v around DCJs, the CC path P _{cc }(y) satisfies conditions in Problem 1.

[0033]

Some embodiments are based on realization that disjunctions can be resolved using clothoid and/or clothoid turns. A clothoid can be a curve satisfying (1), and its curvature κ varies linearly with respect to its arc length s , i.e.,

K S) = <J _{m!0l }s , where r _{max } is a sharpness of the clothoid. In one embodiment, configurations of a clothoid, defined over 0 < s≤ K _{mSK }l cr _{max } , can be determined by integrating system (1) along the clothoid arc

^{s } 7T 2 f ^{5 } ft 2

where C _{f }(s) = ] _{q } cos— τ dr and S _{f }(s) = J _{q } sin— τ άτ are the Fresnel cosine and sine integrals, respectively. The clothoid for a given cr _{msx } can be pre-computed and is plotted in Figure ID as an arc connecting q _{s } and q _{x } .

[0034]

A deflection of the clothoid represents a difference in orientations of q _{s } and q _{x } . For fixed * ^{" } _{max } and cr _{max } , the clothoid has important properties

1. a constant deflection S _{c } -- ^ _{3Χ } (2σ _{π13Χ })

2. a constant arc length s = K _{MA }J cr _{max } .

3. a uniquely determined solution q{s) .

[0035]

These properties can be pre-computed and stored in memory for

computational efficiency.

[0036]

Figure ID shows an example of a clothoid turn used by some embodiments. Specifically, Figure ID plots the left forward clothoid turn 111 , exemplifying a case 2S _{C }≤ δ < 2S _{C } + π , consists of

1. a clothoid from q _{s } to q _{x } , with a sharpness < _{max } ;

2. a circular arc of radius f ^{" }^ and of angle S _{f } - 2S _{C }, starting from q _{x } and ending at q _{2 }

3. a second clothoid starting from q _{2 } with sharpness - c _{max } and ending at q _{f } .

When δ = 0 , a clothoid turn reduces to a line segment of length 2R _{n } sin μ , where R _{Q } and μ are defined later. For 0 < δ < 2S _{C } , a clothoid turn contains a clothoid of a sharpness σ < _{Μ3χ } and a symmetric clothoid of a sharpness - σ , where the shar ness σ is given by

and the length of each clothoid is . Clothoid turns with

2S _{C } + π≤ S < 2π correspond to left backward motions, and can be treated similarly. It is clear that a clothoid turn can be parameterized by its deflection δ , and a line segment can be parameterized by its length / .

[0037]

Figure IE shows a graph illustrating how curvature changes along a clothoid turn such as the clothoid turn is given by an arc 1 1 1 from q _{s } to in

Figure ID, where the arc represents a left forward clothoid turn, equivalently a left turn and forward movement. Depending on its deflection δ ^{~ } 0 _{S } ^{~ } Of , a clothoid turn can include up to two clothoids and one circular arc.

[0038]

One concept in simplifying the sub-optimal CCPP from geometric perspectives is CC Circle. Figure ID shows four CC circles associated with a configuration q _{s } . Particularly, 121 C fe) , 122 C _{/ } ^{" } (^ _{5 } ) , 123 C ^{+ } (# _{5 }) , and

124 C _{r } (q _{s } ) correspond to the left forward clothoid turn 1 11, a left backward clothoid turn, a right forward clothoid turn, and a right backward clothoid turn, respectively. Note that a clothoid turn, enclosed by its corresponding CC circle, always enters and leaves the corresponding CC circle in the direction which forms a μ angle with tangent vectors at the entry and exit points.

[0039]

Notably, each CC circle has exactly the same center as the circular arc of its respective clothoid turn. For instance, Ω , the center of the CC circle 121 Q ^{+ } ( s ) · _{> } coincides with the center of the circular arc which is tangential to the vehicle orientation at q _{x } . Its coordinates are given by

[0040]

For fixed *r _{max } and a _{mwi } , the CC circle radius and //-angle are given by

where μ is the angle between the vehicle orientation and a tangent vector of the

CC circle 121 at q _{s } .

[0041]

In one embodiment, one can perform the sub-optimal CCPP using μ- tangency conditions, which determine combination of clothoid turns and lines. There are three possible combinations in continuous curvature paths fitting the driving patterns in Figure 1C: SC,CC,C \ C . Accordingly, there are three types of / -tangency conditions.

[0042]

Figure IF illustrates SC -tangency condition between a line 133 and a clothoid turn, which uniquely determines a CC circle 132. The clothoid turn begins with 101, and transits to the line 133 at <Z _{\ } 131. The SC /-tangency condition between the clothoid turn and the line 133 is that an angle between a tangent vector 134 of the CC circle 132 at q _{x } 131 and the line 133 is μ .

[0043]

Figure 1GA illustrates CC -tangency conditions between a left forward clothoid turn and a right forward clothoid turn, corresponding to CC circles 141 and 142, respectively. The CC -tangency conditions specifies

1. the CC circles 141 and 142 have to be tangential;

2. a distance L( _{l 2 } ) between centers 143 and 144 of the CC circles 141 and 142 is equal to 2R _{n } .

[0044]

Figure 1GB illustrates C \ C -tangency conditions between a left forward clothoid turn and a right backward clothoid turn, corresponding to CC circles 145 and 146, respectively. The C \ C -tangency conditions specifies that

1. the CC circles 145 and 146 intersect with each other, and at an intersecting configuration 149, tangent vectors 147 and 148, of the CC circles 145 and 146, respectively, form angles of μ with respect to an orientation 149 of the vehicle at 149;

2. a distance Li L^l^) between centers 151 and 152 of the CC circles 145 and 146 should be 2R _{n } cos μ .

[0045]

Some embodiments perform real-time continuous curvature (CC) path planning for vehicles and other robots considering obstacle avoidance. Some embodiments are motivated by the fact that Reeds-Shepp's (RS) based path planning remains unmatched in terms of computation efficiency and reliability when compared with various CC path planning results that consider obstacle avoidance. However, the RS path is a discontinuous curvature path. To that end, some embodiments disclose methods to achieve almost equivalent computational efficiency and reliability as RS path planning by post-processing RS paths. In some implementations, the post-processing resorts to geometric interpretation of -tangency conditions, and thus not only allows a majority of planning work done offline, but also enables development of convergence guaranteed algorithms.

[0046]

Figure 2A shows a block diagram of a system 200 including a motion planner 250 for determining a motion of a vehicle from an initial state to a target state according to some embodiments. The system includes a path planner 210 to determine a discontinuous curvature path 215 connecting the initial state with the target state by a sequential composition of driving patterns. The path planner 210 uses principles employed by the RS based path planning to produce such a discontinuous curvature path that is collision-free within a tolerance envelope centered on the discontinuous curvature path.

[0047]

Some embodiments are based on recognition that some, but not all patterns of the discontinuous curvature path can be replaced with continuous curvature segments remaining within the tolerance envelope. These embodiments are based on analysis of four types of discontinuous curvature junctions (DCJs):

SC,CS,CC,C I C, with S , C , and | representing straight line segments, circular arcs, and change of velocity direction, respectively. L and R specify left and right turns with + or - denoting forward or backward motion, respectively.

[0048]

Based on this analysis, some embodiments identify, i.e., determined offline, different patterns of elementary paths that can be replaced with corresponding patterns forming continuous curvature segments. That is, finding a continuous curvature path staying within a specified neighborhood of paths in these

predetermined patterns is guaranteed. Such a predetermined pattern of elementary paths that can be replaced with continuous curvature segment is referred herein as

"treatable primitive." For example, a sequence of a straight-line segment connected to a first arc segment connected to a second arc segment forms a treatable primitive that can be replaced with a combination of lines, clothoid, and clothoid turns forming a continuous curvature segment. Similarly, a sequence of a first arc segment connected to a second arc segment connected to a straight-line segment also forms a treatable primitive that can be replaced with a pair of clothoid turns and a line if necessary.

[0049]

The types of the treatable primitives can vary, but in order to speed up the online calculation of the CC path, the embodiments predetermine those treatable primitives and corresponding structure of the continuous curvature segments. For example, a sequence of an even number of circular arc segments connected to the straight-line segment can be replaced with a sequence of the even number of clothoid turns. The parameters of the clothoid turns forming the continuous curvature segment are determined online, but because the structure of the continuous curvature segment is predetermined, some embodiments simplifies the online computation of the parameters.

[0050]

To that end, the system 200 includes a path transformer 220 configured to locate and replace at least one treatable primitive in the discontinuous curvature path with a corresponding continuous curvature segment to form a modified path remaining within the tolerance envelope, and a controller 240 to control the motion of the vehicle according to the modified path 235.

[0051]

In one embodiment, the path transformer converts the treatable primitives forming only part of the discontinuous curvature path 215. Thus, in that embodiment, the modified path 235 is also discontinuous but with less

discontinuous curvature junctions then the path 215. In alternative embodiments, the system 200 aim to produce the discontinuous curvature path forms only by treatable primitives, such that the modified path 235 is a continuous curvature path.

[0052] For example, in one embodiment, the path planner 210 is a basic Reeds-Shepp (RS) path planner forming the discontinuous curvature path using a combination of a straight-line segment and a circular arc tangentially connected at discontinuous curvature junctions. This embodiment simplifies the

implementation of the path planner 210, however, does not guarantee that the discontinuous curvature path 215 includes only treatable primitives. To that end, this embodiment also includes a converter 225 to convert the sequential composition of driving patterns of the discontinuous curvature path 215 into another discontinuous curvature path which is formed by a sequence of treatable primitives 225. For example, the converter 220 can insert a new straight-line segment in the path 215 to form the path 225. In such a manner, after the path transformer 230 transforms path 225, the modified path 235 is a continuous curvature path. Examples of the operations of the converter 220 are provided below.

[0053]

Additionally, or alternatively, the path planner 210 can be a modified RS path planner that produces such a discontinuous curvature path 215 that is already formed only by a sequence of treatable primitives. For example, in one

embodiment, the path planner imposes a constraint on the discontinuous curvature path requiring at least one of an initial and a final segment of the discontinuous curvature path being a straight-line segment.

[0054]

Figure 2B shows a system diagram of the system 200 of Figure 2 A according to some embodiments. The system 200 can includes a memory 201 that stores a finite number of treatable primitives 211, discontinuous curvature junction (DC J) treat schemes 212, driving pattern treat schemes 213, and RS path treat schemes 214. All these schemes are present in the memory in the form of codes. A processor 202 implements the motion planning 250, on the basis of inputs states 221 indicating initial and final configurations of the vehicle, determines an RS path 222, and determine a CC path 223. An actuator 203, e.g., a motor, controls the motion of the vehicle according the CC path.

[0055]

Some embodiments post-processes discontinuous curvature paths to generate a CC path, which not only satisfies the CC property, but also stays inside a neighborhood of the discontinuous curvature path to maintain obstacle clearance so that collision detection for the CC path can be avoided. In order to ensure the CC property. Some embodiments exploit the geometric insights casted by -tangency conditions to process discontinuous curvature paths. To that end, in some embodiments, the path transformer 230 determines parameters of the continuous curvature segment to satisfy a μ-tangency constraint.

[0056]

Figure 3 A shows a structure of discontinuous curvature path 215

determined by path planner according to some embodiments. In those

embodiments, the discontinuous curvature path is an RS path 310 formed by one or a sequence of multiple RS driving patterns 311, 312. Example of RS driving patterns 31 1-312 include patterns 199 shown in Figure 1C. As can be seen on Figure 1C, each RS driving pattern includes a sequence of elementary paths 321-322. Path planner selects the elementary paths from a finite set C,.?}, which represents a circular arc and line segment, respectively.

[0057]

Figure 3B shows a table 330 with examples of different RS driving patterns determined by different embodiments. For example, one embodiment determines basic patterns 332 of the RS path. The patterns 332 does not guarantee that the RS path formed by these patterns includes only treatable primitives. Alternatively, some embodiments produce patterns 334, 336, and 336 by deliberately inserting an extra elementary path formed by a straight-line segment S. For example, one embodiment inserts a straight-line segment S in the front of driving patterns 332 to form the pattern 334. Another embodiment inserts a straight-line segment S at the end of driving patterns 332 to form the pattern 336. Yet another embodiment inserts a straight-line segment S at the end and at the beginning of driving patterns 332 to form the pattern 338. Modified RS path planners with insertion of S may simplify post-processing of resultant RS paths, and produce a continuous curvature path with less maneuvers.

[0058]

Figure 4A shows a table with examples of some predetermined pattern of elementary paths forming the treatable primitives 410. This example discloses a set of treatable primitives 410 {S, SCS, SCC, CCS, SC \ C,C \ CS} , which can be treated by one step of treating schemes used by some embodiments. As can be seen in this Figure 4A, the patterns of the treatable primitives 410 starts and/or ends with a straight-line segment S. To that end, the patterns 334, 336, and 336 simplify locating treatable primitives in the discontinuous curvature path generated by the path planner.

[0059]

Figures 4B and 4C show examples of locating and treating the treatable primitives in the discontinuous curvature path according to different

embodiments. For example, some embodiments generate a CC path by locating a finite set of treatable primitives in the sequence of patterns 420 forming

discontinuous curvature path, which represents shortest RS driving patterns that can be converted into a CC driving pattern by one single step.

[0060]

For example, the embodiment of Figure 4B identifies 430 treatable primitives in each RS driving pattern of the RS path independently from each other. Conversely, the embodiment of Figure 4C can process the entire RS path as a whole in order to locate the treatable primitives. This allows to modify the RS path in order to simplify the selection of treatable primitives. For example, this embodiment can modify 440 the RS path by inserting, splitting, or removing an elementary path from the sequence of patterns 440. Also, this embodiment can partition 450 the modified sequence of RS driving patterns into a new sequence of patterns, wherein each new pattern is a treatable primitive, and process 460 the new sequence to produce the continuous curvature path.

[0061]

Figures 5A-G show different examples of converting a different RS pattern into a sequence of treatable primitives according to some embodiments.

[0062]

Figure 5A1 describes a treatment scheme for the RS driving pattern CSC 501a, which is first augmented to the driving pattern 502a by inserting and/or appending S . The word S _{a } represents a line segment with certain length a . Given 502a, the S between two are split into two S _{a } 503b. It is clear that

2

the driving pattern after treatment S _{a } I S _{a }CS _{i }S _{1 }CS _{a } \ S _{a } consists of two treatable

2 2

primitives S , and another two treatable primitives S _{a }CS _{a } 504a and S _{a }CS _{a }

2 2

505a.

[0063]

Figure 5A2 presents a treatment scheme for the RS driving pattern SCSC 511a, which is first augmented to the driving pattern 512a, by appending a driving word S in the end. The driving pattern 512a is further processed by splitting a second S into two $a to give 513a. It is clear that the modified

2

driving pattern 513a consists of treatable primitives SCS 514a, SCS 515a, and an S .

[0064]

Figure 5A3 teaches a treatment scheme for the RS driving pattern CSCS 521a, which is first augmented to the driving pattern 522a, by appending a driving word S in the beginning. The driving pattern 522a is further processed by splitting a third S into two S _{a } to give 523a. It is clear that the modified

2 driving pattern 523a includes treatable primitives SCS 524a, SCS 525a, and an S .

[0065]

Figure 5A4 depicts a treatment scheme for the RS driving patterr SCSCS 531a, which is processed by splitting a second S into two $a to give 532a. It

2

is clear that the modified driving pattern 532a consists of treatable primitives

SCS 532a and SCS 534a.

[0066]

Figure 5B discloses a treatment scheme for the RS driving pattern

SC I C I CS 501b, which is processed by splitting a second C into two C _{a } to

2 give 502b. It is clear that the modified driving pattern 502b consists of treatable primitives SC \ C 503b and C \ CS 504b.

[0067]

Figure 5C discloses a treatment scheme for the RS driving pattern SCC \ C 501c, which is augmented by inserting 5" behind a first driving word C to give 502c. It clear that the modified driving pattern 502c consists of treatable primitives SCC 503c and SCCS 504c.

[0068]

Figure 5D discloses a treatment scheme for the RS driving pattern CC | CC 50 Id, which is processed by inserting S between a second and a third C to give a modified driving pattern 502d. It is clear that the modified driving pattern 502b consists of treatable primitives CCS 503d and SCC 504d.

[0069]

Figure 5E discloses a treatment scheme for the RS driving pattern

SC \ CC \ CS 50 le, which consists of treatable primitives SC \ C 503 e and C \ CS 504e.

[0070] Figure 5F discloses a treatment scheme for the RS driving pattern C | CSCS 50 If, which is processed by splitting a first S into two S _{a } to give 502f. It is

2

clear that the modified driving pattern 502f consists of treatable primitives C | CS

503 f and SCS 504f.

[0071]

Figure 5G discloses a treatment scheme for the RS driving pattern

C I CSC I C 50 lg, which is processed by splitting a first S into two $a to give

2

502g. It is clear that the modified driving pattern 502g consists of treatable primitives SC | C 504g and C | CS 503g.

[0072]

Figure 6 shows a flow chart of a method for processing a sequence of driving patterns 601 forming a discontinues curvature path according to one embodiment. In this embodiment, the converter 220 iteratively tests each tuple of elementary paths forming the discontinuous curvature path for being a treatable primitive and inserts a new straight-line segment into a tuple not representing the treatable primitive. For example, the driving sequence, always untreated in the beginning, is parsed 602, where a first tuple of driving words in the untreated driving sequence are fetched to test whether it belongs to the set of treatable primitives or not 603. If it does, the rest of the driving sequence will be the new untreated driving sequence, otherwise, the fetched tuple is processed by adding S to the front the tuple or after the second driving word of the tuple to yield a treatable primitive; and the new untreated driving sequence includes the last word of the tuple, followed by the old untreated driving sequence. If the whole driving sequence has been treated, the treating process stops; otherwise, the above steps fetching 602, testing 603, adding 604, and verifying 605 are repeated.

[0073]

Figure 7A shows a schematic of discontinuous curvature junction (DC J) treating schemes used by one embodiment. These schemes lead to identification of treatable primitives. Classes of RS driving patterns in Figure 1C reveal that DCJs occur at connections among straight line segments and circular arcs;

specifically, DCJs can generally be categorized as SC , CS , CC and C \ C patterns. Accordingly, DCJ treat schemes contain the treatment of SC 703a, C | 704a, and CC 705a.

[0074]

In some embodiments, a clothoid satisfying system dynamics (1) for a pair of fixed K ^{" } _{max } and <r _{mgx } , is pre-computed and stored in memory 701a, and can be reused online to treat DCJs. Similarly, centers of clothoid circles and radius of clothoid turns702a corresponding a clothoid starting at q _{s } can also be computed offline and stored in memory for online use.

[0075]

Figures 7B1, 7B2, 7B3, and 7B4 show schematic of treating a

discontinuous curvature junction SC to form a continuous curvature junction according to one embodiment. The discontinuous curvature junction SC

corresponds to a transition between a line and a circular arc. By adopting an appropriate coordinate framework, this corresponds to solve the following problem.

[0076]

Problem 2 Let a simple RS driving pattern connecting an initial configuration o ^{= } (Χ _{0 },0,0, ) 701 b and a final configuration

q _{f } = O _{m }L sin0, i sin.0,0, *r _{max } ) 799b consist of[(i)J

1. A line segment S 71 lb described by -∞< x _{0 } < x < 0 and y = ;

2. A circular arc C 712b, which is tangent to S at the origin (0,0) 713b.

The radius and angle of C are given by and θ , respectively. Construct a continuous curvature path P _{cc } connecting o ^{an }d q _{/ } which remains in B(P^ , έ) while satisfying conditions (i)-(ii) in Problem 1.

[0077]

A simple SC RS driving pattern is depicted in Figure 7B 1. Given a value of ^max such that 0 > 2S _{C } , the following SC transition scheme guarantee to make CC connection between S and C (see Figure 7B-1 and Figure 7B-2 for details)

1. A forward straight line segment that starts from the initial configuration q _{0 } = (x _{0 } ,0,0,0) 701b and ends at the configuration q _{s } = (x ^{* }, 0,0,0) 702b with x _{0 } < x ^{* } < 0 .

2. A forward right clothoid turn (CT) that starts from the configuration q _{s } = (x ^{* }, 0,0,0) 702b and ends at q _{x } = (x _{x },y _{x }, S,G) 703b with a deflection

S < 0 ;

3. A forward left CT with a deflection - S that starts from q _{l } 703b and ends at a configuration q _{2 } - (x _{2 },y _{2 }, ,0) 704b.

4. A forward left clothoid that connects q _{2 } 704b and

q _{3 } = (x _{3 }, y ,S _{c }, K _{max } ) 705b, such that the point (x _{3 } , y _{3 } ) e C .

[0078]

The following procedure can be followed to determine parameters in the SC treating scheme, i.e., q _{s } , q _{2 } -

1. As depicted by Figure 7B-2, both the CC Circle C, ^{+ }fa _{2 }) 717b and the circular arc C 712b share the same center. Hence, Ω _{3 } .

Coordinates of _{2 } are then computed by

2. Considering that the configuration q _{2 } is the end configuration of the forward left CT, one can immediately determine the center Ω _{2 } of the CC Circle C (q _{2 }).

3. For a certain x , the CC circle ^(q _{s }) has a center Qj( ), which is uniquely defined by x . On the other hand, the CC circle C*{q _{s }), corresponding to the first right CT starting from q _{s } , should be tangent to

C (q _{2 }). Therefore, x ^{* } can be determined by solving the following nonlinear equation

ΡΩ _{! }-Ω _{2 }Ρ = 2^ _{Ω }, - <x ^{* }<0. _{(8) }

*

4. With x , Ωι ^{= }(¾ »>Ώ ) and Ω _{2 } = ( ^{χ } _{Ω }^ , y _{a }^ ), we compute the auxiliar angle a as follows: a =

[0079]

The deflection of the two CTs are given by δ and - δ accordingly, where δ = α-π/2+μ. The length of the straight line segment is given by

[0080]

The relationship between the CC Circles 715b, 716b, 717b is illustrated in Figure 7B2. Previously-established computation procedure can be summarized as a CCP-SC method described in Fig. 7B3.

[0081]

The CCP-SC method takes the initial configuration q _{0 } and the final configuration q _{f } , along with the error bound ε , as the inputs. The initial design parameter is selected in such a way that Θ > 2S _{C } is satisfied. First, the

CCP-SC algorithm computes (lines 2 to 9) the appropriate ^max such that

I y _{2 } \≤ ε to guarantee that P _{cc } c £(P^, ε) . After obtaining cr _{max } , the CCP - SC

*

algorithm computes (lines 10 to 14) the corresponding x . As explained in the

*

computation procedure, x should satisfy both Equation (8) and the path constraint x _{0 } < x < U . Some implementations reduce (8) to the equation in Line

*

13 based on simple geometric analysis. Once a proper x is obtained by solving the equation, the iteration loop terminates and the CCP-SC algorithm computes and outputs the CC path P _{cc } accordingly (lines 15 to 16); otherwise, the method increases the value of _{max } and updates the value of x ^{* } .

[0082]

Figure 7B4 plots a solution of the proposed SC treat scheme for a specific example.

[0083]

Figures 7C1, 7C2, 7C3, and 7C4 show schematic of treating discontinuous curvature junction CC to form a continuous curvature junction according to one embodiment. The discontinuous curvature junction CC is formed by juncture between two circular arcs. By adopting an appropriate coordinate framework, this corresponds to solve the following problem.

[0084]

Problem 3 Let be a CC type of RS path starting from the initial configuration q _{0 } - (0,0,0) 701c and ending at q _{f } 799c. P _{RS } consists of

1. A circular arc C, 711c which starts from the origin and ends at A

713c. The radius and angle of C, are given by ^ ^{" } _{max } and _{χ } ,

respectively.

2. A circular arc C _{2 } 712c ending at the final configuration

q _{f } 799c which is tangent to C, at A , moves in the same direction as C _{x } , and has the radius and angle given by ^ ^{" }max and 0 _{2 5 } respectively.

Construct a continuous curvature path pattern P _{cc } to achieve the following requirements: (i) P _{cc } connects the null-curvature augmented configurations

4 _{0 }, = (0,0,0,0) and q _{r } = (x _{f },y _{f } , 0 _{f } ,O) _{; } (ii) P ^ B(P£ , ε) _{; } (Hi) P _{cc }

satisfies conditions (i)-(ii) in Problem 1.

[0085]

For illustration purpose, this invention assumes that ^ ^{G } [®> ^{π }] (counterclockwise) and ^{e π }^ (clockwise). This scenario is depicted in Figure 7C1. The centers of ^ and ^{c }^ are obtained as

[0086]

Based on the coordinates of O, and 0 _{2 } , the final configuration q _{f } of the arc C _{2 } (hence ) is then obtained as

Q/ = sin θ _{χ } + f - _{ax } sm(0 _{2 } - θ _{λ } ),

[0087] The embodiments proceeds to the construction of the desired P _{cc } that connects q _{0 }' and q _{f }' . Solving Problem 3 is similar to computing a continuous curvature path that connects the configurations q _{0 } = (0,0,0,0) and q f (9). To this end, we propose a continuous curvature path pattern S \ CC to solve Problem 3, where S and C denotes line segments and CT, respectively. As shown in Figure 7C2, the proposed CC path is composed of:

1. S : a backward line segment that starts from q _{0 } 701c and ends at a configuration q _{x } = (x ^{* }, 0,0,0) 702c with ε≤ ^{χ* } < 0 ;

2. CT _{X } : a forward left CT with a deflection δ _{χ } that starts from q 702c and ends at a configuration q _{2 } ^{= } ( ^{χ }ι » Ί > δ _{\ } ?0) 703c;

3. CT _{2 } : a forward right CT with a deflection S _{2 } = 9 - θ _{2 }— δ _{λ } that starts from q _{2 } 703c and ends at q _{f } 799c.

Given an appropriate a _{max } , the following method can be employed to compute transition parameters x ^{* } , ^{χ } , y , S and S _{2 } .

1. With CT _{2 } ending at q _{f } , given by (9), we have a center

^{Ω } _{2 } ^{= } (½ _{2 } ' ½ _{2 } ) of a CC Circle C ^{~ }(q _{f } ) 716c.

2. Since CT _{X } starts from q _{x } 702c, a center Q _{j }( ^{* }) of CC Circle Ci(<]i) 715c is a vector- valued function of i .

3. CC Circles Q ^{+ } (q _{x } ) 715c and C ^{~ }(q _{f } ) 716c must be tangent to each other. Towards this end, x' shall satisfy

PQ _{x }(x ^{* }) - Q _{2 }P = 2R _{n }, x ^{* } < 0, ( _{10 }) where R _{n } is given by (6).

* Ω, = (x _{n } , y _{n } )

4. Given ^{x } solved from (10), ^{! 1 } can be determined. Compute the auxiliary angle ^{a } as follows

[0088]

The deflections of CT and CT _{2 } are given by δ _{λ } = a + π/2 - μ, and θ _{ι } - θ _{2 } - δ _{ι } , respectively, while the length of the line segment is obtained x | .

[0089]

Similarly, one embodiment can also use a CC \ S path pattern to treat CC , which includes two CTs and a backward motion in a straight line segment. Due to its symmetry to the aforementioned S | CC path pattern, CC | S can also be computed efficiently in a similar manner.

[0090]

The computation procedure of the S \ CC path that connects q _{0 }' and q _{f }' is summarized as the CCP-CC method in Fig. 7C3.

[0091]

Similar to the CCP-SC method, the CCP-CC method takes q _{0 } , q _{f } and ε as the inputs. Intuitively, the CCP-CC algorithm computes an

*

appropriate x in an iterative manner such that (10) can be satisfied. For such a purpose, at the k -th iteration step, the CCP-CC computes R _{n } and μ with respect to cr*^ = σ° _{3χ } + &Δσ (Lines 3 to 4). Next, the center Ω _{2 } of C _{r } (q _{f }) is obtained (Line 5). On the other hand, from a simple geometric reasoning, the center Ω _{χ } of C q _{x }) can be obtained from x ^{* } and Ω _{0 } of C (q _{Q }) (Lines

6-7). Based on the coordinates of , (10) can then be reduced to the equation presented in Line 8 of the CCP-CC algorithm. Solving this equation yields a * _{* } ι * ι

valuation of x ; and x can be accepted whenever \ x \< ε (Lines 9-12).

Afterwards, the desired P _{cc } can be computed in terms of x accordingly.

[0092]

Figure 7C4 plots a solution of the proposed CC treat scheme for a specific example.

[0093]

Figures 7D1, 7D2, 7D3, and 7D4 show schematic of treating discontinuous curvature junction C|C to form a continuous curvature junction according to one embodiment. The discontinuous curvature junction C\C is formed by juncture between two circular arcs with a reverse movement. By adopting an appropriate coordinate framework, this corresponds to solve the following problem.

[0094]

Problem 4 Let be an RS path starting from the initial configuration q _{0 } = (0,0,0) 701d that consists of

1. A circular arc C _{x } 71 Id that starts from q _{0 } with the radius and angle of C _{x } given by and θ _{λ } , respectively.

2. A circular arc C _{2 } 712d that is tangent to C _{x } at point A 713d. The vehicle moves in the opposite direction as in C, until the final

configuration q _{f } = (x _{f },y _{f },0 _{f }) 799d. The radius and angle of C _{2 } are given by and θ _{2 } , respectively.

Construct a continuous curvature path P _{cc } such that: (i) P _{cc } connects the null-curvature augmented configurations q _{0 }, = (0,0,0,0) and q _{f }< = (x _{f } , y _{f } , 0 _{} } ,0) ;

(ii) P c B(P^ ^{C } , ε) ; (iii) P _{cc } satisfies conditions (i)-(ii) in Problem 1.

[0095]

For clarity, this embodiment assumes without loss of generality that θ _{λ } <Ξ [0,π] (counterclockwise) and θ _{2 }≡[0,π] (counterclockwise). Centers of circular arcs C, 71 Id and C _{2 } 712d, denoted as O _{x } and 0 _{2 } , respectively, are computed identically as the CC DCJ circumstance.

[0096]

As shown in Figure 7D1, the RS path P¾ ^{c } ends at the configuration

<lf ^{= } ( ^{2 }* _{m }L sin 0 _{X } - ΛΓ^ _{Χ } sin(<, + # _{2 } ),

KmL - cos Θ _{λ } + f ^ _{x } cos^ + < _{2 } ), ¾ + ¾ ). _{(11} }

[0097]

One embodiment discloses a CC treat scheme, illustrated by Figure 7D2, that constructs a CC path of the pattern S\C\C to connect q _{o }' = (q _{o },0) and

Qf ~ (<7/Ό) in (11). In particular, the proposed CC path includes

1. S : a backward line segment that starts from q _{0 } 70 Id and ends at a configuration q = {x ^{* } ,0,0,0) 702d, with -ff < * < 0;

2. , : a forward left CT with a deflection that starts from q _{x } 702d and ends at a configuration 2 ^{~ } ( i > > >0) 703d;

3. _{2 } : a backward right CT with a deflection t> _{2 } - # _{1 } + θ _{2 } - δ _{χ } that starts from q _{2 } 703d and ends at q _{f } 799d.

Given a fixed CT _{max } , the following procedure can be employed to determine transition parameters x" , x _{x }, y _{x }, δ and S _{2 }. Please refer to Figure 7D2 for illustration.

1. With CT _{2 } ending at q _{f }, given by (11), we have coordinates of a centerof C _{r } ^{+ }(q _{f }) 716das Ω _{2 } = ( _{¾ }, _{¾ }).

2. With CT _{X } starting at q a center Ω,(χ ^{* }) of CC circle C*{q _{x }) 715d is given by a vector- valued function of x .

3. Based on μ- tangency conditions, CC Circles C^(q _{x } 715dand C _{r }(<lf) 716dmustbe - intersecting. That is, x ^{* } shall ensure ΡΩ. _{λ }(χ ^{* })-Ω. _{2 }Ρ = 2R _{a }G,os μ, -ε≤χ ^{* } <0.

4. Given x ^{* } solved from (12), Q _{j } = (x^ , y _{n } ) can be determined. Compute the auxiliary angle a as follows:

[0098]

The deflections of CT _{X } and CT _{2 } are given by δ _{λ } = a + π/2, and δ _{2 } =θ _{ι }+θ _{2 }-δ _{ι }, respectively; whereas the length of the line segment is given by

Ι*Ί·

[0099]

Similar to the CC treat scheme dealing with a CC DCJ, a continuous curvature C\C\S path pattern can also be employed as an alternation of the previously established S\C\C path to treat a C\C DCJ.

[0100]

The aforementioned computation procedure of the S\C\C can be summarized in a CCP -C\C method in Fig.7D3.

[0101]

The inputs of the CCP-C\C algorithm include q _{0 }, q _{f } and ε.

Intuitively, the CCP-C\C algorithm determines x ^{* } iteratively such that it solves (12). Towards this end, the CCP-C\C first computes R _{n } and μ with respect to cr* _{ax } = σ^ _{χ } + kAa at the k -th iteration step (Lines 3 to 4). Next, the center Ω _{2 } of C _{r } ^{+ }(q _{f }) is obtained (Line 5). In addition, the center Ω, of C^{q _{x }) can be computed based on the center Ω _{0 } of C (q _{0 }) (Lines 6-7).

Facing the coordinates of Ω _{1 } , (12) can then be reduced to the equation presented in Line 8 of the CCP-C | C algorithm. Solving this equation yields a valuation of

X ; and the value of X can be accepted whenever | x |< ε is satisfied before the iteration terminates (Lines 9-12). Afterwards, the desired P _{cc } can be

*

computed in terms of x accordingly.

[0102]

Figure 7D4 plots a solution of C \ C treat scheme for an example.

[0103]

Figure 8A shows a function diagram of an automated parking system according to one embodiment. Environment mapping and localization block 801 constructs or updates a map of a parking space, and determines the current location of the vehicle by sensing the environment and vehicle operation condition. For example, an inertial measurement unit, which may include 3 -axis accelerometer(s), 3 -axis gyroscope(s), and/or magnetometer(s), can be used to sense the vehicle operation. A global positioning system sensor can be used to provide position and velocity of the vehicle. Sensors to sense the environment 500 can be video cameras capturing obstacles including other vehicles,

pedestrians, and buildings, ultrasonic/radar sensors detecting distance between the vehicle and obstacles, etc. Final state selection block 801 selects a final state for a parking spot to park the vehicle by identifying parking lot candidates, and sends the final state to a motion planning block 803. In one embodiment, the available parking spots are tracked by a separate system associated with the management of the parking garage.

[0104]

Additionally, or alternatively, the parking spots can be detected using the sensors 803 of the automated parking system. In one embodiment, the motion planning block checks to determine whether the final state is parkable, i.e., there is a feasible path to the parking spot, or not, and notifies the final state selection block 802 the check result. If the final state is not parkable, the target selection block 802 selects another final state for evaluation. In another embodiment, the final state selection block 501 can also evaluate whether the final state is parkable, and only sends a parkable final state to the motion planning block.

[0105]

If the final state is parkable, then the motion planning 803 initiates a complete motion planning procedure to determine a reference trajectory 841 based on the vehicle models 810, the initial and final states of the vehicle, and the map of the parking space. In one embodiment, the reference trajectory defines profiles of the vehicle velocity and steer angle over time. In another embodiment, the reference trajectory defines the profile of the vehicle state over time.

[0106]

Given the reference trajectory 841, vehicle controllers and actuators 804 determine and exert control commands to enforce the vehicle state track the reference trajectory 841 if the reference trajectory is state profile, or to enforce the vehicle velocity and steer angle track the reference trajectory if the reference trajectory is the vehicle velocity and steer angle profiles. In one embodiment, the control commands could be gas pedal pressure or steering torque. The vehicle controller/actuators may also use signal 843 to determine control commands. The signal 843 can be a measured steering angle, or measured currents of motors moving the steering wheel or the gas pedal.

[0107]

Figure 8B shows a general structure of the motion-planning system 803 according to one embodiment. The motion-planning system 803 includes at least one processor 870 for executing modules of the motion-planning system 803. The processor 870 is connected 871 to a memory 880 that stores geometrical

information 881 such as geometry of the vehicle and a map of the parking space. The memory 880 can also store the models of the vehicle 882 such as a kinematic model of the vehicle and a dynamic model of the vehicle. The memory 880 can also store the internal information 883 of the motion planner, including, but not limited to, an initial state of the vehicle, a final state of parked vehicle, cost function, values of each computed state, pre-computed clothoid and its deflection angle and deflection for a set of pairs of σ , radius of a CC circle, the motion leading up to each state, a kinematic graph, reference trajectory. In some embodiments, the memory 880 can include stored instructions implementing the method for the automated parking, wherein the instructions, when executed by the processor 870 carry out at least some steps of the method.

[0108]

Figure 9 shows a schematic of a system according to one embodiment. The system include a vehicle 901 including a processor 902 configured for performing an automated parking 950. The vehicle also includes at least one sensor, such as a LIDAR 910 and/or a camera 920. The LIDAR sensor 910 is the low-resolution first sensor and the camera 920 is the high-resolution second sensor. The sensor 910 and/or 920 is operatively connected to the processor 902 and is configured for sensing information indicative of the geometry of at least part of the parking space. Using this information, the processor 902 determines and/or updates the map of the parking space 900. To that end, the processor 902 performs the automated parking 950 using the map 900.

[0109]

Figure 10 shows a block diagram of an automated parking system 1000 according to some embodiments. The system 1000 can be implemented internal to the vehicle 901. Additionally, or alternatively, the system 1000 can be communicatively connected to the vehicle 901.

[0110]

The system 1000 can include one or combination of a camera 1010, an inertial measurement unit (IMU) 1030, a processor 1050, a memory 1060, a transceiver 1070, and a display/screen 1080, which can be operatively coupled to other components through connections 1020. The connections 1020 can comprise buses, lines, fibers, links or combination thereof.

[0111]

The transceiver 1070 can, for example, include a transmitter enabled to transmit one or more signals over one or more types of wireless communication networks and a receiver to receive one or more signals transmitted over the one or more types of wireless communication networks. The transceiver 1070 can permit communication with wireless networks based on a variety of technologies such as, but not limited to, femtocells, Wi-Fi networks or Wireless Local Area Networks (WLANs), which may be based on the IEEE 802.11 family of standards, Wireless Personal Area Networks (WPANS) such Bluetooth, Near Field Communication (NFC), networks based on the IEEE 802.15x family of standards, and/or Wireless Wide Area Networks (WWANs) such as LTE, WiMAX, etc. The system 400 can also include one or more ports for communicating over wired networks.

[0112]

In some embodiments, the system 1000 can comprise image sensors such as CCD or CMOS sensors, lasers and/or camera 1010, which are hereinafter referred to as "sensor 1010". For example, the sensor 1010 can convert an optical image into an electronic or digital image and can send acquired images to processor 1050. Additionally, or alternatively, the sensor 1010 can sense the light reflected from a target object in a scene and submit the intensities of the captured light to the processor 1050.

[0113]

For example, the sensor 1010 can include color or grayscale cameras, which provide "color information." The term "color information" as used herein refers to color and/or grayscale information. In general, as used herein, a color image or color information can be viewed as comprising 1 to N channels, where

N is some integer dependent on the color space being used to store the image. For example, an RGB image comprises three channels, with one channel each for Red,

Blue and Green information.

[0114]

For example, the sensor 1010 can include a depth sensor for providing "depth information." The depth information can be acquired in a variety of ways using depth sensors. The term "depth sensor" is used to refer to functional units that may be used to obtain depth information independently and/or in conjunction with some other cameras. For example, in some embodiments, the depth sensor and the optical camera can be part of the sensor 1010. For example, in some embodiments, the sensor 1010 includes RGBD cameras, which may capture per-pixel depth (D) information when the depth sensor is enabled, in addition to color (RGB) images.

[0115]

As another example, in some embodiments, the sensor 1010 can include a 3D Time Of Flight (3DTOF) camera. In embodiments with 3DTOF camera, the depth sensor can take the form of a strobe light coupled to the 3DTOF camera, which can illuminate objects in a scene and reflected light can be captured by a CCD/CMOS sensor in the sensor 1010. Depth information can be obtained by measuring the time that the light pulses take to travel to the objects and back to the sensor.

[0116]

As a further example, the depth sensor can take the form of a light source coupled to the sensor 1010. In one embodiment, the light source projects a structured or textured light pattern, which can include one or more narrow bands of light, onto objects in a scene. Depth information is obtained by exploiting geometrical distortions of the projected pattern caused by the surface shape of the object. One embodiment determines depth information from stereo sensors such as a combination of an infra-red structured light projector and an infra-red camera registered to a RGB camera. [0117]

In some embodiments, the sensor 1010 includes stereoscopic cameras. For example, a depth sensor may form part of a passive stereo vision sensor, which may use two or more cameras to obtain depth information for a scene. The pixel coordinates of points common to both cameras in a captured scene may be used along with camera pose information and/or triangulation techniques to obtain per-pixel depth information.

[0118]

In some embodiments, the system 1000 can be operatively connected to multiple sensors 1010, such as dual front cameras and/or a front and rear-facing cameras, which may also incorporate various sensors. In some embodiments, the sensors 1010 can capture both still and video images. In some embodiments, the sensor 1010 can include RGBD or stereoscopic video cameras capable of capturing images at, e.g., 30 frames per second (fps). In one embodiment, images captured by the sensor 1010 can be in a raw uncompressed format and can be compressed prior to being processed and/or stored in memory 1060. In some embodiments, image compression can be performed by the processor 1050 using lossless or lossy compression techniques.

[0119]

In some embodiments, the processor 1050 can also receive input from IMU 1030. In other embodiments, the IMU 1030 can comprise 3-axis accelerometer(s), 3-axis gyroscope(s), and/or magnetometer(s). The IMU 1030 can provide velocity, orientation, and/or other position related information to the processor 1050. In some embodiments, the IMU 1030 can output measured information in

synchronization with the capture of each image frame by the sensor 1010. In some embodiments, the output of the IMU 1030 is used in part by the processor 1050 to fuse the sensor measurements and/or to further process the fused

measurements.

[0120] The system 1000 can also include a screen or display 1080 rendering images, such as color and/or depth images. In some embodiments, the display 1080 can be used to display live images captured by the sensor 1010, fused images, augmented reality (AR) images, graphical user interfaces (GUIs), and other program outputs. In some embodiments, the display 1080 can include and/or be housed with a touchscreen to permit users to input data via some combination of virtual keyboards, icons, menus, or other GUIs, user gestures and/or input devices such as styli and other writing implements. In some

embodiments, the display 1080 can be implemented using a liquid crystal display (LCD) display or a light emitting diode (LED) display, such as an organic LED (OLED) display. In other embodiments, the display 480 can be a wearable display. In some embodiments, the result of the fusion can be rendered on the display 1080 or submitted to different applications that can be internal or external to the system 1000.

[0121]

Exemplary system 1000 can also be modified in various ways in a manner consistent with the disclosure, such as, by adding, combining, or omitting one or more of the functional blocks shown. For example, in some configurations, the system 1000 does not include the IMU 1030 or the transceiver 1070. Further, in certain example implementations, the system 1000 include a variety of other sensors (not shown) such as an ambient light sensor, microphones, acoustic sensors, ultrasonic sensors, laser range finders, etc. In some embodiments, portions of the system 400 take the form of one or more chipsets, and/or the like.

[0122]

The processor 1050 can be implemented using a combination of hardware, firmware, and software. The processor 1050 can represent one or more circuits configurable to perform at least a portion of a computing procedure or process related to sensor fusion and/or methods for further processing the fused

measurements. The processor 1050 retrieves instructions and/or data from memory 1060. The processor 1050 can be implemented using one or more application specific integrated circuits (ASICs), central and/or graphical processing units (CPUs and/or GPUs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), controllers, micro-controllers,

microprocessors, embedded processor cores, electronic devices, other electronic units designed to perform the functions described herein, or a combination thereof.

[0123]

The memory 1060 can be implemented within the processor 1050 and/or external to the processor 1050. As used herein the term "memory" refers to any type of long term, short term, volatile, nonvolatile, or other memory and is not to be limited to any particular type of memory or number of memories, or type of physical media upon which memory is stored. In some embodiments, the memory 1060 holds program codes that facilitate the automated parking.

[0124]

For example, the memory 1060 can store the measurements of the sensors, such as still images, depth information, video frames, program results, as well as data provided by the IMU 1030 and other sensors. The memory 1060 can store a memory storing a geometry of the vehicle, a map of the parking space, a kinematic model of the vehicle, and a dynamic model of the vehicle. In general, the memory 1060 can represent any data storage mechanism. The memory 1060 can include, for example, a primary memory and/or a secondary memory. The primary memory can include, for example, a random access memory, read only memory, etc. While illustrated in Figure 4 as being separate from the processors 1050, it should be understood that all or part of a primary memory can be provided within or otherwise co-located and/or coupled to the processors 1050.

[0125]

Secondary memory can include, for example, the same or similar type of memory as primary memory and/or one or more data storage devices or systems, such as, for example, flash/USB memory drives, memory card drives, disk drives, optical disc drives, tape drives, solid state drives, hybrid drives etc. In certain implementations, secondary memory can be operatively receptive of, or otherwise configurable to a non-transitory computer-readable medium in a removable media drive (not shown). In some embodiments, the non-transitory computer readable medium forms part of the memory 1060 and/or the processor 1050.

[0126]

The above-described embodiments of the present invention can be implemented in any of numerous ways. For example, the embodiments may be implemented using hardware, software or a combination thereof. When implemented in software, the software code can be executed on any suitable processor or collection of processors, whether provided in a single computer or distributed among multiple computers. Such processors may be implemented as integrated circuits, with one or more processors in an integrated circuit

component. Though, a processor may be implemented using circuitry in any suitable format.

[0127]

Also, the embodiments of the invention may be embodied as a method, of which an example has been provided. The acts performed as part of the method may be ordered in any suitable way. Accordingly, embodiments may be constructed in which acts are performed in an order different than illustrated, which may include performing some acts simultaneously, even though shown as sequential acts in illustrative embodiments.

[0128]

Use of ordinal terms such as "first," "second," in the claims to modify a claim element does not by itself connote any priority, precedence, or order of one claim element over another or the temporal order in which acts of a method are performed, but are used merely as labels to distinguish one claim element having a certain name from another element having a same name (but for use of the ordinal term) to distinguish the claim elements.