Measuring Methods of Radius of Curvature and Tread Circle-Fitting Studies for Railway Wheel Proﬁles

: A railway wheel proﬁle consists of short arcs with complex radii of curvature, and wheel wear leads to changes in the proﬁle’s radius of curvature that ultimately affects the dynamic performance of the train. To track the evolution of in-service wheel proﬁle curves, the radii of curvature of new foundry wheel proﬁles need to be measured. This study proposes a series of algorithms and calculation methods for measuring the radius of curvature of wheel proﬁles. Firstly, the curvature was estimated with the U-chord method, and the segment points were located. Secondly, the discrete derivative method and Two-Arcs Tangency Constrain (TATC) method were used to calculate the radius of curvature and the ﬁtting circle radius, respectively. The experimental results of the three types of proﬁles showed that the wheel proﬁle curves were precisely divided according to the estimated curvature method and that the maximum errors of the calculated results compared with standard values by the discrete derivative method and TATC method were 2.50% and 0.42%, respectively. Furthermore, the two measurement methods’ performances and repeated experiments were used to analyze the uncertainty.


Introduction
Rail transit is an important mode of transportation in industrialized countries.Highfrequency and long-distance railway transportation has higher requirements for the safety and stability of trains.Due to the influence of external factors such as carrying capacity, rail terrain, rain erosion, and sudden braking, the situation of wheel wear inevitably occurs during long-term service [1].The wheel profile is an important feature of the wheelset; hence, the timely detection of changes in the wheel profile is the main task of wheelset detection and maintenance.Research in wheel wear has been conducted by some scholars [2,3]; all of them have compared worn wheel profiles to standard ones and neglected errors between the fresh-from-factory forged wheels and standard wheels.Although these errors meet tolerance zone requirements, the manufacturing errors should be analyzed in deeper studies about of the principles of wheel wearing or wheel-rail contact.To conduct studies that explore the actual Hertzian or rolling contact parameters, it is important to measure the true parameters of newly forged wheels rather than using design criterion values directly.
Workers mainly use wheel vernier calipers or laser instruments to measure some horizontal and vertical parameters (i.e., flange height, flange width, rim width, QR), but these parameters cannot manifest changes in the whole wheel curve entirely.The most obvious change for integrated measurements is the radius of curvature.The design of the wheel profile shape is conducted based on many countries' standards, and some additional factors must be considered, such as the degree of wheel-rail matching, vehicle stability, wearing model, and repair strategy.These reasons cause the wheel profile to have a complex curve shape with multiple arcs, as shown in Figure 1.These arcs' radii of curvature influence the train dynamics; some researchers have studied this phenomenon and carried out design optimization research on wheel profile shapes [4,5].Dynamic simulation experiments rely on the input of digital design curves, and the results of simulations present the wheel profile after wearing.When studying the mechanism of wheel wear, the Hertzian or rolling contact theory is often used; the parameters require the measurement of the radius of curvature at the contact point [6].However, the actual wheel profile curve is not exactly consistent with the digital curve, and the curvature of the actual wheel profile curve needs to be detected at all wearing stages, especially from the initial wheel shape.In this way, wear experiments of the actual wheel and simulations can be compared more rigorously.This study laid emphasis on measuring the radius of curvature of a new (before serviced) wheel.This study enriches the original characteristics of wheel profiles and provides a reference for the future analysis of worn wheels.
tional factors must be considered, such as the degree of wheel-rail matching, vehicle stability, wearing model, and repair strategy.These reasons cause the wheel profile to have a complex curve shape with multiple arcs, as shown in Figure 1.These arcs' radii of curvature influence the train dynamics; some researchers have studied this phenomenon and carried out design optimization research on wheel profile shapes [4,5].Dynamic simulation experiments rely on the input of digital design curves, and the results of simulations present the wheel profile after wearing.When studying the mechanism of wheel wear, the Hertzian or rolling contact theory is often used; the parameters require the measurement of the radius of curvature at the contact point [6].However, the actual wheel profile curve is not exactly consistent with the digital curve, and the curvature of the actual wheel profile curve needs to be detected at all wearing stages, especially from the initial wheel shape.In this way, wear experiments of the actual wheel and simulations can be compared more rigorously.This study laid emphasis on measuring the radius of curvature of a new (before serviced) wheel.This study enriches the original characteristics of wheel profiles and provides a reference for the future analysis of worn wheels.The curvature estimation of curves is composed of sequentially sampling discrete points (hereafter referred to as discrete curves) and is commonly used in the field of differential geometry or digital image processing.Mokhtarian et al. [7] proposed the Curvature Scale Space (CSS) algorithm in 1998; the algorithm calculates the curvature based on the contours extracted from an image using the parametric arc length formula for a curve and then measures the degree of inflection response based on the curvature value of each pixel point.The CSS algorithm led to a series of improved algorithms, such as the Curvature Scale Space Corner with Adaptive Threshold (ATCSS) [8], Multi-scale Curvature Product (MSCP) [9], Direct Curvature Scale Space (DCSS) [10], and so on.Zhang et al. [11] employed Chebyshev polynomial fitting to estimate continuous curvatures in corner detection issues.Zhang et al. [12] proposed that variations in the evolution difference in contour curves at different scales could be taken as the corner response function, and the Difference of Gaussian (DoG) detection operator was proposed.Although these methods were based on digital image data, they defined a new curvature expression method to achieve the purpose of feature recognition.However, the true value of curvature was not calculated using these methods.
In the context of how to fit an arc of approximately circular points to an approximate circle, two categories can be divided from recent research [13]: (1) Geometric fitting: the objective is to minimize the orthogonal distance from the fitting circle in order to approximate circular points.(2) Algebraic fitting: the objective is achieved with a series of algebraic formulas.Unfortunately, all of these studies mentioned nothing about when the approximate circular points become short approximate circular arc points (i.e., the approximate circular points are inadequate compared with the whole circle), and the fitting results do not have position constraint requirements.Hu et al. [14] proposed a parameterfree approach for weighted geometric circle fitting, and the corresponding bias was inves- The curvature estimation of curves is composed of sequentially sampling discrete points (hereafter referred to as discrete curves) and is commonly used in the field of differential geometry or digital image processing.Mokhtarian et al. [7] proposed the Curvature Scale Space (CSS) algorithm in 1998; the algorithm calculates the curvature based on the contours extracted from an image using the parametric arc length formula for a curve and then measures the degree of inflection response based on the curvature value of each pixel point.The CSS algorithm led to a series of improved algorithms, such as the Curvature Scale Space Corner with Adaptive Threshold (ATCSS) [8], Multiscale Curvature Product (MSCP) [9], Direct Curvature Scale Space (DCSS) [10], and so on.Zhang et al. [11] employed Chebyshev polynomial fitting to estimate continuous curvatures in corner detection issues.Zhang et al. [12] proposed that variations in the evolution difference in contour curves at different scales could be taken as the corner response function, and the Difference of Gaussian (DoG) detection operator was proposed.Although these methods were based on digital image data, they defined a new curvature expression method to achieve the purpose of feature recognition.However, the true value of curvature was not calculated using these methods.
In the context of how to fit an arc of approximately circular points to an approximate circle, two categories can be divided from recent research [13]: (1) Geometric fitting: the objective is to minimize the orthogonal distance from the fitting circle in order to approximate circular points.(2) Algebraic fitting: the objective is achieved with a series of algebraic formulas.Unfortunately, all of these studies mentioned nothing about when the approximate circular points become short approximate circular arc points (i.e., the approximate circular points are inadequate compared with the whole circle), and the fitting results do not have position constraint requirements.Hu et al. [14] proposed a parameter-free approach for weighted geometric circle fitting, and the corresponding bias was investigated with the classical theory of nonlinear least squares in this study.Liu et al. [15] proposed a radius constraint method to measure the structured light arc of wear of a rail profile.Zhu et al. [16] proposed a center constraint method to measure a rotational workpiece.However, these methods required obtaining prior information about the radius or center of the circle coordinates.In actual on-site measurements, accurate prior information is hard to obtain for railway inspectors.Tao et al. [17] proposed a parameterconstrained fitting method based on a four-parameter circle equation, and the experimental results showed that this method had a high accuracy when measuring short arcs with small central angels without prior information.
In order to compare the evolution of a wheel profile curve before and after wearing, a newly forged wheel profile's radii of curvature need to be measured.This research used a laser profile sensor to collect the wheel profile data.Aiming at the collected complex curve, a series of methods were proposed to measure the radius of curvature, which included demarcating segmented points (positioning tangency points), calculating the radius of curvature, and fitting a circle for small arcs.The remainder of this study consists of three parts.First, in Section 2, the theory of analysis and the calculation processing is proposed.The experimental device and profile data calculation are described in Section 3. In Section 4, the analysis of the uncertainty is summarized.

