[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Loosening Monitoring of a Threaded Pipe Connection Using the Electro-Mechanical Impedance Technique—Experimental and Numerical Studies
Next Article in Special Issue
A Parallel Architecture for the Partitioning around Medoids (PAM) Algorithm for Scalable Multi-Core Processor Implementation with Applications in Healthcare
Previous Article in Journal
Data Collection Based on Opportunistic Node Connections in Wireless Sensor Networks
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient AoA-Based Wireless Indoor Localization for Hospital Outpatients Using Mobile Devices

1
School of Electronic Engineering, Xidian University, Xi’an 710071, China
2
School of Engineering, University of Glasgow, Glasgow G12 8QQ, UK
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(11), 3698; https://doi.org/10.3390/s18113698
Submission received: 27 August 2018 / Revised: 17 October 2018 / Accepted: 26 October 2018 / Published: 30 October 2018

Abstract

:
The motivation of this work is to help outpatients find their corresponding departments or clinics, thus, it needs to provide indoor positioning services with a room-level accuracy. Unlike wireless outdoor localization that is dominated by the global positioning system (GPS), wireless indoor localization is still an open issue. Many different schemes are being developed to meet the increasing demand for indoor localization services. In this paper, we investigated the AoA-based wireless indoor localization for outpatients’ wayfinding in a hospital, where Wi-Fi access points (APs) are deployed, in line, on the ceiling. The target position can be determined by a mobile device, like a smartphone, through an efficient geometric calculation with two known APs coordinates and the angles of the incident radios. All possible positions in which the target may appear have been comprehensively investigated, and the corresponding solutions were proven to be the same. Experimental results show that localization error was less than 2.5 m, about 80% of the time, which can satisfy the outpatients’ requirements for wayfinding.

1. Introduction

With the popularity of mobile devices, such as smartphones, smartwatches, tablets, as well as some wearable devices, wireless localization has spawned numerous location-based services (LBS) to facilitate our daily life in the fields of emergency services, navigation, healthcare, geofence, social networking, entertainment/gaming, asset tracking, etc. [1,2,3,4,5,6]. For example, smartphone users can easily plan their trips using Google Maps and order a nearby taxi. The accurate location is mainly provided by the global positioning system (GPS). However, GPS signal is degraded too severely to be received in indoor environments. Many methods have been developed to meet the increasing demand for indoor LBS applications [7,8,9]. These wireless indoor localization methods and systems can be roughly divided into two main categories: range-based schemes and range-free schemes.
Range-based schemes first measure the distances between the target and at least three reference points, and then determine the target’s position by triangulation or trilateration, using the least square method. The distances are indirectly measured by interpreting some physical parameters obtained from wireless signals, such as time of flight (ToF), time of arrival (ToA), time difference of arrival (TDoA), and received signal strength indication (RSSI) [10,11,12,13,14,15,16]. ToF, ToA, and TDoA are very sensitive to timing error, thus, relying on highly synchronized timers. Furthermore, they usually need work under line of sight (LOS) conditions to get accurate ranging. RSSI is used to figure out the distance between the transmitter and the receiver with log-normal distance path loss (LDPL) model [17]. However, indoor RSSI is subject to temporal and spatial fluctuations because of the multipath effect, usually providing a coarse-grained positioning. The RSSI is not robust to noise because signal attenuates with propagation distance while noise usually does not.
Range-free schemes can be further subdivided into a distance vector hop (DV-hop), proximity detection, fingerprint matching, and angle of arrival (AoA)-based methods [18,19,20,21,22,23,24]. In a wireless sensor network (WSN), hop counts are utilized to estimate the positions of sensors with initially-unknown location information, by using a priori knowledge of the absolute positions of a few sensors (named as anchor nodes) and inter-sensor measurements, such as distance and bearing [25,26,27]. In the proximity detection method, Bluetooth low energy (BLE) beacons are deployed to periodically broadcast advertising messages, which contain their positions and other information [19,28,29]. The nearby mobile devices detect the signal from the beacon and calculate the approximate distance to the beacon. Thus, it needs to deploy a mass of beacons to obtain an accurate location. Wi-Fi can be used in a similar way as the BLE beacons, but requires external power sources and more setup costs. Fingerprint matching is a pattern recognition approach. During the training phase, signal features of all regions of interest are collected to build a fingerprint database [22]. Then, the localization is a process of matching the measured fingerprint, at an unknown location, with those in the database. Finally, the location corresponding to the best-fitted fingerprint is returned. Obviously, fingerprint matching is time-consuming and it is labor-intensive to build and maintain the database. AoA, also known as direction of arrival (DoA), refers to the angle of incidence at which radio signals from a transmitter arrive at the receiver. Given two AoAs and their receivers’ positions, the two-dimensional coordinates of the transmitter can be determined by triangulation. Accurate AoA can only be obtained by directional antennas or an antenna array, which makes the receiver complex and expensive.
Since the RSSI can be easily obtained by common mobile terminals, such as smartphones, Wi-Fi devices, Bluetooth beacons, ZigBee sensors, and so on, it is widely used to study wireless localization in range-based and range-free schemes. The fundamental drawback of RSSI is that it only works well in LOS situations. Unlike the RSSI, which is an aggregated value of all of the subcarriers’ amplitudes, the channel state information (CSI) estimates the channel, on each subcarrier, in the frequency domain. The CSI is to RSSI what a rainbow is to a sunbeam, so it can depict a multipath propagation to some extent and possesses more stable and fine-grained localization capability. The CSI can estimate the AoA and the ToF for multipath components under non line of sight (NLOS) conditions [30]. However, it needs to modify the firmware or update the hardware on available commercial devices to obtain the CSI data, which is difficult to implement on handheld devices, like smartphones. More details about the RSSI and the CSI can be found in Reference [31]. It is worth noting that hybrid methods are increasingly introduced in order to enhance the localization performance, such as the RSSI/DV-hop localization algorithm for the WSN, using fingerprinting method to improve the AoA-based ranging technique, the RSSI/AoA-based target localization and tracking, and so on [27,32,33,34].
In this paper, we propose an efficient AoA-based wireless indoor localization method for the outpatients’ wayfinding, in hospital, using mobile devices, such as tablets, smartphones, smartwatches, and so on. The Wi-Fi access points (APs) with antenna arrays are hung on the ceiling, so that radio signals can reach targets, along straight lines, in most situations. When users inquire about their locations using mobile devices like smartphones, nearby APs will timely return the positions of the APs and the AoAs pointing to the users’ current positions. Then the smartphones efficiently figure out the users’ locations, through geometric calculations. Both the Wi-Fi APs and the smartphones are off-the-shelf. Experimental results show that indoor localization error is less than 2.5 m, about 80%, of the time in LOS environments.
The contributions of this work are: (1) The proposed AoA-based algorithm is to employ only two Aps, at a time, for target localization. (2) The solutions are proven to be in the same form, after all the possible positions where the target may appear are investigated, comprehensively. (3) The calculation only involves a simple tangent function and basic arithmetic operations, which can be efficiently performed on a handheld mobile device.
The rest of this paper is organized as follows. First, Section 2 reviews the related work. Then, Section 3 introduces the design method and localization algorithm, in detail. Section 4 gives the experimental results. Finally, Section 5 concludes this work and discusses future work.

2. Related Work

AoA can be estimated by finding the maximum or minimum signal strength, during the rotation of a directional or a non-ideal omnidirectional monopole antenna [35,36,37,38], or when measuring the phase difference of received signals in an antenna array [39,40,41,42].
For the directional antenna, it needs an additional mechanical unit to rotate, and takes a long period of time to sample the 360° data [43]. So, it is impractical to be equipped in a handheld platform, such as smartphones, smartwatches, tablets, and so on. Zhang et al. proposed an alternative approach to emulate the functionality of a directional antenna, and developed an outdoor AP localization system, called Borealis [44]. In Borealis, a user holds a smartphone in the hand and rotates the body by 360°, a group of RSSI data is recorded and further processed to estimate the AP direction. Borealis produces 30° angular error in simple LOS environment, up to 50° in complex LOS, and 65° in NLOS environments. SpinLoc adopts a similar localization principle as Borealis, except that it replaces the RSSI by the CSI for indoor scenarios [45]. SpinLoc can produce a median error of 20°, offering localization accuracies in order of 6.5 m, with 4 Aps, and up to 5 m, with more APs. Jiang et al. proposed an AoA Localization with RSSI Difference (ALRD) system to estimate the AoA by comparing the RSSI values, about the target, received from two perpendicularly-oriented directional antennas, at the beacon nodes [37]. ALRD has an average localization error of 1.24 m at the CDF (Cumulative Distribution Function) = 50% in a 10 × 10 m indoor area. Although the ALRD need not rotate the directional antennas to collect data, it has to rebuild a new fingerprint database of the RSSI values when the target changes.
Antenna array is another common approach to estimate the AoA, by measuring the phase difference between the received signals of the two antennas. The phase difference is caused by the time difference of the signal that arrives in the different antennas, and is essentially determined by the incident angle of the incoming signal. Amundson et al. employed radio interferometric technique to determine the bearing from an anchor node, with an antenna array, to a distant target node, achieving an accuracy of about 3° [39]. Based on the above radio interferometric AoA estimation, Amundson et al. designed a localization and navigation system, called the TripNav, which could provide an average position and heading accuracy of 0.95 m and 4.75°, when traveling indoor at a speed of 1 m/s [23]. However, it requires an accurate time synchronization between the anchor node and the target node, on the order of microseconds. Chen et al. designed a low-cost and flexible antenna array for the AoA estimation, using modular, commercial-off-the-shelf, software-defined radios (SDRs) [42]. The prototype system still needs an additional clock module to synchronize the connected SDRs, for measuring the phase difference of arrival (PDoA), and offers a localization error less than 3 m, in the open field test. In the indoor location system, called the ArrayTrack, a custom FPGA (Field Programmable Gate Array) hardware with sixteen antennas was designed, and a multi-path suppression algorithm, named MUSIC (Multiple Signal Classification), was applied for the AoA estimation [40]. ArrayTrack can pinpoint forty-one clients, over an indoor office environment, to within 0.23 m median accuracy. Tzur et al. used a commercial Intel 5300 Wi-Fi NIC (Network Interface Card) with two receiving antennas to develop a practical method for identifying the AoA of a Wi-Fi AP [46]. In a realistic indoor environment, they are able to acquire the AoA, with a median error of 8–15°. SpotFi uses a single NIC with three antennas, combining the CSI values across the subcarriers and the antennas to jointly estimate both the AoA and the ToA of each path [47]. The joint estimation procedure enables the SpotFi to provide a similar AoA accuracy that is comparable to systems that require twice as many antennas. Under the LoS conditions, SpotFi can achieve a median error of 5°. Gjengset et al. used multiple NICs and automatic calibration to achieve a median location accuracy of 0.9 m [41].

3. Methodology

In this section, we first introduce the basic idea of our method. And then we give the detailed description of all possible cases in which the target will appear with regards to the Wi-Fi APs. Finally, we show how to deal with the case in which a target can be accessed by several pairs of APs.

3.1. AoA-Based Localization with Two Known APs

We plot Figure 1 to interpret the AoA-based localization with two known APs. For the convenience of description, all positions and angles in this paper are expressed in the two-dimensional Cartesian coordinate system. And the AoA value is defined to be positive if it is counter-clockwise and negative otherwise. In order to avoid confusion, we further define that the AoA specifically refers to the angle of incidence at which the radio signal travels from the transmitter to the receiver, along a straight line.
As shown in Figure 1, two Wi-Fi APs are located at ( x 1 , y 1 ) and ( x 2 , y 2 ) , labelled as A P 1 and A P 2 , respectively. Without the loss of generality, let x 1 x 2 and y 1 y 2 , here. An arbitrary mobile target, labelled as T , is needed to be localized, by determining its coordinates ( x , y ) with the two known APs and the corresponding AoAs. From the standpoint of the target, the AoAs related to A P 1 and A P 2 are θ 1 and θ 2 , respectively, but if the AoAs are measured by A P 1 and A P 2 , they will be α 1 and α 2 .
As stated above, the undetermined coordinates ( x , y ) of target T can be formulated by the following tangent functions
  { y 1 y = ( x 1 x ) tan θ 1 y 2 y = ( x 2 x ) tan θ 2  
The solution to Equation (1) is
  { x = y 2 y 1 + x 1 tan θ 1 x 2 tan θ 2 tan θ 1 tan θ 2                 y = y 2 ( x 2 x 1 ) tan θ 1 ( y 2 y 1 ) tan θ 1 tan θ 2 tan θ 2  
It is easily observed in Figure 1 that two pairs of the AoAs, named as α 1 and α 2 , θ 1 , and θ 2 , have the following relations:
  { θ 1 = α 1 180 ° θ 2 = α 2 180 °  
Considering the periodicity of the tangent function, we can rewrite the solution to Equation (1) by replacing θ 1 and θ 2 with α 1 and α 2
  { x = y 2 y 1 + x 1 tan α 1 x 2 tan α 2 tan α 1 tan α 2                 y = y 2 ( x 2 x 1 ) tan α 1 ( y 2 y 1 ) tan α 1 tan α 2 tan α 2  
If APs are installed in a line, the solution will have a more concise form. For example, when the APs are in the same row, which means y 1 = y 2 , solutions in Equation (4) will be simplified as
  { x = x 1 tan α 1 x 2 tan α 2 tan α 1 tan α 2             y = y 2 ( x 2 x 1 ) tan α 1 tan α 2 tan α 1 tan α 2  
Unlike the least-squares solution, by using a complex matrix equation in Reference [48], the solutions in Equations (2), (4), and (5) show that our method can efficiently figure out the target’s position, through some relatively simple trigonometric calculations.

3.2. All Possible Positions of the Target around the APs

In this section, we will investigate how to figure out all the possible positions of a mobile target, around the Aps, in practice. When a user holds a mobile device (such as smartphone, smartwatch, tablet, etc.) to wander in a shopping mall, workshop, or hospital, he or she will appear in any position randomly. Different positions will get different AoAs. We group these positions into normal and special situations, according to the AoAs. In a normal situation, the tangent functions of the AoAs have non-zero values, while in the special situation they have zero or infinite values.
As shown in Figure 2, the normal situation is summarized into six cases. In order to be more intuitive, let y 1 = y 2 , so that A P 1 and A P 2 are aligned in the same row.
Case 1: α 1 ,   α 2 ( 0 ° , 90 ° ) , as shown in Figure 2a, the equation can be written as
  { y y 1 = ( x x 1 ) tan α 1 y y 2 = ( x x 2 ) tan α 2  
The solution to Equation (6) is
  { x = x 1 tan α 1 x 2 tan α 2 tan α 1 tan α 2             y = y 2 ( x 2 x 1 ) tan α 1 tan α 2 tan α 1 tan α 2  
Considering the geometric relations of θ 1 = α 1 + 180 ° and θ 2 = α 2 + 180 ° , and the periodicity of the tangent function, we can get another form of the solution
  { x = x 1 tan θ 1 x 2 tan θ 2 tan θ 1 tan θ 2             y = y 2 ( x 2 x 1 ) tan θ 1 tan θ 2 tan θ 1 tan θ 2  
Case 2:   α 1 ( 0 ° , 90 ° ) ,   α 2 ( 90 ° , 180 ° ) , as shown in Figure 2b, the equation can be written as
  { y y 1 = ( x x 1 ) tan α 1                 y y 2 = ( x 2 x ) tan ( 180 ° α 2 )  
As tan ( 180 ° α 2 ) = tan α 2 , Equation (9) can be rewritten as
  { y y 1 = ( x x 1 ) tan α 1 y y 2 = ( x x 2 ) tan α 2  
Equation (10) is the same as Equation (6), so their solutions are also the same.
Similarly, since θ 1 = α 1 + 180 ° and θ 2 = α 2 + 180 ° , we can also obtain the same solution as in Equation (8).
Case 3: α 1 ,   α 2 ( 90 ° , 180 ° ) , as shown in Figure 2c, the equation can be written as
  { y y 1 = ( x 1 x ) tan ( 180 ° α 1 ) y y 2 = ( x 2 x ) tan ( 180 ° α 2 )  
As tan ( 180 ° α 1 ) = tan α 1 and tan ( 180 ° α 2 ) = tan α 2 , Equation (11) can be rewritten as
  { y y 1 = ( x x 1 ) tan α 1 y y 2 = ( x x 2 ) tan α 2  
Equation (12) is the same as Equation (6), so their solutions are also the same.
Similarly, since θ 1 = α 1 + 180 ° and θ 2 = α 2 + 180 ° , we can also obtain the same solution as in Equation (8).
Case 4: α 1 ,   α 2 ( 180 ° , 270 ° ) , as shown in Figure 2d, the equation can be written as
  { y 1 y = ( x 1 x ) tan θ 1 y 2 y = ( x 2 x ) tan θ 2  
The solution to Equation (13) is
  { x = x 1 tan θ 1 x 2 tan θ 2 tan θ 1 tan θ 2             y = y 2 ( x 2 x 1 ) tan θ 1 tan θ 2 tan θ 1 tan θ 2  
Considering the geometric relations of θ 1 = α 1 180 ° and θ 2 = α 2 180 ° , and the periodicity of the tangent function, we can get another form of the solution
  { x = x 1 tan α 1 x 2 tan α 2 tan α 1 tan α 2             y = y 2 ( x 2 x 1 ) tan α 1 tan α 2 tan α 1 tan α 2  
Case 5: α 1 ( 270 ° , 360 ° ) ,   α 2 ( 180 ° , 270 ° ) , as shown in Figure 2e, the equation can be written as
  { y 1 y = ( x x 1 ) tan ( 180 ° θ 1 ) y 2 y = ( x 2 x ) tan θ 2                  
As tan ( 180 ° θ 1 ) = tan θ 1 , Equation (16) can be rewritten as
  { y 1 y = ( x 1 x ) tan θ 1 y 2 y = ( x 2 x ) tan θ 2  
Equation (17) is the same as Equation (13), so their solutions are also the same.
Similarly, since θ 1 = α 1 180 ° and θ 2 = α 2 180 ° , we can also obtain the same solution as Equation (15).
Case 6: α 1 ,   α 2 ( 270 ° , 360 ° ) , as shown in Figure 2f, the equation can be written as
  { y 1 y = ( x x 1 ) tan ( 180 ° θ 1 ) y 2 y = ( x x 2 ) tan ( 180 ° θ 2 )  
As tan ( 180 ° θ 1 ) = tan θ 1 and tan ( 180 ° θ 2 ) = tan θ 2 , Equation (18) can be rewritten as
  { y 1 y = ( x 1 x ) tan θ 1 y 2 y = ( x 2 x ) tan θ 2  
Equation (19) is the same as Equation (13), so their solutions are also the same.
Similarly, since θ 1 = α 1 180 ° and θ 2 = α 2 180 ° , we can also obtain the same solution as Equation (15).
So far, we have investigated all the six cases in a normal situation. No matter where the target locates around the two APs and no matter which pairs of the AoAs are used, the results show that the equations are essentially consistent and that solutions are the same.
Compared with the normal situation, the special situation consists of blind and pole points, as shown in Figure 3.
Blind points refer to the positions in which the tangent values of the AoAs are zero. According to the properties of the tangent function, when an angle is 0 ° or integer multiples of 180 ° , its corresponding tangent value is zero. Thus the positions in the same row with the APs, such as Z 1 , Z 2 , and Z 3 , are the blind points.
Pole points refer to the positions in which either of the tangent values of the AoAs is infinity. Mathematically, when an angle is integer multiples of 90 ° , its corresponding tangent value is infinity. Thus, the positions in the same column with either of the APs, such as P 1 , P 2 , P 3 , and P 4 , are pole points.
Obviously, any of blind points fails to constitute a triangle, along with A P 1 and A P 2 . Thus, it cannot figure out the coordinates of the blind points using the tangent functions. Users should avoid these positions when adopting the AoA-based localization algorithm. As for the pole points, it is easier to calculate the target’s coordinates. Taking P 3 for example, we can directly obtain x = x 1 by measuring θ 1 = 90 ° , so that we only need to figure out the vertical coordinate y . Using the tangent function, we have
  y 2 y = ( x 2 x ) tan θ 2  
Thus, y = y 2 ( x 2 x ) tan θ 2 . Along with x = x 1 , we get the position of P 3 as follows
  { x = x 1                                 y = y 2 ( x 2 x 1 ) tan θ 2  

3.3. Roaming Strategy

As shown in Figure 4, the target can be accessed by the four APs. There are six combinations in total for choosing two out of four APs, which are A P 0 and A P 1 , A P 1 and A P 2 , A P 1 and A P 3 , A P 0 and A P 2 , A P 0 and A P 3 , and A P 2 and A P 3 , respectively.
Although it is able to localize the target with the aid of any pair of the APs. For the sake of accuracy and efficiency, we adopt a simple but effective roaming strategy to deal with such uncertainty or collision. When there are more than two APs available for estimating the target’s AoAs, the two with the highest RSSI are used. As shown in Figure 4, the four APs are supposed to be of the same model and the same configuration. It is obvious that A P 1 and A P 2 have the higher RSSI measurements for the radio from target, compared with the other combinations, because the target is much closer in distance to A P 1 and A P 2 .

4. Experiments and Results

In this section, we first introduce the experimental setup and then carry out the performance evaluation.

4.1. Experimental Setup

In our experiment, the AoAs were measured with Cisco AironetTM 4800 Wi-Fi Aps (Cisco Systems Inc., San Jose, CA, USA), and were retrieved offline from the Cisco Connected Mobile Experiences (CMX) software. Aironet 4800 series APs were designed for the multiple-purpose applications, equipped with up to four radios: Flexible 2.4 GHz or 5 GHz for the Wi-Fi, dedicated 5 GHz for Wi-Fi, 2.4 GHz, and 5 GHz for the Hyperlocation, 2.4 GHz for the BLE. Hyperlocation is an enhanced location solution developed by Cisco, consisting of a few hardware instruments and software components. The workflow in Hyperlocation was: 4800 AP listened to the connected mobile devices through a circular array of sixteen antennas and created packets to record the measured phase differences; then, these packets were sent to the CMX, at a rate of about 1 packet per second per AP, via the WLC (Wireless Local Area Network Controller); finally, the CMX extracted the AoA information from the phase differences and performed the localization. It took 10–20 s to complete one localization, based on the AoA information, and provided a localization accuracy of 1–3 m, in the LOS conditions. Hyperlocation could also provide a real-time location service, relying only on the RSSI, if connected mobile device did not support the OFDM (Orthogonal Frequency Division Multiplexing) technique, but at the cost of a much poor accuracy. More details on the Aironet 4800 series AP and the Hyperlocation Solution can be found on the data sheet and the Cisco website [49,50].
Two 4800 APs were hung from a 3 m high ceiling, in an open laboratory, and an Android smartphone (Coolpad C106 Coolpad Group Limited, Shenzhen, China), with dual-band radio transceivers, supporting the IEEE 802.11a/b/g/n/ac standards) was used as the target. As shown in Figure 5, the laboratory was furnished with a few desks, chairs, and PCs to simulate a hospital waiting hall. The due east of the compass in the smartphone was defined as the positive half of the horizontal axis, and the angle was defined to be positive if it was counter-clockwise and negative, otherwise. In order to ensure the signal strength, the smartphone was connected to the 2.4 GHz Wi-Fi. When a user needed to inquire the location, via a smartphone, the user stopped walking and waited to receive the APs’ coordinates and their corresponding AoA measurements. With these parameters, the smartphone pinpointed the user’s current coordinates, as described in Section 3.

4.2. The Impact of the Distance and Direction

In order to investigate the impact of distance and direction on the localization accuracy, we chose a series of positions as the ground truth, which are marked with green squares in Figure 6.
As shown in Figure 6, G 1 to G 7 were located in the same direction, to the center of A P 1 and A P 2 , with a distance increasing from 10 m to 40 m. G 5 , G 8 , G 9 , and G 10 were located in the same distance to the center of A P 1 and A P 2 , and they were in symmetry, with respect to the horizontal or the vertical axis. G 11 to G 15 were random test positions.
The localization results are illustrated in Figure 6. The green squares represent the actual locations of the test positions, and the blue triangles are the estimated locations. Localization was performed for eight times, in every test position, to obtain the statistical outcome, as plotted in Figure 7. One group of the results are listed in Table 1, in which the localization error is the Euclidean distance between the estimated and the actual locations, which was obtained by the following.
Δ d = ( x E x A ) 2 + ( y E y A ) 2  
where, ( x E , y E ) represented the coordinates of the estimated location, and ( x A , y A ) were the coordinates of the actual location.
Figure 7 shows that the localization accuracy decreased as the distance increased, despite the fact that the minimum, maximum, median, and the mean values of the location errors fluctuated. However, it is worth noting that the median and mean location errors did not increase monotonically. For example, most estimations at the distance of 15 m or 20 m, had less location errors than those at the distance of 10 m. This phenomenon was not solely caused by a measurement deviation, but mainly resulted from the ignored fact that the target and the APs were located in a three-dimensional space, rather than on a two-dimensional plane. In other words, the APs provided the AoA measurements in a three-dimensional space, which were mixture of the azimuth and elevation. In our calculations, they were regarded as the azimuth angles. From a mathematical point of view, the farther the distance between the target and APs, the smaller the error will be when the three points of different heights were approximated to be on the XOY plane.
As shown in Table 1, it was obvious that location accuracy depended on the AoAs. Another noteworthy phenomenon was that the G 5 , G 8 , G 9 , and G 10 were located, symmetrically, at the same distance to the center of the A P 1 and A P 2 , but their corresponding AoAs errors and location errors were quite different. This was mainly related to the coverage or directionality of the antenna array. Theoretically, an antenna should be isotropic with a uniform spherical coverage. In fact, no antenna could radiate radio waves, equally, in all directions. In addition, the posture that a user holds the smartphone in, to communicate with the APs and the surroundings, contributed to the deviation in the location errors.

4.3. Comparison with RSSI Range-Based Method

We repeated the localization using the RSSI range-based method in the positions G 1 to G 7 , in order to make a comparison with our proposed AoA-based method. As only two APs were employed, as described in Section 4.1, we redesigned the RSSI range-based localization scheme as follows:
(1)
The user held the smartphone and walked up or back toward the AP, ranging from 0 to 50 m, and recorded the measured RSSI at a step size of 0.5 m.
(2)
The RSSI was fitted a function of the user’s distance to the AP, as illustrated in Figure 8.
(3)
The user held the smartphone and stood facing the two APs in one of the positions G 1 to G 7 , and recorded the RSSI values, corresponding to the A P 1 and A P 2 .
(4)
The distances ( R 1 and R 2 ) to A P 1 and A P 2 were retrieved, based on the function obtained in (2).
(5)
The user’s location could be estimated by solving the following set of equations:
  { ( x x 1 ) 2 + ( y y 1 ) 2 =   R 1 2   ( x x 2 ) 2 + ( y y 2 ) 2 =   R 2 2    
It was obvious that the solutions to Equation (23) were the intersections of the two circles. Two circles were centered on A P 1 and A P 2 , with R 1 and R 2 as the radius, respectively. If the intersection existed, there were usually two solutions to Equation (23), except for the special situation in which there was only one solution, when the two circles were tangent to each other (i.e., ( x 2 x 1 ) 2 + ( y 2 y 1 ) 2 2 = R 1 + R 2 ).
In order to pick the right solution from the two, the blocking effect described in References [44,45] were utilized. As the human body can significantly attenuate Wi-Fi signals (both in 2.4 GHz and 5 GHz frequencies), the RSSI values on the smartphone would vary when the user held smartphone and spun at the location. As shown in Figure 8b, when the user lay between the AP and the smartphone, the direct path was severely blocked, resulting in a large attenuation of the RSSI. With the help of the blocking effect, additional direction information with regard to the APs was obtained, to determine the right solution. Localization was also performed for eight times, in each position, to obtain the statistical result, as plotted in Figure 9, and the comparison with our proposed AoA-based method is given in Table 2.
As shown in Figure 9, a similar phenomenon, as observed in Figure 7, occurred again when the localization accuracy decreased as the distance increased, despite the fluctuating minimum, maximum, median, and mean values of the location errors, which is more intuitional in Table 2. As both localization methods were essentially based on the RSSI information measured in a LOS-dominating environment, it was to be expected that they would have an equivalent location accuracy. However, compared to the solution of the binary quadratic equations in the RSSI range-based method, our proposed AoA-based method was much more efficient, as the calculation only involved a simple tangent function and basic arithmetic operations. Furthermore, the solution was unique in our method, even if only two APs were used for positioning.

4.4. Localization Accuracy

We tested the localization accuracy of our proposed algorithm in a more general environment, where the APs were randomly placed on the ceiling and the user wandered around the APs. The results are shown in Figure 10. Figure 11 plots the cumulative distribution function (CDF) of the median localization errors.
As shown in Figure 11, for over 80% of the test positions, localization error was less than 2.5 m.

5. Conclusions and Discussion

In order to provide room-level localization services for hospital outpatients using their own mobile devices, this paper proposed an efficient AoA-based wireless indoor localization algorithm that uses the coordinates and the angles of the arrival, provided by commercial Wi-Fi access points, to estimate the target’s position. Different from the available triangulation or trilateration localization algorithms, which utilize three or more APs and the complicated least square method, the proposed localization algorithm employs only two Aps, at a time, for positioning. All the possible positions in which the target might appear were comprehensively investigated, and the corresponding solutions were proven to be in the same form. Compared to the solution for the binary quadratic equations in the RSSI range-based method, the calculation in our proposed AoA-based method for positioning, only involves a simple tangent function and basic arithmetic operation, which could be efficiently performed on handheld mobile devices, such as smartphones. Two experiments were conducted to evaluate the feasibility and the localization performance. Experimental results showed that localization error was less than 2.5 m, about 80% of the time, in the LOS-dominating environments, which could satisfy the outpatients’ needs for wayfinding, in the hospital.
However, the drawbacks were also obvious. The main disadvantage was that the infrastructure cost would be much higher to build a wireless local area network, using the Cisco 4800 series Aps, rather than deploying other conventional APs. The localization algorithm is also mainly designed for the LOS-dominating environments, which confines its application to a certain extent, in real practice. Therefore, further work could be done to improve this localization method.
First, hybrid methods could be explored to tackle the positioning in the NLOS environments. For example, the AoA might be used in conjunction with a ToF. To be more specific, the phase difference between a pair of adjacent antennas is marked with a time stamp, by taking into account the relative flight time of their received incident radio waves. The AoA, which corresponds to the minimum ToF, could be considered to be the AoA representing the LOS path.
Second, the time spent in positioning needs to be reduced for the sake of a better user experience. In other word, the user should be able to keep on walking at a normal pace during the positioning. The AoAs measurement and acquisition could be accelerated by developing specific applications to get rid of the manual operation.
Last, the APs’ deployment could be optimized to eliminate the blind points. For example, the APs could be staggered in the corridors or hall. When the target happens to be collinear with the two nearest APs, one of the AP should be replaced by another that is nearby.

Author Contributions

Funding acquisition, Y.H.; Project administration, X.Y.; Validation, Q.H.A.; Writing—original draft, Y.H.

Funding

This research is funded by the National Natural Science Foundation of China under Grant No. 61671349.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

  1. Yim, J.; Ganesan, S.; Kang, B.H. Location-based mobile marketing innovations. Mob. Inf. Syst. 2017, 2017. [Google Scholar] [CrossRef]
  2. Carrizales-Villagómez, S.Y.; Nuño-Maganda, M.A.; Rubio-Loyola, J. A Platform for e-Health Control and Location Services for Wandering Patients. Mob. Inf. Syst. 2018, 2018. [Google Scholar] [CrossRef]
  3. Yang, C.-S.; Tsai, P.-W.; Liao, M.-Y.; Huang, C.-C.; Yeh, C.E. Location-based mobile multimedia push system. In Proceedings of the 2010 International Conference on Cyber-Enabled Distributed Computing and Knowledge Discovery (CyberC), Huangshan, China, 10–12 October 2010; IEEE: Piscataway, NJ, USA, 2010; pp. 181–184. [Google Scholar]
  4. Hinze, A.; Bainbridge, D. Location-triggered mobile access to a digital library of audio books using Tipple. Int. J. Digit. Libr. 2016, 17, 339–365. [Google Scholar] [CrossRef]
  5. Wu, L.; Stilwell, M.A. Exploring the marketing potential of location-based mobile games. J. Res. Interact. Mark. 2018, 12, 22–44. [Google Scholar] [CrossRef]
  6. Huh, J.H.; Seo, K. An indoor location-based control system using bluetooth beacons for IoT systems. Sensors 2017, 17, 2912. [Google Scholar] [CrossRef] [PubMed]
  7. Yassin, A.; Nasser, Y.; Awad, M.; Al-Dubai, A.; Liu, R.; Yuen, C.; Raulefs, R.; Aboutanios, E. Recent Advances in Indoor Localization: A Survey on Theoretical Approaches and Applications. IEEE Commun. Surv. Tutor. 2017, 19, 1327–1346. [Google Scholar] [CrossRef] [Green Version]
  8. Brena, R.F.; García-Vázquez, J.P.; Galván-Tejada, C.E.; Muñoz-Rodriguez, D.; Vargas-Rosales, C.; Fangmeyer, J. Evolution of indoor positioning technologies: A survey. J. Sens. 2017, 2017. [Google Scholar] [CrossRef]
  9. Silva, P.; Kaseva, V.; Lohan, E. Wireless Positioning in IoT: A Look at Current and Future Trends. Sensors 2018, 18, 2470. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, D.; Fattouche, M. OFDM transmission for time-based range estimation. IEEE Signal Process. Lett. 2010, 17, 571–574. [Google Scholar] [CrossRef]
  11. Zhao, G.; Wang, D.; Fattouche, M. Time sum of arrival based BLUE for mobile target positioning. Adv. Sci. Lett. 2011, 4, 165–167. [Google Scholar] [CrossRef]
  12. Tiwari, S.; Wang, D.; Fattouche, M.; Ghannouchi, F. A Hybrid RSS/TOA Method for 3D Positioning in an Indoor Environment. ISRN Signal Process. 2012, 2012. [Google Scholar] [CrossRef]
  13. Podevijn, N.; Plets, D.; Trogh, J.; Martens, L.; Suanet, P.; Hendrikse, K.; Joseph, W. TDoA-Based Outdoor Positioning with Tracking Algorithm in a Public LoRa Network. Wirel. Commun. Mob. Comput. 2018, 2018. [Google Scholar] [CrossRef]
  14. Nurminen, H.; Dashti, M.; Piché, R. A survey on wireless transmitter localization using signal strength measurements. Wirel. Commun. Mob. Comput. 2017, 2017. [Google Scholar] [CrossRef]
  15. Zhu, H.; Alsharari, T. An improved RSSI-based positioning method using sector transmission model and distance optimization technique. Int. J. Distrib. Sens. Netw. 2015, 11, 587195. [Google Scholar] [CrossRef]
  16. Passafiume, M.; Maddio, S.; Cidronali, A. An improved approach for RSSI-based only calibration-free real-time indoor localization on IEEE 802.11 and 802.15.4 wireless networks. Sensors 2017, 17, 717. [Google Scholar] [CrossRef] [PubMed]
  17. Mao, G.; Anderson, B.D.O.; Fidan, B. Path loss exponent estimation for wireless sensor network localization. Comput. Netw. 2007, 51, 2467–2483. [Google Scholar] [CrossRef]
  18. Jia, Z.; Wu, C.; Li, Z.; Zhang, Y.; Guan, B. The Indoor Localization and Tracking Estimation Method of Mobile Targets in Three-Dimensional Wireless Sensor Networks. Sensors 2015, 15, 29661–29684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Zhuang, Y.; Yang, J.; Li, Y.; Qi, L.; El-Sheimy, N. Smartphone-based indoor localization with bluetooth low energy beacons. Sensors 2016, 16, 596. [Google Scholar] [CrossRef] [PubMed]
  20. He, S.; Chan, S.-H.G. Wi-Fi fingerprint-based indoor positioning: Recent advances and comparisons. IEEE Commun. Surv. Tutor. 2016, 18, 466–490. [Google Scholar] [CrossRef]
  21. Wang, X.; Gao, L.; Mao, S.; Pandey, S. CSI-based fingerprinting for indoor localization: A deep learning approach. IEEE Trans. Veh. Technol. 2017, 66, 763–776. [Google Scholar] [CrossRef]
  22. Xiang, P.; Ji, P.; Zhang, D. Enhance RSS-Based Indoor Localization Accuracy by Leveraging Environmental Physical Features. Wirel. Commun. Mob. Comput. 2018, 2018. [Google Scholar] [CrossRef]
  23. Amundson, I.; Sallai, J.; Koutsoukos, X.; Ledeczi, A. Mobile sensor waypoint navigation via RF-based angle of arrival localization. Int. J. Distrib. Sens. Netw. 2012, 2012, 1–15. [Google Scholar] [CrossRef]
  24. Schüssel, M. Angle of Arrival Estimation using WiFi and Smartphones. In Proceedings of the International Conference on Indoor Positioning and Indoor Navigation (IPIN), Alcalá de Henares, Spain, 4–7 October 2016; Volume 4, p. 7. [Google Scholar]
  25. Song, G.; Tam, D. Two Novel DV-Hop Localization Algorithms for Randomly Deployed Wireless Sensor Networks. Int. J. Distrib. Sens. Netw. 2015, 2015. [Google Scholar] [CrossRef]
  26. Xu, Y.; Luo, X.; Wang, W.; Zhao, W. Efficient DV-HOP localization for wireless cyber-physical social sensing system: A correntropy-based neural network learning scheme. Sensors 2017, 17, 135. [Google Scholar] [CrossRef] [PubMed]
  27. Cheikhrouhou, O.; Bhatti, G.M.; Alroobaea, R. A hybrid DV-hop algorithm using RSSI for localization in large-scale wireless sensor networks. Sensors 2018, 18, 1469. [Google Scholar] [CrossRef] [PubMed]
  28. Paek, J.; Ko, J.; Shin, H. A measurement study of BLE iBeacon and geometric adjustment scheme for indoor location-based mobile applications. Mob. Inf. Syst. 2016, 2016. [Google Scholar] [CrossRef]
  29. Kriz, P.; Maly, F.; Kozel, T. Improving indoor localization using bluetooth low energy beacons. Mob. Inf. Syst. 2016, 2016. [Google Scholar] [CrossRef]
  30. Ahmed, A.U.; Arablouei, R.; de Hoog, F.; Kusy, B.; Jurdak, R.; Bergmann, N. Estimating angle-of-arrival and time-of-flight for multipath components using wifi channel state information. Sensors 2018, 18, 1753. [Google Scholar] [CrossRef] [PubMed]
  31. Yang, Z.; Zhou, Z.; Liu, Y. From RSSI to CSI: Indoor localization via channel response. ACM Comput. Surv. 2013, 46, 1–32. [Google Scholar] [CrossRef]
  32. Jachimczyk, B.; Dziak, D.; Kulesza, W.J. Using the fingerprinting method to customize RTLS based on the AoA ranging technique. Sensors 2016, 16, 876. [Google Scholar] [CrossRef] [PubMed]
  33. Tomic, S.; Beko, M.; Dinis, R.; Gomes, J.P. Target tracking with sensor navigation using coupled RSS and AOA measurements. Sensors 2017, 17, 2690. [Google Scholar] [CrossRef] [PubMed]
  34. Tomic, S.; Beko, M.; Dinis, R.; Bernardo, L. On target localization using combined RSS and AoA measurements. Sensors 2018, 18, 1266. [Google Scholar] [CrossRef] [PubMed]
  35. Kim, M.; Chong, N.Y. RFID-based mobile robot guidance to a stationary target. Mechatronics 2007, 17, 217–229. [Google Scholar] [CrossRef]
  36. Kim, M.; Chong, N.Y. Direction Sensing RFID Reader for Mobile Robot Navigation. IEEE Trans. Autom. Sci. Eng. 2009, 6, 44–54. [Google Scholar]
  37. Jiang, J.R.; Lin, C.M.; Lin, F.Y.; Huang, S.T. ALRD: AoA localization with RSSI differences of directional antennas for wireless sensor networks. Int. J. Distrib. Sens. Netw. 2013, 2013, 304–309. [Google Scholar] [CrossRef]
  38. Malajner, M.; Cucej, Z.; Gleich, D. Angle of arrival estimation using a single omnidirectional rotatable antenna. In Proceedings of the 2012 IEEE International Conference on 2012 Wireless Information Technology and Systems (ICWITS), Maui, HI, USA, 11–16 November 2012; Volume 12, pp. 1950–1957. [Google Scholar] [CrossRef]
  39. Amundson, I.; Sallai, J.; Koutsoukos, X.; Ledeczi, A. Radio interferometric angle of arrival estimation. Lect. Notes Comput. Sci. 2010, 5970 LNCS, 1–16. [Google Scholar] [CrossRef]
  40. Xiong, J.; Jamieson, K. ArrayTrack: A fine-grained indoor location system. In Proceedings of the Usenix Conference on Networked Systems Design and Implementation, Lombard, IL, USA, 2–5 April 2013; pp. 71–84. [Google Scholar]
  41. Gjengset, J.; Xiong, J.; McPhillips, G.; Jamieson, K. Phaser: Enabling phased array signal processing on commodity WiFi access points. In Proceedings of the 20th Annual International Conference on Mobile Computing and Networking, Maui, HI, USA, 7–11 September 2014; Volume 19, pp. 153–164. [Google Scholar] [CrossRef]
  42. Chen, H.C.; Lin, T.H.; Kung, H.T.; Lin, C.K.; Gwon, Y. Determining RF angle of arrival using COTS antenna arrays: A field evaluation. In Proceedings of the 2013 IEEE Military Communications Conference, San Diego, CA, USA, 18–20 November 2013; pp. 1–6. [Google Scholar]
  43. Jiang, J.-R.R.; Lin, C.-M.M.; Hsu, Y.-J. Localization with Rotatable Directional Antennas for Wireless Sensor Networks. In Proceedings of the 2010 39th International Conference on IEEE Parallel Processing Workshops (ICPPW), San Diego, CA, USA, 13–16 September 2010; pp. 542–548. [Google Scholar]
  44. Zhang, Z.; Zhou, X.; Zhang, W.; Zhang, Y.; Wang, G.; Zhao, B.Y.; Zheng, H. I am the antenna: Accurate outdoor AP location using smartphones. In Proceedings of the 17th Annual International Conference on Mobile Computing and Networking (MobiCom), Las Vegas, NV, USA, 20–22 September 2011; pp. 109–120. [Google Scholar]
  45. Sen, S.; Choudhury, R.R.; Nelakuditi, S. SpinLoc: Spin once to know your location. In Proceedings of the the Twelfth Workshop on Mobile Computing Systems & Applications, San Diego, CA, USA, 28–29 February 2012; ACM: New York, NY, USA, 2012; p. 12. [Google Scholar]
  46. Tzur, A.; Amrani, O.; Wool, A. Direction Finding of rogue Wi-Fi access points using an off-the-shelf MIMO-OFDM receiver. Phys. Commun. 2015, 17, 149–164. [Google Scholar] [CrossRef]
  47. Kotaru, M.; Joshi, K.; Bharadia, D.; Katti, S. SpotFi: Decimeter level localization using wifi. In ACM SIGCOMM Computer Communication Review; ACM: New York, NY, USA, 2015; Volume 45, pp. 269–282. [Google Scholar]
  48. Sayed, A.H.; Tarighat, A.; Khajehnouri, N. Network-based wireless location: Challenges faced in developing techniques for accurate wireless location information. IEEE Signal Process. Mag. 2005, 22, 24–40. [Google Scholar] [CrossRef]
  49. Cisco Systems. Cisco Aironet 4800 Access Point, Data Sheet; Cisco Systems, Inc.: San Jose, CA, USA, 2018; pp. 1–13. [Google Scholar]
  50. Cisco Hyperlocation Solution. Available online: https://www.cisco.com/c/en/us/solutions/enterprise-networks/hyperlocation-solution/index.html (accessed on 26 August 2018).
Figure 1. The angle of arrival (AoA)-based localization with two known access points (APs).
Figure 1. The angle of arrival (AoA)-based localization with two known access points (APs).
Sensors 18 03698 g001
Figure 2. Six normal cases of the target’s position. (a) Case 1, (b) Case 2, (c) Case 3, (d) Case 4, (e) Case 5, (f) Case 6.
Figure 2. Six normal cases of the target’s position. (a) Case 1, (b) Case 2, (c) Case 3, (d) Case 4, (e) Case 5, (f) Case 6.
Sensors 18 03698 g002aSensors 18 03698 g002b
Figure 3. Special cases of the target’s position.
Figure 3. Special cases of the target’s position.
Sensors 18 03698 g003
Figure 4. The roaming strategy.
Figure 4. The roaming strategy.
Sensors 18 03698 g004
Figure 5. The indoor test environment.
Figure 5. The indoor test environment.
Sensors 18 03698 g005
Figure 6. Indoor localization results.
Figure 6. Indoor localization results.
Sensors 18 03698 g006
Figure 7. Indoor localization errors.
Figure 7. Indoor localization errors.
Sensors 18 03698 g007
Figure 8. The measured RSSI profiles. (a) When the user holds the smartphone and walks; (b) when the user holds smartphone and spins.
Figure 8. The measured RSSI profiles. (a) When the user holds the smartphone and walks; (b) when the user holds smartphone and spins.
Sensors 18 03698 g008
Figure 9. Indoor localization errors using the received signal strength indication (RSSI) range-based method.
Figure 9. Indoor localization errors using the received signal strength indication (RSSI) range-based method.
Sensors 18 03698 g009
Figure 10. Indoor localization results.
Figure 10. Indoor localization results.
Sensors 18 03698 g010
Figure 11. The cumulative distribution function (CDF) of localization error.
Figure 11. The cumulative distribution function (CDF) of localization error.
Sensors 18 03698 g011
Table 1. Indoor localization results.
Table 1. Indoor localization results.
PositionActual AoAsMeasured AoAsLocation Error Δ d  
  α 1   α 2   α 1   α 2 ActualEstimated
  G 1   17.7 °   138.2 °   20 °   140 ° (7.07, 7.07)(5.92, 7.62)1.27
  G 2   22.5 °   112.5 °   25 °   115 ° (10.61, 10.61)(10.65, 9.33)1.28
  G 3   25.9 °   93.4 °   25 °   95 ° (14.14, 14,14)(13.82, 13.44)0.77
  G 4   28.4 °   81.4 °   30 °   80 ° (17.68, 17.68)(18.40, 19.28)1.76
  G 5   30.3 °   73.7 °   33 °   77 ° (21.21, 21.21)(20.29, 22.92)1.94
  G 6   31.9 °   68.5 °   34 °   70 ° (24.75, 24.75)(24.76, 26.82)2.07
  G 7   33.2 °   64.8 °   36 °   68 ° (28.28, 28.28)(27.47, 30.85)2.70
  G 8   106.3 °   149.6 °   105 °   150 ° (−21.21, 21.21)(−20.49, 20.49)1.02
  G 9   253.7 °   210.4 °   256 °   210 ° (−21.21, −21.21)(−20.04, −20.23)1.52
  G 10   329.5 °   286.3 °   330 °   286 ° (21.21, −21.21)(20.95, −20.76)0.52
Table 2. Indoor localization results.
Table 2. Indoor localization results.
PositionLocation Error Using Our Proposed MethodLocation Error Using RSSI Range-Based Method
MaxMinMeanMedianMaxMinMeanMedian
  G 1 1.870.101.271.272.180.381.331.20
  G 2 1.350.261.051.302.850.721.431.26
  G 3 1.760.171.101.212.470.751.671.62
  G 4 2.240.271.491.732.681.001.872.04
  G 5 2.680.521.541.483.241.132.332.23
  G 6 3.330.502.302.695.091.773.163.18
  G 7 3.290.422.402.945.221.423.313.26

Share and Cite

MDPI and ACS Style

Hou, Y.; Yang, X.; Abbasi, Q.H. Efficient AoA-Based Wireless Indoor Localization for Hospital Outpatients Using Mobile Devices. Sensors 2018, 18, 3698. https://doi.org/10.3390/s18113698

AMA Style

Hou Y, Yang X, Abbasi QH. Efficient AoA-Based Wireless Indoor Localization for Hospital Outpatients Using Mobile Devices. Sensors. 2018; 18(11):3698. https://doi.org/10.3390/s18113698

Chicago/Turabian Style

Hou, Yanbin, Xiaodong Yang, and Qammer H. Abbasi. 2018. "Efficient AoA-Based Wireless Indoor Localization for Hospital Outpatients Using Mobile Devices" Sensors 18, no. 11: 3698. https://doi.org/10.3390/s18113698

APA Style

Hou, Y., Yang, X., & Abbasi, Q. H. (2018). Efficient AoA-Based Wireless Indoor Localization for Hospital Outpatients Using Mobile Devices. Sensors, 18(11), 3698. https://doi.org/10.3390/s18113698

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop