Segmented fitting method of line features of digital track map
-
摘要: 基于全球卫星导航系统(GNSS ,Global Navigation Satellite System)的列车定位技术降低了铁路建设与运营的成本,数字轨道地图作为地图匹配方法的基础,通常作为辅助定位手段对卫星定位偏差进行校正,对轨道线形进行准确描述及约简是轨道电子地图构建存储的前提。文章采用方位角近似估计曲率的方法对轨迹平面线形进行初步识别分类,为提高线形的拟合精度,引入自适应权值因子,对最小二乘拟合法进行优化,以消除大偏差观测数据对拟合趋势的影响,从而保证整体轨迹数据的平滑拟合及精度。实测车辆轨迹数据验证表明,该算法可对车辆轨迹线路特征进行有效提取及准确描述。
-
关键词:
- 全球卫星导航系统(GNSS) /
- 数字轨道地图 /
- 平面线形 /
- 分段曲线拟合 /
- 方位角
Abstract: The train positioning technology based on the Global Navigation Satellite System (GNSS) reduces the cost of railway construction and operation. As the basis of map matching method, digital orbit map is usually used as an auxiliary positioning means to correct satellite positioning deviation. Accurate description and reduction of track alignment are the premise of the construction and storage of track electronic map. This paper used the azimuth approximate estimation of curvature to preliminarily identify and classify the track plane alignment. In order to improve the fitting accuracy of the line shape, the paper introduced an adaptive weight factor to optimize the least squares fitting method, and eliminate the influence of large deviation observation data, thereby ensure the smooth fitting and accuracy of the overall trajectory data. The verification of measured vehicle trajectory data shows that the algorithm can effectively extract and accurately describe the characteristics of vehicle trajectory. -
铁路作为国家经济的动脉,承载着推动社会进步与服务经济建设的关键使命。保障铁路安全,是铁路运营的首要职责。依据《铁路安全管理条例》[1],铁路线路两侧必须设立安全保护区,确保铁路安全运营和维护。对于保护区内已存在的威胁铁路运输安全的隐患,必须采取有效防护措施或依法予以清除,以确保铁路运行安全。
彩钢瓦隐患作为铁路周边的主要隐患之一,不容忽视。鉴于我国铁路沿线环境错综复杂,铁路管理部门遵循预防为主的方针,采取定期巡检的举措,旨在及时发现并消除彩钢瓦等安全隐患。当前,铁路周边环境的巡检工作仍主要依赖人力。人工巡检不仅消耗大量的人力资源,且效率低下、周期冗长,无法实现对彩钢瓦隐患的即时发现与迅速应对。
计算机视觉技术的快速发展,使学者们找到了解决上述问题的方向,诞生了大量针对各类隐患的识别模型。黄一昕等人[2]改进了Mask R-CNN实例分割模型,用于提取铁路遥感影像中的房屋,并采用边缘提取算法自动对遥感影像进行矢量化提取,有效生成铁路周边的房屋矢量化影像;吴云鹏等人[3]采用图像增强和缺陷分割技术,通过无人机拍摄轨道图像,提出可突出和均匀化钢轨表面缺陷的图像增强技术,有效检测铁路轨道表面缺陷。
借鉴上述研究成果,本文基于改进YOLO(You Only Look Once)v8模型,结合高分辨率遥感图像技术,设计了铁路沿线彩钢瓦隐患识别方法。针对彩钢瓦结构微小、细节特征难以捕捉的难题,在YOLO(YouOnly Look Once)v8模型中添加区域损失(Area Loss)函数,从而提高彩钢瓦隐患的识别精度。
1 YOLOv8模型改进
YOLO自2015年被提出以来[4],已迭代到了YOLOv8版本。因其兼顾目标检测速度和精度的优点,越来越多研究者使用YOLOv8模型设计目标检测算法[5-7]。由于彩钢瓦在遥感图像中像素过小,本文对YOLOv8模型进行改进,在YOLOv8模型中添加一个Area Loss函数,加强对于小目标的识别性能。改进后的YOLOv8模型(简称:改进模型)网络结构如图1所示,由Backbone、Neck和Head 组成。
1.1 YOLOv8模型改进
1.1.1 Backbone
作为图像特征提取的核心组件,Backbone网络通过运用卷积与反卷积层的组合,并辅以残差连接和瓶颈结构,简化了网络规模,提高了模型的准确性。它基于CSPDarknet53主干网络,保留了跨阶段 (CSP) 模块的精髓,同时用 C2F模块替换了 YOLOv5 中的 C3 模块,使其拥有更加丰富的梯度流信息和更高的识别精度。此外,Backbone还融合了深度可分离卷积和膨胀卷积等改进技术,进一步增强了模型的特征提取能力,为后续处理环节提供了更为精准的特征输入。
1.1.2 Neck
Neck作为图像融合的核心,用于整合Backbone各阶段产出的多尺度特征图,从而显著提升特征的表达力和信息丰富度。它包含了路径聚合网络(PAN,Path Aggregation Network)和特征金字塔网络(FPN,Feature Pyramid Network),旨在克服多级卷积可能导致的特征信息衰减问题,确保模型能够全面捕捉图像细节。FPN通过自上而下与自下而上的双向路径,构建多层次的特征金字塔,有效提取图像的多层次语义信息。同时,引入的横向连接机制,促进了自下而上与自上而下路径间的信息交流,显著增强了跨层特征的融合效果。然而,FPN的融合策略在可能会导致图像特征信息的局部缺失。为此,又引入PAN,通过高效且精准的方式融合高低层次特征,从而有效弥补FPN可能存在的信息缺失问题。
1.1.3 Head
Head负责生成最终的检测结果,由3个独立的检测分支组成,每个分支都能够独立地输出目标分类的概率分布和精确的目标边界框坐标。为进一步优化检测性能并减少超参数依赖,Head采用了Anchor-free机制。该机制摒弃了传统的锚框设定,通过直接预测目标中心点的位置,显著简化了训练流程,并提升了模型的泛化能力。
Head中包含2种损失函数:本文添加的Area Loss函数和原有的分类损失(Cls Loss)函数。Area Loss函数旨在提高改进模型在识别小目标时的准确率;Cls Loss函数用于衡量模型对每个候选框预测的类别标签与真实类别的差异。
1.1.4 损失函数改进
彩钢瓦在遥感图像所占像素较小,相比于小目标物体,YOLOv8模型更容易识别一些大目标物体[8]。为此,本文添加的Area Loss函数,可以使改进模型更关注于小目标特征,从而提升彩钢瓦识别的精度。改进模型的损失函数公式为
Loss=1N(N∑i=1ωiLcls(pi,p∗i)+N∑i=1ωitiLals(gi,g∗i)) (1) 式(1)中,
i 为一张图像上的第i 个彩钢瓦目标;N 为一张图像上的彩钢瓦总数;pi 为物体被预测为是彩钢瓦的概率;p∗i 为物体真实的标签;gi 是包含4个参数的向量,表示预测的彩钢瓦坐标;g∗i 也是包含4个参数的向量,表示真实的彩钢瓦坐标;ti 表示预测彩钢瓦的框和真实框面积的重叠程度,1表示重叠面积大于50%,0表示重叠面积小于50%;Lcls 为Cls Loss函数;Lals 为本文添加的Area Loss函数;ωi 表示锚框的权重系数,为使模型提高小目标识别的准确度,其取值与锚框的面积呈反比,即锚框的面积越小,其取值越大。由于YOLO算法中对锚框的长和宽进行了归一化操作,导致锚框的面积小于1。因此本文提出了一个单调递减函数,使
ωi 的值始终大于1,其公式为ωi=esi+1 (2) 公式(2)中,
si 表示第i个锚框的面积。将公式(2)代入公式(1)可得到
Loss=1N(N∑i=1(e−si+1)Lcls(pi,p∗i)+N∑i=1(e−si+1)tiLlac(gi,g∗i)) (3) Cls Loss函数和Area Loss函数的公式分别为
Lcls=−(p∗i⋅log(pi)+(1−p∗j)⋅log(1−pi)) (4) Lals =scross isi+s∗i−scross i−OC2−c (5) 式(5)中,
Si 表示第i个锚框的真实面积;s∗i 表示预测第i个锚框的面积;scrossi 表示si 与s∗i 相交的面积;O 表示真实的锚框与预测锚框中心点之间的欧氏距离;C 表示真实的锚框与预测锚框之间外接矩形的对角线长度;c 表示真实锚框与预测锚框之间长宽比的相似程度,越相似则值越小。1.2 实验数据处理
1.2.1 数据来源
本文实验采用的遥感图像数据来自于国家遥感数据与应用平台中的高分二号遥感(GF2)图像。共下载卫星遥感图片974景,覆盖铁路的为747景。GF2的空间分辨率为1 m,共有多光谱和全色两种数据。其中,多光谱波段的空间分辨率为4 m,包含蓝光、绿光微米和近红外波段的数据;全色的分辨率为1 m,包含可见光和近红外波段的数据。
1.2.2 数据处理
由于GF2图像为波段影像,其影像数据包括多光谱数据和全色数据,需要对这两种数据进行处理,使其转变为正常的图像数据[9]。数据处理流程如图2所示。
融合GF2的两种数据前需要对其进行辐射定标,以消除传感器本身的误差,确定传感器入口处的准确辐射值;再对多光谱数据进行大气校正,获得地物反射率和辐射率、地表温度等真实物理模型参数,消除大气中的水蒸气、氧气、二氧化碳、甲烷和臭氧对地物反射的影响,同时,消除大气分子和气溶胶散射的影响;将处理后的多光谱数据与全色数据进行融合,得到更清晰、信息更丰富的图像。图3为遥感影像的处理流程示例。
1.2.3 数据标注
数据处理后,将卫星遥感影像裁剪为640×640像素的图片,一景遥感图像可裁出约800张遥感图片,对裁剪完成的图片,使用LabelImg软件进行标注,构建多数据源、多尺度的彩钢瓦数据集,共标注彩钢瓦图片
17816 张。彩钢瓦标注如图4所示。2 实验与结果分析
2.1 评估指标
本文使用平均精度均值(mAP)对彩钢瓦隐患的检测性能进行评估。mAP分为mAP50和mAP50-95两种类型[10]。mAP50表示在模型预测框与真实框交并比(IoU)大于50%时,即检测框与真实框之间的重叠程度大于50%时,检测结果正确的比例。mAP50-95为当IoU阈值从50%到95%时,检测结果正确比例的均值。由于模型准确率会随着IoU阈值的增大而降低,所以mAP50-95指标更为严格[11]。
本文利用的其他常用的评估指标包括准确度(Accuracy)和F1值等[12],计算公式为
Accuracy=TP+TNTP+TN+FP+FN (6) Precision=TPTP+FP (7) Recall=TPTP+FN (8) F1−Score=2(Precision×Recall)Precision+Recall (9) 式(6)~式(9)中,TP表示被模型预测为彩钢瓦隐患且实际是彩钢瓦隐患的个数;TN表示被模型预测为不是彩钢瓦隐患且实际上确实不是彩钢瓦隐患个数;FP表示被模型预测为彩钢瓦隐患但实际上不是彩钢瓦隐患个数;FN表示被模型预测为不是彩钢瓦隐患但实际上是彩钢瓦隐患个数。Accuracy表示模型预测的结果中,预测正确的比例;Precision表示模型预测为彩钢瓦隐患样本中,预测正确的比例;Recall表示模型预测为彩钢瓦隐患的样本,占实际上所有的彩钢瓦隐患的比例。
2.2 实验参数设置
2.2.1 硬件参数与环境配置
本文使用的CPU型号为Inter(R)-Xeon(R)-W1270 3.40 GHz;GPU型号为Nvidia GeForce RTX 3070,显存为8 GB;计算机内存为32 GB;计算机操作系统为Ubuntu20.04;Anaconda版本为23.7.3;Python版本为3.8.16;PyTorch版本为2.0.1。
2.2.2 模型参数
本文共标注遥感图片
17816 张,以8∶2的比例划分训练集和测试集,其中,训练集有15646 张,测试集有2170 张,YOLOv8模型预训练尺寸为YOLOv8x。训练参数设置:epoch设为1000, optimizer设为Adam,初始学习率设为0.001,momentum设为0.937,weight_decay设为0.0005, patience设为100,batch_size设为4。2.3 实验结果
为直观展示本文提出的Area Loss函数的有效性,将未改进的YOLOv8模型的方法和改进模型(YOLOv8+AreaLoss)的方法进行对比实验,方法评估指标如表1所示。
表 1 方法评估指标方法 F1值 mAP50 mAP50-95 YOLOv8 0.786 0.794 0.533 YOLOv8+AreaLoss 0.805 0.820 0.558 如表1所示,与采用未改进的YOLOv8模型相比,采用改进模型的铁路沿线彩钢瓦隐患识别方法,F1值、mAP50和mAP50-95检测指标分别提高了0.019、0.026和0.025。
采用未改进的YOLOv8模型与采用改进模型的实际识别效果如图5所示。图5中,第1列为待检测的原始图像,第2列为采用未改进的YOLOv8模型的识别效果,第3列为采用改进模型的识别效果。
由图5可知,本文设计的方法显著提高了对小目标的检测能力,鲁棒性和准确性得到了显著增强,有效减少了误检和漏检的情况。
3 结束语
本文基于改进YOLOv8模型,设计了铁路沿线彩钢瓦隐患识别方法。改进模型中的Area Loss函数,优化了模型的识别性能,对特征图中面积较小的彩钢瓦隐患给予更多的关注;在自建的彩钢瓦数据集上进行的实验验证了该方法的有效性与优越性;该方法已在中国铁路郑州局集团有限公司铁路工务日常巡检中应用,节约了人力成本和时间成本,为优化铁路领域传统的人工巡检流程提供了技术支撑。
下一步,将考虑改变原有YOLOv8的模型结构,如添加注意力机制等,进一步提升识别的准确性。
-
表 1 铁路线形平面参数特点
平面线形 直线 缓和曲线 圆曲线 半径 / R→∞ R 方位角 定值 二次抛物线变化 线性递增/递减 曲率 恒为0 线性变化 1/R 表 2 算法性能指标对比
性能指标 标准LS 改进LS AVG 1.0421 0.8213 MAX 1.6542 0.9356 VAR 0.3524 0.0572 MSE 0.5446 0.3741 -
[1] OTEGUI J, BAHILLO A, LOPETEGI I, et al. A survey of train positioning solutions [J]. IEEE Sensors Journal, 2017, 17(20): 6788-6797. DOI: 10.1109/JSEN.2017.2747137
[2] 夏 青,杨剑锋. 数字轨道地图线路分段表示方法研究 [J]. 铁路计算机应用,2012,21(4):9-11. DOI: 10.3969/j.issn.1005-8451.2012.04.003 [3] 倪晓林,杨 妮,张 东. 列控数字轨道地图自动生成算法研究 [J]. 铁道运营技术,2010,16(1):24-27. DOI: 10.3969/j.issn.1006-8686.2010.01.009 [4] 左自辉,王开锋,许 聪,等. 生成高精度轨道电子地图的数据处理方法 [J]. 中国铁道科学,2016,37(4):134-138. DOI: 10.3969/j.issn.1001-4632.2016.04.21 [5] 龚 利,沙建领,赵 超,等. 基于北斗定位技术的铁路机车轨道电子地图生成算法 [J]. 计算机系统应用,2019,28(6):172-177. DOI: 10.15888/j.cnki.csa.006945 [6] 刘 江,蔡伯根,唐 涛,等. 基于GPS的列控轨道地图数据生成方法研究 [J]. 测绘学报,2011,40(1):111-117. [7] 陈德旺,蔡伯根,王 剑,等. 轨道交通GPS数据约简的数学模型与算法研究 [J]. 铁道学报,2008(4):116-119. DOI: 10.3321/j.issn:1001-8360.2008.04.023 [8] 曾 强,陈德旺,王丽娟,等. 基于主曲线和自适应半径的多GPS轨迹数据融合算法 [J]. 铁道学报,2015,37(2):46-51. DOI: 10.3969/j.issn.1001-8360.2015.02.007 [9] 陈德旺. 铁路卫星定位数据的智能处理[M]. 北京: 北京交通大学出版社, 2014. [10] Camacho-Torregrosa F J, AM Pérez-Zuriaga, JM Campoy-Ungría, et al. Use of Heading Direction for Recreating the Horizontal Alignment of an Existing Road [J]. Computer-Aided Civil and Infrastructure Engineering, 2015, 30(4): 282-299. DOI: 10.1111/mice.12094
[11] Winter H , Luthardt S , Willert V , et al. Generating Compact Geometric Track-Maps for Train Positioning Applications[C]//2019 IEEE Intelligent Vehicles Symposium (IV), June 09-12, 2019, Paris, France. New York, USA: IEEE, 2019.
[12] 陶维杰,蔡伯根,王 剑,等. 数字轨道地图平面线形特征提取方法研究 [J]. 铁道学报,2019,41(9):81-87. DOI: 10.3969/j.issn.1001-8360.2019.09.011 [13] 陈明晶,方源敏,陈 杰. 最小二乘法和迭代法圆曲线拟合 [J]. 测绘科学,2016,41(1):194-197,202. DOI: 10.16251/j.cnki.1009-2307.2016.01.038