Composition of Wheel Profile Curve and Analysis Interests
During the function of train wheels, the wheel-rail contact point distribution is as shown in Figure 2 [18].The wear of the wheel surface caused by wheel-rail contact mainly occurs in the circular arc AB, the linear section BC, the circular arc CD, the circular arc DE, the circular arc EF, and the linear section FG.There exist sudden changes in the curvature of the wheel profile at the segmented points, which adversely affect the contact point distribution [19].Moreover, if the wheel profile curvature ρ w is close to the rail profile curvature ρ r , there exists a two-point contact or common contact, which aggravates wheel-rail wear [20,21].Therefore, short arcs of the wheel profile curves and segmented points were selected as this study's interests.
is hard to obtain for railway inspectors.Tao et al. [17] proposed a parameter-constrained fitting method based on a four-parameter circle equation, and the experimental results showed that this method had a high accuracy when measuring short arcs with small central angels without prior information.
In order to compare the evolution of a wheel profile curve before and after wearing, a newly forged wheel profile's radii of curvature need to be measured.This research used a laser profile sensor to collect the wheel profile data.Aiming at the collected complex curve, a series of methods were proposed to measure the radius of curvature, which included demarcating segmented points (positioning tangency points), calculating the radius of curvature, and fitting a circle for small arcs.The remainder of this study consists of three parts.First, in Section 2, the theory of analysis and the calculation processing is proposed.The experimental device and profile data calculation are described in Section 3. In Section 4, the analysis of the uncertainty is summarized.

Composition of Wheel Profile Curve and Analysis Interests
During the function of train wheels, the wheel-rail contact point distribution is as shown in Figure 2 [18].The wear of the wheel surface caused by wheel-rail contact mainly occurs in the circular arc AB, the linear section BC, the circular arc CD, the circular arc DE, the circular arc EF, and the linear section FG.There exist sudden changes in the curvature of the wheel profile at the segmented points, which adversely affect the contact point distribution [19].Moreover, if the wheel profile curvature ρw is close to the rail profile curvature ρr, there exists a two-point contact or common contact, which aggravates wheel-rail wear [20,21].Therefore, short arcs of the wheel profile curves and segmented points were selected as this study's interests.The following Figure 3 shows the methodology of the measuring procedure used to process the wheel profile curve.The following Figure 3 shows the methodology of the measuring procedure used to process the wheel profile curve.

Methods of Measuring the Radius of Curvature from the Discrete Curve of a Wheel Profile
The complex curve shape of a wheel profile can be measured by a laser profile sensor, as Figure 4 shows.The laser profile is based on the laser triangle principle [22], and the measured object's profile is collected as a line point cloud, which is composed of a series of sample points with coordinates information.Thus, these sample points comprise a discrete curve of discrete point coordinates (x n , y n ) [23][24][25].

Methods of Measuring the Radius of Curvature from the Discrete Curve of a Wheel Profile
The complex curve shape of a wheel profile can be measured by a laser profile sensor, as Figure 4 shows.The laser profile is based on the laser triangle principle [22], and the measured object's profile is collected as a line point cloud, which is composed of a series of sample points with coordinates information.Thus, these sample points comprise a discrete curve of discrete point coordinates (xn, yn) [23][24][25].

U-Chord Curvature Estimation Method
For the task of wheel profile design, a curve line or a vector function expression of a curve with clear demarcation points can be constructed through computer-aided geometric design (CAGD) software.However, the wheel profile curves collected from on-site measurements contain noisy data points, and the demarcation point of the curve segment is not obvious in the visual field; therefore, it is first necessary to identify the segmented points (tangency points).Segmented points can be discovered by curvature estimation.This research used the U-chord curvature method to estimate the segmented point positions of each short arc and used a refined curve strategy according to the difference between the curve of the discrete points and the digital curve.
The method begins with the following:

Methods of Measuring the Radius of Curvature from the Discrete Curve of a Wheel Profile
The complex curve shape of a wheel profile can be measured by a laser profile sensor, as Figure 4 shows.The laser profile is based on the laser triangle principle [22], and the measured object's profile is collected as a line point cloud, which is composed of a series of sample points with coordinates information.Thus, these sample points comprise a discrete curve of discrete point coordinates (xn, yn) [23][24][25].

U-Chord Curvature Estimation Method
For the task of wheel profile design, a curve line or a vector function expression of a curve with clear demarcation points can be constructed through computer-aided geometric design (CAGD) software.However, the wheel profile curves collected from on-site measurements contain noisy data points, and the demarcation point of the curve segment is not obvious in the visual field; therefore, it is first necessary to identify the segmented points (tangency points).Segmented points can be discovered by curvature estimation.This research used the U-chord curvature method to estimate the segmented point positions of each short arc and used a refined curve strategy according to the difference between the curve of the discrete points and the digital curve.
The method begins with the following: From Equations ( 1) and ( 2), l is defined as a discrete point curve, and Ω is the constraint-supported domain, which is determined by selected parameter U.The parameter

U-Chord Curvature Estimation Method
For the task of wheel profile design, a curve line or a vector function expression of a curve with clear demarcation points can be constructed through computer-aided geometric design (CAGD) software.However, the wheel profile curves collected from on-site measurements contain noisy data points, and the demarcation point of the curve segment is not obvious in the visual field; therefore, it is first necessary to identify the segmented points (tangency points).Segmented points can be discovered by curvature estimation.This research used the U-chord curvature method to estimate the segmented point positions of each short arc and used a refined curve strategy according to the difference between the curve of the discrete points and the digital curve.
The method begins with the following: From Equations ( 1) and ( 2), l is defined as a discrete point curve, and Ω is the constraintsupported domain, which is determined by selected parameter U.The parameter U is subjected to the constraint condition, as shown in Equation ( 3).Since no point of the discrete curve absolutely satisfies Equation (3), the points P i−b and P i+ f that satisfy the support domain are searched by the following algorithm (Algorithm 1), and the theory is shown in Figure 5.
Algorithm 1. U-chord curvature method of discrete curve.

Input:
U: the constraint-supported domain.P i (x i , y i ): points' coordinates of the discrete curve.i max : the maximum number of points in one discrete curve.

Output:
U i : the estimated curvature value of each point in the profile curve.

Output:
Ui: the estimated curvature value of each point in the profile curve.
Machines 2023, 11, 181 6 of 26 support domain are searched by the following alg shown in Figure 5.According to the calculated results of the U i distribution of each point, the segmented point of the wheel profile curve is determined.The coordinates of the segmented points of the profile wheel are determined according to the locations of the U i curve saltation points.

Discrete Derivative Method Analysis for Wheel Profile Discrete Curve
Generally speaking, the connection between the derivative and the radius of curvature results in ways to calculate the radius of curvature.However, the wheel profile curve is discrete, and the function is not written directly.The solution of the discrete derivatives method [26] is based on differential geometry and tangent line estimation.The tangent line l of the point P k (x k , y k ) is set, whose discrete point function is y i = f (x i ) (i∈{1, 2, . . ., n}), and the tangent line l has two attributes: (1) the tangent line passes through a fixed point P k ; and (2) the tangent line is a local line approximation of the fixed point's symmetric neighborhood domain.In other words, the sum of the perpendicular distances from the points from the P k symmetric neighborhood domain to the tangent line is a minimum.According to the above attributes, the relationship between tangent line l and y i = f (x i ) based on the least squares method can be described as follows: The Lagrange multiplier is used to construct the Lagrange function: When setting ∂L ∂a = 0, ∂L ∂b = 0, ∂L ∂λ = 0, the result is as follows: The differential result value f (x k ) of point P k is replaced by the curvature value a, and m 1 is the value of half of the neighborhood domain's distance calculated from the point P k .Furthermore, the differential coefficient result value f (x k ) can be calculated  9), the second derivative value of point P k is calculated with Equation (10).(10) Due to the correlation between the numerical value of the radius of curvature and the derivative value, the radius of the curvature value is calculated with Equation (11).
In spatial analytic geometry, space curves are represented by parameterization, i.e., the coordinates of a point p on the space curve are written as a scalar function of a parameter u, as shown in Equation (12).In differential geometry, these forms are written as row vectors; thus, the normal curve is represented as a vector function about parameter u, as shown in Equation ( 13).
x = x(u), y = y(u), z = z(u) In Equation ( 13), due to the wheel profile in two-dimensional space, only the x and y coordinates are parameterized.In contrast to the normal continuous differentiable parametrized curve, the wheel profile curve is discrete.The purpose of parameterization is to assign corresponding parameters to the curve points and, in turn, build an approximate function expression to approach the discrete curve.The common parameterization methods include uniform, cumulative chord length, endocentric parameters, Foley, and so on.In this research, the cumulative chord length method was selected to assign parameter u, and the u i of each point P i (x i , y i ) was calculated as shown in Equation (14).
The wheel profile curve can be written as follows: According to the definition of partial derivatives and discrete differentials, the specific first derivative is shown in Equation ( 16), and the second derivative has a similar form.
Naturally, the radius of curvature value of Equation ( 11) can be re-written as follows:

The Circle-Fitting Method of TATC
The sectors of a connection tread circle and central tread circle captured from a wheel profile have flat and large curvature values.Moreover, these two arcs comprise a large proportion of the whole profile curve; thus, the quantity of noise points (outliers) caused by uneven light reflection or overexposure is greater than those for other arcs.The discrete derivative method is defined in the local neighborhood approximation domain, and if a set of noise points are mixed in the neighborhood domain, the calculation result is widely affected.Therefore, the fitting circle method is applied in the circle of the connection tread and the circle of the central tread, as well as according to the result of the circle fitting in order to estimate the radius of the osculating fitting circle.
A wheel profile consists of continuous short arcs, that is, there exists a tangency relation between every two arcs.The circle of connection tread and the circle of the central tread are no exception.In this study, we proposed a two-arcs tangency constraint (hereinafter referred as to TATC) method to fit a circle via a new objective function.The basic idea of TATC is to fix a circle and establish another circle that fits the objective function.Then, the LM (Levenberg-Marquardt) algorithm and IRLS (iteratively reweighted least squares) method were used to solve the objective function.Since the radius value of the former circle is included in the latter circle's parameters, the tangency relationship can be guaranteed.This study selected the circular arc of DE as the fixed circle O 2 , using the TATC method in detail as follows.
The set of points from the circular arc of EF were defined as establishing the fitting circle objective function of EF: In Equation ( 18), (a 1 , b 1 ) is the center coordinate of fitting circle O 1 , n is the number of points from the arc of EF, and R 1 is the radius of fitting circle O 1 .On account of the tangency relation there exists the following relationship between circles O 1 and O 2 : In Equation ( 19), (a 2 , b 2 ) is the center coordinate of fitting circle O 2 , and R 2 is the radius of fitting circle O 2 , which requires R 1 > R 2 .a 2 , b 2 , and R 2 could be obtained from previous studies [27].Equation ( 19) is substituted into Equation (18), and the objective function is changed, as the following shows: In Equations ( 20) and ( 21), ε i is the error of fitting precision, ρ(•) represents the loss function, and Cauchy's function [28] was chosen to construct the loss function in this study, as shown in Equation (22).Equation ( 23) is the derivation result of Equation ( 22), and Equation ( 23) was used to calculate the Jacobian matrix of the objective function.Equation ( 24) is the ratio of Equation ( 23)/Equation (21) and is a defined weight function.
Machines 2023, 11, 181 9 of 26 In the equations, σ is a variable of scale parameter, and the variable is updated by Equations ( 26)- (28).ε i,k−1 represents the fitting error value, which is calculated by inputting k−1th times the fitting circle results into the calculation at point P i .
In Equation ( 26), M (•) is the median function, and ε i,0 is the fitting error value of the first iterative calculation.σ 0 is the scale parameter of the first iterative calculation, where the value of σ 0 is determined by Equation ( 26), and the latter value is determined by Equation (27).In Equation ( 27), ξ is a convergent parameter, and e is the natural logarithm.n s is the number of iterations.τ is an accelerated parameter of convergence, as shown in Equation (28).In Equation (28), δ is the noise intensity value of measurement, and it is determined by the sensor's precision.Through Equations ( 26)-( 28), the value of σ in Equations ( 22)-( 24) is decided.
The nonlinear objective function, as shown in Equation (25), is difficult to solve with the least squares method.Common solution methods include the Gauss-Newton iteration method and the gradient descent method.The Gauss-Newton method is sensitive to the initial values, which is in contrast to the gradient descent method.Therefore, in order to preserve solution procedure accuracy, the LM (Levenberg-Marquardt) algorithm was chosen in this research, and is written concretely as follows: In Equation (29), µ is the damping coefficient, h is the direction of iteration, J is the Jacobian matrix of the objective function, I is the identity matrix, and f is the objective function.The trust region method was applied in the LM algorithm, which can update the damping coefficient µ.Next, in order to reduce the influence of noise points when calculating a fitting circle, the idea of IRLS was introduced: when the kth calculation is completed, the calculation results (i.e., the coordinate of the center of circle ) of k times are input into Equation ( 23); each point's error (i.e., ε i,k−1 ) in arc EF is calculated.Then, Equations ( 24)-( 29) are utilized to update the objective function Equation (25).Then, the objective function is solved using Equation ( 29 [27].w 0 : initial weight matrix.q max : the maximum number of iteration times.∆ q : the minimum precision.

Experimental Equipment and Test Objectives
A laser profile sensor was used to scan the wheel profile, and the sensor equipment and measurement procedure are shown in Figure 6.The detailed key performance parameters of the sensor are shown in Table 1.The three types of scanned wheel profiles included LMA, LMC, and JM, all of which are common types of wheels in service train bogies, and the standard values of each profile are shown in Table 2 [29].

Experimental Equipment and Test Objectives
A laser profile sensor was used to scan the wheel profile, and the sensor equipment and measurement procedure are shown in Figure 6.The detailed key performance parameters of the sensor are shown in Table 1.The three types of scanned wheel profiles included LM A , LM C , and JM, all of which are common types of wheels in service train bogies, and the standard values of each profile are shown in Table 2 [29].
eters of the sensor are shown in Table 1.The three type cluded LMA, LMC, and JM, all of which are common types and the standard values of each profile are shown in Tab

Segmentation Point Experiments Based on U-Chord Curvature Estimation
The wheel profile collected by the sensor is shown in Figure 7.The LM A -type profile had 1588 points, and the LM C -and JM-type profiles had 1540 and 1539 points, respectively.Thus, the three types of wheel profiles had similar shapes but different radius curvature arcs, as Table 2 shows.

Segmentation Point Experiments Based on U-Chord Curvature Estimation
The wheel profile collected by the sensor is shown in Figure 7.The LMA-type profile had 1588 points, and the LMCand JM-type profiles had 1540 and 1539 points, respectively.Thus, the three types of wheel profiles had similar shapes but different radius curvature arcs, as Table 2 shows.When inputting the discrete curve of the wheel profile, Algorithm 1 was chosen to estimate the U-chord curvature method.Then, the segmented points' coordinates of each arc were determined.The U-chord curvature method's result was affected by input parameter U; thus, inputting different parameters and comparing the results was used to find the best distinct curvature distribution.Taking the LMA-type profile as an example, the results are shown in Figure 8.When inputting the discrete curve of the wheel profile, Algorithm 1 was chosen to estimate the U-chord curvature method.Then, the segmented points' coordinates of each arc were determined.The U-chord curvature method's result was affected by input parameter U; thus, inputting different parameters and comparing the results was used to find the best distinct curvature distribution.Taking the LM A -type profile as an example, the results are shown in Figure 8. Figure 8a-f represent the U-chord curvature results from selecting different values of U (i.e., U = 1.2, …, 2.2.); in each parameter case, the curvature of each point was estimated, and the wheel profile curve was drawn on the same sheet, which facilitated the observation of the curvature distribution of each point on the profile curve.In fact, the U-chord curvature method is also a class of curvature scale space, and different U values present different parameters of scale space.It is necessary to find the most distinctive features in a scale space.The demarcation points of each segment arcs can be determined by the local peak point or bottom point.In Figure 8a-c, the wave of the estimated curvature was too Figure 8a-f represent the U-chord curvature results from selecting different values of U (i.e., U = 1.2, . . ., 2.2.); in each parameter case, the curvature of each point was estimated, and the wheel profile curve was drawn on the same sheet, which facilitated the observation of the curvature distribution of each point on the profile curve.In fact, the U-chord curvature method is also a class of curvature scale space, and different U values present different parameters of scale space.It is necessary to find the most distinctive features in a scale space.The demarcation points of each segment arcs can be determined by the local peak point or bottom point.In Figure 8a-c, the wave of the estimated curvature was too disordered.In Figure 8e, point E and point F were not clearly recognized.In Figure 8d, when the value of U was equal to 1.8, the segmented points A~F were distinguishable, and the segmented points were marked by recognition.

R REVIEW 12 of 24
In order to verify the accuracy of the segmented point coordinates, a collected curve of a wheel profile was registered to the wheel profile of the Chinese national standard [29,30], as shown in Figure 9.According to the corresponding points of the two different source curves, the selected segmented points' indexes and coordinates from the scanned curve were compared to position accuracy.Figure 9 shows the coordinate transformation results of the scanned LM A -type profile, which were registered to standard values plotted by the CAD software.The points' indexes were recorded before and after registration, as shown in Table 3.The maximum error was 0.790 mm; thus, the selected points could be thought of as arc demarcation points.According to the results of Section 2.2, we adopted the discrete derivative method to calculate the radius of curvature of the wheel profile.The discrete derivative needs to determine the value of the neighborhood domain m 1 , m 2 .Hence, a pre-experiment was performed before the formal experiment, whose purpose was to select a suitable neighborhood domain to adapt different curvatures.
The wheel profile curve was divided into two classifications due to the morphology of the curve, as shown in Figure 10.On the basis of the wheel profile types in service from China's national standard, the value of the radius of curvature of the convex arcs included 15, 16, 18 mm, and so on; the radius of the concave arcs included 13 and 14 mm.Thus, a radius calibration block was used in the pre-experimentation, and a suitable neighborhood domain value was determined.The pre-experiment process is shown in Figure 11.

