Evaluation of Deployment Challenges of Wireless Sensor Networks at Signalized Intersections

With the growing demand of Intelligent Transportation Systems (ITS) for safer and more efficient transportation, research on and development of such vehicular communication systems have increased considerably in the last years. The use of wireless networks in vehicular environments has grown exponentially. However, it is highly important to analyze radio propagation prior to the deployment of a wireless sensor network in such complex scenarios. In this work, the radio wave characterization for ISM 2.4 GHz and 5 GHz Wireless Sensor Networks (WSNs) deployed taking advantage of the existence of traffic light infrastructure has been assessed. By means of an in-house developed 3D ray launching algorithm, the impact of topology as well as urban morphology of the environment has been analyzed, emulating the realistic operation in the framework of the scenario. The complexity of the scenario, which is an intersection city area with traffic lights, vehicles, people, buildings, vegetation and urban environment, makes necessary the channel characterization with accurate models before the deployment of wireless networks. A measurement campaign has been conducted emulating the interaction of the system, in the vicinity of pedestrians as well as nearby vehicles. A real time interactive application has been developed and tested in order to visualize and monitor traffic as well as pedestrian user location and behavior. Results show that the use of deterministic tools in WSN deployment can aid in providing optimal layouts in terms of coverage, capacity and energy efficiency of the network.


Introduction
In the past two decades, the number of vehicles circulating in cities has increased considerably. This fact, together with the tendency of citizens to cluster has caused congestion of roads and junctions, especially in cities with large populations. The existing infrastructures have become inefficient to optimize traffic. The traditional traffic control system based on traffic lights depending on preset date and time patterns is not efficient in reducing waiting times and energy consumption. In fact, simulations show that adaptive traffic control based on information from Wireless Sensor Networks (WSNs) can be improved by 65% in normal traffic conditions compared to traffic management by preset patterns [1]. Since the nineties, there has been extensive research looking for different solutions in Intelligent Transportation Systems (ITS), with the objective of improving and performing traffic simulation procedure for the channel characterization of the considered scenarios and the scenarios description. Section 3 shows simulation results such as bi-dimensional planes of received power, power delay profiles and delay spread, among others. Section 4 presents the experimental-setup in the real scenario and the comparison between simulation and measurements. Section 5 describes the radio planning analysis for different technologies (ZigBee, Bluetooth and the 802.11p radio standard). Section 6 describes the proposed application within the considered scenario, followed by the last section presenting concluding remarks.

Ray Launching Technique
A deterministic method based on an in-house developed 3-D RL code has been used to analyze the radio electric behavior of the considered scenario. The 3D RL algorithm is based on Geometrical Optics (GO) and Geometrical Theory of Diffraction (GTD). It has been verified in the literature with different applications, such as the analysis of wireless propagation in closed environments [27][28][29][30], interference analysis [31] or electromagnetic dosimetry evaluation in wireless systems [32]. The main principle of the RL techniques is to identify a single point on the wave front of the radiated wave with a ray that propagates in the space following a combination of optic and electromagnetic theories. This phenomenon is illustrated in Figure 1 for the considered scenario. Each ray propagates in the space as a single optic ray. The incident electric field (E i ) created by an antenna at a distance r in the free space can be calculated by [33]: where P rad is the radiated power with a directivity D t pθ t , φ t q, where the sub-index t refers to the transmitted angle, and polarization ratio´X K , X ¯. β 0 " 2π f c ? ε 0 µ 0 , ε 0 " 8.854ˆ10´1 2 , µ 0 " 4πˆ10´7 and η 0 " 120π. L K are the path loss coefficients for each polarization. The parameter j in Equation (1) refers to the complex number.
Sensors 2016, 16, 1140 4 of 25 the real scenario and the comparison between simulation and measurements. Section 5 describes the radio planning analysis for different technologies (ZigBee, Bluetooth and the 802.11p radio standard). Section 6 describes the proposed application within the considered scenario, followed by the last section presenting concluding remarks.

Ray Launching Technique
A deterministic method based on an in-house developed 3-D RL code has been used to analyze the radio electric behavior of the considered scenario. The 3D RL algorithm is based on Geometrical Optics (GO) and Geometrical Theory of Diffraction (GTD). It has been verified in the literature with different applications, such as the analysis of wireless propagation in closed environments [27][28][29][30], interference analysis [31] or electromagnetic dosimetry evaluation in wireless systems [32]. The main principle of the RL techniques is to identify a single point on the wave front of the radiated wave with a ray that propagates in the space following a combination of optic and electromagnetic theories. This phenomenon is illustrated in Figure 1 for the considered scenario. Each ray propagates in the space as a single optic ray. The incident electric field (Ei) created by an antenna at a distance r in the free space can be calculated by [33]: where is the radiated power with a directivity ( , ), where the sub-index t refers to the transmitted angle, and polarization ratio ( , ∥ ) .  When the rays impinge with an obstacle in their path, a reflected and a transmitted ray are created with new angles provided by Snell's law [34]. It is important to take into account that the rays considered in the GO only approach are only direct, reflected and refracted rays, leading to the existence of abrupt areas, which correspond to the boundaries of the regions where these rays exist. Because of that, the diffracted rays are introduced with the GTD and its uniform extension, the Uniform GTD (UTD). The main purpose of these diffracted rays is to remove field discontinuities and to introduce proper field corrections, especially in the zero-field regions predicted by GO. The diffracted field is calculated by [35]: When the rays impinge with an obstacle in their path, a reflected and a transmitted ray are created with new angles provided by Snell's law [34]. It is important to take into account that the rays considered in the GO only approach are only direct, reflected and refracted rays, leading to the existence of abrupt areas, which correspond to the boundaries of the regions where these rays exist. Because of Sensors 2016, 16, 1140 5 of 25 that, the diffracted rays are introduced with the GTD and its uniform extension, the Uniform GTD (UTD). The main purpose of these diffracted rays is to remove field discontinuities and to introduce proper field corrections, especially in the zero-field regions predicted by GO. The diffracted field is calculated by [35]: where s 1 , s 2 are the distances from the source to the edge and from the edge to the receiver point, respectively. D K are the diffraction coefficients given by [35,36]. The RL algorithm is performed three-dimensionally, with angular resolution (horizontal and vertical planes) in a predefined solid angle that considers the radiation diagram of the transceivers sources. Spatial resolution is also defined by a uniform hexahedral mesh. Parameters such as frequency of operation, radiation patterns of the antennas, number of multipath reflections, separation angle between rays, and cuboid dimension can be taken into account. Besides, all the material properties for all the elements within the scenario can also be considered, given the dielectric constant and the loss tangent at the frequency range of operation of the system under analysis. When a ray impacts with an obstacle, reflection, refraction and diffraction will occur, depending on the geometry and the electric properties of the object, as is depicted in Figure 2.
where , are the distances from the source to the edge and from the edge to the receiver point, respectively.
The RL algorithm is performed three-dimensionally, with angular resolution (horizontal and vertical planes) in a predefined solid angle that considers the radiation diagram of the transceivers sources. Spatial resolution is also defined by a uniform hexahedral mesh. Parameters such as frequency of operation, radiation patterns of the antennas, number of multipath reflections, separation angle between rays, and cuboid dimension can be taken into account. Besides, all the material properties for all the elements within the scenario can also be considered, given the dielectric constant and the loss tangent at the frequency range of operation of the system under analysis. When a ray impacts with an obstacle, reflection, refraction and diffraction will occur, depending on the geometry and the electric properties of the object, as is depicted in Figure 2.

Simulation Scenario Description
The simulation scenario implemented for calculation by means of the in-house developed 3D RL code consists of an outdoor environment of an urban city area, placed at the old town of Pamplona (Navarre, Spain). The considered scenario consists of an outdoor environment with trees of different heights, grass, junctions of roads, cars, people and different types of buildings. This environment can be considered as really complex, in terms of radio wave propagation, because it has a broad range of elements with different material properties within it. The dimensions of the scenario are (150 m, 130 m, 30 m). Figure 3 presents a view of the real and schematic considered scenarios.
All the material properties for all the elements within the scenarios have been considered, given the dielectric constant and the loss tangent at the frequency range of operation of the system under analysis. Taking into account that radio wave propagation in this scenario could change depending of the weather, as it is an outdoor environment, it is relevant to consider different conditions for the material properties of the vegetation. For that purpose, the values obtained in [37] for the material properties of the wood and the foliage of the vegetation have been used. A human body model [32], which take into account all the elements within the human body, and a car model [38], which also consider all the elements within a car, has also been used in the simulations. These models have been created to be implemented in the in-house 3D RL algorithm. The material parameters used in the simulation are defined in Table 1 [37][38][39][40].

Simulation Scenario Description
The simulation scenario implemented for calculation by means of the in-house developed 3D RL code consists of an outdoor environment of an urban city area, placed at the old town of Pamplona (Navarre, Spain). The considered scenario consists of an outdoor environment with trees of different heights, grass, junctions of roads, cars, people and different types of buildings. This environment can be considered as really complex, in terms of radio wave propagation, because it has a broad range of elements with different material properties within it. The dimensions of the scenario are (150 m, 130 m, 30 m). Figure 3 presents a view of the real and schematic considered scenarios.
All the material properties for all the elements within the scenarios have been considered, given the dielectric constant and the loss tangent at the frequency range of operation of the system under analysis. Taking into account that radio wave propagation in this scenario could change depending of the weather, as it is an outdoor environment, it is relevant to consider different conditions for the material properties of the vegetation. For that purpose, the values obtained in [37] for the material properties of the wood and the foliage of the vegetation have been used. A human body model [32], which take into account all the elements within the human body, and a car model [38], which also consider all the elements within a car, has also been used in the simulations. These models have been created to be Sensors 2016, 16, 1140 6 of 25 implemented in the in-house 3D RL algorithm. The material parameters used in the simulation are defined in Table 1 [37][38][39][40].

Simulation Results
Once the simulation scenario has been defined, simulation results can be obtained. These results have been validated with real measurements which have been made in a campaign of measurements in the same scenario, which is shown in Section 4.
For the simulations, 16 antennas have been placed in the considered scenario, one for each traffic light. The specific position of the antennas in the scenario is shown in Figure 4. To reduce interferences with people and vehicles, each antenna has been placed 10 cm above each traffic light. The radiating element is a wireless ZigBee mote which has been configured as a dipole, transmitting 18 dBm at 2.41 GHz. Simulation parameters are shown in Table 2.

Simulation Results
Once the simulation scenario has been defined, simulation results can be obtained. These results have been validated with real measurements which have been made in a campaign of measurements in the same scenario, which is shown in Section 4.
For the simulations, 16 antennas have been placed in the considered scenario, one for each traffic light. The specific position of the antennas in the scenario is shown in Figure 4. To reduce interferences with people and vehicles, each antenna has been placed 10 cm above each traffic light. The radiating element is a wireless ZigBee mote which has been configured as a dipole, transmitting 18 dBm at 2.41 GHz. Simulation parameters are shown in Table 2.  Figures 5 and 6 shows the power distribution within the considered scenario for different heights and planes, and for different sensors placed in different traffic lights. As it can be seen, the influence of the obstacles (like the trees, streetlights, buildings, etc.) can be easily appreciated. It is shown that the morphology as well as the topology of the considered scenario has a noticeable impact on radio wave propagation.  Figures 5 and 6 shows the power distribution within the considered scenario for different heights and planes, and for different sensors placed in different traffic lights. As it can be seen, the influence of the obstacles (like the trees, streetlights, buildings, etc.) can be easily appreciated. It is shown that the morphology as well as the topology of the considered scenario has a noticeable impact on radio wave propagation.  Figures 5 and 6 shows the power distribution within the considered scenario for different heights and planes, and for different sensors placed in different traffic lights. As it can be seen, the influence of the obstacles (like the trees, streetlights, buildings, etc.) can be easily appreciated. It is shown that the morphology as well as the topology of the considered scenario has a noticeable impact on radio wave propagation.    Figure 7 shows the different zones of estimated received power extracted from the data of Figure 5a,c. It can be seen that received power results are accurate but not homogeneous, and because of that, the radio planning task has an important role. Different values of the mean received power are obtained depending of the obstacle density of the different zones. Table 3 show the different values of power mean with its standard deviation and the obstacle density of the zone considered. It can be seen that the obstacle density has great influence in received power, but it is also relevant to consider the area of the considered zone and the distance from the transmitter of each zone. From these results, it is shown that received power in different points of the scenarios depends greatly of the position of the transmitter, the obstacle density of the different zones and the different material properties of the obstacles.
An important radioelectric phenomenon in this type of environment is given by multipath propagation. To illustrate this fact, the delay spread for the whole scenario and the power delay profile for the central location of the scenario has been obtained and it is shown in Figures 8 and 9, respectively. The delay spread has been defined as the time delay between the first and the last ray which arrives to each spatial point. It has been calculated using as threshold the noise floor, with a value of −100 dBm.   Figure 5a,c. It can be seen that received power results are accurate but not homogeneous, and because of that, the radio planning task has an important role. Different values of the mean received power are obtained depending of the obstacle density of the different zones. Table 3 show the different values of power mean with its standard deviation and the obstacle density of the zone considered. It can be seen that the obstacle density has great influence in received power, but it is also relevant to consider the area of the considered zone and the distance from the transmitter of each zone. From these results, it is shown that received power in different points of the scenarios depends greatly of the position of the transmitter, the obstacle density of the different zones and the different material properties of the obstacles.
An important radioelectric phenomenon in this type of environment is given by multipath propagation. To illustrate this fact, the delay spread for the whole scenario and the power delay profile for the central location of the scenario has been obtained and it is shown in Figures 8 and 9, respectively. The delay spread has been defined as the time delay between the first and the last ray which arrives to each spatial point. It has been calculated using as threshold the noise floor, with a value of´100 dBm.  Figure 7 shows the different zones of estimated received power extracted from the data of Figure 5a,c. It can be seen that received power results are accurate but not homogeneous, and because of that, the radio planning task has an important role. Different values of the mean received power are obtained depending of the obstacle density of the different zones. Table 3 show the different values of power mean with its standard deviation and the obstacle density of the zone considered. It can be seen that the obstacle density has great influence in received power, but it is also relevant to consider the area of the considered zone and the distance from the transmitter of each zone. From these results, it is shown that received power in different points of the scenarios depends greatly of the position of the transmitter, the obstacle density of the different zones and the different material properties of the obstacles.
An important radioelectric phenomenon in this type of environment is given by multipath propagation. To illustrate this fact, the delay spread for the whole scenario and the power delay profile for the central location of the scenario has been obtained and it is shown in Figures 8 and 9, respectively. The delay spread has been defined as the time delay between the first and the last ray which arrives to each spatial point. It has been calculated using as threshold the noise floor, with a value of −100 dBm.  From Figure 8 it can be seen that multipath phenomena have a lot of influence in the environment. It is observed that the delay spread is higher in the areas closest to the transmitter antenna and it is lower in other areas because the radiation power in these parts is smaller due to the larger distance. Each spatial sample of the delay spread is associated with a power delay profile, shown in Figure 9 for the central position of the scenario. As it is observed, the scenario is really complex and there are several echoes in the scenario due to this multipath channel behavior.  From Figure 8 it can be seen that multipath phenomena have a lot of influence in the environment. It is observed that the delay spread is higher in the areas closest to the transmitter antenna and it is lower in other areas because the radiation power in these parts is smaller due to the larger distance. Each spatial sample of the delay spread is associated with a power delay profile, shown in Figure 9 for the central position of the scenario. As it is observed, the scenario is really complex and there are several echoes in the scenario due to this multipath channel behavior.  As previously shown, a clear dependence on the topology and morphology of the scenario as well as on the system parameters is observed and must be carefully considered in order to assess the performance of multiple systems within a particular scenario. These results allow the optimization of the design of the WSN since the designer can select the minimum number of nodes required to grant a certain communication level, and the optimal emplacements of these nodes. It is important to consider also that a larger density of nodes can lead to increased interference levels, which could degrade system performance. This comment also applies in the case of changes in the transmitted service, associated with a certain allocated bandwidth and hence, changes in sensitivity levels.
The computational time is also highly important when analyzing these complex scenarios. In this case, simulations have been performed in an Intel Xeon CPU X5650 @ 6.67 GHz and 2.66 GHz, and the simulation computational time has been of 43,537 s using the software Matlab from Mathworks (Natick, MA, USA). It must be pointed out that the input parameters have been chosen according to a convergence analysis of the algorithm previously done for different types of scenarios [41].
A key parameter that is used in assessing systems that transmit digital data from one location to another is the Bit Error Rate (BER). The BER expression for QPSK modulation can be calculated by:  As previously shown, a clear dependence on the topology and morphology of the scenario as well as on the system parameters is observed and must be carefully considered in order to assess the performance of multiple systems within a particular scenario. These results allow the optimization of the design of the WSN since the designer can select the minimum number of nodes required to grant a certain communication level, and the optimal emplacements of these nodes. It is important to consider also that a larger density of nodes can lead to increased interference levels, which could degrade system performance. This comment also applies in the case of changes in the transmitted service, associated with a certain allocated bandwidth and hence, changes in sensitivity levels.
The computational time is also highly important when analyzing these complex scenarios. In this case, simulations have been performed in an Intel Xeon CPU X5650 @ 6.67 GHz and 2.66 GHz, and the simulation computational time has been of 43,537 s using the software Matlab from Mathworks (Natick, MA, USA). It must be pointed out that the input parameters have been chosen according to a convergence analysis of the algorithm previously done for different types of scenarios [41].
A key parameter that is used in assessing systems that transmit digital data from one location to another is the Bit Error Rate (BER). The BER expression for QPSK modulation can be calculated by: As previously shown, a clear dependence on the topology and morphology of the scenario as well as on the system parameters is observed and must be carefully considered in order to assess the performance of multiple systems within a particular scenario. These results allow the optimization of the design of the WSN since the designer can select the minimum number of nodes required to grant a certain communication level, and the optimal emplacements of these nodes. It is important to consider also that a larger density of nodes can lead to increased interference levels, which could degrade system performance. This comment also applies in the case of changes in the transmitted service, associated with a certain allocated bandwidth and hence, changes in sensitivity levels.
The computational time is also highly important when analyzing these complex scenarios. In this case, simulations have been performed in an Intel Xeon CPU X5650 @ 6.67 GHz and 2.66 GHz, and the simulation computational time has been of 43,537 s using the software Matlab from Mathworks (Natick, MA, USA). It must be pointed out that the input parameters have been chosen according to a convergence analysis of the algorithm previously done for different types of scenarios [41].
A key parameter that is used in assessing systems that transmit digital data from one location to another is the Bit Error Rate (BER). The BER expression for QPSK modulation can be calculated by: where E b " P RX {R b and N 0 is the noise value. The received power P RX has been calculated with the 3D RL algorithm for each spatial sample in the considered scenario. With these values of P RX , the BER has been calculated and it is shown in Figures 10 and 11, for different values of data rate (R b ) and different values of N 0 . It can be seen the high variability between the different cases, with higher values of BER when N 0 is higher. In addition, differences between the different data rates considered are observed, leading to a lower BER for all cases with the lowest data rate considered. These results can be really helpful in order to optimize the design and deployment of a WSN depending on the modulation, the used data rate and the level of N 0 . where = ⁄ and N0 is the noise value. The received power has been calculated with the 3D RL algorithm for each spatial sample in the considered scenario. With these values of , the BER has been calculated and it is shown in Figures 10 and 11, for different values of data rate ( ) and different values of . It can be seen the high variability between the different cases, with higher values of BER when is higher. In addition, differences between the different data rates considered are observed, leading to a lower BER for all cases with the lowest data rate considered. These results can be really helpful in order to optimize the design and deployment of a WSN depending on the modulation, the used data rate and the level of .

