CN103093088B - The safety evaluation method of abrupt slope and crankcase ventilaton - Google Patents

The safety evaluation method of abrupt slope and crankcase ventilaton Download PDF

Info

Publication number
CN103093088B
CN103093088B CN201310008140.9A CN201310008140A CN103093088B CN 103093088 B CN103093088 B CN 103093088B CN 201310008140 A CN201310008140 A CN 201310008140A CN 103093088 B CN103093088 B CN 103093088B
Authority
CN
China
Prior art keywords
section
road
coordinate system
segmentation
gps
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201310008140.9A
Other languages
Chinese (zh)
Other versions
CN103093088A (en
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.)
Xian Feisida Automation Engineering Co Ltd
Original Assignee
Xian Feisida Automation Engineering Co Ltd
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 Xian Feisida Automation Engineering Co Ltd filed Critical Xian Feisida Automation Engineering Co Ltd
Priority to CN201310008140.9A priority Critical patent/CN103093088B/en
Publication of CN103093088A publication Critical patent/CN103093088A/en
Application granted granted Critical
Publication of CN103093088B publication Critical patent/CN103093088B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses the safety evaluation method of a kind of abrupt slope and crankcase ventilaton, the technical matters cannot measured during for solving the comparatively large and shortage of data of stochastic error and cumulative errors in existing road slope measuring, belongs to road driving security fields.The safety evaluation method of abrupt slope and crankcase ventilaton in the present invention, metrical information according to vehicle GPS carries out segmentation description and piecewise fitting to detection section, when GPS record is not enough to smoothing process or GPS record lacks as travelled in cavern, spline interpolation and smoothing processing are carried out in the section less to GPS survey record in conjunction with IMU metrical information, obtain the fitting surface equation of road, and then calculate the radius of turn of road, the gradient, then utilize recurrence estimation by road forward piecemeal, the impact of cumulative errors in measurement is eliminated by restricted memory length, improve road grade, the computational accuracy of the geometric parameters such as camber, solve the technical matters that existing hill path is difficult to measure, for vehicle supervision department carries out pavement of road monitoring and inspect for acceptance providing foundation, in this section safety traffic, load-carrying is provided to vehicle simultaneously, rate request.

Description