The Radius of Curvature Calculation Based on Discrete Deri
According to the results of Section 2.2, we adopted the calculate the radius of curvature of the wheel profile.The determine the value of the neighborhood domain m1, m2.performed before the formal experiment, whose purpose w borhood domain to adapt different curvatures.
The wheel profile curve was divided into two classific of the curve, as shown in Figure 10.On the basis of the whee China's national standard, the value of the radius of curvatu 15, 16, 18mm, and so on; the radius of the concave arcs inc radius calibration block was used in the pre-experimentat hood domain value was determined.The pre-experiment p    According to the results of the pre-experiment, and combined with the segmented points, the results of the discrete derivative experiment are shown in Figure 12; each segment arc's radii of curvature values are shown in Tables 4-6.The Huang method [31] and the Yang method [32] were selected for comparison.In the figures and tables, the discrete derivatives are written as DD; the different methods' relative errors are plotted in Figure 13.ment arc's radii of curvature values are shown in Tables 4-6.The Huang method [31] and the Yang method [32] were selected for comparison.In the figures and tables, the discrete derivatives are written as DD; the different methods' relative errors are plotted in Figure 13.As shown in Figure 12, to the left side of point D, or in other words, in the small radius of curvature segments, the discrete derivative method had high-precision results.

The Fitting Circle Experiment of the Circle of the Connection Tread and the Circle of Central Tread
Using the Pratt fitting circle method to fit the data points from arc DE, the result of the radius was 90.3969 mm, and the relative error of the standard val 0.44%.According to Algorithm 2, the fitting circle result of arc EF was calcula shown in Figure 14.The result of the Pratt method used to fit the points of arc 481.869 mm; this gave a large discrepancy compared with the standard value.Th because arc EF had noise data.The arc was short and had a small central angle, m As shown in Figure 12, to the left side of point D, or in other words, in the small radius of curvature segments, the discrete derivative method had high-precision results.According to Table 4, the relative error was less than 3.05%.The Huang and Yang methods are based on calculating the tangency circle in the radius curvature definition; however, these methods had large errors in their high-density line point cloud data.In the arcs of DE and EF, with the increase in the radius of curvature, the noise points caused by the laser source had a great effect on the calculation results, and the discrete derivative method was no longer applicable.The same situation was observed in the measurements of the other two profiles.
As shown in Tables 5 and 6, similar situations also occurred in the LM C -and JM-type wheel profiles.All of the results in the small-radius curvatures had a high accuracy, and the maximum relative errors in the two types were less than 3.23% and 3.58%, respectively.Furthermore, both the DE and EF arcs' measurement errors severely deviated.The test environment and the associated definition of the discrete derivative, as well as the existence of noise data points, contributed to the approximate calculation errors.Moreover, the noise points affected the symmetric neighborhood domain, and a tiny change in the neighborhood domain influenced the calculation results.