Experimental Setup
In order to validate the results previously obtained, real measurements in the same scenario shown in Figure 3 have been performed. A transmitter antenna, connected to a signal generator at 2.41 GHz and 5.9 GHz has been located at the coordinates (X = 31 m, Y = 76.8 m, Z = 2.05 m), just above the traffic light #1. The employed signal generator is a portable N1996A (Agilent Technologies, Santa Clara, CA, USA) unit and the spectrum analyzer is an Agilent N9912 Field Fox. The antennas used are ACA-4HSRPP-2458 from Zentri (Los Gatos, CA, USA), both omnidirectional antennas, with a gain of 0.3 dB for 2.41 GHz and a gain of 3.74 dB for 5.9 GHz. Measurements have been performed along the measurement points depicted in Figure 12   where = ⁄ and N0 is the noise value. The received power has been calculated with the 3D RL algorithm for each spatial sample in the considered scenario. With these values of , the BER has been calculated and it is shown in Figures 10 and 11, for different values of data rate ( ) and different values of . It can be seen the high variability between the different cases, with higher values of BER when is higher. In addition, differences between the different data rates considered are observed, leading to a lower BER for all cases with the lowest data rate considered. These results can be really helpful in order to optimize the design and deployment of a WSN depending on the modulation, the used data rate and the level of .

Experimental Setup
In order to validate the results previously obtained, real measurements in the same scenario shown in Figure 3

Experimental Setup
In order to validate the results previously obtained, real measurements in the same scenario shown in Figure 3  To visualize the frequency variation with time and the possible influence of pre-existent interference in the scenario, a spectrogram has been measured before starting the measurement campaign. Figure 13 shows the measured spectrogram for both frequencies in Max Hold mode, for the transmitter position represented in Figure 12 in the considered scenario, with the aid of the Agilent N9912 Field Fox portable spectrum analyzer. It can be seen that the measured values can be considered as noise and there is not a significant pre-existent interference at the frequency bands of interest. It is also observed that the noise is To visualize the frequency variation with time and the possible influence of pre-existent interference in the scenario, a spectrogram has been measured before starting the measurement campaign. Figure 13 shows the measured spectrogram for both frequencies in Max Hold mode, for the transmitter position represented in Figure 12 in the considered scenario, with the aid of the Agilent N9912 Field Fox portable spectrum analyzer. To visualize the frequency variation with time and the possible influence of pre-existent interference in the scenario, a spectrogram has been measured before starting the measurement campaign. Figure 13 shows the measured spectrogram for both frequencies in Max Hold mode, for the transmitter position represented in Figure 12 in the considered scenario, with the aid of the Agilent N9912 Field Fox portable spectrum analyzer. It can be seen that the measured values can be considered as noise and there is not a significant pre-existent interference at the frequency bands of interest. It is also observed that the noise is It can be seen that the measured values can be considered as noise and there is not a significant pre-existent interference at the frequency bands of interest. It is also observed that the noise is slightly higher for 5.9 GHz frequency. This is due to the fact the gain of the antenna in 5.9 GHz is 3.44 dB higher than at 2.41 GHz frequency. Figure 14 shows the comparison between simulation and measurement results for both frequencies. Measurements were performed with 100 MHz bandwidth at 2.41 GHz and 5.9 GHz frequency. The measurement time at each point was 60 s, and the value of received power represented by each point is the highest peak of power shown by the spectrum analyzer for the considered bandwidth (MaxHold function in the spectrum analyzer of Agilent). The received power values estimated by simulation have been obtained for the same spatial samples as the real measurements, considering the corresponding cuboid in the three-dimensional mesh of cuboids in which the scenario have been divided. It is observed a good agreement between simulation and measurement results, with a mean error of 0.27 dBm with a standard deviation of 1.13 dB for 2.41 GHz frequency, and 0.66 dBm with a standard deviation of 1.46 dB for 5.9 GHz frequency. The resulting error means are very low, indicating that the proposed 3D RL simulation method works properly, validating in the same way the simulation results shown in the previous sections of this work. slightly higher for 5.9 GHz frequency. This is due to the fact the gain of the antenna in 5.9 GHz is 3.44 dB higher than at 2.41 GHz frequency. Figure 14 shows the comparison between simulation and measurement results for both frequencies. Measurements were performed with 100 MHz bandwidth at 2.41 GHz and 5.9 GHz frequency. The measurement time at each point was 60 s, and the value of received power represented by each point is the highest peak of power shown by the spectrum analyzer for the considered bandwidth (MaxHold function in the spectrum analyzer of Agilent). The received power values estimated by simulation have been obtained for the same spatial samples as the real measurements, considering the corresponding cuboid in the three-dimensional mesh of cuboids in which the scenario have been divided. It is observed a good agreement between simulation and measurement results, with a mean error of 0.27 dBm with a standard deviation of 1.13 dB for 2.41 GHz frequency, and 0.66 dBm with a standard deviation of 1.46 dB for 5.9 GHz frequency. The resulting error means are very low, indicating that the proposed 3D RL simulation method works properly, validating in the same way the simulation results shown in the previous sections of this work.

Radioplanning Analysis
The deployment of three different communication systems has been analyzed within the considered scenario. Firstly, an analysis of the received power and the receiver's sensitivity for the different technologies has been analyzed, leading to the most efficient coverage ratio for the deployment of the different systems. The first technology evaluated has been ZigBee technology. Specifically, the wireless devices used for simulation have been both, the XBee Pro and XBee motes from Digi International Inc. (Minnetonka, MN, USA). The main difference between them is the larger transmitter power of the XBee Pro, leading to longer range distance. In order to analyze the worst case conditions, the transmitted power levels considered has been reduced to the minimum default value. After that, the second technology evaluated has been Bluetooth, specifically Bluetooth Low Energy (BLE) and classic Bluetooth. BLE has been designed as a low-power solution for control and monitoring applications, in contrast with classic Bluetooth. Finally, the third technology evaluated has been the 802.11p radio, which is the standard to add wireless access in vehicular environments (WAVE). It must be pointed out that for this technology we have performed new simulations in the 5.9 GHz frequency. Three different modulations with three different data rates have been analyzed for the 802.11p radio standard. The transmitted power and sensitivity for the different communication systems considered is shown in Table 4.

Radioplanning Analysis
The deployment of three different communication systems has been analyzed within the considered scenario. Firstly, an analysis of the received power and the receiver's sensitivity for the different technologies has been analyzed, leading to the most efficient coverage ratio for the deployment of the different systems. The first technology evaluated has been ZigBee technology. Specifically, the wireless devices used for simulation have been both, the XBee Pro and XBee motes from Digi International Inc. (Minnetonka, MN, USA). The main difference between them is the larger transmitter power of the XBee Pro, leading to longer range distance. In order to analyze the worst case conditions, the transmitted power levels considered has been reduced to the minimum default value. After that, the second technology evaluated has been Bluetooth, specifically Bluetooth Low Energy (BLE) and classic Bluetooth. BLE has been designed as a low-power solution for control and monitoring applications, in contrast with classic Bluetooth. Finally, the third technology evaluated has been the 802.11p radio, which is the standard to add wireless access in vehicular environments (WAVE). It must be pointed out that for this technology we have performed new simulations in the 5.9 GHz frequency. Three different modulations with three different data rates have been analyzed for the 802.11p radio standard. The transmitted power and sensitivity for the different communication systems considered is shown in Table 4.  Figure 15 shows the linear distributions of received power along the X-axis for different values of Y for different heights in the case of ZigBee and 802.11p radio and for 1.4 m height in the case of Bluetooth. In order to clarify which distribution lines are those represented in Figure 15 in the real scenario, the aerial view of these distribution lines are shown in Figure 16. It is observed that the distribution lines of received power, which corresponds with the transmitter antenna placed at the traffic light #1 are represented. These simulations results are shown in a bi-dimensional plane in Figure 5a.   Figure 15 shows the linear distributions of received power along the X-axis for different values of Y for different heights in the case of ZigBee and 802.11p radio and for 1.4 m height in the case of Bluetooth. In order to clarify which distribution lines are those represented in Figure 15 in the real scenario, the aerial view of these distribution lines are shown in Figure 16. It is observed that the distribution lines of received power, which corresponds with the transmitter antenna placed at the traffic light #1 are represented. These simulations results are shown in a bi-dimensional plane in Figure 5a.   Figure 15c shows the comparison between BLE with the maximum and minimum transmitted power and the higher and lower receiver sensitivity depending of which type of BLE has being used. Figure 15d shows the comparison between Class 1, Class 2, and Class 3 of classic Bluetooth. Figure 15e,f shows the comparison between different devices with different modulations and data rates of the 802.11p standard, for different values of Y. It can be seen that there are some points where the signal goes down below the sensitivity level. The comparison of the linear distribution of received power along the X-axis and Y-axis for different lines in the considered scenario with the different sensitivity levels, leads to radio planning engineers to know the radio coverage of the different systems. These results, along with the different planes of received power achieved for the whole scenario with the 3D Ray Launching Code, gives precise results but no homogeneous. As stated before and it is represented in Figure 7, the received   Figure 15c shows the comparison between BLE with the maximum and minimum transmitted power and the higher and lower receiver sensitivity depending of which type of BLE has being used. Figure 15d shows the comparison between Class 1, Class 2, and Class 3 of classic Bluetooth. Figure 15e,f shows the comparison between different devices with different modulations and data rates of the 802.11p standard, for different values of Y. It can be seen that there are some points where the signal goes down below the sensitivity level.   Figure 15c shows the comparison between BLE with the maximum and minimum transmitted power and the higher and lower receiver sensitivity depending of which type of BLE has being used. Figure 15d shows the comparison between Class 1, Class 2, and Class 3 of classic Bluetooth. Figure 15e,f shows the comparison between different devices with different modulations and data rates of the 802.11p standard, for different values of Y. It can be seen that there are some points where the signal goes down below the sensitivity level. The comparison of the linear distribution of received power along the X-axis and Y-axis for different lines in the considered scenario with the different sensitivity levels, leads to radio planning engineers to know the radio coverage of the different systems. These results, along with the different planes of received power achieved for the whole scenario with the 3D Ray Launching Code, gives precise results but no homogeneous. As stated before and it is represented in Figure 7, the received  Figure 15, along the X-axis and Y-axis.
The comparison of the linear distribution of received power along the X-axis and Y-axis for different lines in the considered scenario with the different sensitivity levels, leads to radio planning engineers to know the radio coverage of the different systems. These results, along with the different planes of received power achieved for the whole scenario with the 3D Ray Launching Code, gives precise results but no homogeneous. As stated before and it is represented in Figure 7, the received power can be approximated for regular coverage zones. Depending of these results, and considering also the capacity of the system, an optimal design of the wireless sensor network coverage can be obtained. In order to represent a practical example, a schematic view of a final implementation with omnidirectional antennas for the different systems analyzed in the considered scenario is shown in Figure 17. It must be pointed out that the different radio coverage of the different circles has been obtained according to simulation results represented in Figure 15, and also taking into account the complete planes of received power, which have been obtained for the whole scenario and several examples have been represented in Figure 5 (XY planes), Figure 6 (YZ planes) and Figure 7 (XY planes divided by zones). power can be approximated for regular coverage zones. Depending of these results, and considering also the capacity of the system, an optimal design of the wireless sensor network coverage can be obtained. In order to represent a practical example, a schematic view of a final implementation with omnidirectional antennas for the different systems analyzed in the considered scenario is shown in Figure 17. It must be pointed out that the different radio coverage of the different circles has been obtained according to simulation results represented in Figure 15, and also taking into account the complete planes of received power, which have been obtained for the whole scenario and several examples have been represented in Figure 5 (XY planes), Figure 6 (YZ planes) and Figure 7 (XY planes divided by zones).

BLE with Maximum Transmitted Power
Bluetooth Class 1 Another important parameter which is important when designing a WSN is capacity. Once coverage levels are satisfied, we must have adequate capacity in order to have a good quality of service in the communication link. The channel capacity depends on the number of users who are connected to the same communication link at the same time, the data rate of the transceivers and the number of gateways in which the information is gathered. Figures 18 and 19 show the channel capacity vs. the number of users for different data rates considered and for different number of Another important parameter which is important when designing a WSN is capacity. Once coverage levels are satisfied, we must have adequate capacity in order to have a good quality of service in the communication link. The channel capacity depends on the number of users who are connected to the same communication link at the same time, the data rate of the transceivers and the number of gateways in which the information is gathered. Figures 18 and 19 show the channel capacity vs. the number of users for different data rates considered and for different number of gateways. It is observed as expected that the channel capacity increases for the same number of users, for every data rate considered, if the information is gathered with a higher number of gateways. gateways. It is observed as expected that the channel capacity increases for the same number of users, for every data rate considered, if the information is gathered with a higher number of gateways.  The great difference in the channel capacity depending on the data rate considered is also very significant. This leads us to conclude that it is highly important to adequately fix the number of gateways in the design phase, depending of the data rate of the transceivers and the expected number of users. gateways. It is observed as expected that the channel capacity increases for the same number of users, for every data rate considered, if the information is gathered with a higher number of gateways.  The great difference in the channel capacity depending on the data rate considered is also very significant. This leads us to conclude that it is highly important to adequately fix the number of gateways in the design phase, depending of the data rate of the transceivers and the expected number of users. The great difference in the channel capacity depending on the data rate considered is also very significant. This leads us to conclude that it is highly important to adequately fix the number of gateways in the design phase, depending of the data rate of the transceivers and the expected number of users.

