# Evaluation of a Novel Radar Based Scanning Method

Authors: Paul Fritsche; Bernardo Wagner
Journal: Journal of Sensors
 (2016)

---

## Abstract

The following paper introduces a novel scanning
method for mapping and localization purposes in mobile robotics.
Our method is based on a rotating monostatic radar network,
which determines the positions of objects around the scanner via
a continuously running lateration algorithm. The estimation of
surfaces with ultrawideband radar networks has been studied
experimentally in lab environments, especially with lateration,
envelopes of spheres, and SEABED algorithms. But we do not
see a link to the field of mapping and localization of mobile
robots, where laser scanners are dominating. Indeed, only few
research groups use radars for mapping and localization, but
their applied sensor principle is based on a rotating focused radar
beam. Consequently, only 2D radar scanners are known inside
the robotic world and methods for 3D scanning with radars need
to be investigated. This paper will derive the theoretical background of the
sensor principle, which is based on a radar network on a rotating
joint, and discuss its erroneous influences. We were performing
first scans of standard geometries and deriving a model in
order to compare theoretical and experimental measurement
results. Furthermore, we present first mapping approaches and
a simulation of a scanner with multiple sensors.


---

## Body

## 1. Introduction

Obtaining all-around range information of an environment is essential in many areas of mobile robotics. Commonly, popular sensors like laser scanners, sonar sensors, and stereocameras have established themselves as state-of-the-art for most tasks in mobile robotics. Nevertheless, radar sensors frequently appear in field robotics but are seldom used to perform tasks like mapping and localization. Radar can penetrate certain materials, basically nonconductors, which provides advantages in dusty, foggy, rainy, or other harsh environments. But limited resolution, noisy data, and influence of optical effects like refraction, reflection, and absorption make the application in mobile robotics challenging.The use of radar sensors in mobile robotics is challenging but not impossible. The first appearance of radar sensors in the robotic community is traced back to the Australian Centre for Field Robotics in the early nineties, where fundamental work on probabilistic SLAM algorithms in combination with radar was developed (Clark and Whyte [1]). Because of their limited resolution and other aforementioned drawbacks, radar sensors are not very suitable to use in indoor environments. Nevertheless, Detlefsen et al. [2] were investigating the use of radar sensors in an industrial environment and Marck et al. [3] in an office. As far as we can see, all radar sensor principles in mobile robotics are based on mechanical beam-forming. Usually, the radar beam is focussed via a parabolic antenna and panned mechanically over the environment. Electrical beam-forming through phased array antennas is not seen very often in mobile robotics, but rather in automotive systems of the car industry.Besides beam-forming techniques, position estimation can be achieved through lateration, which is a common technique in radar networks for aircraft surveillance. Lateration is a measurement method, where the position of a point target is calculated through distance information fromn sensors with known location. The term trilateration refers to the measurement of three distances to define the position of an object (in contrast to triangulation, where three angles are used to calculate an object’s position). There exist two types of radar networks. In case of a monostatic radar network, the transmitter and receiver of the radar signal are collocated at the same location and can only receive signals that have been emitted by themselves. Multistatic radar networks consist of transmitters and receivers at different locations and can receive other sensors’ signal after it has been reflected from an object.In this paper, we introduce a scanning method, which is based on a rotating monostatic radar network. We use frequency modulated continuous wave (FMCW) radar sensors, which provide distance but no angle information of objects inside the observation area. The sensors work in 24 GHz ISM band and accordingly are limited in Germany to a bandwidth (B) of 250 MHz, which corresponds to a theoretical distance resolution of 0.6 m (see (1), c
o—speed of light). But the real distance resolution of most radar sensors is by the factor two or three larger. The availability of sensors with a high resolution depends on national and international bandwidth regulations. An ultrawideband (UWB) channel between 22 GHz and 26,65 GHz has been closed in 2013, but it is moved to 79 GHz for automotive purposes recently ([4, p. 20]):(1)
Δ
d
=
c
0
2
B
.The resolutionΔ
d of a radar sensor is equal to its minimum detection range. A radar’s resolution is its capability to distinguish objects. If the difference between the radial distances of two or more objects to the sensor is less than its resolution, then the sensor merges the two or more distance information into one. Additionally, the detection of objects depends on their radar cross-section (RCS) and the background noise of the environment.This paper is organized as follows. In Section2.1, we present a short overview on how position estimation via lateration in radar networks is commonly solved. In Section 2.2, we describe a problem that arises from data association and how to apply the bottom-up data association method of Fölster and Rohling [5] to resolve this problem. Section 2.3 describes the influences of errors in a radar network. In our first experiment, which will be described in Section 3, we were performing first scans, whose results will be presented and discussed in Section 4.


## 2. Materials and Methods

Estimating the position of an object with a radar network can be solved by standard lateration methods. For example, in order to define an object’s position in two-dimensional space, at least two sensors are necessary. Two radii can break down the object’s position to two possible locations. Usually, only one location is plausible due to the antenna’s direction. Geometrically, ghost objects can appear in lateration networks, which represents a wrong data association of the sensor’s object lists. A precise derivation, how ghost objects appear, is given by Rabe et al. [6].Besides lateration techniques, the envelopes of spheres and SEABED method exist. This method has been studied by Sakamoto and Kidera et al. [7, 8], who use it as a surface estimation technique for three-dimensional imaging with a UWB pulse radar system with high resolution. Like Sakamoto and Kidera, we assume boundary scattering on metallic surfaces to describe the model of our scanner.


### 2.1. Principle and Derivation