The Fitting Circle Experiment of the Circle of the Connection Tread and the Circle of the Central Tread
Using the Pratt fitting circle method to fit the data points from arc DE, the fitting result of the radius was 90.3969 mm, and the relative error of the standard value was 0.44%.According to Algorithm 2, the fitting circle result of arc EF was calculated, as shown in Figure 14.The result of the Pratt method used to fit the points of arc EF was 481.869 mm; this gave a large discrepancy compared with the standard value.This was because arc EF had noise data.The arc was short and had a small central angle, making the circle shape  The TATC method considered the link of the parameters of arc DE when fitting the circle to arc EF, and the result was 451.820 mm, and the relative error was 0.40%.This was because the TATC method's objective function was based on the tangency relation of arc DE, and the radius was constrained by arc DE's parameters.In addition, the distance of the two circles' tangency point to the segmented point E was 0.719 mm; this result showed   The TATC method considered the link of the parameters of arc circle to arc EF, and the result was 451.820 mm, and the relative error The TATC method considered the link of the parameters of arc DE when fitting the circle to arc EF, and the result was 451.820 mm, and the relative error was 0.40%.This was because the TATC method's objective function was based on the tangency relation of arc DE, and the radius was constrained by arc DE's parameters.In addition, the distance of the two circles' tangency point to the segmented point E was 0.719 mm; this result showed that the tangency relation was met from another point of view.
In addition, the wheel profile of the JM type was also tested.The general result of fitting the circle of arc EF is shown in Figure 16a.Similarly to the calculating process of the LM A -type curve, the Pratt fitting result of arc DE's radius was 100.512 mm, and the arc EF was 470.451 mm.The local results of the TATC method are shown in Figure 16bd.The radius of curvature value calculated by the TATC method was 492.461 mm, and the relative error was 1.51%.It was found that the cambered curve position of fitting the circle by the TATC method was closer to the standard circle than that of the Pratt method, as shown in Figure 16b,c, and the value showed a high precision.Figure 16d shows the noise points of the scanned curve segments, and the Pratt method was effectively influenced by these undulating points.The IRLS method, which was used in the TATC calculating process, played a role in facing noise points.Despite the Pratt method fully taking into consideration the minimization of the fitting errors, the noise points' influence affected the measurement accuracy.Thus, the procedure of reweighting with Equation ( 24

The Measurement Analysis of Uncertainty Evaluation
As a result of the existence of measuring errors and calculation errors, the true value was difficult to obtain.Generally, the uncertainty of measurement is used to reflect the

The Measurement Analysis of Uncertainty Evaluation
As a result of the existence of measuring errors and calculation errors, the true value was difficult to obtain.Generally, the uncertainty of measurement is used to reflect the quality and credibility of the determined values.Through a discrete derivative method and the TATC method, the LM A type was used as an example, and arc AB and arc EF were selected as the objects of the study.
In the following chapters, we kept four decimal places during the display of the calculated processing (the data figures collected by the laser profile sensor were in units of microns).In the final results, we consulted the relevant regulations of significant figures and the round-off method, and we kept two significant figures in the uncertainty values to ensure that the decimal point had the same placeholder of digits.7 shows three sets of experiments; every set was scanned at different times (different environmental sunlight), different temperatures, and different operators, but all of them were in the same ten positions of the train wheel profile.Bayes information fusion, which is based on both prior information and present sample information, was selected to analyze the results of the wheel profile radius of curvature [33].Figure 17 shows the calculation procedure for arc AB.Hence, the best estimated value of the radius of curvature from the three sets was 15.0984 mm, and the type-A standard uncertainty u a was 0.0211 mm.The degrees of freedom in the type-A evaluation of measurement uncertainty was equal to the degrees of freedom of the standard deviation, i.e., 3 * (10 − 1) = 27.

Type-B Uncertainty of arc AB
According to Table 1, the precision parameters of the laser sensor included linearity and repeatability.In the case of arc AB, the accuracy of the uncertainty results of the radius of curvature largely depended on the ordinate axes of the area of the arc, i.e., the sag across the arc.For these purposes, the arch height and chord length of arc AB were considered in the uncertainty analysis, and the uncertainty model was as follows: As shown in Figure 18, the relationship of the radius and the straight-line distance parameters (arch height h and chord length S) is shown in Equation (32).
The sensitivity coefficients were calculated as follows:

Type-B Uncertainty of arc AB
According to Table 1, the precision parameters of the laser sensor included linearity and repeatability.In the case of arc AB, the accuracy of the uncertainty results of the radius of curvature largely depended on the ordinate axes of the area of the arc, i.e., the sag across the arc.For these purposes, the arch height and chord length of arc AB were considered in the uncertainty analysis, and the uncertainty model was as follows: As shown in Figure 18, the relationship of the radius and the straight-line distance parameters (arch height h and chord length S) is shown in Equation (32).
FOR PEER REVIEW The actual measurement results show and 3.0820 mm, respectively.Table 8 show arc AB.The degrees of freedom of the typecalculated with the following equation:  The sensitivity coefficients were calculated as follows: The actual measurement results showed that the values of S and h were 15.3247 mm and 3.0820 mm, respectively.Table 8 shows the budget of calculating the uncertainty of arc AB.The degrees of freedom of the type-B evaluation of measurement uncertainty were calculated with the following equation: The combine standard uncertainty The actual calculation results showed that ∆u(x i )/u(x i ) ≈ 0 and that the value of v i was too large.This meant that the probability of the measured value falling outside the interval was low.Finally, the coefficient value k was 1.96 when setting the confidence interval as 95%.Consequently, the expand uncertainty U = k * u c ≈ 0.1803 mm.

Type-A Uncertainty Analysis of arc EF
Similarly, aiming to assess arc EF by means of using Equations ( 30) and (31), the best estimated value of the radius of curvature from the three sets was 450.8590 mm, and the type-A standard uncertainty u a was 1.0945 mm.In the same way, the degrees of freedom in the type-A evaluation of measurement uncertainty was equal to the degrees of freedom of standard deviation, i.e., 3 * (10 − 1) = 27.

Type-B Uncertainty Analysis of arc EF
There was a difference with arc AB, because the arc of EF was measured by the fitting circle method; the computational solution function can be expressed approximately as a least squares model: In Equation (36), (x k , y k ) is the coordinate point of fitting points, and M is the number of points.R is the value of the fitting circle's radius.According to Equation (37), the value of the radius is affected by every point in the arc curve, and the errors of those points are measured independently.Thus, the type-B standard uncertainty u B was calculated as follows: In Equation (37), u x i and u y i are the errors of the x-axis and z-axis, respectively.C x i and C y i are sensitivity coefficients, which can be calculated by taking Equation (37) to the partial derivatives of x k and y k .
The detailed calculation parameters are shown in Table 9. Referring to Equation (35), the degrees of freedom of the type-B standard uncertainty was also too large.This was subjective for the estimation of the value u b to a certain extent.Finally, the combined standard uncertainty u c was calculated as u c = (u a ) 2 + (u b ) 2 = 1.1200 mm.The expand uncertainty U = k * u c ≈ 2.1952 mm.From the abovementioned results, and according to the conservative rule of rounding off, with a 95% confidence interval, the final result of the radius curvature of arc AB was 15.10 ± 0.19 mm, the radius curvature of arc EF was 450.9 ± 2.2 mm, for which two significant figures were kept in the uncertainty approximation, and the point was kept to the hundredths place.Compared with standard values 15 mm and 450 mm, these errors kept 0.60~1.93%and 0.07~0.28% in each arc, respectively.

Conclusions
Aiming at complex multi-segmented short arcs in a new wheel profile, this study completed the radius of curvature and fitting circle measurement in a wheel profile curve using a laser profile sensor.The measurement methods of three types were verified by comparing them to standard design values.This research made the following major observations: (1) The U-chord curvature method fairly met the demand of the precision of segmented points, and the segmented points' position errors were less than 0.790 mm in the LM A type. (2) By the pre-experiment and the discrete derivative calculations in the small radius of curvature segments, the radii of curvature of the three wheel profile types were measured.In arc AB, the relative errors were 1.83%, 1.85%, and 0.60% for the LM A , LM c , and JM types, respectively.In arc CD, the relative errors were 2.19%, 2.5%, and 2.22%, respectively.(3) In the large radius of curvature segments, the discrete derivative calculation results occurred with large errors; this was because the discrete derivative method, which was based on a limiting approach method, did not work, and the noise points affected this method's accuracy.For the continuous two short arcs EF of the wheel profile, the proposed TATC method and solution method accurately calculated the radius of curvature of arc EF, and the maximum relative errors were 0.40% and 1.51% for the LM A and JM types, respectively.(4) With the LM A -type profile as an example, the uncertainty evaluation results of arc AB and arc EF showed that the radius curvature confidence intervals were ±0.19 mm and ±2.1 mm, respectively, with a 95% confidence probability.However, the algorithm of measuring radius of curvature of the wheel profile was restricted by using newly made wheels, and these curves had clearly comparable parameters.However, in the situation of worn wheels, the worn wheel profile curves are random, and the algorithm procedure was not well-suited and needs to be improved.Moreover, a worn wheel profile's radius of curvature is complex.Further research requires more field measurement data to attain a more comprehensive understanding of the relationship between wheel wear and profile curvature.

Figure 1 .
Figure 1.The complex curve of a wheel profile (unit: mm).

Figure 1 .
Figure 1.The complex curve of a wheel profile (unit: mm).

Figure 2 .
Figure 2. The measured curve in this study.(a) The distribution of wheel-rail contact points.(b) Each circular arc of the standard wheel profile.

Figure 2 .
Figure 2. The measured curve in this study.(a) The distribution of wheel-rail contact points.(b) Each circular arc of the standard wheel profile.

Figure 3 .
Figure 3.The procedure for measuring wheel profile radius of curvature.

Figure 4 .
Figure 4. Using a laser profile to collect the wheel profile curve.(a) Using a laser profile to measure the wheel profile.(b) Collecting the wheel profile as a discrete curve.

Figure 3 .
Figure 3.The procedure for measuring wheel profile radius of curvature.

Figure 3 .
Figure 3.The procedure for measuring wheel profile radius of curvature.

Figure 4 .
Figure 4. Using a laser profile to collect the wheel profile curve.(a) Using a laser profile to measure the wheel profile.(b) Collecting the wheel profile as a discrete curve.

Figure 4 .
Figure 4. Using a laser profile to collect the wheel profile curve.(a) Using a laser profile to measure the wheel profile.(b) Collecting the wheel profile as a discrete curve.

Figure 5 .
Figure 5.The supported domain of U-chord curvature.

Figure 5 .Figure 5 .
Figure 5.The supported domain of U-chord curvature

Figure 6 .
Figure 6.Using laser profile sensor to scan the wheel profile.

Table 1 .
The key performance parameters of the laser sensor.Parameters ValueMaximum sample points in one line point cloud 1920 Repeatability of x-axis 40 µm Repeatability of z-axis 20 µm Linearity ±0.01 %FS Temperature drift ±0.01 %FS• • C −1

Figure 7 .
Figure 7. Different types of wheel profiles.

Figure 7 .
Figure 7. Different types of wheel profiles.

Figure 8 .
Figure 8. Curvature estimation based on the U-chord method.

Figure 8 .
Figure 8. Curvature estimation based on the U-chord method.

Figure 9 .
Figure 9. Registration of measured curve to standard profile curve.

Figure 10 .
Figure 10.The classification of wheel profile curve based on morp

Figure 10 .
Figure 10.The classification of wheel profile curve based on morphology.

Figure 11 .
Figure 11.The pre-experiment of selecting su

Figure 11 .
Figure 11.The pre-experiment of selecting suitable neighborhood domain values.

Figure 12 .
Figure 12.Absolute errors of radii of curvature for all points in the LMA-type wheel profile.

Figure 12 .
Figure 12.Absolute errors of radii of curvature for all points in the LM A -type wheel profile.

Figure 13 .
Figure 13.The relative errors of measuring different types of wheel profiles' radii of curva LMA-type profile.(b) LMC-type profile.(c) JM-type profile.

Figure 13 .
Figure 13.The relative errors of measuring different types of wheel profiles' radii of curvature.(a) LM A -type profile.(b) LM C -type profile.(c) JM-type profile.

Figure 14 .
Figure 14.The result of fitting circle to arc EF of LMA type.

Figure 15 .
Figure 15.The contrast between measured and standard arc EF.

Figure 14 .
Figure 14.The result of fitting circle to arc EF of LM A type.

Figure 14 .
Figure 14.The result of fitting circle to arc EF of LMA type.

Figure 15 .
Figure 15.The contrast between measured and standard arc EF.

Figure 15 .
Figure 15.The contrast between measured and standard arc EF.

Figure 16 .
Figure 16.The result of fitting circle to arc EF of JM type.

Figure 16 .
Figure 16.The result of fitting circle to arc EF of JM type.

4. 1 .
The Uncertainty Analysis in arc AB 4.1.1.Type-A Standard Uncertainty of arc AB Table

Figure 17 .
Figure 17.Type-A standard uncertainty calculation procedure based on Bayes information fusion.

Figure 17 .
Figure 17.Type-A standard uncertainty calculation procedure based on Bayes information fusion.

Figure 18 .
Figure 18.Estimation of the radius of curvature

Figure 18 .
Figure 18.Estimation of the radius of curvature by arch height and chord length.

Table 2 .
The segmented parameters of different wheel profile types.

Table 3 .
Comparison of the estimated and the standard segmented points' coordinates according to registration.

Table 4 .
The results of calculating LM A -type profile's radius of curvature.

Table 5 .
The results of calculating LM c -type profile's radius of curvature.

Table 6 .
The results of calculating JM-type profile's radius of curvature.

Table 6 .
The results of calculating JM-type profile's radius of curvature.

Table 7 .
Three sets of wheel profile curve measuring results.

Table 8 .
The procedure budget about the uncer

Table 8 .
The procedure budget about the uncertainty of arc AB.

Table 9 .
The procedure budget about uncertainty of arc AB.