CN108204985A - A kind of shield tunnel mileage localization method - Google Patents

A kind of shield tunnel mileage localization method Download PDF

Info

Publication number
CN108204985A
CN108204985A CN201810038976.6A CN201810038976A CN108204985A CN 108204985 A CN108204985 A CN 108204985A CN 201810038976 A CN201810038976 A CN 201810038976A CN 108204985 A CN108204985 A CN 108204985A
Authority
CN
China
Prior art keywords
row
shield tunnel
denoted
circumferential weld
gray scale
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
CN201810038976.6A
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.)
Shanghai Same Rock Civil Engineering Polytron Technologies Inc
Original Assignee
Shanghai Same Rock Civil Engineering Polytron Technologies Inc
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 Shanghai Same Rock Civil Engineering Polytron Technologies Inc filed Critical Shanghai Same Rock Civil Engineering Polytron Technologies Inc
Priority to CN201810038976.6A priority Critical patent/CN108204985A/en
Publication of CN108204985A publication Critical patent/CN108204985A/en
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06Systems determining position data of a target
    • G01S17/08Systems determining position data of a target for measuring distance only
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497Means for monitoring or calibrating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • G01N21/8851Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
    • G01N2021/8854Grading and classifying of flaws
    • G01N2021/888Marking defects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • G01N21/8851Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
    • G01N2021/8887Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges based on image processing techniques

Abstract

The present invention relates to a kind of shield tunnel mileage localization methods, mainly solve the technical issues of existing tunnel mileage positioning method accuracy is not high, technical scheme of the present invention includes the following steps, the section of jurisdiction circumferential weld in shield tunnel surface image is identified using image-recognizing method, and using the position correction rangefinder mileage location information of each section of jurisdiction circumferential weld, obtain final tunnel mileage location information.Specific steps:S1, shield tunnel surface image and rangefinder mileage location information are obtained by mobile detection apparatus;S2, the shield tunnel surface image according to acquisition identify the section of jurisdiction circumferential weld in image, and export the position of each section of jurisdiction circumferential weld;S3, rangefinder mileage location information is corrected according to the judging result of duct pieces of shield tunnel circumferential weld position, exports final tunnel mileage location information.

Description

A kind of shield tunnel mileage localization method
Technical field
The present invention relates to the positioning of shield tunnel mileage, and a kind of shield tunnel mileage localization method is disclosed, applied to underground Engineering field.
Background technology
With the development of science and technology, the construction scale of shield tunnel project is growing day by day, greatly facilitate people trip and Life, and the shield tunnel project being largely completed has been enter into the maintenance stage, therefore, as shield tunnel quantity continues Increase, the operation state of shield tunnel construction also becomes particularly important with Defect inspection.Traditional shield tunnel surface disease inspection Work is surveyed mainly to detect by artificial detection or installation retainer instrument equipment(Such as total powerstation lays convergence gauge, optical fiber) Come what is realized, but since shield tunnel is more and more, detection portfolio sharply increases, and traditional method cannot expire in efficiency The work requirements of sufficient shield tunnel surface Defect inspection.Relative to traditional method, the shield tunnel based on Digital Image Processing Surface Defect inspection method has the advantages that efficient, precision is high, and stability is high, it has also become shield tunnel surface Defect inspection industry The main stream approach of business.Shield tunnel surface Defect inspection method based on Digital Image Processing is broadly divided into three parts, and A part is digital image acquisition, is generally completed by the mobile equipment such as Tunnel testing vehicle, and second part is Digital Image Processing, Disease recognition carries out the image collected by designed image Processing Algorithm, Part III is that defect information is integrated, according to institute The disease recognition result and corresponding mileage information for having image carry out disease positioning.Wherein, mileage location information directly affects How the precision of disease positioning, therefore, obtain the problem of accurate mileage location information is one very crucial.In recent years, Research is expanded to tunnel mileage localization method there are many scholar.
The patent document of Publication No. CN101943577B discloses a kind of metro tunnel fracture surface deformation detection system, should Profile scanning device, mileage positioning device, data processing system and other auxiliary devices are installed on subway detection vehicle by method On, the positioning of tunnel cross-section information is realized using mileage positioning device.Wherein mileage positioning device is mainly by measuring wheel and rotation Encoder is formed, measuring wheel and rail contact, and rotary encoder connect the two synchronous rotary with measuring wheel by shaft, real Existing tunnel mileage positioning.
The patent document of Publication No. CN105548205A discloses tunnel surface defect positioning method and alignment system, This method is multiple for marking the kilometer stones of tunnel mileage by being set in advance in tunnel surface, then along tunnel extending direction into Row Image Acquisition, then the image based on acquisition identify the kilometer stone in image, obtain the mileage location information in tunnel.
Mileage localization method in tunnel disclosed in document above, tunnel mileage location information source is single, places one's entire reliance upon Measuring wheel or image identification, due to Tunnel testing vehicle attitudes vibration, skidding etc., the mileage information and tunnel of measuring wheel are practical Mileage information is often not consistent, therefore precision is not high.And the kilometer stone in image is identified by image procossing, it needs pre- Kilometer stone is first laid inside tunnel, increases manpower and time cost, also there is some potential safety problems, and precision easily by Picture quality influences.Therefore, it is necessary to study a kind of stability and high efficiency, the tunnel mileage localization methods of high accurancy and precision.
Invention content
The purpose of the invention is to overcome the shortcomings of existing method, it is proposed that a kind of shield tunnel mileage localization method, To improve the precision of the shield tunnel surface disease positioning based on Digital Image Processing.The present invention thinking be:Use image Section of jurisdiction circumferential weld in recognition methods identification shield tunnel surface image, and utilize the position correction rangefinder mileage of each section of jurisdiction circumferential weld Location information obtains final tunnel mileage location information.
The technical scheme is that:A kind of shield tunnel mileage localization method, which is characterized in that include the following steps:
S1, shield tunnel surface image and mileage location information are obtained by mobile detection apparatus;
S2, the shield tunnel surface image according to acquisition identify the section of jurisdiction circumferential weld in image, and export the position of each section of jurisdiction circumferential weld It puts;
S3, it is corrected according to the judging result of duct pieces of shield tunnel circumferential weld position, final tunnel mileage location information is exported, It is denoted as
The particular content of the S1 steps is as follows:
(1)It is single channel gray level image to acquire image by the laser scanner of mobile detection apparatus, and the required precision of image is extremely Reach 3mm/pixel less;
(2)By the preliminary tunnel mileage location information of the rangefinder of mobile detection apparatus, encoder output, it is denoted as
The concrete operations of the step S2 are as follows:
(1)Contrast enhancement processing is carried out to original tunnel image using histogram equalization method, original tunnel image is denoted as, the enhanced image of contrast is denoted as I, and the width of image is high to be denoted as W, H respectively.
(2)According to the image I after contrast enhancement processing, statistics wherein each column gray scale accumulated value.It is as follows:
(a)The initial range of each column gray scale accumulated value in statistical picture I is set, is denoted as, wherein,
(b)According to the statistics initial range of each column gray scale accumulated value of setting, valid data statistical regions are obtained, are denoted as, system The row gray scale accumulated value in valid data statistical regions in image I is counted, is denoted as,, computational methods are as follows:
(1)
WhereinExist for image IThe gray value of position, and,
(3)Processing is zoomed in and out to the row gray scale accumulated value after statistics, scales it setting range.Concrete operations are such as Under:
(a)The maximum value of row gray scale accumulated value after setting scaling, is denoted as
(b)The row gray scale accumulated value after scaling is calculated, is denoted as, the scaling computational methods of row gray scale accumulated value are as follows:
(2)
WhereinFor the maximum value in the row gray scale accumulated value before scaling,
(4)Duct pieces of shield tunnel circumferential weld decision threshold is calculated according to the row gray scale accumulated value after scaling, is denoted as;Specific behaviour Make as follows:
(a)Row gray scale accumulated value after scaling needs to cut into ranks, cuts columns and is denoted as, in order to eliminate image The zero interference on boundary, the column position index after row are cut are denoted as,
(b)The maximum value and minimum value in the row gray scale accumulated value after row are cut are calculated, is denoted as respectively, Computational methods are as follows:
(3)
(4)
(c)The maximum value and minimum value in row gray scale accumulated value after being cut according to row determines that duct pieces of shield tunnel circumferential weld judges Threshold value is denoted as, computational methods are as follows:
(5)
(6)
WhereinIt is punishment parameter.
(5)Traverse all positions effectively arranged, and according toWhether the position for judging the currently active row is shield tunnel Position where piece circumferential weld;Concrete operations are as follows:
(a)The row gray scales for choosing several sampling row at left and right sides of the currently active row add up and data, the position effectively arranged are denoted as,, wherein N is sampling columns,,It is that the sampling that sampling arranges at left and right sides of the currently active row is inclined Shifting amount.Starting sample column position at left and right sides of the currently active row can be inclined according to the sampling that sampling arranges at left and right sides of the currently active row Shifting amount determines that the starting sample column position at left and right sides of the currently active row is denoted as respectively, wherein,,For the position of the currently active row,.The row gray scale of sampling row at left and right sides of the currently active row adds up Value is denoted as respectively, wherein,
(b)Statistics is located at the average value of the row gray scale accumulated value of several samplings row at left and right sides of the currently active row, is denoted as respectively,, computational methods are as follows:
(7)
(8)
(c)Calculate the currently active row row gray scale add up and at left and right sides of it several samplings row row gray scale accumulated value it is equal The difference of value, is denoted as,, computational methods are as follows:
(9)
(10)
(d)If the mean value of the row gray scale accumulated value of the currently active row and the row gray scale accumulated value of several sampling row in its left and right sides Difference be more than duct pieces of shield tunnel circumferential weld decision threshold, then where the position of the currently active row is duct pieces of shield tunnel circumferential weld Position, record this position;Conversely, the position of the currently active row is not the position where duct pieces of shield tunnel circumferential weld, do not record This position;
(e)The duct pieces of shield tunnel circumferential weld position tentatively recorded is post-processed, post-processing approach is:One search model of setting It encloses, is denoted as SR, in the range of SR, only retain an optimal duct pieces of shield tunnel circumferential weld position;
(f)It repeats above-mentioned(a)~(e)Step is completed until the positions effectively arranged all in image traverse.
(6)The judging result of final duct pieces of shield tunnel circumferential weld position is exported, is denoted as
The concrete operations of the step S3 are as follows:
(1)InitializationFor
(2)TraversalThe position of middle duct pieces of shield tunnel circumferential weld and, obtain the position pair of duct pieces of shield tunnel circumferential weld The mileage location information answered.
(3)According to the mileage information of duct pieces of shield tunnel length correction, shield tunnel leaf length is denoted as, meter Calculation method is as follows:
Wherein,, N is the duct pieces of shield tunnel circumferential weld number detected,For i-th of section of jurisdiction circumferential weld position Corresponding mileage location information.
(4)It repeats above-mentioned(2)~(3)Step, untilIn all duct pieces of shield tunnel circumferential welds position traversal complete.
Compared with prior art, the present invention has the following advantages:(1)Present invention incorporates the high efficiency and image of rangefinder The high accurancy and precision of identification realizes the mileage positioning of shield tunnel, and comprehensive performance is stronger.(2)Mathematics shape based on image gradation data State feature identifies that, without complicated image processing operations, computational efficiency is high to do duct pieces of shield tunnel circumferential weld.(3)Profit of the invention With the tunnel mileage location information of duct pieces of shield tunnel circumferential weld recognition methods correction rangefinder, robustness higher.
Description of the drawings
Fig. 1 is the system flow schematic diagram of the method for the present invention.
Fig. 2 is the section of jurisdiction circumferential weld identification process schematic diagram of the method for the present invention.
Fig. 3 is the circumferential weld resolution principle schematic diagram of the method for the present invention.
Fig. 4 is the tunnel mileage location information bearing calibration schematic diagram of the method for the present invention.
Specific embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings.
With reference to the accompanying drawings 1, the present invention is a kind of shield tunnel mileage localization method, is included the following steps:
S1, shield tunnel surface image and mileage location information are obtained by mobile detection apparatus;
(1)Mobile detection apparatus mainly includes image collecting device, mileage positioning device;
(2)Image collecting device use is single channel ash inspiring confidence in scientific and technological 9012 laser scanners of Z+F PROFILER, acquisition image Image is spent, the required precision of image is at least up to 3mm/pixel;
(3)Mileage positioning device mainly includes rangefinder, encoder, exports preliminary tunnel mileage location information, is denoted as; Wherein rangefinder is laser range finder, and encoder is incremental optical-electricity encoder.
S2, the shield tunnel surface image according to acquisition identify the section of jurisdiction circumferential weld in image, and export each section of jurisdiction circumferential weld Position;
With reference to the accompanying drawings 2, the concrete operations of the step S2 are as follows:
(1)Contrast enhancement processing is carried out to original tunnel image using histogram equalization method, original tunnel image is denoted as, the enhanced image of contrast is denoted as I, and the width of image is high to be denoted as W, H respectively.
(2)According to the image I after contrast enhancement processing, statistics wherein each column gray scale accumulated value.It is as follows:
(a)The initial range of each column gray scale accumulated value in statistical picture I is set, is denoted as, wherein,
(b)According to the statistics initial range of each column gray scale accumulated value of setting, valid data statistical regions are obtained, are denoted as, system Count imageRow gray scale accumulated value in middle valid data statistical regions, is denoted as,, computational methods are as follows:
(1)
WhereinFor image The gray value of position, and,
Valid data statistical regions control parameter can be defined as follows:
struct GapRecRegCtlPara
{
int nVStart;// setting data extraction scope:Vertical direction starting point
int nVEnd;// setting data extraction scope:Vertical direction terminal
}。
(3)Processing is zoomed in and out to the row gray scale accumulated value after statistics, scales it setting range.Concrete operations are such as Under:
(a)The maximum value of row gray scale accumulated value after setting scaling, is denoted as
(b)The row gray scale accumulated value after scaling is calculated, is denoted as, the scaling computational methods of row gray scale accumulated value are as follows:
)
WhereinFor the maximum value in the row gray scale accumulated value before scaling,
(4)With reference to the accompanying drawings 3, duct pieces of shield tunnel circumferential weld decision threshold, note are calculated according to the row gray scale accumulated value after scaling For.Concrete operations are as follows:
(a)Row gray scale accumulated value after scaling needs to cut into ranks, cuts columns and is denoted as, in order to eliminate image The zero interference on boundary, the column position index after row are cut are denoted as,
(b)The maximum value and minimum value in the row gray scale accumulated value after row are cut are calculated, is denoted as respectively, Computational methods are as follows:
(3)
(4)
(c)The maximum value and minimum value in row gray scale accumulated value after being cut according to row determines that duct pieces of shield tunnel circumferential weld judges Threshold value is denoted as, computational methods are as follows:
(5)
(6)
WhereinIt is punishment parameter, to control the size of duct pieces of shield tunnel circumferential weld decision threshold, which can generally be set as 0.5;
The control parameter for calculating duct pieces of shield tunnel circumferential weld decision threshold can be defined as follows:
struct GapRecThCtlPara
{
double dPenaltyCoef;// duct pieces of shield tunnel circumferential weld decision threshold
int nBoundaryClip;// image boundary cuts range
}。
(5)Traverse all positions effectively arranged, and according toWhether the position for judging the currently active row is shield tunnel Position where piece circumferential weld.Concrete operations are as follows:
(a)The row gray scales for choosing several sampling row at left and right sides of the currently active row add up and data, the position effectively arranged are denoted as,, wherein N is sampling columns,,It is that the sampling that sampling arranges at left and right sides of the currently active row is inclined Shifting amount.Starting sample column position at left and right sides of the currently active row can be inclined according to the sampling that sampling arranges at left and right sides of the currently active row Shifting amount determines that the starting sample column position at left and right sides of the currently active row is denoted as respectively, wherein,,For the position of the currently active row,.The row gray scale of sampling row at left and right sides of the currently active row adds up Value is denoted as respectively, wherein,
(b)Statistics is located at the average value of the row gray scale accumulated value of several samplings row at left and right sides of the currently active row, is denoted as respectively,, computational methods are as follows:
(7)
(8)
(c)Calculate the currently active row row gray scale add up and at left and right sides of it several samplings row row gray scale accumulated value it is equal The difference of value, is denoted as,, computational methods are as follows:
(9)
(10)
(d)If the mean value of the row gray scale accumulated value of the currently active row and the row gray scale accumulated value of several sampling row in its left and right sides Difference be more than duct pieces of shield tunnel circumferential weld decision threshold, then where the position of the currently active row is duct pieces of shield tunnel circumferential weld Position, record this position;Conversely, the position of the currently active row is not the position where duct pieces of shield tunnel circumferential weld, do not record This position;
(e)The duct pieces of shield tunnel circumferential weld position tentatively recorded is post-processed, post-processing approach is:One search model of setting It encloses, is denoted as SR, in the range of SR, only retain an optimal duct pieces of shield tunnel circumferential weld position;
(f)It repeats above-mentioned(a)~(e)Step is completed until the positions effectively arranged all in image traverse.
(6)The judging result of final duct pieces of shield tunnel circumferential weld position is exported, is denoted as
The control parameter of duct pieces of shield tunnel circumferential weld judgement can be defined as follows:
struct GapRecJudgeCtlPara
{
int nShiftLR;// sampling pixel points offset
int nPointNum;// sampling pixel points quantity
}。
S3, it is corrected according to the judging result of duct pieces of shield tunnel circumferential weld position, export final tunnel mileage positioning Information is denoted as
With reference to the accompanying drawings 4, the concrete operations of the step S3 are as follows:
(1)InitializationFor
(2)TraversalThe position of middle duct pieces of shield tunnel circumferential weld and, the position for obtaining duct pieces of shield tunnel circumferential weld corresponds to Mileage location information;
(3)According to the mileage information of duct pieces of shield tunnel length correction, shield tunnel leaf length is denoted as, which is Standard setting, generally 1.5m are built according to shield tunnel.Computational methods are as follows:
Wherein,,For the duct pieces of shield tunnel circumferential weld number detected,It isA section of jurisdiction circumferential weld position Corresponding mileage location information;
(4)It repeats above-mentioned(2)~(3)Step, untilIn all duct pieces of shield tunnel circumferential welds position traversal complete.
The preferred embodiment of the present invention described in detail above.It should be appreciated that those of ordinary skill in the art without Creative work is needed according to the present invention can to conceive and makes many modifications and variations.Therefore, all technologies in the art Personnel are available by logical analysis, reasoning, or a limited experiment on the basis of existing technology under this invention's idea Technical solution, all should be in the protection domain being defined in the patent claims.

Claims (5)

1. a kind of shield tunnel mileage localization method, which is characterized in that identify shield tunnel exterior view using image-recognizing method Section of jurisdiction circumferential weld as in, and using the position correction rangefinder mileage location information of each section of jurisdiction circumferential weld, obtain in final tunnel Journey location information.
2. a kind of shield tunnel mileage localization method according to claim 1, which is characterized in that include the following steps:
S1, shield tunnel surface image and rangefinder mileage location information are obtained by mobile detection apparatus;S2, according to acquisition Shield tunnel surface image identifies the section of jurisdiction circumferential weld in image, and exports the position of each section of jurisdiction circumferential weld;
S3, rangefinder mileage location information is corrected according to the judging result of duct pieces of shield tunnel circumferential weld position, exports final tunnel Road mileage location information.
3. a kind of shield tunnel mileage localization method according to claim 2, which is characterized in that the S1 steps it is specific Content is as follows:
(1)It is single channel gray level image to acquire image by the laser scanner of mobile detection apparatus, and the required precision of image is extremely Reach 3mm/pixel less;
(2)By the preliminary tunnel mileage location information of the rangefinder of mobile detection apparatus, encoder output, it is denoted as
A kind of 4. shield tunnel mileage localization method according to claim 2, which is characterized in that the concrete operations of the S2 It is as follows:
(1)Contrast enhancement processing is carried out to original tunnel image using histogram equalization method, original tunnel image is denoted as, the enhanced image of contrast is denoted as, the width of image is high to be denoted as respectively,
(2)According to the image after contrast enhancement processing, statistics wherein each column gray scale accumulated value;It is as follows:
(a)Set statistical pictureThe initial range of middle each column gray scale accumulated value, is denoted as, wherein,
(b)According to the statistics initial range of each column gray scale accumulated value of setting, valid data statistical regions are obtained, are denoted as, system Count imageRow gray scale accumulated value in middle valid data statistical regions, is denoted as,, computational methods are as follows:
WhereinFor image The gray value of position, and,
(3)Processing is zoomed in and out to the row gray scale accumulated value after statistics, scales it setting range;It is as follows:
(a)The maximum value of row gray scale accumulated value after setting scaling, is denoted as
(b)The row gray scale accumulated value after scaling is calculated, is denoted as, the scaling computational methods of row gray scale accumulated value are as follows:
WhereinFor the maximum value in the row gray scale accumulated value before scaling,
(4)Duct pieces of shield tunnel circumferential weld decision threshold is calculated according to the row gray scale accumulated value after scaling, is denoted as;Specific steps are such as Under:
(a)Row gray scale accumulated value after scaling needs to cut into ranks, cuts columns and is denoted as, in order to eliminate image side The zero interference on boundary, the column position index after row are cut are denoted as,
(b)The maximum value and minimum value in the row gray scale accumulated value after row are cut are calculated, is denoted as respectively, meter Calculation method is as follows:
(c)The maximum value and minimum value in row gray scale accumulated value after being cut according to row determines that duct pieces of shield tunnel circumferential weld judges Threshold value is denoted as, computational methods are as follows:
WhereinIt is punishment parameter;
(5)Traverse all positions effectively arranged, and according toWhether the position for judging the currently active row is duct pieces of shield tunnel circumferential weld The position at place;It is as follows:
(a)The row gray scales for choosing several sampling row at left and right sides of the currently active row add up and data, the position effectively arranged are denoted as,, wherein N is sampling columns,,It is that the sampling that sampling arranges at left and right sides of the currently active row is inclined Shifting amount;Starting sample column position at left and right sides of the currently active row can be inclined according to the sampling that sampling arranges at left and right sides of the currently active row Shifting amount determines that the starting sample column position at left and right sides of the currently active row is denoted as respectively, wherein,,For the position of the currently active row,;The row gray scale of sampling row at left and right sides of the currently active row adds up Value is denoted as respectively, wherein,
(b)Statistics is located at the average value of the row gray scale accumulated value of several samplings row at left and right sides of the currently active row, is denoted as respectively,, computational methods are as follows:
(c)Calculate the currently active row row gray scale add up and at left and right sides of it several samplings row row gray scale accumulated value it is equal The difference of value, is denoted as,, computational methods are as follows:
(d)If the mean value of the row gray scale accumulated value of the currently active row and the row gray scale accumulated value of several sampling row in its left and right sides Difference be more than duct pieces of shield tunnel circumferential weld decision threshold, then where the position of the currently active row is duct pieces of shield tunnel circumferential weld Position, record this position;Conversely, the position of the currently active row is not the position where duct pieces of shield tunnel circumferential weld, do not record This position;
(e)The duct pieces of shield tunnel circumferential weld position tentatively recorded is post-processed, post-processing approach is:One search model of setting It encloses, is denoted as, in the range of SR, only retain an optimal duct pieces of shield tunnel circumferential weld position;
(f)It repeats above-mentioned(a)~(e)Step is completed until the positions effectively arranged all in image traverse;
(6)The judging result of final duct pieces of shield tunnel circumferential weld position is exported, is denoted as
A kind of 5. shield tunnel mileage localization method according to claim 2, which is characterized in that the specific steps of the S3 It is as follows:
(1)Initialize final tunnel mileage location informationFor rangefinder mileage location information
(2)TraversalThe position of middle duct pieces of shield tunnel circumferential weld and, the position for obtaining duct pieces of shield tunnel circumferential weld is corresponding Mileage location information;
(3)According to the mileage information of duct pieces of shield tunnel length correction, shield tunnel leaf length is denoted as, calculating side Method is as follows:
Wherein,,For the duct pieces of shield tunnel circumferential weld number detected,It isA section of jurisdiction circumferential weld position pair The mileage location information answered;
(4)It repeats above-mentioned(2)~(3)Step, untilIn all duct pieces of shield tunnel circumferential welds position traversal complete.
CN201810038976.6A 2018-01-16 2018-01-16 A kind of shield tunnel mileage localization method Withdrawn CN108204985A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810038976.6A CN108204985A (en) 2018-01-16 2018-01-16 A kind of shield tunnel mileage localization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810038976.6A CN108204985A (en) 2018-01-16 2018-01-16 A kind of shield tunnel mileage localization method

Publications (1)

Publication Number Publication Date
CN108204985A true CN108204985A (en) 2018-06-26

Family

ID=62606259

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810038976.6A Withdrawn CN108204985A (en) 2018-01-16 2018-01-16 A kind of shield tunnel mileage localization method

Country Status (1)

Country Link
CN (1) CN108204985A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109029373A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 The synthesized positioning method and device of vehicle are detected in a kind of tunnel
CN109029372A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 A kind of detection vehicle localization method and device
CN109029374A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 It is a kind of to analyze the method and device positioned to measuring car using visual pattern
CN109029377A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 It is a kind of using visual analysis to detection car weight positioning square law device and system
CN110161043A (en) * 2019-05-10 2019-08-23 同济大学 A kind of subway tunnel structure synthetic detection vehicle
CN112819070A (en) * 2021-01-29 2021-05-18 中国电建集团华东勘测设计研究院有限公司 Tunnel segment circular seam linearization identification and mileage ring number marking method
DE102022103565A1 (en) 2022-02-16 2023-08-17 Bayerische Motoren Werke Aktiengesellschaft Longitudinal positioning of a vehicle

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101323971B1 (en) * 2012-09-03 2013-10-31 현대엠엔소프트 주식회사 A method for automatic generation of tunnel information using a mobile mapping system
CN104266592A (en) * 2014-09-29 2015-01-07 樊晓莉 Sleeper-based displacement measuring method for mobile device inside tunnel
CN104266591A (en) * 2014-09-29 2015-01-07 樊晓莉 Displacement detection method for moving device in tunnel
CN105737849A (en) * 2015-12-24 2016-07-06 武汉武大卓越科技有限责任公司 Calibration method of relative position between laser scanner and camera on tunnel car
CN107392216A (en) * 2017-07-27 2017-11-24 上海同岩土木工程科技股份有限公司 A kind of duct pieces of shield tunnel circumferential weld method for quickly identifying based on gradation data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101323971B1 (en) * 2012-09-03 2013-10-31 현대엠엔소프트 주식회사 A method for automatic generation of tunnel information using a mobile mapping system
CN104266592A (en) * 2014-09-29 2015-01-07 樊晓莉 Sleeper-based displacement measuring method for mobile device inside tunnel
CN104266591A (en) * 2014-09-29 2015-01-07 樊晓莉 Displacement detection method for moving device in tunnel
CN105737849A (en) * 2015-12-24 2016-07-06 武汉武大卓越科技有限责任公司 Calibration method of relative position between laser scanner and camera on tunnel car
CN107392216A (en) * 2017-07-27 2017-11-24 上海同岩土木工程科技股份有限公司 A kind of duct pieces of shield tunnel circumferential weld method for quickly identifying based on gradation data

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109029373A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 The synthesized positioning method and device of vehicle are detected in a kind of tunnel
CN109029372A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 A kind of detection vehicle localization method and device
CN109029374A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 It is a kind of to analyze the method and device positioned to measuring car using visual pattern
CN109029377A (en) * 2018-07-16 2018-12-18 银河水滴科技(北京)有限公司 It is a kind of using visual analysis to detection car weight positioning square law device and system
CN110161043A (en) * 2019-05-10 2019-08-23 同济大学 A kind of subway tunnel structure synthetic detection vehicle
CN110161043B (en) * 2019-05-10 2021-03-26 同济大学 Subway tunnel structure comprehensive detection vehicle
CN112819070A (en) * 2021-01-29 2021-05-18 中国电建集团华东勘测设计研究院有限公司 Tunnel segment circular seam linearization identification and mileage ring number marking method
DE102022103565A1 (en) 2022-02-16 2023-08-17 Bayerische Motoren Werke Aktiengesellschaft Longitudinal positioning of a vehicle