For terminology, we define the number of sensors in a radar network asn and the sensor-index i. Accordingly, S
i are the sensors with the coordinates x
S
i and y
S
i. Every sensor outputs an object-list O
L
i
=
[
d
i
1
,
d
i
2
,
…
,
d
i
m
], which contains distance information d
i
j from the sensor S
i to an object O
j.For a radar network in two-dimensional space and two sensors (see Figure1), we obtain two equations:(2)
x
S
1
-
x
O
1
2
+
y
S
1
-
y
O
1
2
=
d
11
2
x
S
2
-
x
O
1
2
+
y
S
2
-
y
O
1
2
=
d
21
2
.Figure 1
Our FMCW radar sensors measure radial distances to objects inside the observation area. The position of an objectO
j can be estimated through the intersection of two radii.In case the number of sensors is higher than the dimension of the space plus one (e.g., three sensors in two-dimensional space), then the system of equations becomes overdetermined. Due to errors in each sensor measurement, the overdetermined system of equations has no exact solution; hence a regression has to be found. The common way to solve this problem is the minimum mean square method. A detailed derivation and example of the method can be seen in the dissertation of Schneider [9, p. 11–14]. A general representation of the system of equations for a sensor network in two-dimensional space is given in (3). An expansion to three dimensions is self-explanatory. In order to estimate the object’s location, the system needs to be solved for o
→:(3)
A
=
2
·
x
S
1
-
x
S
n
2
·
y
S
1
-
y
S
n
2
·
x
S
1
-
x
S
n
2
·
y
S
1
-
y
S
n
⋮
2
·
x
S
n
-
1
-
x
S
n
2
·
y
S
n
-
1
-
y
S
n
o
→
j
=
x
O
j
y
O
j
d
→
j
=
x
S
1
2
-
x
S
n
2
+
y
S
1
2
-
y
S
n
2
-
d
1
j
2
-
d
n
j
2
x
S
2
2
-
x
S
n
2
+
y
S
2
2
-
y
S
n
2
-
d
2
j
2
-
d
n
j
2
⋮
x
S
n
-
1
2
-
x
S
n
2
+
y
S
n
-
1
2
-
y
S
n
2
-
d
n
-
1
;
j
2
-
d
n
j
2
A
·
o
→
j
=
d
→
j
.


### 2.2. Ghost Objects: A Data Association Problem

If more than one object is located inside the observation area of the radar network and the distance of the sensors to each other is larger than two times the sensor’s resolution, then the so-called ghost objects can appear (see Figure2). Ghost objects represent a wrong data association of the distance information from the object lists due to geometrical ambiguity. For our experiment, the ghost object issue is not relevant because we chose a distance between the sensors d
A
n
t smaller than 1,2 m. But if the experiments were performed with UWB radars, then ghost object would be an issue.Figure 2
Two objects (blue) are located in front of three sensors. In this case, due to ambiguity, six ghost objects appear. The observation area is discretized into grid cells. By calculating for every cell the error according to (4), an optimal data association can be performed.Fölster and Rohling [5] present a method called bottom-up data association. In two-dimensional space, at least three sensors are necessary for this method. In order to distinguish ghost objects from real objects, the observation area in front of the sensor network is discretized into cells, which represent a finite set of possible object positions. Then, a simple minimum distance calculation is done. Each cell contains an error value E
(
x
,
y
), which represents the square of the minimum distance of the cell to the sensor S
i minus the distance d
i
j between object O
j and sensor S
i, summarized over all n sensors (see (4)). This calculation results in the lowest error values in cells that are closest to the real objects:(4)
E
x
,
y
=
∑
i
=
1
n
min
d
ij
∈
OL
i
⁡
d
i
j
-
d
x
,
y
2
.The grid can be represented in Cartesian or polar coordinates. The grid size should be chosen under the consideration of the radar network’s range resolution and azimuth-angle accuracy [5].


### 2.3. Considerations on Erroneous Influences on Position Estimation in Monostatic Radar Networks

