Next Article in Journal
Automatic Detection Method for Cancer Cell Nucleus Image Based on Deep-Learning Analysis and Color Layer Signature Analysis Algorithm
Next Article in Special Issue
Vibroarthrographic Signal Spectral Features in 5-Class Knee Joint Classification
Previous Article in Journal
Transduction Mechanisms, Micro-Structuring Techniques, and Applications of Electronic Skin Pressure Sensors: A Review of Recent Advances
Previous Article in Special Issue
Distribution Transformer Parameters Detection Based on Low-Frequency Noise, Machine Learning Methods, and Evolutionary Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Possibility of Identifying a Complex Polygonal Tram Track Layout Using Multiple Satellite Measurements

1
Department of Electrified Transportation, Gdańsk University of Technology, Gabriela Narutowicza 11-12, 80-233 Gdańsk, Poland
2
Department of Geodesy and Oceanography, Gdynia Maritime University, Morska 81-87, 81-225 Gdynia, Poland
3
Department of Rail Transportation and Bridges, Gdańsk University of Technology, Gabriela Narutowicza 11-12, 80-233 Gdańsk, Poland
4
Department of Transport and Logistics, Gdynia Maritime University, Morska 81-87, 81-225 Gdynia, Poland
5
Department of Geodesy, Gdańsk University of Technology, Gabriela Narutowicza 11-12, 80-233 Gdańsk, Poland
6
Department of Navigation and Marine Hydrography, Polish Naval Academy, Inżyniera Jana Śmidowicza 69, 81-103 Gdynia, Poland
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(16), 4408; https://doi.org/10.3390/s20164408
Submission received: 7 July 2020 / Revised: 28 July 2020 / Accepted: 5 August 2020 / Published: 7 August 2020

Abstract

:
We present the main assumptions about the algorithmization of the analysis of measurement data recorded in mobile satellite measurements. The research team from the Gdańsk University of Technology and the Maritime University in Gdynia, as part of a research project conducted in cooperation with PKP PLK (Polish Railway Infrastructure Manager), developed algorithms supporting the identification and assessment of track axis layout. This article presents selected issues concerning the identification of a tramway line’s axis system. For this purpose, the supporting algorithm was developed and measurement data recorded using Global Navigation Satellite System (GNSS) techniques was evaluated and analyzed. The discussed algorithm identifies main track directions from multi-device data and repeated position recordings. In order to observe the influence of crucial factors, the investigated route was carefully selected. The chosen tramway track was characterized by its location in various field conditions and a diversified and complex geometric layout. The analysis of the obtained results was focused on the assessment of the signal’s dispersion and repeatability using residuals in relation to the estimated track’s direction. The presented methodology is intended to support railway infrastructure management processes, mainly in planning and maintenance through an efficient inventory of the infrastructure in service.

1. Introduction

In recent years, the managers of rail infrastructure have been increasingly recreating the course of railroads, as part of the overall inventory, in the global coordinate system. Initially, the geometric system in the global respect—the concept of Absolute Track Geometry (ATG)—was applied only to selected lines, which were characterized by specific operating conditions, such as railway lines of high operational priority (high speed rail lines) or sections coexisting with other structures of transport infrastructure, such as bridges, tunnels or squares. Other lines were measured using conventional measurement methods [1,2]. This was due to the relatively large labor input associated with the measurement and reproduction of the coordinates of the track courses in the global systems of spatial references. At present, in the processes of the modernization of railway lines, railway geodetic control networks are mounted more and more often, which is related to the state of geodetic networks. They allow the effective reconstruction of the track axis in the global system by tacheometric measurement (backward cutting to several reference points).
The real revolution in the field of track axis reconstruction was brought about by mobile surveying methods, which allow a continuous measurement at a higher speed than in the case of static or quasi static tacheometry. Due to techniques such as Global Navigation Satellite System (GNSS), Inertial Navigation System (INS) and Mobile Laser Scanning (MLS), the effective inventory of all lines at the railway network scale is now possible [3,4,5,6,7,8]. Moreover, these techniques enable a mobile measurement of coordinates with a high density of data recording at a relatively high speed of the measurement process. In extreme cases, it can even be the measurement with the speed of regular trains [9]. Contemporary methods of data processing are also worth mentioning. With frequencies of data recording of tens of hertz, it is necessary to use efficient algorithms for processing the measurement signal. The signal achieved is the so-called point cloud, which without further interpretation does not contribute valuable information. Hence, the signals should be interpreted in such a way that the data set can be converted into the model representation of geometric systems of the track, mainly in the horizontal plane (straight sections, transition curves and circular arcs), and also in the vertical alignment (sections of a constant slope and sections in vertical arcs) [10,11,12]. Because of the large number of collected coordinates, the effective analysis must include methods of automatic selection of particular sections of geometric systems [13,14,15,16].
Until recently, a questionable issue was the accuracy of GNSS measurements. The methods using tacheometers and railway control networks guarantee an accuracy of up to several millimeters. On the other hand, until recently, such accuracy in GNSS measurements was difficult to achieve. At present, it has been reported that an accuracy of about a 1 cm level is achievable also in this survey technique. Generally, the positions along the track are obtained based on the conversion of measurement signals from the different devices, based on the data fusion algorithms [17,18]. Furthermore, the accuracy often analyzed relates to certain specific measurement conditions, but railroads’ conditions are varied.
For more than 10 years, the authors have been developing the technique of satellite measurements in order to evaluate implementation possibilities of this method [19,20]. At present, as part of the project InnoSatTrack shared with the Polish Railways, current possibilities of using GNSS, INS and MLS systems are subject to assessment for the purposes of track axis inventory (position and shape) [21,22] as well as for train positioning possibilities [23]. In the works under way, the biggest accent is put on GNSS technology, because the positioning of devices with high accuracy is the key to reconstruct the track axis in the global coordinates system [24]. In turn, measurements using additional devices, such as marine satellite compasses, inertial devices or laser scanners perform auxiliary functions, in particular for the purpose of determining corrections on the basis of additional information about the movement of measurement platform relative to the track [25,26].
In this paper we present the assessment of the possibility of determining the main directions of a polygon system, based on the example of a complex geometric system of the track. The main directions of the route are defined as tangents to horizontal curves in their boundaries with straight sections. It could be assumed that in specific cases the main direction covers a straight track section or a tangent between reverse curves. According to the above definition, the polygon of the route is understood as a system of main directions defined in the Cartesian coordinate system.
Geometric layouts in railway lines are usually characterized by mild changes of curvature (long straight lines and transition curves, large radii). For this reason, the authors tested both the measurement methods and data analysis methodology in relation to layouts of more diverse geometry. The examples of such systems are railway lines of agglomeration importance and mainly tramway lines, which in conditions of urban development require the use of short geometric sections [27,28]. In such situations, the exact reconstruction of the route’s polygon is a prerequisite for the proper reconstruction of the system of non-linear sections by using analytical methods of design [29,30,31,32].
An additional advantage of the analysis of the system on the tramway network is a varied degree of the occurrence of field screens. In such a situation the assessment of the accuracy of the reconstruction of the polygon is more reliable and indicates the potential possibilities for planning measurement campaigns on the lines of rail transport.

2. Materials and Methods

2.1. Characteristics of the Measured Tramway Route

For the analysis of the main directions, a tramway line including a street terminus in the district of Gdańsk Brzeźno (coastal town in northern Poland) was selected. Choosing the route, the following assumptions were taken into account:
  • The route should be characterized by diversified visibility conditions of GNSS satellites, in order to verify and assess the impact of the number of visible satellites (their distribution on the horizon) on the possibility of the reconstruction of track geometry;
  • The measuring section should have a diversified condition of track (structural and geometrical) so as to assess the impact of the state on the reconstruction quality;
  • The measuring section should allow multiple measurements, in order to assess the reproducibility of the results obtained related to the parameters of the identified geometric system of the axis of track.
The measurement campaign was carried out as a part of the project InnoSatTrack from 28.11.2018 until 29.11.2018 on a section of the tramway network in Gdańsk [22,24]. A large street terminus in Brzeźno district was selected (shown in Figure 1). The length of the test section was approximately 3 km, allowing multiple measurements of the same section overnight. Additionally, the system of the route in the form of a street terminus assured that there was no need to perform any maneuvers to prolong the duration of measurement sessions. Finally, the selected street terminus was characterized by diverse pavement conditions and the presence of field obstacles along particular sections.
Detailed characteristics of sections of the selected terminus:
Section a:
Section of a length of about 850 m, starting at the terminus Brzeźno Plaża. Double-track tramway route is located on a strip dividing four-lane alleys. The section is characterized by a small number of field screens (single trees near the track at the terminus) and a very good technical condition of the track’s structure. The section was chosen as the reference section in relation to others.
Section b:
Section of a length of about 400 m. One-track tramway route is located asymmetrically on the north side of a four-lane avenue. The fragment is characterized by a small number of field screens (single trees near the track), and a very good technical condition of the track’s structure. This section was also chosen as a reference section.
Section c:
Section of a length of about 1000 m. One-track tramway route is located asymmetrically on the west side of a two-lane street. This section is characterized by a larger number of field screens than in previous sections (groups of trees and buildings located along the track) and very poor technical condition of the pavement (high values of twist and irregularities of the track in both the horizontal and vertical planes). This section was selected as the section to study the effect of track conditions on the accuracy of measurement.
Section d:
This section is the tramway route carried out independently of the road system with a length of approximately 150 m. This track is characterized by a complex geometry and the presence of the most extensive field screens (trees with crowns located over the streetcar track), and at the same time with a very good quality of the track. It was chosen to assess the impact of field screens on the possibility to reconstruct a complex geometric system of the track (field obstacles are shown in Figure 2). This section includes a diverging track to the Brzeźno terminus.
Section e:
Section of a length of about 330 m. One-track tramway route located asymmetrically on the south side of an access one-lane road. It is characterized by the appearance of very extensive field screens (tall buildings adjacent to the track) and the good condition of the track. It was chosen to assess the impact of field screens on the possibility of the reconstruction of a typical geometric system of the track (see Figure 3)
Section f:
This section is independent of the road system, and its length is about 270 m. It is characterized by an average number of field screens and good track condition. This section was chosen as a reference for section c.
Figure 2 shows the track in sections a and b, while Figure 3 shows examples of locations on sections c, d, e and f. Sections a and b are characterized by rather good conditions due to the visibility to satellites, while sections cf are located in close proximity to high-rise buildings or trees. Here, disturbances in satellite signals from the sections of the northern edge of the terminus are expected.

2.2. Characteristics of Measurement Unit

To carry out the measurement, 10 GNSS receivers installed on a measurement platform were used. For the investigation, Leica GS18 and Trimble R10 receivers were selected. The corrections for real time kinematic measurement were obtained from the closest reference stations of the SmartNet and VRSNet networks for the Leica and Trimble receivers, respectively.
As measurement cars, the bogies of the GWF 300 car series were adopted. They are characterized by favorable conditions for matching their ride axes into the track axis (no extensive amortization, no clearances) and a relatively low height of the measurement deck. On each of both bogies, five receivers were installed in a rectangular formation. In the centers of both bogies, directly above the bogie’s pivot, central antennas were placed. A diagram of the layout of the receivers is shown in Figure 4.
As the tramway bogies do not have their own drive device, there is a problem concerning the proper selection of a motor vehicle. The vehicle should meet the following requirements:
  • The vehicle should provide both constant and low velocity of about 30 km/h. Only relatively modern tramways, equipped with power electronics drive control fulfill this requirement;
  • The vehicle should provide a 230 V AC power supply to measuring devices and recording computers from the onboard network (preferably equipped with a standard 230 V socket).
The vehicle that met the above-mentioned requirements proved to be a Bombardier NGT-6/2. Finally, three measurement bogies were attached to the engine car. On the first bogie (the closest to tramway) the devices recording additional parameters of a run (such as acceleration, inclination and direction of the run) were installed. On the other two, a set of 10 GNSS receivers were placed. To avoid a negative influence from naturally occurring obstacles, i.e., tramway cabs, an additional bogie was located between the tramway and the bogies with GNSS antennas. A diagram of the measuring train is shown in Figure 5.
In the present measurement methodology, the fact is taken into account, that the position of the axis of the track relative to the rail vehicle is not unambiguously determined. One must keep in mind that each rail vehicle has a small degree of freedom in the lateral direction due to the bogie structural clearances and clearances between the rails and wheels. These clearances ensure the correct fitting of the bogie’s axles into the track. Moreover, besides of the constructional clearances of the bogies, displacements due to the spring amortization of the rolling stock occur. Taking into account also conicity in a rail and wheel interface, some kind of typical oscillation (so-called hunting oscillation) should be expected. However, this kind of movement is hard to assess without additional reference measurements (for example vision methods) and is not discussed in detail in this article.
It is therefore assumed, that the most reliable results with respect to the track axis are obtained by situating the antenna in the axis of symmetry with respect to the construction of the wheel set. This point for the construction of the bogie is in the pivot axis. This operation is difficult due to the limited access to the structural elements of pivot from the outside. Thus, there is the first factor related to the uncertainty of measurement, i.e., determination of the precise location of the base receiver on the vehicle’s platform. Due to the above problem, the positioning of the GNSS receivers’ centers was performed in conditions that minimize positioning errors. Installation of equipment was carried out in the hall of the tram depot. The process of measuring devices installation is shown in Figure 6. Special support structures enabled the precise installation of each tribrach according to the scheme assumed. Figure 7 presents the bogies instrumented with controllers and receivers during measurements.
During the measurement campaign, six passages on the specified street terminus were performed between 0:09 and 3:37. At this time, 10 receivers operated in continuous recording mode at a frequency of 20 Hz. During the study, a total of 1,186,518 positions were collected. All of the records met the assumed accuracy of designation (uncertainty radius not greater than 0.2 m).

2.3. Methodology for Determining the Main Directions of the Routes

The axis of the railway or tramway track physically does not have any representation. Considering a typical cross-section of a track, its axis is located at a point distant by a length of half the width (track gauge) between the rail courses at a height of the line connecting running surfaces of rail heads. In general, the axis of the track does not have to coincide with the central axis (so-called centerline) of a railway track. This is the case in curves with small radii, in which a widening of the track is often made. Figure 8 shows the definition of the axis of a railway track and its centerline.
Reconstruction of the track position in straight sections does not require additional adjustments related to the transverse or longitudinal inclinations. It is assumed that the track in the straight section is designed and built as a track without a lateral inclination (cant). Therefore, any found deformations resulting from the technical degradation of the track will be removed from the signal if necessary. Thus, the main directions will be assessed based on the analysis of the recorded positions on the straight sections of the track in particular passing measurements.
Prior to the analysis, the following assumptions were made:
  • Multiple repeated measurement of a precisely indicated track section will enable to evaluate the error of the identification of the main directions. The directions will be defined on the grid of horizontal coordinates;
  • Generally, the main directions will be designated by antennas positioned in the pivot axis of the measurement bogies;
  • The geometric layout of the selected test section will be identified as a system of straight lines, regardless of the nature of deformation occurring in the track;
  • Obvious deformations observed in the measurement signal, such as the so-called broken straights, will be removed, as they distort the assessment of the main directions;
  • It was assumed that the assessment of the accuracy of the identification in the horizontal coordinate system will be made on the basis of the analysis of the reproduction of straight sections signals.

2.3.1. Data

The selected straights located in the designated sections of the tramway terminus were analyzed. Figure 9 shows the location of individual straights in the form of coordinates on the grid of the PL-2000 system in force in Poland (Gauss–Krüger projection for Geodetic Reference System ‘80). These straights were separated from GNSS measurement data based on the visual assessment of the results. Table 1 shows the signatures of straight lines in the corresponding sections of the tested track. Data in the form of coordinates and additional information, such as the accuracy recorded by the particular receiver as well as the number of satellites used to calculate the position at a given moment of time, were used in the analysis. Table 2 shows an example of the scope of measurement data.

2.3.2. Main Directions

The main direction of the existing railway track is defined as a straight line that describes the rectilinear section of the track best or as a tangent to the reverse curves system at the point of their contact. The main directions form a system called a route’s polygon.
As defined, the main direction is described by the parameters of a linear model equation of the straight line. In this analysis, these parameters were estimated by the least squares method. Suitable sets of points were selected for the analysis with Y X coordinates in the PL-2000 system, which are the measured positions of GNSS receiver.
Signals from 10 receivers recorded during six measurement runs were analyzed. It can therefore be assumed, that each of the 60 signals is an independent description of the location of the analyzed section of track. Obviously, the most reliable signals are the signals recorded by antennas situated over the axis. In these places, transverse shifts of the tram bogies in relation to the track should be minimal. The antennas located at the corners of measuring platforms should also draw straight lines with the same direction coefficient. However, these positions can be subject to increased uncertainty due to the shifts and vibrations of the platform during test runs.
The signals recorded by individual antennas in each of the runs were analyzed in the local systems of coordinates, of which the abscissa axis is in accordance with the least squares line. This methodology allows for a clear interpretation of the differences between the signals recorded. A relatively large number of points used in the regression calculation makes it so that the straight lines are determined at a very low matching error, especially in situations of long straight lines.
The methodology of transformation that uses homogeneous coordinates was accepted. The transformation of the Y X coordinates of the PL-2000 system to homogeneous coordinates can be written as:
( Y , X ) [ Y X 1 ]
Thus, transformations of coordinates can be described with a matrix equation in the form:
[ Y 1 X 1 1 Y 2 X 2 1 Y n X n 1 ] = [ Y 1 X 1 1 Y 2 X 2 1 Y n X n 1 ] · M
where M is a matrix of transformation in the form:
M = [ m 11 m 12 1 m 21 m 22 1 m 31 m 32 1 ]
The elements of matrix M ( m ij , i = 1 , 2 , 3 , j = 1 , 2 ) depend on a character of transformation and in specific cases of a single transformation this matrix takes the form of a translation matrix T ( t h , t v ) , with rotation R ( φ ) or conversion of symmetry F performed with respect to the O x or O y axes of coordinates:
T ( t h , t v ) = [ 1 0 0 0 1 0 t h t v 1 ] ;   R ( φ ) = [ cos φ sin ( φ ) 0 sin ( φ ) cos φ 0 0 0 1 ] ;   F = [ ± 1 0 0 0 ± 1 0 0 0 1 ]
where:   t h , t v are the distance from the origin to a translated point along the horizontal and vertical directions, and φ is the rotation angle of the vector of the transformed point with respect to the origin of the coordinate system.

2.4. Automatic Identification of the Polygon-Algorithm for Identifying Main Directions

