轨道曲线参数检测算法优化.pdf

第 63 卷 第 10 期2023 年10 月铁道建筑Railway EngineeringVol.63 No.10October 2023文章编号：10031995（2023）10000105轨道曲线参数检测算法优化李颖1 程朝阳1 陆永彪2 陈仕明1 余宁1 秦哲1 那强21.中国铁道科学研究院集团有限公司 基础设施检测研究所，北京 100081；2.北京铁科英迈技术有限公司，北京 100081摘要 现有的轨道检测系统对曲线相关参数检测准确性不高，尤其是5 000 m以上大半径曲线，导致曲线数据报表利用率较低。针对这一问题，提出在使用经典DP（DouglasPeucker）算法提取粗略曲线线型变化特征点的基础上，在该特征点周围窗口内使用角度阈值法提取角点，将角点连线成走向段，通过顺次寻找极值走向段的方法，确定曲线线型变化的直缓点、缓圆点、圆缓点、缓直点四大特征点，从而准确计算曲线平均半径、长度等参数，提高系统曲线检测准确性。将本文算法应用到车载轨道检测系统曲线参数检测中，并与现有算法及只采用DP算法提取的特征点进行对比。结果表明，使用本文算法对大半径和小半径曲线参数检测准确率均能达到95%以上，本文算法具有良好的适应性。关键词 轨道几何；曲线参数；特征点提取；DP算法；角点检测；极值段中图分类号 U216.3 文献标识码 A DOI：10.3969/j.issn.10031995.2023.10.01引用格式：李颖，程朝阳，陆永彪，等.轨道曲线参数检测算法优化 J.铁道建筑，2023，63（10）：15.曲线地段轨道在列车离心力的作用下，比直线地段受力大，轨道几何尺寸不易保持，是铁路线路的薄弱环节。轨道曲线参数的平顺性关系到列车运行的平稳性与稳定性，影响列车通过曲线地段容许的最高运行速度。为使曲线保持圆顺，维护良好状态，工务部门根据轨检车提供的周期性检测数据对线路进行调整，因此对轨道曲线参数准确检测提出了要求1。目前国内广泛使用的是GJ6型轨道检测系统以及最新研制的数字式轨道检测系统，能够检测的曲线半径为R=150 12 000 m，但是对R=5 000 8 000 m的曲线线路检测准确性不高，对于提速线路和高速线路上R 8 000 m的大半径曲线线路的识别和检测准确性较低，从而导致系统提供的轨道不平顺超限大值所在的曲线线型（缓和曲线、圆曲线、直线）往往存在误判。同时，对曲线的起始位置、平均半径、长度等参数检测的准确性低、偏差大，使得系统提供的曲线摘要报表和数据利用率较低。国内在役的轨道检测系统使用步长跟踪曲线线型的方法计算曲线相关参数，但是该方法受跟踪步长和设定阈值参数的影响，对不同大小的半径不能自适应选择参数，检测准确性达不到要求。侯广东2提出拟合迭代算法识别曲线关键参数，该方法受迭代终止阈值的影响，线型分段自动化能力有限。熊丽娟等3提出使用霍夫变换、稳健估计和迭代法三种基本算法对曲线关键参数先进行识别，再进行拟合与迭代，使用轨道检查仪采集数据进行后处理，实现较为复杂，不适用于车载轨道检测实时系统。曲线参数检测时，只要确定曲线的直缓点、缓圆点、圆缓点、缓直点四大特征点的准确位置，就能准确计算曲线平均半径、长度、线型等参数。轨道检测系统曲线检测数据具有明显的时间和空间特征，空间上无交叉和重复，与地理测绘领域中运动轨迹数据特征提取类似4。轨迹特征提取算法有两种：线性分段法5，算法简单且计算复杂度低；非线性拟合法6，其结果更接近真实轨迹，但算法复杂，计算量大。在线性分段轨迹特征提取的经典算法中，DP算法7最具代表性。该算法是一种线性压缩算法，能有效提取曲线轨迹特征点，运算复杂度较低。其主要思想是将复杂线段进行抽稀和压缩8，从而寻找出线型变化处的特征点，但是DP算法在局部细节上不能很好地反映曲线走向特征。本文根据轨道检测系统曲线检测数据特点，提出在经典DP算法粗略确定曲线特征点的基础上，在特征点周围窗口内使用角度阈值法提取角点，再将各个角点逐一连线成为走向段，设定阈值寻找极值走向段，根据极值段筛选出曲线线型变化处的全局特征点，实现对曲线相关参数的准确检测。收稿日期：20230321；修回日期：20230719基金项目：中国国家铁路集团有限公司科技研究开发计划（P2021G014）第一作者：李颖(1979 )，女，研究员，博士。E-mail：铁道建筑第 63 卷1 轨检系统曲线检测数据特点 车载轨道检测系统通过安装在车体或检测梁内的惯性和位移传感器来测量曲线相关参数，要求传感器信号带宽不能太高，以免噪声信号叠加到线型数据中，干扰曲线特征点的准确获取。目前在役的轨道检测系统能输出曲率、曲率变化率、曲线正矢等曲线相关检测波形，在利用这些数据之前需要对信号进行预处理，降低噪声信号的影响。在多种信号滤波方法中，选取能有效去除噪声又能保持曲线线型不失真的线性高斯滤波算法对传感器信号进行预处理滤波。高斯滤波窗口长度不宜设置过大，按照轨检系统空间0.25 m采样间隔，一般选取100 120 m窗口长度，以防止曲线特征点被滤除，同时可以满足检测系统延时输出最大长度限制。本文选取窗口长度120 m进行高斯滤波（图1），可以看出轨道曲线线型变化的特点，整体形状类似于互相间隔的梯形，存在变形并含有噪声，使得底边与斜边的转折点不易准确确定。2 DP算法提取特征点 车载轨道检测系统虽然是在线实时检测系统，但内部需要根据数学模型作计算和滤波处理，运算需要积累一定长度的区段数据，而且曲线摘要报表的输出是按各铁路局或者工务段管辖区段进行统计，因此可以应用DP算法对曲线线型数据进行特征点提取。DP算法的基本思路8为：设置一定长度（L）区段，连接区段的首尾两点形成一条直线，计算区段内所有点到直线的距离，将最大距离（Dmax）与设定的阈值（D）进行比较，若Dmax D则保留此坐标点，并以该点为界，将区段分为两部分，对这两部分递归使用该方法。原理见图2。该算法能够删除曲线上的非特征点，寻找最有代表性的特征点，不足之处是在局部细节上不能很好地反映曲线的走向特征。如图3所示，对一铁路轨道部分曲线数据使用纯DP算法进行特征点提取。该段线路设计半径5 000 12 000 m。首先设定距离阈值（D），使用经典DP算法得到曲线的粗略特征点和轨迹，D可以稍大，使DP算法迭代次数减少，提高计算速度。L根据需要统计或检测的里程确定，D根据曲率数值选取。在动态轨道检测系统中曲率的单位为 rad/km，对于半径 5 000 12 000 m的区段，可以取D=0.02 km。由图3可以看出，使用DP算法能够提取特征点并确定曲线线型变化轨迹，但直缓点、缓圆点、圆缓点、缓直点四个线型变化特征点的位置存在偏差。特征点位置偏差不仅使检测系统提供的曲线起始位置、曲线长度、平均半径、限速半径等曲线参数的检测结果产生较大误差，还会对高低、轨向、轨距、水平、三角坑等轨道几何参数超限峰值所在的线型出现误报，因此须要对算法进行改进。3 曲线线型变化特征点提取算法 使用DP算法能够从全局出发确定轨迹线的特征点，保留曲线数据的整体形态，但是对四大特征点位置的确定不够精准，须要进一步采取局部特征提取方法。因此，在使用DP算法确定粗略特征点的一定宽度窗口内，依据角度阈值法筛选角点9，再将所有角点连线成为走向段，设定阈值寻找极值走向段，从而准确得到曲线线型变化处四大特征点。图1曲线检测数据预处理图2DP算法提取特征点原理图3使用DP算法提取曲线轨迹与特征点2第 10 期李颖等：轨道曲线参数检测算法优化3.1角点提取角点是计算机视觉领域提出的概念，可应用于轮廓曲线的特征点提取，其定义为曲线上局部曲率最大值点10。本文对DP算法确定的粗略特征点窗口内角点提取方法为：设定角度阈值（），依次顺序获取窗口内曲线数据的3个点，计算相邻两点连线的夹角i1和i2；根据式（1），将i与进行比较，若i 则将该点放入角点集合（S）中。i=|i1+|i2 i1i2 0|i1-i2 i1i2 0（1）设定曲线上窗口内数据点之间角度变化值作为选择角点的依据。用绝对角度来度量曲线轨迹变化，规避了丢失曲率变化特征点的问题，同时也避免了使用斜率差提取角点方法带来的非线性弊端。3.2曲线线型变化特征点提取轨道检测系统的曲线检测数据虽然通过多种信号处理，但是仍然存在噪声，而且线路上曲线的圆顺性并不理想，存在大值波动。噪声和大值波动的存在容易使曲线趋势特征被破坏并形成多个特征点11。通过角度阈值法筛选的曲线角点集合（S)，可能由局部特征点、转折点、极值点、波动较大的非极值点等组成的，因此要再次对这些角点进行评价和筛选，确定准确的曲线四大特征点。将S中相邻的两个角点逐一连线成为走向段，使寻找表征曲线线型变化特征点的问题转化为寻找极值走向段的问题，采用顺序对比走向段直至确定极值段12，该极值段对应的角点即为曲线线型变化四大特征点。将角点集合S=s1，s2，sN中两个相邻角点连线，以其平均值tj=(sj+sj+1)/2代表该走向段，表征该段的趋势程度，tj T，其中 T 为走向段集合，T=t1，t2，tM。当轨道曲线线型转变时会存在一个极值段，设定阈值（)对极值段进行筛选。筛选条件为：tj tj-1 tj tj-1 tj tj+1，|tj-tj-1|或|tj-tj+1|。当极值段为极大值时，极值段角点大的点为曲线线型变化特征点；当极值段为极小值时，极值段角点小的点为曲线线型变化特征点。曲线线型变化四大特征点提取算法步骤：设置DP算法距离阈值，粗略特征点周围窗口宽度、角度阈值、走向段极值阈值。轨道曲线检测原始信号数据采用高斯滤波进行预处理。使用DP算法确定粗略特征点。由于此步为粗筛，对距离阈值的选取可以稍微增大，以降低DP算法的迭代次数，提高运算时间。在已确定的窗口宽度内提取角点。根据式（1）依次计算窗口内曲线相邻数据点的角度变化值，当角度变化大于角度阈值时，把该点加入角点集合。将相邻角点连线为走向段，以角点平均值建立走向段集合。确定极值走向段并确定曲线线型变化全局特征点。根据筛选条件删除走向段集合内不满足条件的区段，得到极值走向段。极值段为极大值或极小值时，取相应极值段角点大或小的点为曲线线型变化最终特征点。由于DP算法能够确定曲线特征点的粗略位置，结合TB 100982017 铁路线路设计规范 对5 000 m以上大半径曲线线路中缓和曲线的设计要求，窗口宽度取50 m较为合适；角度阈值的选取根据自适应阈值法，依次计算窗口宽度范围内相邻两点连线夹角的平均值作为角度阈值；走向段极值阈值的选取可以稍大，避免特征点具有局部性，同时减少计算量。算法流程见图4。4 算法验证 将本文提出的算法应用到车载轨道检测系统曲线参数检测中，与只采用DP算法提取的特征点进行对比，见图5。可知：纯DP算法提取的特征点往往受噪声和大值波动的干扰，线型变化特征点位置不够准确；本文算法提取的线型变化特征点在DP算法的基础上进行了修正，准确性得到提高。对于图3曲线数据，使用本文算法对四大特征点位置进行修正（图6），可以明确看出曲线线型变化处，圆缓点、缓直点位置更加准确。为了进一步验证算法的准确性，使用本文提出的图4算法流程3铁道建筑第 63 卷算法与在役轨道检测系统现有算法、纯DP算法进行对比，对某高速铁路和普速铁路线路的曲线进行检测，部分曲线参数检测数据见表1。其中台账为对应线路上曲线标准设计半径和长度参数。以台账数据为参考标准值，根据表1数据进行准确率统计对比，见表2。本文不考虑检测数据和台账之间的里程偏差。由表1和表2可知：对半径5 000 12 000 m的高速铁路线路大半径曲线，轨道检测系统曲线检测现有算法无法检测出数据，使用纯DP算法对曲线半径和长度参数检测的平均准确率为93.5%，使用本文算法可达95.8%；对半径5 000 m以下的普速铁路线路，现有算法对曲线半径和长度检测的平均准确率为91.1%，使用纯DP算法为93.7%，使用本文算法可达95.8%。经过大量的数据测试和验证，本文算法对曲线参数检测的平均准确率能达到95%以上。图5部分曲线数据特征点提取图6本文算法提取曲线轨迹与特征点表1三种算法曲线参数检测结果对比m线路高速铁路线路普速铁路线路曲线台账半径6 00011 00012 0008 00512 0008 0009 9957 0009 9958 0052 5303 0001 5405 0008131 1004 000500450长度7872 0314061 8961 2301 8571 49724761 3722 3475369233416349839769400844现有算法半径2 8573 2261 6395 8828471 0874 762493450长度5298323093849136857410855纯DP算法半径5 62110 01312 6978 26612 7428 13811 1766 78010 9438 3502 7022 7891 4365 4657711 0224 139495459长度8112 2454382 0141 2941 9251 3992 2141 4332 5955591 0263836051741977430867本文算法半径5 73711 65211 5347 88712 4617 89310 6246 7869 7958 2672 6932 9031 5925 4738071 1264 118495458长度7632 1174221 7761 2851 7931 3782 2221 3422 3195018863576046638665421867表2三种算法曲线参数检测准确率对比曲线半径/m5 000 12 000 5 000算法现有算法纯DP算法本文算法现有算法纯DP算法本文算法准确率/%半径无法检出93.996.492.294.696.5长度无法检出93.195.390.092.995.2平均93.595.891.193.795.84第 10 期李颖等：轨道曲线参数检测算法优化5 结论 本文在轨道检测系统曲线参数检测中引入经典DP算法提取曲线线型变换粗略特征点，在粗略特征点一定窗口宽度内使用角度阈值法提取角点，将角点连线成走向段，通过寻找极值走向段的方法确定曲线线型变化处准确特征点，根据准确的四大特征点位置来计算曲线平均半径、长度等曲线参数。本文算法能够准确确定轨道曲线线型变换四大特征点的位置，对大半径和小半径曲线的检测准确率均达到95%以上，具有良好的适应性，能够更好地为铁路工务部门服务。参考文献1 戚妍娟.高速铁路大半径曲线参数检测方法研究 J.铁道建筑，2019，59（6）：109-111.2 侯广东.铁路曲线特征点提取算法研究 D.成都：西南交通大学，2015.3 熊丽娟，朱洪涛，魏晖，等.基于轨道相对测量数据识别曲线关键参数的研究 J.铁道学报，2018，40（4）：127-134.4 张达夫，张昕明.基于时空特性的GPS轨迹数据压缩算法J.交通信息与安全，2013（3）：6-9.5 梁明，陈文静，段平，等.轨迹压缩的典型方法评价 J.测绘通报，2019（4）：60-646 李卫东，高樱莺.列车轨道曲线拟合算法改进研究 J.计算机与数字工程，2019，47（6）：1326-1330.7 DOUGLAS D H，PEUCKER T K.Algorithms for the Deduction of the Number of Points Required to Represent a Digitized Line or its CaricatureJ.Canadian Journal of Cardiology，1973，10（2）：112-1228 韩晓霞，崔浩，孙钰珊.顾及曲线走向及局部面积特征的矢量数据压缩算法 J.北京测绘，2017，31（6）：6-9.9 周腾，李晓妍.改进的道格拉斯-普克算法在电信代维系统的应用 J.计算机与数字工程，2020，48（7）：1576-1579.10 刘晓红，李树军.矢量数据压缩的角度分段道格拉斯算法研究 J.四川测绘，2005，28（2）：51-52.11 朱思聪，周德龙.角点检测技术综述 J.计算机系统应用，2020，29（1）：22-28.12 王大维，汪浩.一种剔除GPS坐标时间序列非线性趋势项的方法 J.测绘与空间地理信息，2019，42（11）：142-144.Optimization of Algorithm for Detecting Track Curve ParametersLI Ying1，CHENG Zhaoyang1，LU Yongbiao2，CHEN Shiming1，YU Ning1，QIN Zhe1，NA Qiang21.Infrastructure Inspection Research Institute，CARS，Beijing 100081，China；2.Beijing IMAP Technology Co.Ltd.，Beijing 100081，ChinaAbstract The existing track detection systems have low accuracy in detecting curve related parameters，especially for large radius curves above 5 000 m，which results in low utilization of curve data reports.In response to this issue，it was proposed to use the classic DP（Douglas-Peucker）algorithm to extract rough curve line transformation feature points，and then use the angle threshold method to extract corner points within the window around the feature points.The corner points were connected into a trend segment，and the method of sequentially searching for the extreme value trend segment was used to determine the four major feature points of curve line transformation，between straight，gentle circle，circular，and gentle straight sections，which could accurately calculate the average radius，length，and other parameters of the curve，and improve the accuracy of system curve detection.The algorithm in this article was applied to the curve parameter detection of the vehicle mounted track detection system，which was compared with existing algorithms and feature points extracted only using DP algorithm.The results show that the accuracy of using the algorithm presented in this paper in detecting parameters of large and small radius curves can reach over 95%，and the algorithm in this paper has good adaptability.Key words track geometry；curve parameters；feature point extraction；DP algorithm；corner point detection；extreme value sectionCitation format：LI Ying，CHENG Zhaoyang，LU Yongbiao，et al.Optimization of Algorithm for Detecting Track Curve Parameters J.Railway Engineering，2023，63（10）：15.（编辑：苗蕾 校对：李付军）5