Understanding and analysis of errors that influence our sensor principle are essential to evaluate and discuss the feasibility and experimental results. Estimating the position of objects with lateration requires sensors with very high range accuracy. Nevertheless, every sensor has a measurement error. In case of lateration, the maximum position measurement errorσ
P
D can be approximated by the maximal range measurement error σ
R
m
a
x of all sensors and the angle α (see (5)). Figure 3 clarifies the relation between the range measurement error σ
R, the angle α, and the position measurement error σ
P
D graphically. This figure displays the assumed case that σ
R is constant. A closer look on σ
R is done in the following and shows that it can not be assumed to be constant; indeed, it is impossible to predict. From (5) and Figure 3, it can be seen that the accuracy is getting very bad at the sides of the sensor network, where α is approaching zero. An enhancement of the accuracy can be achieved by a larger distance between the sensors:(5)
σ
PD
≈
σ
R
max
⁡
sin
⁡
α
.Figure 3
(a) Presentation of the distribution of the position estimation errorσ
P
D over an observation area with two sensors at S
1
(
12.3,0
) and S
2
(
12.7,0
), in case of a constant σ
R of = 0.03 m. Although the assumption of a constant σ
R is not correct, this diagram aims to clarify the effect that σ
P
D gets very high at the sides of the radar network; hence the accuracy of our sensor principle is bad, if an object is entering from the side of the observation area, which is usually the case in our rotating sensor network. (b) Area of ambiguity, which is caused by the two rms errors σ
R
1 and σ
R
2.
(a)
(b)The range measurement accuracy is defined by the root-mean-square (rms) measurement errorσ
R [10, p. 167]. According to (6), the range measurement error is formed of the root-sum-square of three error components. The dominating component is the S/N-dependent random range measurement error σ
R
N. Its standard deviation is given by (7). The fixed random range error σ
R
F is the remaining range error, which remains if the S/N ratio gets very high. It represents the remaining error, caused by the sensor and sensor electronics architecture in case of a perfect S/N ratio. Range bias errors σ
R
B are constant for n
M
C measurements. In case of our sensor principle, the bias error does not affect the function or needs to be considered in lateration algorithms because it only results in a scaling that can be calibrated. In general, the accuracy of the sensor’s range measurement can be enhanced by increasing the number of measurements n
M
C. If measurements, which are recorded during a scan of one object, can be averaged will be discussed in Section 4. In Section 2.3, it is explained that the root standard deviation σ
R is not constant and depends on unpredictable influences. Mathematically, if different groups of random variables with different standard deviations have the same expected value, then together they still have the same expected value and a weighted standard deviation. This rule makes it legitimate to average over measurement cycles if radar network and object are static. In case of a dynamic radar network, we need to consider the non-point-target-case (NPTC), in order to clarify the characteristics of the expected value in dependency of the angle of the radar network. The NPTC is based on the fact that a radar is not measuring the distance to the same point of an object, if the sensor is placed at different locations. Besides obtaining n
M
C measurements from different measurement cycles, an averaging can be achieved as well by increasing the number of sensors; hence it was shown in Section 2.1 that overdetermined systems of equations are solved by a regression. A normal distribution of the range measurement values of our sensors has been investigated and confirmed, but will not be explained in further details in this paper:(6)
σ
R
=
σ
RN
2
n
MC
+
σ
RF
2
n
MC
+
σ
RB
2
.The dominating component of (6) is the S/N-dependent random range error σ
R
N. Interestingly, the sensor’s bandwidth has influence on the sensor’s resolution and accuracy as well. So far, estimating the radar’s accuracy seems to be feasible. But a closer look at the signal-to-noise ratio is leading to a complex relation, which is impossible to resolve in our experimental setup. Without going into details of the derivation of the S/N ratio, the proportionality in (8) is sufficient to know for our experiment:(7)
σ
RN
=
Δ
R
2
·
S
/
N
=
c
0
2
B
2
S
/
N
.From the proportionality in (8), we can summarize the following facts. First, the S/N ratio is higher and results accordingly in a better accuracy, if the RCS (σ) of an object is high. Consequently, our radar principle results in better position estimations for objects with high RCS. But objects with a high RCS are entering the observation area earlier from the sides than objects with a lower RCS; hence an object with high RCS suffers more from the geometrical issue described in Figure 3(a). Another problem is that the RCS cannot be assumed to be constant for an object. Even for standard geometries, for example, corner reflectors, σ differs from the aspect of view. From a historical point, research on the RCS and its dependence on the aspect angle has been already performed in the year 1946 by Robertson [11]. But there is not only a change of the RCS because of its geometric shape. There is a fluctuation of the RCS, which can be explained by the Swerling models. Ludloff explains in [12, p. 3–14] how the fluctuation can be modelled. The model is based on the idea that one radar target consists of multiple reflector elements, which are distributed over the volume of the target. The model assumes the reflector elements to be isotropic and with the same RCS and neglects the effects of reflection or shadows among themselves. Through overlapping of reflected radar waves on this multiple isotropic reflector elements, phase differences result in complex interferences. This model explains the appearance of high fluctuation of the RCS, even if the aspect angle is changed only slightly. To sum it up, an exact estimation of the RCS, even of standard geometries, is not possible in the real world and consequently needs to be represented by a probability. In addition to other influencing parameters, like the wave length λ of the radar’s center frequency, the object’s distance to the sensor R, and the thermal noise (see (8)), it leads us to the conclusion that it is not possible to give an exact estimation for the range measurement error σ
R:(8)
S
N
~
σ
λ
2
R
4
BT
S
.After having clarified that the position estimation errorσ
P
D depends on an unpredictable component, namely, the objects in the sensor’s observation area, we want to call attention to another source of error, which is caused by the sensor principle itself. We can see from the approximation in (5) that it is possible to lower the position estimation error σ
P
D by increasing the distance between the sensors. But, with a larger sensor distance, the influence of the NPTC raises.The NPTC occurs due to the fact that in most environments we cannot assume to have only point objects, whose dimensions are much smaller than the radar’s resolution. Normally, objects exceed the resolution cells of the sensors; in other words, we need to be able to handle surfaces with our sensor principle. That boundary surfaces can be estimated with UWB radars in combination with lateration algorithms as was shown by Mirbach and Menzel [13]. But, as can be seen in Figure 4, the estimation of the object position is wrong, if sensors measure the distance to different points of the object.Figure 4
Measuring the position of an object via lateration is based on the assumption of detecting with different sensors at different locations always the distances to the same point. This assumption is only acceptable for certain geometries or if the objects dimension is much smaller as the sensors’ resolution. In a real world environment, this becomes a surface estimation issue.


## 2.1. Principle and Derivation

For terminology, we define the number of sensors in a radar network asn and the sensor-index i. Accordingly, S
i are the sensors with the coordinates x
S
i and y
S
i. Every sensor outputs an object-list O
L
i
=
[
d
i
1
,
d
i
2
,
…
,
d
i
m
], which contains distance information d
i
j from the sensor S
i to an object O
j.For a radar network in two-dimensional space and two sensors (see Figure1), we obtain two equations:(2)
x
S
1
-
x
O
1
2
+
y
S
1
-
y
O
1
2
=
d
11
2
x
S
2
-
x
O
1
2
+
y
S
2
-
y
O
1
2
=
d
21
2
.Figure 1
Our FMCW radar sensors measure radial distances to objects inside the observation area. The position of an objectO
j can be estimated through the intersection of two radii.In case the number of sensors is higher than the dimension of the space plus one (e.g., three sensors in two-dimensional space), then the system of equations becomes overdetermined. Due to errors in each sensor measurement, the overdetermined system of equations has no exact solution; hence a regression has to be found. The common way to solve this problem is the minimum mean square method. A detailed derivation and example of the method can be seen in the dissertation of Schneider [9, p. 11–14]. A general representation of the system of equations for a sensor network in two-dimensional space is given in (3). An expansion to three dimensions is self-explanatory. In order to estimate the object’s location, the system needs to be solved for o
→:(3)
A
=
2
·
x
S
1
-
x
S
n
2
·
y
S
1
-
y
S
n
2
·
x
S
1
-
x
S
n
2
·
y
S
1
-
y
S
n
⋮
2
·
x
S
n
-
1
-
x
S
n
2
·
y
S
n
-
1
-
y
S
n
o
→
j
=
x
O
j
y
O
j
d
→
j
=
x
S
1
2
-
x
S
n
2
+
y
S
1
2
-
y
S
n
2
-
d
1
j
2
-
d
n
j
2
x
S
2
2
-
x
S
n
2
+
y
S
2
2
-
y
S
n
2
-
d
2
j
2
-
d
n
j
2
⋮
x
S
n
-
1
2
-
x
S
n
2
+
y
S
n
-
1
2
-
y
S
n
2
-
d
n
-
1
;
j
2
-
d
n
j
2
A
·
o
→
j
=
d
→
j
.