The geometric layout of the railway in the horizontal plane consists of straight-line sections and the sections arranged in curves, consisting of circular arcs and transition curves with variable curvature along the length of the elements. The analysis presented in this paper concerns the polygon, which is formed by the main directions in accordance with the straight sections of the track or which are tangent to reverse curves. It is also assumed that two successive directions in places of their intersections are formed by the vertices of horizontal curves. The angle through which the direction changes from direction i to i + 1 is the intersection angle.
For the description of polygon system, the equations of the straights of the Y X coordinates (in accordance with the PL-2000 system) as X = a m d · Y + b m d was used. Parameters of the equation ( a m d , b m d ) refer to the main directions of the analyzed route.
The basic steps of the algorithm supporting the identification of main directions’ system are described below.
(1) The first step is the discretization of the data to a finite number of sections. Fragmentation of the measurement signal took place on the basis of assigning each coordinate to intervals of a fixed length. At this stage the geometric system is not known, so the road is calculated by summing up distances from point to point. This is the first estimate of the distance along the route, which—when carried out in an appropriate way (filtering out concentrations of points recorded at standstills etc.)—will enable further analyses to identify the polygon. Figure 10 shows a schematic diagram of the route discretization, while Equation (5) shows a basic method of counting the distance along the route. The length of each section depends on the density of records and it is generally set in the phase of tuning the algorithm to the needs arising from the characteristics of the measurement signal.
D i = D i 1 + d i = D i 1 + ( Y i + 1 Y i ) 2 + ( X i + 1 X i ) 2
where: D i —distance from the starting point to the current one; Y i , X i —flat coordinates of measuring points p i in the local system with directions in accordance with the PL-2000.
In each of the sections, there are a number of points with coordinates Y X . When testing the algorithm, the minimum number of points in a particular interval was fixed; it should be remembered, that in continuous measurements there are places with no records of a measurement signal. Then, for each of the intervals the analysis of regression is made between the Y and X coordinates, due to which the estimate of the directions characterizing individual intervals of discretized route takes place, as shown in Figure 11. Due to the auxiliary character of estimated directions, the calculations were performed according to the least squares algorithm method for linear correlation. The method in this case involves the determination of the parameters of simple regression by minimizing the deviations of the measurement values of vertical coordinate X of the reference system PL-2000 (or the local system in accordance with the directions compatible with PL-2000) with respect to the linear model. The method accepted can be described with the equation:
f ( a j , b j ) = i = 1 m ( X i ( a j Y i + b j ) ) 2 m i n
where: a j , b j —parameters of linear model in section j , and Y i , X i —coordinates of measurement points p i in the local system with directions in accordance with the PL-2000.
Each section has its own estimator of direction, which can be the coefficient of the least squares straight slope or the angle to the axis of abscissae Y . It could also be the azimuth of the designated direction.
(2) With estimators of the direction in the individual sections of the route, variation of this direction is analyzed in the passages from one interval to another, i.e., in the series of consecutive sections. It is assumed that the angle representing the direction in each section maps the variation of the shape of the axis of the route analyzed. On the straight sections, a slight oscillation of the estimator value is expected, and in the areas of curves beside the oscillation a clear increasing trend of this value should appear. In order to examine these trends, the algorithm groups the adjacent sections with one another (of which each has a calculated estimator of direction in the previous step) in a minimum number and carries out the next analysis on the sets thus designated. The number of sections in each group depends on the geometric system (length of homogeneous sections of the route) and is determined by the analysis of Pearson’s correlation coefficient. This time, the analysis of regression concerns the values of the angles of individual sections. If we assume the angle as azimuth ( A z ) then the analysis of the grouped regressions of section can be written as follows:
f ( a g , b g ) = j ( A z j ( a g D m i d j + b g ) ) 2 m i n
where:
  • a g , b g —parameters of a linear model of the group of adjacent sections;
  • A z j —angles characterizing the individual sections in the group;
  • D m i d j —middle values of the distance covered in each section (equivalent to mileage).
Schematic diagram of the analysis of the regression of the sequence of angles indicated in the individual sections of discretized route is shown in Figure 12.
According to the assumptions, if the absolute value of the direction coefficient a g does not exceed the threshold value assumed in a particular group of sections, the algorithm assigns to a particular group (at this distance) the attribute of a straight section. The total length of the straight line is calculated as the sum of the lengths of the next (consecutive) groups of sections characterized by the stability of angle, and thus the direction coefficient of regression did not exceed a certain value close to zero. This approach makes it possible to find the main direction as well in combination with reverse curves. Then, the location with the value of regression close to zero will appear in the analysis.
(3) Once the records belonging to the straight sections are identified, the total, i.e., the linear regression defined on the whole set of points, is determined. Due to the fact that the straight sections can vary in length and for longer sections precision parameters of determining the position can undergo dynamic changes, the authors used a weighted regression analysis, in which the weights are defined by measurement uncertainties recorded for flat coordinates by GNSS receivers. The method implemented in this way maps the main directions on the sets of points with the differentiated distribution of position error more accurately than the simple linear regression (unweighted).
(4) The coefficients of the equations of straights constitute the basis for building the polygon, in which the coordinates of vertices are calculated as points of the intersection of adjacent (by indexation) straight lines. The coordinates of the points of the start and end points of the polygon are calculated based on measured extreme coordinates. A schematic diagram of the reconstructed system of vertices is shown in Figure 13.
(5) The last step of algorithm is the analysis of the slopes of the adjacent straight lines. In this step, it is possible to merge the adjacent straight lines with one another. This is the moment in which the two adjacent straight lines with different direction slopes should be accepted (or not) as the area of the change of main directions. The adjacent two straight lines can therefore be designated as two tangents of a horizontal curve or connected together in one tangent (main direction) of the polygon.

3. Results

For detailed analyses, the data obtained from receivers number 5 and 10 were selected. These receivers were mounted on the pivots of the tramway bogies, which means that the trajectory recorded by them in the horizontal plane maps the axis of the track. The exceptions are sections located in the curves with a cant, but the test section was selected in such a way that the impact of the cant is as small as possible. These receivers recorded a total number of 253,435 coordinates during the passage through the test section.
Then, the analysis of the resulting measurement signal was made, according to the algorithm shown in Section 2.4. For input data in the form of the measurement signal recorded by two receivers (with a frequency of 20 Hz) during the six passages of the route, the algorithm detected the presence of 536 straight sections, giving an average of about 45 straight sections making up the whole geometrical system. Figure 14, Figure 15 and Figure 16 show the calculation results of the algorithm. They represent scalable windows of a program for data analysis, implemented in a Scilab environment. The automatically developed computer program defined the beginnings and ends of the straight sections, and in half the reconstructed section a label with the next number of section was pasted. Since the algorithm operated simultaneously on all 12 measuring signals, these fragments, which a number of straight lines overlap, can be treated—as the way of visual assessment of the result in the graph—as repeatedly reconstructed. If the numbers of straight lines are offset with respect to one another, it means that the algorithm was not able to reconstruct the beginnings and ends of straight or divided sections into shorter ones at a high level of repeatability, or indicated that the straight lines as false positive, recognizing, for example, the distortion of a signal or inequality of a track in a straight section as a curvilinear section.
It can be noted that on sections a and b, the algorithm identified the vertices most accurately. Figure 15 shows the vertices determined on the basis of identified straight lines in sections b and c. Section c compared to section b featured less precise characteristics of the straight lines reconstructed. The result was a clear dispersion of the designations of vertex in the whole series of measurements. In addition, Figure 14 shows that on the other sections of the test track (mostly c, d, e and f) a series of very short straight sections were measured, hindering the overall analysis (expressed as a loss of clarity of description at the scale assumed). In addition, in section c there were fragments of the track with very short reverse curves (shown in Figure 16), which constituted an additional problem for the automatic identification of the algorithm.

4. Discussion

4.1. Analysis of Main Directions

Bearing in mind that in the case of tramway lines, deformations of the track axis in the horizontal plane can be of similar length to the short straight sections—full automation turned out to be very difficult. In this situation, for the further analysis of the main directions only those sections were selected for which there was little doubt that they constitute the main directions of the route.
Among the sections identified by the algorithm as straight sections, a, b and c were selected for further analysis of eight straight of sections. Figure 17 shows the final selected straights and straights selected initially (on the basis of a visual evaluation of the measurement data). To distinguish the straight lines from the earlier selected sections (Figure 9), for the first ones Str’ markings were introduced.
Comparing both sets of straight lines, it can be concluded that there was a large convergence of both methods. There were significant differences in the division of Str 6 to Str’ 5 and Str’ 6 and division of Str 7 to Str’ 7 and Str’ 8. It is worth noting that the reconstructions of straight sections of relatively short length and the longest ones were different. In addition, existing GNSS measurement uncertainties affect the repeatability of the reconstruction of very short sections. However, the uncertainty of designating the angle in the case of short sections will cause small deviations of the model course of the straights with respect to the actual course. However, in the case of long sections, even a small error in the designation of angle can have a significant effect on the results obtained at the ends of the reproduced geometric layout. That is why the indicator taking into account the length of a section was introduced. This indicator, designated as σ A z ( r e l ) , is the standard deviation of the reconstructed azimuth multiplied by the length of section. The obtained parameters of the sections examined are shown in Table 3 and Table 4.
For each measurement, a set of coordinates was determined. Then, the mean and standard deviation of the corresponding values were calculated. Based on the coordinates of vertices the distances | W | between these vertices were calculated (of the individual measurements i = 1 , , n ) as well as the averaged vertex X a v ,   Y a v . The measure of uncertainty determined the vertex based on the straights reconstructed were therefore the distance | W | a v and its standard deviation σ | W | . All parameters describing the vertices, together with intersection angles α , are shown in Table 4.
The values in Table 1 and Table 2 clearly show that Section 5 and the associated vertex W5–6 were reproduced with the highest uncertainty. Section 5 is a straight insert between the crossover needle (the measurement train moved in the reverse crossover direction) and the curve behind the crossover. This section is very short, and the intersection angle at this point is very small (less than 1 degree). The measurement on such a short section is burdened with considerable uncertainty to designate its direction, and also coordinates of vertices W4–5 and W5–6. Additionally, the vertices analyzed featured by a very small intersection angle, comparable to the value with the uncertainty to determine the angle. The coincidence of a short length of Section 5 and the small intersection angle in vertex W5–6 generated a very high level of uncertainty (mean value of error and standard deviation at 1.3 m). However, this was not the result of an incorrectly accepted measurement technique, but of geometric relationships in the measured geometric system. The coordinates of vertices with the favorable system of tangents were determined with very high reproducibility.

4.2. Accuracy of Measurement Based on the Repeatability of Signals

All originally selected straight sections were analyzed. The existence of numerous unfavorable deformations of the track and the field conditions for GNSS measurement should appear in the analysis of signal reproducibility. Figure 18 shows the results of all measurement series—for each straight section—recorded by one of the antennas positioned on the pivot (receiver 5). In most cases, a very good repeatability of the trajectory of this device was observed. It is clear that in cases of bends on the straight lines all runs recorded the same qualitative character of deformation. This observation makes it possible to trim signals so that they do not distort the precision assessment. Of course, cutting out only the deformed parts from the straight line disturbs the actual direction of the straight line; however, to assess the accuracy of the measurement itself such an operation is necessary. Figure 19 shows the results recorded by all receivers during the measurements of section Str 9. This section features the poor track conditions. Numerous irregularities and deformations caused the occurrence of vibrations on the measurement vehicle. Graphs of recorded signals indicated a large number of bends (deformations) along the section of this track. This is one of the sections where the algorithm of the detection of main directions selected a series of short straights. The signals recorded by all receivers revealed the same deformation of the track.
Figure 20 and Figure 21 show the results of the analysis of the reproducibility of the data recorded by receivers 5 and 10 for sections Str 12 and Str 6. As seen in Table 1, these sections were characterized by a very unfavorable environment in the form of buildings and trees. This resulted in a discontinuity in data recording during part of the particular run. The analysis was done for a short segment that was recorded during all six crossings. As regards the accuracy in the reconstruction of direction, the standard deviation was at a level of 1–2 cm, depending on the recorded segment of a particular section, but reliability of measurement under these conditions was relatively low—every crossing showed disturbances of signal loss by different receivers.
The example of section Str 6 (Figure 21) shows that in the accepted methodology of measurement, a very high reproducibility of independent measurements can be achieved, which directly proves a high precision of the reconstruction of the main directions of the tested route. However, the overall assessment ultimately depends on the field conditions in which the measurement is made. Systems based on several receivers operating independently increase the reliability. Therefore, even in difficult measurement conditions, data collection is still possible. However, the scope can be limited to certain sections of the route.
Figure 22 shows a graph of least squares lines in the local coordinate system for section Str 6. Also in this case, reference straights were brought to the local system, which facilitated the assessment of mutual relationships with individual straights. The black line represents the estimator of the expected value, that is, the real position of the track (reduced to level 0 of the local system). The lines in other colors are six estimated straight lines from the six individual runs. The graphs show that the maximum values of the straight lines most distant from one another did not exceed 1 cm.

5. Conclusions

The selected test track met the assumptions as regards different conditions of satellite measurements. The analysis of the results clearly showed the impact of urban infrastructure on the possibilities of the reproduction of main directions along tramway lines. The algorithm presented using the regression analysis of measured points (representing the track axis) with the appropriate approach allowed efficient determination of the position of straight sections in a global system of references. However, the specification of geometric systems of tramway tracks (short geometric elements) led to the analysis of the main directions requiring a sufficiently high density of data. In the case of a large number of points, which are burdened with a high level of uncertainty, there is a risk of incorrect estimation of the identified parameters of the route directions, which in turn leads to the incorrect reconstruction of the polygon. The biggest problems with the polygon reconstruction are the intersection angles with very short tangents. The occurrence of small values of intersection angles (about 3 degrees or less) generates a relatively large uncertainty of the reproduction of the intersection of tangents (vertices). Finally, the coincidence of small intersection angles with short tangents makes that even a small deviation of the direction of the tangent causes very significant differences in the coordinates of the polygon’s vertices.
As shown, the tramway infrastructure is characterized by an extremely unfavorable shape from the point of view of the reproduction of geometric systems in the global system. It is dominated by very short-length straight sections and curves with a complex system. In the measurements made, the shortest uniquely located straight section had a length of 0.7 m. It was on the extension of the main track in the crossing (straight section inserted within the turnout construction). Since railway systems have much longer straight sections, a successful measurement made on the tramway tracks is a positive prognostic for further research related essentially to the inventory of railway lines.
The analysis also showed that the geometric condition of track (rails’ deformations and geometric imperfections) has a significant impact on the identification of the main directions of the route. One of the straight lines in the street terminus analyzed featured a high degree of geometric deformation, with the result that the algorithm of the identification of the main directions did not distinguish real directions (design) from the short straight lines resulting from the deformations present. The automatic detection of straight sections in such conditions turned out to be very difficult, requiring more advanced computational algorithms. Sections characterized by a good geometric condition were identified properly and the polygon designated on this basis was characterized by a relatively small error of the reproduction of vertices, at the level of a few centimeters.

Author Contributions

Conceptualization, A.W., C.S., W.K., K.K., J.S. (Jacek Szmagliński), and P.C.; data curation, C.S., J.S. (Jacek Szmagliński), P.C., and P.D.; formal analysis, J.S. (Jacek Szmagliński). and P.C.; investigation, A.W., C.S., W.K., K.K., J.S. (Jacek Skibicki), J.S. (Jacek Szmagliński), P.C., P.D., M.S. (Mariusz Specht), M.Z., and S.J., M.S. (Marcin Skóra) and S.G.; methodology, A.W., C.S., W.K., K.K., J.S. (Jacek Szmagliński), P.C., P.D., and M.Z.; resources, A.W., C.S., and K.K.; software, J.S. (Jacek Szmagliński) and P.C.; supervision, W.K.; validation, J.S. (Jacek Skibicki); visualization J.S. (Jacek Szmagliński) and P.C.; writing—original draft, W.K., J.S. (Jacek Szmagliński), and P.C.; writing—review and editing, W.K. and P.C. All authors have read and agreed to the published version of the manuscript.

Funding

Work were carried out as part of the research project “Development of innovative methods for determining the precise trajectory of rail vehicle” InnoSatTrack (POIR.04.01.01-00-0017/17), financed by NCBiR and PKP PLK SA.

Acknowledgments

The authors wish to thank the Company Gdańskie Autobusy i Tramwaje for making this measurement campaign possible.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. NR (Network Rail). NR/L3/TRK/0030 Iss. 1, Reinstatement of Absolute Track Geometry (WCRL Routes), Network Rail 2008. Available online: https://standards.globalspec.com/std/10216330/nr-l3-trk-0030 (accessed on 5 August 2020).
  2. CEN (European Committee for Standardization). BS EN 13803:2017 Railway Applications–Track–Track Alignment Design Parameters–Track Gauges 1435 mm and Wider. CEN. 2017. Available online: https://shop.bsigroup.com/ProductDetail/?pid=000000000030279796 (accessed on 5 August 2020).
  3. Szwilski, T.B.; Begley, R.; Dailey, P.; Sheng, Z. Determining rail track movement trajectories and alignment using HADGPS. In Proceedings of the AREMA 2003 Annual Conference, Chicago, IL, USA, 5–8 October 2003; pp. 1–10. [Google Scholar]
  4. Szwilski, A.B.; Dailey, P.; Sheng, Z.; Begley, R.D. Employing HADGPS to survey track and monitor movement at curves. In Proceedings of the 8th International Conference “Railway Engineering 2005”, London, UK, 29–30 June 2005. [Google Scholar]
  5. Glaus, R. The Swiss Trolley- A Modular System for Track Surveying. Geodätisch-Geophys. Arb. in der Schweiz. 2006, p. 184. Available online: https://www.sgc.ethz.ch/publications.html (accessed on 5 August 2020).
  6. Akpinar, B.; Gülal, E. Multisensor railway track geometry surveying system. IEEE Trans. Instrum. Meas. 2012, 61, 190–197. [Google Scholar] [CrossRef]
  7. Stein, D.; Spindler, M.; Kuper, J. Rail detection using lidar sensors. Int. J. Sustain. Dev. Plan. 2016, 11, 65–78. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, Q.; Niu, X.; Zuo, L.; Zhang, T.; Xiao, F.; Liu, Y.; Liu, J. A Railway Track Geometry Measuring Trolley System Based on Aided INS. Sensors 2018, 18, 538. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. RILA Track. Available online: https://www.fugro.com/our-services/asset-integrity/raildata/rila-track-rila-360 (accessed on 14 May 2020).
  10. Gikas, V.; Daskalakis, S. Determining rail track axis geometry using satellite and terrestrial geodetic data. Surv. Rev. 2008, 40, 392–405. [Google Scholar] [CrossRef]
  11. Lenda, G. Determining the geometrical parameters of exploited rail track using approximating spline functions. Arch. Civ. Eng. 2014, 60, 295–305. [Google Scholar] [CrossRef] [Green Version]
  12. Li, W.; Pu, H.; Schonfeld, P.; Song, Z.; Zhang, H.; Wang, L.; Wang, J.; Peng, X.; Peng, L. A method for automatically recreating the horizontal alignment geometry of existing railways. Comput.-Aided. Civ. Inf. 2018, 34, 71–94. [Google Scholar] [CrossRef] [Green Version]
  13. Elberink, S.O.; Khoshelham, K. Automatic Extraction of Railroad Centerlines from Mobile Laser Scanning Data. Remote Sens. 2015, 7, 5565–5583. [Google Scholar] [CrossRef] [Green Version]
  14. Singh, K.; Swarup, A.; Agarwal, A.; Singh, D. Vision based rail track extraction and monitoring through drone imagery. ICT Express 2019, 5, 250–255. [Google Scholar] [CrossRef]
  15. Xiong, Z.; Li, Q.; Mao, Q.; Zou, Q.A. 3D Laser Profiling System for Rail Surface Defect Detection. Sensors 2017, 17, 1791. [Google Scholar] [CrossRef] [Green Version]
  16. Yang, B.; Fang, L. Automated extraction of 3-d railway tracks from mobile laser scanning point clouds. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 4750–4761. [Google Scholar] [CrossRef]
  17. Akpinar, B.; Gülal, E. Railway track geometry determination using adaptive Kalman filtering model. Measurement 2013, 46, 639–645. [Google Scholar] [CrossRef]
  18. Zhou, Y.; Chen, Q.; Niu, X. Kinematic Measurement of the Railway Track Centerline Position by GNSS/INS/Odometer Integration. IEEE Access 2019, 7, 157241–157253. [Google Scholar] [CrossRef]
  19. Koc, W.; Specht, C. Selected problems of determining the course of railway routes by use of GPS network solution. Arch. Transp. 2011, 23, 303–320. [Google Scholar] [CrossRef]
  20. Specht, C.; Koc, W.; Chrostowski, P.; Szmaglinski, J. Accuracy Assessment of Mobile Satellite Measurements Relation to the Geometrical Layout of Rail Tracks. Metrol. Meas. Syst. 2019, 26, 309–321. [Google Scholar] [CrossRef]
  21. Wilk, A.; Specht, C.; Koc, W.; Karwowski, K.; Chrostowski, P.; Szmagliński, J.; Dąbrowski, P.; Specht, M.; Judek, S.; Skibicki, J.; et al. Research project BRIK: Development of an innovative method for determining the precise trajectory of a railway vehicle. Transp. Overv. 2019, 7, 32–47. [Google Scholar] [CrossRef]
  22. Dabrowski, P.; Specht, C.; Koc, W.; Wilk, A.; Czaplewski, K.; Karwowski, K.; Specht, M.; Chrostowski, P.; Szmaglinski, J.; Grulkowski, S. Installation of GNSS receivers on a mobile platform—methodology and measurement aspects. Sci. J. Marit. Univ. Szczec. 2019, 60, 18–26. [Google Scholar] [CrossRef]
  23. Specht, M.; Szmagliński, J.; Specht, C.; Koc, W.; Wilk, A.; Czaplewski, K.; Karwowski, K.; Dąbrowski, P.S.; Chrostowski, P.; Grulkowski, S. Analysis of Positioning Methods Using Global Navigation Satellite Systems (GNSS) in Polish State Railways (PKP). Sci. J. Marit. Univ. Szczec. 2020, 62, 26–35. [Google Scholar] [CrossRef]
  24. Specht, M.; Specht, C.; Wilk, A.; Koc, W.; Smolarek, L.; Czaplewski, K.; Karwowski, K.; Dąbrowski, P.S.; Skibicki, J.; Chrostowski, P.; et al. Testing the Positioning Accuracy of GNSS Solutions during the Tramway Track Mobile Satellite Measurements in Diverse Urban Signal Reception Conditions. Energies 2020, 13, 3646. [Google Scholar] [CrossRef]
  25. Koc, W.; Specht, C.; Szmaglinski, J.; Chrostowski, P. A method for determination and compensation of a cant influence in a track centerline identification using GNSS methods and inertial measurement. Appl. Sci. 2019, 9, 4347. [Google Scholar] [CrossRef] [Green Version]
  26. Dąbrowski, P.S.; Specht, C.; Felski, A.; Koc, W.; Wilk, A.; Czaplewski, K.; Karwowski, K.; Jaskólski, K.; Specht, M.; Chrostowski, P.; et al. The Accuracy of a Marine Satellite Compass under Terrestrial Urban Conditions. J. Mar. Sci. 2020, 8, 18. [Google Scholar] [CrossRef] [Green Version]
  27. Specht, C.; Chrostowski, P.; Dąbrowski, P.; Szmagliński, J.; Specht, M.; Dera, M.; Koc, W.; Skóra, M. Mobile satellite measurements on the Pomeranian Metropolitan Railway. Transp. Overv. 2016, 5, 24–35. [Google Scholar] [CrossRef]
  28. Specht, C.; Koc, W.; Chrostowski, P.; Szmaglinski, J. The Analysis of Tram Tracks Geometric Layout Based on Mobile Satellite Measurements. Urban Rail Transit 2017, 3, 214–226. [Google Scholar] [CrossRef] [Green Version]
  29. Koc, W. Design of rail-track geometric systems by satellite measurement. J. Transp. Eng. 2012, 138, 114–122. [Google Scholar] [CrossRef]
  30. Koc, W. Analytical method of modelling the geometric system of communication route. Math. Probl. Eng. 2014, 2014, 1–13. [Google Scholar] [CrossRef]
  31. Koc, W.; Chrostowski, P. Computer-aided design of railroad horizontal arc areas in adapting to satellite measurements. J. Transp. Eng. 2014, 140, 1–8. [Google Scholar] [CrossRef]
  32. Chrostowski, P.; Koc, W.; Palikowska, K. Prospects in elongation of railway transition curves. P I Civil Eng-Transp. 2017, 1–12. [Google Scholar] [CrossRef]
