WO2016056569A1 - 路面状態判別システム - Google Patents
路面状態判別システム Download PDFInfo
- Publication number
- WO2016056569A1 WO2016056569A1 PCT/JP2015/078398 JP2015078398W WO2016056569A1 WO 2016056569 A1 WO2016056569 A1 WO 2016056569A1 JP 2015078398 W JP2015078398 W JP 2015078398W WO 2016056569 A1 WO2016056569 A1 WO 2016056569A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- road surface
- surface state
- value
- feature
- vehicle
- Prior art date
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W40/00—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
- B60W40/02—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to ambient conditions
- B60W40/06—Road conditions
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60C—VEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
- B60C19/00—Tyre parts or constructions not otherwise provided for
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60T—VEHICLE BRAKE CONTROL SYSTEMS OR PARTS THEREOF; BRAKE CONTROL SYSTEMS OR PARTS THEREOF, IN GENERAL; ARRANGEMENT OF BRAKING ELEMENTS ON VEHICLES IN GENERAL; PORTABLE DEVICES FOR PREVENTING UNWANTED MOVEMENT OF VEHICLES; VEHICLE MODIFICATIONS TO FACILITATE COOLING OF BRAKES
- B60T8/00—Arrangements for adjusting wheel-braking force to meet varying vehicular or ground-surface conditions, e.g. limiting or varying distribution of braking force
- B60T8/17—Using electrical or electronic regulation means to control braking
- B60T8/172—Determining control parameters used in the regulation, e.g. by calculations involving measured or detected parameters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
-
- G—PHYSICS
- G07—CHECKING-DEVICES
- G07C—TIME OR ATTENDANCE REGISTERS; REGISTERING OR INDICATING THE WORKING OF MACHINES; GENERATING RANDOM NUMBERS; VOTING OR LOTTERY APPARATUS; ARRANGEMENTS, SYSTEMS OR APPARATUS FOR CHECKING NOT PROVIDED FOR ELSEWHERE
- G07C5/00—Registering or indicating the working of vehicles
- G07C5/008—Registering or indicating the working of vehicles communicating information to a remotely located station
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
- G08G1/0112—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
- G08G1/0141—Measuring and analyzing of parameters relative to traffic conditions for specific applications for traffic information dissemination
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60C—VEHICLE TYRES; TYRE INFLATION; TYRE CHANGING; CONNECTING VALVES TO INFLATABLE ELASTIC BODIES IN GENERAL; DEVICES OR ARRANGEMENTS RELATED TO TYRES
- B60C19/00—Tyre parts or constructions not otherwise provided for
- B60C2019/004—Tyre sensors other than for detecting tyre pressure
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60T—VEHICLE BRAKE CONTROL SYSTEMS OR PARTS THEREOF; BRAKE CONTROL SYSTEMS OR PARTS THEREOF, IN GENERAL; ARRANGEMENT OF BRAKING ELEMENTS ON VEHICLES IN GENERAL; PORTABLE DEVICES FOR PREVENTING UNWANTED MOVEMENT OF VEHICLES; VEHICLE MODIFICATIONS TO FACILITATE COOLING OF BRAKES
- B60T2210/00—Detection or estimation of road or environment conditions; Detection or estimation of road shapes
- B60T2210/10—Detection or estimation of road conditions
- B60T2210/12—Friction
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60T—VEHICLE BRAKE CONTROL SYSTEMS OR PARTS THEREOF; BRAKE CONTROL SYSTEMS OR PARTS THEREOF, IN GENERAL; ARRANGEMENT OF BRAKING ELEMENTS ON VEHICLES IN GENERAL; PORTABLE DEVICES FOR PREVENTING UNWANTED MOVEMENT OF VEHICLES; VEHICLE MODIFICATIONS TO FACILITATE COOLING OF BRAKES
- B60T2210/00—Detection or estimation of road or environment conditions; Detection or estimation of road shapes
- B60T2210/30—Environment conditions or position therewithin
- B60T2210/36—Global Positioning System [GPS]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2422/00—Indexing codes relating to the special location or mounting of sensors
- B60W2422/70—Indexing codes relating to the special location or mounting of sensors on the wheel or the tire
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/02—Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
Definitions
- the present invention relates to a system for determining the state of a road surface during traveling, and in particular, a road surface state determination system for determining a road surface state using vehicle information of a plurality of vehicles traveling in the same place within a predetermined time. About.
- the vibration of the tire during traveling is detected, and the road surface state is estimated from the time series waveform of the detected tire vibration (for example, Patent Documents 1 to 4). 3) or a method of detecting a tire generated sound generated from a tire and estimating a road surface state from a sound pressure level of the detected tire generated sound (see, for example, Patent Document 4).
- JP 2011-242303 A JP 2005-59800 A JP 2014-35279 A JP-A-8-261993
- the present invention has been made in view of the conventional problems, and an object thereof is to provide a road surface state determination system capable of improving the road surface state determination accuracy.
- a road surface state determination system includes an in-vehicle sensor that acquires vehicle information that is information on behavior of a vehicle during travel, and a road surface state (hereinafter referred to as the vehicle surface that the vehicle travels from vehicle information acquired by the in-vehicle sensor).
- FIG. It is a figure which shows the structure of the road surface state determination system which concerns on this Embodiment 1.
- FIG. It is a figure which shows the example of arrangement
- FIG. It is a figure which shows the structure of the road surface state determination apparatus which concerns on this Embodiment 2.
- FIG. It is a figure which shows the structure of the road surface state determination system which concerns on this Embodiment 3.
- FIG. It is a figure which shows the structure of the road surface state determination apparatus which concerns on this Embodiment 4.
- FIG. It is a figure which shows the structure of the road surface state determination system which concerns on this Embodiment 1.
- FIG. 1 is a functional block diagram showing the configuration of the road surface state determination system 1 according to the first embodiment.
- the vehicles W 1 to W n are vehicles that pass through a place within a predetermined range within a predetermined time.
- Reference numeral 20 denotes a server including a receiver 21, data storage means 22, and a transmitter 23.
- Reference numeral 30 denotes a road surface state determination device that determines the state of a road surface within a predetermined time and within a predetermined range.
- the server 20 and the road surface state determination device 30 are provided in the road surface state management center 2.
- the feature quantity calculation means 13 and the road surface state determination device 30 are configured by software of a computer, for example.
- the acceleration sensor 11 is disposed at a substantially central portion on the tire air chamber 42 side of the inner liner portion 41 of the tire 40, and inputs vibration (tire vibration) input to the tread 43 of the tire 40 from the road surface.
- the wheel speed sensor 12 detects the rotational speed (wheel speed) of the wheel.
- a well-known electromagnetic induction type wheel speed sensor that detects the rotational angle of the wheel can be used.
- Feature calculating unit 13, a vibration waveform detecting unit 13a, a region signal extraction unit 13b, and a feature amount calculation unit 13c the time-series waveform of the detected tire vibration by the acceleration sensor 11, the travel of the vehicle W i
- a feature amount for estimating a road surface state (road surface state) to be calculated is calculated.
- the vibration waveform detector 13 a detects an acceleration waveform for one rotation of the tire using the wheel speed detected by the wheel speed sensor 12.
- the region signal extraction unit 13b generates an acceleration waveform for one rotation of the tire, a step region R1 before the stepping side peak that appears at the stepping end, a stepping region R2 that forms the stepping side peak, and the stepping step.
- a kick-out region R3 between a peak on the side and a kick-out peak appearing at the kick-out end, a kick-out region R4 forming the kick-out peak, and a rear of the kick-out peak Are divided into the post-kick region R5, and the time-series waveforms of the vibration levels in the regions R1 to R5 are extracted.
- the feature amount calculation unit 13c passes each time-series waveform of each of the regions R1 to R5 through a bandpass filter, and sets m band values that are magnitudes of vibration components in a predetermined region and a predetermined frequency region.
- the acceleration waveform for one rotation of the tire may be divided into three areas, a pre-depression area, a ground contact area, and a post-kick area.
- the vibration level in the frequency band, the band value of the fluctuation spectrum of the air pressure detected by the pressure sensor installed in the tire chamber, or the band power of the tire generated sound collected by the microphone installed in the tire or in the front frame of the rear wheel The value etc. can be mentioned.
- the vehicle sensor, bidirectionally vibration sensor is provided to detect the vibration of the vibration and the tire radial direction of the tire circumferential direction, may be a specific frequency component of the vibration waveform as the feature quantity X ij of the vehicle W i.
- the feature vector X i has two components. Further, as the feature amount, tread shear stress arranged in the tread block, circumferential tire inner surface strain detected by a strain sensor provided on the inner surface of the tire, time-series data of wheel speed detected by the wheel speed sensor, and the like are used. be able to.
- GPS device 14 is provided with a GPS antenna (not shown) and a receiver, and acquires the position data of the vehicle W i.
- time data an extraction time of an acceleration time series waveform, an acquisition time of position data, a data transmission time, or the like may be used. Since these times are almost the same time, there is no problem even if any time is used as the time data.
- the transmitter 16 transmits m feature quantity X ij data collected by the vehicle information collecting means 15, vehicle position data, and time data together with a vehicle ID for identifying the vehicle, not shown. It transmits to the server 20 of the road surface state management center 2 from the antenna.
- the receiver 17 receives the road surface condition data at a place within a predetermined time and within a predetermined range, which is determined by the road surface state determination device 30 of the road surface state management center 2 and transmitted from the transmitter 23 of the server 20.
- the received road surface state data is displayed on a monitor provided in the vehicle, so that the driver can know the road surface state in the predetermined time and within the predetermined range. Moreover, if the road surface state data is fed back to the vehicle control, the traveling safety of the vehicle can be improved.
- Data of road surface conditions at a place within a predetermined time and within a predetermined range determined by the state determination device 30 is transmitted to the registered vehicle.
- Data of the feature quantity X ij of each vehicle W i is classified and stored in the data storage means 22.
- the predetermined time indicates a range of 1 to 5 minutes including a preset determination time
- the predetermined range indicates, for example, a range including a preset place on a road map.
- the predetermined range for example, by using grids having a predetermined size frame on the road map for navigation, because it is easy to ensure the position consistency between each vehicle Wi and the server 20.
- the vehicles registered not only the vehicle W i for transmitting data of the feature quantity X ij, are also included vehicles are connected by a communication with the server 20.
- Road surface condition determination apparatus 30 includes a feature quantity calculation value calculating means 31, and the feature classifier value calculating means 32, a road surface condition judging means 33, the data of the feature quantity X ij transmitted from each vehicle W i It is used to determine the road surface condition within a predetermined time and within a predetermined range.
- the road surface state R is classified into a DRY road surface, a WET road surface, a SNOW road surface, and an ICE road surface.
- Feature calculation value calculating means 31 calculates a calculation value X j of the feature X ij of the vehicle W i at the location of the range and the predetermined given time.
- X j (X 1j + X 2j +... + X nj ) / n, that is, the average value of the feature amount X ij
- the number of operation values Xj is m, which is the same as the number of feature values Xij .
- the calculated value may be the median instead of the average value.
- the feature quantity identification value calculation means 32 inputs the calculated value X j to a road surface identification function F R (R; road surface state) set in advance for each road surface state, and is a feature that is an identification value for each feature amount X j.
- a quantity identification value Y R is calculated.
- the road surface identification function Y R (X) separates the road surface state R and the other road surface state R ′, and is created for each road surface state.
- the coefficient a R is a vector having m components like the calculation value X j, and b R is a constant.
- a RK is obtained in advance by running the test vehicle on the road surfaces of DRY, WET, SNOW, and ICE.
- the feature amount identification value Y D on the DRY road surface can be expressed by the following equation (1).
- Y D a D1 ⁇ X 1 + a D2 ⁇ X 2 + whil + a Dm ⁇ X m + b D (1)
- the feature quantity identification value Y W on the WET road surface, the feature quantity identification value Y S on the SNOW road surface, and the feature quantity identification value Y I on the ICE road surface are respectively expressed by the following equations (2) to (4). I can express.
- Y W a W1 ⁇ X 1 + a W2 ⁇ X 2 + whil + a Wm ⁇ X m + b W (2)
- Y S a S1 ⁇ X 1 + a S2 ⁇ X 2 + whil + a Sm ⁇ X m + b S (3)
- Y I a I1 ⁇ X 1 + a I2 ⁇ X 2 + whil + a Im ⁇ X m + b I (4)
- the road surface state determination means 33 determines the road surface state R corresponding to the feature value identification value Y R taking the maximum value as the road surface state at a place within a predetermined time and within a predetermined range.
- the road surface identification function Y R (X) a known identification method such as a kernel method, SVM, least square method, or neural network may be used.
- the calculated value X j of the feature value X ij of the vehicle W i at a place within a predetermined time and within a predetermined range is calculated, and the calculated calculated value X j is calculated.
- each is input to a road surface identification function Y R (X) set in advance for each road surface state, and feature value identification values Y D , Y W , Y S , Y I for each road surface state R are calculated. Since the road surface state R corresponding to the feature amount identification value Y R to be taken is the road surface state at a place within a predetermined time and within a predetermined range, the road surface state can be accurately determined.
- Embodiment 2 the feature value identification value Y R obtained by inputting the calculated value X j of the feature value X ij of the vehicle W i to the road surface identification function Y R (X) set in advance for each road surface state.
- Y W, Y S has been determined road surface state R with Y I, instead of the road surface condition determination unit 30, provided with a configuration of a road surface state detecting apparatus 30A shown in FIG. 4, wherein for each vehicle W i
- the quantity identification value Y R (X ij ) may be obtained, and the road surface state R may be determined using these feature quantity identification values Y R (X ij ).
- the road surface state determination device 30 ⁇ / b> A includes an identification function value calculation unit 34, an identification function calculation value calculation unit 35, and a road surface state determination unit 33.
- the discriminant function value Y iR F R (X ij ) of the feature value X ij of the vehicle W i at the location is calculated.
- Feature classifier value Y iD in DRY road surface the vehicle W i can be represented by the following equation (5).
- Y iD a D1 ⁇ X i1 + a D2 ⁇ X i2 + whil + a Dm ⁇ X im + b D (5)
- the feature quantity identification value Y iW on the WET road surface, the feature quantity identification value Y iS on the SNOW road surface, and the feature quantity identification value Y iI on the ICE road surface are respectively expressed by the following equations (6) to (8). I can express.
- Y iW a W1 ⁇ X i1 + a W2 ⁇ X i2 + whil + a Wm ⁇ X im + b W (6)
- Y iS a S1 ⁇ X i1 + a S2 ⁇ X i2 + whil + a Sm ⁇ X im + b S (7)
- Y iI a I1 ⁇ X i1 + a I2 ⁇ X i2 + whereas + a Im ⁇ X im + b I (8)
- the discriminant function calculation value Z R which is the calculated value is calculated.
- Z R (Y 1R + Y 2R +... + Y nR ) / n, that is, the average value of the discrimination function values Y iR of the respective vehicles W i can be used.
- the discriminant function calculation value Z R may be a median value instead of the average value.
- the road surface state determination means 33 determines the road surface state R corresponding to the discriminant function calculation value Z R taking the maximum value as the road surface state in a place within a predetermined time and within a predetermined range. If the determination of the road surface state is performed for each region including a plurality of preset locations at predetermined times, the road surface state at each time on the road map can be accurately determined.
- Embodiment 3 Wherein in the second embodiment, calculated after obtaining the n base vehicle W i of the feature quantity X ij classification function value Y iR respectively, the identification function computation value Z R is a calculated value using the discriminant function value Y iR Then, the road surface state is discriminated using this discriminant function calculation value Z R , but the discriminant function value distribution state of the discriminant function value Y iR and the distribution state of discriminant function value C R for each road surface state obtained in advance. a similarity between calculated identification value similarity H R, from the identification value similarity H R, if so to determine the road surface condition at the location of the range and the predetermined given time, the road surface condition The discrimination accuracy can be further improved.
- the road surface state determination system 3 includes an acceleration sensor 11, a wheel speed sensor 12, a feature amount calculation unit 13, and the like.
- N vehicles W i (i 1 to n) having a GPS device 14, vehicle information collecting means 15, a transmitter 16, and a receiver 17, a server 20 provided in the road surface state management center 2, and And a road surface state determination device 30B.
- the road surface state determination device 30B will be described.
- the road surface state determination device 30 ⁇ / b> B includes an identification function value calculation unit 34, an identification value similarity calculation unit 36, and a road surface state determination unit 33.
- the identification function value calculation unit 34 using the classification function value calculating means 34 in the second embodiment, the identification function of the feature quantity X ij of the vehicle W i at the location of the range and the predetermined given time
- the calculation formula of the feature amount identification value Y iR on the road surface R in the vehicle W i is as follows.
- the identification value similarity H R can be obtained by adding the similarity h iR for each vehicle, which is a value obtained by substituting the identification function value Y i of each vehicle W i into the variable x in the above formula [Equation 1]. it can.
- the similarity h iR for each vehicle is a value corresponding to the distance between the identification road surface model M R and the identification function value vector Y i of each vehicle W i , and the identification value similarity H R is the sum thereof.
- Road surface condition judging means 33 a road surface state R corresponding to the identification value similarity H R taking the maximum value, it is determined that the road surface condition at the location of the range and the predetermined given time.
- the road surface state R is obtained by examining which road surface model in the road surface model M R constituting the mixed Gaussian model G R (x) is similar to the distribution state of the discriminant function value vector Y i. Since the determination is made, it is possible to determine the road surface condition R at a place within a predetermined time and within a predetermined range with higher accuracy.
- Embodiment 4 it has been determined road surface state R from the distribution of the classification function value Y i of the distribution and identification road model M R of the vehicle W i, in place of the road surface condition determination unit 30B, FIG. providing the structure of a road surface state detecting apparatus 30C shown in 7, wherein the amount of the road surface model which represents the distribution of the feature vectors P i constituting the feature quantity X ij of each vehicle W i, a distribution of the feature vectors of each road surface it may be determined road surface state R from the distribution of the C R.
- the road surface state determination device 30 ⁇ / b> C includes a feature amount similarity calculation unit 37 and a road surface state determination unit 33.
- the feature quantity road surface model C R is 4 (number of road surface states) ⁇ N feature quantity vectors P i obtained by running a test vehicle equipped with a standard tire many times (N times) at a constant speed on the R road surface. Is the distribution state.
- the degree of similarity between the distribution state of the feature vector P i and the feature road model C R is a mixed Gaussian model G R (x) as described in the above equation [Equation 1]. It can be calculated using a statistical model such as
- the similarity S R is also larger as the distance between the feature amount road surface model C R and the feature amount X ij of each vehicle W i is shorter in the feature space.
- the road surface state determination means 33 determines the road surface state R corresponding to the similarity S R having the maximum value as the road surface state at a place within a predetermined time and within a predetermined range.
- the first to fourth embodiments there are four types of road surface conditions: DRY road surface, WET road surface, SNOW road surface, and ICE road surface. It may be divided into “ ⁇ path ( ⁇ ⁇ 0.7)”, “medium ⁇ path (0.3 ⁇ ⁇ 0.7)”, “low ⁇ path ( ⁇ ⁇ 0.3)”, and the like. Further, in the first to fourth embodiments have been judged road surface condition from the collected vehicle information for each vehicle W i, to each vehicle W i, Web information acquisition for acquiring outside information and connect to an unillustrated Internet Provide means, and acquire weather information such as weather, temperature, precipitation, wind speed, sunshine hours and road information such as terrain information, traffic information, road management information, etc.
- weather information such as weather, temperature, precipitation, wind speed, sunshine hours
- road information such as terrain information, traffic information, road management information, etc.
- the road surface identification function F R (X), the identification road surface model M R , or the feature amount road surface model C R is corrected, or the road surface identification function F R (X), the identification road surface model is obtained for each vehicle outside information. If M R or the feature amount road surface model CR is set or created, the road surface state determination accuracy can be further improved.
- Web information acquisition means may be provided in the server 20, and the server 20 may acquire weather information and road information.
- a road surface state determination system includes an in-vehicle sensor that acquires vehicle information that is information on behavior of a vehicle during travel, and a road surface state (hereinafter referred to as the vehicle surface that the vehicle travels from vehicle information acquired by the in-vehicle sensor).
- the road state determination device uses the feature amounts of the plurality of vehicles, characterized in that to determine the road surface condition at the location of the range and the predetermined said predetermined time.
- the road surface state is determined using a plurality of feature amounts obtained from the vehicle information of the plurality of vehicles, the road surface state determination accuracy can be improved.
- the “feature value” means the magnitude (band value) of vibration in a plurality of specific frequency bands of the time series waveform of the acceleration sensor provided in the tire, wheel, or knuckle, or provided in the tire chamber.
- a physical quantity used for estimation of road surface conditions such as a band value of an air pressure fluctuation signal detected by a detected pressure sensor or a band value of a sound pressure signal detected by a microphone provided in a tire or a front wheel front frame.
- the “predetermined time” can be changed, for example, by setting it longer or shorter according to the section in which the vehicle travels and the time at which the vehicle travels.
- the “predetermined time” may be a time required for a predetermined number (for example, five) of vehicles to pass through a place within a predetermined range. Also in this case, the number of vehicles set in advance can be changed according to the traveling section and the traveling time.
- the vehicle-mounted sensor is an acceleration sensor disposed in the tire, the difference in tire vibration due to the road surface condition can be accurately grasped, and the road surface condition determination accuracy can be further improved.
- the road surface state determination device has a feature of the plurality of vehicles at a place within a predetermined time and within a predetermined range for each of the plurality of feature amounts.
- a feature classifier value calculating means for calculating an amount identification value Y R for each said road identification function, using the feature classifier value Y R for each road identification function calculated, the predetermined time and within a predetermined range Road surface state determining means for determining the road surface state R at a place inside It is characterized by.
- the discriminant value the road surface state is discriminated using all the data of the feature amounts X ij of a plurality of vehicles, so that the road surface state discriminating accuracy can be improved.
- Discrimination function value calculation means for calculating for each discrimination function; discrimination function calculation value calculation means for calculating an identification function calculation value Z R that is a calculation value using the calculated discrimination function value Y iR for each road surface discrimination function; And a road surface state determining means for determining a road surface state at a location within the predetermined time and within a predetermined range using the calculated discriminant function calculation value Z R.
- the discrimination function value Y iR is obtained for each vehicle, and the discrimination function calculation value Z R that is the calculation value of these discrimination function values Y iR is used as the discrimination value of the road surface state R, thereby improving the discrimination accuracy of the road surface state. Can be made.
- the road surface state determination device has the road surface identification function Y in which the plurality of feature amounts X ij are set in advance for each road surface state for each of the plurality of vehicles.
- R F R (X) is input respectively, and an identification function value Y iR that is an identification value for each vehicle is calculated for each road surface identification function, and for each calculated road surface identification function of the distribution of classification function value Y iR, advance distribution of discriminant function values for each road surface condition R which has been determined (hereinafter, referred to as identification road model M R) which is a similarity identification value similarity H R and calculated, from said calculated identification value similarity H R, characterized in that it comprises a road surface condition judging means for judging road surface condition at the location of the range and the predetermined said predetermined time.
- the classification function value Y iR of distribution obtained from the feature X ij is, for example, a vehicle on DRY road
- the identified road surface model M R ′ that is close to the identified road surface model M D obtained from the feature value obtained when traveling or is obtained from the feature value obtained when traveling on another road surface R ′. Since the determination criterion for road surface discrimination is used as the determination criterion, the road surface state determination accuracy can be further improved.
- the road surface state determination device determines the distribution state of the feature amount X ij for each of the plurality of vehicles for each road surface state.
- feature amount of distribution hereinafter
- feature amount similarity is similarity between S R and the feature amount similarity calculating means for calculating the feature quantity similarity S of the vehicle as the calculated From R
- road surface state determining means for determining a road surface state at a place within a predetermined time and within a predetermined range is provided.
- the feature amount road surface model CR is created for each road surface state in advance, and the distribution state of the actually measured feature amount X ij is obtained, for example, when the vehicle travels on the DRY road surface. was close to or characteristic quantity road model C D obtained from the feature, or the road surface or near the 'feature quantity road model C R calculated from the obtained feature amount when is traveling on' another road R.
- the discrimination accuracy of the road surface condition can be further improved as a discrimination criterion.
- the road surface state determination system is provided with a means for acquiring vehicle outside information, and using the acquired vehicle outside information, the road surface identification function F R (X) or the road surface. It is characterized in that the distribution state of the identification function value for each state (identification road surface model M R ) or the distribution of the feature amount for each road surface state (feature amount road surface model C R ) is corrected.
- the information outside the vehicle refers to weather information such as weather, temperature, precipitation, wind speed, and sunshine hours acquired through the Internet, or road information such as topographic information, traffic information, and road management information. Further, information used for the correction includes an estimation result of a road surface state at a location within another predetermined range.
- the road surface state estimation results of the other location B within the predetermined range correlated with the location A within the predetermined range are used.
- the estimation result at the location A within the predetermined range can be corrected.
- the road surface state in a place within a predetermined time and within a predetermined range can be determined with higher accuracy.
- the present invention it is possible to accurately determine the road surface state in a place within the predetermined time and within a predetermined range. If the vehicle to be notified is notified, the traveling safety of the vehicle can be improved.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Environmental & Geological Engineering (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Transportation (AREA)
- Biodiversity & Conservation Biology (AREA)
- Environmental Sciences (AREA)
- Ecology (AREA)
- Atmospheric Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Automation & Control Theory (AREA)
- Mathematical Physics (AREA)
- Traffic Control Systems (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Tires In General (AREA)
Abstract
路面状態の判別精度を向上させた路面状態判別システムを提供するために、車両の走行時の挙動の情報である車両情報を取得する車載センサ11,12、車載センサ11,12により取得した車両情報から車両の走行する路面の状態を推定するための複数の特徴量を算出する特徴量算出手段13、及び、特徴量を当該車両の外部に送信する送信機16を備えた複数台の車両Wiと、複数台の車両の複数の特徴量を蓄積するデータ記憶手段22を有するサーバー20と、蓄積された特徴量を用いて路面状態の判別を行う路面状態判別装置30とを備え、所定の時間内かつ所定の範囲内の場所における路面状態を判別する。
Description
本発明は、走行中の路面の状態を判別するシステムに関するもので、特に、所定の時間内に同じ箇所を走行する複数台の車両の車両情報を用いて、路面状態を判別する路面状態判別システムに関する。
自動車の走行安全性を高めるため、走行中の路面の状態を精度良く推定し、車両制御へフィードバックすることが求められている。走行中の路面の状態を推定することができれば、制駆動や操舵といった危険回避の操作を起こす前に、例えば、ABSブレーキのより高度な制御等が可能になり、安全性が一段と高まることが予想される。
走行中の路面の状態を推定する方法としては、例えば、走行中のタイヤの振動を検出し、この検出されたタイヤの振動の時系列波形から路面状態を推定する方法(例えば、特許文献1~3参照)や、タイヤから発生するタイヤ発生音を検出し、この検出されたタイヤ発生音の音圧レベルから路面状態を推定する方法が提案されている(例えば、特許文献4参照)。
しかしながら、従来は、1台の車両の情報から走行中の路面状態を推定しているため、路面状態の判別精度が必ずしも十分とはいえなかった。
本発明は、従来の問題点に鑑みてなされたもので、路面状態の判別精度を向上させることのできる路面状態判別システムを提供することを目的とする。
本発明の一態様にかかる路面状態判別システムは、車両の走行時の挙動の情報である車両情報を取得する車載センサ、前記車載センサにより取得した車両情報から前記車両の走行する路面の状態(以下、路面状態という)を推定するための複数の特徴量を算出する特徴量算出手段、及び、前記特徴量を当該車両の外部に送信する送信手段を備えた複数台の車両と、前記複数台の車両Wiの送信手段からそれぞれ送信される特徴量を蓄積するサーバーと、前記蓄積された複数台の車両の複数の特徴量Xij(i=1~n:nは車両の数、j=1~m:mは特徴量の数)を用いて路面状態の判別を行う路面状態判別装置とを備えた路面状態判別システムであって、前記送信される特徴量が、所定の時間内かつ所定の範囲内の場所における前記特徴量であり、前記路面状態判別装置は、前記複数台の車両の特徴量を用いて、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別することを特徴とする。
実施の形態1.
図1は、本実施の形態1に係る路面状態判別システム1の構成を示す機能ブロック図で、W1~Wnは車両で、各車両Wi(i=1~n)は、加速度センサ11と、車輪速センサ12と、特徴量算出手段13と、GPS装置14と、車両情報収集手段15と、送信機16と、受信機17とを備える。
なお、車両W1~Wnは、所定の時間内に所定の範囲内の場所を通過する車両である。
20は受信機21と、データ記憶手段22と、送信機23とを備えたサーバーで、30は所定の時間内かつ所定の範囲内の場所における路面の状態を判別する路面状態判別装置である。
サーバー20と路面状態判別装置30とは、路面状態管理センター2に設けられる。
特徴量算出手段13、及び、路面状態判別装置30は、例えば、コンピュータのソフトウェアにより構成される。
加速度センサ11は、図2に示すように、タイヤ40のインナーライナー部41のタイヤ気室42側のほぼ中央部に配置されて、路面からタイヤ40のトレッド43に入力する振動(タイヤ振動)を検出する。
車輪速センサ12は、車輪の回転速度(車輪速)を検出するもので、例えば、車輪の回転角度を検出する周知の電磁誘導型の車輪速センサを用いることができる。
図1は、本実施の形態1に係る路面状態判別システム1の構成を示す機能ブロック図で、W1~Wnは車両で、各車両Wi(i=1~n)は、加速度センサ11と、車輪速センサ12と、特徴量算出手段13と、GPS装置14と、車両情報収集手段15と、送信機16と、受信機17とを備える。
なお、車両W1~Wnは、所定の時間内に所定の範囲内の場所を通過する車両である。
20は受信機21と、データ記憶手段22と、送信機23とを備えたサーバーで、30は所定の時間内かつ所定の範囲内の場所における路面の状態を判別する路面状態判別装置である。
サーバー20と路面状態判別装置30とは、路面状態管理センター2に設けられる。
特徴量算出手段13、及び、路面状態判別装置30は、例えば、コンピュータのソフトウェアにより構成される。
加速度センサ11は、図2に示すように、タイヤ40のインナーライナー部41のタイヤ気室42側のほぼ中央部に配置されて、路面からタイヤ40のトレッド43に入力する振動(タイヤ振動)を検出する。
車輪速センサ12は、車輪の回転速度(車輪速)を検出するもので、例えば、車輪の回転角度を検出する周知の電磁誘導型の車輪速センサを用いることができる。
特徴量算出手段13は、振動波形検出部13aと、領域信号抽出部13bと、特徴量算出部13cとを備え、加速度センサ11で検出したタイヤ振動の時系列波形から、当該車両Wiの走行する路面の状態(路面状態)を推定するための特徴量を算出する。
振動波形検出部13aは、図3に示すように、車輪速センサ12で検出した車輪速度を用いて、タイヤ一回転分の加速度波形を検出する。
領域信号抽出部13bは、タイヤ一回転分の加速度波形を、踏み込み端に出現する踏み込み側のピークよりも前の踏み込み前領域R1と、前記踏み込み側のピークを形成する踏み込み領域R2と、前記踏み込み側のピークと蹴り出し端に出現する蹴り出し側のピークとの間の蹴り出し前領域R3と、前記蹴り出し側のピークを形成する蹴り出し領域R4と、前記蹴り出し側のピークよりも後の蹴り出し後領域R5とに分割し、各領域R1~R5における振動レベルの時系列波形を抽出する。
特徴量算出部13cは、各領域R1~R5のそれぞれの時系列波形をバンドパスフィルターに通して、所定の領域で、かつ、所定の周波数領域の振動成分の大きさである帯域値をm個算出し、これらの帯域値を車両Wi(i=1~n)の特徴量Xij(j=1~m)として、送信機16から路面状態管理センター2に送信する。なお、各領域R1~R5からの帯域値は複数であってもよいし、特徴量Xijを算出しない領域があってもよい。また、タイヤ一回転分の加速度波形を、踏み込み前領域、接地領域、蹴り出し後領域の3つに分けるなどしてもよい。
換言すれば、車両Wi(i=1~n)からは、成分が(Xi1,Xi2,……,Xim)であるm次元の特徴ベクトルXiが路面状態管理センター2に送信される。
振動波形検出部13aは、図3に示すように、車輪速センサ12で検出した車輪速度を用いて、タイヤ一回転分の加速度波形を検出する。
領域信号抽出部13bは、タイヤ一回転分の加速度波形を、踏み込み端に出現する踏み込み側のピークよりも前の踏み込み前領域R1と、前記踏み込み側のピークを形成する踏み込み領域R2と、前記踏み込み側のピークと蹴り出し端に出現する蹴り出し側のピークとの間の蹴り出し前領域R3と、前記蹴り出し側のピークを形成する蹴り出し領域R4と、前記蹴り出し側のピークよりも後の蹴り出し後領域R5とに分割し、各領域R1~R5における振動レベルの時系列波形を抽出する。
特徴量算出部13cは、各領域R1~R5のそれぞれの時系列波形をバンドパスフィルターに通して、所定の領域で、かつ、所定の周波数領域の振動成分の大きさである帯域値をm個算出し、これらの帯域値を車両Wi(i=1~n)の特徴量Xij(j=1~m)として、送信機16から路面状態管理センター2に送信する。なお、各領域R1~R5からの帯域値は複数であってもよいし、特徴量Xijを算出しない領域があってもよい。また、タイヤ一回転分の加速度波形を、踏み込み前領域、接地領域、蹴り出し後領域の3つに分けるなどしてもよい。
換言すれば、車両Wi(i=1~n)からは、成分が(Xi1,Xi2,……,Xim)であるm次元の特徴ベクトルXiが路面状態管理センター2に送信される。
なお、車両Wiの特徴量Xijとしては、タイヤ内,ホイール、もしくは、ナックルに設けられた加速度センサにより取得した加速度の時系列波形の周波数スペクトルから得られた、路面状態により変化する複数の周波数帯域の振動レベルや、タイヤ気室内に設けられた圧力センサで検出した空気圧の変動スペクトルの帯域値、あるいは、タイヤ内や後輪前方フレームに設けられたマイクロフォンで採取したタイヤ発生音のバンドパワー値などを挙げることができる。
あるいは、車載センサとして、二方向振動センサを設けてタイヤ周方向の振動とタイヤ径方向の振動とを検出し、これらの振動波形の特定周波数成分を車両Wiの特徴量Xijとしてもよい。この場合には、特徴ベクトルXiの成分は2つである。
更に、特徴量としては、トレッドブロック内に配置したトレッド剪断応力や、タイヤ内面に設けられた歪センサで検出した周方向タイヤ内面歪、車輪速センサで検出した車輪速の時系列データなどを用いることができる。
あるいは、車載センサとして、二方向振動センサを設けてタイヤ周方向の振動とタイヤ径方向の振動とを検出し、これらの振動波形の特定周波数成分を車両Wiの特徴量Xijとしてもよい。この場合には、特徴ベクトルXiの成分は2つである。
更に、特徴量としては、トレッドブロック内に配置したトレッド剪断応力や、タイヤ内面に設けられた歪センサで検出した周方向タイヤ内面歪、車輪速センサで検出した車輪速の時系列データなどを用いることができる。
GPS装置14は、図示しないGPSアンテナと受信機とを備え、当該車両Wiの位置データを取得する。
車両情報収集手段15は、特徴量算出手段13で算出された車両Wiの特徴量Xijと、GPS装置14で取得した車両の位置データと、車両情報の取得時刻のデータ(時刻データ)とを収集する。なお、時刻データとしては、加速度の時系列波形の抽出時刻,位置データの取得時刻、あるいは、データ送信時刻などを用いればよい。これらの時刻は、ほぼ同時刻なので、どの時刻を時刻データとしても問題はない。
送信機16は、車両情報収集手段15により収集されたm個の特徴量Xijのデータと、車両の位置データと、時刻データとを、当該車両を識別するための車両IDとともに、図示しない送信アンテナから、路面状態管理センター2のサーバー20に送信する。
受信機17は、路面状態管理センター2の路面状態判別装置30で判別され、サーバー20の送信機23から送信された所定の時間内かつ所定の範囲内の場所における路面状態のデータを受信する。なお、この受信された路面状態のデータを、車内に設けられたモニター上に表示することで、運転者は、前記所定の時間内かつ所定の範囲内の場所における路面状態を知ることができる。
また、路面状態のデータを車両制御へフィードバックすれば、車両の走行安全性を向上させることができる。
車両情報収集手段15は、特徴量算出手段13で算出された車両Wiの特徴量Xijと、GPS装置14で取得した車両の位置データと、車両情報の取得時刻のデータ(時刻データ)とを収集する。なお、時刻データとしては、加速度の時系列波形の抽出時刻,位置データの取得時刻、あるいは、データ送信時刻などを用いればよい。これらの時刻は、ほぼ同時刻なので、どの時刻を時刻データとしても問題はない。
送信機16は、車両情報収集手段15により収集されたm個の特徴量Xijのデータと、車両の位置データと、時刻データとを、当該車両を識別するための車両IDとともに、図示しない送信アンテナから、路面状態管理センター2のサーバー20に送信する。
受信機17は、路面状態管理センター2の路面状態判別装置30で判別され、サーバー20の送信機23から送信された所定の時間内かつ所定の範囲内の場所における路面状態のデータを受信する。なお、この受信された路面状態のデータを、車内に設けられたモニター上に表示することで、運転者は、前記所定の時間内かつ所定の範囲内の場所における路面状態を知ることができる。
また、路面状態のデータを車両制御へフィードバックすれば、車両の走行安全性を向上させることができる。
サーバー20は、各車両Wi(i=1~n)から送られてくる車両情報を受信機21で受信し、これらのデータを、データ記憶手段22に保存するとともに、送信機23から、路面状態判別装置30で判別された所定の時間内かつ所定の範囲内の場所における路面状態のデータを登録車両に送信する。
各車両Wiの特徴量Xijのデータは、データ記憶手段22に分類・保存される。
ここで、所定の時間は、予め設定された判別時刻を含む1分~5分の範囲を指し、所定の範囲とは、例えば、予め設定された道路地図上の場所を含む範囲を指す。また、所定の範囲の指定は、例えば、ナビゲーション用の道路地図上の所定の大きさ枠のグリットを用いれば、各車両Wiとサーバー20とでの位置の整合性を取り易いので好ましい。
なお、登録車両とは、特徴量Xijのデータを送信する車両Wiだけではなく、サーバー20と通信により接続されている車両も含まれる。
各車両Wiの特徴量Xijのデータは、データ記憶手段22に分類・保存される。
ここで、所定の時間は、予め設定された判別時刻を含む1分~5分の範囲を指し、所定の範囲とは、例えば、予め設定された道路地図上の場所を含む範囲を指す。また、所定の範囲の指定は、例えば、ナビゲーション用の道路地図上の所定の大きさ枠のグリットを用いれば、各車両Wiとサーバー20とでの位置の整合性を取り易いので好ましい。
なお、登録車両とは、特徴量Xijのデータを送信する車両Wiだけではなく、サーバー20と通信により接続されている車両も含まれる。
路面状態判別装置30は、特徴量演算値算出手段31と、特徴量識別値算出手段32と、路面状態判別手段33とを備え、各車両Wiから送られてきた特徴量Xijのデータを用いて、所定の時間内かつ所定の範囲内の場所における路面状態の判別を行う。なお、本例では、路面状態Rを、DRY路面,WET路面,SNOW路面,ICE路面に分類した。
特徴量演算値算出手段31は、所定の時間内かつ所定の範囲内の場所における車両Wiの特徴量Xijの演算値Xjを算出する。
演算値Xjとしては、例えば、Xj=(X1j+X2j+……+Xnj)/n、すなわち、特徴量Xijの平均値を用いることができる。演算値Xjの個数はm個で、特徴量Xijの個数と同じである。なお、演算値としては、平均値に代えて、中央値としてもよい。
特徴量識別値算出手段32は、演算値Xjを、予め路面状態毎に設定した路面識別関数FR(R;路面状態)にそれぞれ入力して、特徴量Xj毎の識別値である特徴量識別値YRを算出する。路面識別関数YR(X)は、路面状態Rと他の路面状態R’とを分けるもので、路面状態毎に作成される。
本例では、路面識別関数FR(X)を、YR=FR(X)=aR・X+bRとした。
係数aRは、演算値Xjと同じくm個の成分を持つベクトルで、bRは定数である。ここで、aR=(aR1,aR2,……,aRm)とすると、ベクトルの各成分aRK(k=1~m)は、特徴量XKに対する「重み」で、路面状態Rにより異なる値をとる。aRKは、予め、試験車両を,DRY,WET,SNOW,ICEの各路面上を走行させて求められる。
DRY路面における特徴量識別値YDは、以下の式(1)で表せる。
YD=aD1・X1+aD2・X2+……+aDm・Xm +bD……(1)
同様に、WET路面における特徴量識別値YW、SNOW路面における特徴量識別値YS、及び、ICE路面における特徴量識別値YIは、それぞれ以下の式(2)~(4)のように表せる。
YW=aW1・X1+aW2・X2+……+aWm・Xm +bW……(2)
YS=aS1・X1+aS2・X2+……+aSm・Xm +bS……(3)
YI=aI1・X1+aI2・X2+……+aIm・Xm +bI……(4)
路面状態判別手段33は、最大値を取る特徴量識別値YRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
この路面状態の判別を、予め設定された時刻毎に、予め設定された複数の箇所を含む領域毎に行えば、道路地図上における各時刻における路面状態を精度よく推定することができる。
なお、路面識別関数YR(X)としては、カーネル法、SVM、最小二乗法、ニューラルネットなどの既知の識別手法を用いてもよい。
特徴量演算値算出手段31は、所定の時間内かつ所定の範囲内の場所における車両Wiの特徴量Xijの演算値Xjを算出する。
演算値Xjとしては、例えば、Xj=(X1j+X2j+……+Xnj)/n、すなわち、特徴量Xijの平均値を用いることができる。演算値Xjの個数はm個で、特徴量Xijの個数と同じである。なお、演算値としては、平均値に代えて、中央値としてもよい。
特徴量識別値算出手段32は、演算値Xjを、予め路面状態毎に設定した路面識別関数FR(R;路面状態)にそれぞれ入力して、特徴量Xj毎の識別値である特徴量識別値YRを算出する。路面識別関数YR(X)は、路面状態Rと他の路面状態R’とを分けるもので、路面状態毎に作成される。
本例では、路面識別関数FR(X)を、YR=FR(X)=aR・X+bRとした。
係数aRは、演算値Xjと同じくm個の成分を持つベクトルで、bRは定数である。ここで、aR=(aR1,aR2,……,aRm)とすると、ベクトルの各成分aRK(k=1~m)は、特徴量XKに対する「重み」で、路面状態Rにより異なる値をとる。aRKは、予め、試験車両を,DRY,WET,SNOW,ICEの各路面上を走行させて求められる。
DRY路面における特徴量識別値YDは、以下の式(1)で表せる。
YD=aD1・X1+aD2・X2+……+aDm・Xm +bD……(1)
同様に、WET路面における特徴量識別値YW、SNOW路面における特徴量識別値YS、及び、ICE路面における特徴量識別値YIは、それぞれ以下の式(2)~(4)のように表せる。
YW=aW1・X1+aW2・X2+……+aWm・Xm +bW……(2)
YS=aS1・X1+aS2・X2+……+aSm・Xm +bS……(3)
YI=aI1・X1+aI2・X2+……+aIm・Xm +bI……(4)
路面状態判別手段33は、最大値を取る特徴量識別値YRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
この路面状態の判別を、予め設定された時刻毎に、予め設定された複数の箇所を含む領域毎に行えば、道路地図上における各時刻における路面状態を精度よく推定することができる。
なお、路面識別関数YR(X)としては、カーネル法、SVM、最小二乗法、ニューラルネットなどの既知の識別手法を用いてもよい。
このように、本実施の形態1では、所定の時間内かつ所定の範囲内の場所における車両Wiの特徴量Xijの演算値Xjをそれぞれ算出し、この算出された演算値Xjを、予め路面状態毎に設定した路面識別関数YR(X)にそれぞれ入力して、路面状態R毎の特徴量識別値YD,YW,YS,YIを算出した後、最大値をとる特徴量識別値YRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態としたので、路面状態を精度よく判別することができる。
実施の形態2.
前記実施の形態1では、車両Wiの特徴量Xijの演算値Xjを予め路面状態毎に設定した路面識別関数YR(X)にそれぞれ入力して得られた特徴量識別値YR,YW,YS,YIを用いて路面状態Rを判定したが、路面状態判別装置30に代えて、図4に示す構成の路面状態判別装置30Aを設けて、車両Wi毎に特徴量識別値YR(Xij)を求め、これらの特徴量識別値YR(Xij)を用いて路面状態Rを判定するようにしてもよい。
路面状態判別装置30Aは、識別関数値算出手段34と、識別関数演算値算出手段35と、路面状態判別手段33とを備える。
識別関数値算出手段34は、実施の形態1と同様の、路面状態毎に設定された路面識別関数FR(X)=aR・X+bRを用いて、所定の時間内かつ所定の範囲内の場所における車両Wiの特徴量Xijの識別関数値YiR=FR(Xij)を算出する。
車両WiにおけるDRY路面における特徴量識別値YiDは、以下の式(5)で表せる。
YiD=aD1・Xi1+aD2・Xi2+……+aDm・Xim +bD……(5)
同様に、WET路面における特徴量識別値YiW、SNOW路面における特徴量識別値YiS、及び、ICE路面における特徴量識別値YiIは、それぞれ以下の式(6)~(8)のように表せる。
YiW=aW1・Xi1+aW2・Xi2+……+aWm・Xim +bW……(6)
YiS=aS1・Xi1+aS2・Xi2+……+aSm・Xim +bS……(7)
YiI=aI1・Xi1+aI2・Xi2+……+aIm・Xim +bI……(8)
識別関数演算値算出手段35は、識別関数値算出手段34で算出されたn台の車両Wiの識別関数値YiR(i=1~n、R;DRY,WET,SNOW,ICE)を用いた演算値である識別関数演算値ZRを算出する。
識別関数演算値ZRとしては、例えば、ZR=(Y1R+Y2R+……+YnR)/n、すなわち、各車両Wiの識別関数値YiRの平均値を用いることができる。
なお、識別関数演算値ZRとしては、平均値に代えて、中央値としてもよい。
路面状態判別手段33は、最大値を取る識別関数演算値ZRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
この路面状態の判別を、予め設定された時刻毎に、予め設定された複数の箇所を含む領域毎に行えば、道路地図上における各時刻における路面状態を精度よく判別することができる。
前記実施の形態1では、車両Wiの特徴量Xijの演算値Xjを予め路面状態毎に設定した路面識別関数YR(X)にそれぞれ入力して得られた特徴量識別値YR,YW,YS,YIを用いて路面状態Rを判定したが、路面状態判別装置30に代えて、図4に示す構成の路面状態判別装置30Aを設けて、車両Wi毎に特徴量識別値YR(Xij)を求め、これらの特徴量識別値YR(Xij)を用いて路面状態Rを判定するようにしてもよい。
路面状態判別装置30Aは、識別関数値算出手段34と、識別関数演算値算出手段35と、路面状態判別手段33とを備える。
識別関数値算出手段34は、実施の形態1と同様の、路面状態毎に設定された路面識別関数FR(X)=aR・X+bRを用いて、所定の時間内かつ所定の範囲内の場所における車両Wiの特徴量Xijの識別関数値YiR=FR(Xij)を算出する。
車両WiにおけるDRY路面における特徴量識別値YiDは、以下の式(5)で表せる。
YiD=aD1・Xi1+aD2・Xi2+……+aDm・Xim +bD……(5)
同様に、WET路面における特徴量識別値YiW、SNOW路面における特徴量識別値YiS、及び、ICE路面における特徴量識別値YiIは、それぞれ以下の式(6)~(8)のように表せる。
YiW=aW1・Xi1+aW2・Xi2+……+aWm・Xim +bW……(6)
YiS=aS1・Xi1+aS2・Xi2+……+aSm・Xim +bS……(7)
YiI=aI1・Xi1+aI2・Xi2+……+aIm・Xim +bI……(8)
識別関数演算値算出手段35は、識別関数値算出手段34で算出されたn台の車両Wiの識別関数値YiR(i=1~n、R;DRY,WET,SNOW,ICE)を用いた演算値である識別関数演算値ZRを算出する。
識別関数演算値ZRとしては、例えば、ZR=(Y1R+Y2R+……+YnR)/n、すなわち、各車両Wiの識別関数値YiRの平均値を用いることができる。
なお、識別関数演算値ZRとしては、平均値に代えて、中央値としてもよい。
路面状態判別手段33は、最大値を取る識別関数演算値ZRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
この路面状態の判別を、予め設定された時刻毎に、予め設定された複数の箇所を含む領域毎に行えば、道路地図上における各時刻における路面状態を精度よく判別することができる。
実施の形態3.
前記実施の形態2では、n台の車両Wiの特徴量Xijの識別関数値YiRをそれぞれ求めた後、識別関数値YiRを用いた演算値である識別関数演算値ZRを算出し、この識別関数演算値ZRを用いて路面状態を判別したが、識別関数値YiRの識別関数値の分布状態と、予め求めておいた路面状態毎の識別関数値CRの分布状態との類似度である識別値類似度HRを計算し、この識別値類似度HRから、所定の時間内かつ所定の範囲内の場所における路面状態を判別するようにすれば、路面状態の判別精度を更に向上させることができる。
図5は、本実施の形態3に係る路面状態判別システム3の構成を示す機能ブロック図で、路面状態判別システム3は、加速度センサ11と、車輪速センサ12と、特徴量算出手段13と、GPS装置14と、車両情報収集手段15と、送信機16と、受信機17とを備えたn台の車両Wi(i=1~n)と、路面状態管理センター2に設けられるサーバー20及び路面状態判別装置30Bとを備える。
以下に、路面状態判別装置30Bについて説明する。なお、車両Wi、及び、サーバー20の構成及び動作については、前記実施の形態1と同様であるので、説明を省略する。
前記実施の形態2では、n台の車両Wiの特徴量Xijの識別関数値YiRをそれぞれ求めた後、識別関数値YiRを用いた演算値である識別関数演算値ZRを算出し、この識別関数演算値ZRを用いて路面状態を判別したが、識別関数値YiRの識別関数値の分布状態と、予め求めておいた路面状態毎の識別関数値CRの分布状態との類似度である識別値類似度HRを計算し、この識別値類似度HRから、所定の時間内かつ所定の範囲内の場所における路面状態を判別するようにすれば、路面状態の判別精度を更に向上させることができる。
図5は、本実施の形態3に係る路面状態判別システム3の構成を示す機能ブロック図で、路面状態判別システム3は、加速度センサ11と、車輪速センサ12と、特徴量算出手段13と、GPS装置14と、車両情報収集手段15と、送信機16と、受信機17とを備えたn台の車両Wi(i=1~n)と、路面状態管理センター2に設けられるサーバー20及び路面状態判別装置30Bとを備える。
以下に、路面状態判別装置30Bについて説明する。なお、車両Wi、及び、サーバー20の構成及び動作については、前記実施の形態1と同様であるので、説明を省略する。
路面状態判別装置30Bは、識別関数値算出手段34と、識別値類似度算出手段36と、路面状態判別手段33とを備える。
本例では、識別関数値算出手段34として、実施の形態2の識別関数値算出手段34を用いて、所定の時間内かつ所定の範囲内の場所における車両Wiの特徴量Xijの識別関数値YiR=FR(Xij)を算出した。
車両Wiにおける路面Rにおける特徴量識別値YiRの算出式は以下の通りである。
YiR=aR1・Xi1+aR2・Xi2+……+aRm・Xim+bR
識別値類似度算出手段36は、識別関数値算出手段34で算出されたn台の車両Wiの識別関数値YiR(i=1~n、R;DRY,WET,SNOW,ICE)の分布状態と、予め路面状態R毎に求めておいた識別関数の分布(以下、識別路面モデルMRという)との識別値類似度HRを算出する。
識別路面モデルMRは、標準タイヤを装着した試験車両をR路面にて一定速度で多数回(N回)走行させて得られた4(路面状態数)×N個の特徴量XRj(j=1~m)をそれぞれ路面識別関数FRに代入して求めた4(路面状態数)×N個の識別関数値YiR=FR(XRj)の分布状態である。
識別関数値算出手段34を用いて算出された特徴量識別値YiRの分布状態と、識別路面モデルMRとの類似度は、例えば、混合ガウスモデルなどの統計モデルを用いて算出することができる。
混合ガウスモデルGR(x)は路面状態R毎(ここでは、4つ)に作成されるもので、以下の式[数1]で表せる。すなわち、[数1]の式は、識別路面モデルMRを混合ガウスモデルとしたときの式である。
Kは混合数で、dは路面状態Rの数(ここでは、d=4)である。
また、混合比πR、混合要素の共分散行列ΣR、混合要素の平均ベクトルμRなどのパラメータは、EMアルゴリズムなどの既知の手法を用いて求めることができる。
本例では、識別関数値算出手段34として、実施の形態2の識別関数値算出手段34を用いて、所定の時間内かつ所定の範囲内の場所における車両Wiの特徴量Xijの識別関数値YiR=FR(Xij)を算出した。
車両Wiにおける路面Rにおける特徴量識別値YiRの算出式は以下の通りである。
YiR=aR1・Xi1+aR2・Xi2+……+aRm・Xim+bR
識別値類似度算出手段36は、識別関数値算出手段34で算出されたn台の車両Wiの識別関数値YiR(i=1~n、R;DRY,WET,SNOW,ICE)の分布状態と、予め路面状態R毎に求めておいた識別関数の分布(以下、識別路面モデルMRという)との識別値類似度HRを算出する。
識別路面モデルMRは、標準タイヤを装着した試験車両をR路面にて一定速度で多数回(N回)走行させて得られた4(路面状態数)×N個の特徴量XRj(j=1~m)をそれぞれ路面識別関数FRに代入して求めた4(路面状態数)×N個の識別関数値YiR=FR(XRj)の分布状態である。
識別関数値算出手段34を用いて算出された特徴量識別値YiRの分布状態と、識別路面モデルMRとの類似度は、例えば、混合ガウスモデルなどの統計モデルを用いて算出することができる。
混合ガウスモデルGR(x)は路面状態R毎(ここでは、4つ)に作成されるもので、以下の式[数1]で表せる。すなわち、[数1]の式は、識別路面モデルMRを混合ガウスモデルとしたときの式である。
また、混合比πR、混合要素の共分散行列ΣR、混合要素の平均ベクトルμRなどのパラメータは、EMアルゴリズムなどの既知の手法を用いて求めることができる。
識別値類似度HRは、上記の式[数1]の変数xに各車両Wiの識別関数値Yiを代入した値である車両毎の類似度hiRを加算することで求めることができる。
識別関数値Yiはベクトルで、Yi=(YiD,YiW,YiS,YiI)である。
すなわち、DRY路面の類似度HDは、HD=ΣiGD(Yi)により算出され、WET路面の類似度HWは、HW=ΣiGW(Yi)により算出される。また、SNOW路面の類似度HSは、HS=ΣiGS(Yi)により算出され、ICE路面の類似度HIは、HI=ΣiGI(Yi)により算出される(i=1~n)。
車両毎の類似度hiRは、識別路面モデルMRと各車両Wiの識別関数値ベクトルYiとの距離に対応する値で、識別値類似度HRはその総和である。
路面状態判別手段33は、最大値を取る識別値類似度HRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
このように、識別関数値ベクトルYiの分布状態が、混合ガウスモデルGR(x)を構成する路面モデルMRのなかのどの路面モデルと類似しているかを調べることで、路面状態Rを判別するようにしたので、所定の時間内かつ所定の範囲内の場所における路面状態Rを更に精度よく判別することができる。
識別関数値Yiはベクトルで、Yi=(YiD,YiW,YiS,YiI)である。
すなわち、DRY路面の類似度HDは、HD=ΣiGD(Yi)により算出され、WET路面の類似度HWは、HW=ΣiGW(Yi)により算出される。また、SNOW路面の類似度HSは、HS=ΣiGS(Yi)により算出され、ICE路面の類似度HIは、HI=ΣiGI(Yi)により算出される(i=1~n)。
車両毎の類似度hiRは、識別路面モデルMRと各車両Wiの識別関数値ベクトルYiとの距離に対応する値で、識別値類似度HRはその総和である。
路面状態判別手段33は、最大値を取る識別値類似度HRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
このように、識別関数値ベクトルYiの分布状態が、混合ガウスモデルGR(x)を構成する路面モデルMRのなかのどの路面モデルと類似しているかを調べることで、路面状態Rを判別するようにしたので、所定の時間内かつ所定の範囲内の場所における路面状態Rを更に精度よく判別することができる。
実施の形態4.
前記実施の形態3では、各車両Wiの識別関数値Yiの分布状態と識別路面モデルMRとの分布状態とから路面状態Rを判定したが、路面状態判別装置30Bに代えて、図7に示す構成の路面状態判別装置30Cを設けて、車両Wi毎の特徴量Xijを構成する特徴量ベクトルPiの分布状態と、路面ごとの特徴量ベクトルの分布を表す特徴量路面モデルCRとの分布状態とから路面状態Rを判定するようにしてもよい。
路面状態判別装置30Cは、特徴量類似度算出手段37と、路面状態判別手段33とを備える。
特徴量類似度算出手段37は、特徴量算出手段13で算出された特徴量Xij(i=1~n、j=1~m)から成る車両Wi毎の特徴量ベクトルPi=(Xi1,Xi2,……,Xim)の分布状態と、予め路面状態R毎に求めておいた特徴量の分布状態(以下、特徴量路面モデルCRという)との類似度SRを算出する。
特徴量路面モデルCRは、標準タイヤを装着した試験車両をR路面にて一定速度で多数回(N回)走行させて得られた4(路面状態数)×N個の特徴量ベクトルPiの分布状態である。
特徴量ベクトルPiの分布状態と、特徴量路面モデルCRとの類似度は、実施の形態3と同様に、前記式[数1]に記載したような、混合ガウスモデルGR(x)などの統計モデルを用いて算出することができる。
類似度SRは、上記の式[数1]の変数xに各車両Wiの特徴量ベクトルPiを代入した値である。
すなわち、DRY路面の類似度SDは、SD=ΣiGD(Pi)により算出され、WET路面の類似度SWは、SW=ΣiGW(Pi)により算出される。また、SNOW路面の類似度SSは、SS=ΣiGS(Pi)により算出され、ICE路面の類似度SIは、SI=ΣiGI(Pi)により算出される(i=1~n)。
類似度SRも、特徴空間において、特徴量路面モデルCRと各車両Wiの特徴量Xijとの距離が近いほど大きい。
路面状態判別手段33は、最大値を取る類似度SRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
このように、特徴量ベクトルPiが、混合ガウスモデルGR(x)を構成する特徴量路面モデルMXのなかのどの路面モデルと類似しているかを調べるようにしても、所定の時間内かつ所定の範囲内の場所における路面状態Rを更に精度よく判別することができる。
前記実施の形態3では、各車両Wiの識別関数値Yiの分布状態と識別路面モデルMRとの分布状態とから路面状態Rを判定したが、路面状態判別装置30Bに代えて、図7に示す構成の路面状態判別装置30Cを設けて、車両Wi毎の特徴量Xijを構成する特徴量ベクトルPiの分布状態と、路面ごとの特徴量ベクトルの分布を表す特徴量路面モデルCRとの分布状態とから路面状態Rを判定するようにしてもよい。
路面状態判別装置30Cは、特徴量類似度算出手段37と、路面状態判別手段33とを備える。
特徴量類似度算出手段37は、特徴量算出手段13で算出された特徴量Xij(i=1~n、j=1~m)から成る車両Wi毎の特徴量ベクトルPi=(Xi1,Xi2,……,Xim)の分布状態と、予め路面状態R毎に求めておいた特徴量の分布状態(以下、特徴量路面モデルCRという)との類似度SRを算出する。
特徴量路面モデルCRは、標準タイヤを装着した試験車両をR路面にて一定速度で多数回(N回)走行させて得られた4(路面状態数)×N個の特徴量ベクトルPiの分布状態である。
特徴量ベクトルPiの分布状態と、特徴量路面モデルCRとの類似度は、実施の形態3と同様に、前記式[数1]に記載したような、混合ガウスモデルGR(x)などの統計モデルを用いて算出することができる。
類似度SRは、上記の式[数1]の変数xに各車両Wiの特徴量ベクトルPiを代入した値である。
すなわち、DRY路面の類似度SDは、SD=ΣiGD(Pi)により算出され、WET路面の類似度SWは、SW=ΣiGW(Pi)により算出される。また、SNOW路面の類似度SSは、SS=ΣiGS(Pi)により算出され、ICE路面の類似度SIは、SI=ΣiGI(Pi)により算出される(i=1~n)。
類似度SRも、特徴空間において、特徴量路面モデルCRと各車両Wiの特徴量Xijとの距離が近いほど大きい。
路面状態判別手段33は、最大値を取る類似度SRに対応する路面状態Rを、所定の時間内かつ所定の範囲内の場所における路面状態と判別する。
このように、特徴量ベクトルPiが、混合ガウスモデルGR(x)を構成する特徴量路面モデルMXのなかのどの路面モデルと類似しているかを調べるようにしても、所定の時間内かつ所定の範囲内の場所における路面状態Rを更に精度よく判別することができる。
以上、本発明を実施の形態を用いて説明したが、本発明の技術的範囲は前記実施の形態に記載の範囲には限定されない。前記実施の形態に、多様な変更または改良を加えることが可能であることが当業者にも明らかである。そのような変更または改良を加えた形態も本発明の技術的範囲に含まれ得ることが、特許請求の範囲から明らかである。
例えば、前記実施の形態1~4では、路面状態をDRY路面,WET路面,SNOW路面,ICE路面の4種類としたが、これに限るものではなく、路面摩擦係数μによりに分けたり、「高μ路(μ≧0.7)」、「中μ路(0.3<μ<0.7)」、「低μ路(μ≦0.3)」などに分けてもよい。
また、前記実施の形態1~4では、各車両Wiで採取した車両情報から路面状態を判別したが、各車両Wiに、図外のインターネットに接続して車外情報を取得するWeb情報取得手段を設けて、天候,気温,降水量,風速,日照時間などの気象情報や地形情報,交通情報,道路管理情報などの道路情報などを取得するとともに、取得された車両外情報を用いて、路面識別関数FR(X)、または、識別路面モデルMR、または、特徴量路面モデルCRを修正するか、もしくは、車両外情報毎に、路面識別関数FR(X)、識別路面モデルMR、または、特徴量路面モデルCRを設定もしくは作成しておけば、路面状態の判別精度を更に向上させることができる。
なお、Web情報取得手段をサーバー20に設け、サーバー20にて気象情報や道路情報を取得するようにしてもよい。
また、前記実施の形態1~4では、各車両Wiで採取した車両情報から路面状態を判別したが、各車両Wiに、図外のインターネットに接続して車外情報を取得するWeb情報取得手段を設けて、天候,気温,降水量,風速,日照時間などの気象情報や地形情報,交通情報,道路管理情報などの道路情報などを取得するとともに、取得された車両外情報を用いて、路面識別関数FR(X)、または、識別路面モデルMR、または、特徴量路面モデルCRを修正するか、もしくは、車両外情報毎に、路面識別関数FR(X)、識別路面モデルMR、または、特徴量路面モデルCRを設定もしくは作成しておけば、路面状態の判別精度を更に向上させることができる。
なお、Web情報取得手段をサーバー20に設け、サーバー20にて気象情報や道路情報を取得するようにしてもよい。
本発明の一態様にかかる路面状態判別システムは、車両の走行時の挙動の情報である車両情報を取得する車載センサ、前記車載センサにより取得した車両情報から前記車両の走行する路面の状態(以下、路面状態という)を推定するための複数の特徴量を算出する特徴量算出手段、及び、前記特徴量を当該車両の外部に送信する送信手段を備えた複数台の車両と、前記複数台の車両Wiの送信手段からそれぞれ送信される特徴量を蓄積するサーバーと、前記蓄積された複数台の車両の複数の特徴量Xij(i=1~n:nは車両の数、j=1~m:mは特徴量の数)を用いて路面状態の判別を行う路面状態判別装置とを備えた路面状態判別システムであって、前記送信される特徴量が、所定の時間内かつ所定の範囲内の場所における前記特徴量であり、前記路面状態判別装置は、前記複数台の車両の特徴量を用いて、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別することを特徴とする。
このように、複数台の車両の車両情報から得られた複数の特徴量を用いて路面状態の判別を行うようにしたので、路面状態の判別精度を向上させることができる。
ここで、「特徴量」とは、タイヤ内、ホイール、もしくは、ナックルに設けられた加速度センサの時系列波形の複数の特定周波数帯域における振動の大きさ(帯域値)や、タイヤ気室内に設けられた圧力センサで検出した空気圧変動信号の帯域値、あるいは、タイヤ内や後輪前方フレームに設けられたマイクロフォンで検出した音圧信号の帯域値などの、路面状態の推定に用いられる物理量を指す。
また、「所定の時間」は、車両の走行する区間や走行している時刻に応じて、長く設定したり、短く設定したりするなど、変更可能である。
また、「所定の時間」を、所定の範囲内の場所を、予め設定された台数(例えば、5台)の車両が通過する時間としてもよい。この場合も、予め設定する車両の台数は、走行する区間や走行している時刻に応じて、変更可能とする。
このように、複数台の車両の車両情報から得られた複数の特徴量を用いて路面状態の判別を行うようにしたので、路面状態の判別精度を向上させることができる。
ここで、「特徴量」とは、タイヤ内、ホイール、もしくは、ナックルに設けられた加速度センサの時系列波形の複数の特定周波数帯域における振動の大きさ(帯域値)や、タイヤ気室内に設けられた圧力センサで検出した空気圧変動信号の帯域値、あるいは、タイヤ内や後輪前方フレームに設けられたマイクロフォンで検出した音圧信号の帯域値などの、路面状態の推定に用いられる物理量を指す。
また、「所定の時間」は、車両の走行する区間や走行している時刻に応じて、長く設定したり、短く設定したりするなど、変更可能である。
また、「所定の時間」を、所定の範囲内の場所を、予め設定された台数(例えば、5台)の車両が通過する時間としてもよい。この場合も、予め設定する車両の台数は、走行する区間や走行している時刻に応じて、変更可能とする。
また、前記車載センサを、タイヤ内に配置された加速度センサとしたので、路面状態によるタイヤ振動の違いを的確に把握でき、路面状態の判別精度を更に向上させることができる。
また、本発明の一態様にかかる路面状態判別システムは、前記路面状態判別装置が、前記複数の特徴量のそれぞれについて、所定の時間内かつ所定の範囲内の場所における前記複数台の車両の特徴量Xij(i=1~n:n車両の数、jは特定の特徴量の番号)を用いた演算値である特徴量演算値Xjを算出する特徴量演算値算出手段と、前記複数の特徴量演算値Xj(j=1~m;mは特徴量の数)を、予め路面状態R毎に設定した路面識別関数YR=FR(X)にそれぞれ入力して得られる特徴量識別値YRを前記路面識別関数毎に算出する特徴量識別値算出手段と、前記算出された路面識別関数毎の特徴量識別値YRを用いて、前記所定の時間内かつ所定の範囲内の場所における路面状態Rを判別する路面状態判別手段とを備えることを特徴とする。
このように、複数台の車両の特徴量を用いた特徴量演算値Xjを路面状態毎に設定した識別関数YR=FR(X)に代入して、得られた値を路面状態Rの判別値とすることで、複数台の車両の特徴量Xijのデータ全てを用いて路面状態の判別を行うようにしたので、路面状態の判別精度を向上させることができる。
このように、複数台の車両の特徴量を用いた特徴量演算値Xjを路面状態毎に設定した識別関数YR=FR(X)に代入して、得られた値を路面状態Rの判別値とすることで、複数台の車両の特徴量Xijのデータ全てを用いて路面状態の判別を行うようにしたので、路面状態の判別精度を向上させることができる。
また、本発明の一態様にかかる路面状態判別システムは、前記路面状態判別装置が、前記複数台の車両毎に、前記複数の特徴量Xij(iは車両の番号、j=1~m:mは特徴量の数)を、予め路面状態R毎に設定した路面識別関数YR=FR(X)にそれぞれ入力して、車両毎の識別値である識別関数値YiRを、前記路面識別関数毎に算出する識別関数値算出手段と、前記算出された路面識別関数毎の識別関数値YiRを用いた演算値である識別関数演算値ZRを算出する識別関数演算値算出手段と、前記算出された識別関数演算値ZRを用いて、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別する路面状態判別手段とを備えることを特徴とする。
このように、車両毎に識別関数値YiRを求め、これらの識別関数値YiRの演算値である識別関数演算値ZRを路面状態Rの判別値としても、路面状態の判別精度を向上させることができる。
このように、車両毎に識別関数値YiRを求め、これらの識別関数値YiRの演算値である識別関数演算値ZRを路面状態Rの判別値としても、路面状態の判別精度を向上させることができる。
また、本発明の一態様にかかる路面状態判別システムは、前記路面状態判別装置が、前記複数台の車両毎に、前記複数の特徴量Xijを、予め路面状態毎に設定した路面識別関数YR=FR(X)にそれぞれ入力して、車両毎の識別値である識別関数値YiRを、前記路面識別関数毎に算出する識別関数値算出手段と、前記算出された路面識別関数毎の識別関数値YiRの分布状態と、予め求めておいた路面状態R毎の識別関数値の分布状態(以下、識別路面モデルMRという)との類似度である識別値類似度HRを計算し、前記計算された識別値類似度HRから、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別する路面状態判別手段とを備えることを特徴とする。
このように、予め、路面状態毎に識別路面モデルMRを作成しておき、実際に計測した特徴量Xijから求められる識別関数値YiRの分布状態が、例えば、車両をDRY路面上で走行させたときに得られた特徴量から求めた識別路面モデルMDに近いか、あるいは、他の路面R’上を走行させたときに得られた特徴量から求めた識別路面モデルMR’に近いかを路面判別の判別基準としたので、路面状態の判別精度を一層向上させることができる。
このように、予め、路面状態毎に識別路面モデルMRを作成しておき、実際に計測した特徴量Xijから求められる識別関数値YiRの分布状態が、例えば、車両をDRY路面上で走行させたときに得られた特徴量から求めた識別路面モデルMDに近いか、あるいは、他の路面R’上を走行させたときに得られた特徴量から求めた識別路面モデルMR’に近いかを路面判別の判別基準としたので、路面状態の判別精度を一層向上させることができる。
また、本発明の一態様にかかる路面状態判別システムは、前記路面状態判別装置が、前記複数台の車両のそれぞれについて、前記特徴量Xijの分布状態を、予め求めておいた路面状態毎の特徴量の分布状態(以下、特徴量路面モデルCRという)との類似度である特徴量類似度SR計算する特徴量類似度算出手段と、前記計算された各車両の特徴量類似度SRから、所定の時間内かつ所定の範囲内の場所における路面状態を判別する路面状態判別手段とを備えることを特徴とする。
このように、予め、路面状態毎に特徴量路面モデルCRを作成しておき、実際に計測した特徴量Xijの分布状態が、例えば、車両をDRY路面上を走行させたときに得られた特徴量から求めた特徴量路面モデルCDに近いか、あるいは、他の路面R’上を走行させたときに得られた特徴量から求めた特徴量路面モデルCR’に近いかを路面判別の判別基準としても、路面状態の判別精度を一層向上させることができる。
このように、予め、路面状態毎に特徴量路面モデルCRを作成しておき、実際に計測した特徴量Xijの分布状態が、例えば、車両をDRY路面上を走行させたときに得られた特徴量から求めた特徴量路面モデルCDに近いか、あるいは、他の路面R’上を走行させたときに得られた特徴量から求めた特徴量路面モデルCR’に近いかを路面判別の判別基準としても、路面状態の判別精度を一層向上させることができる。
また、本発明の一態様にかかる路面状態判別システムは、車両外情報を取得する手段を設けるとともに、前記取得された車両外情報を用いて、前記路面識別関数FR(X)、または、路面状態毎の識別関数値の分布状態(識別路面モデルMR)、または、路面状態毎の特徴量の分布状態(特徴量路面モデルCR)を修正することを特徴とする。
ここで、車両外情報とは、インターネットなどにより取得した、天候,気温,降水量,風速,日照時間などの気象情報、もしくは、地形情報,交通情報,道路管理情報などの道路情報などを指す。
また、修正に用いられる情報としては、他の所定の範囲内の場所の路面状態の推定結果などがある。具体的には、過去の条件(天候、交通量、路面状態の推定結果)において、所定の範囲内の場所Aと相関がある他の所定の範囲内の場所Bの路面状態の推定結果を用いて、所定の範囲内の場所Aでの推定結果を補正することができる。
これにより、所定の時間内かつ所定の範囲内の場所における路面状態を更に精度よく判別することができる。
ここで、車両外情報とは、インターネットなどにより取得した、天候,気温,降水量,風速,日照時間などの気象情報、もしくは、地形情報,交通情報,道路管理情報などの道路情報などを指す。
また、修正に用いられる情報としては、他の所定の範囲内の場所の路面状態の推定結果などがある。具体的には、過去の条件(天候、交通量、路面状態の推定結果)において、所定の範囲内の場所Aと相関がある他の所定の範囲内の場所Bの路面状態の推定結果を用いて、所定の範囲内の場所Aでの推定結果を補正することができる。
これにより、所定の時間内かつ所定の範囲内の場所における路面状態を更に精度よく判別することができる。
なお、上記発明の概要は、本発明の必要な全ての特徴を列挙したものではなく、これらの特徴群のサブコンビネーションもまた、発明となり得る。
以上説明したように、本発明によれば、前記所定の時間内かつ所定の範囲内の場所における路面状態を精度よく判別することができるので、この判別結果を、所定の範囲内の場所を走行する車両に通知するなどすれば、車両の走行安全性を向上させることができる。
1;3 路面状態判別システム、2 路面状態管理センター、
11 加速度センサ、12 車輪速センサ、13 特徴量算出手段、
14 GPS装置、15 車両情報収集手段、16 送信機、
17 受信機、20 サーバー、21 受信機、22 データ記憶手段、
23 送信機、30,30A~30C 路面状態判別装置、
31 特徴量演算値算出手段、32 特徴量識別値算出手段、
33 路面状態判別手段、34 識別関数値算出手段、
35 識別関数演算値算出手段、36 識別値類似度算出手段、
37 特徴量類似度算出手段、40 タイヤ、41 インナーライナー部、
42 タイヤ気室、43 トレッド、44 リム、W1~Wn 車両。
11 加速度センサ、12 車輪速センサ、13 特徴量算出手段、
14 GPS装置、15 車両情報収集手段、16 送信機、
17 受信機、20 サーバー、21 受信機、22 データ記憶手段、
23 送信機、30,30A~30C 路面状態判別装置、
31 特徴量演算値算出手段、32 特徴量識別値算出手段、
33 路面状態判別手段、34 識別関数値算出手段、
35 識別関数演算値算出手段、36 識別値類似度算出手段、
37 特徴量類似度算出手段、40 タイヤ、41 インナーライナー部、
42 タイヤ気室、43 トレッド、44 リム、W1~Wn 車両。
Claims (7)
- 車両の走行時の挙動の情報である車両情報を取得する車載センサ、前記車載センサにより取得した車両情報から前記車両の走行する路面の状態を推定するための複数の特徴量を算出する特徴量算出手段、及び、前記特徴量を当該車両の外部に送信する送信手段を備えた複数台の車両と、
前記複数台の車両の送信手段からそれぞれ送信される特徴量を蓄積するサーバーと、
前記蓄積された複数台の車両の特徴量を用いて路面状態の判別を行う路面状態判別装置とを備えた路面状態判別システムであって、
前記送信される特徴量が、所定の時間内かつ所定の範囲内の場所における前記特徴量であり、
前記路面状態判別装置は、
前記複数台の車両の特徴量を用いて、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別することを特徴とする路面状態判別システム。 - 前記車載センサが、タイヤ内に配置された加速度センサであることを特徴とする請求項1に記載の路面状態判別システム。
- 前記路面状態判別装置は、
前記複数の特徴量のそれぞれについて、所定の時間内かつ所定の範囲内の場所における前記複数台の車両の特徴量を用いた演算値である特徴量演算値を算出する特徴量演算値算出手段と、
前記複数の特徴量演算値を、予め路面状態毎に設定した路面識別関数にそれぞれ入力して得られる特徴量識別値を、前記路面識別関数毎に算出する特徴量識別値算出手段と、
前記算出された路面識別関数毎の特徴量識別値を用いて、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別する路面状態判別手段とを備えることを特徴とする請求項1または請求項2に記載の路面状態判別システム。 - 前記路面状態判別装置は、
前記複数台の車両毎に、前記複数の特徴量を、予め路面状態毎に設定した路面識別関数にそれぞれ入力して、車両毎の識別値である識別関数値を、前記路面識別関数毎に算出する識別関数値算出手段と、
前記算出された路面識別関数毎の識別関数値を用いた演算値である識別関数演算値を算出する識別関数演算値算出手段と、
前記算出された識別関数演算値を用いて、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別する路面状態判別手段とを備えることを特徴とする請求項1または請求項2に記載の路面状態判別システム。 - 前記路面状態判別装置は、
前記複数台の車両毎に、前記複数の特徴量を、予め路面状態毎に設定した路面識別関数にそれぞれ入力して、車両毎の識別値である識別関数値を、前記路面識別関数毎に算出する識別関数値算出手段と、
前記算出された路面識別関数毎の識別関数値の分布状態と、予め求めておいた路面状態毎の識別関数値の分布状態との類似度である識別値類似度を計算し、前記計算された識別値類似度から、前記所定の時間内かつ所定の範囲内の場所における路面状態を判別する路面状態判別手段とを備えることを特徴とする請求項1または請求項2に記載の路面状態判別システム。 - 前記路面状態判別装置は、
前記複数台の車両のそれぞれについて、前記特徴量の分布状態を、予め求めておいた路面状態毎の特徴量の分布状態との類似度である特徴量類似度を計算する特徴量類似度算出手段と、
前記計算された各車両の特徴量類似度から、所定の時間内かつ所定の範囲内の場所における路面状態を判別する路面状態判別手段とを備えることを特徴とする請求項1または請求項2に記載の路面状態判別システム。 - 車両外情報を取得する手段を設けるとともに、
前記取得された車両外情報を用いて、前記路面識別関数、または、路面状態毎の識別関数値の分布状態、または、路面状態毎の特徴量の分布状態を修正することを特徴とする請求項3~請求項6のいずれかに記載の路面状態判別システム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201580066169.XA CN107003429B (zh) | 2014-10-06 | 2015-10-06 | 路面状态判断系统 |
US15/516,970 US10377385B2 (en) | 2014-10-06 | 2015-10-06 | Road surface condition determining system |
EP15848631.6A EP3206056B1 (en) | 2014-10-06 | 2015-10-06 | Road surface conditions distinguishing system |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014-206017 | 2014-10-06 | ||
JP2014206017A JP6408852B2 (ja) | 2014-10-06 | 2014-10-06 | 路面状態判別システム |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2016056569A1 true WO2016056569A1 (ja) | 2016-04-14 |
Family
ID=55653180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2015/078398 WO2016056569A1 (ja) | 2014-10-06 | 2015-10-06 | 路面状態判別システム |
Country Status (5)
Country | Link |
---|---|
US (1) | US10377385B2 (ja) |
EP (1) | EP3206056B1 (ja) |
JP (1) | JP6408852B2 (ja) |
CN (1) | CN107003429B (ja) |
WO (1) | WO2016056569A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018054910A3 (en) * | 2016-09-22 | 2018-05-03 | Omniklima Ab | Method and arrangement for determining a condition of a road surface |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6544302B2 (ja) * | 2016-06-22 | 2019-07-17 | 株式会社Soken | 路面状況推定装置 |
JP6624152B2 (ja) * | 2017-04-26 | 2019-12-25 | 株式会社Soken | タイヤ側装置およびそれを含むタイヤ装置 |
JP6783184B2 (ja) * | 2017-05-12 | 2020-11-11 | 株式会社ブリヂストン | 路面状態判別方法及び路面状態判別装置 |
JP6930355B2 (ja) * | 2017-10-11 | 2021-09-01 | 株式会社Soken | 路面状態判別装置およびそれを備えたタイヤシステム |
DE102017219519B4 (de) * | 2017-11-02 | 2019-06-27 | Ford Global Technologies, Llc | Luftreifen für ein Fahrzeug und Verfahren zu seiner Herstellung |
NO20171853A1 (en) * | 2017-11-21 | 2019-05-22 | Fm Equipment As | Assembly and method for a vehicle |
WO2019113022A1 (en) * | 2017-12-04 | 2019-06-13 | University Of Massachusetts | Method to measure road roughness characteristics and pavement induced vehicle fuel consumption |
WO2019131568A1 (ja) * | 2017-12-27 | 2019-07-04 | 株式会社デンソー | 路面状態判別装置 |
JP6969399B2 (ja) * | 2018-01-19 | 2021-11-24 | 株式会社Soken | タイヤシステム |
JP2019185647A (ja) * | 2018-04-17 | 2019-10-24 | スズキ株式会社 | 路面環境情報収集システム |
CN110263844B (zh) * | 2019-06-18 | 2021-04-06 | 北京中科原动力科技有限公司 | 一种在线学习和实时估计路面状态的方法 |
JP7306943B2 (ja) * | 2019-09-30 | 2023-07-11 | Toyo Tire株式会社 | 車両走行条件評価方法及びシステム |
JP7211350B2 (ja) * | 2019-11-29 | 2023-01-24 | トヨタ自動車株式会社 | 路面損傷検出装置、路面損傷検出方法、プログラム |
JP7211349B2 (ja) * | 2019-11-29 | 2023-01-24 | トヨタ自動車株式会社 | 路面損傷検出装置、路面損傷検出方法、プログラム |
US11436877B2 (en) * | 2019-12-26 | 2022-09-06 | Globiz Co., Ltd. | Apparatus for determining overworked consumable for vehicle |
JP7347219B2 (ja) * | 2020-01-06 | 2023-09-20 | トヨタ自動車株式会社 | 道路診断システム |
EP4425463A2 (en) * | 2020-02-25 | 2024-09-04 | Huawei Technologies Co., Ltd. | Method and apparatus for recognizing special road condition, electronic device, and storage medium |
CN111452800B (zh) * | 2020-03-06 | 2021-05-04 | 武汉理工大学 | 水陆两栖车的参数输出方法、装置和存储介质 |
CN111806166B (zh) * | 2020-06-24 | 2021-08-31 | 浙江大学 | 一种基于摩擦起电材料的地形监测轮胎 |
US11878695B2 (en) * | 2021-01-26 | 2024-01-23 | Motional Ad Llc | Surface guided vehicle behavior |
EP4140782A1 (en) * | 2021-08-27 | 2023-03-01 | Bridgestone Europe NV/SA | A method to synchronize the time of a tire-mounted sensor to the road impact and measure the contact patch duration and amplitude |
KR102414233B1 (ko) * | 2022-01-18 | 2022-06-29 | (주)에스알디코리아 | 주행차량의 노면인지에 의한 제어 시스템 및 그 방법 |
CN115272930A (zh) * | 2022-07-28 | 2022-11-01 | 广西北投交通养护科技集团有限公司 | 一种基于探地雷达的路面状态评估方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002008198A (ja) * | 2000-04-21 | 2002-01-11 | Sumitomo Rubber Ind Ltd | 路面情報配信システム、車両の情報集約配信システムおよび車両の情報発信装置、ならびに車両制御プログラム |
US6919821B1 (en) * | 2000-05-19 | 2005-07-19 | Navteq North America, Llc | Method and system for collecting meteorological data using in-vehicle systems |
JP2008297872A (ja) * | 2007-06-04 | 2008-12-11 | Japan Radio Co Ltd | 路面性状検知装置における路面性状検知方法 |
US20090312943A1 (en) * | 2008-02-04 | 2009-12-17 | Kelly Terence F | System and Method for Weather Mapping to Road Segments |
JP2011242303A (ja) * | 2010-05-19 | 2011-12-01 | Bridgestone Corp | 路面状態推定方法 |
US20140081573A1 (en) * | 2012-09-20 | 2014-03-20 | Google Inc. | Detecting road weather conditions |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3473158B2 (ja) | 1995-03-22 | 2003-12-02 | 住友電気工業株式会社 | 路面状態検出装置 |
US6577943B2 (en) * | 2000-04-21 | 2003-06-10 | Sumitomo Rubber Industries, Ltd. | System for distributing road surface information, system for collecting and distributing vehicle information, device for transmitting vehicle information and program for controlling vehicle |
US6853894B1 (en) * | 2000-04-24 | 2005-02-08 | Usa Technologies, Inc. | Global network based vehicle safety and security telematics |
JP4554176B2 (ja) | 2003-08-19 | 2010-09-29 | 株式会社ブリヂストン | 路面状態推定方法 |
US7355509B2 (en) | 2005-02-25 | 2008-04-08 | Iwapi Inc. | Smart modem device for vehicular and roadside applications |
WO2006135090A1 (ja) * | 2005-06-17 | 2006-12-21 | Kabushiki Kaisha Bridgestone | 路面状態推定方法、路面状態推定用タイヤ、路面状態推定装置、及び、車両制御装置 |
JP4427759B2 (ja) * | 2007-06-29 | 2010-03-10 | アイシン・エィ・ダブリュ株式会社 | 車両挙動学習装置及び車両挙動学習プログラム |
US8451140B2 (en) * | 2007-12-20 | 2013-05-28 | International Business Machines Corporation | Monitoring road surface conditions |
JP5837341B2 (ja) * | 2011-06-24 | 2015-12-24 | 株式会社ブリヂストン | 路面状態判定方法とその装置 |
US8942861B2 (en) * | 2011-07-20 | 2015-01-27 | Bridgestone Corporation | Road surface condition estimation method, and road surface condition estimation apparatus |
JP5937921B2 (ja) * | 2012-08-09 | 2016-06-22 | 株式会社ブリヂストン | 路面状態判別方法とその装置 |
DE102012219631A1 (de) * | 2012-10-26 | 2014-04-30 | Robert Bosch Gmbh | Verfahren und Vorrichtung zum Detektieren von zumindest einer Fahrbahnunebenheit |
CN104755345B (zh) * | 2014-03-19 | 2019-04-19 | 株式会社小松制作所 | 路面状态判定方法及装置、路面状态输出方法及装置 |
-
2014
- 2014-10-06 JP JP2014206017A patent/JP6408852B2/ja active Active
-
2015
- 2015-10-06 EP EP15848631.6A patent/EP3206056B1/en active Active
- 2015-10-06 CN CN201580066169.XA patent/CN107003429B/zh active Active
- 2015-10-06 US US15/516,970 patent/US10377385B2/en active Active
- 2015-10-06 WO PCT/JP2015/078398 patent/WO2016056569A1/ja active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002008198A (ja) * | 2000-04-21 | 2002-01-11 | Sumitomo Rubber Ind Ltd | 路面情報配信システム、車両の情報集約配信システムおよび車両の情報発信装置、ならびに車両制御プログラム |
US6919821B1 (en) * | 2000-05-19 | 2005-07-19 | Navteq North America, Llc | Method and system for collecting meteorological data using in-vehicle systems |
JP2008297872A (ja) * | 2007-06-04 | 2008-12-11 | Japan Radio Co Ltd | 路面性状検知装置における路面性状検知方法 |
US20090312943A1 (en) * | 2008-02-04 | 2009-12-17 | Kelly Terence F | System and Method for Weather Mapping to Road Segments |
JP2011242303A (ja) * | 2010-05-19 | 2011-12-01 | Bridgestone Corp | 路面状態推定方法 |
US20140081573A1 (en) * | 2012-09-20 | 2014-03-20 | Google Inc. | Detecting road weather conditions |
Non-Patent Citations (1)
Title |
---|
See also references of EP3206056A4 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018054910A3 (en) * | 2016-09-22 | 2018-05-03 | Omniklima Ab | Method and arrangement for determining a condition of a road surface |
Also Published As
Publication number | Publication date |
---|---|
US10377385B2 (en) | 2019-08-13 |
CN107003429B (zh) | 2020-02-28 |
JP6408852B2 (ja) | 2018-10-17 |
US20170297580A1 (en) | 2017-10-19 |
CN107003429A (zh) | 2017-08-01 |
EP3206056A1 (en) | 2017-08-16 |
EP3206056B1 (en) | 2019-01-23 |
JP2016075575A (ja) | 2016-05-12 |
EP3206056A4 (en) | 2017-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6408852B2 (ja) | 路面状態判別システム | |
US11993260B2 (en) | System and method for predicting tire traction capabilities and active safety applications | |
JP6408854B2 (ja) | 路面状態予測方法 | |
US8942861B2 (en) | Road surface condition estimation method, and road surface condition estimation apparatus | |
JP5657917B2 (ja) | 路面状態推定方法 | |
JP6392623B2 (ja) | 路面状態判別システム | |
JP5878612B2 (ja) | 路面状態推定方法 | |
CN105793687B (zh) | 轮胎异常磨损的估计方法和估计设备 | |
US11505015B2 (en) | Determining a tire pressure status in a vehicle | |
EP3800100A1 (en) | Road surface condition estimation apparatus and road surface condition estimation method using the same | |
CN110997441B (zh) | 用于求取驾驶行为的方法和设备 | |
WO2020031513A1 (ja) | タイヤ種判別方法及びタイヤ種判別装置 | |
JP5908953B2 (ja) | 車両制御方法及び路面状態警告方法 | |
JP7551934B2 (ja) | 音響フットプリント分析を使用してタイヤ摩耗を推定するためのシステム及び方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 15848631 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 15516970 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REEP | Request for entry into the european phase |
Ref document number: 2015848631 Country of ref document: EP |