CN105575120A - Floating car data parking behavior mode cleaning method specific to road real time speed calculation - Google Patents

Floating car data parking behavior mode cleaning method specific to road real time speed calculation Download PDF

Info

Publication number
CN105575120A
CN105575120A CN201511031346.9A CN201511031346A CN105575120A CN 105575120 A CN105575120 A CN 105575120A CN 201511031346 A CN201511031346 A CN 201511031346A CN 105575120 A CN105575120 A CN 105575120A
Authority
CN
China
Prior art keywords
speed
ldid
gps
behavior
fdc
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.)
Granted
Application number
CN201511031346.9A
Other languages
Chinese (zh)
Other versions
CN105575120B (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.)
Yinjiang Technology Co.,Ltd.
Original Assignee
Enjoyor 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 Enjoyor Co Ltd filed Critical Enjoyor Co Ltd
Priority to CN201511031346.9A priority Critical patent/CN105575120B/en
Publication of CN105575120A publication Critical patent/CN105575120A/en
Application granted granted Critical
Publication of CN105575120B publication Critical patent/CN105575120B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0108Measuring and analyzing of parameters relative to traffic conditions based on the source of data
    • G08G1/0112Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/052Detecting movement of traffic to be counted or controlled with provision for determining speed or overspeed

Abstract

Provided is a floating car data parking behavior mode cleaning method specific to road real time speed calculation, comprising the steps of: (1) reading map data which comprises road section numbers and road section point location information, and gridding a map to obtain a ''grid-road section number'' mapping table; (2) reading in an ascending sequence the floating car data of a preset time slice before the current time according to ''license plate number'' and ''creation time''; (3) calculating the vehicle distance and time difference between twice continuous recordings, calculating the speed according to the distance and time difference, and adding the speed to the road section number where the vehicle point location information of secondary recording is located according to the ''grid-road section number'' mapping table; (4) obtaining the speed list of each vehicle at each road section of the time slice; and (5) identifying and rejecting abnormal stop points in the speed list caused by ''pseudo stop behaviors'', and returning to the speed list. The floating car data parking behavior mode cleaning method can effectively identify floating car parking behaviors.

Description