## 2.2. Ghost Objects: A Data Association Problem

If more than one object is located inside the observation area of the radar network and the distance of the sensors to each other is larger than two times the sensor’s resolution, then the so-called ghost objects can appear (see Figure2). Ghost objects represent a wrong data association of the distance information from the object lists due to geometrical ambiguity. For our experiment, the ghost object issue is not relevant because we chose a distance between the sensors d
A
n
t smaller than 1,2 m. But if the experiments were performed with UWB radars, then ghost object would be an issue.Figure 2
Two objects (blue) are located in front of three sensors. In this case, due to ambiguity, six ghost objects appear. The observation area is discretized into grid cells. By calculating for every cell the error according to (4), an optimal data association can be performed.Fölster and Rohling [5] present a method called bottom-up data association. In two-dimensional space, at least three sensors are necessary for this method. In order to distinguish ghost objects from real objects, the observation area in front of the sensor network is discretized into cells, which represent a finite set of possible object positions. Then, a simple minimum distance calculation is done. Each cell contains an error value E
(
x
,
y
), which represents the square of the minimum distance of the cell to the sensor S
i minus the distance d
i
j between object O
j and sensor S
i, summarized over all n sensors (see (4)). This calculation results in the lowest error values in cells that are closest to the real objects:(4)
E
x
,
y
=
∑
i
=
1
n
min
d
ij
∈
OL
i
⁡
d
i
j
-
d
x
,
y
2
.The grid can be represented in Cartesian or polar coordinates. The grid size should be chosen under the consideration of the radar network’s range resolution and azimuth-angle accuracy [5].


## 2.3. Considerations on Erroneous Influences on Position Estimation in Monostatic Radar Networks

Understanding and analysis of errors that influence our sensor principle are essential to evaluate and discuss the feasibility and experimental results. Estimating the position of objects with lateration requires sensors with very high range accuracy. Nevertheless, every sensor has a measurement error. In case of lateration, the maximum position measurement errorσ
P
D can be approximated by the maximal range measurement error σ
R
m
a
x of all sensors and the angle α (see (5)). Figure 3 clarifies the relation between the range measurement error σ
R, the angle α, and the position measurement error σ
P
D graphically. This figure displays the assumed case that σ
R is constant. A closer look on σ
R is done in the following and shows that it can not be assumed to be constant; indeed, it is impossible to predict. From (5) and Figure 3, it can be seen that the accuracy is getting very bad at the sides of the sensor network, where α is approaching zero. An enhancement of the accuracy can be achieved by a larger distance between the sensors:(5)
σ
PD
≈
σ
R
max
⁡
sin
⁡
α
.Figure 3
(a) Presentation of the distribution of the position estimation errorσ
P
D over an observation area with two sensors at S
1
(
12.3,0
) and S
2
(
12.7,0
), in case of a constant σ
R of = 0.03 m. Although the assumption of a constant σ
R is not correct, this diagram aims to clarify the effect that σ
P
D gets very high at the sides of the radar network; hence the accuracy of our sensor principle is bad, if an object is entering from the side of the observation area, which is usually the case in our rotating sensor network. (b) Area of ambiguity, which is caused by the two rms errors σ
R
1 and σ
R
2.
(a)
(b)The range measurement accuracy is defined by the root-mean-square (rms) measurement errorσ
R [10, p. 167]. According to (6), the range measurement error is formed of the root-sum-square of three error components. The dominating component is the S/N-dependent random range measurement error σ
R
N. Its standard deviation is given by (7). The fixed random range error σ
R
F is the remaining range error, which remains if the S/N ratio gets very high. It represents the remaining error, caused by the sensor and sensor electronics architecture in case of a perfect S/N ratio. Range bias errors σ
R
B are constant for n
M
C measurements. In case of our sensor principle, the bias error does not affect the function or needs to be considered in lateration algorithms because it only results in a scaling that can be calibrated. In general, the accuracy of the sensor’s range measurement can be enhanced by increasing the number of measurements n
M
C. If measurements, which are recorded during a scan of one object, can be averaged will be discussed in Section 4. In Section 2.3, it is explained that the root standard deviation σ
R is not constant and depends on unpredictable influences. Mathematically, if different groups of random variables with different standard deviations have the same expected value, then together they still have the same expected value and a weighted standard deviation. This rule makes it legitimate to average over measurement cycles if radar network and object are static. In case of a dynamic radar network, we need to consider the non-point-target-case (NPTC), in order to clarify the characteristics of the expected value in dependency of the angle of the radar network. The NPTC is based on the fact that a radar is not measuring the distance to the same point of an object, if the sensor is placed at different locations. Besides obtaining n
M
C measurements from different measurement cycles, an averaging can be achieved as well by increasing the number of sensors; hence it was shown in Section 2.1 that overdetermined systems of equations are solved by a regression. A normal distribution of the range measurement values of our sensors has been investigated and confirmed, but will not be explained in further details in this paper:(6)
σ
R
=
σ
RN
2
n
MC
+
σ
RF
2
n
MC
+
σ
RB
2
.The dominating component of (6) is the S/N-dependent random range error σ
R
N. Interestingly, the sensor’s bandwidth has influence on the sensor’s resolution and accuracy as well. So far, estimating the radar’s accuracy seems to be feasible. But a closer look at the signal-to-noise ratio is leading to a complex relation, which is impossible to resolve in our experimental setup. Without going into details of the derivation of the S/N ratio, the proportionality in (8) is sufficient to know for our experiment:(7)
σ
RN
=
Δ
R
2
·
S
/
N
=
c
0
2
B
2
S
/
N
.From the proportionality in (8), we can summarize the following facts. First, the S/N ratio is higher and results accordingly in a better accuracy, if the RCS (σ) of an object is high. Consequently, our radar principle results in better position estimations for objects with high RCS. But objects with a high RCS are entering the observation area earlier from the sides than objects with a lower RCS; hence an object with high RCS suffers more from the geometrical issue described in Figure 3(a). Another problem is that the RCS cannot be assumed to be constant for an object. Even for standard geometries, for example, corner reflectors, σ differs from the aspect of view. From a historical point, research on the RCS and its dependence on the aspect angle has been already performed in the year 1946 by Robertson [11]. But there is not only a change of the RCS because of its geometric shape. There is a fluctuation of the RCS, which can be explained by the Swerling models. Ludloff explains in [12, p. 3–14] how the fluctuation can be modelled. The model is based on the idea that one radar target consists of multiple reflector elements, which are distributed over the volume of the target. The model assumes the reflector elements to be isotropic and with the same RCS and neglects the effects of reflection or shadows among themselves. Through overlapping of reflected radar waves on this multiple isotropic reflector elements, phase differences result in complex interferences. This model explains the appearance of high fluctuation of the RCS, even if the aspect angle is changed only slightly. To sum it up, an exact estimation of the RCS, even of standard geometries, is not possible in the real world and consequently needs to be represented by a probability. In addition to other influencing parameters, like the wave length λ of the radar’s center frequency, the object’s distance to the sensor R, and the thermal noise (see (8)), it leads us to the conclusion that it is not possible to give an exact estimation for the range measurement error σ
R:(8)
S
N
~
σ
λ
2
R
4
BT
S
.After having clarified that the position estimation errorσ
P
D depends on an unpredictable component, namely, the objects in the sensor’s observation area, we want to call attention to another source of error, which is caused by the sensor principle itself. We can see from the approximation in (5) that it is possible to lower the position estimation error σ
P
D by increasing the distance between the sensors. But, with a larger sensor distance, the influence of the NPTC raises.The NPTC occurs due to the fact that in most environments we cannot assume to have only point objects, whose dimensions are much smaller than the radar’s resolution. Normally, objects exceed the resolution cells of the sensors; in other words, we need to be able to handle surfaces with our sensor principle. That boundary surfaces can be estimated with UWB radars in combination with lateration algorithms as was shown by Mirbach and Menzel [13]. But, as can be seen in Figure 4, the estimation of the object position is wrong, if sensors measure the distance to different points of the object.Figure 4
Measuring the position of an object via lateration is based on the assumption of detecting with different sensors at different locations always the distances to the same point. This assumption is only acceptable for certain geometries or if the objects dimension is much smaller as the sensors’ resolution. In a real world environment, this becomes a surface estimation issue.