The safety evaluation method of abrupt slope and crankcase ventilaton
Technical field
The present invention relates to a kind of road safety appraisal procedure, particularly the safety evaluation method of a kind of abrupt slope and crankcase ventilaton Dangerous Area, belong to road driving security fields.
Background technology
The crankcase ventilaton on abrupt slope is had to be the multiple area of traffic hazard, higher road accident rate and zig zag, the road alignment of continuous turning is relevant, and the traffic hazard of bend is even more serious than the traffic hazard result of straight-line segment, in addition, on continuous print is long, the road accident rate in descending section is apparently higher than the road accident rate of level road, visible, unreasonable and the bad linear combination of road geometric element design, it is the immediate cause causing traffic hazard to occur, the linear feature of real road is evaluated, in conjunction with the feature of road to the load-carrying of driving vehicle, speed limits, contribute to the generation reducing traffic hazard, improve the safety traffic performance of road.
Relate generally to the measurement to road grade and camber to the linear feature evaluation of real road, the measuring method of the existing pavement of road gradient and camber mainly contains mapping method, based on digital elevation model (DEM) graticule mesh computing method, the computing method etc. measured based on onboard sensor.Mapping method according to the principle of spirit-leveling instrument, the point-to-point measurement of artificial stepping, and then calculate the gradient on road surface; Obtaining the angle of gradient on road surface based on the graticule mesh method of digital elevation model by calculating on road surface certain any gradient, needing the surface equation knowing road surface, or utilizing the coordinate that road surface is put to carry out matching to realize to local surface; The computing method measured based on onboard sensor mainly install different sensors on vehicle, as GPS, gyroscope, laser sensor, displacement transducer etc., according to the metrical information of sensor and the grade information of geometric relationship reckoning road.
Document " based on any direction gradient computing method of DEM, Liu Xuejun, Ren Zhifeng, Wang Yanfang, Jin Bei, Areal Research and Development, 2009,28(4): 139-141 " gives gradient computation model and the calculation process of any direction case site in Grid DEM; Patent 201110125126.2, date of publication 2011.11.23, disclose a kind of car-mounted road surface gradient measuring method based on laser sensor, by installing gyroscope in vehicle centroid position, in the left and right sides symmetry of vehicle body front and back end, laser displacement sensor is installed, utilizes the geometric relationship of road to obtain the road surface horizontal wall inscription gradient and the longitudinal gradient gradient of road according to metrical information.
Need when utilizing mapping method to measure to consume a large amount of manpower and materials, be difficult to carry out, and have cumulative errors in measuring, measurement efficiency is low, DEM graticule mesh computing method, the surface equation of road pavement is needed to carry out matching, during the computing method that onboard sensor carries out measuring mainly travel according to vehicle, the prompting message of measured single locus calculates, because the stochastic error in measuring is inevitable, single-point information is utilized merely to calculate the gradient not necessarily continuous and derivable of the road obtained, and the impact that when not considering in existing Algorithm of road slope estimation to measure, stochastic error and cumulative errors cause, computing method when sensor failure does not have a measurement data situation are not considered yet, the technical matters cannot measured when cannot solve the comparatively large and shortage of data of stochastic error and cumulative errors in measurement.
Summary of the invention
In order to overcome the technological deficiency that in existing road slope measuring method, stochastic error and cumulative errors cannot be measured more greatly and when data is missing, the present invention proposes the safety evaluation method of a kind of abrupt slope and crankcase ventilaton, metrical information according to GPS carries out segmentation description and piecewise fitting to detected section, eliminate the impact of spot measurement stochastic error, when GPS record is not enough to smoothing process or GPS record lacks as travelled in cavern, spline interpolation and smoothing processing are carried out in the section less to GPS survey record in conjunction with IMU metrical information, obtain the fitting surface equation of road, and then calculate the radius of turn of road, the gradient, then utilize recurrence estimation by road forward piecemeal, the impact of cumulative errors in measurement is eliminated by restricted memory length, improve road grade, the computational accuracy of the geometric parameters such as camber, solve the technical matters that existing Mountain Pass Geomtry feature is difficult to measure, for vehicle supervision department carries out pavement of road monitoring and inspect for acceptance providing foundation, in this section safety traffic, load-carrying is provided to vehicle simultaneously, rate request.
The technical solution adopted for the present invention to solve the technical problems is, the safety evaluation method of a kind of abrupt slope and crankcase ventilaton, and its feature comprises the following steps:
1) according to the linear feature of road, segmentation description and modeling are carried out to road: with the initial point measured for true origin sets up geographical rectangular coordinate system, and the change in coordinate axis direction of set up rectangular coordinate system is calibrated, make it consistent with the direction of terrestrial coordinate system, under the measurement data of sensor record is converted to geographical rectangular coordinate system relation, according to the plots changes in detected section, in the following manner segmentation description and modeling are carried out to section:
In formula, represent with the coordinate measured under geographical rectangular coordinate system that initial point is true origin, and be respectively waypoint position on direction, for the number in divided section, the individual ( ) functional relation in section;
When carrying out segmentation to detected section, follow following principle: 1. each segmentation section only has a turning, and be corresponding height single-valued function; 2. turning waits Dangerous Area section length less, and linear section section length can be relatively large;
2) the rectangular coordinate system state equation that segmentation section GPS expresses is set up, for different segmentation sections, because its length is different, the record number of sensor gained measured value is also not quite similar, when GPS record value is not enough to fitting surface equation, or gps data lack as vehicle travel in cavern time, the metrical information of IMU is utilized to carry out spline interpolation and level and smooth estimation to state equation, each segmentation section obtains enough multiple multiple measurement estimated value, realize the spline smooth process to segmentation section, choose for state variable, for input variable, set up the state equation in segmentation section:
In formula, for the coordinate components under geographical rectangular coordinate system, be respectively the speed in lower three directions of vehicle carrier coordinate system, be respectively the longitude of vehicle position, latitude and height, for the angular velocity of rotation in earth relative inertness space , for the radius-of-curvature of prime vertical plane, for the radius-of-curvature of meridian ellipse, for the first excentricity of the earth, for major semi-axis is long, for the length of minor semi-axis, for gravity angular velocity, for the overload in carrier coordinate system three directions, for systematic error, average is zero, and variance is respectively ;
3) according to matching least square index:
According to the following formula piecewise polynomial fitting is carried out to the multiple measurement estimated value adopting GPS/IMU jointly to obtain section rectangular coordinate system:
In formula, for front the number of the multiple measurement estimated value adopted during individual section matching, the multiple measurement estimated value number in order to smooth curve translation backward, be individual multiple measurement estimated value, be the estimated value that section is obtained by model, for the coefficient to be determined of model;
With the function expression in individual segmentation section for example illustrates coefficient determination, if for this section of matching multiple measurement estimated value have individual, in order to represent convenient, be expressed as ,
In formula, for noise vector;
4) is calculated according to the following formula the gradient in individual section and radius of turn
5) according to the design standards of road geometric linear, for certain road grade angle, radius of turn, driving vehicle speed, load-carrying quantitative limitation requirement is provided.
The invention has the beneficial effects as follows, when GPS records disappearance in cavern or record number is not enough, the measurement of another kind of sensor as IMU is utilized to supplement it, multiple onboard sensor is utilized to measure road information, by to the segmentation description of road and piecewise fitting, eliminate the impact of spot measurement on geometric parameter enchancement factors such as road grades, recursion is utilized smoothly to estimate to eliminate the impact of cumulative errors, by the process to a mobile unit metrical information, accurately can obtain the geometric parameter information of road, Dangerous Area is judged according to constructing road standard, for vehicle supervision department is to the pavement detection of road and inspect for acceptance and provide foundation, simultaneously to the load-carrying of Dangerous Area driving vehicle, speed claims, reduce the generation of traffic hazard, improve the safety traffic performance of road.
Below in conjunction with drawings and Examples, invention is elaborated.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of safety evaluation method;
Fig. 2 is the plots changes figure in section in specific embodiment;
Fig. 3 is the longitudinal gradient gradient in section in specific embodiment;
Fig. 4 is the radius of turn in section in specific embodiment.
Embodiment
With reference to accompanying drawing.
1) according to the linear feature of road, segmentation description and modeling are carried out to road: with the initial point measured for true origin sets up geographical rectangular coordinate system, and the change in coordinate axis direction of set up rectangular coordinate system is calibrated, make it consistent with the direction of terrestrial coordinate system, under the measurement data of sensor record is converted to geographical rectangular coordinate system relation, draw out the changing trend diagram in detected section, see accompanying drawing 2, as follows segmentation description and modeling are carried out to section:
In formula, represent with the coordinate measured under geographical rectangular coordinate system that initial point is true origin, and be respectively waypoint position on direction, for the number in divided section, the individual ( ) functional relation in section;
When carrying out segmentation to detected section, follow following principle: 1. each segmentation section only has a turning, and be corresponding height single-valued function; 2. turning waits Dangerous Area section length less, and linear section section length can be relatively large;
2) the rectangular coordinate system state equation that the GPS setting up segmentation section expresses, for different segmentation sections, because its length is different, the record number of sensor gained measured value is also not quite similar, when GPS record value is not enough to matched curve, or gps data lack as vehicle travel in cavern time, the metrical information of IMU is utilized to carry out spline interpolation and level and smooth estimation to state equation, each segmentation section obtains enough multiple multiple measurement estimated value, realize the spline smooth process to segmentation section
Choose for state variable, for input variable, set up the state equation in segmentation section:
In formula, for the coordinate components under geographic coordinate system, be respectively the speed in lower three directions of carrier coordinate system, be respectively the longitude of target position, latitude and height, for the angular velocity of rotation in earth relative inertness space , for the radius-of-curvature of prime vertical plane, for the radius-of-curvature of meridian ellipse, for the first excentricity of the earth, for major semi-axis is long, for the length of minor semi-axis, for gravity angular velocity, for the overload in carrier coordinate system three directions, for systematic error, average is zero, and variance is respectively ;
3) according to matching least square index:
According to the following formula piecewise polynomial fitting is carried out to the multiple measurement estimated value adopting GPS/IMU jointly to obtain section rectangular coordinate system:
In formula, for front the number of the multiple measurement estimated value adopted during individual section matching, the multiple measurement estimated value number in order to smooth curve translation backward, be individual multiple measurement estimated value, be the estimated value that section is obtained by model, for the coefficient to be determined of model;
With the function expression in individual segmentation section for example illustrates coefficient determination, if for this section of matching multiple measurement estimated value have individual, in order to represent convenient, be expressed as ,
In formula, for noise vector;
By carrying out fitting of a polynomial to the data under road section selected rectangular coordinate system, obtain the coefficient of polynomial fitting in section,
4) is calculated according to the following formula the gradient in individual section and radius of turn
Known by calculating, the slope change curve in this section is shown in accompanying drawing 3, and the corresponding radius of turn in section is shown in accompanying drawing 4, visible, and the ruling grade in this section is 1.1428%, and minimal curve radius is 8.1786 meters;
5) according to the design standards of road geometric linear, for certain road grade angle, radius of turn, driving vehicle speed, load-carrying quantitative limitation requirement is provided.