Towards the floating car data parking behavior pattern cleaning method that road real-time speed calculates
Technical field
The present invention relates in data mining technical field, be specifically related to a kind of floating car data parking behavior pattern cleaning method calculated towards road real-time speed.
Background technology
Road speeds is point duty, navigation, induced travel, block up administer etc. macroscopic view, microcosmic application prerequisite, be important basic data.Therefore obtain more close to real data
Floating car data, as one of the important component part of traffic data, has calculated at Real-time Road passage rate and has applied widely.Can find during existing floating car data research experiment, the data that the Floating Car behavior that abends produces are applied to the accuracy that can have a strong impact on its speed during Real-time Road passage rate calculates.Because floating car data amount is large, and the behavior individual behaviour and irregularly to follow of abending, artificial being difficult to finds timely and rejects.
Summary of the invention
In order to overcome the existing deficiency that effectively cannot identify parking behavior, the invention provides a kind of floating car data parking behavior pattern cleaning method calculated towards road real-time speed of effective identification Floating Car parking behavior.
The technical solution adopted for the present invention to solve the technical problems is:
Towards the floating car data parking behavior pattern cleaning method that road real-time speed calculates, described cleaning method comprises the steps:
Step (1). read map datum, map datum comprises section numbering (LDID), some position, section information, map grid is obtained " grid-section numbering " mapping table;
Step (2). according to the floating car data of setting-up time sheet before " license plate number ", " creation-time " ascending order reading current time, floating car data comprises number-plate number CPHM, some position information GPS_FDC and creation-time CJSJ;
Step (3). calculate the Distance geometry mistiming of the double record of vehicle, then by Distance geometry mistiming computing velocity, and according to " grid-section numbering " mapping table, this speed is added in the section numbering at the information place, vehicle point position of second time record;
Step (4). the speed list of each vehicle on each section of this timeslice is obtained according to above step;
Step (5). identify and reject abending a little of causing due to " pseudo-stopping behavior " in speed list, and return speed list, process is as follows:
In the data structure LDID_SpeedMap that 5-1. obtains in step (4), searching loop obtains CPHM in LDID kspeed list speedList j, searching loop speedList jin all speed v iif, v i=0 and v i+s≠ 0,0<s<=n-i, s are the number of continuous halt, then the car speed record number δ=s+1 of " pseudo-stopping behavior " process;
5-2. calculates " pseudo-stopping behavior " car speed sudden change threshold value SPSSC according to formula (2); When i ≠ 1, the formula left side is the average velocity of all speed records before this speed list first halt, as i=1, with the average velocity v in a moment on this section tas average velocity threshold value;
S P S S C = &Sigma; f = 1 i - 1 v f i - 1 * &delta; , i &NotEqual; 1 v t * &delta; , i = 1 - - - ( 2 )
If 5-3. is v i+s>=SPSSC, then this speed record is the mutating speed that " pseudo-stopping behavior " causes, and rejects speed record and the mutating speed speed record of continuous halt from speed list;
5-4. completes " pseudo-stopping behavior " identifies and returns a new LDID_SpeedMap after rejecting.
Further, described cleaning method also comprises the steps:
Step (6). identify and reject abending a little of causing due to " really stopping behavior " in speed list, and return speed list, process is as follows:
6-1. after step (5) from new LDID_SpeedMap searching loop obtain CPHM in LDID kspeed list speedList j, searching loop speedList jin all speed v iif, v i=0 and v i+s≠ 0,0<s<=n-i, s are the number of continuous halt; If s>=3, think that this car has occurred that " stopping behavior " skips to 6-2 and judge whether to belong to " crossing stopping behavior ";
6-2. obtains the v of this stopping behavior icorresponding LDID and GPS_FDC i, traversal intersectionMap, if there is this LDID in intersectionMap, illustrates that this LDID is section, crossing, enters 6-3 if do not exist; From intersectionMap, obtain crossing point (GPS_QD, GPS_ZD) that this LDID is corresponding, calculate GPS_FDC by formula (3) iget its minimum value L with the distance at crossing, if L<L_min, then think that this " stopping behavior " belongs to " crossing stopping behavior " and do not reject, otherwise enter 6-3 and judge whether to belong to " collective stops behavior ";
L=min(abs(GPS_FDC i-GPS_QD),abs(GPS_FDC i-GPS_ZD)(3)
6-3. obtains the v of this stopping behavior icorresponding CJSJ iand GPS_FDC iif the speed list of other vehicles occurs that " stopping behavior " obtains the GPS_FDC of its corresponding halt under this LDID of searching loop i'and CJSJ i', calculate Rule of judgment according to formula (4), formula (5),
abs(GPS_FDC i-GPS_FDC i')<Δs_min,i≠i'(4)
abs(CJSJ i-CJSJ i')<Δt_min,i≠i'(5)
If other vehicles all meet above Rule of judgment under this LDID, " collective stops behavior " does not reject to think that this stopping behavior belonging to, otherwise this stopping behavior belonging to " really stopping behavior ", and rejects the speed record of continuous halt from speed list;
6-4. completes " really stopping behavior " identifies and returns a new LDID_SpeedMap after rejecting.
Further again, described cleaning method also comprises the steps:
Step (7). after above step, calculate the average velocity of all speed records corresponding to LDID and be stored in database real-time speed table, update time, sheet repeated step (2) to step (6).
Further, in described step (3), velocity computation process is as follows:
3-1. obtains all records of each vehicle under this timeslice according to step (2), is designated as { CPHM k} _ { GPS_FDC i, CJSJ i, k ∈ Ki ∈ n, n is the number of all records, and K is the sum of this timeslice vehicle;
3-2. calculates car speed v according to formula (1.1) or (1.2) i, obtain car speed list { CPHM k} _ { v i, GPS_FDC i, CJSJ i.
What adopt is the absolute distance of the GPS point of double record, as the operating range of vehicle,
1) when not across using formula (1.1) when section
v i = a b s ( G P S _ FDC i - G P S _ FDC i - 1 ) CJSJ i - CJSJ i - 1 - - - ( 1.1 )
2) when across using formula (1.2) when section
v i = a b s ( G P S _ FDC i - G P S _ JD i ) + a b s ( G P S _ JD i - G P S _ FDC i - 1 ) CJSJ i - CJSJ i - 1 - - - ( 1.2 )
GPS_JD ifor the GPS point of double record, when across section, the some position information of the intersection point in two sections.
In described step (7), obtain the speed record of all car speed lists under LDID according to LDID_SpeedMap and calculate average velocity corresponding to LDID according to formula (6), and applying for other stored in database Real-time Road velometer.
v a v g = &Sigma; k = 1 K &Sigma; i = 1 n k v i K &CenterDot; n k - - - ( 6 )
K is the vehicle number under this LDID, n kfor the speed record number of a kth car.
Technical conceive of the present invention is: according to real-time floating car data, calculates floating vehicle travelling speed and adds in corresponding section, obtains the speed list of each vehicle in section and identifies the behavior that abends and rejecting abnormalities speed.
Described Floating Car refers to taxi in urban public transport and bus.Floating Car (can be generally 30 seconds) at set intervals and upload a secondary data in vehicle travel process, and data comprise: Floating Car license plate number, Floating Car positional information (GPS), data creation time etc.; " stopping behavior " comprises " normally stopping behavior " and " abend behavior "." normally stop behavior " to comprise:
" crossing stopping behavior ": vehicle is parked in crossing (in most cases waiting for that red light belongs to the normal phenomenon in traffic process)
" collective stop behavior ": refer to the behavior (there will be the normal phenomenon also belonged in traffic process when running into and blocking up) that stops of same place at one time of all Floating Car.
" abend behavior " comprising:
" really stop behavior ": vehicle abnormality stops (individual vehicle is parked in a point for a long time, and in traffic process, belong to individual behaviour is abnormal occurrence)
" pseudo-stopping behavior ": the abnormal stopping caused of data upload, shows as and occur that speed is the point of 0 continuously, and velocity jump after this point.(as shown in Figure 1).
" really stopping behavior ", " pseudo-stopping behavior " is distinguished from the feature of data itself." really stop behavior " and comprise " true stopping _ individual behavior " and " true stopping _ group behavior ".This programme identification be " true stopping _ individual behavior ", the individual behavior of this stopping is disturbing factor when calculating road speeds, will not count sample data when calculating road real-time speed.For " pseudo-stopping behavior " that this patent is set forth, mainly refer to that those are because the looking of causing such as equipment failure, transmission network fault or database Write fault is as the data sample stopping behavior.Although these stopping behaviors showing as stopping on data mode, the possibility that data are produced and transmission link mistake causes is larger, and in other words, we think that it is " pseudo-stopping behavior ".This " pseudo-stopping behavior " is also got rid of outside the sample data calculating road real-time speed.
Beneficial effect of the present invention is mainly manifested in: two kinds of stopping behaviors at big data quantity and in irregular governed floating car data fast and accurately in identification data are also rejected, while the real-time ensureing data, significantly increase the accuracy of real-time speed.
Accompanying drawing explanation
Fig. 1 is " pseudo-stopping behavior " key diagram.
Fig. 2 is the operating range calculating schematic diagram across vehicle when section.
Fig. 3 is the process flow diagram of the floating car data parking behavior pattern cleaning method towards the calculating of road real-time speed.
Fig. 4 is LDID_SpeedMap data structure diagram.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described.
With reference to Fig. 1 ~ Fig. 4, a kind of floating car data parking behavior pattern cleaning method calculated towards road real-time speed, definition related symbol is as follows:
V i: car speed, unit km/h, i ∈ n, n are the speed record number of certain car on same LDID under single timeslice.
V avg: Road average-speed, unit km/h.
δ: the car speed record number of " pseudo-stopping behavior " process (speed be continuously 0 number of times add 1 velocity jump).
SPSSC (StopPointSpeedSuddenChange): " pseudo-stopping behavior " car speed sudden change threshold value.
Δ t_min: the speed record mistiming threshold value (establishing value 30 seconds) of same section different vehicle.
Δ s_min: the speed record location interval threshold value (arranging 50 meters) of same section different vehicle
L_min: judge whether vehicle is parked in the distance threshold (arranging 50 meters) at crossing.
Described floating car data parking behavior pattern cleaning method comprises the steps:
Step (1). map datum pre-service:
1-1. reads map datum, map datum comprises section numbering (LDID), starting point point position information (GPS_QD), terminal point position information (GPS_ZD), mid point point position information (GPS_MD), map grid (grid length of side 100-300 rice) is obtained " grid-section numbering " mapping table (as shown in table 1).
Table 1
1-2. is because the connectivity in section, and between section 1 with section 2, terminal is connected with starting point (being same point), and in statistical map data, the starting point in all sections and terminal, obtain " point-section is numbered " mapping table (as shown in table 2).
Table 2
1-3. finds the point of fork in the road or crossroad and LDID number to be more than or equal to the point of 3 according to " point-section numbering " mapping table, form new " point-section numbering " mapping table to be called " crossing point mapping table ", " crossing point mapping table " is carried out upset formation " section numbering-point " mapping table (as shown in table 3, note: what record in Value is the some position information that this section belongs to fork in the road or crossroad, if the starting point in certain section and terminal all belong to fork in the road or four crossway stomion, in Value, record starting point point position information and the terminal point position information in this section, if only have wherein some points to be fork in the road or four crossway stomion, then only record the some position information of this point) be called that " section, crossing mapping table " remembers intersectionMap.
Table 3
Step (2). read the floating car data of current time toward previous timeslice (5min) according to " number-plate number ", " creation-time " ascending order, floating car data comprises the number-plate number (CPHM), some position information (GPS_FDC), creation-time (CJSJ).
Step (3). calculate road speeds:
3-1. obtains all records of each vehicle under this timeslice according to step (2), is designated as { CPHM k} _ { GPS_FDC i, CJSJ i, k ∈ Ki ∈ n, n is the number of all records, and K is the sum of this timeslice vehicle.
3-2. calculates car speed v according to formula (1.1) or (1.2) i, obtain car speed list { CPHM k} _ { v i, GPS_FDC i, CJSJ i.
What adopt is the absolute distance of the GPS point of double record, as the operating range of vehicle.
1) when not across using formula (1.1) when section
v i = a b s ( G P S _ FDC i - G P S _ FDC i - 1 ) CJSJ i - CJSJ i - 1 - - - ( 1.1 )
2) when across using formula (1.2) when section
v i = a b s ( G P S _ FDC i - G P S _ JD i ) + a b s ( G P S _ JD i - G P S _ FDC i - 1 ) CJSJ i - CJSJ i - 1 - - - ( 1.2 )
GPS_JD ifor the GPS point of double record, when across section, the some position information of the intersection point in two sections.
Step (4). according to the GPS_FDC in speed list icoupling " grid-section numbering " mapping table, obtains GPS_FDC iall LDID of affiliated grid, calculate GPS_FDC icalculate distance with the GPS_MD of each LDID, get apart from the section of minimum LDID as this Point matching, and this speed list is added in LDID, obtain data structure LDID_SpeedMap (as shown in Figure 3).Step (5). identify and reject " pseudo-stopping behavior ":
5-1. searching loop in data structure LDID_SpeedMap obtains CPHM in LDID kspeed list speedList j, searching loop speedList jin all speed v iif, v i=0 (i.e. halt) and v i+s≠ 0 (0<s<=n-i, s are the number of continuous halt), then the car speed record number δ=s+1 of " pseudo-stopping behavior " process.
5-2. calculates " pseudo-stopping behavior " car speed sudden change threshold value SPSSC according to formula (2).When i ≠ 1, the formula left side is the average velocity of all speed records before this speed list first halt, as i=1, with the average velocity v in a moment on this section tas average velocity threshold value.
S P S S C = &Sigma; f = 1 i - 1 v f i - 1 * &delta; , i &NotEqual; 1 v t * &delta; , i = 1 - - - ( 2 )
If 5-3. is v i+s>=SPSSC, then this speed record is the mutating speed that " pseudo-stopping behavior " causes, and rejects speed record and the mutating speed speed record of continuous halt from speed list.
5-4. completes " pseudo-stopping behavior " identifies and returns a new LDID_SpeedMap after rejecting.
Citing:
In data structure LDID_SpeedMap, obtain the speed list of certain section car, be (35.52,27.89,0,0,131.07,32.12,9.1,33.17) in chronological order, unit is km/h.Now we can find that continuous halt number is 2 then " pseudo-stopping behavior " SPSSC (velocity jump threshold value)=(35.52+27.89)/2* (2+1)=31.7*3=95.1; 95.1<131.07 then thinks that the behavior is " pseudo-stopping behavior ".
Step (6). identify and reject " really stopping behavior "
Citing: " really stopping behavior ", as shown in the table: in same section, sheet only has a car to occur parking behavior at one time.Namely the behavior is " really stopping behavior "
Zhejiang CT1**0 Speed Zhejiang CT3**3 Speed Zhejiang CT3**7 Speed Zhejiang CT5**6 Speed
8:55:45 12.29 8:55:24 0 8:55:35 34.4 8:55:49 55.5
8:56:03 42.26 8:55:36 0 8:56:07 24.33 8:56:15 28.56
8:56:25 27.62 8:55:57 4.93 8:56:31 85.7 8:56:51 0
8:56:44 8.29 8:56:18 11.14 8:57:02 25.41 8:57:14 36.2
8:56:44 8.29 8:56:37 19.19 8:57:32 75.01 8:57:50 32.5
8:57:04 47.48 8:56:55 8 8:58:03 21.49 8:58:14 47.82
8:57:25 50.36 8:56:55 8 8:58:32 0 8:58:50 22.25
8:57:43 50.5 8:57:17 0 8:58:56 11.25 8:59:15 28.98
8:58:04 21 8:57:36 18.24 8:59:23 54.17 8:59:48 18
8:58:24 27 8:57:55 0 8:59:23 54.46
8:58:45 25.07 8:58:17 0
8:59:04 22.03 8:58:37 2.7
8:59:25 4.29 8:58:56 0
8:59:45 19.35 8:59:18 12.27
6-1. after step (5) from new LDID_SpeedMap searching loop obtain CPHM in LDID kspeed list speedList j, searching loop speedList jin all speed v iif, v i=0 (i.e. halt) and v i+s≠ 0 (0<s<=n-i, s are the number of continuous halt); If s>=3, think that this car has occurred that " stopping behavior " skips to 6-2 and judge whether to belong to " crossing stopping behavior ".
6-2. obtains the v of this stopping behavior icorresponding LDID and GPS_FDC i.Traversal intersectionMap, if there is this LDID in intersectionMap, illustrates that this LDID is section, crossing, enters 6-3 if do not exist.From intersectionMap, obtain crossing point (GPS_QD, GPS_ZD) that this LDID is corresponding, calculate GPS_FDC by formula (3) iget its minimum value L with the distance at crossing, if L<L_min, then think that this " stopping behavior " belongs to " crossing stopping behavior " and do not reject, otherwise enter 6-3 and judge whether to belong to " collective stops behavior ".
L=min(abs(GPS_FDC i-GPS_QD),abs(GPS_FDC i-GPS_ZD)(3)
6-3. obtains the v of this stopping behavior icorresponding CJSJ iand GPS_FDC iif the speed list of other vehicles occurs that " stopping behavior " obtains the GPS_FDC of its corresponding halt under this LDID of searching loop i'and CJSJ i', calculate Rule of judgment according to formula (4), formula (5),
abs(GPS_FDC i-GPS_FDC i')<Δs_min,i≠i'(4)
Abs (CJSJ i-CJSJ i') < Δ t_min, i ≠ i'(5) " collective stops behavior " does not reject if other vehicles all meet above Rule of judgment under this LDID, to think that this stopping behavior belonging to, otherwise this stopping behavior belonging to " really stopping behavior ", and rejects the speed record of continuous halt from speed list.
6-4. completes " really stopping behavior " identifies and returns a new LDID_SpeedMap after rejecting.
Step (7). after processing all data of this timeslice, obtain the speed record of all car speed lists under LDID according to LDID_SpeedMap and calculate average velocity corresponding to LDID according to formula (6), and applying for other stored in database Real-time Road velometer.
v a v g = &Sigma; k = 1 K &Sigma; i = 1 n k v i K &CenterDot; n k - - - ( 6 )
K is the vehicle number under this LDID, n kfor the speed record number of a kth car.
Update time, sheet, repeated step (2) to step (6).

Claims (5)

1., towards the floating car data parking behavior pattern cleaning method that road real-time speed calculates, it is characterized in that: described cleaning method comprises the steps:
Step (1). read map datum, map datum comprises section numbering LDID, some position, section information, map grid is obtained " grid-section numbering " mapping table;
Step (2). according to the floating car data of setting-up time sheet before " license plate number ", " creation-time " ascending order reading current time, floating car data comprises number-plate number CPHM, some position information GPS_FDC and creation-time CJSJ;
Step (3). calculate the Distance geometry mistiming of the double record of vehicle, then by Distance geometry mistiming computing velocity, and according to " grid-section numbering " mapping table, this speed is added in the section numbering at the information place, vehicle point position of second time record;
Step (4). the speed list of each vehicle on each section of this timeslice is obtained according to above step;
Step (5). identify and reject abending a little of causing due to " pseudo-stopping behavior " in speed list, and return speed list, process is as follows:
In the data structure LDID_SpeedMap that 5-1. obtains in step (4), searching loop obtains CPHM in LDID kspeed list speedList j, searching loop speedList jin all speed v iif, v i=0 and v i+s≠ 0,0<s<=n-i, s are the number of continuous halt, then the car speed record number δ=s+1 of " pseudo-stopping behavior " process;
5-2. calculates " pseudo-stopping behavior " car speed sudden change threshold value SPSSC according to formula (2); When i ≠ 1, the formula left side is the average velocity of all speed records before this speed list first halt, as i=1, with the average velocity v in a moment on this section tas average velocity threshold value;
S P S S C = &Sigma; f = 1 i - 1 v f i - 1 * &delta; , i &NotEqual; 1 v t * &delta; , i = 1 - - - ( 2 )
If 5-3. is v i+s>=SPSSC, then this speed record is the mutating speed that " pseudo-stopping behavior " causes, and rejects speed record and the mutating speed speed record of continuous halt from speed list;
5-4. completes " pseudo-stopping behavior " identifies and returns a new LDID_SpeedMap after rejecting.
2. as claimed in claim 1 a kind of towards road real-time speed calculate floating car data parking behavior pattern cleaning method, it is characterized in that: described cleaning method also comprises the steps:
Step (6). identify and reject abending a little of causing due to " really stopping behavior " in speed list, and return speed list, process is as follows:
6-1. after step (5) from new LDID_SpeedMap searching loop obtain CPHM in LDID kspeed list speedList j, searching loop speedList jin all speed v iif, v i=0 and v i+s≠ 0,0<s<=n-i, s are the number of continuous halt; If s>=3, think that this car has occurred that " stopping behavior " skips to 6-2 and judge whether to belong to " crossing stopping behavior ";
6-2. obtains the v of this stopping behavior icorresponding LDID and GPS_FDC i, traversal intersectionMap, if there is this LDID in intersectionMap, illustrates that this LDID is section, crossing, enters 6-3 if do not exist; From intersectionMap, obtain crossing point (GPS_QD, GPS_ZD) that this LDID is corresponding, calculate GPS_FDC by formula (3) iget its minimum value L with the distance at crossing, if L<L_min, then think that this " stopping behavior " belongs to " crossing stopping behavior " and do not reject, otherwise enter 6-3 and judge whether to belong to " collective stops behavior ";
L=min(abs(GPS_FDC i-GPS_QD),abs(GPS_FDC i-GPS_ZD)(3)
6-3. obtains the v of this stopping behavior icorresponding CJSJ iand GPS_FDC iif the speed list of other vehicles occurs that " stopping behavior " obtains the GPS_FDC of its corresponding halt under this LDID of searching loop i'and CJSJ i', calculate Rule of judgment according to formula (4), formula (5),
abs(GPS_FDC i-GPS_FDC i')<Δs_min,i≠i'(4)
abs(CJSJ i-CJSJ i')<Δt_min,i≠i'(5)
If other vehicles all meet above Rule of judgment under this LDID, " collective stops behavior " does not reject to think that this stopping behavior belonging to, otherwise this stopping behavior belonging to " really stopping behavior ", and rejects the speed record of continuous halt from speed list;
6-4. completes " really stopping behavior " identifies and returns a new LDID_SpeedMap after rejecting.
3. as claimed in claim 2 a kind of towards road real-time speed calculate floating car data parking behavior pattern cleaning method, it is characterized in that: described cleaning method also comprises the steps:
Step (7). after above step, calculate the average velocity of all speed records corresponding to LDID and be stored in database real-time speed table, update time, sheet repeated step (2) to step (6).
4. a kind of floating car data parking behavior pattern cleaning method calculated towards road real-time speed as described in one of claims 1 to 3, is characterized in that: in described step (3), velocity computation process is as follows:
3-1. obtain all records of each vehicle under this timeslice according to step (2), be designated as { CPHM k} _ { GPS_FDC i, CJSJ i, k ∈ Ki ∈ n, n is the number of all records, and K is the sum of this timeslice vehicle;
3-2. calculates car speed v according to formula (1.1) or (1.2) i, obtain car speed list { CPHM k} _ { v i, GPS_FDC i, CJSJ i.
What adopt is the absolute distance of the GPS point of double record, as the operating range of vehicle,
1) when not across using formula (1.1) when section
v i = a b s ( G P S _ FDC i - G P S _ FDC i - 1 ) CJSJ i - CJSJ i - 1 - - - ( 1.1 )
2) when across using formula (1.2) when section
v i = a b s ( G P S _ FDC i - G P S _ JD i ) + a b s ( G P S _ JD i - G P S _ FDC i - 1 ) CJSJ i - CJSJ i - 1 - - - ( 1.2 )
GPS_JD ifor the GPS point of double record, when across section, the some position information of the intersection point in two sections.
5. a kind of floating car data parking behavior pattern cleaning method calculated towards road real-time speed as described in one of claims 1 to 3, it is characterized in that: in described step (7), obtain the speed record of all car speed lists under LDID according to LDID_SpeedMap and calculate average velocity corresponding to LDID according to formula (6), and applying for other stored in database Real-time Road velometer.
v a v g = &Sigma; k = 1 K &Sigma; i = 1 n k v i K &CenterDot; n k - - - ( 6 )
K is the vehicle number under this LDID, n kfor the speed record number of a kth car.
CN201511031346.9A 2015-12-31 2015-12-31 The floating car data parking behavior pattern cleaning method calculated towards road real-time speed Active CN105575120B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511031346.9A CN105575120B (en) 2015-12-31 2015-12-31 The floating car data parking behavior pattern cleaning method calculated towards road real-time speed

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511031346.9A CN105575120B (en) 2015-12-31 2015-12-31 The floating car data parking behavior pattern cleaning method calculated towards road real-time speed

Publications (2)

Publication Number Publication Date
CN105575120A true CN105575120A (en) 2016-05-11
CN105575120B CN105575120B (en) 2017-11-07

Family

ID=55885195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511031346.9A Active CN105575120B (en) 2015-12-31 2015-12-31 The floating car data parking behavior pattern cleaning method calculated towards road real-time speed

Country Status (1)

Country Link
CN (1) CN105575120B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106548444A (en) * 2016-11-03 2017-03-29 杭州电子科技大学 Towards the Floating Car carrying behavior pattern cleaning method of road real-time speed estimation
CN106875680A (en) * 2017-03-21 2017-06-20 杭州电子科技大学 Crossing average latency computational methods based on big data analysis
CN108922193A (en) * 2018-08-03 2018-11-30 北京航空航天大学 A kind of intersection signal phase estimate method based on Floating Car track data
CN110766937A (en) * 2019-05-22 2020-02-07 北京嘀嘀无限科技发展有限公司 Parking spot identification method and device, electronic equipment and readable storage medium
CN114822016A (en) * 2022-03-11 2022-07-29 上海工程技术大学 Road section uniform speed estimation system based on GPS data flow

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001188986A (en) * 1999-12-30 2001-07-10 Hiroshi Tanaka Method for measuring vehicle travel speed and traffic volume on road by aerial photography
JP4809321B2 (en) * 2007-12-07 2011-11-09 住友電気工業株式会社 Traffic information calculation device, traffic information calculation program, and traffic information calculation method
CN102982678A (en) * 2012-11-28 2013-03-20 青岛海信网络科技股份有限公司 Traffic data information service system and method for realizing traffic data information service
CN103927873A (en) * 2014-04-28 2014-07-16 中国航天系统工程有限公司 Matching method for probe car and road section and method for obtaining real-time traffic status in parallel
CN103971516A (en) * 2014-05-13 2014-08-06 紫光捷通科技股份有限公司 Traffic data preprocessing method and road condition detection method
CN104282151A (en) * 2014-09-30 2015-01-14 北京交通发展研究中心 Real-time floating vehicle path matching method based on high-frequency satellite positioning data
CN105139645A (en) * 2015-07-23 2015-12-09 合肥革绿信息科技有限公司 Urban regional road network operation index assessment method based on floating car technology

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001188986A (en) * 1999-12-30 2001-07-10 Hiroshi Tanaka Method for measuring vehicle travel speed and traffic volume on road by aerial photography
JP4809321B2 (en) * 2007-12-07 2011-11-09 住友電気工業株式会社 Traffic information calculation device, traffic information calculation program, and traffic information calculation method
CN102982678A (en) * 2012-11-28 2013-03-20 青岛海信网络科技股份有限公司 Traffic data information service system and method for realizing traffic data information service
CN103927873A (en) * 2014-04-28 2014-07-16 中国航天系统工程有限公司 Matching method for probe car and road section and method for obtaining real-time traffic status in parallel
CN103971516A (en) * 2014-05-13 2014-08-06 紫光捷通科技股份有限公司 Traffic data preprocessing method and road condition detection method
CN104282151A (en) * 2014-09-30 2015-01-14 北京交通发展研究中心 Real-time floating vehicle path matching method based on high-frequency satellite positioning data
CN105139645A (en) * 2015-07-23 2015-12-09 合肥革绿信息科技有限公司 Urban regional road network operation index assessment method based on floating car technology

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106548444A (en) * 2016-11-03 2017-03-29 杭州电子科技大学 Towards the Floating Car carrying behavior pattern cleaning method of road real-time speed estimation
CN106875680A (en) * 2017-03-21 2017-06-20 杭州电子科技大学 Crossing average latency computational methods based on big data analysis
CN108922193A (en) * 2018-08-03 2018-11-30 北京航空航天大学 A kind of intersection signal phase estimate method based on Floating Car track data
CN108922193B (en) * 2018-08-03 2019-06-04 北京航空航天大学 A kind of intersection signal phase estimate method based on Floating Car track data
CN110766937A (en) * 2019-05-22 2020-02-07 北京嘀嘀无限科技发展有限公司 Parking spot identification method and device, electronic equipment and readable storage medium
CN110766937B (en) * 2019-05-22 2020-10-20 滴图(北京)科技有限公司 Parking spot identification method and device, electronic equipment and readable storage medium
CN114822016A (en) * 2022-03-11 2022-07-29 上海工程技术大学 Road section uniform speed estimation system based on GPS data flow
CN114822016B (en) * 2022-03-11 2023-07-25 上海工程技术大学 Road section uniform speed estimation system based on GPS data flow

Also Published As

Publication number Publication date
CN105575120B (en) 2017-11-07

Similar Documents

Publication Publication Date Title
CN102136190B (en) Dispatching management system and method for event emergency response of urban bus passenger transport
CN105575120A (en) Floating car data parking behavior mode cleaning method specific to road real time speed calculation
CN104933869B (en) Anti-noise type satellite positioning data-based method for identifying vehicle violation behaviors
CN108345666B (en) Vehicle abnormal track detection method based on time-space isolated points
CN109191849B (en) Traffic jam duration prediction method based on multi-source data feature extraction
CN104700646B (en) A kind of taxi exception track real-time detection method based on online gps data
CN113870570B (en) ETC-based road network operation state method, system and storage medium
CN103578295B (en) A kind of tail of expressway queue rear-end collision risk real-time estimate and prevention method
CN103646542B (en) A kind of Forecasting Methodology of Evaluating traffic impact area and device
CN106781468B (en) Link Travel Time Estimation method based on built environment and low frequency floating car data
CN103440768B (en) Dynamic-correction-based real-time bus arrival time predicting method
CN104240499A (en) Abnormal congestion point judgment method based on microwave data
CN104408920B (en) The method judging long-distance passenger transportation vehicular traffic violation behavior based on bayonet socket traffic information
CN104408925A (en) Array radar based intersection running state evaluation method
CN103971523A (en) Mountainous road traffic safety dynamic early-warning system
Goyani et al. Investigation of crossing conflicts by vehicle type at unsignalized T-intersections under varying roadway and traffic conditions in India
CN100470603C (en) Urban road traffic congestion detecting alarm system and operation method thereof
CN106297304A (en) A kind of based on MapReduce towards the fake-licensed car recognition methods of extensive bayonet socket data
CN103632541B (en) Traffic incident road chain detection and data filling method
CN102426783B (en) Low flow road traffic incident detection method based on vehicle tracking
CN109493449A (en) A kind of lorry loading method for estimating state based on lorry GPS track data and high speed transaction data
Singh et al. Driver behaviour modelling of vehicles at signalized intersection with heterogeneous traffic
CN105551241B (en) A kind of real-time jamming analysis method based on FCD and EP multi-source datas
CN105261218A (en) Floating car accompany behavior mode digging method based on large data analysis
CN105046958A (en) Highway traffic information acquisition node nonequidistance optimized layout method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 310012 floor 1, building 1, No. 223, Yile Road, Xihu District, Hangzhou City, Zhejiang Province

Patentee after: Yinjiang Technology Co.,Ltd.

Address before: 310012 floor 1, building 1, No. 223, Yile Road, Xihu District, Hangzhou City, Zhejiang Province

Patentee before: ENJOYOR Co.,Ltd.

CP01 Change in the name or title of a patent holder