## 3. Experiments

In order to evaluate the sensor principle, we were performing first scans of standard objects in an indoor environment. The goal of the first experiment is to find out about error influences in our sensor principle. As mentioned before, a limited resolution can be problematic in an indoor environment, for example, an office. There might be metallic radiators, steal-beams behind the walls, computer towers, and many other objects that can have a RCS huge enough of being detected by our radar sensor.Accordingly, the probability of occurrence of two objects, with a smaller difference of their radial distances to the sensor than the radar’s resolution, is high; hence we can rarely trust our scan results, if performed in an indoor environment. For fundamental research, our radar sensors with an resolution of approximately 0.6 m are sufficient. Of course, by using better sensors, a better result can be achieved.For our first experiment, we were performing two test series, where we placed a standard geometry, far away (at least two times the resolution) from disturbing objects as mentioned before. During experiment A, we were using a planar metallic board, and during experiment B, we were using a corner reflector. Both experiments will be explained in the following section.


### 3.1. Scan in front of a Metallic Planar Board and a Corner Reflector

During experiment A, we were scanning an indoor environment with our sensor unit. We were placing a planar metallic surface perpendicular to the sensor unit. We placed the board at the same height like the sensor unit in order to perform a 2D experiment in a 3D environment. Each measurement contains the accumulation of ten 360° degree scans with a step size of 0.7° degrees. Not every measurement cycle leads to a successful position estimation. A successful position estimation can be processed if both sensors detect an object. Experiment B is performed according to experiment A with a corner reflector. All relevant details of the experiments can be seen in Table1.Table 1
Experiment A.
Experiment
A
B
Object
Metallic planar surface (100
cm
×
45
cm)
Aluminium corner reflector
Antenna distance
d
Ant
=
40 cm
d
Ant
=
40 cm
Distance between sensor unit and object
d
0
≈
198 cm
d
0
≈
215 cmIn Figure6, the results of experiment A and experiment B are displayed. Further discussion of the scan results will be given in Section 4. A picture of the experimental environment is shown in Figure 5.Figure 5
In the front, the scanning unit is visible. In the back, the metallic planar surface and the corner reflector can be seen.Figure 6
Both diagrams present the results of scans inside an indoor environment. (a) A planar metallic surface (A) and (b) a corner reflector (B) have been positioned inside the scan area of the sensor. Both scans have been postprocessed with a Gaussian kernel density estimator in order to represent the probability of a detected objectp
(
o
).
(a)
(b)


### 3.2. Scan of a Hallway

In the second experiment, we want to find out if our sensor principle is suitable for robotic mapping. Therefore, we were recording raw data for a map of our hallway. To avoid influences of control and odometry errors of our robots, we were performing scans at known poses (see Figure7). The walls of the corridor have a distance of approximately 2 m.Figure 7
Ground truth of our office environment with scan positions of the sensor unit.


## 3.1. Scan in front of a Metallic Planar Board and a Corner Reflector