Application
Once we have modeled the scenario in terms of potential restrictions given by wireless channel operation and the required radio-planning phase, the main goal of the system is to produce a visualization platform presenting a city's traffic status on real time. The main stakeholder is the municipality. It is worth noting that the results previously obtained in Section 3 provide insight into the characteristics of wireless propagation for the scenario under analysis. In this way, estimation of path loss is derived, a required information in order to perform the radio-planning analysis which is described in Section 5. These results provide the density of nodes as a function of the employed transceiver. This information in turn is used in order to adequately dimension the applications which can be executed, as well as the location of the nodes within the context of the specific items of the application, such as pedestrian monitoring. The proposed application that has been implemented, its characteristics and the functional/system architecture will be described in this section.
The traffic in a city, as many other parameters, is very time/day dependent. Mobility patterns highly depend on the day and time of the day. Therefore, traffic light policy is a very complex task and difficult to execute in real time. Additionally, each decision at one place in the city impacts directly or indirectly on other areas of the city. It is a very complex task to maintain a fluid traffic flow in the whole city. In fact, it could be impossible in some situations like blocked streets due to meteorological issues.
In order to analyze the traffic in a city, usually the city is logically divided in areas. Each area defines a conflict region where cars, pedestrian and bicycles are inputs, which should be processed to get out of the conflict region (some of them could remain inside). The conflict area is modeled using queuing theory. Each cross is a join of streets at a common point. Then, each street is modeled as a queue and server, with a processing policy. Mainly, this policy is first-in-first-out (FIFO). Traffic lights are the components of the system that regulate the order to process inputs on queues (server order). In this model, the arrival process of the inputs is a key factor to analyze each conflict area. Such a model is applied mainly for worst and best case modeling. This approach models different conflict areas independently, and not the whole city.
In order to make decisions about traffic light policies, it is important to monitor the traffic, and to store historic data. Such data will be used to adjust the model: determine the arrival process, and the server processing time. However, such simulations are particular cases and while useful to analyze worst case, and average case, it is difficult to optimize general parameters like average time to cross the city.
In our case, we propose another option to analyze a city's traffic by gathering traffic information on-line and presenting it on a visual interface. In this case, traffic managers could analyze in real time the situation of the traffic in the city. The collected information could be used to adjust the previous model (queue network model to analyze the traffic by conflict areas).
In our case, we gather such information based on a WSN. Each sensor is located at a traffic light. Each one collects information about the number of cars and the time interval between each pair of cars. Periodically, we integrate this information in a single packet, and transfer it to a central control, where data is stored and processed. In the initial test, we have deployed 20 nodes within the scenario under analysis, although there was no limitation to increase the network size, scalable as a function of the system requirements. In relation with the sensors deployed, we have focused mainly on weather, pollution (CO 2 and NO) and presence. The employed nodes allow monitoring other items, such as multiple polluting substances, ambient conditions, liquid flow or existence of electric currents, to name a few.
The graphical interface has been implemented in order to illustrate the city's traffic state with the aid of lines and different colors. We have parameterized our sensors to evaluate levels of congestion. We use two parameters to establish line's color, and type (points or partial lines): arrival rate (AR), to count the number of cars arriving at the traffic light each time unit; and processing rate (PR), to measure the time the car requires to pass next to the traffic light.
Both parameters are measured at the traffic light, and integrated periodically to evaluate the average value of the above parameters. Based on them, we evaluate the Congestion Threshold (CTh, as the ratio between AR and PR), which represents the level of congestion of a traffic lane. With blue lines we identify fluid traffic (CTh = 1, orange is used to denote partial congestion (1 < CTh < 5), and red means dense traffic (5 < CTh). The style of the line (points or small lines) denotes the arrival rate, using points to identify the case of low arrival rate. Figure 20 illustrates the traffic congestion at two different locations at the city of Pamplona.
Sensors 2016, 16, 1140 19 of 25 Both parameters are measured at the traffic light, and integrated periodically to evaluate the average value of the above parameters. Based on them, we evaluate the Congestion Threshold (CTh, as the ratio between AR and PR), which represents the level of congestion of a traffic lane. With blue lines we identify fluid traffic (CTh = 1, orange is used to denote partial congestion (1 < CTh < 5), and red means dense traffic (5 < CTh). The style of the line (points or small lines) denotes the arrival rate, using points to identify the case of low arrival rate. Figure 20 illustrates the traffic congestion at two different locations at the city of Pamplona.

Application Architecture
In our scenario we have deployed numerous sensors. Each one has common software to execute. Each node communicates its information to a common server. The server collects the information, and present the information to the users. As we have similar tasks at different locations, their results are collected in a common point, and it is not required any temporal restriction. We have deployed a message bus software architecture. The main advantage is the scalability of the approach, and minimization of application coupling. Each sensor can work independently. Figure 21 depicts the message bus architecture, which uses a software system that can receive and send messages using one or more communication channels, so that applications can interact without needing to know specific details about each other. Interaction between applications is accomplished by passing messages (asynchronously) over a common bus. We implement the message bus architecture using a messaging router, and it is implemented using message queuing.
The motes placed around the scenario collect traffic information from sensors (microwave radar, image processing, infrared and ultrasonic in our case) and introduce it into the system by means of the message bus. The information sensed is persistently stored into a database. Data collected is analyzed in order to detect traffic congestion levels. Such congestion levels are presented to experts (traffic managers) in order to allow them to take decisions within the existing traffic control systems. The congestion level service may be used to provide real-time information to experts or either it could be embedded as part of a traffic control system. Figure 20 shows the traffic monitoring tool that allows experts to visualize the traffic incidences observed by the sensor network. Making use of the publish/subscribe paradigm, the congestion level service pushes the traffic flows to those monitoring tools that have subscribed the service. In the same way, the traffic control system may also subscribe this service, so that the traffic control system has a new source of information to be considered in decision algorithms and internal processes.

Application Architecture
In our scenario we have deployed numerous sensors. Each one has common software to execute. Each node communicates its information to a common server. The server collects the information, and present the information to the users. As we have similar tasks at different locations, their results are collected in a common point, and it is not required any temporal restriction. We have deployed a message bus software architecture. The main advantage is the scalability of the approach, and minimization of application coupling. Each sensor can work independently. Figure 21 depicts the message bus architecture, which uses a software system that can receive and send messages using one or more communication channels, so that applications can interact without needing to know specific details about each other. Interaction between applications is accomplished by passing messages (asynchronously) over a common bus. We implement the message bus architecture using a messaging router, and it is implemented using message queuing.
The motes placed around the scenario collect traffic information from sensors (microwave radar, image processing, infrared and ultrasonic in our case) and introduce it into the system by means of the message bus. The information sensed is persistently stored into a database. Data collected is analyzed in order to detect traffic congestion levels. Such congestion levels are presented to experts (traffic managers) in order to allow them to take decisions within the existing traffic control systems. The congestion level service may be used to provide real-time information to experts or either it could be embedded as part of a traffic control system. Figure 20 shows the traffic monitoring tool that allows experts to visualize the traffic incidences observed by the sensor network. Making use of the publish/subscribe paradigm, the congestion level service pushes the traffic flows to those monitoring tools that have subscribed the service. In the same way, the traffic control system may also subscribe this service, so that the traffic control system has a new source of information to be considered in decision algorithms and internal processes. Data integration is performed at motes. Each mote collects the information from sensors and aggregates this information in order to minimize the number and length of traffic messages. Each mote communicates the information, using a ZigBee interface, to a gateway. The gateway inserts the information on the message bus where a set of parallel process manage the existing queues and store the information into a database. The message bus includes a set of message queues geographically distributed in order to collect information according to the geographical location of the sensors. The queue hierarchy allows the arranged gathering of the traffic information, avoiding duplicated or incomplete information. A pool of database connections allows maximizing the data insertion rate into the database. Figure 22 depicts the gathering process. The message bus uses services for communication between the bus and components (sensors) attached to the bus. In addition, the bus provides services that transform message from one format to another, allowing sensors using incompatible message formats to communicate with each other. In Figure 19 we can see the capacity of the system in terms of number of users. This means applications can be hosted in the cloud instead of on an enterprise network. The cloud improves the scalability of Data integration is performed at motes. Each mote collects the information from sensors and aggregates this information in order to minimize the number and length of traffic messages. Each mote communicates the information, using a ZigBee interface, to a gateway. The gateway inserts the information on the message bus where a set of parallel process manage the existing queues and store the information into a database. The message bus includes a set of message queues geographically distributed in order to collect information according to the geographical location of the sensors. The queue hierarchy allows the arranged gathering of the traffic information, avoiding duplicated or incomplete information. A pool of database connections allows maximizing the data insertion rate into the database. Figure 22 depicts the gathering process. Data integration is performed at motes. Each mote collects the information from sensors and aggregates this information in order to minimize the number and length of traffic messages. Each mote communicates the information, using a ZigBee interface, to a gateway. The gateway inserts the information on the message bus where a set of parallel process manage the existing queues and store the information into a database. The message bus includes a set of message queues geographically distributed in order to collect information according to the geographical location of the sensors. The queue hierarchy allows the arranged gathering of the traffic information, avoiding duplicated or incomplete information. A pool of database connections allows maximizing the data insertion rate into the database. Figure 22 depicts the gathering process. The message bus uses services for communication between the bus and components (sensors) attached to the bus. In addition, the bus provides services that transform message from one format to another, allowing sensors using incompatible message formats to communicate with each other. In Figure 19 we can see the capacity of the system in terms of number of users. This means applications The message bus uses services for communication between the bus and components (sensors) attached to the bus. In addition, the bus provides services that transform message from one format to another, allowing sensors using incompatible message formats to communicate with each other. In Figure 19 we can see the capacity of the system in terms of number of users. This means applications can be hosted in the cloud instead of on an enterprise network. The cloud improves the scalability of the system, because multiple instances of motes can be attached to the bus in order to handle multiple requests at the same time.

Congestion Level Service
The analysis of the information stored at the database is performed using a fuzzy algorithm developed with the aid of Matlab. We have trained an Adaptive Neuro-fuzzy Inference System (ANFIS) according to data collected from sensors and stored into the database.
The fuzzy system has two input and one output variables. The input variables are the traffic density observed during a given period (traversingVehiclesRate, number of vehicles observed from a traffic light during the green period); and the number of vehicles awaiting the next traffic light cycle (#waitingGreen). The congestion level is the output variable. The output values vary on the range [1,10] as described above. The membership degree of each variable is given by three fuzzy sets (low, medium and high) provided by the ANFIS. Table 5 briefly summarizes the decision rules. In this paper, in order to illustrate congestion level service operation we have deployed a simplified interface, where data collected is used to present the current status of the traffic in the city. However, such information is processed to augment the information at each point taking into account previous analysis.
In our case, a first approach has been the evaluation of the status of each conflict area considering the congestion level at different days. We integrate on a graphical view the level of congestion during different periods of time. We group the week on labor days (LD) and holidays (HD). We average the congestion level for both groups, and provide it to the user (traffic manager). In Figure 20 it corresponds to the pie chart next to some traffic lights. The colors of the pie chart correspond to the congestion level at each location. Figure 23 shows some of the parameters and indicators considered (temperature, wind orientation and intensity, pollution degree . . . ). The average, maximum and minimum temperatures observed during the three first weeks of June 2016 are depicted in Figure 23b, while Figure 23a shows the variation of the temperature each hour during the same period. The graphs on the bottom show the pollution degree (CO, SO 2 and NO 2 ) and the wind direction (bottom right). Traffic managers can use such information in order to temporarily modify traffic light policies, or to send human resources (police) to regulate the traffic. Currently we are evaluating a few parameters, but the sensors deployed at traffic light have an enormous potential to provide better services for citizens. Nowadays, we obtain situational awareness from temperature and pollution sensors and also from passive infrared sensor (PIR), which aid to detect the presence of pedestrians and vehicles, allowing the estimation of the degree of occupancy of both lanes and crosswalks. Figure 24a shows the relative humidity observed, while the graph on the right shows the percentage of operating time during which the infrared sensor detects pedestrians. The system here proposed provides complete and adequate situational awareness to municipal experts in order to allow the modification of the traffic light cycles according to the specific needs of each instant considering multiple environmental factors in addition to the traditional vehicle counting.

Conclusions
In this work, the demands for modeling the radio channel in an urban city area for intelligent traffic light control systems are presented. The obtained results show the complexity of this type of scenario, with traffic lights, vehicles, people, buildings, vegetation and urban environment. It can be concluded that the morphology as well as the topology of the scenario play a key role in the estimation of radio signal propagation, due to the strong impact of multipath components in the overall loss mechanism of the propagating radio waves. The validation of the 3D Ray Launching algorithm has been done with a good match in the comparison with real measurements. The use of simulation tools, such as the 3D Ray Launching code, is a great aid for radioplanning analysis with the aim of deploying WSN within this type of environments. Moreover, the context aware scenario implementation using different wireless technologies can be properly designed with a detailed characterization and analysis of the radio wave propagation for different frequencies. The system

Conclusions
In this work, the demands for modeling the radio channel in an urban city area for intelligent traffic light control systems are presented. The obtained results show the complexity of this type of scenario, with traffic lights, vehicles, people, buildings, vegetation and urban environment. It can be concluded that the morphology as well as the topology of the scenario play a key role in the estimation of radio signal propagation, due to the strong impact of multipath components in the overall loss mechanism of the propagating radio waves. The validation of the 3D Ray Launching algorithm has been done with a good match in the comparison with real measurements. The use of simulation tools, such as the 3D Ray Launching code, is a great aid for radioplanning analysis with the aim of deploying WSN within this type of environments. Moreover, the context aware scenario implementation using different wireless technologies can be properly designed with a detailed characterization and analysis of the radio wave propagation for different frequencies. The system

Conclusions
In this work, the demands for modeling the radio channel in an urban city area for intelligent traffic light control systems are presented. The obtained results show the complexity of this type of scenario, with traffic lights, vehicles, people, buildings, vegetation and urban environment. It can be concluded that the morphology as well as the topology of the scenario play a key role in the estimation of radio signal propagation, due to the strong impact of multipath components in the overall loss mechanism of the propagating radio waves. The validation of the 3D Ray Launching algorithm has been done with a good match in the comparison with real measurements. The use of simulation tools, such as the 3D Ray Launching code, is a great aid for radioplanning analysis with the aim of deploying WSN within this type of environments. Moreover, the context aware scenario implementation using different wireless technologies can be properly designed with a detailed characterization and analysis of the radio wave propagation for different frequencies. The system has been tested by means of an ad hoc application and its corresponding platform, providing real time traffic monitoring within the scenario under analysis. Results show that the use of deterministic tools in WSN deployment can aid in providing optimal layouts in terms of coverage, capacity and energy efficiency of the network.

Acknowledgments:
The authors wish to acknowledge the financial support of project TEC2013-45585-C2-1-R, funded by the Spanish Ministry of Economy and Competiveness.
Author Contributions: Leyre Azpilicueta, Peio López-Iturri, Erik Aguirre, Carlos Martínez and Francisco Falcone conducted the characterization of wireless propagation mechanisms in the outdoor scenario, as well as the WSN campaign of measurements. José Javier Astrain and Jesús Villadangos participated in the development of the WSN and the application for such scenario.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: