An Analysis of the Impact of LED Tilt on Visible Light Positioning Accuracy

Whereas the impact of photodiode noise and reflections is heavily studied in Visible Light Positioning (VLP), an often underestimated deterioration of VLP accuracy is caused by tilt of the Light Emitting Diodes (LEDs). Small LED tilts may be hard to avoid and can have a significant impact on the claimed centimeter-accuracy of VLP systems. This paper presents a Monte-Carlo-based simulation study of the impact of LED tilt on the accuracy of Received Signal Strength (RSS)-based VLP for different localization approaches. Results show that trilateration performs worse than (normalized) Least Squares algorithms, but mainly outside the LED square. Moreover, depending on inter-LED distance and LED height, median tilt-induced errors are in the range between 1 and 6 cm for small LED tilts, with errors scaling linearly with the LED tilt severity. Two methods are proposed to estimate and correct for LED tilts and their performance is compared.


Introduction on Visible Light Positioning
The introduction of Light Emitting Diodes (LEDs) for traditional lighting applications has also led to new research lines in other application domains, as LEDs can be modulated to transmit information over the visible light propagation channel.Besides large interest in high-speed Visible Light Communication (VLC), another promising application is Visible Light Positioning (VLP) [1,2], where e.g., the location of a photo diode (PD) is estimated.When using Received Signal Strength Indicator (RSSI)-based positioning [3], VLP has an advantage over well-known Radio-Frequency (RF) solutions, thanks to the absence of small-scale fading effects.Research has shed light on the impact of noise on positioning performance [4], the impact of reflections [5,6], or the impact of LED power uncertainty [7].Also the impact of receiver tilt has been characterized [8][9][10][11], and ways to compensate for this, but it remains unclear to what extent random small tilts of the LED at the ceiling impact positioning performance.

Expected Issues Considering LED Tilt
Besides tilt induced via the LED die placement and packaging, another important source of tilt is the operator suspending the LED.As suspending a LED is typically done while looking up and standing on a ladder or a hydraulic platform lift, it becomes harder for the operator to maintain a good feeling of horizontality, compared to when standing on the ground.Intuitively, one would expect that the impact of LED tilt on the positioning error strongly relates to the lateral deviation in the receiver plane, defined as the distance between the intersect of the untilted LED normal and the receiver plane (i.e., right below the LED), and the intersect of the tilted LED normal and the receiver plane (i.e., where the LED beam is directed to).This lateral deviation is equal to the tangent of the tilt angle (in radians) multiplied by the LED-PD height difference, or, for small tilt angles, the product of the LED-PD height difference and the tilt angle (in radians) itself.Figure 1a shows the value of this deviation in the receiver plane due to a tilted LED, as a function of the LED height (i.e., the height difference between the LED and the PD).For real-life industrial VLP applications, LED heights of 7 m or more are not uncommon and lead to deviations between 12 cm (for a tilt of 1°) and 61 cm (for a tilt of 5°).This paper will quantitatively and qualitatively investigate the validity of this intuitive assumption.Furthermore, it needs investigation to what extent induced errors of multiple LEDs, each with a certain (unknown) tilt, accumulate or compensate each other.

Paper Content and Structure
This paper will characterize LED-tilt impact within a typical square VLP configuration by means of a Monte-Carlo simulation, in which it is assumed that each of the four deployed LEDs has a certain unknown (small) horizontal tilt, and a random azimuthal rotation.This way, the cumulative distribution function (cdf) of the positioning error will be constructed for different locations within the test site, and for different localization approaches.Also, the impact of LED tilt will be related with the LED height and inter-LED distance.Finally, a brief exploration is presented towards estimating LED tilt by investigating two tilt estimation methods, which can then be used to adjust model-based RSS-VLP fingerprinting maps, and eventually, reduce positioning errors.The remainder of this paper is structured as follows.Section 2 will present the visible light channel model that will be used in the positioning algorithm, the simulation configuration, the positioning algorithms, and the proposed LED tilt estimation methods.In Section 3, the results will be presented, after which the main findings of this work will be discussed in Section 4, together with related future research work.

Channel Model
In this work, only the Line-of-Sight (LoS) path between the LED transmitters and the PD receivers are accounted for.We will not consider the impact of reflections as a source of "noise", in order to be able to unambiguously assess the effect of LED tilt only.For the same reason, no shot noise or thermal noise will be considered in this study.It should be noted that even when reflections are accounted for via a model-based fingerprinting approach [6], LED tilt will have an impact since the reflected contributions could either increase or decrease, depending on the LED location and tilt, the PD location, and the location of the reflective surface.This is considered to be future work.Figure 2 and Table 1 define the model parameters of the visible light channel.The power P R received at the PD is calculated according to the channel model used in [12]: with P E the emitted optical power by the LED.h LoS is the channel gain along the direct link and can be described as follows, when assuming a Lambertian radiator: with m the order of the Lambertian emitter, and φ the angle of irradiance (i.e., the angle between the LED normal and the vector v LED2PD from the LED to the PD, with φ equal to either ϕ tilt or ϕ in Figure 2, depending on the LED being tilted or not).T R (ψ) and G R (ψ) are the optical filter's gain and the optical concentrator's gain at the receiver, respectively, with ψ the angle of incidence (i.e., the angle between the PD normal n PD and the vector from the PD to the LED).The LEDs will be assumed to be within the field-of-view (FOV) of the PD and T R (ψ) and G R (ψ) are assumed equal to 1. d is the distance between the LED and the PD, and A R the actual PD area, here assumed to be 1 cm 2 .The PD will be assumed to be horizontally oriented, so that cos(ψ) reduces to h/d, with h the height difference between the LED and the PD.When the LED is not tilted (i.e., horizontally oriented), the angle of irradiance φ = ϕ = ψ (see Figure 2) with cos(ϕ) also equal to h/d.However, in case of LED tilt, the angle of irradiance φ is equal to ϕ tilt (see Figure 2), with cos(ϕ tilt ) as follows: with n tilt LED the (tilted) LED normal (compared to n no tilt LED , being the untilted LED normal).The LED is assumed to be tilted over an angle θ (determining the severity of the tilt), and rotated over an angle α (determining where the LED is tilted to), as shown in Figure 2. Since | n tilt LED | = 1 and | v LED2PD | = d, cos(ϕ tilt ) is obtained as follows: with (x PD , y PD , 0) and (x LED , y LED , h) the coordinates of the PD and the LED respectively.Please note that the assumption of the PD being located in the xy-plane, does not retract from the generality of the work, as only the LED-PD height difference h matters.As such, the power P tilted R received from a tilted LED on a horizontally oriented PD is calculated as follows:

Simulation Configuration
Simulations will be executed for the room depicted in Figure 3.The dimensions of the room are 7 m × 7 m, with a ceiling at a height difference h above the PD, with h = 2.5 m (office environment) or h = 6 m (industrial environment).The Lambertian order m of the LEDs is equal to 1, for all four LEDs.In this scenario, we assume that the receiver height is fixed and known (e.g., a PD attached to the top of a cart), so the evaluation of the receiver location is reduced to a plane.A receiver grid of 5 mm will be considered here, meaning that the PD center can be located at N L = 1401 2 candidate locations.Furthermore, we assume that the receiver hardware is able to demultiplex the contributions of the different LED sources [13].Four LEDs are attached to the ceiling, at the locations indicated in Figure 3.We here assume the frequency division multiple access (FDMA) scheme proposed in [13] and further investigated in [14], which combines the transmission of square waves with 'power spectrum identification' based on the Fast Fourier Transform of the incident signal.The authors make use of the even harmonics of square waves being zero, to separate the different contributions of each LED transmitter at the receiver side.The modulation frequencies are typically in the range of 1-100 kHz, which is sufficiently below the LED and LED driver bandwidth (>1 MHz) to not have a detrimental impact on the positioning system.Although the exact optical power of the LEDs does not impact the findings, as no receiver noise is added, we here assume an optical power of 10 W, in line with common assumptions in literature [15,16].
The parameter θ is assumed to be the absolute value of a normally distributed variable θ N , and α to be uniformly distributed: . Two values of σ tilt will be mainly considered in this work: 1°and 2°, meaning that 95% of the LEDs will be mounted with a tilt below 2°and 4°respectively.
The positioning error due to tilt will be evaluated at 100 locations in a quarter of the receiver plane, i.e., at the blue locations indicated in Figure 3. Thanks to the symmetry of the setup, the resulting error distributions can be extrapolated to the three other zones, since the distribution of the positioning errors will repeat themselves at the corresponding locations of the different parts of the 7 m × 7 m area.It should be noted though that each single simulation will result in an asymmetric setup, since for each simulation, the tilt of each LED will be randomly and independently chosen.However, the resulting cdf at corresponding locations will be the same when a sufficient number of simulations is considered, due to the fact that the statistical distribution of the LED tilt is assumed to be the same for each of the LEDs.In this work, a position estimation for each of the 100 positions will be executed for 10,000 random sets of LED tilts in a Monte-Carlo simulation.Each position estimation will be done according to the algorithms presented in Section 2.3.

Positioning Algorithms
Three different positioning algorithms will be compared: a traditional trilateration method, and two model-based fingerprinting methods using a Least-Squares Estimation (LSE) and a normalized LSE (nLSE) respectively.

Trilateration
Given a measured power P meas Ri from LED i (i = 1..N) and assuming that cos(φ) = cos(ψ) = h d i for horizontally oriented LEDs and PD (see Figure 2), Equation ( 2) can be rewritten to allow the calculation of the estimated distance d i between LED i and PD: The real squared horizontal distance d 2 i between LED i and the PD is given by: with (x, y) the PD location, and (x i , y i ) the coordinates of LED i .After eliminating the quadratic terms in x 2 and y 2 by subtracting d 2 N from d 2 i , N-1 equations are obtained (i = 1, . . ., N − 1): These linear equations in x and y can be written as x y can then be estimated as x y using the Moore-Penrose pseudo-inverse of M: In case the LEDs are tilted, the assumption of cos(φ) = cos(ψ) = h d i will no longer hold, and errors will be introduced.

Least-Squares Estimation
The following two methods are based on the comparison of the set of so-called measured received powers P meas Ri from each (tilted) LED i (i = 1, . . ., N) at the unknown PD location (x PD , y PD ), with the set of fingerprinted PD powers P L,model Ri from LED i at all (1401 2 ) locations L in the grid.For the construction of the fingerprinting database of the P L,model Ri values, each LED is assumed to be untilted, as it is its most probable position.The set of so-called measurements (P meas R1 , P meas R2 , . . ., P meas RN ) represent the observed values in the realistic setup investigated here.They are obtained from (P E1 , P E2 , . . ., P EN ), where P meas Ri values are obtained from Equation ( 5) with θ and α values as samples from their respective statistical distributions.The larger the tilt of the LEDs (larger σ 2 tilt values), the larger the positioning errors will be.The algorithm estimates the unknown location with coordinates (x PD , y PD ) to be at the location L where the Least-Squares cost function C L LSE has a minimum [6]:

Normalized Least-Squares Estimation
The normalized Least-Squares Estimation algorithm was shown to perform better than the LSE algorithm when there is uncertainty on the respective LED powers [17].The unknown location (x PD , y PD ) is estimated at the the location L where the cost function C L nLSE has a minimum: For the LSE and nLSE methods, each position estimation thus consists of a comparison of the set of measurements (P meas R1 , P meas R2 , . . ., P meas RN ) against all (P L,model

R1
, P L,model R2 , . . ., P L,model RN ) sets that are stored in the database.In total, N L sets of N values of P L,model Ri are precalculated and stored in a fingerprinting database, i.e., the received power at N L locations from each of the N LEDs, according to the LoS channel model from Section 2.1.For the configuration under test, N = 4 and N L = 1401 2 = 1,962,801, meaning that 7,851,204 values are stored.

LED Tilt Estimation Methods
Based on the configuration and the algorithms described in the previous sections, the impact of LED tilt will be assessed.In this section, two methods for estimating the LED tilt are presented, i.e., the angles θ and α from Figure 2. Once this tilt is known, the fingerprinting map with powers (P L,model Ri values of Sections 2.3.2 and 2.3.3)can be adjusted, in order to reduce aforementioned LED tilt impact.

Exhaustive Search
The first method encompasses an exhaustive search over all possible (x, y) locations in the receiver plane.As Equation (2) can be rewritten as the location (x L ,y L ) in the receiver plane which is pointed to by the LED normal, can be determined as the (x,y) location for which the angle of irradiance φ is minimal or cos m (φ)(x, y) is maximal.For a horizontal PD, and when removing all (x,y)-independent factors, (x L ,y L ) is thus found at the location (x,y) where is maximal.Especially for small LED tilt values, it is fair to assume that the optical filter's gain and the optical concentrator's gain at the receiver are also independent of the receiver location, so that the LED tilt is determined by the location (x L ,y L ) where the product of the measured received power P R (x, y) and the cube of the LED-PD distance becomes maximal, i.e., From the location (x L ,y L ), the tilt values (α, θ) are easily derived using basic trigonometry as (arctan ).Although this method in principle delivers the exact LED tilt, it is very cumbersome, as it requires the execution of a large set of measurements, whereby the exact ground truth of the measurement location has to be known.The measurement campaign should be executed with an automated positioning system, e.g., with a stepping motor [18].

Quick Search
Unlike in the case of the first method, the second method, denoted here as 'quick search', uses a power measurement at only a limited set of N M locations under the considered LED.Based on these power measurements P meas i , a minimum-search in the (θ,α) space is executed for finding the most likely values for these tilt parameters.We propose a cost function using power ratios at the different locations instead of absolute power values, since the expected powers P model i are sometimes not known without doing a calibration phase like e.g., in [19].This is due to e.g., unknown deviations of the assumed LED power, PD responsivity,. . .When working with power ratios and thus excluding these unknown factors to a certain extent, the cost function is expected to be better suited to identify power differences that are solely due to the LED tilt.As such, the proposed cost function is the following: with N M the number of measurement locations in the receiver plane, and P meas i the measured power at location i underneath the considered LED.The power values P i , i=1..N M , are first sorted in descending order.P model i (α, θ) is the modeled received power at location i, according to Equation (5) for a LED tilt determined by (α, θ).The cost function iterates over α between 0°and 360°, with a 1°resolution.The parameter θ is iterated between 0°and θ max = 5, with 0.1°resolution.C tilt will reach its minimal value for (α, θ) values equal to the actual LED tilt values, although the outcome will be corrupted by external factors such as noise, unaccounted PD tilt, deviations from the tabulated Lambertian order.

LED Tilt Estimation Scenario
The performance of the described methods will be evaluated for a 10 W LED at a height of 2.5 m above the receiver plane.The LED is assumed to have a horizontal tilt of θ = 1.6°and an azimuthal rotation α of 230°.For the 'quick search' method, four measurement locations (N M = 4) in the receiver plane will be considered, with (quite randomly chosen) coordinates (1,0), (−0.5, ± √ 2 ), and (1, −1).Finally, a PD area of A R = 1 cm 2 is assumed.

Assessment of LED Tilt Impact on Positioning Accuracy for Typical Configurations Using Different Metrics
Figure 4 shows the cdf of the positioning errors of the 10 4 simulations at each of the 100 evaluation locations for (a) a LED-PD height difference of 2.5 m and σ tilt = 1°(denoted as normal office placement); (b) a LED-PD height difference of 2.5 m and σ tilt = 2°(denoted as sloppy office placement); (c) a LED-PD height difference of 6 m and σ tilt = 1°(denoted as normal industrial placement); and (d) a LED-PD height difference of 6 m and σ tilt = 2°(denoted as sloppy industrial placement).In each plot, the three localisation metrics from Section 2.3 are compared.Table 2 lists the median (p 50 ) and 95%-percentile (p 95 ) values of the errors for the different configurations.It is observed that for a normal office placement (see Figure 4a), median errors are around 2-3 cm and 95-percentile errors between 6 and 10 cm.The LSE metric slightly outperforms the nLSE metric, whereas traditional trilateration produces errors that are almost 40% higher on average than for LSE.For a sloppy industrial placement (see Figure 4b), median (p 50 ) and maximal (p 95 ) errors more or less double, compared to the normal placement, for all metrics.
When increasing the height to 6 m (industrial), errors increase, but not significantly.For a normal industrial placement (see Figure 4c), although median errors increase between 8 (trilateration) and 24% (LSE) compared to the normal office placement, maximal errors even see a slight decrease for the LSE and trilateration approaches.The next section will elaborate on this impact of LED height more thoroughly.Finally, for a sloppy industrial placement (see Figure 4d), the errors again more or less double compared to the normal industrial placement, suggesting that doubling σ tilt for each of the LEDs also doubles the resulting positioning error.
To give more insight on the spatial distribution of the postion estimates for the different localization approaches, Figure 5a-c show a scatter plot of the estimated positions for a LED-PD height difference of 2.5 m, and a σ tilt value of 2°, for LSE, nLSE, and trilateration, respectively.For each of the 100 evaluation points, 250 estimations are displayed.These 100 evaluation points correspond to the blue dots in Figure 3.The square that is formed by the LEDs is partly shown with orange lines.These scatter plots correspond to the values in Table 2 for a 'sloppy office' deployment.The dominant LED at (1.5, 1.5) causes a more circular pattern due its LED tilt for the LSE metric than for nLSE, showing that nLSE indeed reduces the impact of the dominant LED.Trilateration performs worse than the other two metrics, but especially outside the LED square (i.e., for x or y coordinates <1.5 m).Inside the LED square, the performance of the algorithms is comparable.

LED Tilt Impact for Different Inter-LED Distances, LED Height and Tilt
In this section, we further characterize the errors due to random LED tilt, using the common trilateration method described in Section 2.3.1.Different configurations with four LEDs in a square configuration are investigated.We vary the side of the square (i.e., the inter-LED distance) between 2 and 8 m in steps of 1 m, the LED height between 2 and 8 m in steps of 1 m, and σ tilt between 1 and 3°in steps of 1°.For each of these 7 × 7 × 3 =147 combinations, each of the 4 LEDs are randomly tilted and the localisation accuracy is evaluated on a uniform 11 × 11 grid in between the LED locations, meaning that the accuracy is evaluated only inside the LED square, for which was observed that all three localization approaches yielded comparable estimation clouds (see Figure 5).These 121 evaluation points correspond to the orange dots in Figure 3. Please note that Figure 3 represents a configuration with square side length S = 4.For each of the 147 combinations, 5000 LED tilt settings are generated, with the LEDs tilted according to the statistical tilt distributions presented in Section 2.2.As (thermal or shot) noise is not considered here, the value of the transmit power P E or the PD area A R has no influence on the positioning accuracy, so conclusions are also valid for other LED optical powers or PD areas.
Figure 6a shows the p 50 error due to LED tilt for 7 different square side lengths (2 to 8 m), as a function of the LED height, for a σ tilt value of 1°.The figure shows that irrespective of the LED height, the impact of tilt is the lowest for smaller LED squares.e.g., for a LED height of 2 m, the median induced positioning error increases from 1 cm for a 2 × 2 m LED square to 6.5 cm for an 8 × 8 m LED square.As the LED height increases, the influence of the square side length reduces, with median errors between 2.4 cm (2 × 2 m square) and 3.8 cm (8 × 8 m square).Furthermore, it is interesting to note that the impact of LED tilt on positioning accuracy in the receiver plane does not necessarily increase with LED height (i.e., higher LED-PD height differences), while intuitively (and from Figure 1a), one might expect a monotonic increase of the induced error with LED height.However, for a given LED square size S, there appears to be an optimal LED height with respect to the impact of LED tilt.e.g., for an 8 × 8 m LED square, this optimal LED height is 6 m, decreasing to 4 m for a 5 × 5 m LED square, and to smaller than 2 m for a 2 × 2 m LED square.
We explain this using Figure 1b, where a 1D-simplified side view of the considered configuration is shown, with h the LED-PD height difference and S the LED square side length.With γ between 0 and 1, γS indicates the range of possible lateral displacements L = h • tanφ of the PD in the receiver plane, with respect to perpendicular irradiance (γ = 0 is below the LED, γ = 1 is below an adjacent LED of the LED square).φ is the angle of irradiance for an untilted LED.The variable of interest here is ∂L ∂φ , indicating the additional lateral displacement in the receiver plane, due to a LED tilt h .The result comprises two opposed phenomena.One the one hand, we observe that larger LED heights h correspond to larger additional displacements in the receiver plane.This corresponds to the phenomenon described in Section 1.2 and Figure 1, and is is reflected by the first factor h of ∂L ∂φ , meaning that the tilt impact is smaller for smaller LED heights.On the other hand, we see that for a given (and fixed) LED square size S, relatively more receiver locations will have a large angle of irradiance φ when h is smaller (see also Figure 1b), leading to a larger additional lateral displacement.This is reflected by the second factor γ 2 S 2 h of ∂L ∂φ .The combination of these two opposed phenomena causes, depending on the size of the LED square S, the minimal tilt impact to be observed at different height differences h.When the LED-PD height difference is equal to the inter-LED distance (h = S), forming a cube, the error sees a perfect linear increase with h, as indicated by the red curve in Figure 6a, again corresponding to the intuitive assumption described in Section 1.2 ( ∂L ∂φ reduces to h(1 + γ 2 ) for this simplified 1D-case).
Finally, it was again observed that p 95 errors scale approximately to values equal to three times the p 50 values.Similarly, the p 50 errors scale also linearly with the σ tilt value, as is shown in Figure 6b, depicting the induced median error for a 4 × 4 m LED square, as a function of LED height for three σ tilt values.

Evaluation of LED Tilt Estimation Methods
In this section, the performance of the LED tilt estimation methods presented in Section 2.4 are presented.Figure 7a,b show the P R • d 3 values in the receiver plane without receiver noise and with added noise to the received powers (σ noise = 10 −7 W) respectively ('exhaustive search method').The figures indicate the LED location (green asterisk in (0,0)), the real (black asterisk) and estimated (red asterisk) intersect of the LED normal with the receiver plane, and the location receiving the maximal power P R (blue asterisk).It should be noted that the direction of the tilted LED normal is indeed not only determined by the maximal received power P R (blue asterisk), as this maximal-power location is determined by a tradeoff between being located along the LED normal (pulling towards the black asterisk) and having the shortest distance to the LED (pulling the location back towards right underneath the LED, green asterisk).This is also why the tilted LED normal intersect is not found at the location with a maximal P R , but instead at the location with a maximal P R • d 3 .Figure 7a shows that under the absence of noise, the LED tilt can be estimated exactly (red asterisk co-located with black asterisk), indicating the correctness of the proposed method.When noise is added, θ and α are estimated at 2.1°and 201°respectively (compared to the real values θ = 1.6°and α = 230°).To reduce or avoid the possible effect of obtaining an erroneous noise-induced maximum of P R • d 3 , an alternative approach could be to calculate the centre of the set of the X locations with the highest measured P R • d 3 values.The value of X should be determined based on a tradeoff between a sufficient tilt estimation precision and a sufficient spatial averaging of possible noise.16) without noise and with added noise to the considered measurement points (σ noise = 10 −7 W). Figure 7c indicates the correct operation of the quick search method (green asterisk co-located with red asterisk).Figure 7d shows a slight deviation of the estimation due to the added noise (error in θ and α of 0.1°and 4.5°respectively).Comparison of the LED tilt estimations of the two methods under the presence of noise, shows that the quick search using only four measurement points typically performs better than the exhaustive search, despite the fact that fewer measurements were used.This is explained by the fact that the intersects of the tilted LED normals are usually very close to the untilted LED normal (a few cm, see Figure 7a,b).Due to the low gradient of the received power in the area underneath the LED (m = 1), differences in P R (or in P R • d 3 ) are very small, making them susceptible to noise.The quick search method uses measurement points that are further from the intersect of the untilted LED normal with the receiver plane.Thanks to the larger power gradient there, this method appears to be less susceptible to noise.Or differently stated, LED tilt is better noticed at larger angles of irradiance, the phenomenon also discussed in the previous section and noticed in Figure 1b.

Discussion
In this paper, we investigated to what extent the uncertainty on the LED tilt impacts RSS-based VLP accuracies.For a 7 m × 7 m room with four LEDs placed in a 4 m × 4 m square, positioning errors at 100 locations in the receiver plane were compared for three localization approaches, based on a Monte Carlo simulation consisting of 10 4 simulations.It was shown the model-based fingerprinting methods (Least-Squares Estimation and normalized Least-Squares Estimation) performed slightly better than a traditional trilateration, mainly for locations outside the LED square.We observed that tilt-induced errors are in the order of centimeters, depending on the LED configuration (LED height and inter-LED distance), and the severity of the tilt.The errors scale linearly with the severity of the LED tilt.
We found that for a σ tilt value of 1°, median errors for a LED height of 7 m were between 2 and 4 cm (see Figure 6a, depending on inter-LED distance), which indicates a more limited impact than intuitively assumed from Figure 1a, where a lateral deviation of 12.2 cm was found.This indicates that the lateral deviation of the LED normal in the receiver plane corresponds to a significant overestimation of tilt-induced positioning errors, suggesting that induced errors of multiple LEDs compensate each other.It has also been shown that increasing the LED height does not necessarily increase the tilt-induced errors: depending on the inter-LED distance, there is an LED height optimum where LED tilt impact is minimal.
We now compare these errors with errors induced by noise, and errors due to deviations on the tabulated optical LED power.In [6], the performance of VLP under the presence of reflections was investigated for h = 1.65 m and S = 2.5 m and a wall reflectance factor of 0.3.Median errors between 6.7 and 8.7 cm were found, depending on the metric used.Based on the curves presented in Figure 6a, a median error between 1 and 1.5 cm due to tilt can be expected for a configuration with these h and S values, when σ tilt = 1°.In [17], it was investigated to what extent deviations from the tabulated transmitted optical LED power impacted positioning accuracy.Again, for h = 1.65 m and S = 2.5 m, median and maximal (95% percentile) errors of 4.21 and 10.75 cm were found for a LED tolerance (i.e., 3-sigma value) of 10%.For the same configuration, median errors ranging from about 1 mm to 10 cm were obtained for standard deviations of the observed noise power ranging between 10 −8 and 10 −6 W [17].
Knowing that the tilt-induced errors investigated here, add up with those introduced by thermal noise and shot noise, reflections, receiver tilt, and imperfections in the LED radiation pattern, LED tilt is one of the crucial aspects to consider and compensate for, since the other aspects are often harder to compensate for: receiver tilt might be variable while moving, noise is a random process, and reflections and LED radiations patterns are not easy to model well, or require the execution of a large measurement campaign.Therefore, estimating LED tilt and compensating for it might consistently reduce VLP errors by a few centimeters.As a first exploration, two methods to estimate the LED tilt have been presented.Although intuitively it could be assumed that an exhaustive analysis of the (maximal) powers in the area right underneath the tilted LED would be the best approach to accurately estimate LED tilts, it seems that a limited set of measurements at different, more distant points would deliver better estimates.It remains to be (experimentally) investigated how many points would ideally be required and which locations would be optimal.This is considered as future work.The presented method shows to have the potential to estimate LED tilt, which in turn allows adjusting model-based fingerprinting maps for an improved RSS-based VLP performance.Another interesting future research track is a full sensitivity analysis, including the impact of noise, LED tilt, PD tilt, and wall reflections.

Figure 1 .
Lateral displacements in receiver plane in case of LED tilt.(a) Maximal lateral displacement of LED normal in the receiver plane due to LED tilt, as a function of LED-PD height difference; (b) Additional displacement ∂L in receiver plane due to LED tilt ∂φ as a function of angle of irradiance φ (S = square side length, h = LED-PD height difference).

Figure 2 .
Figure 2. Overview of visible light channel.

Figure 3 .
Figure 3. Simulation configuration (A,B,C,D indicate LED locations, blue dots indicate scatter plot ground truth locations, see Section 3.1, orange dots indicate the evaluation points used in Section 3.2).

Figure 5 .
Figure 5. Scatter plot of 250 location estimations for each of 100 locations (indicated by white-edged dots) in top left part of the 7 m by 7 m square of Figure 3, for a LED height of 2.5 m and a σ tilt = 2°.The LED square is indicated with an orange line.

Figure 6 .
Induced median error due to LED tilt for different LED square sizes, heights, and tilts.(a) Induced median error for square side lengths between 2 and 8 m, as a function of LED height, for a σ tilt value of 1°; (b) Induced median error for a 4 × 4 m LED square, as a function of LED height for three σ tilt values.

Figure
Figure 7c,d show the cost function from Equation (16) without noise and with added noise to the considered measurement points (σ noise = 10 −7 W).Figure7cindicates the correct operation of the quick search method (green asterisk co-located with red asterisk).Figure7dshows a slight deviation of the estimation due to the added noise (error in θ and α of 0.1°and 4.5°respectively).Comparison of the LED tilt estimations of the two methods under the presence of noise, shows that the quick search using

Table 1 .
Summary of the parameters defined in Figure2.

Table 2 .
Median and maximal positioning errors for the three localization approaches, for four typical LED deployments.