During experiment A, we were scanning an indoor environment with our sensor unit. We were placing a planar metallic surface perpendicular to the sensor unit. We placed the board at the same height like the sensor unit in order to perform a 2D experiment in a 3D environment. Each measurement contains the accumulation of ten 360° degree scans with a step size of 0.7° degrees. Not every measurement cycle leads to a successful position estimation. A successful position estimation can be processed if both sensors detect an object. Experiment B is performed according to experiment A with a corner reflector. All relevant details of the experiments can be seen in Table1.Table 1
Experiment A.
Experiment
A
B
Object
Metallic planar surface (100
cm
×
45
cm)
Aluminium corner reflector
Antenna distance
d
Ant
=
40 cm
d
Ant
=
40 cm
Distance between sensor unit and object
d
0
≈
198 cm
d
0
≈
215 cmIn Figure6, the results of experiment A and experiment B are displayed. Further discussion of the scan results will be given in Section 4. A picture of the experimental environment is shown in Figure 5.Figure 5
In the front, the scanning unit is visible. In the back, the metallic planar surface and the corner reflector can be seen.Figure 6
Both diagrams present the results of scans inside an indoor environment. (a) A planar metallic surface (A) and (b) a corner reflector (B) have been positioned inside the scan area of the sensor. Both scans have been postprocessed with a Gaussian kernel density estimator in order to represent the probability of a detected objectp
(
o
).
(a)
(b)


## 3.2. Scan of a Hallway

In the second experiment, we want to find out if our sensor principle is suitable for robotic mapping. Therefore, we were recording raw data for a map of our hallway. To avoid influences of control and odometry errors of our robots, we were performing scans at known poses (see Figure7). The walls of the corridor have a distance of approximately 2 m.Figure 7
Ground truth of our office environment with scan positions of the sensor unit.


## 4. Results and Discussion

In order to make a correct interpretation of our scan results, two models for the case of the board and the corner reflector have been derived. A model can be helpful when designing parameters of a system or understanding effects.


### 4.1. Model of the Scanned Objects

As in most physical models, we assume simplifications in order to break down the complexity of the problem. We assume that all radar waves are reflected from the object’s boundary surface; hence a penetration of the material is neglected. Furthermore, we assume that a single sensor is only measuring distances to surfaces, which are inside its observation area and perpendicular to the sensor itself. Double reflections, like those can appear between two parallel walls, are neglected. In case of orthogonal lines or areas, we assume a distance measuring to the corner point.Figures8 and 9 demonstrate a simplified understanding of the experimental setup. The value d
0 represents the distance between the center of rotation of the sensor unit and the object. In case of the planar metallic surface (Figure 8), we cannot expect that each sensor is measuring the distance to the same point due to a NPTC. According to our model, the distances can be calculated with (9), where X
i represents a random variable with normal distribution and standard deviation σ
R which is added to the geometrical part of the equation:(9)
d
1
A
=
d
0
-
d
Ant
2
sin
⁡
β
+
X
1
d
2
A
=
d
0
+
d
Ant
2
sin
⁡
β
+
X
2
.Figure 8
Each sensor measures the distance to a perpendicular surface. In case of a planar surface, the sensors of a radar network measure the distances to different points.Figure 9
Due to the reflective property of a corner reflector, the sensors of the radar network measure always the distance to the corner point; hence no NPTC occurs and it is legitimate to average measurement values of a scan.Due to the influence of the NPTC, it is not legitimate to average the measurement points of the planar metallic surface because the expected value, resulting from the lateration algorithm, is not constant and depends on the aspect angle.In case of a corner reflector, we can expect a distance measuring always to the same point (the corner point). Consequently, an averaging is legitimate. A measurement cycle in front of a corner reflector can be represented with the following equation:(10)
d
1
B
=
d
0
-
d
Ant
2
sin
⁡
β
2
+
d
Ant
2
cos
⁡
β
2
+
X
1
d
2
B
=
d
0
+
d
Ant
2
sin
⁡
β
2
+
d
Ant
2
cos
⁡
β
2
+
X
2
.Figure10 presents a comparison between the model of the planar surface and the measurements of the experiment A. We cannot confirm our model to be correct because of side effects and a nonperfectly planar surface. Nevertheless, we learn that the hugest error influence of our sensor unit is not the range measurement error σ
R. The biggest problem of our principle is the NPTC, which results in a wide spread of the measurement values. Figure 6 displays scan results with a KDE postprocessing, which result in a very high probability p
(
o
) for our scanned object’s locations. Accordingly, we obtain probabilities which describe areas, where an object could be located.Figure 10
(a) This diagram displays the characteristics of the distance measurementd
1
A (green) and d
2
A (red) of object A. The blue points represent the distance measurements of the model according to (9). It can be seen that model and experiment do not fit very well. Side effect from the edges of the metallic surface, the metallic stand, and a nonperfectly planar surface are reasons for the nonoptimal characteristics. (b) The diagram presents the distribution of the position estimation after the lateration algorithm. Again, the blue points represent the modelled result and the green ones the measurement results. Obviously, the NPTC situation is leading to the spread results, which can be seen already in Figure 6.
(a)
(b)Figure11 represents the measurement results of the corner reflector, which fit very well with our model. We can assume to have almost no influence of the NPTC situation. The remaining spread of the measurement values is caused by the range measurement error σ
R of the sensors.Figure 11
(a) The diagram represents the characteristics ofd
1
B (green) and d
2
B (red) of object B. The blue points represent the distance measurements of the model according to (11). Corner reflectors are very suitable for our scanning principle and can function as landmarks because they do not suffer under the NPTC. (b) Consequently, the experimental and modelled values result in very similar position estimations.
(a)
(b)


### 4.2. Mapping with Known Poses

There exist several mapping algorithms. An overview is given by Thrun in [14, p. 7]. Thrun introduces algorithms, which are suitable for mapping with unknown robot poses, which is named simultaneous localization and mapping (SLAM). In this paper, we focus on mapping with known poses, which is simpler. But mapping with known poses is leading to more promising results because odometry and control errors do not influence the map. Occupancy grid mapping with Bayes filter might be the most popular probabilistic representation of a map and will be the next step to combine with our sensor principle.Figure12 presents two simple approaches for grid mapping with our sensor principle. Figure 12(a) presents the measurement values of the hallway in a two-dimensional histogram. In general, double reflection between parallel walls or from the ground can occur and cause wrong position estimations. Due to the fact that we can expect more measurements of true objects than from wrong detection, a histogram is accumulating accordingly true object locations; hence the hallway’s shape is visible. The histogram representation can be fine-tuned via changing the bin size. Our second approach for further grid mapping algorithm is shown in Figure 12(b). We applied on every 360° scan a Gaussian KDE and normalized the values. Afterwards, we summed all normalized kernel estimated scans to obtain the map. This method can be fine-tuned with the kernel size of the estimator. In general, this approach is leading to a blurry representation of our map. We realized that we can build a dynamic inverse sensor model with both approaches, which will be the focus of further investigation.Figure 12
(a) Representation of the measurement points in a two-dimensional histogram. (b) A normalized Gaussian kernel density estimation for each 360° scan.
(a)
(b)


### 4.3.
3D Scan with Multiple Sensors

We simulated a scenario with six point targets and a scanner which is equipped with nine sensors. We used the simulation software V-REP from Coppelia Robotics, which allows the simulation of sonar sensor, which have a similar behaviour to radar sensors. We added to every distance measurement a Gaussian noise and solved the system of equation from (3) via the minimum mean square method. Furthermore, the simulation assumes sensors with a very high resolution. Accordingly, the simulation is used only to demonstrate the influence of the NPTC effect without influences of low resolution, when detecting multiple point targets with one sensor. The setup of the scenario can be seen in Figure 13.Figure 13
(a) This picture shows the simulated scene with six objects. Spheres are suitable as a point target. The scanner is equipped with 9 radar sensors and panning over all six objects from the left to the right. (b) Representation of the measurement points which are calculated via minimum mean square method. As it can be seen in the picture above, wrong position estimation of the spheres is the consequence of not detecting the same object with all sensors. (c) The filter criterion (see (11)) lowers the occurrence of wrong position estimations due to the NPTC. The red triangle in both diagrams represents the scanning unit.
(a)
(b)
(c)Due to NPTC, which can be seen in the upper picture of Figure13, the lateration algorithm results in many wrong position estimations, which can be seen as a spread of the estimated points. Therefore, we defined a criterion to filter the measurement. d
j is the distance between the average position of all sensors and the object position o
→:(11)
-
d
Ant
max
⁡
2
+
1
n
∑
i
=
1
n
d
i
j
<
d
j
<
d
Ant
max
⁡
2
+
1
n
∑
i
=
1
n
d
i
j
.The filter criterion enhances the result significantly. A comparison of the result with and without filter is shown at the bottom of Figure13.


## 4.1. Model of the Scanned Objects

As in most physical models, we assume simplifications in order to break down the complexity of the problem. We assume that all radar waves are reflected from the object’s boundary surface; hence a penetration of the material is neglected. Furthermore, we assume that a single sensor is only measuring distances to surfaces, which are inside its observation area and perpendicular to the sensor itself. Double reflections, like those can appear between two parallel walls, are neglected. In case of orthogonal lines or areas, we assume a distance measuring to the corner point.Figures8 and 9 demonstrate a simplified understanding of the experimental setup. The value d
0 represents the distance between the center of rotation of the sensor unit and the object. In case of the planar metallic surface (Figure 8), we cannot expect that each sensor is measuring the distance to the same point due to a NPTC. According to our model, the distances can be calculated with (9), where X
i represents a random variable with normal distribution and standard deviation σ
R which is added to the geometrical part of the equation:(9)
d
1
A
=
d
0
-
d
Ant
2
sin
⁡
β
+
X
1
d
2
A
=
d
0
+
d
Ant
2
sin
⁡
β
+
X
2
.Figure 8
Each sensor measures the distance to a perpendicular surface. In case of a planar surface, the sensors of a radar network measure the distances to different points.Figure 9
Due to the reflective property of a corner reflector, the sensors of the radar network measure always the distance to the corner point; hence no NPTC occurs and it is legitimate to average measurement values of a scan.Due to the influence of the NPTC, it is not legitimate to average the measurement points of the planar metallic surface because the expected value, resulting from the lateration algorithm, is not constant and depends on the aspect angle.In case of a corner reflector, we can expect a distance measuring always to the same point (the corner point). Consequently, an averaging is legitimate. A measurement cycle in front of a corner reflector can be represented with the following equation:(10)
d
1
B
=
d
0
-
d
Ant
2
sin
⁡
β
2
+
d
Ant
2
cos
⁡
β
2
+
X
1
d
2
B
=
d
0
+
d
Ant
2
sin
⁡
β
2
+
d
Ant
2
cos
⁡
β
2
+
X
2
.Figure10 presents a comparison between the model of the planar surface and the measurements of the experiment A. We cannot confirm our model to be correct because of side effects and a nonperfectly planar surface. Nevertheless, we learn that the hugest error influence of our sensor unit is not the range measurement error σ
R. The biggest problem of our principle is the NPTC, which results in a wide spread of the measurement values. Figure 6 displays scan results with a KDE postprocessing, which result in a very high probability p
(
o
) for our scanned object’s locations. Accordingly, we obtain probabilities which describe areas, where an object could be located.Figure 10
(a) This diagram displays the characteristics of the distance measurementd
1
A (green) and d
2
A (red) of object A. The blue points represent the distance measurements of the model according to (9). It can be seen that model and experiment do not fit very well. Side effect from the edges of the metallic surface, the metallic stand, and a nonperfectly planar surface are reasons for the nonoptimal characteristics. (b) The diagram presents the distribution of the position estimation after the lateration algorithm. Again, the blue points represent the modelled result and the green ones the measurement results. Obviously, the NPTC situation is leading to the spread results, which can be seen already in Figure 6.
(a)
(b)Figure11 represents the measurement results of the corner reflector, which fit very well with our model. We can assume to have almost no influence of the NPTC situation. The remaining spread of the measurement values is caused by the range measurement error σ
R of the sensors.Figure 11
(a) The diagram represents the characteristics ofd
1
B (green) and d
2
B (red) of object B. The blue points represent the distance measurements of the model according to (11). Corner reflectors are very suitable for our scanning principle and can function as landmarks because they do not suffer under the NPTC. (b) Consequently, the experimental and modelled values result in very similar position estimations.
(a)
(b)


## 4.2. Mapping with Known Poses

There exist several mapping algorithms. An overview is given by Thrun in [14, p. 7]. Thrun introduces algorithms, which are suitable for mapping with unknown robot poses, which is named simultaneous localization and mapping (SLAM). In this paper, we focus on mapping with known poses, which is simpler. But mapping with known poses is leading to more promising results because odometry and control errors do not influence the map. Occupancy grid mapping with Bayes filter might be the most popular probabilistic representation of a map and will be the next step to combine with our sensor principle.Figure12 presents two simple approaches for grid mapping with our sensor principle. Figure 12(a) presents the measurement values of the hallway in a two-dimensional histogram. In general, double reflection between parallel walls or from the ground can occur and cause wrong position estimations. Due to the fact that we can expect more measurements of true objects than from wrong detection, a histogram is accumulating accordingly true object locations; hence the hallway’s shape is visible. The histogram representation can be fine-tuned via changing the bin size. Our second approach for further grid mapping algorithm is shown in Figure 12(b). We applied on every 360° scan a Gaussian KDE and normalized the values. Afterwards, we summed all normalized kernel estimated scans to obtain the map. This method can be fine-tuned with the kernel size of the estimator. In general, this approach is leading to a blurry representation of our map. We realized that we can build a dynamic inverse sensor model with both approaches, which will be the focus of further investigation.Figure 12
(a) Representation of the measurement points in a two-dimensional histogram. (b) A normalized Gaussian kernel density estimation for each 360° scan.
(a)
(b)


## 4.3.
3D Scan with Multiple Sensors

We simulated a scenario with six point targets and a scanner which is equipped with nine sensors. We used the simulation software V-REP from Coppelia Robotics, which allows the simulation of sonar sensor, which have a similar behaviour to radar sensors. We added to every distance measurement a Gaussian noise and solved the system of equation from (3) via the minimum mean square method. Furthermore, the simulation assumes sensors with a very high resolution. Accordingly, the simulation is used only to demonstrate the influence of the NPTC effect without influences of low resolution, when detecting multiple point targets with one sensor. The setup of the scenario can be seen in Figure 13.Figure 13
(a) This picture shows the simulated scene with six objects. Spheres are suitable as a point target. The scanner is equipped with 9 radar sensors and panning over all six objects from the left to the right. (b) Representation of the measurement points which are calculated via minimum mean square method. As it can be seen in the picture above, wrong position estimation of the spheres is the consequence of not detecting the same object with all sensors. (c) The filter criterion (see (11)) lowers the occurrence of wrong position estimations due to the NPTC. The red triangle in both diagrams represents the scanning unit.
(a)
(b)
(c)Due to NPTC, which can be seen in the upper picture of Figure13, the lateration algorithm results in many wrong position estimations, which can be seen as a spread of the estimated points. Therefore, we defined a criterion to filter the measurement. d
j is the distance between the average position of all sensors and the object position o
→:(11)
-
d
Ant
max
⁡
2
+
1
n
∑
i
=
1
n
d
i
j
<
d
j
<
d
Ant
max
⁡
2
+
1
n
∑
i
=
1
n
d
i
j
.The filter criterion enhances the result significantly. A comparison of the result with and without filter is shown at the bottom of Figure13.


## 5. Conclusion

Robust localization and navigation in hazardous and tough environments are still a difficult issue in field robotics research. Dust, rain, fog, and inadequate illumination are conditions, which make popular sensors, such as laser scanners or cameras, not suitable. Radar overcomes the aforementioned difficulties.In this paper, we were investigating a new scanning method and took a closer look at failure influences in order to judge the principle’s feasibility. We were focusing on two influences: first, the range measurement error of the sensor itself and second, the influence of wrong position estimation due to the lateration principle (NPTC). We proved that the influence of the NPTC needs to be considered. Only objects, which are similar to corner reflectors, do not suffer under this effect.The lateration technique focusses on the estimation of points inside the observation area. A comparison with surface estimation algorithms like envelopes of spheres or SEABED might be an interesting point of investigation.We proved that our sensor principle is suitable for robotic mapping. An investigation on an dynamic inverse sensor model, which is obtained from kernel density estimations, will be the focus of our future research.To sum it up, our proposed principle is an alternative to standard radar based scanning-methods. Mechanical beam-forming techniques require an antenna and electric beam-forming techniques need phase array radars, which are commonly more expansive. Although no antenna construction is required, our principle needs more than one sensor. From one single 360° scan that is obtained through mechanical beam-forming, we can expect from each measurement step at an incremental angle distance information, which results in more measurement points than our principle, which detects only the objects with the highest RCS due to a nonfocused observation area. But our principle is recording more than one measurement of an object during one scan rotation, which raises the possibility of a correct detection of an object. An advantage over traditional rotating mechanical beam-forming techniques is the possibility to perform 3D scans as well, which would be mechanically complicated in case of mechanical beam-forming and is only known in combination with electrical beam-forming radars. But our principle suffers more from the NPTC, bad accuracy and resolution, and wrong calibration or asynchronism of measurements than traditional techniques and is accordingly limited.



---

*Source: 6952075-2015-11-22.xml*