Figure 1. Measuring testing ground based on Google Maps.
Figure 1. Measuring testing ground based on Google Maps.
Sensors 20 04408 g001
Figure 2. View of the streetcar track in the area of sections a and b.
Figure 2. View of the streetcar track in the area of sections a and b.
Sensors 20 04408 g002
Figure 3. View of the tramway track in the area of sections c, d, e and f.
Figure 3. View of the tramway track in the area of sections c, d, e and f.
Sensors 20 04408 g003
Figure 4. Diagram of the layout of GNSS receivers on the bogies.
Figure 4. Diagram of the layout of GNSS receivers on the bogies.
Sensors 20 04408 g004
Figure 5. Diagram of the configuration of the measuring train.
Figure 5. Diagram of the configuration of the measuring train.
Sensors 20 04408 g005
Figure 6. Setting out the position of receivers in relation to the bogie’s structure and the tramway track in the hall of the Wrzeszcz tram depot.
Figure 6. Setting out the position of receivers in relation to the bogie’s structure and the tramway track in the hall of the Wrzeszcz tram depot.
Sensors 20 04408 g006
Figure 7. The measuring set used in the mobile satellite measurements duringmeasurements in Gdańsk.
Figure 7. The measuring set used in the mobile satellite measurements duringmeasurements in Gdańsk.
Sensors 20 04408 g007
Figure 8. A schema of the location of the track axis and the track centerline.
Figure 8. A schema of the location of the track axis and the track centerline.
Sensors 20 04408 g008
Figure 9. Measurement of the coordinates of the route in the PL-2000 system with the designation of individual straight sections.
Figure 9. Measurement of the coordinates of the route in the PL-2000 system with the designation of individual straight sections.
Sensors 20 04408 g009
Figure 10. Schematic diagram of the route discretization into sections; p i —measured points, d i —distance between points.
Figure 10. Schematic diagram of the route discretization into sections; p i —measured points, d i —distance between points.
Sensors 20 04408 g010
Figure 11. Schematic diagram of determining directions characterizing individual sections.
Figure 11. Schematic diagram of determining directions characterizing individual sections.
Sensors 20 04408 g011
Figure 12. Schematic diagram of the analysis of the variability of azimuths along the discretized route.
Figure 12. Schematic diagram of the analysis of the variability of azimuths along the discretized route.
Sensors 20 04408 g012
Figure 13. Schematic diagram of the reconstructed vertices system.
Figure 13. Schematic diagram of the reconstructed vertices system.
Sensors 20 04408 g013
Figure 14. Automatic division of the measuring route into straight sections.
Figure 14. Automatic division of the measuring route into straight sections.
Sensors 20 04408 g014
Figure 15. Result of the multiple reconstruction of the selected vertex (intersection) of the polygon.
Figure 15. Result of the multiple reconstruction of the selected vertex (intersection) of the polygon.
Sensors 20 04408 g015
Figure 16. Close-up of the system with reverse curves.
Figure 16. Close-up of the system with reverse curves.
Sensors 20 04408 g016
Figure 17. Comparison of the set of selected straights through a visual assessment (left) to a set of straights identified automatically (right).
Figure 17. Comparison of the set of selected straights through a visual assessment (left) to a set of straights identified automatically (right).
Sensors 20 04408 g017
Figure 18. The results of a series of six measuring sections Str 1Str 12 recorded by receiver no. 5.
Figure 18. The results of a series of six measuring sections Str 1Str 12 recorded by receiver no. 5.
Sensors 20 04408 g018
Figure 19. The results of a series of six measuring section of Str 9 characterized by poor track conditions, recorded by all receivers of the measurement system.
Figure 19. The results of a series of six measuring section of Str 9 characterized by poor track conditions, recorded by all receivers of the measurement system.
Sensors 20 04408 g019
Figure 20. Analysis of the accuracy of the reproduction of the main direction in section Str 12 featuring extremely unfavorable field screens.
Figure 20. Analysis of the accuracy of the reproduction of the main direction in section Str 12 featuring extremely unfavorable field screens.
Sensors 20 04408 g020
Figure 21. Analysis of the accuracy of the reproduction of the main direction in section Str 6 featuring very good visibility of satellites.
Figure 21. Analysis of the accuracy of the reproduction of the main direction in section Str 6 featuring very good visibility of satellites.
Sensors 20 04408 g021
Figure 22. The relative position of straight line Str 6 estimated from individual runs around the least squares straight line—antenna no. 5 (left) and no. 10 (right).
Figure 22. The relative position of straight line Str 6 estimated from individual runs around the least squares straight line—antenna no. 5 (left) and no. 10 (right).
Sensors 20 04408 g022
Table 1. Selected straights in the scope of the af sections of the geometric layout analyzed.
Table 1. Selected straights in the scope of the af sections of the geometric layout analyzed.
Straight LineLocation in the RouteAssessment of Track Condition/QualityAssessment of Visibility
Str 1Sect. aFairGood
Str 2Sect. aExcellentFair
Str 3Sect. aExcellentGood
Str 4Sect. aExcellentExcellent
Str 5Sect. aExcellentGood
Str 6Sect. bExcellentExcellent
Str 7Sect. cBadPoor
Str 8Sect. cBadPoor
Str 9Sect. cBadFair
Str 10Sect. cBadPoor
Str 11Sect. dGoodBad
Str 12Sect. eGoodPoor
Table 2. Input data for the analysis of the main directions for a selected straight.
Table 2. Input data for the analysis of the main directions for a selected straight.
Record No.GPS Time (s)Horizontal Coordinate Y (m)Vertical Coordinate X (m)Two Dimensional Position Error (m)Number of Satellites
3586702,029,347,009.006,540,504.81386,030,362.96990.04008
3586712,029,347,009.056,540,504.93456,030,362.84420.04008
3586722,029,347,009.106,540,505.05616,030,362.71320.04008
3586732,029,347,009.156,540,505.18876,030,362.57970.04008
Table 3. Comparison of the reconstructed parameters of straight sections.
Table 3. Comparison of the reconstructed parameters of straight sections.
Section NumberNumber of Measurements n Average   Length   L a v   [ m ] Average   Azimuth   A z a v   [ deg ] Standard   Deviation   σ A z   [ deg ] Indicator   σ A z ( r e l )   [ deg * m ]
11219.1196.1390.0260.500
21271.5186.5360.0020.116
311194.9201.4270.0010.176
412187.1202.1960.0020.332
5120.7135.7920.3460.233
612274.5136.4150.0020.670
71250.944.0510.0120.603
8737.341.7060.0080.313
Table 4. Comparison of the reconstructed parameters of the polygon’s vertices.
Table 4. Comparison of the reconstructed parameters of the polygon’s vertices.
Marking of Vertex Angle   of   Return   α a v   [ o ] Number   of   Measurements   n Vertical   Coordinate   X a v   [ m ] Standard   Deviation σ X   [ m ] Horizontal   Coordinate Y a v   [ m ] Standard   Deviation   σ Y   [ m ] Distance | W | a v   [ m ] Standard   Deviation   σ | W |   [ m ]
W1-29.60312−69.6560.011−62.4220.0640.0490.041
W2-314.89111−86.5350.007−209.7460.0250.0230.011
W3-40.76911−181.8610.152−452.6480.3910.3370.226
W4-566.40412−271.7960.098−673.0740.2470.1500.215
W5-60.62312−240.9441.284−704.8021.3391.2611.306
W6-792.36412−5.1700.017−952.5180.0200.0240.008
W7-82.345786.0520.126−858.2190.1250.1370.099

Share and Cite

MDPI and ACS Style

Wilk, A.; Specht, C.; Koc, W.; Karwowski, K.; Skibicki, J.; Szmagliński, J.; Chrostowski, P.; Dabrowski, P.; Specht, M.; Zienkiewicz, M.; et al. Evaluation of the Possibility of Identifying a Complex Polygonal Tram Track Layout Using Multiple Satellite Measurements. Sensors 2020, 20, 4408. https://doi.org/10.3390/s20164408

AMA Style

Wilk A, Specht C, Koc W, Karwowski K, Skibicki J, Szmagliński J, Chrostowski P, Dabrowski P, Specht M, Zienkiewicz M, et al. Evaluation of the Possibility of Identifying a Complex Polygonal Tram Track Layout Using Multiple Satellite Measurements. Sensors. 2020; 20(16):4408. https://doi.org/10.3390/s20164408

Chicago/Turabian Style

Wilk, Andrzej, Cezary Specht, Wladyslaw Koc, Krzysztof Karwowski, Jacek Skibicki, Jacek Szmagliński, Piotr Chrostowski, Pawel Dabrowski, Mariusz Specht, Marek Zienkiewicz, and et al. 2020. "Evaluation of the Possibility of Identifying a Complex Polygonal Tram Track Layout Using Multiple Satellite Measurements" Sensors 20, no. 16: 4408. https://doi.org/10.3390/s20164408

APA Style

Wilk, A., Specht, C., Koc, W., Karwowski, K., Skibicki, J., Szmagliński, J., Chrostowski, P., Dabrowski, P., Specht, M., Zienkiewicz, M., Judek, S., Skóra, M., & Grulkowski, S. (2020). Evaluation of the Possibility of Identifying a Complex Polygonal Tram Track Layout Using Multiple Satellite Measurements. Sensors, 20(16), 4408. https://doi.org/10.3390/s20164408

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop