Diffraction Loss Prediction of Multiple Edges Using Bullington Method with Neural Network in Mountainous Areas

This paper proposes a neural network approach to improve the Bullington method by using parameters obtained from ignored obstacles in mountainous areas. Measurements were performed in mountainous areas to compare the prediction accuracy of propagation loss. And the measured data were used for neural network training. A detailed description of the input parameters of the proposed neural network is presented. The prediction performances were improved by up to 3.20 dB in the average error and 2.11 dB in the standard deviation of errors by the proposed method when compared to traditional diffraction methods.


Introduction
Propagation loss prediction is one of the main problems in planning of radio communication links.And it is also very important to design mobile radio systems [1].In Korea, the mountainous area is about 70% of the whole land area and the hilly terrain can cause multiple obstacles to the signal path.In hilly terrains, line-of-sight propagation is not possible at all the places and the diffraction becomes dominant on propagation [2].However, accurate prediction of the diffraction losses is still challenging for realistic propagation environments such as mountainous regions [3].Practically simplified models are used such as Bullington, Deygout, Causebrook, and Giovanelli for computational efficiency [3,7].Those models use knife edges to replace the mountain peaks and ridges for computing diffraction losses [2,3].In contrast, more sophisticated techniques such as UTD (uniform theory of diffraction) can be used to improve accuracy but need much more computation time and detailed information of terrains and obstacles such as conductivity of the materials [4].
Here, we propose a neural network with Bullington model in the presence of three or more edges.To improve the prediction accuracy of diffraction loss, the proposed method uses ignored edges in the Bullington method to calculate additive diffraction losses.Parameters obtained from terrain elevation information were used as inputs to the neural network.The details are described in Section 3.
Measurements were conducted in mountainous areas.The measured data were used to evaluate the diffraction loss prediction and train the neural network using training algorithm.
The structure of this paper is as follows.Section 2 introduces the traditional diffraction models.And the proposed method and training are described in Sections 3 and 4, respectively.Section 5 illustrates the measurement campaign.Section 6 gives the evaluation of the prediction performance.Conclusions are drawn in Section 7.

Diffraction Methods
2.1.Bullington Method.In the Bullington method, the real terrain is reduced to a single equivalent knife edge.The location of the equivalent knife edge is the point at which the extended lines joining the transmitter and receiver to their respective dominant (the greatest angle of elevation as viewed from transmitter or receiver) obstacles meet as shown in Figure 1 [5].Then the diffraction loss is computed using (1) and ( 2) [4,6].This method has the advantage of simplicity but obstacles below the paths of the horizon rays can be ignored.And this may cause large prediction errors.In general, it underestimates path losses [4,5].
The geometrical parameters are combined in a dimensionless parameter denoted by ν as shown in [6]: where h is the height of the top of the obstacle above (or below) the straight path line and d 1 and d 2 are the distances of the path ends from the top of the obstacle.If the obstacle is below the straight path line, then h is negative.The diffraction loss as a function of ν is presented in (2) for ν greater than −0.78 [6].It should be noted that the diffraction loss can be avoided for ν ≤ −0.78 [6].
The approximated diffraction loss, L B , due to the equivalent Bullington edge is equal to (2), and the total path loss using the Bullington method is [4] where L fs is free space loss [16].

Deygout Method.
The Deygout method is drawn in Figure 2 for a path with three obstacles.The first step is to compute ν parameter using (1) for each edge alone, as if all other edges were absent, that is, all ν parameters are calculated for the paths Tx-A-Rx, Tx-B-Rx, and Tx-C-Rx in Figure 2. If edge B is the main edge in Figure 2, then the diffraction losses, which are J(ν D_A ) and J(ν D_C ) , for edge A and edge C are found with respect to a line joining the main edge to the Tx and Rx.And those are added to the main edge loss (J(ν D_B )) to obtain a total approximated diffraction loss (L D ), by Deygout method [4,7].The diffraction loss at each edge is calculated using (2).The total path loss using the Deygout method is presented in (4) [4].This procedure can be repeated until all the edges have been considered for more than three edges [4,8].However, it is common to compute the total loss as the sum of one main edge and two subsidiary main edges on either side in practical applications [4,8].
where L fs is free space loss [16].