Similar Documents

Publication Publication Date Title
CN108204985A (en) A kind of shield tunnel mileage localization method
CN111899288B (en) Tunnel leakage water area detection and identification method based on infrared and visible light image fusion
CN109708615B (en) Subway tunnel clearance dynamic detection method based on laser scanning
CN102297660B (en) Measuring method of shield tunnel lining segment seam open width and apparatus thereof
CN100504288C (en) Article geometrical size measuring device and method based on multi-source image fusion
CN107632022A (en) A kind of detection method of surface flaw of steel rail and device based on data processing
CN105091922A (en) Virtual watch plate based pointer reading identifying method
CN112433203A (en) Lane alignment detection method based on millimeter wave radar data
CN103778681A (en) Vehicle-mounted high-speed road inspection system and data acquisition and processing method
CN102768022A (en) Tunnel surrounding rock deformation detection method adopting digital camera technique
CN103486971B (en) A kind of subway tunnel fracture width detection and correcting algorithm
CN105890634B (en) A kind of absolute type encoder and its decoding measuring method
CN106918598B (en) Bridge pavement strain and crack detection analysis system and method based on digital picture
CN105387801A (en) Subway tunnel segment dislocation quantity detection method
CN102305795A (en) Method for positioning tiny crack on surface of concrete
CN103837087A (en) Automatic pantograph detection method based on active shape model
CN111539436B (en) Rail fastener positioning method based on straight template matching
Huang et al. Rapid surface damage detection equipment for subway tunnels based on machine vision system
CN107392216B (en) Method for quickly identifying circumferential seams of shield tunnel segments based on gray data
CN105674893B (en) Absolute grating scale based on cmos image sensor and its measurement method
CN114459539B (en) Double-line tunnel excavation safety on-line monitoring analysis system based on artificial intelligence
CN110793501A (en) Subway tunnel clearance detection method
CN103810676A (en) Method for monitoring running speed of steel pipe
CN107022933B (en) For the orbit adjusting method during high-speed railway combined test or during operation
CN108596968B (en) Sleeper counting method based on track 3D depth image

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WW01 Invention patent application withdrawn after publication

Application publication date: 20180626

WW01 Invention patent application withdrawn after publication