CN104296682A - Contour registration method used in train guide rail contour measurement based on machine vision - Google Patents

Contour registration method used in train guide rail contour measurement based on machine vision Download PDF

Info

Publication number
CN104296682A
CN104296682A CN201410588757.7A CN201410588757A CN104296682A CN 104296682 A CN104296682 A CN 104296682A CN 201410588757 A CN201410588757 A CN 201410588757A CN 104296682 A CN104296682 A CN 104296682A
Authority
CN
China
Prior art keywords
image
profile
contour line
rail
line
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.)
Withdrawn
Application number
CN201410588757.7A
Other languages
Chinese (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201410588757.7A priority Critical patent/CN104296682A/en
Publication of CN104296682A publication Critical patent/CN104296682A/en
Priority to CN201510037287.XA priority patent/CN104655043B/en
Withdrawn legal-status Critical Current

Links

Abstract

The invention discloses a contour registration method used in train guide rail contour measurement based on machine vision. Contour measurement comprises the first step of collecting steel rail images, the second step of obtaining a single-pixel curved line of a steel rail profile through computer processing, the third step of carrying out calibration restoring and the fourth step of carrying out contour registration. The contour registration method comprises the steps of enabling an actual contour line and a standard contour line to be unified in a coordinate system and utilizing asymptotic lines of the two contour lines to obtain an included angle between the actual contour line and the standard contour line, rotating the actual contour line and enabling the asymptotic lines of the two contour lines to be paralleled, and horizontally moving the actual contour line and enabling ends of the actual contour line to be coincident with ends of the standard contour line. According to the contour registration method, the actual contour line and the standard contour line are unified into the coordinate system and the actual contour line is processed, the asymptotic lines of the actual contour line and the standard contour line are paralleled and coincide at ends, and therefore the registration result is accurate, observation and calculation are carried out conveniently and reliable abrasion loss data can be achieved.

Description

A kind of based on profile method for registering in the measurement of machine vision train guide rail profile
Technical field
The invention belongs to the contactless automatic fields of measurement based on image, relate to image procossing, laser and optics, computing machine, automation field, specifically refer to a kind of based on profile method for registering in the measurement of machine vision train guide rail profile.
Background technology
Rail, as the important carrier in transportation by railroad, is one of most base unit in transportation by railroad, and its state quality directly affects the security of operation of train.Especially, when train runs under high speed, heavy duty, high density, rail wear strengthens, and the lost of life, causes high risks to train safe.The Ministry of Railways is in rail track " dimension rule ", and, severe injury light for rail defines the limit value of rail head abrasion, and requires will carry out in time processing or changing, to ensure the security of operation of train after rail head of rail abrasion exceed limit value when detecting.
Along with the lifting of transportation by railroad rate of traffic flow, being increased sharply of travelling speed, the abrasion of rail are more and more serious, especially the sharp radius curve Innerside wear at rail head place is particularly serious, have impact on the even running of train, the geometrical deviation of small track all likely causes train (especially high ferro) rocks and various vibration, impact is had, even because path wear increases the weight of to cause the skew laterally run on the comfort level of train driving and traveller's.When taking turns transverse force and increasing, just likely can there is the situations such as derail, bring great harm to passenger and traffic safety.The degree of wear of rail determines the serviceable life of rail, when attrition value arrives certain limit, needs to change in time.The measurement data of rail wear is that work business system determines that rail is the need of the important evidence changed, and therefore, detects the abrasion value size of track exactly, and to analyze each section abrasion situation be a very important job.
Mechanical Contact mensuration is that China makes a kind of method with the most use at present, slide calliper rule are mainly adopted to carry out hand dipping, this method is rely on manually completely, the shortcoming of contact type measurement is that speed is too slow, error is large, human factor impact is large, efficiency is low, and is not easy to ex-post analysis statistics and files.Another branch of contact type measurement is exactly the rail wear detector of machinery and electronics, the technological synthesis application and development such as able to programme.Also have photoelectric coding method, displacement sensing method and the non-contact type measuring method based on image procossing in addition, this method obtains the contour curve of rail profile indirectly by sensor, and the curve of drafting and the curve comparison of standard are obtained wearing valve.Photoelectric encoder and displacement transducer can be divided into again according to the difference of sensor.The MINIPROF system of greenwood company of Denmark development, it is made up of a notebook, special expansion box and wheel track measurement mechanism.Contact tested Rail Surface by the little magnet-wheel on measurement mechanism, the little magnet-wheel of manual control, obtains the polar data of sensor by computing machine, calculate profile after being converted to Cartesian coordinates.The track of magnet-wheel core can go on record, and computing machine just can calculate the actual profile of section.Domestic photoelectric coding checkout equipment has the WRS2000 shaped steel rail rail head appearance instrument for measuring and recording by railway institute of Tongji University and Jingjiang Hua Xing company cooperation research and development, and it is realized by high-precision photoelectric encoder.Tongji University still further developed a kind of wave abrasion measuring device of low-arranged sliding type railway steel rail in addition, this device adopts simple mechanical component to form longitudinal measuring basis and guide rail, and adopt USB interface and control technology and computer radio transmission measurement result, can immediately result be presented on screen.
A kind of noncontact measurement based on machine vision profilometry, have speed fast, the advantages such as profile overall picture, data management be convenient can be obtained.In recent years due to the high speed development of Digital image technology, laser technology and computer technology, the train guide rail profile measuring technique based on machine vision is greatly developed.A kind of noncontact measurement based on machine vision profilometry, have speed fast, the advantages such as profile overall picture, data management be convenient can be obtained.
Measure based on machine vision train guide rail profile, wherein the link of most critical is by obtained guide rail profile and real profile comparison, and could calculate and carry out abrasion calculating, the process of this comparison is called profile registration.Profile registration not pair time, the abrasion size calculated has very large, full of prunes deviation, and only have profile registration correct, the abrasion size calculated is only trust data within the scope of certain error.
Summary of the invention
The object of the invention is to overcome the deficiencies in the prior art, profile method for registering in providing a kind of contactless train guide rail profile to measure, by treated guide rail real contour line and the unification of standard wheels profile in same coordinate system, and real contour line is rotated, the process such as translation, make real contour line parallel with the asymptotic line of standard wheels profile and end points coincidence, the relative position of two outline lines is determined with this, registration result is accurate, clear and easy to understand, observe and convenience of calculation, can obtain reliable abrasion data.
The object of the invention is to be achieved through the following technical solutions: a kind of measure based on machine vision train guide rail profile in profile method for registering, wherein, what profile measurement adopted is non-contact measurement, and it can improve the work efficiency that on-the-spot railway detects department greatly.The cardinal principle measured obtains rail profile contour curve by laser light cutting method, simultaneously, CCD camera gathers laser image, and image information is transferred in the process chip of the computing machine be loaded on dolly, through measurement of rail wear based program, a series of process is carried out to the image collected, obtain single pixel curve of rail profile.By the matching and correlation to curve image, detect with railway the rail wear point that department defines and compare, draw the value of corresponding side wear and vertical abrasion.
Profile measurement specifically comprises the following steps: S1, gathers rail image; S2, processes the rail image transmitting collected to computing machine, and obtains single pixel curve of rail profile; S3, carries out demarcation reduction to single pixel curve, obtains real contour line; S4, carries out profile registration;
Profile registration in described step S4 comprises following sub-step:
S41: using the most end end points on guide rail inoperative limit as the scaling point of outline line, i.e. initial point, set up cartesian coordinate system, by the real contour line of guide rail and the unification of standard wheels profile in cartesian coordinate system, get the asymptotic line of two outline lines respectively, and utilize the asymptotic line of two outline lines to obtain angle theta between the two;
S42: rotate real contour line with twiddle factor R (θ), makes the asymptotic line of two outline lines be parallel to each other, obtains postrotational real contour line;
S43: postrotational real contour line is carried out translation, makes this real profile line endpoints overlap with the end points of standard wheels profile, completes the registration of outline line, and now, the real contour line relative position after standard wheels profile and registration is determined.
Further, the method gathering rail image in described step S1 is: obtain rail profile contour curve by laser light cutting method, and adopt CCD camera to gather laser image, CCD camera collection be the scattered light of profiled outline, the scattered light of this profile is two dimensional image model.
Further, demarcation in step S3 reduction refers to the two dimensional image model for collecting, and by calculating the information of actual three dimensions epigraph of reducing, is reconstructed the process of the corresponding relation between two dimensional image and three dimensional practicality by geometric model.Described step S3 specifically comprises following sub-step:
S31: by the change of camera imaging and pattern distortion, if its hidden parameter relational expression is:
x = au 2 + buv + cv 2 + du + ev + f y = gu 2 + huv + iv 2 + ju + kv + l - - - ( 1 )
Wherein, if P (x, y) represents is world coordinate system, P (u, what v) represent is image space coordinate system, parameter a, b, c, d, e, f, g, h, i, j, k, l is undetermined coefficient, and the object of demarcation is exactly will obtain stealthy undetermined parameter by the relation of fixing world coordinate system and image space coordinate, then obtains corresponding world coordinates by the parameter of stealth and image space coordinate points in actual applications and fastens object coordinates point;
S32: in order to eliminate the impact of indivedual points on result, adopt least square method to ask for undetermined parameter: the optimum matching being asked for function parameter by the quadratic sum asking for least error, adopts 20 groups of corresponding coordinates to ask for stealthy parameter;
S33: adopt semi-automation to determine calibration method, by presetting the coordinate points of the point in calibration plate, i.e. object coordinates, gathering image space coordinate figure corresponding in image space plane by the method for image procossing, formula (1) being separated into the form of matrix:
C=[a b c d e f g h i j k l] T (4)
C=A -1B (5)
Wherein, formula (2) is object coordinates system matrix, and formula (3) is corresponding image space coordinate points matrix, and formula (4) is matrix of coefficients, and through type (5) solves matrix of coefficients totally 12 parameters;
S34: in measuring system, is transformed into formula (5):
A=BC -1 (6)
Test matrix of coefficients C is known, then carrys out by the coordinate B matrix of image mid point the object coordinates that reverse goes out correspondence, the namely real contour line of rail profile section;
S35: in the process of image calibration reduction, because the image gathered there will be geometry deformation, so need to carry out geometric transformation to the image obtained, geometric transformation is in another image by the pixel-map in image, comprise the mapping of horizontal direction and the mapping of vertical direction, geometric transformation comprises the content of two aspects: pixel coordinate conversion and brightness interpolating:
S351: pixel coordinate conversion is mapping between output image and input picture, the plane of delineation (x, y) is be through geometric transformation to converted relational expression to the plane of delineation (X, Y):
X = Σ r = 0 m Σ k = 0 m - r a rk x r y k Y = Σ r = 0 m Σ k = 0 m - r b rk x r y k - - - ( 7 )
Like this, by the coefficient a of conversion rkand b rkjust can determine the relation of geometric transformation, whole coordinate system may be changed when image procossing is carried out in geometric transformation, be determined the information of coordinate system transformation by Jacobi coefficient:
J = | ∂ ( X , Y ) ∂ ( x , y ) | = | ∂ X ∂ x ∂ X ∂ y ∂ Y ∂ x ∂ Y ∂ y | - - - ( 8 )
As J=0, conversion is irreversible, work as J=1, represent that image area has unchangeability when converting, when measuring system gathers image, the most incident image geometry distortion is exactly the rotation of image, in the geometric transformation of image, rotating R is the most frequently used one conversion, and transformation relation is such as formula (9):
So, as long as obtain the angle θ rotated, just by formula (9), image is carried out rotational transform;
S352: the object of brightness interpolating is the gray-scale value of the pixel determined after conversion, the new picture point obtained after conversion may be a non-integral coordinate, need to process to the pixel value of image the gray scale obtaining this point, it has direct impact to the quality of image after geometric transformation, and the method for interpolation comprises arest neighbors, linear, bicubic three kinds.
Particularly, described arest neighbors Outlet approach concrete operations are: establish the result of gray-level interpolation to be expressed as Convolution Formula:
f n ( X , Y ) = Σ l = - ∞ ∞ Σ k = - ∞ ∞ g s ( lΔX , kΔY ) h n ( X - lΔX , Y - kΔY ) - - - ( 10 )
Wherein, h nwhat represent is interpolation kernel, represents very little field;
Arest neighbors interpolation is that its interpolation formula is expressed as by point (X, Y) with the gray-scale value put to replace this from the gray-scale value of its nearest point:
f 1(X,Y)=g s[round(x),round(y)] (11)。
Described linear interpolation considers four consecutive point of point (X, Y), and the gray-scale value of these four points is represented the gray-scale value of this point, interpolation formula is:
f 2(X,Y)=(1-a)(1-b)g s(l,k)+a(1-b)g s(l+1,k)
+b(1-a)g s(l,k+1)+abg s(l+1,k+1) (12)。
The invention has the beneficial effects as follows:
1, the profile method for registering in providing a kind of contactless train guide rail profile to measure, by treated guide rail real contour line and the unification of standard wheels profile in same coordinate system, and real contour line is rotated, the process such as translation, make real contour line parallel with the asymptotic line of standard wheels profile and end points coincidence, the relative position of two outline lines is determined with this, registration result is accurate, clear and easy to understand, observe and convenience of calculation, then the calculating of relevant abrasion is carried out in the abrasion definition specified according to the Ministry of Railways, can obtain reliable abrasion data;
2, adopt contactless measurement, adopt laser light cutting method to obtain rail profile contour curve, and adopt CCD camera to gather laser image, and the image transmitting collected is processed to computing machine; There is image acquisition speed fast, the advantage such as profile overall picture and data management convenience can be obtained;
3, by rail profile line that laser line generator and CCD camera collect, outline line relatively is reliably obtained through the image processing means such as noise processed, connection of broken lines, and then suitable demarcation reduction treatment is carried out to outline line, two dimensional image model CCD camera collected is by calculating the information of actual three dimensions epigraph of reducing, can by elimination of the distortion, rediscover outline line, makes result of calculation more accurately and reliably.
Accompanying drawing explanation
Fig. 1 is profile measurement schematic diagram of the present invention;
Fig. 2 is rail profile skeleton pattern figure of the present invention;
Fig. 3 is outline line registration process schematic diagram of the present invention;
Fig. 4 is the square trrellis diagram of standard;
Fig. 5 is the image formed after two kinds typically distortion;
Fig. 6 is the image rotation schematic diagram due to guide rail, the generation of CCD attitudes vibration of the present invention;
Fig. 7 is linear gradation interpolation schematic diagram.
Embodiment
Further illustrate technical scheme of the present invention below in conjunction with accompanying drawing, but the content that the present invention protects is not limited to the following stated.
Profile measurement principle as shown in Figure 1, employing be non-contact measurement, it greatly can improve on-the-spot railway and detect the work efficiency of department.The cardinal principle measured obtains rail profile contour curve by laser light cutting method, adopt semiconductor county laser instrument Emission Lasers, simultaneously, CCD camera gathers laser image, and image information is transferred in the process chip of the computing machine be loaded on dolly, through measurement of rail wear based program, a series of process is carried out to the image collected, obtain single pixel curve of rail profile.By the matching and correlation to curve image, detect with railway the rail wear point that department defines and compare, draw the value of corresponding side wear and vertical abrasion.
In measurement of rail wear based process, detect department and fix, have corresponding definition to rail head of rail Detection of content according to rail, main measuring content comprises side wear, vertical abrasion and always wears away, and Figure 2 shows that rail wear section and calibration schematic diagram.In figure, what W1 represented is vertical abrasion, and its definition is the thickness of 2/3 place's vertical direction abrasion of the rail head width on distance rail inoperative limit; What W2 represented is side wear, and its definition is that the definition of total abrasion is that vertical wearing valve adds 1/2 times of horizontal wearing valve from the rail head top down width that wears away in the horizontal direction of 16mm.
The rail profile line collected by laser line generator and ccd array, is obtained outline line relatively reliably through image processing meanses such as noise processed (filtering etc.), connection of broken lines, and then carries out suitable calibration process to outline line.But the plane formed due to laser line generator is not orthogonal with contoured surface, therefore the outline line obtained not is real outline line, that is there is distortion, so need elimination of the distortion, with rediscover outline line.
A kind of based on profile method for registering in the measurement of machine vision train guide rail profile, wherein, profile measurement specifically comprises the following steps: S1, gathers rail image; S2, processes the rail image transmitting collected to computing machine, and obtains single pixel curve of rail profile; S3, carries out demarcation reduction to single pixel curve, obtains real contour line; S4, carries out profile registration;
Profile registration in described step S4 comprises following sub-step:
S41: as shown in Fig. 3 (a), using the most end end points P on guide rail inoperative limit as the scaling point of outline line, i.e. initial point, set up cartesian coordinate system, by the real contour line of guide rail and the unification of standard wheels profile in cartesian coordinate system, get the asymptotic line of two outline lines respectively, and utilize the asymptotic line of two outline lines to obtain angle theta between the two;
S42: rotate the real contour line in Fig. 3 (a) with twiddle factor R (θ), makes the asymptotic line of two outline lines be parallel to each other, obtains postrotational real contour line;
S43: the postrotational real contour line in Fig. 3 (b) is carried out translation, this real profile line endpoints is overlapped with the end points of standard wheels profile, as shown in Fig. 3 (c), complete the registration of outline line, now, the real contour line relative position after standard wheels profile and registration is determined.
Further, the method gathering rail image in described step S1 is: obtain rail profile contour curve by laser light cutting method, and adopt CCD camera to gather laser image, CCD camera collection be the scattered light of profiled outline, the scattered light of this profile is two dimensional image model.
When obtaining rail profile image based on light cross-section method, due to CCD camera collection is the scattered light of profiled outline, CCD camera imaging plane, laser plane are at an angle, based in the non-contact measurement system of machine vision, all can be designed into the demarcation reduction of image.This is because, for the two dimensional image model that CCD camera collects, need, by calculating the information of actual three dimensions epigraph of reducing, to reconstruct the corresponding relation between two dimensional image and three dimensional practicality by the geometric model of CCD imaging.The process obtaining the parameter of this geometric model is exactly the process that CCD camera is demarcated.
Ideally, video camera imaging is a kind of linear pin-hole imaging model.But in practice, the camera lens that video camera is installed or lens group, due to the reason of processing or other some reasons, there will be the phenomenon of distortion in the place away from optical axis center during imaging.Be illustrated in figure 4 the square trrellis diagram of standard, shown in Fig. 5 (a) He Fig. 5 (b) is the image that Fig. 4 is formed after two kinds typically distortion.After CCD imaging, the distortion type of appearance is primarily of two classes: the barrel distortion shown in the pincushion type distortion shown in Fig. 5 (a) and Fig. 5 (b).Now, camera imaging model is a kind of nonlinear, so the process of image rectification also should comprise the impact of distortion correction on image.
Camera plane is relative with the position of laser plane is fixing, so, only need surely show once for the mechanical hook-up fixed.But when actual measurement, because dolly shake and other unknown causes are on the impact of machinery, in order to ensure the accuracy measured, still need regularly to demarcate camera.Demarcation reduction in described step S3 specifically comprises following sub-step:
S31: the model adopting hidden parameter to calibrate reacts the corresponding relation between two dimensional image coordinate and three dimensional practicality, these parameters calculated do not have concrete physical significance, just characterize corresponding variation relation; By the change of camera imaging and pattern distortion, if its hidden parameter relational expression is:
x = au 2 + buv + cv 2 + du + ev + f y = gu 2 + huv + iv 2 + ju + kv + l - - - ( 1 )
Wherein, if P (x, y) represents is world coordinate system, P (u, what v) represent is image space coordinate system, parameter a, b, c, d, e, f, g, h, i, j, k, l is undetermined coefficient, and the object of demarcation is exactly will obtain stealthy undetermined parameter by the relation of fixing world coordinate system and image space coordinate, then obtains corresponding world coordinates by the parameter of stealth and image space coordinate points in actual applications and fastens object coordinates point;
S32: the stealthy parameter asked in theory in formula (1) needs 6 groups of corresponding image space points and object space point, but in order to eliminate the impact of indivedual points on result, adopt least square method to ask for undetermined parameter: the optimum matching being asked for function parameter by the quadratic sum asking for least error, adopt 20 groups of corresponding coordinates to ask for stealthy parameter;
S33: adopt semi-automation to determine calibration method, by presetting the coordinate points of the point in calibration plate, i.e. object coordinates, gathering image space coordinate figure corresponding in image space plane by the method for image procossing, formula (1) being separated into the form of matrix:
C=[a b c d e f g h i j k l] T (4)
C=A -1B (5)
Wherein, formula (2) is object coordinates system matrix, and formula (3) is corresponding image space coordinate points matrix, and formula (4) is matrix of coefficients, and through type (5) solves matrix of coefficients totally 12 parameters;
S34: in measuring system, is transformed into formula (5):
A=BC -1 (6)
Test matrix of coefficients C is known, then carrys out by the coordinate B matrix of image mid point the object coordinates that reverse goes out correspondence, the namely real contour line of rail profile section;
S35: in the process of image calibration reduction, the visual position of laser plane and camera is fixing, and the scaling parameter obtained by formula (5) records in the lab, and during calibration, dolly lies on the rail level of rail head of rail.In actual measurement, the position of dolly can not lie on rail level completely, inevitably has inclination, shake in the process of promotion.Illustrate for the some sections on rail, when dolly tilts, just in the section surface thereof perpendicular to Rail Surface, now the position of laser plane still matches with the section of this point, do not change, but the position due to camera is the bottom being fixed on dolly, the rotation of dolly has just driven the rotation of camera position, and the position be equivalent to as plane and object plane there occurs small change.Now, the stealthy parameter that through type (1) obtains is to when having the rail profile curve of certain inclination to calibrate, and the image obtained just has certain inclination.When graphical analysis, because the image gathered there will be geometry deformation, so need to carry out geometric transformation to the image obtained, the impact brought to graphical analysis could be out of shape by removal of images like this.Geometric transformation be by the pixel-map in image in another image, comprise the mapping of horizontal direction and the mapping of vertical direction.When carrying out geometric transformation, mainly comprise two aspects, i.e. pixel coordinate conversion and brightness interpolating:
S351: pixel coordinate conversion is mapping between output image and input picture, the plane of delineation (x, y) is be through geometric transformation to converted relational expression to the plane of delineation (X, Y):
X = Σ r = 0 m Σ k = 0 m - r a rk x r y k Y = Σ r = 0 m Σ k = 0 m - r b rk x r y k - - - ( 7 )
Like this, by the coefficient a of conversion rkand b rkjust can determine the relation of geometric transformation, whole coordinate system may be changed when image procossing is carried out in geometric transformation, be determined the information of coordinate system transformation by Jacobi coefficient:
J = | ∂ ( X , Y ) ∂ ( x , y ) | = | ∂ X ∂ x ∂ X ∂ y ∂ Y ∂ x ∂ Y ∂ y | - - - ( 8 )
As J=0, conversion is irreversible, works as J=1, represent that image area has unchangeability when converting, when measuring system gathers image, the most incident image geometry distortion is exactly the rotation of image, and Fig. 6 is the image rotation schematic diagram for producing due to guide rail, CCD attitudes vibration; Solid line shown in figure is the actual curve obtained, the rail image collected when dotted line is and keeps flat, and the degree of tilt of two curves is θ.In the geometric transformation of image, rotating R is the conversion of the most frequently used one, transformation relation such as formula for:
So, as long as obtain the angle θ rotated, just by formula (9), image is carried out rotational transform;
S352: the object of brightness interpolating is the gray-scale value of the pixel determined after conversion, the new picture point obtained after conversion may be a non-integral coordinate, need to process to the pixel value of image the gray scale obtaining this point, it has direct impact to the quality of image after geometric transformation, and the method for interpolation comprises arest neighbors, linear, bicubic three kinds.
Particularly, described arest neighbors Outlet approach concrete operations are: establish the result of gray-level interpolation to be expressed as Convolution Formula:
f n ( X , Y ) = Σ l = - ∞ ∞ Σ k = - ∞ ∞ g s ( lΔX , kΔY ) h n ( X - lΔX , Y - kΔY ) - - - ( 10 )
Wherein, h nwhat represent is interpolation kernel, represents very little field;
Arest neighbors interpolation is that its interpolation formula is expressed as by point (X, Y) with the gray-scale value put to replace this from the gray-scale value of its nearest point:
f 1(X,Y)=g s[round(x),round(y)] (11)。
Particularly, described linear interpolation considers four consecutive point of point (X, Y), and the gray-scale value of these four points is represented the gray-scale value of this point, interpolation formula is:
f 2(X,Y)=(1-a)(1-b)g s(l,k)+a(1-b)g s(l+1,k)
+b(1-a)g s(l,k+1)+abg s(l+1,k+1) (12)。
Fig. 7 is linear interpolation schematic diagram, and wherein solid line represents input picture, and dotted line represents the image grid that output image obtains through inverse transformation.
The geometric transformation problem of image when solving measurement by pixel transform and brightness interpolating.But, above rotation is only applicable to the low-angle variation that dolly rotates, here set the maximum anglec of rotation is by experiment 5 °, when the angle rotated is excessive, owing to also needing the factor considered to be the aberration problems of image during image restoring, between three dimensional practicality image to two-dimensional image plane, be transformed to nonlinear relation.At this moment, rotate the attitudes vibration caused just there will be larger error by reacting dolly as the change of plane.So when designing dolly, having a scroll wheel be fixed and clamped to the position of dolly on rail, allowing little truck position lie in rail rail level as far as possible.In the process promoted, we also will allow dolly keep equilibrium position.Like this, the impact of nonlinear transformation on result can be reduced.
Technique effect of the present invention is further illustrated below in conjunction with experiment.
When certain rail welding plant is tested, the rail type of experiment is 60 rails.The flow process of whole measuring system according to setting is gathered rail profile curve and preserve data.By automatic and manual mode image data respectively, system can normally be run.Following table 1 is depicted as the contrast of system for field measurement data and kind of calliper data, reliability and stability when data carry out evaluating system field survey by experiment.
The reliability and stability of table 1 field survey verification msg
Table 1 (a) field survey different measuring point experimental data
The same measurement point experimental data of table 1 (b) field survey
Standard deviation in table is the dispersion degree of data deviation average, and above two list datas can be found out, system still has good Stability and veracity for in-site measurement.
When experienced Chongqing line is tested, it is very serious that guide rail has all worn away, and the rail head of rail has back edge phenomenon, and it is very serious that rail limit to be measured all wears away.Following table 2 is depicted as the one piece of data wherein tested, and by obtaining the information of rail wear after data analysis software process, and is printed by printing reports.
Table 2 field survey wears away comparatively heavy rail experimental data
In form, what first print is the essential information measured.First form is the abrasion statistics exceeded standard, the measurement data of a certain section that what table 2 represented is.Wherein, the font that have passed different colours for the point beyond slight wound or severely injured standard marks.In in-site measurement, the sampling interval of setting is 3m, and the spacing that dolly is passed by is 250m.System is for manually and automatically sampling and can normally work.A contrast has roughly been passed through with the carrying out of kind of calliper (in measuring at the scene for the data collected, kind of calliper also has some enchancement factors and occurs), wherein vertically wear away error roughly at ± 0.5mm, level abrasion error is ± 1.0mm.
Experiment determines rationality and the correctness of matching process of the present invention, and this measuring method can stable operation obtain good measurement effect, meets the measurement requirement that the Ministry of Railways specifies.
Those of ordinary skill in the art will appreciate that, embodiment described here is to help reader understanding's principle of the present invention, should be understood to that protection scope of the present invention is not limited to so special statement and embodiment.Those of ordinary skill in the art can make various other various concrete distortion and combination of not departing from essence of the present invention according to these technology enlightenment disclosed by the invention, and these distortion and combination are still in protection scope of the present invention.

Claims (6)

1. a profile method for registering in measuring based on machine vision train guide rail profile, it is characterized in that, profile measurement comprises the following steps: S1, gathers rail image; S2, processes the rail image transmitting collected to computing machine, and obtains single pixel curve of rail profile; S3, carries out demarcation reduction to single pixel curve, obtains real contour line; S4, carries out profile registration;
Profile registration in described step S4 comprises following sub-step:
S41: using the most end end points on guide rail inoperative limit as the scaling point of outline line, i.e. initial point, set up cartesian coordinate system, by the real contour line of guide rail and the unification of standard wheels profile in cartesian coordinate system, get the asymptotic line of two outline lines respectively, and utilize the asymptotic line of two outline lines to obtain angle theta between the two;
S42: rotate real contour line with twiddle factor R (θ), makes the asymptotic line of two outline lines be parallel to each other, obtains postrotational real contour line;
S43: postrotational real contour line is carried out translation, makes this real profile line endpoints overlap with the end points of standard wheels profile, completes the registration of outline line, and now, the real contour line relative position after standard wheels profile and registration is determined.
2. profile method for registering according to claim 1, it is characterized in that, the method gathering rail image in described step S1 is: obtain rail profile contour curve by laser light cutting method, and adopt CCD camera to gather laser image, CCD camera collection be the scattered light of profiled outline, the scattered light of this profile is two dimensional image model.
3. profile method for registering according to claim 1, it is characterized in that, demarcation reduction in described step S3 refers to the two dimensional image model for collecting, by calculating the information of actual three dimensions epigraph of reducing, reconstructed the process of the corresponding relation between two dimensional image and three dimensional practicality by geometric model.
4. profile method for registering according to claim 3, is characterized in that, described step S3 specifically comprises following sub-step:
S31: by the change of camera imaging and pattern distortion, if its hidden parameter relational expression is:
Wherein, if P (x, y) represents is world coordinate system, P (u, what v) represent is image space coordinate system, parameter a, b, c, d, e, f, g, h, i, j, k, l is undetermined coefficient, and the object of demarcation is exactly will obtain stealthy undetermined parameter by the relation of fixing world coordinate system and image space coordinate, then obtains corresponding world coordinates by the parameter of stealth and image space coordinate points in actual applications and fastens object coordinates point;
S32: in order to eliminate the impact of indivedual points on result, adopt least square method to ask for undetermined parameter: the optimum matching being asked for function parameter by the quadratic sum asking for least error, adopts 20 groups of corresponding coordinates to ask for stealthy parameter;
S33: adopt semi-automation to determine calibration method, by presetting the coordinate points of the point in calibration plate, i.e. object coordinates, gathering image space coordinate figure corresponding in image space plane by the method for image procossing, formula (1) being separated into the form of matrix:
B=[x 1 x 2....x 20 y 1 y 2.....y 20] T (3)
C=[a b c d e f g h i j k l] T (4)
C=A -1B (5)
Wherein, formula (2) is object coordinates system matrix, and formula (3) is corresponding image space coordinate points matrix, and formula (4) is matrix of coefficients, and through type (5) solves matrix of coefficients totally 12 parameters;
S34: in measuring system, is transformed into formula (5):
A=BC -1 (6)
Test matrix of coefficients C is known, then carrys out by the coordinate B matrix of image mid point the object coordinates that reverse goes out correspondence, the namely real contour line of rail profile section;
S35: in the process of image calibration reduction, because the image gathered there will be geometry deformation, so need to carry out geometric transformation to the image obtained, geometric transformation is in another image by the pixel-map in image, comprise the mapping of horizontal direction and the mapping of vertical direction, geometric transformation comprises the content of two aspects: pixel coordinate conversion and brightness interpolating:
S351: pixel coordinate conversion is mapping between output image and input picture, the plane of delineation (x, y) is be through geometric transformation to converted relational expression to the plane of delineation (X, Y):
Like this, by the coefficient a of conversion rkand b rkjust can determine the relation of geometric transformation, whole coordinate system may be changed when image procossing is carried out in geometric transformation, be determined the information of coordinate system transformation by Jacobi coefficient:
As J=0, conversion is irreversible, work as J=1, represent that image area has unchangeability when converting, when measuring system gathers image, the most incident image geometry distortion is exactly the rotation of image, in the geometric transformation of image, rotating R is the most frequently used one conversion, and transformation relation is such as formula (9):
So, as long as obtain the angle θ rotated, just by formula (9), image is carried out rotational transform;
S352: the object of brightness interpolating is the gray-scale value of the pixel determined after conversion, the new picture point obtained after conversion may be a non-integral coordinate, need to process to the pixel value of image the gray scale obtaining this point, it has direct impact to the quality of image after geometric transformation, and the method for interpolation comprises arest neighbors, linear, bicubic three kinds.
5. profile method for registering according to claim 4, is characterized in that, described arest neighbors Outlet approach concrete operations are: establish the result of gray-level interpolation to be expressed as Convolution Formula:
Wherein, h nwhat represent is interpolation kernel, represents very little field;
Arest neighbors interpolation is that its interpolation formula is expressed as by point (X, Y) with the gray-scale value put to replace this from the gray-scale value of its nearest point:
f 1(X,Y)=g s[round(x),round(y)] (11)。
6. profile method for registering according to claim 4, is characterized in that, described linear interpolation considers four consecutive point of point (X, Y), and the gray-scale value of these four points is represented the gray-scale value of this point, interpolation formula is:
f 2(X,Y)=(1-a)(1-b)g s(l,k)+a(1-b)g s(l+1,k)
+b(1-a)g s(l,k+1)+abg s(l+1,k+1) (12)。
CN201410588757.7A 2014-10-28 2014-10-28 Contour registration method used in train guide rail contour measurement based on machine vision Withdrawn CN104296682A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201410588757.7A CN104296682A (en) 2014-10-28 2014-10-28 Contour registration method used in train guide rail contour measurement based on machine vision
CN201510037287.XA CN104655043B (en) 2014-10-28 2015-01-26 One kind is based on contour registration method in the measurement of machine vision train guide rail profile

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410588757.7A CN104296682A (en) 2014-10-28 2014-10-28 Contour registration method used in train guide rail contour measurement based on machine vision

Publications (1)

Publication Number Publication Date
CN104296682A true CN104296682A (en) 2015-01-21

Family

ID=52316518

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201410588757.7A Withdrawn CN104296682A (en) 2014-10-28 2014-10-28 Contour registration method used in train guide rail contour measurement based on machine vision
CN201510037287.XA Expired - Fee Related CN104655043B (en) 2014-10-28 2015-01-26 One kind is based on contour registration method in the measurement of machine vision train guide rail profile

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN201510037287.XA Expired - Fee Related CN104655043B (en) 2014-10-28 2015-01-26 One kind is based on contour registration method in the measurement of machine vision train guide rail profile

Country Status (1)

Country Link
CN (2) CN104296682A (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634273A (en) * 2015-01-08 2015-05-20 成都多极子科技有限公司 Contour registration method based on machine vision train guide rail contour measurement
CN105151083A (en) * 2015-09-28 2015-12-16 成都多极子科技有限公司 Full-automatic train guide rail loss image acquisition device
CN105841633A (en) * 2016-05-23 2016-08-10 电子科技大学 Large-area optical profile measurement calibration method based on double-wave-surface interference fringe array
CN105865369A (en) * 2016-05-23 2016-08-17 电子科技大学 Large-area optical profile measurement device and method based on dual-wave interference fringe array
CN106080661A (en) * 2016-05-27 2016-11-09 电子科技大学 Train guide rail profile undulatory wear measuring method
CN106093071A (en) * 2016-08-22 2016-11-09 合肥德泰科通测控技术有限公司 The three-dimensional population parameter intellectualized detection method of rail
CN106643545A (en) * 2016-10-09 2017-05-10 湖南大学 Calibration method for steel rail profile measured by adopting laser displacement technology
CN107401979A (en) * 2017-08-01 2017-11-28 中国铁道科学研究院 A kind of body oscillating displacement compensation device and method for Contact Line Detection
CN108230394A (en) * 2016-12-14 2018-06-29 中南大学 A kind of orbital image auto-correction method
CN111824206A (en) * 2019-04-14 2020-10-27 宋夏子 In-service steel rail detection device for suburban railway
CN113463720A (en) * 2021-06-30 2021-10-01 广西柳工机械股份有限公司 System and method for identifying contact material of loader bucket
CN114396873A (en) * 2021-12-30 2022-04-26 郑州大学 Method and device for measuring geometric dimension of trapezoidal steel wire
CN114663438A (en) * 2022-05-26 2022-06-24 浙江银轮智能装备有限公司 Track detection method, system, apparatus, storage medium and computer program product
CN114858092A (en) * 2022-05-12 2022-08-05 万岩铁路装备(成都)有限责任公司 Calibration system and calibration method for railway steel rail profile measurement

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105043292B (en) * 2015-06-24 2017-10-27 湘潭大学 A kind of image capturing device analyzed for the quality testing of machine tool guideway face
CN105004280A (en) * 2015-07-13 2015-10-28 成都多极子科技有限公司 Image restoring method in train guiderail contour measurement based on machine vision
CN106020125B (en) * 2016-06-23 2018-10-02 埃夫特智能装备股份有限公司 One kind being based on the continuous industrial robot kinematics algorithm of speed
CN108128323B (en) * 2016-08-30 2019-09-17 大连民族大学 The selection method of laser image characteristic quantity relevant to rail wear amount
JP6920158B2 (en) * 2017-09-29 2021-08-18 株式会社日立ハイテクファインシステムズ Inspection device and inspection method for traffic road equipment
CN108896000A (en) * 2018-07-03 2018-11-27 中国铁道科学研究院集团有限公司 Profile of steel rail alignment schemes and device, rail wear calculation method and device
CN108921904B (en) * 2018-07-13 2021-10-08 云南大学 Method for calibrating pinhole camera by using properties of single ball and asymptote
CN112577445B (en) * 2020-12-18 2021-08-03 深圳市泽峰光电科技有限公司 Log end face profile extraction method of rotary grating

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4082805B2 (en) * 1998-11-04 2008-04-30 株式会社トキメック Rail cross-section measuring device

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634273A (en) * 2015-01-08 2015-05-20 成都多极子科技有限公司 Contour registration method based on machine vision train guide rail contour measurement
CN105151083A (en) * 2015-09-28 2015-12-16 成都多极子科技有限公司 Full-automatic train guide rail loss image acquisition device
CN105841633B (en) * 2016-05-23 2018-03-27 电子科技大学 Large area optical profilometry scaling method based on double wave face interference fringe array
CN105841633A (en) * 2016-05-23 2016-08-10 电子科技大学 Large-area optical profile measurement calibration method based on double-wave-surface interference fringe array
CN105865369A (en) * 2016-05-23 2016-08-17 电子科技大学 Large-area optical profile measurement device and method based on dual-wave interference fringe array
CN105865369B (en) * 2016-05-23 2019-03-15 电子科技大学 Based on double wave face interference fringe array large area optical profilometry device and method
CN106080661A (en) * 2016-05-27 2016-11-09 电子科技大学 Train guide rail profile undulatory wear measuring method
CN106093071A (en) * 2016-08-22 2016-11-09 合肥德泰科通测控技术有限公司 The three-dimensional population parameter intellectualized detection method of rail
CN106643545B (en) * 2016-10-09 2018-10-26 湖南大学 Laser displacement technology measures the calibration method of rail profile
CN106643545A (en) * 2016-10-09 2017-05-10 湖南大学 Calibration method for steel rail profile measured by adopting laser displacement technology
CN108230394A (en) * 2016-12-14 2018-06-29 中南大学 A kind of orbital image auto-correction method
CN107401979A (en) * 2017-08-01 2017-11-28 中国铁道科学研究院 A kind of body oscillating displacement compensation device and method for Contact Line Detection
CN107401979B (en) * 2017-08-01 2020-06-05 中国铁道科学研究院集团有限公司 Vehicle body vibration displacement compensation device and method for catenary detection
CN111824206A (en) * 2019-04-14 2020-10-27 宋夏子 In-service steel rail detection device for suburban railway
CN111824206B (en) * 2019-04-14 2023-01-13 宋夏子 In-service steel rail detection device for suburban railway
CN113463720A (en) * 2021-06-30 2021-10-01 广西柳工机械股份有限公司 System and method for identifying contact material of loader bucket
CN113463720B (en) * 2021-06-30 2023-02-17 广西柳工机械股份有限公司 System and method for identifying contact material of loader bucket
CN114396873A (en) * 2021-12-30 2022-04-26 郑州大学 Method and device for measuring geometric dimension of trapezoidal steel wire
CN114858092A (en) * 2022-05-12 2022-08-05 万岩铁路装备(成都)有限责任公司 Calibration system and calibration method for railway steel rail profile measurement
CN114663438A (en) * 2022-05-26 2022-06-24 浙江银轮智能装备有限公司 Track detection method, system, apparatus, storage medium and computer program product

Also Published As

Publication number Publication date
CN104655043B (en) 2017-09-22
CN104655043A (en) 2015-05-27

Similar Documents

Publication Publication Date Title
CN104296682A (en) Contour registration method used in train guide rail contour measurement based on machine vision
CN104634273A (en) Contour registration method based on machine vision train guide rail contour measurement
CN105004280A (en) Image restoring method in train guiderail contour measurement based on machine vision
CN103217111B (en) A kind of non-contact contact line geometric parameter detection method
CN100480627C (en) Steel rail wearing integrative parameter vehicle-mounted dynamic measuring device and method
CN103557788B (en) A kind of high ferro contact net connects geometric parameter and detects non-contact compensation and Kalman filtering modification method
CN106595630B (en) It is a kind of that drawing system and method are built based on laser navigation Intelligent Mobile Robot
CN101666716B (en) Railway locomotive running attitude measuring method
US20040095585A1 (en) Wheel profile inspection apparatus and method
CN110874596B (en) Intelligent recognition early warning method and system for instrument
CN109060820B (en) Tunnel disease detection method and tunnel disease detection device based on laser detection
CN102441581A (en) Machine vision-based device and method for online detection of structural steel section size
CN107817044B (en) Device and method for measuring plate vibration based on machine vision
CN102589530A (en) Method for measuring position and gesture of non-cooperative target based on fusion of two dimension camera and three dimension camera
CN106080661A (en) Train guide rail profile undulatory wear measuring method
CN104236498A (en) Method for measuring end size of switch rail
AU2020386611A1 (en) Pavement flatness measurement method and system
CN101915852A (en) Velocity measurement method based on stereoscopic vision
CN104897269A (en) Overhead wire vibration frequency measuring system based on linear scanning imaging
CN104239904A (en) Non-contact detection method for external outline of railway vehicle
CN204286406U (en) Use the portable railway Geometrical Parameter Testing system of MEMS gyro instrument
CN110132990A (en) A kind of profile of steel rail detection method based on image recognition
CN103810676A (en) Method for monitoring running speed of steel pipe
KR20100073880A (en) System for measuring railway abrasion using wavelet, and method thereof
CN104878679A (en) Vehicle-mounted rut real-time detection system and detection method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C04 Withdrawal of patent application after publication (patent law 2001)
WW01 Invention patent application withdrawn after publication

Application publication date: 20150121