Claims (1)

1. a safety evaluation method for abrupt slope and crankcase ventilaton, its feature comprises the following steps:
1) according to the linear feature of road, segmentation description and modeling are carried out to road: with the initial point measured for true origin sets up geographical rectangular coordinate system, and the change in coordinate axis direction of set up rectangular coordinate system is calibrated, make it consistent with the direction of terrestrial coordinate system, under the measurement data of sensor record is converted to geographical rectangular coordinate system relation, according to the plots changes in detected section, in the following manner segmentation description and modeling are carried out to section:
In formula, represent with the coordinate measured under geographical rectangular coordinate system that initial point is true origin, and be respectively waypoint position on direction, for the number in divided section, the individual ( ) functional relation in section;
When carrying out segmentation to detected section, follow following principle: 1. each segmentation section only has a turning, and be corresponding height single-valued function; 2. turning waits Dangerous Area section length less, and linear section section length can be relatively large;
2) the rectangular coordinate system state equation that segmentation section GPS expresses is set up, for different segmentation sections, because its length is different, the record number of sensor gained measured value is also not quite similar, when GPS record value is not enough to fitting surface equation, or gps data lack as vehicle travel in cavern time, the metrical information of IMU is utilized to carry out spline interpolation and level and smooth estimation to state equation, each segmentation section obtains enough multiple multiple measurement estimated value, realize the spline smooth process to segmentation section, choose for state variable, for input variable, set up the state equation in segmentation section:
In formula, for the coordinate components under geographical rectangular coordinate system, be respectively the speed in lower three directions of vehicle carrier coordinate system, be respectively the longitude of vehicle position, latitude and height, for the angular velocity of rotation in earth relative inertness space , for the radius-of-curvature of prime vertical plane, for the radius-of-curvature of meridian ellipse, for the first excentricity of the earth, for major semi-axis is long, for the length of minor semi-axis, for gravity angular velocity, for the overload in carrier coordinate system three directions, for systematic error, average is zero, and variance is respectively ;
3) according to matching least square index:
According to the following formula piecewise polynomial fitting is carried out to the multiple measurement estimated value adopting GPS/IMU jointly to obtain section rectangular coordinate system:
In formula, for front the number of the multiple measurement estimated value adopted during individual section matching, the multiple measurement estimated value number in order to smooth curve translation backward, be individual multiple measurement estimated value, be the estimated value that section is obtained by model, for the coefficient to be determined of model;
With the function expression in individual segmentation section for example illustrates coefficient determination, if for this section of matching multiple measurement estimated value have individual, in order to represent convenient, be expressed as ,
In formula, for noise vector;
4) is calculated according to the following formula the gradient in individual section and radius of turn
5) according to the design standards of road geometric linear, for certain road grade angle, radius of turn, driving vehicle speed, load-carrying quantitative limitation requirement is provided.
CN201310008140.9A 2013-01-09 2013-01-09 The safety evaluation method of abrupt slope and crankcase ventilaton Expired - Fee Related CN103093088B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310008140.9A CN103093088B (en) 2013-01-09 2013-01-09 The safety evaluation method of abrupt slope and crankcase ventilaton

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310008140.9A CN103093088B (en) 2013-01-09 2013-01-09 The safety evaluation method of abrupt slope and crankcase ventilaton

Publications (2)

Publication Number Publication Date
CN103093088A CN103093088A (en) 2013-05-08
CN103093088B true CN103093088B (en) 2015-08-19

Family

ID=48205647

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310008140.9A Expired - Fee Related CN103093088B (en) 2013-01-09 2013-01-09 The safety evaluation method of abrupt slope and crankcase ventilaton

Country Status (1)

Country Link
CN (1) CN103093088B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105975782A (en) * 2016-05-11 2016-09-28 辽宁工程技术大学 Composite soft fluctuating base waste dump slope form optimizing method

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103760569B (en) * 2013-12-31 2016-03-30 西安交通大学 A kind of drivable region detection method based on laser radar
CN105138733B (en) * 2015-07-30 2018-01-23 河北工业大学 Two-lane highway Traffic safety evaluation method based on driver comfort
CN107544983A (en) * 2016-06-24 2018-01-05 高德信息技术有限公司 A kind of road method for drafting, apparatus and system
US10495493B2 (en) * 2017-06-28 2019-12-03 GM Global Technology Operations LLC Systems and methods for controlling sensing device field of view
CN107351916A (en) * 2017-07-20 2017-11-17 合肥顺顺信息咨询有限公司 A kind of anticipation turning path system based on the steering wheel anglec of rotation
CN108711281A (en) * 2018-03-23 2018-10-26 武汉理工大学 A method of with section actual driving speed data evaluation Highway Safety
CN110103823B (en) * 2019-05-21 2021-06-11 东南大学 Vehicle rollover advance early warning method based on enhanced digital map
CN110186428B (en) * 2019-05-27 2022-03-18 北京泓达九通科技发展有限公司 Road slope determination method and system based on big data curve fitting
CN111693299B (en) * 2020-06-16 2021-07-20 安徽江淮汽车集团股份有限公司 Method and device for formulating driving condition of power test, terminal equipment and storage medium
CN113051994B (en) * 2020-11-18 2021-12-17 庄景江 Real-time lifting system of grip ability based on environmental monitoring
CN113341966B (en) * 2021-05-31 2022-05-17 广州文远知行科技有限公司 Ramp vehicle control method and device, vehicle and storage medium
CN114141011B (en) * 2021-11-23 2022-09-27 北京安融畅信息技术有限公司 Steep slope and continuous downhill road section identification method and potential safety hazard troubleshooting method
CN115424441A (en) * 2022-08-31 2022-12-02 南京慧尔视智能科技有限公司 Road curve optimization method, device, equipment and medium based on microwave radar
CN116578891B (en) * 2023-07-14 2023-10-03 天津所托瑞安汽车科技有限公司 Road information reconstruction method, terminal and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101059851A (en) * 2007-06-05 2007-10-24 天津市市政工程设计研究院 Highway route evaluation method
CN102402644A (en) * 2011-08-11 2012-04-04 西北工业大学 Dynamical model modeling method of vehicle driven on mountainous road

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1764580B1 (en) * 2005-09-14 2008-07-30 C.R.F. Società Consortile per Azioni Method and system for recognizing the sign of the velocity of a vehicle and for estimating the road slope

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101059851A (en) * 2007-06-05 2007-10-24 天津市市政工程设计研究院 Highway route evaluation method
CN102402644A (en) * 2011-08-11 2012-04-04 西北工业大学 Dynamical model modeling method of vehicle driven on mountainous road

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
利用GPS实现大跨桥梁的实时安全监测;过静珺等;《全球定位系统》;20011231;第26卷(第4期);2-8 *
基于GIS的交通安全工作实地调查与评估方法;李平凡等;《交通信息与安全》;20090820;第27卷(第4期);101-103,117 *
基于GPS/IMU的山区道路参数拟合与行车安全分析;王慧丽等;《中国公路学报》;20140630;第27卷(第6期);27-33 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105975782A (en) * 2016-05-11 2016-09-28 辽宁工程技术大学 Composite soft fluctuating base waste dump slope form optimizing method

Also Published As

Publication number Publication date
CN103093088A (en) 2013-05-08

Similar Documents

Publication Publication Date Title
CN103093088B (en) The safety evaluation method of abrupt slope and crankcase ventilaton
CN104061899B (en) A kind of vehicle side inclination angle based on Kalman filtering and angle of pitch method of estimation
Tao et al. Lane marking aided vehicle localization
Jo et al. Road slope aided vehicle position estimation system based on sensor fusion of GPS and automotive onboard sensors
EP3719449A1 (en) Driving condition specific sensor quality index
CN111868798B (en) Generation and update of a lane network map model
JP6260114B2 (en) Traveling route information generation device
CN103162689B (en) The assisted location method of auxiliary vehicle positioning system and vehicle
CN104240536A (en) Lane monitoring method with electronic horizon
Jiang et al. Carloc: Precise positioning of automobiles
CN102226700B (en) Method for matching electronic map of flyover road network
Ben-Arieh et al. Geometric modeling of highways using global positioning system data and B-spline approximation
CN110864696A (en) Three-dimensional high-precision map drawing method based on vehicle-mounted laser inertial navigation data
US10907972B2 (en) 3D localization device
CN111831960B (en) Dynamic load measuring and calculating method of internet wagon based on gradient disturbance identification and elimination
CN108303043A (en) Plant leaf area index detection method and system combined of multi-sensor information
CN103400416A (en) City environment robot navigation method based on multi-layer probabilistic terrain
CN112629530B (en) Vehicle positioning method, device, equipment and storage medium
CN107274721A (en) Many vehicle cooperative localization methods in a kind of intelligent transportation system
CN104216405B (en) The air navigation aid and equipment of field robot
Gupta et al. Road grade estimation using crowd-sourced smartphone data
US10697775B2 (en) Travel distance calculation device, charging system, travel distance calculation method, program, and storage medium
AU2010200766A1 (en) Database for efficient storage of track geometry and feature locations
CN101739840A (en) Poly GPS/INS and transportation image fusion and positioning method
CN202614240U (en) Dead reckoning device combined with vehicle diagnosis system

Legal Events

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

Granted publication date: 20150819