Efficient 3D trilateration algorithm for visible light positioning
David Plets,Yousef Almadani,Sander Bastiaens,Muhammad Ijaz,Luc Martens,Wout Joseph
DOI: https://doi.org/10.1088/2040-8986/ab1389
IF: 2.1
2019-04-17
Journal of Optics
Abstract:This letter presents an efficient algorithm for estimating the three-dimensional (3D) location of a photodiode (PD) receiver via visible light positioning. It solely works on measured powers from different light-emitting diode (LED) sources and does not require any prior knowledge of the PD receiver height. It is found that four LEDs are required that are not on the same circle, in order to unambiguously determine the 3D location. The algorithm is optimized towards a minimized calculation time in view of real-time operation on energy-constrained lightweight and mobile devices such as drones.The advent of light-emitting diodes (LEDs) has sparked a large research interest in visible light positioning (VLP), whereby the location of a photodiode (PD) receiver is being estimated based on its received powers from different LEDs. For many applications, the PD is expected to maintain a fixed and known height, reducing the localization problem to a two-dimensional (2D) problem. This is either solved via (manual or model-based) fingerprinting maps or via a classic trilateration method. For three-dimensional (3D) problems however, a fingerprinting approach with cm-level granularity becomes unfeasible due to the requirement of large memory and computation time to iterate over the 3D map. Alternatively, 3D trilateration can be applied, for which a plethora of algorithms is available. In VLP however, the distance cannot be estimated directly from the received power, as this estimation also depends on the height difference between the LED and the PD. Therefore, some workarounds have been sought, e.g. one can assume that the PD height is either known or estimated from on-board sensors, leading to a '2.5D solution'. In [<a href="#joptab1389bib1">1</a>], a 3D algorithm is presented, using multiple receivers instead of multiple transmitters. In [<a href="#joptab1389bib2">2</a>], a full 3D VLP algorithm is presented, estimating the 2D position based on the previous height estimate, followed by a height adjustment of this 2D location. However, this introduces errors as the best 2D estimate at the previous height will be different from the one at the new height. In [<a href="#joptab1389bib3">3</a>], another full 3D approach is presented, but the method requires integral calculations and is limited to three sources, which will be shown further to not suffice for an unambiguous location estimation. The solution of [<a href="#joptab1389bib4">4</a>] is also limited to three sources, and is mathematically relatively complex. Compared to available research, this letter presents a novel simple and error-free method for estimating the full 3D (<em>x</em>, <em>y</em>, <em>z</em>) location of a single PD, i.e. without any prior knowledge of the PD receiver height. It is shown that three LEDs, or even four LEDs on a square's corner, do not allow an unambiguous location estimation. Moreover, thanks to its low complexity and an additional fast search optimization, the method is suited for real-time 3D VLP operation.We assume a generic configuration with <em>N</em> LEDs, each horizontally placed at a fixed height <em>h</em><sub>LED</sub>, denoted with LED<sub>i</sub> and with coordinates (<em>x</em><sub><em>i</em></sub><em>, y</em><sub><em>i</em></sub><em>, h</em><sub>LED</sub>), <em>i</em> = 1 <em>.. N</em>. The Lambertian mode of LED<sub>i</sub> is denoted as <em>m</em><sub><em>i</em></sub> and its transmitted power as <em>P</em><sub><em>Ti</em></sub>. The PD, with area <em>A</em>, is also horizontally oriented and located at the unknown location <em>(x</em>, <em>y</em>, <em>z)</em> at an unknown distance <em>d</em><sub><em>i</em></sub> from LED<sub>i</sub>. The angle of irradiance is defined as <em>θ</em><sub><em>i</em></sub> and the angle of incidence as as shown in figure <a href="#joptab1389f1">1</a>(a). A setup with four LEDs in a 5 <b>×</b> 5 <b>×</b> 5 m<sup>3</sup> volume is considered; see figure <a href="#joptab1389f1">1</a>(b). According to the channel model of [<a href="#joptab1389bib5">5</a>], the power <em>P</em><sub><em>Ri</em></sub> received from LED<sub>i</sub>, at <em>(x</em>, <em>y</em>, <em>z)</em> is given by:<strong>Figure 1.</strong> (a) Definition of configuration parameters and (b) considered setup with indication of the involved LED locations and evaluation grid.Download figure:Inversely, given a measured power <em>P</em><sub><em>Ri</em></sub> received at an unknown location, and knowing for horizontally oriented LEDs and PD (see figure <a href="#joptab1389f1">1</a>(a)), the distance <em>d</em><sub><em>i</em></sub> between the PD and LED<sub>i</sub> can be estimated as based on the power <em>P</em><sub><em>Ri</em></sub> received from LED<sub>i</sub> (<em>i</em> = 1 <em>.. N</em>):with Δ<em>h</em> = <em>h</em><sub>LED</sub>−<em>z</em> = the (unknown) height difference between LED<sub>i</sub> and the PD. As such, cannot be directly determined from <em>P</em> <p>-Abstract Truncated-</p>
optics