
II. RELATED WORK
The current state-of-the-art formulation of RA-SLAM is
through nonlinear least-squares (NLS) optimization based
upon MAP inference [12]. This work specifically focuses on
initialization for this MAP formulation. Though we focus
on MAP estimation, there are many important previous
formulations of RA-SLAM using: extended Kalman filters
[3], [13]–[15], particle filters [16], [17], mixture models [18],
and NLS optimization [1], [2], [19].
We first survey previous initialization strategies for RA-
SLAM problems and then expand the scope to cover initial-
ization strategies for general pose-graph SLAM problems.
As our initialization approach is a convex relaxation of the
RA-SLAM problem, we discuss related convex relaxations
in robotic state-estimation. Finally, as the key challenge
that differentiates RA-SLAM from pose-graph SLAM is
the nature of range measurements, we discuss the closely
related field of sensor network localization (SNL) and convex
relaxations developed specifically for SNL problems.
A. Initialization in RA-SLAM
Notable work in the single-robot RA-SLAM case used
spectral graph clustering to initially estimate beacon loca-
tions [20]. Similarly, [21] used spectral decomposition to
jointly estimate robot poses and beacon locations. However,
these approaches [20], [21] do not account for sensor noise
models and do not readily extend to multi-robot scenar-
ios. Other initializations for multi-robot RA-SLAM used
coordinated movement patterns to find the (single) relative
transform between robots [22], [23].
Finally, a series of works developed increasingly sophis-
ticated methods to compute relative transforms between two
robots using range and odometry measurements [24]–[30].
While these approaches represent great progress, they are
only capable of solving for a single inter-robot relative
transform and rely on (noisy) odometric pose composition
to initialize all other poses. In fact, [28]–[30] utilize convex
relaxation to solve the relative transformation problem.
B. Initialization and Convex Relaxation in State-Estimation
In pose-graph SLAM many existing initializations solve
approximations of the MAP problem. [7] leveraged the
relationships between rotation and translation measurements
in SLAM and explored the use of several rotation averaging
algorithms [31]–[34] to obtain initial estimates. The same
relationships were used to approximate the SLAM problem
as two sequential linear estimation problems [35]. However,
these approaches do not consider range measurements and
thus do not extend to RA-SLAM problems.
With exception to [34], these initialization procedures rep-
resent convex relaxations of the pose-graph SLAM problem.
Subsequent works in pose-graph SLAM developed convex
relaxations which obtained exact solutions to the non-convex
MAP problem [8], [11], [36]–[38]. Other works developed
exact convex relaxations for the problems of extrinsic cali-
bration [9], two-view relative pose estimation [39], [40], and
spline-based trajectory estimation from range measurements
with known beacon locations [41]. With exception to [41],
these convex relaxations also do not extend to range mea-
surements. However, [41] required known beacon locations,
allowed only for range measurements, and only considered
single-robot localization. SCORE allows for multi-robot RA-
SLAM, and combines measurements of rigid transformations
and ranges without necessitating any information a priori.
A number of convex relaxations exist for sensor network
localization (SNL). As SNL centers around point-to-point
distance measurements, these works are closely related to
SCORE. Previous works developed and analyzed semidef-
inite program (SDP) [10], [42], [43] and SOCP [44]–[47]
relaxations of the SNL problem. Whereas these relaxations
only consider range measurements and can only estimate
the Euclidean position of points, SCORE allows for rigid
transformation measurements and estimation of poses. Ad-
ditionally, as SCORE is a SOCP, it is substantially more
scalable than these SDP relaxations.
C. Novelty of Our Approach
SCORE is similar to [31] in its relaxation of rotations and
to [45] in its relaxation of range measurements. However,
it differs from these works in that it jointly considers both
relative pose and range measurements. Notably, SCORE
generalizes the chordal-initialization of [31] and the SOCP
relaxation of [45].
Additionally, there is much commonality to [28]–[30]
in that the convex relaxation of SCORE considers both
range and relative transformation measurements. Unlike
these works, SCORE generalizes to arbitrary dimensions
(e.g. 2- or 3-D) and multi-robot cases. Furthermore, SCORE
differs in that the convex relaxation is derived from the full
MAP problem, and thus utilizes all measurements and jointly
solves for an initialization of all RA-SLAM variables.
III. NOTATION AND MATHEMATICAL PRELIMINARIES
Symbols: This work generalizes across dimensionalities
(e.g., 2-D or 3-D). As such the mathematical presentation
will assume to be working in d-dimensional space. We
denote the d-dimensional identity matrix as Idand the special
orthogonal group as SO(d). We refer to the set of real values
as Rand the set of nonnegative reals as R+. Furthermore,
we indicate vectors and scalars with lowercase symbols (e.g.,
t) and matrices with uppercase symbols (e.g., R). Noisy
measurements are indicated with a tilde (e.g., ˜
R). The true
(latent) value of a quantity will be underlined (e.g., R)
Mathematical Operators: We denote the matrix trace as
tr(·). A diagonal matrix with diagonal elements [a1, . . . , an]
is indicated as diag([a1, . . . , an]). The vector 2-norm is
indicated as k·k2and the Frobenius norm is indicated
as k · kF. We note the following relationship between the
Frobenius norm and trace operators: kAk2
F= tr(AA>) =
tr(A>A).
Probability Distributions: We indicate a Gaussian prob-
ability distribution over Euclidean space, with mean µ∈
Rdand variance Σ∈Rd×das N(µ, Σ). Similarly, we