CN100462678C - 钢轨磨耗激光视觉动态测量方法 - Google Patents

钢轨磨耗激光视觉动态测量方法 Download PDF

Info

Publication number
CN100462678C
CN100462678C CNB2005101237250A CN200510123725A CN100462678C CN 100462678 C CN100462678 C CN 100462678C CN B2005101237250 A CNB2005101237250 A CN B2005101237250A CN 200510123725 A CN200510123725 A CN 200510123725A CN 100462678 C CN100462678 C CN 100462678C
Authority
CN
China
Prior art keywords
rail
measurement
profile
coordinate
laser
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CNB2005101237250A
Other languages
English (en)
Other versions
CN1776364A (zh
Inventor
周富强
张广军
朱奎义
闫伟
江洁
魏振忠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CNB2005101237250A priority Critical patent/CN100462678C/zh
Publication of CN1776364A publication Critical patent/CN1776364A/zh
Application granted granted Critical
Publication of CN100462678C publication Critical patent/CN100462678C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明属于测量技术领域,涉及一种钢轨磨耗综合参数的非接触动态测量测量装置及测量方法。本发明测量装置由激光视觉传感器7、计算机、高速图像采集卡、标定靶标和测量软件等组成。本发明测量方法由激光投射器3投射出光平面,经钢轨断面成像***6成像后,由高速图像采集卡采集到计算机中,对采集图像进行处理,提取出断面特征轮廓的图像坐标,根据测量模型,计算出垂直磨耗、侧边磨耗和总磨耗。本发明首次实现了单激光视觉传感器同时测量单侧钢轨断面轮廓、垂直磨耗、侧面磨耗和总磨耗参数;实现了动态条件下对钢轨磨耗的现场测量。

Description

钢轨磨耗激光视觉动态测量方法
技术领域
本发明属于测量技术领域,涉及一种钢轨磨耗综合参数的非接触动态测量方法。
背景技术
铁路的状况与安全性密切相关,特别是对于高速列车更加重要。钢轨的定期检测对于合理计划和保持低的维护成本是非常重要的。在钢轨磨损和变形早期进行测量有助于制定更好的维护时间表,以避免等观测到危险状况时再作紧急处理,另一方面,进行提前维护使有限的人力和仪器设备资源得到更好的利用,从而降低了成本。此外,正确的维护减少了因车轮与钢轨之间不良接触所产生的噪声污染。近年来我国铁路提速和重载运输的大力发展、客货运量和行车密度大幅度增长,导致钢轨磨耗日益严重。如何高效而精确地测量钢轨磨耗值一直是国内外铁道部门的一个重要研究课题。
为了保证铁路运行安全,要求在不影响火车运行的条件下,每个月要对钢轨至少进行一次检修。九五期间,我国铁路钢轨里程达到68000多公里,至今没有合适的测量手段。长期以来,国内对钢轨磨耗的测量都是采用手工测量,用专用卡尺人工抽样检测其断面损耗。这种方式效率低下,耗费大量人力物力,而且处于高度危险中。由于在测量中不可避免地引入了测量者的人为因素,直接影响了测量的精度和可靠性。手工测量已经不能适应我国铁路事业的发展。因此,研制一种快速准确的钢轨磨耗自动测量装置具有很大的现实意义。
国内外许多研究机构和学者对钢轨磨耗的检测方法进行了研究,并成功研制了各种测量装置,根据检测方式的不同,大致可以分为接触式和非接触式检测。美国GreenWood公司研制的MiniProf专用量规钢轨磨耗测量装置,采用接触式方式逐点测量钢轨轮廓,然后与标准轮廓进行比对。Wroclaw科技大学研制了类似的接触式高精度钢轨断面磨耗测量***。在国内,哈尔滨铁路计量局研制了一种类似国外的卡尺式钢轨磨耗测量装置。接触式钢轨磨耗测量***测量精度高,但只适合静态测量,操作复杂,测量效率低,主要用于实验室研究钢轨耐磨性能时使用,不适合在线测量。
北美和欧洲发达国家从20世纪80年代初期就开始利用计算机进行自动、精确钢轨磨耗测量的研究。八十年代澳大利亚率先研制了LitesLice固体传感器钢轨断面轮廓自动测量***。1984年后相继在澳大利亚、美国、加拿大和英国等投入使用,并不断完善,逐渐成为世界上使用最为普遍的钢轨断面磨耗自动测量***。该方法的缺陷在于,采用图像上边缘和左边缘为定位基准,因此要求测量传感器不能有较大的振动存在,因而比较适合于静态测量。实际上,安装在机车上的测量传感器必然存在较大振动,由此影响到测量结果的精度和可靠性。此外,国外有采用两套测量传感器分别测量同一钢轨的两侧,采用没有磨耗的一侧测量结果作为基准,该方法虽然解决振动对测量结果的影响,但成本高,且需要两套测量传感器坐标系的统一,操作复杂。目前,国内研究的钢轨磨耗检测***基本上与国外类似,而且由于测量坐标系到标准(设计)坐标系的转换环节没有很好的解决,因此实际研究的装置仅仅用于静态实验研究,而没有能够解决动态条件下钢轨磨耗的自动测量。
发明内容
本发明所要解决的技术问题是:提供一种精度较高,测量速度快,单一传感器能够同时测量单侧钢轨垂直磨耗、侧面磨耗以及总磨耗等综合参数的在线动态实时测量方法及***装置。以减小测量设备的体积,降低测量设备成本,提高测量效率,改善其工程化应用的可操作性和便捷性。
本发明的技术解决方案是:钢轨磨耗激光视觉动态测量方法,使用钢轨磨耗激光视觉动态测量装置进行综合参数的动态测量,所说的钢轨磨耗激光视觉动态测量装置由激光投射器3和钢轨断面成像***6组成的激光视觉传感器7、计算机、高速图像采集卡、标定靶标和测量软件组成;激光投射器3由半导体激光器1和光平面发生光学***2组成;钢轨断面成像***6由CCD摄像机4和滤光***5组成,钢轨断面成像***6的光轴和激光投射器3的光轴交汇角为45°~60°;激光视觉传感器7固定在火车底部车体上,与被测钢轨10的距离为300~700mm;激光投射器3投射的光平面8垂直于被测钢轨10并与钢轨断面平行;钢轨断面成像***6的成像***能够观测到投射光平面8与被测钢轨10表面相交形成的轨头和轨腰特征轮廓9;所说的标定靶标为一个二维平面,靶标平面上有预先设置的特征点,在靶标平面上有成矩阵排列的黑色方块,方块的数量为4~100个,方块的边长为3~50mm,其精度为0.01~0.05mm,方块间距为3~50mm,其精度为0.01~0.05mm,选取每个方块的顶点为特征点;其特征在于,测量的具体步骤是:
1、建立激光视觉传感器7的摄像机模型和测量模型,标定摄像机参数和传感器结构参数;
1.1、在光平面上建立二维测量坐标系om-xmym,原点om为光平面上任意点,xm轴和ym轴相互垂直且位于光平面内;
1.2、标定摄像机的内部参数和激光视觉传感器的结构参数并保存到计算机;
1.3、建立标准钢轨二维设计坐标系on-xnyn,以钢轨断面对称轴为yn轴,轨底为xn轴,选取标准钢轨上轨腰部分大圆和小圆的两个圆心以及圆心连线与小圆的交点为对准特征点,分别计算出对准特征点的设计坐标,并保存到计算机,以备测量阶段调用;
2、测量钢轨磨耗;
2.1、将里程表通过串口连接到计算机上;
2.2、设定标准轮廓与测量轮廓对准参数,设定钢轨轮廓最小长度阈值TL,它由轮廓特征点数表示,取值范围是40~100;轨头轮廓与轨腰轮廓间距最小阈值TD,它由轮廓图像的像素坐标表示,取值范围为15~30像素,典型值20像素;保存设定参数,以备后续测量调用;
2.3、测量程序监视串口数据状态,当接收到串口发送的启动命令和里程数据时,保存里程数据并控制摄像机拍摄一幅包括钢轨断面特征轮廓的图像,称为测量图像;
2.4、提取测量图像中钢轨断面特征轮廓的图像坐标,若未提取到两条满足长度大于TL的特征轮廓,或者两条特征轮廓的间距小于TD,则返回到步骤2.3,重新根据串口状态,控制采集钢轨测量图像;若找到两条正确的钢轨特征轮廓,则根据激光视觉传感器7的测量模型,利用特征轮廓的图像坐标计算钢轨断面特征轮廓的测量坐标;
2.5、根据钢轨断面特征轮廓的测量坐标,将测量轮廓自动分割成轨头、轨腰大圆和轨腰小圆三部分。由分割的轨腰大圆和轨腰小圆部分的特征轮廓,采用半径约束方法分别拟合轨腰大圆和轨腰小圆,小圆半径为20mm,大圆半径取400mm,称为60kg对准模式,大圆半径取350mm,称为50kg对准模式;
2.6、分别采用60kg对准模式和50kg对准模式,计算对准特征点的测量坐标,由其中至少2个对准特征点的测量坐标及设计坐标,计算测量坐标系到设计坐标系的二维变换;
2.7、根据在两种模式下求出的测量坐标系到设计坐标系的变换,分别将钢轨断面轨头部分的测量坐标转换到设计坐标系中,直接和标准钢轨比对,计算得到钢轨的垂直磨耗、侧面磨耗和总磨耗;
2.8、比较两种对准模式下的3个磨耗值,选取较小值作为最后的测量结果,并根据所采用的对准模式得到被测钢轨的轨型,将测量结果、轨型以及对应的里程数保存到测量结果文件中;
2.9、向串口发送测量完成信号;
2.10、重复2.3~2.9步,进行下一次测量。
本发明的优点是:
第一、首次实现了单激光视觉传感器同时测量单侧钢轨轮廓以及垂直磨耗、侧面磨耗和总磨耗多参数,降低了设备成本;
第二、自动分割测量钢轨断面轮廓的轨头、轨腰大圆及轨腰小圆,实现了动态条件下对钢轨磨耗的现场测量;
第三、以钢轨断面轨腰大圆及小圆圆心及其衍生特征为对准特征点,自动获得测量基准,实现了钢轨测量轮廓与设计轮廓的自动对准,避免了火车运行过程中振动对测量的影响;
第四、实现了被测钢轨型号自动识别及磨耗实际检测点的里程记录;
第五、测量***装置结构简单、成本低、自动化程度高,测量精度高、速度快,测量过程简单。
第六、钢轨磨耗激光视觉测量***得到工程化应用。
附图说明
图1是钢轨磨耗激光视觉动态测量***基本原理示意图。图1中,1是半导体激光器,2是光平面发生光学***,3是由1和2组成的激光投射器,4是CCD摄像机,5是滤光***,6是由4和5组成的钢轨断面成像***,7是由3和6构成的激光视觉传感器,8是光平面,10是被测钢轨,9是8与10相交形成的钢轨断面特征轮廓9。
图2是激光视觉传感器的测量数学模型。4是CCD摄像机,3是激光投射器,7是由4和3组成的激光视觉传感器,204是图像平面,8是光平面,206是被测物体。
图3是钢轨设计二维坐标系的建立及对准特征点和磨耗测量点的选择示意图。
图4是投射到钢轨表面的特征光条及采集的包含钢轨断面特征轮廓图像。
图5是测量坐标系到设计坐标系的变换计算。
图6是“钢轨磨耗激光视觉动态检测***”软件界面。
图7是50kg右钢轨测量轮廓图像及提取特征轮廓。
图8是50kg右钢轨轮廓在测量坐标系和设计坐标系中的坐标。
图9是测量轮廓与设计轮廓对准后及磨耗计算结果。
图10是同一钢轨磨耗连续540次测量误差曲线图。
图11是测量软件流程图。
具体实施方式
下面对本发明做进一步详细说明。本发明将现代激光技术和计算机视觉技术相结合,使用单激光视觉传感器,在动态条件下,实现了单侧钢轨断面轮廓、垂直磨耗、侧面磨耗和总磨耗多参数的现场自动测量。
本发明钢轨磨耗激光视觉动态测量装置***基本原理如图1所示,它由激光投射器3和钢轨断面成像***6组成的激光视觉传感器7、计算机、高速图像采集卡和测量软件等组成。激光投射器3由半导体激光器1和光平面发生光学***2组成;钢轨断面成像***6由CCD摄像机4和滤光***5组成;钢轨断面成像***6的光轴和激光投射器3的光轴交汇角为45°~60°;激光视觉传感器7固定在火车底部车体上,与被测钢轨10的距离为300~700mm。固定激光视觉传感器7的装置为现有的固定装置或者类似的简单固定机构。
调整激光投射器3,使得投射光平面8垂直于被测钢轨10并与钢轨断面平行;调整钢轨断面成像***6,保证成像***能够观测到投射光平面8与被测钢轨10表面相交形成的轨头和轨腰特征轮廓9;分别固紧激光投射器3和成像***6。激光投射器3和钢轨断面成像***6的调整机构可以采用通用的各种二维机械调整和固紧机构。
所说的标定靶标为一个二维平面,靶标平面上有预先设置的特征点。在靶标平面上有成矩阵排列的黑色方块,方块的数量为4~100个,方块的边长为3~50mm,其精度为0.01~0.05mm,方块间距为3~50mm,其精度为0.01~0.05mm,选取每个方块的顶点为特征点。
本发明测量装置一式两套安装在火车底部,分别测量左右侧两条钢轨的磨耗。
钢轨磨耗激光视觉动态检测的工作原理为:由高强度的激光投射器3投射出光平面8,光平面8对准被测钢轨10的内侧进行照射并垂直于钢轨。光平面8与钢轨10相交,形成钢轨横断面特征轮廓9,经钢轨断面成像***6成像后,由高速图像采集卡采集到计算机内存,对图像进行处理分析得到钢轨断面特征轮廓的图像坐标,根据测量传感器的数学模型,求得钢轨断面的实际轮廓在测量坐标系中的二维坐标。将钢轨轮廓的测量坐标与设计坐标自动对准,根据磨耗测量的规定,计算垂直磨耗、侧边磨耗以及总磨耗。
参见图11,“钢轨磨耗激光视觉动态检测***”软件的主要功能包括:钢轨断面特征轮廓图像的实时采集;特征轮廓图像坐标的自动“亚像素”精度提取;摄像机镜头畸变校正;钢轨表面轮廓测量;测量***的参数标定;轨型(左右钢轨以及规格50kg和60kg)自动识别;测量基准自动建立、钢轨表面测量轮廓与标准轮廓的自动对准;垂直磨耗和侧面磨耗的自动计算;测量磨耗与里程对应等功能。
激光视觉传感器数学模型。
由CCD摄像机4和激光投射器3组成的激光视觉传感器7的数学模型如图2所示。摄像机三维坐标系为oc-xcyczc,图像坐平面204的图像坐标系为ou-xuyu。在光平面8上以一点om为原点,建立三维参考坐标系om-xmymzm,其中om-xmym为二维测量坐标系,光平面8在参考坐标系下的方程为:
zm=0             [1]
设光平面上任意一点P的三维参考坐标为(xm,ym,zm),相应二维测量坐标系为(xm,ym),点P在图像平面204上的透视投影点为p,p的图像坐标为(xu,yu),归一化图像坐标为(xn,yn),则摄像机模型可以表为:
ρ x n y n 1 = R c m T c m x m y m z m 1 = r 1 r 2 r 3 t x r 4 r 5 r 6 t y r 7 r 8 r 9 t z x m y m z m 1 ρ ≠ 0 - - - [ 2 ]
λ x n y n 1 = f x 0 u 0 0 f y v 0 0 0 1 - 1 x u y u 1 = A - 1 x u y u 1 λ ≠ 0 - - - [ 3 ]
其中ρ,λ为比例因子,A为摄像机302的内部参数矩阵,fx,fy为x、y方向上的有效焦距,(u0,v0)为主点坐标。
Figure C200510123725D00103
为3×3正交旋转矩阵,
Figure C200510123725D00104
为3×1平移矢量。
Figure C200510123725D00105
Figure C200510123725D00106
表示了参考坐标系om-xmymzm到摄像机坐标系Oc-xcyczc的变换。
由公式[1]~[3]可得,归一化图像坐标与测量坐标的变换为:
ρ 1 x n y n 1 = r 1 r 2 t x r 4 r 5 t y r 6 r 7 t z x m y m 1 = H x m y m 1 ρ 1 ≠ 0 - - - [ 4 ]
公式[4]表明光平面与图像平面之间的关系可以通过一个3×3矩阵表示。公式[1]~[4]构成了欧式空间上的激光视觉传感器的数学模型。根据公式[3]由光平面上点的图像坐标可以计算归一化图像坐标,然后根据公式[4],计算其对应的测量坐标。
激光视觉传感器模型参数的标定方法参见周富强等的发明专利“一种基于平面靶标的结构光视觉传感器标定方法”,申请号:03142658.1,公开日:2005年1月19日。
根据激光视觉传感器7的测量模型和本发明的测量原理,钢轨磨耗激光视觉动态测量分为测量***模型参数标定和根据模型进行测量两个阶段。
测量***模型参数的标定具体步骤如下:
1、设定激光视觉传感器。激光视觉传感器7为由钢轨断面成像***6和激光投射器3组成。钢轨断面成像***6的光轴和激光投射器3的光轴按成45°~60°角度设置好后,其相对位置一经固紧后保持不变。激光投射器3由半导体激光器1和光平面发生光学***2构成,半导体激光器的激光波长为红外激光或者可见激光,光平面发生光学***2由准直光路和柱面镜构成。激光投射器可以采用商业产品,如西安华科光电有限公司生产的一字激光器等。钢轨断面成像***由CCD摄像机4和滤光***5组成,滤光***的峰值波长与激光波长相同。激光视觉传感器7固定在火车底部车体上,与被测钢轨10的距离为300~700mm。
2、调整激光投射器3,使得投射光平面8垂直于被测钢轨10并与钢轨断面平行;调整钢轨断面成像***6,保证成像***能够观测到投射光平面8与被测钢轨10表面相交形成的轨头和轨腰特征轮廓9;分别固紧激光投射器3和成像***6。
3、所说的标定靶标为一个二维平面,靶标平面上有预先设置的特征点。在靶标平面上有成矩阵排列的黑色方块,方块的数量为4~100个,方块的边长为3~50mm,其精度为0.01~0.05mm,方块间距为3~50mm,其精度为0.01~0.05mm,选取每个方块的顶点为特征点。
4、在光平面上建立二维测量坐标系om-xmym,原点om为光平面上任意点,xm轴和ym轴相互垂直且位于光平面内;
5、打开传感器的电源,启动“钢轨磨耗激光视觉动态检测***”软件,初始化图像采集卡。
6、标定激光视觉传感器的摄像机参数和传感器结构参数,标定方法和步骤参见周富强等的发明专利“一种基于平面靶标的结构光视觉传感器标定方法”,申请号:03142658.1,公开日:2005年1月19日。将标定参数保存到计算机,以备测量阶段调用。
7、如图3所示,建立标准钢轨二维设计坐标系onxnyn,以钢轨对称轴为yn轴,轨底为xn轴。AB段轮廓表示标准钢轨上轨腰部分大圆(50kg钢轨的半径为350mm,60kg钢轨的半径为400mm),称为轨腰大圆轮廓段。BC段轮廓表示轨腰部分小圆(半径为20mm),称为轨腰小圆轮廓段。EF段轮廓为轨头部分,称为轨头轮廓段。G点为侧面磨耗测量点,H点为垂直磨耗测量点。根据“铁道线路维修规则”的规定,G点位于标准断面距离钢轨顶面往下16mm处,H点位于钢轨顶面宽1/3处(距标准作用边)。选取轨腰大圆和轨腰小圆的2个圆心(用D1和D2表示)以及圆心连线D1D2与小圆的交点(用D3和D4)为对准特征点,分别计算出50kg和60kg钢轨的对准特征点的设计坐标,并保存到计算机中,以备测量阶段调用。
测量***的参数只需要标定一次,标定好测量***后,就可以进行钢轨磨耗动态检测。钢轨磨耗的具体测量步骤如下:
8、将里程表固定在火车适当位置,并通过串口连接到计算机上。
9、在火车运动或者静止的两种情况下,都可以启动测量软件的“磨耗动态检测”功能。
10、设定钢轨设计轮廓与测量轮廓对准参数。设定钢轨轮廓最小长度(由轮廓特征点数表示)阈值TL,取值范围为40~100,典型值为50左右。轨头轮廓与轨腰轮廓间距(由轮廓图像的像素坐标表示)最小阈值TD,取值范围为15~30 pixels,典型值20 pixels。设定轨腰大圆的最小长度为TBCL(由轮廓特征点数表示,保证比实际图像处理得到的长度小,典型值为70),轨腰小圆的最大长度为TSCL(由轮廓特征点数表示,典型值为30)。以上对准参数可以根据实际拍摄的钢轨图像,通过图像处理得到,设定时,比实际得到数据小5左右。设定小圆关键点距离参数T14,T14=‖D1D4‖-rBC,rBC为轨腰大圆半径。对于50kg钢轨,T14=3.51mm,对于60kg钢轨,T14=3.07mm。保存设定参数到计算机,以备后续测量调用。
11、测量程序监视串口数据状态,当接收到串口发送的启动命令和里程数据时,测量程序保存里程数据并同时控制摄像机拍摄一幅包括钢轨断面特征轮廓的图像,称为测量图像。
12、提取测量图像中钢轨断面特征轮廓的图像坐标,根据传感器测量模型,利用特征轮廓的图像坐标计算特征轮廓在测量坐标系下的二维坐标,称为特征轮廓的测量坐标。特征轮廓的图像坐标提取算法参见Carsten Steger的文章“一种对称曲线结构探测器”[An UnbiasedDetector of Curvilinear Structures],IEEE期刊《模式分析及机器智能》,20(2),第113~125页,1998年。[IEEE Transactions on PatternAnalysis and Machine Intelligence,20(2),February1998]。如图4所示,激光投射器投射光平面在钢轨表面形成的特征轮廓有两条,分别对应图4的AC段轮廓和EF段轮廓,若AC和EF轮廓的长度TL,或者A点与C点对应的图像坐标之间的距离小于TD,则表明图像中无钢轨轮廓特征,返回步骤11,重新根据串口状态,控制采集钢轨测量图像。正常情况,得到两条满足设定长度阈值以及距离阈值的轮廓特征,分别将两条轮廓的图像坐标按行坐标排序,轮廓开始点的行坐标较小的称为轨头轮廓,其长度为Nhead,另一条称为轨腰轮廓,长度为Nwaist。比较轨头轮廓的开始点和最后点的列坐标,如果开始点列坐标大于最后点,则钢轨为右钢轨,否则为左钢轨。
13、由轨腰轮廓的图像坐标,根据公式[3]和[4]表示的传感器的测量模型,计算得到轨腰轮廓的二维测量坐标。
14、根据轨腰轮廓的二维测量坐标,将轨腰轮廓自动分割成轨腰大圆和轨腰小圆两部分,提取对准特征点D1~D4的二维测量坐标。
轨腰大圆轮廓数据占整个大圆的比例不到1/3,通过无约束拟合圆,获得圆半径和中心坐标,误差大。因此,本发明假定圆半径为设计时的标准值,采用半径约束非线性优化方法拟合圆,主要的目的是获得圆的圆心坐标。非线性优化的目标函数为:
f ( x m 0 , y m 0 ) = Σ i = 1 N ( ( x mi - x m 0 ) 2 + ( y mi - y m 0 ) 2 - r ) 2 - - - [ 5 ]
(xm0,ym0)为拟合圆心坐标,r为已知半径,N为参加拟合圆的轮廓点数量。
由TBCL个轨腰轮廓点的二维测量坐标拟合轨腰大圆,得到大圆圆心D1的二维测量坐标。从轨腰的第TBCL+1个点开始,按下式计算:
d i = ( x mi - x m 0 ) 2 + ( y mi - y m 0 ) 2 - r BC   i=TBCL+1~(Nwaist-TSCL/2)     [6]
当满足(Nwaist-TSCL/2)>di>T14时,停止计算。称此时的轮廓点为小圆关键特征点,其坐标为第i个轨腰轮廓点。若未满足条件,则返回到第12步,重新根据串口状态,控制采集钢轨测量图像。
以第i-TSCL/2到第i+TSCL/2个轨腰轮廓点的二维测量坐标,采用约束半径拟合轨腰小圆,求得小圆圆心D2的二维测量坐标。由直线D1D2和小圆的方程,求得D3和D4的二维测量坐标。
15、由轨腰大圆和轨腰小圆的部分轮廓特征点二维测量坐标,按步骤14所述方法,计算测量坐标系到标准设计坐标的变换。小圆半径为20mm,大圆半径取400mm,称为60kg对准模式,大圆半径取350mm,称为50kg对准模式。分别采用60kg对准模式和50kg对准模式,计算出对准特征点D1~D4的二维测量坐标,由其中至少两个对准特征点,计算出测量坐标系到设计坐标系的二维变换。
如图5所示,测量坐标系omxmym到标准设计坐标系onxnyn之间的变换包括旋转和平移。大圆圆心D1和小圆圆心D2在标准坐标系下的坐标已知,测量坐标系下的坐标在步骤19求出,则测量坐标系中向量与设计坐标系中对应向量
Figure C200510123725D00141
之间的夹角θ表示两坐标系之间的旋转。即
θ = arctan ( y 1 m - y 2 m x 1 m - x 2 m ) - arctan ( y 1 n - y 2 n x 1 n - x 2 n ) - - - [ 7 ]
测量坐标系到设计坐标系的变换为:
x n y n = cos θ - sin θ sin θ cos θ x m y m + t x t y - - - [ 8 ]
由4对对准特征点,根据式[8],采用线性最小二乘方法求得两坐标系之间的平移。由此得到两坐标系之间的变换参数θ、tx和ty
16、根据在两种模式下求出的测量坐标系到设计坐标系的变换,分别将钢轨断面轨头部分的测量坐标转换到设计坐标系中,直接和标准钢轨比对,求出钢轨的垂直磨耗、侧面磨耗和总磨耗。
根据“铁道线路维修规则”的规定,钢轨的垂直磨耗Wv在钢轨顶面1/3处(距标准作用边)测量,侧面磨耗Wh在标准踏面(按标准断面)下16mm处测量。
总磨耗为:Ws=Wv+Wh/2。
17、比较两种对准模式下的3个磨耗值,取较小值作为最后的测量结果,并根据所采用的对准模式自动识别轨型,将测量结果、轨型以及目前所测钢轨对应的里程数保存到文件中。
18、向串口发送测量完成信号。
19、重复11~18步,进行下一次测量。
实施例
根据本发明提出测量装置和测量方法,设计了实际激光视觉传感器,对钢轨进行了现场动态测量,测量***的软件界面如图6所示。
图7给出了一幅包括50kg钢轨的测量图像及提取出的轮廓特征图像。图8为在测量坐标系的钢轨断面轮廓及变换到设计坐标系下的轮廓;图9为对准后,测量轮廓和设计轮廓。
最后测量磨耗分别为:垂直磨耗为1.9mm,侧边磨耗为-0.1mm,总磨耗为1.8mm,轨型为50R(50kg右钢轨)。在PIV3.0计算机,单次测量时间为25.5ms。
图10为火车停止时,对同一固定钢轨表面,连续进行540次测量,得到的钢轨表面磨耗重复性误差曲线。RMS误差分别为:0.10mm,0.04mm,0.10mm。

Claims (1)

1、钢轨磨耗激光视觉动态测量方法,使用钢轨磨耗激光视觉动态测量装置进行综合参数的动态测量,所说的钢轨磨耗激光视觉动态测量装置由激光投射器(3)和钢轨断面成像***(6)组成的激光视觉传感器(7)、计算机、高速图像采集卡、标定靶标和测量软件组成;激光投射器(3)由半导体激光器(1)和光平面发生光学***(2)组成;钢轨断面成像***(6)由CCD摄像机(4)和滤光***(5)组成,钢轨断面成像***(6)的光轴和激光投射器(3)的光轴交汇角为45°~60°;激光视觉传感器(7)固定在火车底部车体上,与被测钢轨(10)的距离为300~700mm;激光投射器(3)投射的光平面(8)垂直于被测钢轨(10)并与钢轨断面平行;钢轨断面成像***(6)的成像***能够观测到投射光平面(8)与被测钢轨(10)表面相交形成的轨头和轨腰特征轮廓(9);所说的标定靶标为一个二维平面,靶标平面上有预先设置的特征点,在靶标平面上有成矩阵排列的黑色方块,方块的数量为4~100个,方块的边长为3~50mm,其精度为0.01~0.05mm,方块间距为3~50mm,其精度为0.01~0.05mm,选取每个方块的顶点为特征点;其特征在于,测量的具体步骤是:
1.1、建立激光视觉传感器(7)的摄像机模型和测量模型,标定摄像机参数和传感器结构参数;
1.1.1、在光平面上建立二维测量坐标系om-xmym,原点om为光平面上任意点,xm轴和ym轴相互垂直且位于光平面内;
1.1.2、标定摄像机的内部参数和激光视觉传感器的结构参数并保存到计算机;
1.1.3、建立标准钢轨二维设计坐标系on-xnyn,以钢轨断面对称轴为yn轴,轨底为xn轴,选取标准钢轨上轨腰部分大圆和小圆的两个圆心以及圆心连线与小圆的交点为对准特征点,分别计算出对准特征点的设计坐标,并保存到计算机,以备测量阶段调用;
1.2、测量钢轨磨耗;
1.2.1、将里程表通过串口连接到计算机上;
1.2.2、设定标准轮廓与测量轮廓对准参数,设定钢轨轮廓最小长度阈值TL,它由轮廓特征点数表示,取值范围是40~100;轨头轮廓与轨腰轮廓间距最小阈值TD,它由轮廓图像的像素坐标表示,取值范围为15~30像素,典型值20像素;保存设定参数,以备后续测量调用;
1.2.3、测量程序监视串口数据状态,当接收到串口发送的启动命令和里程数据时,保存里程数据并控制摄像机拍摄一幅包括钢轨断面特征轮廓的图像,称为测量图像;
1.2.4、提取测量图像中钢轨断面特征轮廓的图像坐标,若未提取到两条满足长度大于TL的特征轮廓,或者两条特征轮廓的间距小于TD,则返回到步骤1.2.3,重新根据串口状态,控制采集钢轨测量图像;若找到两条正确的钢轨特征轮廓,则根据激光视觉传感器(7)的测量模型,利用特征轮廓的图像坐标计算钢轨断面特征轮廓的测量坐标;
1.2.5、根据钢轨断面特征轮廓的测量坐标,将测量轮廓自动分割成轨头、轨腰大圆和轨腰小圆三部分,由分割的轨腰大圆和轨腰小圆部分的特征轮廓,采用半径约束方法分别拟合轨腰大圆和轨腰小圆,小圆半径为20mm,大圆半径取400mm,称为60kg对准模式,大圆半径取350mm,称为50kg对准模式;
1.2.6、分别采用60kg对准模式和50kg对准模式,计算对准特征点的测量坐标,由其中至少2个对准特征点的测量坐标及设计坐标,计算测量坐标系到设计坐标系的二维变换;
1.2.7、根据在两种模式下求出的测量坐标系到设计坐标系的变换,分别将钢轨断面轨头部分的测量坐标转换到设计坐标系中,直接和标准钢轨比对,计算得到钢轨的垂直磨耗、侧面磨耗和总磨耗;
1.2.8、比较两种对准模式下的3个磨耗值,选取较小值作为最后的测量结果,并根据所采用的对准模式得到被测钢轨的轨型,将测量结果、轨型以及对应的里程数保存到测量结果文件中;
1.2.9、向串口发送测量完成信号;
1.2.10、重复1.2.3~1.2.9步,进行下一次测量。
CNB2005101237250A 2005-11-22 2005-11-22 钢轨磨耗激光视觉动态测量方法 Expired - Fee Related CN100462678C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005101237250A CN100462678C (zh) 2005-11-22 2005-11-22 钢轨磨耗激光视觉动态测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005101237250A CN100462678C (zh) 2005-11-22 2005-11-22 钢轨磨耗激光视觉动态测量方法

Publications (2)

Publication Number Publication Date
CN1776364A CN1776364A (zh) 2006-05-24
CN100462678C true CN100462678C (zh) 2009-02-18

Family

ID=36766004

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005101237250A Expired - Fee Related CN100462678C (zh) 2005-11-22 2005-11-22 钢轨磨耗激光视觉动态测量方法

Country Status (1)

Country Link
CN (1) CN100462678C (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109059775A (zh) * 2016-08-30 2018-12-21 大连民族大学 图像边缘提取的方法
ES2749734A1 (es) * 2018-09-21 2020-03-23 Univ Madrid Carlos Iii Dispositivo móvil de generación de imágenes 3D mediante haz láser plano

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100480627C (zh) 2007-10-26 2009-04-22 北京航空航天大学 一种钢轨磨耗综合参数车载动态测量装置及方法
CN101949116B (zh) * 2010-08-06 2012-02-08 中国科学院长春光学精密机械与物理研究所 一种采用激光修补钢轨的方法
CN102310869A (zh) * 2011-06-01 2012-01-11 王平 钢轨廓形测试小车
CN104180767B (zh) * 2013-05-24 2017-05-10 武汉铁路局武汉大型养路机械运用检修段 基于机器视觉的钢轨磨损测量装置及其打磨策略方法
CN104417583A (zh) * 2013-08-29 2015-03-18 中国铁道科学研究院铁道建筑研究所 一种实时检测钢轨铣磨车作业精度的装置
CN103630088B (zh) * 2013-11-06 2017-01-04 北京市地铁运营有限公司 基于双激光带的高精度隧道断面检测方法及装置
CN104179101A (zh) * 2014-07-22 2014-12-03 上海工程技术大学 一种精确修复地铁弯道钢轨波浪形磨损的方法
CN104359444B (zh) * 2014-11-10 2017-01-18 上海同儒机电科技有限公司 钢轨轮廓检测方法及检测***
CN104359421B (zh) * 2014-11-10 2017-05-24 上海同儒机电科技有限公司 钢轨轮廓检测***及检测方法
CN105675272A (zh) * 2014-11-20 2016-06-15 中国航空工业第六八研究所 一种判断摆式零件动态特性优略的检测方法
CN104442917A (zh) * 2014-12-31 2015-03-25 北京二七轨道交通装备有限责任公司 轨道检测方法
CN105091778B (zh) * 2015-04-28 2017-10-24 长春机械科学研究院有限公司 基于激光线结构光的异形导轨单面特征检测方法及装置
CN105172839A (zh) * 2015-09-28 2015-12-23 成都多极子科技有限公司 一种全自动火车导轨轮廓检测***
CN105423924B (zh) * 2015-11-17 2018-07-27 北京交通大学 声屏障的在线测量方法和***
CN105480256B (zh) * 2015-11-20 2018-06-22 武汉滨湖电子有限责任公司 一种高铁轨道测量装置及测量方法
CN106080661A (zh) * 2016-05-27 2016-11-09 电子科技大学 火车导轨轮廓波浪磨耗测量方法
CN108128323B (zh) * 2016-08-30 2019-09-17 大连民族大学 与钢轨磨耗量相关的激光图像特征量的选择方法
CN106643545B (zh) * 2016-10-09 2018-10-26 湖南大学 激光位移技术测量钢轨轮廓的校准方法
CN108195291B (zh) * 2018-01-03 2020-05-05 中山大学 一种基于差分光斑的运动车辆三维检测方法及检测装置
CN108917591A (zh) * 2018-04-27 2018-11-30 中国铁道科学研究院集团有限公司基础设施检测研究所 一种动态行车环境下的钢轨轮廓自动配准方法及装置
CN109017867B (zh) * 2018-08-01 2021-05-25 湖南大学 钢轨波磨动态测量方法
CN109540027A (zh) * 2018-10-11 2019-03-29 伍波 一种快速3d视觉激光测量定位方法
CN109798844B (zh) * 2018-12-28 2020-07-28 北京交通大学 一种钢轨磨耗动态测量中的轨廓自动配准方法
CN110645908A (zh) * 2019-08-16 2020-01-03 广州瑞松北斗汽车装备有限公司 车身外观缺陷检测方法及检测***
CN110926373A (zh) * 2019-12-10 2020-03-27 中南大学 铁路异物检测场景下结构光光平面标定方法及其***
CN111003018B (zh) * 2019-12-27 2021-02-26 郑州铁路职业技术学院 一种高速铁路轨道状况动态检测***及方法
CN111882599A (zh) * 2020-07-31 2020-11-03 北斗国铁(武汉)智能科技有限公司 一种钢轨磨损的数据处理方法、装置及存储介质
CN112414318B (zh) * 2020-11-10 2022-09-30 杭州申昊科技股份有限公司 基于结构光钢轨磨耗快速测量算法
CN112683195B (zh) * 2020-12-07 2022-10-25 中国铁道科学研究院集团有限公司基础设施检测研究所 钢轨纵向标定方法及装置
CN112902869B (zh) * 2021-01-21 2022-10-25 中国铁道科学研究院集团有限公司基础设施检测研究所 轨廓测量***激光平面调整方法及装置
CN113432630B (zh) * 2021-05-13 2023-02-17 国网河北省电力有限公司衡水供电分公司 分布式光纤传感监测***
CN114046744A (zh) * 2021-11-19 2022-02-15 武汉利德测控技术有限公司 一种2d钢轨轮廓智能检测***

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001304824A (ja) * 2000-04-19 2001-10-31 Tokimec Inc 携帯式レール断面測定装置
CN1566906A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种基于平面靶标的结构光视觉传感器标定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001304824A (ja) * 2000-04-19 2001-10-31 Tokimec Inc 携帯式レール断面測定装置
CN1566906A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种基于平面靶标的结构光视觉传感器标定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
An Unbiased Detector of Curvilinear Structures. Carsten Strger.IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE,Vol.20 No.2. 1998
An Unbiased Detector of Curvilinear Structures. Carsten Strger.IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE,Vol.20 No.2. 1998 *
实时动态检测钢轨磨耗***研究. 纪淑波,蒋本和,靳文瑞.光电技术应用,第19卷第4期. 2004
实时动态检测钢轨磨耗***研究. 纪淑波,蒋本和,靳文瑞.光电技术应用,第19卷第4期. 2004 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109059775A (zh) * 2016-08-30 2018-12-21 大连民族大学 图像边缘提取的方法
CN109059775B (zh) * 2016-08-30 2020-11-24 大连民族大学 具有图像边缘提取步骤的钢轨磨耗检测方法
ES2749734A1 (es) * 2018-09-21 2020-03-23 Univ Madrid Carlos Iii Dispositivo móvil de generación de imágenes 3D mediante haz láser plano
WO2020058557A1 (es) * 2018-09-21 2020-03-26 Universidad Carlos Iii De Madrid Dispositivo móvil de generación de imágenes 3d mediante haz láser plano

Also Published As

Publication number Publication date
CN1776364A (zh) 2006-05-24

Similar Documents

Publication Publication Date Title
CN100462678C (zh) 钢轨磨耗激光视觉动态测量方法
CN101532827B (zh) 一种激光视觉钢轨磨耗测量的纠偏方法
US6909514B2 (en) Wheel profile inspection apparatus and method
CN100480627C (zh) 一种钢轨磨耗综合参数车载动态测量装置及方法
CN107607044B (zh) 一种基于激光位移传感器的车轮踏面磨耗检测方法
Sun et al. Motion deviation rectifying method of dynamically measuring rail wear based on multi-line structured-light vision
CN102749061B (zh) 基于动态模板的钢轨磨耗测量方法
CN102901457B (zh) 一种列车车轮直径动态测量方法及***
CN107421502B (zh) 一种铁路货车限界自动化测量方法
CN101957182B (zh) 一种大口径高陡度光学镜面在线测量***
CN105292180A (zh) 基于多种传感器的非接触式轮对尺寸在线检测方法和装置
CN105133471B (zh) 一种基于线结构光路面检测***的路面深度图像生产方法
CN103499302A (zh) 基于结构光视觉成像***的凸轮轴直径尺寸在线测量方法
CN101694084A (zh) 地面车载移动检测***
CN109017867B (zh) 钢轨波磨动态测量方法
CN100999220A (zh) 车辆轮对直径在线检测方法及装置
CN103075976A (zh) 一种高速列车动态包络线测量方法
CN102564335A (zh) 一种大型隧道变形测量方法
CN110375722A (zh) 一种基于点云数据提取盾构隧道管片张开量的方法
CN104648440B (zh) 一种轮对几何参数在线测量***及其方法
CN102997863A (zh) 一种全口径光学非球面面形误差直接检测***
CN109443245A (zh) 一种基于单应性矩阵的多线结构光视觉测量方法
CN105292181A (zh) 一种基于两种传感器的轮对尺寸在线检测方法及装置
CN102944188A (zh) 一种点扫描三维形貌测量***标定方法
CN110243293A (zh) 基于结构光和机器视觉的管片错台快速检测装置与方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090218

Termination date: 20121122