Narrow-Band Radio Propagation Prediction Based on a Highly Accurate Three-Dimensional Railway Environment Model

In this paper, radio propagation predictions were performed by a self-developed ray tracing simulator based on highly accurate three-dimensional (3D) environment models obtained by oblique aerial photography (OAP). In order to support ray tracing computations, the highly complex environment models were managed by three steps, including model separation, electromagnetic parameter assignment, and model lightweighting. Simulations and real-world experiments were conducted to test the proposed method. The results show that the root mean square errors (RMSE) in the predicted received signal strength curves were achieved 4 dB with comparably low computational complexity in the selected railway scenarios.


Introduction
Nowadays, the Global System for Mobile Communications-Railway (GSM-R) networks in China mainly transmit the train-to-trackside safety-related information for China Train Control System Level 3 (CTCS-3), which makes the wireless communication networks vital for high-speed railway controlling [1]. During the railway wireless network design, it is a basic step to locate the communication base stations (BS). If the network planning is conducted based on inaccurate radio propagation prediction results, it is highly potential that the deployment of BSs is not optimized and the telecom infrastructure needs to modify during constructions, which will surely add to the cost. Since the costs of the facilities supporting the BSs, including the antenna towers, backhaul equipment, communication cables, electric power supply, and housing, are much higher than the expense of the wireless system devices, the BS site locations affect not only the quality of the railway wireless communication networks but also the overall investment of the railway constructions. As a result, accurate radio propagation predictions are significant in railway wireless system design [2,3]. The GSM-R downlink frequency band is 930-934 MHz in China. The existing channel models for radio propagation prediction in this specific band are mainly classified into two categories: empirical channel models and deterministic channel models. Empirical channel models, such as the Hata-Okumura model [4,5], the Stanford University Interim (SUI) model [6], and the Cost 231 model [7] are widely used owing to their high execution speed. However, the empirical channel models are usually inaccurate before measurement-based model corrections [8,9]. Moreover, even if the model corrections are performed, the reliance on the measurements in the selected environment makes the empirical channel models less reliable when the environment changes. In order to increase the robustness and accuracy of the radio propagation predictions, deterministic channel models are developed [10][11][12]. By literally computing different radio propagation factors, the deterministic channel models obtain radio propagation prediction results with higher accuracy [13]. Among the deterministic channel models, ray tracing is a widely accepted tool. Ray tracing is a three-dimensional (3D) ray-optic technique computing the influence of multiple environmental factors that contribute to the signal propagation processes [14]. Direct, transmitted, reflected, and diffracted electromagnetic fields propagating in the environment are coherently summed to simulate the realworld radio propagation. With the fast development of computer CPUs and 3D reconnaissance maps, ray tracing has already become a vital tool in wireless communication system design [15]. Ray tracing is a powerful tool in accurately describing multipath propagation, especially in complex railway scenarios. For high-speed railways, it is required to reduce the bumpiness and ensure safe operation for high-speed trains. As a result, characteristic railway environments, such as cuttings, tunnels, flyovers, and viaducts, are widely implemented in railway constructions to straighten the track, which put together may occupy over 70% of the whole railway line when the line frequently encounters obstructions or waters. The specificity of railway scenarios adds to the difficulty in network planning. Over the past decades, efforts had been made on accurately characterizing the propagation models in those specific railway environments using ray tracing [16][17][18][19]. Commercial tools like CloudRT [15], Winprop, Ranplan, and EDX were also developed to support railway scenarios. Most existing ray tracing methods calculate the radio propagation based on digital elevation model (DEM) or manual 3D modeling. The DEM describes the natural terrain in high resolution, but it is not able to fully present the railway constructions and the surrounding buildings in the environment. To make reasonable environment models, the constructions should be manually added to the 3D environment model, thus causing uncertainty to the simulations. Though the 3D environment models are imperfect, it should be emphasized that ray tracing still performs much more accurate radio propagation prediction than the empirical channel models. To further improve the performance of ray tracing and achieve highly accurate radio propagation predictions, realistic 3D environment descriptions are needed. Light Detection and Ranging (LiDAR) and Oblique Aerial Photography (OAP) are currently two representative 3D modeling technologies that are used to obtain large-scale open-area 3D environment models in high efficiency. LiDAR obtains high precision point clouds in occlusion areas, which describes the natural terrain efficiently [20]. However, the problem with LiDAR is that it is costly to fully present railway constructions and the texture of the environment models should be obtained by other technology so that the model can be applied to ray tracing. In order to produce environment models which are highly consistent with the realistic railway environments, OAP is used in this paper. OAP is a highly precise realistic 3D model restoring technology, which was widely used for producing 3D building textures over the last decades [21,22]. With the development of computer technology and unmanned aerial vehicles (UAV), OAP attracts massive attention in geological surveys [23]. OAP obtains aerial images from different vertical and oblique angles via oblique photography platforms. Most commercial oblique photography platforms are of the "1 + 4" system structure, where one vertical and four oblique cameras are integrated on a single UAV platform. By implementing the inertial naviga-tion system, the obtained high-resolution aerial images are arranged with their corresponding geographic locations [24,25]. After image matching and 3D modeling, the obtained images are put into a consistent coordinate system to produce highly realistic 3D environment models. OAP obtains 3D environment models with high accuracy, which is promising of supporting highly accurate radio propagation predictions in railway scenarios. However, the challenge is that large-scale OAP environment models produced by high-resolution aerial images contain large amounts of data, which result in excessive computational complexity when performing ray tracing computations [26]. In this paper, we proposed a method to manage the OAP environment models before implementing them in a self-developed ray tracing simulator. We took advantage of OAP to accurately establish the environment models and exported them to light-weighted environment models. The proposed method was performed in two typical railway scenarios and propagation prediction results were obtained. The results show that a 4 dB root mean square error (RMSE) was achieved with comparably low computational complexity.
The rest of this paper is organized as follows. Section 2 presents the simulation methods including the ray tracing simulator and environment model management. Section 3 describes and discusses the comparative experiment results of the ray tracing simulation and the field experiment. Finally, Section 4 makes conclusions on the modeling and simulation processes.

Materials and Methods
2.1. Ray Tracing Simulator. In order to make accurate radio propagation predictions in high-speed railway scenarios, a ray tracing simulator was developed. The existing schemes for ray tracing include the shooting and bouncing rays (SBR) method [27] and the image method (IM) [28]. Since the railway scenarios are usually distributed on a large physical scale and the radio propagation along the railway track is the primary consideration for most railway applications, only the propagation paths received by the selected observation points are necessarily considered. Therefore, we adopted the IM scheme for the ray tracing simulator to link each observation point to their corresponding source signal before calculating the radio propagation concerning the facets and edges of the obstacles in the environment model. The received power for each communication link at the observation point coherently sums up all multipath components and can be expressed in [15] where Gðθ i , φ i Þ is the receiving antenna gain, λ is the radio working wavelength, η 0 is the characteristic impedance of the vacuum, θ i and φ i are the arrival angles of the ray i to the observation point at E i which is the associated receiving The electric field is in the form of the matrix considering different polarization components of the transmitting antenna, where P r is the radiation power of the transmitting antenna, G M is the transmitting antenna gain, r is the radiated distance between the transmitting antenna and the observation point, E θ i ðθ s , φ s Þ and E φ i ðθ s , φ s Þ are the vertical and horizontal components of the normalized antenna patterns of the transmitting antenna, and T θ i and T φ i represent the propagation path loss caused by the multipath components in the two corresponding polarization directions.
The multipath components include the direct, transmitted, reflected, and diffracted field components with respect to the obstacles in the environment model. For the direct, transmitted, and reflected field components, the ray tracing algorithm went through every triangular facet in the 3D environment model to find the propagation paths before calculating the path loss, implementing Fresnel's equation when necessary. Meanwhile, for the diffracted field components, each of the edges in the 3D environment model was tested using Uniform Theory of Diffraction (UTD) [15]. Because the algorithm computational complexity increases exponentially when there are higher-order reflections and diffractions, which seldom occur in wide-open scenarios, we consider only the first and second orders of reflected and diffracted components in the simulator.
Before conducting the simulations, the antenna patterns for both transmitting and receiving antennas were modeled, which are demonstrated in Figure 1. The 3D antenna pattern was obtained by interpolating the 2D patterns. The transmitting antenna was a directional dual-polarization planar antenna (YX0809X6517-DT0, Dongzhou, China) with 17 dBi gain and 65°horizontal and 9°vertical beamwidth. Meanwhile, the receiving antenna was a shark fin antenna (84002818, HUBER+SUHNER, German) located at the train rooftop. The simulation iteratively traces the rays starting from the transmitting antenna to the receiving antenna to obtain the radio propagation simulation results. The antenna gains of the propagation rays at the transmitting or receiving antennas were determined by their corresponding azimuth and elevation angles with respect to the antenna patterns.

3D Environment Model Management.
Ray tracing simulations were performed in realistic railway environments along the Zhangjiajie-Jishou-Huaihua (ZJH) high-speed railway. We built realistic 3D railway environment models using OAP. During the modeling process, the aerial images were obtained by commercial OAP devices (DG3 Pro, Rainpoo, China) and processed to highly accurate 3D railway environment models using the software ContextCapture. Figures 2(a)-2(c) present the environment models of the three selected railway scenarios, respectively. The models for all scenarios were reconstructed by 5 cm resolution aerial images. Figure 2(a) shows the OAP environment model of scenario 1, which is a typical wide-open railway district with an area of 400 × 1200 m 2 . Typical elements for wide-open railway sections, such as viaduct, cuttings, and tunnel entrances, were found in this specified railway area. Figure 2(b) shows the OAP environment model of scenario 2, which is a railway traffic junction area linking to the Huaihua South railway station in Hunan Province in China with an area of 500 × 2600 m 2 . The surrounding environment in scenario 2 is more complex than that in scenario 1 because of the urban structures and the overpassed railway tracks. Figure 2(c) shows the OAP environment model of scenario 3, which is a valley area with an area of 400 × 2600 m 2 . The bending railway track makes the radio propagation mostly obstructed by the surrounding terrain, which makes a good NLOS (non-line-of-sight) situation for railway wireless communications.
After we obtained the original OAP 3D railway environment models, the models were managed to meet the demands of the ray tracing simulator in three steps, which were the model separation, electromagnetic parameters assignment, and model lightweighting.
First of all, we separated different environment model elements into different submodels from the OAP environment model. Specifically, all environment elements were classified into six categories: vegetation, concrete, water, soil, metal, and cutting slopes. Amount the six categories of elements, the vegetation elements contained the vegetated areas including forests, lawns, and farmland. The concrete elements covered most of the artificial constructions, including the viaduct, surrounding buildings, tunnel entrances, and solid grounds. The water elements included all ponds, streams, and rivers. Meanwhile, all bare ground areas were classified into soil elements and all metal-made structures were classified into metal elements. There was a special case for cutting slopes since the cuttings along the railways contributed massively to the ray tracing simulation results and were often combined material covered with bushes. In addition to the submodel separations, we also separated the artificial constructions from the natural terrain area since the artificial constructions along the railway tracks are usually close to the receiving antennas on the operating railway train and have more influence on the radio propagation. The above-mentioned model separations were conducted manually with the OAP 3D modeling software DpModeler (Wuhan Tianjihang Information Technology Co., Ltd., Wuhan, China).
After we finished the environment model separation, the electromagnetic parameters were assigned to the 3D environment model. To differentiate all six categories of elements in the 3D environment model, different colors, such as green, grey, baby blue, brown, red, and navy blue, were implemented to represent the vegetation, concrete, water, soil, metal, and cutting slope elements,    Wireless Communications and Mobile Computing At last, the environment models were light-weighted before being implemented in the ray tracing simulator. The model lightweighting aimed to reduce the data volume of the environment models and decrease the computational complexity while maintaining the high prediction accuracy of the simulations. The model lightweighting was performed using the ProOptimizer tool in the software 3dsmax (Autodesk, California, USA). The amounts of 3D triangular facets in the original environment models were over 11 million for scenario 1 (Figure 2(a)), 74 million for scenario 2 (Figure 2(b)), and 45 million for scenario 3 (Figure 2(c)), respectively. In this paper, we used different strategies for different submodels. Specifically, artificial constructions, such as cutting slopes, viaducts, tunnel entrances, and   (Figures 2(d)-2(f)), while the rest of the submodel in the environment model, including all surrounding vegetation, soil, and water elements as well as the remaining concrete and metal elements, were considered contributing lighter to the simulation results and were classified as the class B elements (Figures 2(g)-2(i)). During the model lightweighting, the class B elements remained only 0.1% of the model vertices, while the class A elements experienced a relatively lower lightweighting level and 1% of the model vertices remained. Figures 2(j)-2(l) demonstrated the simplified environment models. The overall numbers of triangular facets in the three environment models were reduced to 20k, 84k, and 53k, respectively.
Apart from the 3D environment models for simulation, accurately locating the transmitting antennas is also an advantage of implementing highly accurate 3D railway environment models. According to the GSM-R system design documents along the ZJH high-speed railway, the transmitting antennas in all selected scenarios are shown in Figure 3. Figure 3(a) depicts the location of the transmitting antennas in scenario 1. Antenna A and antenna B were located at the two opposite tunnel entrances. Meanwhile, the electromagnetic environment in scenario 2 is much more complex compared with that in scenario 1; hence, we discussed only the railway area before entering the railway station to simplify the simulation processes and only antenna C was concerned in this area (Figure 3(b)). Figure 3(c) shows the location of the transmitting antennas in scenario 3. Antenna D was located 74 m away from the railway track at the foot of the earth mound beside the tunnel entrance. Antenna E was located close to the railway track at the opposite end of the district. All the antenna parameters are listed in Table 2.

Results and Discussion
3.1. Field Experiment. In order to obtain the real-world radio propagation in the two selected railway scenarios, field experiments were conducted. The received signal strength indicator (RSSI) curves along the railway track were obtained by conducting measurements on operating railway trains. Figure 4 demonstrates the field test environment. The transmitting antenna was installed on the installation platform on the BS tower beside the railway track (Figure 4(a)), and the receiving antenna was the shark fin antenna located at the train rooftop (Figure 4(b)).
The RSSI curves along the railway tracks were obtained by the downlink received signal levels (RXLev) exported from the Abis interface of the GSM-R communication   In compliance with the CTCS-3 standards for the ZJH high-speed railway, coverage of interleaved stations is demanded in the GSM-R wireless system design. As a result, all scenarios are covered with two sets of transmitting systems. Scenario 1 is an intercell handover area with the two remote unit (RU) devices belonging to different cells; hence, we resolved the RSSI curves for each of the transmitting antennas, while scenario 2 performs single-sided wireless coverage with the RU device acting in an active-standby scheme using a single transmitting antenna (antenna C). Scenario 3 is also covered with two sets of transmitting systems. However, the RU devices supporting the two antennas share the same BTS and the transmitted frequency from both antennas are the same. As a result, the receiving antenna is not able to distinguish the downlink signal from each of the two transmitting antennas and the received signal was a compound signal from both two signal sources. All field experiment results shown in Figure 5 were used for verifying the ray tracing simulation results.

Comparative
Results. Radio propagation predictions implementing ray tracing were performed in all selected scenarios. Figure 6 demonstrates the comparisons between the simulation results and the experimental results. The observation points in the simulation were all chosen as 1 m evenly spaced points which were 4.5 m above the railway tracks. The height of the observation point was chosen for the reason that the shark fin antennas on the train rooftop are usually considered about 4.5 m above the railway track. Since the altitude of the railway track along a single railway line is not completely flat, the observation points were set strictly matching the railway track altitudes but not on the same 3D plane. The field test measurement results in all scenarios were obtained with the train operating at a speed about 120 km/h. So in a SACCH multiframe period (480 ms), the train rooftop receiving antenna is about to move 16 m distance. To physically match the simulation results with the obtained measurement results exported from the GSM-R Abis interface, we averaged every 16 adjacent observation points in the simulation.
In scenario 1 and scenario 2, the three antennas in the two scenarios were simulated separately. In Figure 6(a), the received radio signal for antenna 1 was dominated by the direct and first-order reflected field components because of the comparably simple railway scenario with no obstructions. Specifically, the reflections were mainly caused by the cutting slopes beside the railway track. Antenna 2 was also in scenario 1, which was located at the opposite side of the area against to antenna 1. In this situation, there was a trackside building in front of the BS tower. Although the building was not tall enough to obstruct the direct rays, it was located within the Fresnel zone of the transmitting antenna, hence causing some diffraction. As a result, the RSSI curve in Figure 6(b) shows some fluctuations compared with the antenna 1 results. Figure 6(c) demonstrates the simulated results for antenna C in scenario 2 with its corresponding measurement results. Scenario 2 was much more complex compared to scenario 1, which caused more violent signal fluctuations.
In scenario 3, the two antennas transmit radio signals from the RUs sharing the same BTS. As a result, the measured data cannot distinguish the received signal with the transmitting antenna and a compound signal was obtained. In order to make a fair comparison, we also compounded the received signal from both antennas in the simulation. There is an NLOS area for both two transmitting antennas in this scenario, so an obvious signal loss can be found at the distance of 1400 m to 1800 m in Figure 6(d).
To statistically analyze all three simulation results with their corresponding experimental results, the RMSE and the mean absolute error (MAE) of the differences between the measurement results and prediction results are listed in Table 3. All analyses considered only the far-field results regardless of those near field values topped at -47 dBm. The simulation results obtained by the proposed method are consistent with the measurement data, with the RMSE around 4 dB and the MAE around 3 dB. Comparing the results with the existing published works for high-speed railway ray tracing applications based on DEM maps or manual 3D models, which obtained the simulation with RMSE at  [2] and MAE above 5 dB in railway tunnels [19], our proposed method implementing highly accurate 3D railway environment model promotes the accuracy of radio propagation prediction.
By analyzing the result, two main aspects may give rise to RMSE and MAE. The first aspect lies in the ray tracing method. Since ray tracing predicts the received power under ideal conditions, regardless of other external factors, such as    Wireless Communications and Mobile Computing the influence of humidity on vegetation and soil. In our proposed method, we regard all soil material and vegetation material as uniformed materials. However, according to ref. [31], the humidity is highly relevant to the relative permittivity of the material, especially for soil. So the inaccuracy in material parameters causes inaccurate predictions. On the other hand, the second aspect affecting the predicted RMSE lies in the environment model. The real-world open-area environment is not static. Since the model acquirement and the experiment was not simultaneously conducted, the environmental changes may make the 3D environment model different from reality when conducting field test in the environment. To further improve the accuracy of ray tracing simulations, an accurate relative permittivity extraction method together with a timely updated 3D environment model is needed.

Model Lightweighting and Computational Complexity.
OAP 3D environment models consist of large amounts of data, which massively adds to the ray tracing computational complexity. The problem is particularly serious when dealing with higher-order field components. Therefore, the environment model should be light-weighted before implementing in the ray tracing simulator.
To evaluate the model lightweighting, we first governed the computational complexity of the ray tracing algorithm. The ray tracing simulator calculates the RSSI of the propaga-tion paths in two steps. Firstly, the paths from the transmitting to the receiving antennas were literally traced. Then the attenuation suffered by the obstacles in each path was calculated. For large-scale ray tracing applications implementing the IM scheme, every facet or edge in the environment model should be tested. However, only a small amount of the environment model elements actually react to the ray propagation processes. As a result, most of the computational complexity lies in pathfinding rather than path loss calculations. We analyzed the computational complexity for each multipath component by multiplying the number of observation points with their corresponding interacting environment model elements. Assuming that the computational complexities in finding all subpaths are the same as calculating the direct ray, all results refer to the computational complexity of a single direct path as the reference unit (U LOS ). The computational complexities for each reflected and diffracted component were calculated using where C R i and C D i are the computational complexity for the i-th order reflected and diffracted field components, respectively, N P is the number of observation points along the railway track, N f is the total number of triangular facets, and N e is the total number of edges in the environment model. The computational complexity for the direct field (C Direct ) and transmitted field (C Trans ) are both equal to N P since no extra steps should be taken when finding the propagation paths. The calculation results are demonstrated in Table 4. The statistical analysis shows that the    Although the proposed method largely reduces the data volume and promotes calculation efficiency, the implementation of OAP still puts limits on the processed data volume in the light-weighted environment models. OAP established environment models out of a large number of aerial photographs, which described the natural terrain with lower efficiency than the regular artificial constructions. Especially when dealing with vegetated areas, the bushes with leaves often made the 3D models distorted with overlapped triangular facets. For most ray tracing algorithms, vegetated areas are often regarded as simplified 3D surface models with a given empirical parameter that considers the influence of the vegetation [33]. As a result, methods should be taken to the model to extremely simplify those areas before the ray tracing calculations. In this paper, we used the 3D modeling software to manage the model lightweighting and applied higher levels of lightweighting for the natural terrain. However, over-simplification of the environment models led to severe model distortions, which lowered the accuracy of the simulations. Especially when the signal propagation was obstructed by vegetated natural terrain, the model distortion and breaking points led to severe errors in calculating diffracted field components. In order to avoid serious simulation degradation, the natural terrain model was described using excessive triangular facets than necessary after the lightweighting. That was the reason the lightweighted model still contains over 0.1% of the model vertices to describe the above-mentioned scenarios.
Compounding different modeling techniques to establish the environment model and choosing a more efficient way to express the natural terrain is a promising way to solve the problem and accelerate the ray tracing process one step further. For example, we can implement OAP to describe the artificial constructions and LiDAR to describe the natural terrain. Such compounding techniques take advantage of different reconnaissance technologies and further lightweight the environment model for ray tracing usage.
3.4. Ray Tracing Strategy. Outdoor railway scenarios are often open-air applications with the transmitting antenna above the roof of trains. In our selected railway scenarios, the observation points along the railway tracks in scenarios 1 and 2 are mostly in good line-of-sight (LOS) with the signal towers. As a result, higher-order multipath components seldom appear, and propagation predictions considering merely low-order multipath components should achieve fairly high accuracy. Scenario 3 is a NLOS dominating area. In such areas, natural terrain structures are usually large in size, which makes the higher-order diffractions experience large path loss comparing with the low-order components. As a result, high-order multipath components do not contribute much to the RSSI results in railway scenarios. Therefore, it is a straightforward way to ignore the second-order calculations and accelerate the ray tracing algorithm with slightly compromised prediction accuracy.
To support the above-mentioned ray tracing strategy, we summed the absolute value of the RSSIs for each separated path with their corresponding multipath components in all observation points, the results of all three simulation processes are demonstrated in Table 5. Antenna A and antenna B were in scenario 1, which was a LOS dominating scenario. According to the statistical analysis results, the direct field and the first-order reflected field consisted of most of the received signal. In this scenario, the first-order diffraction  only occurred occasionally and the second-order components did not exist. Meanwhile, the situation for antenna C (Scenario 2) was much more complex. Since there were many non-line-of-sight observation points, the proportion of the transmitted and diffracted field components increased dramatically with appearing second-order components. However, the second-order components added together contributed to 0.04% of the total received signal. Scenario 3 is a special case for railways. The overall scenario is still an open area, which is similar to Scenario 1. However, the high-speed railway track experienced a curved line amount the valley, which resulted in a high NLOS coverage rate for wireless communication since the LOS transmission for the wireless signal was blocked by the surrounding environment to a great extent. According to the result, 61.29% and 78.98% of the wireless signal LOS transmission was blocked, respectively, for the two antennas in this scenario, and the two antennas showed different multipath component distributions. Antenna D was located close to the railway track, so the summed directed field was strong even under the circumstances that most areas in this scenario were NLOS areas. In this situation, the transmitted field, as well as the first-order reflected field, dominated the NLOS areas. In the meantime, antenna E shows a different situation. When the antenna is distant from the railway track, most NLOS areas were covered with the diffracted field. Judging from all the above results, the sum of the direct, transmitted, first-order reflected, and first-order diffracted field components account for over 99% of the total received signal in this railway scenario. Furthermore, the MSE and MAE comparison results for all radio coverages when included and excluded the second-order components in the simulations are demonstrated in Table 6. The errors for all curves shows only slightly differences when excluding the 2 nd -order calculations, indicating that implementing loworder multipath components in ray tracing can still accurately describe the radio propagation in the selected railway scenarios.

Conclusions
In this paper, we proposed a method to implement highly accurate 3D railway environment models on the ray tracing simulator and make radio propagation predictions along the railways. To adapt the 3D environment model to the ray tracing algorithm, we decomposed the environment models obtained by OAP into submodels for artificial construction and natural terrain, respectively. Different simplification levels were implemented on different submodels to lightweight environment models before they were imported to the ray tracing simulator. Simulations and real-world experiments were conducted to test the proposed method. The results show that a 4 dB RMSE in the selected railway scenarios was achieved with comparably low computational complexity. It was also found that the low-order multipath components account for most of the received radio signals in the selected railway scenarios. By focusing ray tracing on those low-order components regardless of the highorder ones, the algorithm can be further accelerated.

Data Availability
The datasets used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.