Causebrook Method (Correction).
To reduce an overestimation problem of the Deygout method, Causebrook proposed an approximate correction derived from the exact analysis of the two-edge solution [24].The corrected form is given in where L D is the diffraction loss from (4).The correction factors C 1 and C 2 are presented below [24]: where cos  2 International Journal of Antennas and Propagation 2.4.Giovanelli Method.Another development of the Deygout method has been proposed [24,25].The concept for the Giovanelli method is shown in Figure 3.It is assumed that edge A is the main edge in this case.Then a reference field point F ′ is found by projecting AB onto F ′ F ″ and h 1 ″ is defined in [24] h where H = h 2 + md 3 and m = (h 2 − h 1 )/d 2 .The effective height for the secondary edge B is then given by [24] The total diffraction loss is now given in [24] where J(•) is the loss from ( 2).This method is conveniently extended to more than two edges by recursively applying the above procedure [24].

Proposed Method
As mentioned in Section 2, there may be ignored edges in the Bullington method, which generally cause underestimation.Therefore, for multiple edges on the path, we propose a method to obtain the total diffraction loss by adding the loss from the Bullington method and the loss from the ignored edges.Where the diffraction loss from the ignored edges is calculated using a neural network as shown in Figure 4.There are 15 inputs of the proposed neural network, all of which can be obtained from DTM (digital terrain model).Inputs 1 through 9 are consisted of the heights of Tx and Rx, the three largest ν of ignored obstacles, the distances between Tx and three ignored obstacles, and the distance of the total path.Inputs 10 through 15 are approximate front and rear slopes of ignored edges.An approximate slope is calculated as the height difference between the edge and the point spaced by 30 m from the edge as shown in Figure 5.The use of slopes is intended to improve the prediction performance by considering the obstacle radius of curvature [9,10].This is an approximate application for the obstacle radius of curvature in the diffraction loss calculation.The proposed method is applied when there are at least three edges to include an ignored edge on the paths.For paths with three or more edges, the prediction of the proposed method is compared with the two traditional diffraction methods in Section 2, and the results are presented in Section 6.

Training
Figure 6 shows a simple neuron.The neuron is presented with inputs as shown in and an output value of a neuron is produced as shown in [11] net = w T x, 11 where (•) T denotes the transpose, and the weights of a neuron w are presented in The activation function has been selected to be the commonly used hyperbolic tangent sigmoid transfer function [11].Figure 7 shows the architecture of the neural network with inputs of the proposed method.The inputs in Figure 7 correspond to those shown in Figure 5. Prior to training, data with three or more obstacles should be extracted from all of the measurement data.And it is necessary to divide all the measurement data with three or more edges into each of the two disjoint sets which are training and testing sets.The training sets are only used in the training procedure, and the testing sets are used for the evaluation of prediction performances.The training data are randomly chosen to be about 50% of the measurement data with three or more edges in this paper.The input vectors shown in Figures 5 and 7 can be obtained from the DTM corresponding to the transmitting and the receiving positions as described in Section 3. The receiving power may be expressed in (13) in wireless environments [12,13].
where P t is the transmitting power, G t and G r are the gains of the Tx and Rx antennas.L p is the path loss.Bullington proposed that the diffraction loss has to be added to the free space loss when the path is obstructed [4,14,15].In order to find target values, ( 14) can be obtained from (13) by using the Bullington method as shown in (3).
where L fs is the free space loss and L B is the Bullington diffraction loss.L t is the target value, which is desired loss due to ignored edges for training.The target value can be obtained using (14) where P r is measured power in dBm, G t and G r are 2.15 dBi gain, and P t is 36 dBm as drawn in Figure 8.And the free space loss L fs is computed using [16] L fs = 32 4 + 20log f + 20log d dB , 15 3 International Journal of Antennas and Propagation where f is the frequency (MHz) and d is distance (km).
Here, the error is defined for the feed-forward neural network as shown in where o(j) is the output of the neural network corresponding to an input vector, L t (j) is a target value, and j is the index number of training data.In the backpropagation procedure, the weights of the neural network are adjusted by minimizing the mean square error as follows.
The Levenberg-Marquardt algorithm [17] is used to update the weights of the proposed neural network.The    The transmitting signal is 1399 MHz continuous wave.This frequency was temporarily allowed in UHF band for the measurement campaign.The transmitting antenna is fixed on tripod of 6 m height, and the receiving antenna is placed on the roof of a car at a height of 1.8 meter above ground.The power of the transmitter is 36 dBm with a 2.15 dBi gain antenna and a cable of 2.71 dB loss.The measurement data are obtained from a spectrum analyzer connected to a 2.15 dBi gain antenna and a 40 dB gain low-noise-amplifier.A measurement configuration diagram is shown in Figure 8.The measurement data with three or more obstacles from the real mountain areas are used to produce train sets for the proposed method and later evaluate the performances of Bullington, Deygout, and the proposed method.The numbers of data with three or more obstacles are as follows: 1676, 5863, and 1147 for the west side of Mt.Maebong, Mt.Ongma, and Mt.Hwaak, respectively.And the numbers of training data are as follows: 838 and 2932 for the west side of Mt.Maebong and Mt.Ongma.And all training data consist of 70% for training and 30% for validation.

Comparison of Prediction Performance
The proposed method is applied to real mountainous environments.Two scenarios are discussed for the comparisons of diffraction methods in mountainous areas as shown in Figures 12-14.The details of the scenarios are drawn in Sections 6.1 and 6.2.In Section 6.1, the area-specific performance of the proposed method is investigated as shown in Figures 12 and  14.In Section 6.2, we try to find out the performance of the proposed method in another area as shown in Figure 14.7 International Journal of Antennas and Propagation 6.1.Prediction Results of the Proposed Method with Traditional Methods.The proposed method is evaluated by comparing with the four diffraction methods described in Section 2. As described in Section 4, the testing sets are applied to the proposed method, Bullington method, Deygout method, Causebrook method, and Giovanelli method for the evaluation of prediction performances.The error is defined as shown in [2,20,21] e i = p i − m i dB , 18 where p is the prediction of received power, m is the measured value, and i is the point number of measurement data.
The prediction values are obtained by using (14) with outputs of the neural network.The prediction values are calculated for each method as follows: (1) Prediction values for the proposed method using ( 14) (2) Prediction values for the Bullington method using (3) and ( 13) (3) Prediction values for the Deygout method using ( 4) and ( 13) (4) Prediction values for the Causebrook method using (5) and ( 13) (5) Prediction values for the Giovanelli method using ( 9) and ( 13) The average error and the standard deviation of errors from (18) are considered to represent the prediction performance [22,23].
As described in Table 1 and Figures 15 and 16, the results show that prediction performances are improved when using the proposed method in mountain areas with three or more obstacles on the paths.For the paths with three or more obstacles, neural network has been trained with all of the training data from the two areas in the previous subsection.The calculation for the prediction of each method uses (19), ( 20), ( 21), (22), and ( 23) as in Section 6.1.As shown in Table 2 and Figure 17, the comparison results show that prediction performance is improved by using the proposed method (pretrained neural network) in the area where data are not used in training.

Conclusion
A neural network-based diffraction method is proposed, and the method can be applied with the DTM in mountainous areas which have three or more obstacles on the paths.With a set of measured data of real mountain areas, the proposed method was evaluated with Bullington, Deygout, Causebrook, and Giovanelli methods and the results showed that the prediction performances were improved by up to 3.20 dB in the average error and 2.11 dB in the standard deviation of errors by the proposed method.The remaining results are shown in Table 1.
In addition, the prediction performance of the proposed method is validated by applying the proposed method which was trained with data from two areas to the area where data are not used in training.The results show that prediction performance is improved as described in Table 2.We used 838 data as minimum train set within about 10 km between Tx and Rx.To use the proposed method at longer distances, it is recommended to train the neural network with data from corresponding distances.
It is expected that the proposed method based on neural network would show better prediction results with more measurement data in hilly mountainous areas with multiple obstacles.
and L 1 and L 3 are the losses due to edges A and C, as if they existed on their own between the Tx and Rx.And d 1 , d 2 , d 3 , and d 4 are shown in Figure 2.

Figure 5 :
Figure 5: Slopes of an ignored edge.
Measurements have been performed in the three mountain areas as shown in Figures 9-11.The heights of the selected mountains for measurement purposes are as follows: 746.7 m, 533.1 m, and 1420.2 m for the west side of Mt.Maebong in Inje, Mt.Ongma in Boryeong, and Mt.Hwaak in Gapyeong, respectively, in Figures 9-11.The selected mountains are typical hilly mountains in the Republic of Korea.In these mountains, it was convenient and safe to conduct measurements using a vehicle.The red triangle marks in Figures 9-11 are the transmitting positions, and the circles are measuring points.The total number of measurement data in Figures 9-11 is as follows: 80,163 for the west side of Mt.Maebong, 77,215 for Mt.Ongma, and 35,057 for Mt.Hwaak.

Figure 9 :
Figure 9: Measurements in the west side of Mt.Maebong.

6. 2 .
Prediction Results of the Trained Neural Network in the Area Where Data Are Not Used in Training.The proposed

Figure 15 :
Figure 15: Comparison of the measurements and predictions in the west side of Mt.Maebong.

Figure 16 :
Figure 16: Comparison of the measurements and predictions in Mt.Ongma.

Figure 17 :
Figure 17: Comparison in the area where data are not used in training (Mt.Hwaak).

Table 1 :
Comparison of prediction performance.

Table 2 :
Comparison in the area where data are not used in training (Mt.Hwaak).