<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Indoor Positioning based on Active Radar Sensing and Passive Reflectors: Concepts &amp; Initial Results</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Pascal Schlachter</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Zhibin Yu</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Naveed Iqbal</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Xiaofeng Wu</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Sven Hinderer</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Bin Yang</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Institute of Signal Processing and System Theory, University of Stuttgart</institution>
          ,
          <addr-line>Pfafenwaldring 47, 70569 Stuttgart</addr-line>
          ,
          <country country="DE">Germany</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Munich Research Center, Huawei Technologies Duesseldorf GmbH</institution>
          ,
          <addr-line>Riesstraße 25, 80992 Munich</addr-line>
          ,
          <country country="DE">Germany</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>To navigate reliably in indoor environments, an industrial autonomous vehicle must know its position. However, current indoor vehicle positioning technologies either lack accuracy, usability or are too expensive. Thus, we propose a novel concept called local reference point assisted active radar positioning, which is able to overcome these drawbacks. It is based on distributing passive retroreflectors in the indoor environment such that each position of the vehicle can be identified by a unique reflection characteristic regarding the reflectors. To observe these characteristics, the autonomous vehicle is equipped with an active radar system. On one hand, this paper presents the basic idea and concept of our new approach towards indoor vehicle positioning and especially focuses on the crucial placement of the reflectors. On the other hand, it also provides a proof of concept by conducting a full system simulation including the placement of the local reference points, the radar-based distance estimation and the comparison of two diferent positioning methods. It successfully demonstrates the feasibility of our proposed approach.</p>
      </abstract>
      <kwd-group>
        <kwd>Indoor positioning</kwd>
        <kwd>active radar sensing</kwd>
        <kwd>passive local reference point</kwd>
        <kwd>autonomous vehicle system</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>
        Introduced in 1953, automated guided vehicles (AGVs) have been used extensively in industrial
environments during the past decades. Recently, they are being replaced by autonomous
mobile robots (AMRs). In contrast to their predecessors, AMRs are able to independently move
through their environment without relying on tracks, predefined paths or supervision by human
operators. In this way, they can increase flexibility, diversify applications and simplify the way
of life of human beings within a broad field of diferent branches [
        <xref ref-type="bibr" rid="ref1 ref2">1, 2</xref>
        ]. In order to reliably
navigate towards their destination and to safely avoid collisions with known obstacles, walls
or other AMRs, they must determine their position. Hence, a dependable and accurate indoor
vehicle positioning system is required. Therefore, in this paper we propose local reference point
(LRP) assisted active radar positioning. It is a novel concept towards indoor positioning which
is specifically tailored to the application on autonomous vehicles.
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. Related Work</title>
      <p>All common approaches towards indoor vehicle positioning can be classified into three
categories, namely fingerprint-based, geometry-based or SLAM-based methods.</p>
      <p>
        First, fingerprint-based methods [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ] try to identify each position in an indoor environment
by a unique set of significant features (a so-called fingerprint) of a wireless transmission
channel. Such methods can be built upon existing communication systems (e.g. Wi-Fi) by using
features that are already available like the channel state information (CSI) or the received signal
strength indicator (RSSI). However, a large dataset of fingerprints must be collected by taking
measurements on a dense grid. Moreover, this dataset must be updated each time anything in
the environment changes.
      </p>
      <p>
        Second, geometry-based methods determine the distances and/or angles to multiple active
reference points and use the result to perform positioning via trilateration or triangulation.
Most commonly, this is either done using the Bluetooth Low Energy (BLE) standard [
        <xref ref-type="bibr" rid="ref4 ref5">4, 5</xref>
        ] or
the ultra-wideband (UWB) technology [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. While BLE-based systems rely on measurements
of the RSSI and/or the angles of arrival, UWB systems mostly use time-based measurements.
Nevertheless, also geometry-based methods sufer from drawbacks for the application on AMRs
such as a lack of accuracy (BLE) or the need of costly hardware and challenging calibration
(UWB).
      </p>
      <p>
        Third, simultaneous localization and mapping (SLAM) [
        <xref ref-type="bibr" rid="ref7 ref8">7, 8</xref>
        ] can be applied. Its idea is to
incrementally build a map of the unknown environment and use it for positioning at the same
time. Thereby, the vehicle only relies on its onboard sensors scanning the environment and
does not need an external reference system. However, like fingerprint-based methods, also
SLAM-based methods struggle with non-stationary and large environments. Moreover, they
can only estimate the relative position w.r.t. the simultaneously learned map.
      </p>
      <p>
        Finally, sensor fusion, meaning the combination of multiple diferent positioning methods, can
be applied in order to increase the accuracy and robustness. Especially, an iterative combination
of the positioning results with the observations of an inertial measurement unit (IMU) is
considered promising [
        <xref ref-type="bibr" rid="ref10 ref9">9, 10</xref>
        ]. Thereby, the IMU captures the odometry data of the vehicle,
meaning its change of position over time. Nevertheless, an IMU is not suited as a standalone
positioning technique because it can only provide a relative position and its results typically
drift away from the ground truth over time.
      </p>
      <p>In conclusion, there is a lack of low-cost and easy to use but still highly accurate indoor
positioning techniques. This especially holds for industrial applications like AMRs, since they
typically require highly accurate estimates of the absolute position while operating in a regularly
changing environment.</p>
    </sec>
    <sec id="sec-3">
      <title>3. Basic Idea</title>
      <p>The basic idea of our novel approach towards indoor positioning is to equip the environment with
multiple low-cost passive retroreflectors such that a radar system mounted on an AMR observes
a unique reflection characteristic at each position. In general, this reflection characteristic can
include the distances, angles and relative velocities between the vehicle and the reflectors as
well as the radar cross sections (RCSs) and polarimetrics of the corresponding reflectors. In the
remainder of this paper, the reflectors will be denoted as local reference points (LRPs) and we
will refer to the observed reflection characteristics as the LRP fingerprints.</p>
      <p>Accordingly, our method combines the ideas of geometry- and fingerprint-based methods
to get the best of both worlds. By using a well-designed radar system, our method is not only
able to provide robust and highly accurate distance measurements to the LRPs compared to
RSSI based estimates but also only requires cheap and simple passive instead of active reference
points. Moreover, in contrast to existing fingerprint-based methods, the LRP fingerprint of
each position can be calculated from the well-known positions and properties of the LRPs and
does not change as long as the positions of the LRPs remain constant. In this way, our method
does not need the extensive collection of a dataset and is also applicable in regularly changing
environments.</p>
      <p>In this initial work, we do not want to limit our approach to only work with high-end radar
technologies. Instead, we propose and investigate a baseline which only requires a simple
lowcost radar system with a single non-polarimetric TX and RX antenna, respectively. Showing
that it can already provide reliable position estimates allows us to build onto this foundation and
further develop our method in future work e.g. by additionally using angular or polarimetric
information to further improve its accuracy.</p>
      <p>In the following, we describe the four subtasks of the iterative positioning and navigation of
an AMR in detail.</p>
      <sec id="sec-3-1">
        <title>3.1. LRP placement</title>
        <p>The first subtask is to place the LRPs into the environment. The resulting LRP layout is desired
to provide a unique LRP fingerprint for each position. In this work, we restrict ourselves to
the simplest setup, namely using only the set of distances to the LRPs as the LRP fingerprint
for each position. By using a more advanced radar system, also other measurements can be
incorporated into the LRP fingerprint.</p>
        <p>Note that, in contrast to existing geometry-based indoor positioning methods or GNSS, the
individual distances cannot be assigned to the respective LRP they belong to. Concretely, due to
the application of passive instead of active reference points, the only way to distinguish between
LRPs is to use diferent reflectors which difer in their RCSs. However, this is challenging in
practice because of the large multi-path fading efects within indoor environments. Accordingly,
we assume that at most two diferent types of LRPs can be distinguished. As a result, positioning
via trilateration is not possible.</p>
        <p>In the following section 4, we analyze which constraints an LRP layout must fulfill to ensure
a bijective mapping between the LRP fingerprints and the positions. Thereby, we must deal
with a so far unsolved mathematical ambiguity problem. From all layouts which provide unique
LRP fingerprints for all positions, we can then select the one which best fulfills secondary
requirements like having large distances among the LRPs or keeping a distance to all walls.</p>
      </sec>
      <sec id="sec-3-2">
        <title>3.2. Radar transmission</title>
        <p>If the environment is equipped with LRPs, the iterative positioning and navigation of an AMR
can start. The first step is the radar transmission. Although, we only assume a simple radar
hardware which cannot capture angular or polarimetric information, measuring the relative
velocity and RCS of targets besides their distances comes without additional cost. Both are
helpful, for instance, to distinguish LRPs from other targets as described in more detail below.
Moreover, the simple radar hardware does not limit the measurement accuracy. Instead, the
resolution and maximal detectable range of the distance and relative velocity are defined by the
choice of the radar system parameters. Thus, they must be designed carefully as they directly
influence the achievable positioning accuracy.</p>
        <p>We decide to apply a chirp sequence frequency modulated continuous wave (FMCW) radar
system. It is state-of-the-art and, in contrast to other radar technologies, enables a joint
estimation of the distances, relative velocities and RCSs of multiple targets at the same time.
Obviously, this is required here as we intend to detect multiple LRPs simultaneously.</p>
        <p>
          If multiple vehicles equipped with LRP assisted active radar positioning operate in the same
environment, their respective radar systems interfere with each other. This can lead to ghost
targets and complicates the detection of the LRPs. To mitigate the interference, well-known
interference cancellation and coordination techniques like multiplexing, listen-before-talk or
radar-specific approaches either based on classical signal processing [
          <xref ref-type="bibr" rid="ref11">11</xref>
          ] or deep learning [
          <xref ref-type="bibr" rid="ref12">12</xref>
          ]
can be applied.
        </p>
      </sec>
      <sec id="sec-3-3">
        <title>3.3. Radar signal processing</title>
        <p>The radar signal processing estimates the distances, relative velocities and RCSs for all LRPs
from the radar RX signal. Therefore, first the peaks within the range-Doppler spectrum are
detected using a constant false alarm rate (CFAR) detector. Thereby, the position of each peak
determines the distance and relative velocity of the corresponding target while the RCS is
obtained from its magnitude. Subsequently, the detected targets must be classified to determine
whether they correspond to an LRP or to other objects in the environment. Most easily, this
can be done by comparing the observed RCSs to the well-known properties of the reflectors.
Additionally, to increase robustness, we can track the distances of the targets over time by
applying a motion model. Concretely, we use the relative velocities to calculate the changes
of the distances assuming a straight movement with constant velocity. This assumption is
suficient if the time between two measurements is small compared to the velocity of the AMR.
In this way, we can predict the distance to each target for the next iteration based on the current
measurements. These predictions help to identify the LRPs within the target list. Moreover, by
applying a Kalman filter, the predictions are combined with the actual distance measurements to
smooth the results. In conclusion, tracking the targets over time increases the robustness against
measurement, detection and classification errors. The resulting estimated LRP fingerprints are
then used for positioning.</p>
      </sec>
      <sec id="sec-3-4">
        <title>3.4. Positioning</title>
        <p>In this work, we present two diferent approaches to find the position of the AMR based on the
set of smoothed distance measurements to the LRPs (i.e. the estimated LRP fingerprint) output
by the Kalman filter, namely using a look up table and AMCL. Both approaches are described in
detail below.</p>
        <sec id="sec-3-4-1">
          <title>3.4.1. Look up table</title>
          <p>Most straightforwardly, the positioning can be achieved by using a look up table. Since the LRP
layout is well-known, for each position the corresponding LRP fingerprint can be calculated
easily. Now the inverse mapping is desired, meaning getting the position based on the fingerprint.
Nevertheless, the latter can be obtained using the first mapping. Concretely, for each position
the distances to the LRPs are calculated and the resulting fingerprints are stored in a look up
table. Then, for an estimated distance set, we search for the corresponding entry in the look
up table and finally get the position. However, note that it is not possible to store all infinitely
precise positions into a look up table. Instead, the entries must be limited to a finite number.
For instance, this can be achieved by quantizing the distances to a certain resolution. Then we
get a finite number of grids which can be identified by their quantized LRP fingerprints. These
ifngerprints are then mapped to the center points of the corresponding grids to get a finite look
up table. Thereby, obviously, the resolution of the quantization directly determines the size of
the look up table. It must be chosen as a trade-of between the memory consumption and the
positioning accuracy.</p>
        </sec>
        <sec id="sec-3-4-2">
          <title>3.4.2. Adaptive Monte Carlo localization</title>
          <p>While the look up table approach is straightforward, it is sensitive to erroneous or ambiguous LRP
ifngerprint estimates. Concretely, if one or multiple elements of the estimated LRP fingerprint
are still inaccurate, the look up table directly leads to a wrong position estimate. To gain
robustness, the idea is to use the odometry data of the AMR, meaning its movement information
between two measurements, to additionally track the position over time. The odometry data
can be either taken from the motion control commands which are passed to the drive unit of
the AMR or measured by an IMU.</p>
          <p>
            The position tracking is done by applying a particle filter based approach called Adaptive
Monte Carlo Localization (AMCL) [
            <xref ref-type="bibr" rid="ref13 ref14">13, 14</xref>
            ]. It recursively estimates the three-dimensional state
vector  = [, ,  ] , i.e. the position and orientation of the AMR at time instance k.
Therefore, the posterior density (|) of the current state conditioned on all previous LRP
^
ifngerprint estimates  = { ,  = 1, . . . , } is recursively approximated by a weighted set
of  random particles  = {{, },  = 1, . . . , } at each iteration . Thereby, each
particle  describes one possible state vector and its corresponding weight  denotes the
probability that this particle describes the true current state vector .
          </p>
          <p>After drawing an initial population of 0 particles from a prior distribution (0), each
particle is weighted according to the likelihood  = (^|). Hence, the weight of each
^
particle  denotes the likelihood that the LRP fingerprint estimate  is observed given that
the AMR is located at this respective particle . Finally, the weights are normalized such
that they sum to one. Then they are used to do the positioning by calculating the weighted
mean of the particles. Subsequently, we do a resampling by duplicating particles with large
weights and discarding particles with small weights. This prevents a depletion of the population
of particles since over time many particles drift far enough for their weight to become too
small to contribute to the posterior. Furthermore, the number of particles  can be adapted
to the variance of the previous population of particles to achieve a good trade-of between
accuracy and computational complexity. Finally, the movement of the AMR between two radar
measurements is modeled by propagating the particles according to the odometry data. Thereby,
we account for possible inaccuracies by adding a Gaussian distributed motion uncertainty.</p>
          <p>The previously described propagation of the particles, update of the weights, position
estimation and resampling is done in each iteration of the positioning and navigation process as shown
in Algorithm 1. Thereby, the population of particles more and more converges towards the true
position of the AMR with each new iteration. After convergence, since the position estimate
can only range within the set of particles, the positioning will still be suficiently accurate even
in the case of completely erroneous LRP fingerprint estimates. Moreover, in this way possible
ambiguities, meaning diferent positions that observe the same LRP fingerprint, can be resolved.
This enables the application of simpler LRP layouts.</p>
        </sec>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Placement of the Local Reference Points</title>
      <p>Now, we focus on one of the most crucial steps, namely the placement of the LRPs. Thereby, for
a given indoor environment, the goal is to find a placement of  LRPs
ℛ = {(1, 1, 1, 1) , . . . , ( ,  ,  ,  ) }
(1)
such that each position  = (, ) of an AMR can be characterized by a unique set of distances
to the LRPs (LRP fingerprint) . In other words, the desired solution is an LRP layout ℛ such
that the mapping  ↦→  is bijective. Thereby, (, , ) denotes the position of the -th LRP
while  describes its type. To the best of our knowledge this corresponds to a novel and so
y</p>
      <p>LRP 3
LRP 1
d2</p>
      <p>LRP 2
far unsolved mathematical ambiguity problem, since (in contrast to existing geometry-based
methods or GNSS) we cannot uniquely assign the distances to their corresponding LRPs.</p>
      <p>For our initial study in this work, we restrict ourselves to small-scaled rooms and make the
following assumptions:
• At most two types of LRPs can be distinguished reliably. Hence,  ∈ {0, 1},  =
1, . . . ,  .
• Diferent LRPs of the same type cannot be distinguished.
• All LRPs are detectable from all positions of the room. This combines two assumptions.</p>
      <p>First, we do not consider shadowing, meaning that the direct view of the radar towards
an LRP can be blocked. Second, we only consider environments where at each position
all LRPs are within the maximal detectable distance of the radar system.
• All LRPs are mounted at the same height, meaning 1 = . . . =  = . Therefore it can
be neglected and the problem can be reduced to the -plane.</p>
      <p>Under these assumptions, we systematically investigated the following research questions
(RQs) 1:</p>
      <p>RQ1 Does an LRP layout exist which provides unique LRP fingerprints for all positions?
RQ2 If yes: How can we systematically find such a layout, i.e. which constraints must it
fulfill and how many LRPs are necessary at least?
Concretely, we started with one LRP and analyzed step by step how each additional LRP can
help to avoid ambiguities, meaning diferent AMR positions with the same LRP fingerprint.
Thereby, we found that in general two types of ambiguities are present. We call them systematic
and random ambiguities. In the following we summarize our results. For more insights and
detailed proofs, we refer to [16] due to the limited space in this paper.</p>
      <sec id="sec-4-1">
        <title>4.1. Systematic ambiguities</title>
        <p>Systematic ambiguities denote all ambiguities which we can describe systematically. The
foundation for all systematic ambiguities is the case of using two LRPs of one type. Thereby,
intuitively, for each position in the room three other positions exist which have the same LRP
ifngerprint. These ambiguous positions are symmetric according to two symmetry axes. The
ifrst one is the line through the two LRPs and the second one is orthogonal to the first axis
passing through the middle between the LRPs. In this way, the ambiguities are also point
symmetric to the intersection point of the two axes. The symmetry axes (blue dotted lines)
and one arbitrary set of ambiguities having the same LRP fingerprint {1, 2} (red crosses) are
visualized in Figure 1a for an exemplary LRP layout. We can now use this knowledge to look at
three and four LRPs of the same type as three and six combinations of two diferent LRPs each,
respectively.</p>
        <p>Adding a third LRP of the same type cannot resolve all these ambiguities. Concretely, for
each combination of two LRPs there is a line containing point symmetric ambiguities which
also have the same distance to the respective third LRP. As shown in Figure 1b for an exemplary
LRP layout, each of these ambiguity lines (red, green and blue lines) passes through the center
point between two LRPs (i.e. the symmetry point) and is orthogonal to the line connecting this
point to the respective third LRP.</p>
        <p>Finally, using a fourth LRP can resolve all systematic ambiguities. This is the case if for
all six combinations of two LRPs, the symmetry axes do not cross the symmetry point of the
respective complementary combination. Otherwise, as shown in Figure 1c, we again get (at
least) a line of point symmetric ambiguities (red dotted line) which passes the symmetry point
and is orthogonal to the corresponding symmetry axis.
1A further analysis of these RQs under relaxed assumptions is done in [15] based on the results of this paper.</p>
        <p>Introducing a second type of LRP does not change much. When using three LRPs, one
distance can now be assigned uniquely to the one LRP being of a diferent type. Nevertheless,
still the ambiguity line corresponding to the combination of the two LRPs of the same type is
left. When using four LRPs where two of each are of the same type, still no ambiguities are
present if the previously stated condition holds. However, now we only need to check this
condition for the two combinations of two LRPs of the same type.</p>
        <p>In summary, a layout with three or less LRPs always contains systematic ambiguities while
a layout with four LRPs does not contain any systematic ambiguities as long as for all
combinations of two LRPs, the symmetry axes do not cross the symmetry point of the respective
complementary combination.</p>
      </sec>
      <sec id="sec-4-2">
        <title>4.2. Random ambiguities</title>
        <p>For layouts with three or four LRPs of one type, we find that there are additional ambiguities
which we cannot describe systematically. We denote them as random ambiguities. So far, when
searching for systematic ambiguities, for each combination of two LRPs and each position we
only considered the three ambiguous positions w.r.t. this exact combination of two LRPs and
checked under which circumstances they are also ambiguous w.r.t. the remaining LRP(s). E.g. in
Figure 1d, when considering the combination of the LRPs 1 and 2 and searching for ambiguities
of the arbitrary position indicated by a red cross, we only considered the positions indicated by
blue crosses. However, there are additional ambiguity candidates having the exact same two
distances in common with the considered position but w.r.t. one of the other combinations of
two LRPs. Accordingly, in Figure 1d the green and purple crosses show the positions whose
LRP fingerprints also contain the two distances 1 and 2 but w.r.t. the combination of the LRPs
1 and 3 or 2 and 3, respectively. Since we cannot assign the distances to the LRPs when using
only one LRP type, these candidates may also result in ambiguities.</p>
        <p>In contrast to the systematic ambiguities, here introducing a second type of LRP makes a
diference. When using four LRPs where two of each are of the same type, the only ambiguity
can stem from common ambiguities of the two complementary combinations of two similar LRPs.
This corresponds to the case we already studied when searching for systematic ambiguities.</p>
        <p>As a result, there are no random ambiguities when using layouts with two diferent types
of LRPs. By combining this result with our findings for the systematic ambiguities we can
answer RQ1 and RQ2 as follows: A layout using four LRPs of two types provides unique LRP
ifngerprints for all positions if the symmetry axes of one combination of LRPs of the same type
do not cross the symmetry point of the other combination and vice versa. A detailed proof can
be found in [16].</p>
      </sec>
    </sec>
    <sec id="sec-5">
      <title>5. Proof of Concept</title>
      <p>5.1. Setup
Finally, we conduct a system simulation to provide a proof of concept for our novel approach
towards indoor positioning under idealized conditions of the environment and the LRPs.
Concretely, we estimate the positions of an arbitrary path of an AMR within an empty room of</p>
      <p>8
15 14 13 12
x
9
11
10
size 5 m × 5 m × 4 m. The path consists of 15 locations which are 0.5 m apart from each other,
respectively. Furthermore, we assume the AMR to move linearly between two such locations
with a constant velocity of 2 m/s. Hence, the position is estimated every 0.25 s.</p>
      <p>The room is equipped with four LRPs of one type which are all at the same height of 3 m. We
decide to use only one LRP type because distinguishing between diferent types of LRPs which
only difer in their RCSs is challenging and error-prone in practice. Nonetheless, we desire
a layout with as few ambiguities as possible. Accordingly, we generate 500 random layouts,
which fulfill the previously derived constraint to avoid systematic ambiguities and select the
one with the smallest amount of random ambiguities. We expect the few remaining ambiguities
to be resolvable by applying AMCL. The LRP layout and corresponding ambiguities together
with the exemplary path of the AMR are shown in Figure 2.</p>
      <p>The radar system parameters are designed such that it can measure distances up to
max = 19.125 m with a resolution of Δ = 0.075 m and relative velocities in the range
[− 5.6816 m/s, 5.6816 m/s) with a resolution of Δ = 0.3551 m/s. The carrier frequency is
selected to  = 60 GHz. We simulate the radar transmission at the 15 exemplary positions
by applying a ray tracer in a 3D model of the environment. Thereby, for the proof of concept,
we neglect noise and difused scattering efects while reflections and thus multipath fading
are taken into account in the simulation. For each position, the ray tracer outputs the channel
impulse response (CIR) and the Doppler shift for each channel tap, which we use to calculate
the RX signal and finally the range-Doppler spectrum.</p>
      <p>Subsequently, we find the peaks within the range-Doppler spectrum by applying a CFAR
detector. Among the detected targets, we identify the four that correspond to the LRPs by a
combination of comparing the estimated RCS of each target to the known value of the LRPs
and tracking the targets over time. Concretely, the tracking is done by applying a linear motion
model using the distance and relative velocity measurements of the previous position to predict
the distances to the LRPs for the current position. The estimated LRP fingerprint then consists
of the four most likely targets regarding both the RCS estimates and the distance predictions.
Look up table
AMCL
1
2
3
4
5
Fig. 3: Errors between the estimated and true positions for all samples of the exemplary path using the
look up table and AMCL, respectively. Thereby, for AMCL the mean and standard deviation of ten runs
is reported. In contrast, the look up table approach is deterministic.</p>
      <p>Finally, using the estimated LRP fingerprints, we compare both the look up table and the
AMCL based position estimates with the true positions. Thereby, we choose the quantization
resolution of the look up table entries equally to the distance resolution of the radar system,
namely to 0.075 m, resulting in 6574 entries. Each estimated LRP fingerprint is then mapped
to the entry with the smallest Euclidean distance. Regarding AMCL, we initially start with
0 = 10000 uniform distributed particles. Concerning the propagation of the particles, we
add Gaussian distributed motion uncertainties to both the turned angle and traveled distance
contained in the odometry data. They have zero mean and the standard deviations   = 5∘
and   = 5 cm, respectively. For weighting the particles, (^|) is modeled as a Gaussian
likelihood whose mean is the calculated LRP fingerprint corresponding to the respective particle
. Moreover, its covariance matrix is chosen to  = (0.075 m)2, where  denotes the identity
matrix. Finally, to estimate the position of the AMR, we calculate the weighted mean of the 20
best particles that are within a range of 15 cm to the particle with the largest weight.</p>
      <sec id="sec-5-1">
        <title>5.2. Results</title>
        <p>The errors between the position estimates and the true positions for both the look up table
approach and AMCL are shown in Figure 3. Thereby, the look up table provides deterministic
position estimates while the results of the AMCL are subject to random deviations. Therefore,
AMCL is performed ten times and we report the mean and standard deviation for each position.</p>
        <p>The results prove two points. First and foremost, the system simulation successfully
demonstrates the feasibility of our new approach towards indoor vehicle positioning and therefore
provides a proof of concept. When using AMCL for positioning, the average of the mean
errors for the 15 exemplary positions is at about 0.097 m only marginally larger than the range
resolution of the radar system. Hence, our novel method proves to be able to reliably provide
accurate position estimates.</p>
        <p>Second, as expected, we observe AMCL to be significantly more robust than the look up
table. As long as the estimated LRP fingerprints are correct and unambiguous, both positioning
approaches achieve comparable accuracies. However, if they are erroneous or ambiguous
the look up table approach fails completely while AMCL can still provide accurate position
estimates. Concretely, the errors for the ninth and twelfth position using the look up table are
around 3.5 m and 1 m, respectively, while AMCL leads to errors of only around 0.1 m for both
positions. Nevertheless, the full truth also includes the fact that this finding only holds once the
whole set of particles  has converged tightly enough towards the true position of the AMR.
This typically takes a few iterations (here around six) which explains why for the positions
two and four both positioning techniques perform equally poorly. Accordingly, in practice
we recommend to initialize 0 already as accurate as possible or to perform some “warm-up”
iterations.</p>
      </sec>
    </sec>
    <sec id="sec-6">
      <title>6. Discussion about practical aspects</title>
      <sec id="sec-6-1">
        <title>6.1. Robustness against multipath fading</title>
        <p>In indoor environments, the efect of multipath fading is generally stronger than outdoors since
they are closed and narrow. While fingerprinting takes advantage of this, it is a challenge for
most other radio transmission based indoor positioning methods in practice. Also the radar
transmission of our proposed approach needs to cope with it. Accordingly, in the following, we
will discuss the robustness of our method against multipath fading.</p>
        <p>In terms of our application, multipath fading describes the efect that there are not only the
desired direct paths from the radar TX antenna to the reflectors and back to the RX antenna.
Instead, there are also countless other disturbing paths from the TX to the RX antenna like
reflections from objects, the walls, the ceiling or even multi-reflections. First, we take a closer
look at the direct paths. Their received power strictly follows the radar equation given as
(2)
(3)
 =  · (4 )342 ·  ,
where  and  are the transmitted and received powers, respectively,  and  denote
the antenna gains of the TX and RX antennas and   is the carrier wavelength of the radar
system. Moreover,  describes the distance between the radar and the reflector and  denotes
the RCS of the reflector. Thereby, note that  ,  ,  and   are all constant parameters of
the radar system. Other single-reflection paths also follow this equation but of course the RCS
 then is defined by the corresponding reflection surface/object while  describes the distance
to it. Multi-reflection paths, however, do not follow equation 2. Instead, they obviously need to
consider multiple RCSs, one for each reflection, and also difer in their path loss which depends
on the length of the path. Accordingly, the corresponding received power of a path experiencing
 reflections is calculated as</p>
        <p>2  1  2
 =  · (4 )22 · 4 21 · 4 22 · · ·
0</p>
        <p>4 2 ,
where 0 is the length of the path between the TX antenna and the first reflection, 1, . . . , − 1
denote the lengths of the paths between the individual reflections and  is the distance from
the final reflection to the RX antenna. Furthermore,  1, . . . ,   are the RCSs corresponding to
the  reflections.</p>
        <p>In order for our method to be robust against multipath fading, we face two challenges. First,
the received powers of the desired direct LRP-paths must stand out enough compared to other
paths with a similar travel time and to the noise floor such that they are detected as targets by
the CFAR detector. Second, among all the paths that are detected as targets we must find the
ones corresponding to the LRPs.</p>
        <p>The keys to the robustness against multipath fading of our proposed radar-based approach,
solving these two challenges, lie in the RCS of the reflectors and the RCS based LRP identification.
Using appropriate and carefully designed reflectors with the largest possible RCS solves the
ifrst challenge. In this case, we expect the corresponding direct LRP-paths to be much stronger
than reflections of random surfaces/objects and even more compared to multi-reflections.
Accordingly, the following section 6.2 provides a discussion about the reflector choice. A further
remedy can be to focus the radar transmission to a cone facing upwards. In this way, the number
of disturbing paths is reduced by, among others, eliminating the possibly strong horizontal
single-reflection paths. The second challenge is solved by the RCS based LRP identification
together with the Kalman filter based LRP tracking as already introduced in section 3.3. Since
an FMCW radar measures the distance  of each target based on the frequency shift, their RCSs
can be calculated from the observed received power  of the corresponding path according to
equation 2. These estimations can then be compared to the well known RCS of the reflectors to
identify the targets corresponding to the LRPs. Even in the unlikely cases that either another
path which does not correspond to an LRP has the same estimated RCS as the reflectors or LRP
targets are not detected by the CFAR detector, the Kalman filter based LRP tracking will help to
still identify/recover the desired LRP targets.</p>
        <p>The previously presented proof of concept simulation already showed that in this way our
method is indeed robust against multipath fading. Thereby, the applied ray tracer simulated all
paths between an omni-directional TX and the RX antenna within the simulation environment.
Even in this worst-case-scenario, meaning a completely unfocused radar transmission in a
narrow environment, the targets corresponding to the LRPs could be successfully detected and
identified for nearly all positions. This can be seen from the good positioning accuracies for
most positions using the look up table shown in Figure 3, since an erroneous LRP detection or
identification directly leads to large positioning errors when using the look up table. Moreover,
as already discussed in the previous section 5.2, the simulations also showed that AMCL can
still provide accurate position estimates for wrong LRP fingerprint estimates and in this way
even further improves the robustness against the efects of multipath fading.</p>
      </sec>
      <sec id="sec-6-2">
        <title>6.2. Reflector choice</title>
        <p>In order to allow a reliable detection of the LRPs in a wide area, the retroreflectors are desired
to have a large RCS and incidence angle. A first intuitive choice would be a spherical reflector
due to its unlimited incidence angle. Its RCS is given as
 spherical =  ·
︂( 0 )︂ 2
2
,
(4)
where 0 denotes its diameter [17]. Accordingly, for a suficient RCS large reflectors are required
which is undesirable in practice (e.g. for  spherical = 0 dbm2 already a diameter of 0 = 1.13 m
is required).</p>
        <p>Alternatively, trihedral corner reflectors are simple yet efective and therefore widely used.
Their RCS is given as
 corner =
4 4
3 2 ,
(5)
where  is the edge length of the trihedral and   denotes the carrier bandwidth of the radar
system. Consequently, already a small reflector can provide a large RCS (e.g. assuming   =
5 mm like in the simulation setup, then for  corner = 0 dbm2 only a edge length of  = 5 cm is
required). Although their incidence angle is limited they work well for angles within almost
± 40∘ [18]. Furthermore, multiple such reflectors can be assembled to further increase the
covered area. Nonetheless, [15] shows that our novel approach towards indoor positioning still
works well with LRPs having limited incidence angles by further analyzing the LRP positioning
problem under this constraint based on the results of this paper.</p>
        <p>Besides trihedral corner reflectors, promising choices combining suficiently large RCSs and
incidence angles are the Luneburg lens and the Van-Atta array [19]. In further studies, we
will experimentally evaluate which of the introduced retroreflectors is the best choice for our
application.</p>
      </sec>
    </sec>
    <sec id="sec-7">
      <title>7. Conclusion and Outlook</title>
      <p>In conclusion, we proposed a novel approach towards indoor vehicle positioning. Its basic idea
is to identify each position by the unique reflection characteristic (the so-called LRP fingerprint)
an active radar system on top of an autonomous vehicle observes regarding a predefined layout
of passive LRPs. To avoid having to make limiting assumptions regarding the radar system
and to ensure robustness for further developments, we focused in this work on the set of
distances to the LRPs as LRP fingerprints and assumed that at most two types of LRPs can be
distinguished. Since this restriction does not allow the use of trilateration, we studied how
to find an LRP layout which provides a unique LRP fingerprint for each position. Thereby,
we found that a layout using four LRPs does not provide any systematic ambiguities if for all
combinations of two LRPs, the symmetry axes do not cross the symmetry point of the respective
complementary combination. Moreover, we showed that random ambiguities can be prevented
by using two types of LRPs instead of only one. Finally, we provided a successful proof of
concept by conducting a system simulation which demonstrated the feasibility of our proposed
approach.</p>
      <p>Further studies will be conducted in diferent directions. The first is the mitigation of the
assumption that all LRPs are visible from all positions. This especially includes the extension to
larger rooms, the application of reflectors with limited incident angles and the investigation
of shadowing (some of this is examined in [15]). The second is to investigate the potential
of additionally using angular and/or polarimetric informations. Are their benefits worth the
additional costs? Finally, a hardware implementation can be helpful to enable experiments in a
real-world scenario.
2004. URL: http://www.cim.mcgill.ca/~yiannis/particletutorial.pdf.
[15] S. Hinderer, P. Schlachter, Z. Yu, X. Wu, B. Yang, Indoor positioning based on active radar
sensing and passive reflectors: Reflector placement optimization (2023).
[16] P. Schlachter, Indoor Positioning based on Active Radar Sensing and Passive Reference</p>
      <p>Points, Master’s thesis, University of Stuttgart, Stuttgart, Germany, 2021.
[17] P. Rajyalakshmi, G. Raju, Characteristics of radar cross section with diferent objects,</p>
      <p>International Journal of Electronics and Communication Engineering 4 (2011) 205–216.
[18] C. Li, J. Zhao, J. Yin, G. Zhang, X. Shan, Analysis of rcs characteristic of dihedral corner and
triangular trihedral corner reflectors, in: 2010 5th International Conference on Computer
Science &amp; Education, 2010, pp. 40–43. doi:10.1109/ICCSE.2010.5593647.
[19] D. Bird, Design and manufacture of a low-profile radar retro-reflector, in: NATO RTO SCI
Symposium on Sensors and Sensor Denial by Camouflage, Concealment and Deception
(RTO-MP-SCI-145). Brussels, Belgium, 2004.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>E. A.</given-names>
            <surname>Oyekanlu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. C.</given-names>
            <surname>Smith</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W. P.</given-names>
            <surname>Thomas</surname>
          </string-name>
          ,
          <string-name>
            <given-names>G.</given-names>
            <surname>Mulroy</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Hitesh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Ramsey</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D. J.</given-names>
            <surname>Kuhn</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. D.</given-names>
            <surname>Mcghinnis</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S. C.</given-names>
            <surname>Buonavita</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N. A.</given-names>
            <surname>Looper</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Ng</surname>
          </string-name>
          ,
          <string-name>
            <surname>A</surname>
          </string-name>
          . Ng'oma, W. Liu,
          <string-name>
            <given-names>P. G.</given-names>
            <surname>Mcbride</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M. G.</given-names>
            <surname>Shultz</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>Cerasi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Sun</surname>
          </string-name>
          ,
          <article-title>A review of recent advances in automated guided vehicle technologies: Integration challenges and research areas for 5g-based smart manufacturing applications</article-title>
          ,
          <source>IEEE Access 8</source>
          (
          <year>2020</year>
          )
          <fpage>202312</fpage>
          -
          <lpage>202353</lpage>
          . doi:
          <volume>10</volume>
          .1109/ACCESS.
          <year>2020</year>
          .
          <volume>3035729</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>M. B.</given-names>
            <surname>Alatise</surname>
          </string-name>
          ,
          <string-name>
            <given-names>G. P.</given-names>
            <surname>Hancke</surname>
          </string-name>
          ,
          <article-title>A review on challenges of autonomous mobile robot and sensor fusion methods</article-title>
          ,
          <source>IEEE Access 8</source>
          (
          <year>2020</year>
          )
          <fpage>39830</fpage>
          -
          <lpage>39846</lpage>
          . doi:
          <volume>10</volume>
          .1109/ACCESS.
          <year>2020</year>
          .
          <volume>2975643</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>A.</given-names>
            <surname>Khalajmehrabadi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.</given-names>
            <surname>Gatsis</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Akopian</surname>
          </string-name>
          ,
          <article-title>Modern wlan fingerprinting indoor positioning methods and deployment challenges</article-title>
          ,
          <source>IEEE Communications Surveys &amp; Tutorials</source>
          <volume>19</volume>
          (
          <year>2017</year>
          )
          <fpage>1974</fpage>
          -
          <lpage>2002</lpage>
          . doi:
          <volume>10</volume>
          .1109/COMST.
          <year>2017</year>
          .
          <volume>2671454</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>Y.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>X.</given-names>
            <surname>Yang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Zhao</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Liu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Cuthbert</surname>
          </string-name>
          ,
          <article-title>Bluetooth positioning using rssi and triangulation methods</article-title>
          ,
          <source>in: 2013 IEEE 10th Consumer Communications and Networking Conference (CCNC)</source>
          ,
          <year>2013</year>
          , pp.
          <fpage>837</fpage>
          -
          <lpage>842</lpage>
          . doi:
          <volume>10</volume>
          .1109/CCNC.
          <year>2013</year>
          .
          <volume>6488558</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>G.</given-names>
            <surname>Pau</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Arena</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y. E.</given-names>
            <surname>Gebremariam</surname>
          </string-name>
          ,
          <string-name>
            <surname>I. You</surname>
          </string-name>
          ,
          <article-title>Bluetooth 5.1: An analysis of direction ifnding capability for high-precision location services</article-title>
          ,
          <source>Sensors</source>
          <volume>21</volume>
          (
          <year>2021</year>
          ). URL: https: //www.mdpi.com/1424-8220/21/11/3589. doi:
          <volume>10</volume>
          .3390/s21113589.
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>A.</given-names>
            <surname>Alarifi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Al-Salman</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Alsaleh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Alnafessah</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Al-Hadhrami</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M. A.</given-names>
            <surname>Al-Ammar</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H. S.</given-names>
            <surname>Al-Khalifa</surname>
          </string-name>
          ,
          <article-title>Ultra wideband indoor positioning technologies: Analysis and recent advances</article-title>
          ,
          <source>Sensors</source>
          <volume>16</volume>
          (
          <year>2016</year>
          ). URL: https://www.mdpi.com/1424-8220/16/5/707. doi:
          <volume>10</volume>
          . 3390/s16050707.
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>H.</given-names>
            <surname>Durrant-Whyte</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            <surname>Bailey</surname>
          </string-name>
          ,
          <article-title>Simultaneous localization and mapping: part i</article-title>
          ,
          <source>IEEE Robotics Automation Magazine</source>
          <volume>13</volume>
          (
          <year>2006</year>
          )
          <fpage>99</fpage>
          -
          <lpage>110</lpage>
          . doi:
          <volume>10</volume>
          .1109/MRA.
          <year>2006</year>
          .
          <volume>1638022</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>J.</given-names>
            <surname>Fuentes-Pacheco</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Ruiz-Ascencio</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. M.</given-names>
            <surname>Rendón-Mancha</surname>
          </string-name>
          ,
          <article-title>Visual simultaneous localization and mapping: a survey</article-title>
          ,
          <source>Artificial intelligence review 43</source>
          (
          <year>2015</year>
          )
          <fpage>55</fpage>
          -
          <lpage>81</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>D.</given-names>
            <surname>Feng</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>He</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Zhuang</surname>
          </string-name>
          , X.-G. Xia,
          <article-title>Kalman-filter-based integration of imu and uwb for high-accuracy indoor positioning and navigation</article-title>
          ,
          <source>IEEE Internet of Things Journal</source>
          <volume>7</volume>
          (
          <year>2020</year>
          )
          <fpage>3133</fpage>
          -
          <lpage>3146</lpage>
          . doi:
          <volume>10</volume>
          .1109/JIOT.
          <year>2020</year>
          .
          <volume>2965115</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>A.</given-names>
            <surname>Poulose</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Kim</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D. S.</given-names>
            <surname>Han</surname>
          </string-name>
          ,
          <article-title>A sensor fusion framework for indoor localization using smartphone sensors and wi-fi rssi measurements</article-title>
          ,
          <source>Applied Sciences</source>
          <volume>9</volume>
          (
          <year>2019</year>
          ). URL: https: //www.mdpi.com/2076-3417/9/20/4379. doi:
          <volume>10</volume>
          .3390/app9204379.
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>G.</given-names>
            <surname>Hakobyan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Yang</surname>
          </string-name>
          ,
          <article-title>High-performance automotive radar: A review of signal processing algorithms and modulation schemes</article-title>
          ,
          <source>IEEE Signal Processing Magazine</source>
          <volume>36</volume>
          (
          <year>2019</year>
          )
          <fpage>32</fpage>
          -
          <lpage>44</lpage>
          . doi:
          <volume>10</volume>
          .1109/MSP.
          <year>2019</year>
          .
          <volume>2911722</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>Z.</given-names>
            <surname>Geng</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H.</given-names>
            <surname>Yan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Zhang</surname>
          </string-name>
          , D. Zhu,
          <article-title>Deep-learning for radar: A survey</article-title>
          ,
          <source>IEEE Access 9</source>
          (
          <year>2021</year>
          )
          <fpage>141800</fpage>
          -
          <lpage>141818</lpage>
          . doi:
          <volume>10</volume>
          .1109/ACCESS.
          <year>2021</year>
          .
          <volume>3119561</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>F.</given-names>
            <surname>Dellaert</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Fox</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.</given-names>
            <surname>Burgard</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Thrun</surname>
          </string-name>
          ,
          <article-title>Monte carlo localization for mobile robots</article-title>
          ,
          <source>in: Proceedings 1999 IEEE International Conference on Robotics and Automation (Cat. No.99CH36288C)</source>
          , volume
          <volume>2</volume>
          ,
          <year>1999</year>
          , pp.
          <fpage>1322</fpage>
          -
          <lpage>1328</lpage>
          . URL: https://ieeexplore.ieee.org/ document/772544. doi:
          <volume>10</volume>
          .1109/ROBOT.
          <year>1999</year>
          .
          <volume>772544</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <given-names>I.</given-names>
            <surname>Rekleitis</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A Particle</given-names>
            <surname>Filter</surname>
          </string-name>
          <article-title>Tutorial for Mobile Robot Localization</article-title>
          ,
          <source>Technical Report TRCIM-04-02</source>
          ,
          <article-title>Centre for Intelligent Machines</article-title>
          , McGill University, Montreal, Québec, Canada,
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>