CN106781479B - A method of highway operating status is obtained based on mobile phone signaling data in real time - Google Patents
A method of highway operating status is obtained based on mobile phone signaling data in real time Download PDFInfo
- Publication number
- CN106781479B CN106781479B CN201611205475.XA CN201611205475A CN106781479B CN 106781479 B CN106781479 B CN 106781479B CN 201611205475 A CN201611205475 A CN 201611205475A CN 106781479 B CN106781479 B CN 106781479B
- Authority
- CN
- China
- Prior art keywords
- user
- speed
- grid
- base station
- highway
- 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.)
- Active
Links
Classifications
-
- 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/012—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from other sources than vehicle or roadside beacons, e.g. mobile networks
-
- 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/0125—Traffic data processing
Abstract
The present invention relates to a kind of method for obtaining highway operating status in real time based on mobile phone signaling data, this method delimit highw ay m onitoring region first, and carry out rasterizing processing to monitoring region；Raster series in grid system are mapped to according to established grid system acquisition highway；Base station information in monitoring region, grid where obtaining each base station according to the location information of base station and grid system are then obtained, and the road grid that base station maps arrive is obtained according to base station range；Establish base station information table；Then the signaling data generated in monitoring region in monitoring time section is obtained, edit is carried out to signaling data and obtains the track sets of each user；Judge whether it travels on expressway in monitoring time section according to the user trajectory of monitoring time section and its historical track；The user for travelling and generating new track sets data in monitoring time section on expressway is obtained accordingly；Road grid real time execution speed is calculated followed by the signaling data that these users update.The present invention is not necessarily to additional hardware device, and its implementation is simple, and accuracy is higher.
Description
Technical field
The invention belongs to highway real-time road monitoring technical fields, are related to one kind and are obtained in real time based on mobile phone signaling data
The method for taking highway operating status.
Background technique
With being continuously increased for China's vehicle fleet size, highway vehicle flowrate constantly increases, and leads to expressway traffic accident frequency
Hair, congestion occurrence frequency are higher and higher.The acquisition of traffic state information is mainly using fixed sensor at present.As tachymeter,
Vehicle checker, camera etc..And the cost of these hardware facilities, installation and the expense of maintenance are relatively high, and its monitoring range has
Limit, cannot cover whole highway.Therefore its urgent need one kind can monitor entire highway network, and lower-cost
Method.
Mobile phone signaling be it is a kind of result from mobile communications network and can reflect that this communication behavior occurs for mobile phone user where
The data of geographical location and time.Such data contain mobile phone user's travel behaviour feature, can be not only used for analysis and use
The trip mode and trip preference at family, moreover it is possible to by analyzing user's sample of group, obtain the real-time running state of highway.
Summary of the invention
Highway operation is obtained based on mobile phone signaling data in real time in view of this, the purpose of the present invention is to provide one kind
The method of state, this method analyze the signaling data acquired in real time using the relevant technologies, identify the user travelled on expressway；
And the signaling data analytical calculation generated in monitoring cycle using user obtains express highway section average running speed, according to
The road-section average speed of service obtains highway real-time road.
In order to achieve the above objectives, the invention provides the following technical scheme:
A method of highway operating status being obtained based on mobile phone signaling data in real time, this method includes following step
It is rapid:
Step 1: according to highway actual path, delimiting rectangle and monitor region；It is the square grid of L using length and width
Rasterizing processing is carried out to monitoring region；To monitor the region lower left corner as starting point, using geographical horizontal line as horizontal axis, vertical line is the longitudinal axis
Divide grid；Grid is expressed as B [x, y], wherein x line number where grid thus, y columns where grid thus；If highway
A certain driving direction, if the grid containing express highway section is road grid, thus obtains highway forward direction road to be positive
Road raster series L_HB={ B [h1,w1],B[h2,w2],B[h3,w3],…,B[hi,wi],…,B[hn,wn], wherein B [hi,wi]
For a road grid, hiLine number where grid thus, wiColumns where grid thus；
Step 2: obtaining the base station data in monitoring region, and the grating map determined according to step 1, establish base station information
Table；Base station information table (is base comprising base station location area number LAC, base station cell number CellID, base station geographic position information
Stand latitude and longitude coordinates, referred to as LOC), affiliated grid (referred to as CB), whether cover highway (referred to as IOW), be corresponding
Road grid (referred to as CWB)；
Step 3: set T moment user information collection and be combined into S_U, wherein user information include user's unique identification UID and this
Track sets of the user before the T moment；It obtains and monitors all signaling datas that region generates in T to T+t monitoring time section,
Wherein signaling data includes user's unique identification UID, signaling data own base station position area number LAC, base station cell number
CellID and time stamp T imeStamp field；Then by signaling data by user classify and to the data of each user temporally
Stamp TimeStamp is ranked up, and filters pingpang handoff, and user's data co-located in continuous time period are carried out
Merge, to construct track sets of the user in monitoring time section；User information set S_U is then updated according to UID, if S_
There are this users in U, then update the track sets of this user；If this user is not present in S_U, this user is added into S_U
Thus information obtains the user information set S_U1 that track data updates in monitoring time section；
Step 4: being located at the user travelled on expressway is highway user, obtains T moment highway user set
S_HU (wherein S_HU is contained in S_U)；Effectively trip section segmentation is carried out to the track sets of user each in S_U1, and analyzes it most
The track sets of latter section of section of effectively going on a journey, judge whether it travels on expressway, and updating maintenance highway is used accordingly
Family set S_HU and generate calculate each road grid real-time speed of highway needed for positive track sets set S_PE and negative
To track sets set S_NE；
Step 5: extracting each track sets in S_PE and S_NE, the positive negative sense that analytical calculation obtains each road grid is fast in real time
Degree；Wherein road grid real-time speed is the express highway section real-time speed that road grid includes, wherein positive negative direction road
Road grid real-time speed calculation method is consistent.
Further, in step 2, the affiliated grid CB in base station is determined according to base station LOC first, model is then covered according to base station
It encloses, determines the grid set of base station covering；If containing road grid in the grid set of base station covering, this base station IOW is belonged to
Property value is set as iw, this IOW attribute value is otherwise set as ow；Wherein iw indicates that this base station range contains highway, ow table
Show that this base station cannot cover highway；If base station IOW attribute value is iw, the road grid of calculation base station covering and its place
The distance of grid CB, taking apart from the smallest road grid is the grid CWB in base station maps to road；If base station IOW attribute value
For ow, then its CWB value is null；If grid where certain base station is B [a, b], the road grid set of covering is S_CHB=
{B[c1,d1],B[c2,d2],B[c3,d3],…,B[ci,di],…,B[cm,dm]}；Calculate each grid and base station place in S_CHB
Distance S_D={ the D of grid B [a, b]c1d1,Dc2d2,Dc3d3,…,Dcidi,…,Dcmdm(wherein DcidiFor B [ci,di] and B [a, b]
The distance between, Dcidi=((ci-a)2+(di-b)2)(1/2)), choosing the smallest grid of S_D intermediate value is that base station maps arrive road
Grid.
Further, in step 4, it specifically includes:
Step 4.1: extracting the track sets of each user in S_U1, and carry out effectively trip section segmentation；It extracts and is used in S_U1
Family final stage is effectively gone on a journey section, is effectively jumped to step 4.2 if trip section signaling number is greater than certain threshold value if this, is otherwise judged
It whether there is this user in S_HU, this user if it exists, then according to user's driving direction by its track in monitoring time section
Sequence is added in corresponding track sets set；
Step 4.2: obtaining effectively trip section track sets L_P={ P1,P2,P3,…,Pi,…,Pu(wherein PiIndicate rail
I-th signaling data in mark sequence, in chronological sequence sequence arranges between each signaling), it is obtained according to the base station information table of step 2
Thus the corresponding road grid CWB in the base station Pi obtains the corresponding road raster series L_PCWB={ CWB of track sets1,CWB2,
CWB3,…,CWBi,…,CWBu(wherein CWBiFor PiCorresponding CWB, CWBiIt may be null)；
Step 4.3: calculating the accounting of the CWB of non-null in L_PCWB sequence；Obtaining value in L_PCWB sequence is null
CWB number N1 and the non-null of value CWB number N2；Calculate the accounting R1=N2/ (N1+ of the CWB of non-null
N2)；Step 4.4 is gone to if R1 is greater than certain threshold value, this user non-freeway user is otherwise assert, gos to step 4.5；
Step 4.4: extracting value in L_PCWB is that the CWB of non-null forms the road grid of the corresponding non-empty of track sets P
Lattice sequence L_PCWB2={ CWBg1,CWBg2,CWBg3,…,CWBgi,…,CWBgq(wherein L_PCWB2 is contained in L_PCWB, gi∈
[1,u])；Extract CWBgiSerial number O in Expressway Road raster series L_HBgi, then the corresponding sequence number sequence of L_PCWB2
For L_O={ Og1,Og2,Og3,…,Ogi,…,Ogq}；Calculate the average value Agi=(O of two neighboring element in L_Ogi+Og(i+1))/
2, then the sequence after L_O smoothing processing is L_OA={ Ag1,Ag2,Ag3,…,Agi,…,Ag(q-1)}；If there are bright for L_OA sequence
Aobvious growth trend then assert that the moving direction of this user is matched with highway and it is highway user；Become if increasing
Gesture is positive, then it is assumed that this user's moving direction is highway positive direction, and its track sets in monitoring time section is added
It is added in S_PE；If growth trend is negative, then it is assumed that this user's moving direction is highway negative direction, and by it in monitoring
Between track sets in section be added in S_NE；
Step 4.5: updating highway user set S_HU；It obtains T moment highway user collection and is combined into S_HU, if this
User is highway user and this user is not present in S_HU, then adds this user into S_HU；If this user is judged as non-height
There are this users in fast highway user and S_HU, then this user is deleted from S_HU.
Further, in steps of 5, by taking positive direction speed is sought as an example, step includes:
Step 5.1: corresponding speed container being set for road grid each in L_HB, for storing mobile phone in each track sets
Movement speed of the user between two tracing points；Speed is divided into Z grade according to the path length between two tracing points, then often
It include Z speed container in the speed container of a road grid, each speed container stores corresponding speed data；
Step 5.2: movement speed of the user between two tracing points in each track sets in calculating S_PE, according to two tracks
Distance between point carries out grade classification to the speed being calculated, and puts it into the corresponding speed of corresponding road grid and hold
In device；
Step 5.3: the speed data in each speed container of road grid each in L_HB being gathered using clustering algorithm
Alanysis obtains the corresponding speed of each road grid difference distance grade；If road grid B [hi,wi] corresponding speed sets of containers
It is combined into S_B [hi,wi] _ Con={ B [hi,wi]_C1,B[hi,wi]_C2,B[hi,wi]_C3,…,B[hi,wi]_Cj,…,B[hi,
wi]_Cz, B [h is calculated using focusing solutions analysisi,wi]_CjIt is Con_ that middle speed data, which obtains the corresponding speed of grade j,
Speedj, then road grid B [hi,wi] different brackets sets of speeds be S_Con_Speed={ Con_Speed1,Con_Speed2,
Con_Speed3,…,Con_Speedj,…,Con_Speedz}；Speed data sample number in each speed container is obtained, then road
Grid B [hi,wi] different brackets speed sample number collection is combined into S_Con_NS={ NS1,NS2,NS3,…,NSj,…NSz}；Then road
Grid B [hi,wi] speed total sample number B [hi,wi] _ SpeedSwatch_Sum=NS1+NS2+NS3+…+NSj+…+NSz；It utilizes
Following formula calculates road grid B [hi,wi] real-time speed B [hi,wi] _ RTSpeed:
Step 5.4: the traffic behavior of express highway section therein is determined according to the real time execution speed of road grid.
Further, in the step 5.2, the specific steps are as follows:
Step 5.2.1: a certain track sets L_EP={ R in S_PE is obtained1,R2,R3,…,Ri,…,Rk(wherein RiIt indicates
I-th of tracing point in this track sets)；Each tracing point in order traversal L_EP, calculate mobile phone user itself and subsequent each point it
Between movement speed, and the speed being calculated is included into the speed container of corresponding road grid；
Step 5.2.2: tracing point R is seti-1And RiCorresponding road grid is respectively B [e(i-1),f(i-1)] and B [ei,fi],
Obtain B [e(i-1),f(i-1)] and B [ei,fi] serial number in L_HB, if the difference of two serial numbers and the moving direction of this user are positive together
Or with being negative, then jump to 5.2.3；
Step 5.2.3: B [e is calculated(i-1),f(i-1)] and B [ei,fi] between distance Dis(i-1)i；Obtain B in L_HB
[e(i-1),f(i-1)] and B [ei,fi] between subsequence L_EHB={ B [s1,l1],B[s2,l2],…,B[si,li],…,B[sv,
lv], then B [e(i-1),f(i-1)] and B [ei,fi] between distance Dis(i-1)iIt can be calculated with following formula；
Step 5.2.4: tracing point R is seti-1And RiCorresponding time stamp T imeStamp(i-1)And TimeStampi, then user
In Ri-1And RiBetween movement speed Speed(i-1)i=Dis(i-1)i/(TimeStampi-TimeStamp(i-1))；
Step 5.2.5: the distance Dis obtained according to step 5.2.3(i-1)iBy Speed(i-1)iIt is put into B [e(i-1),f(i-1)]
With B [ei,fi] between subsequence L_EHB in the corresponding speed container of each road grid in；If Speed(i-1)iIt is corresponding
Grade is Level_Speed(i-1)i, rate sequence is L_Level={ Level1,Level2,Level3,…,Levelj,…,
Levelz(wherein LeveljIndicate j-th of speed class)；Then Level_Speed can be determined by following formula(i-1)i；If being
Dis(i-1)iMore than or equal to DisNumjAnd it is less than DisNumj+1, then Level_Speed(i-1)iFor Levelj；Wherein DisNumj+1Greatly
In DisNumjIf Dis(i-1)iNot in any section, then this speed data is given up；
The beneficial effects of the present invention are: the signaling data that the present invention is acquired in real time by analysis identifies on expressway
The user of traveling；And the signaling data analytical calculation generated in monitoring cycle using user is obtained express highway section and averagely transported
Scanning frequency degree obtains highway real-time road according to the road-section average speed of service；The present invention is not necessarily to additional hardware device, and its
Implementation method is simple, and accuracy is higher.
Detailed description of the invention
In order to keep the purpose of the present invention, technical scheme and beneficial effects clearer, the present invention provides following attached drawing and carries out
Illustrate:
Fig. 1 is a kind of method for obtaining highway operating status in real time based on mobile phone signaling data provided by the invention
Process.
Specific embodiment
Below in conjunction with attached drawing, a preferred embodiment of the present invention will be described in detail.
Fig. 1 is a kind of method for obtaining highway operating status in real time based on mobile phone signaling data provided by the invention
Process, as shown, this method specifically includes the following steps:
Step 1: according to highway actual path, delimiting rectangle and monitor region.It is the square grid of L using length and width
Rasterizing processing is carried out to monitoring region.To monitor the region lower left corner as starting point, using geographical horizontal line as horizontal axis, vertical line is the longitudinal axis
Divide grid.Grid is expressed as B [x, y], wherein x line number where grid thus, y columns where grid thus.If highway
A certain driving direction, if the grid containing express highway section is road grid, thus obtains highway forward direction road to be positive
Road raster series L_HB={ B [h1,w1],B[h2,w2],B[h3,w3],…,B[hi,wi],…,B[hn,wn], wherein B [hi,wi]
For a road grid, hiLine number where grid thus, wiColumns where grid thus.
Step 2: obtaining the base station data in monitoring region, and the grating map determined according to step 1, establish base station information
Table.Base station information table (is base comprising base station location area number LAC, base station cell number CellID, base station geographic position information
Stand latitude and longitude coordinates, referred to as LOC), affiliated grid (referred to as CB), whether cover highway (referred to as IOW), be corresponding
Road grid (referred to as CWB).The affiliated grid CB in base station is determined according to base station LOC first, then according to base station range, really
Determine the grid set of base station covering.If containing road grid in the grid set of base station covering, this base station IOW attribute value is set
For iw, this IOW attribute value is otherwise set as ow.Wherein iw indicates that this base station range contains highway, and ow indicates this base
Highway cannot be covered by standing.If base station IOW attribute value is iw, the road grid and grid CB where it of calculation base station covering
Distance, take apart from the smallest road grid be base station maps to road on grid CWB.If base station IOW attribute value is ow,
Its CWB value is null.If grid where certain base station is B [a, b], the road grid set of covering is S_CHB={ B [c1,
d1],B[c2,d2],B[c3,d3],…,B[ci,di],…,B[cm,dm]}.Calculate each grid and base station place grid B in S_CHB
Distance S_D={ the D of [a, b]c1d1,Dc2d2,Dc3d3,…,Dcidi,…,Dcmdm(wherein DcidiFor B [ci,di] with B [a, b] between
Distance, Dcidi=((ci-a)2+(di-b)2)(1/2)), choosing the smallest grid of S_D intermediate value is grid of the base station maps to road
Lattice.
Step 3: set T moment user information collection and be combined into S_U, wherein user information include user's unique identification UID and this
Track sets of the user before the T moment.It obtains and monitors all signaling datas that region generates in T to T+t monitoring time section,
Wherein signaling data includes user's unique identification UID, signaling data own base station position area number LAC, base station cell number
CellID and time stamp T imeStamp field.Then by signaling data by user classify and to the data of each user temporally
Stamp TimeStamp is ranked up, and filters pingpang handoff, and user's data co-located in continuous time period are carried out
Merge, to construct track sets of the user in monitoring time section.User information set S_U is then updated according to UID, if S_
There are this users in U, then update the track sets of this user.If this user is not present in S_U, this user is added into S_U
Information.Thus the user information set S_U1 that track data updates in monitoring time section is obtained.
Step 4: being located at the user travelled on expressway is highway user, obtains T moment highway user set
S_HU (wherein S_HU is contained in S_U).Effectively trip section segmentation is carried out to the track sets of user each in S_U1, and analyzes it most
The track sets of latter section of section of effectively going on a journey, judge whether it travels on expressway, and updating maintenance highway is used accordingly
Family set S_HU and generate calculate each road grid real-time speed of highway needed for positive track sets set S_PE and negative
To track sets set S_NE.Specific step is as follows:
Step 4.1: extracting the track sets of each user in S_U1, and carry out effectively trip section segmentation.It extracts and is used in S_U1
Family final stage is effectively gone on a journey section, is effectively jumped to step 4.2 if trip section signaling number is greater than certain threshold value if this, is otherwise judged
It whether there is this user in S_HU, this user if it exists, then according to user's driving direction by its track in monitoring time section
Sequence is added in corresponding track sets set.
Step 4.2: obtaining effectively trip section track sets L_P={ P1,P2,P3,…,Pi,…,Pu(wherein PiIndicate rail
I-th signaling data in mark sequence, in chronological sequence sequence arranges between each signaling), it is obtained according to the base station information table of step 2
Thus the corresponding road grid CWB in the base station Pi obtains the corresponding road raster series L_PCWB={ CWB of track sets1,CWB2,
CWB3,…,CWBi,…,CWBu(wherein CWBiFor PiCorresponding CWB, CWBiIt may be null).
Step 4.3: calculating the accounting of the CWB of non-null in L_PCWB sequence.Obtaining value in L_PCWB sequence is null
CWB number N1 and the non-null of value CWB number N2.Calculate the accounting R1=N2/ (N1+ of the CWB of non-null
N2).Step 4.4 is gone to if R1 is greater than certain threshold value, this user non-freeway user is otherwise assert, gos to step 4.5.
Step 4.4: extracting value in L_PCWB is that the CWB of non-null forms the road grid of the corresponding non-empty of track sets P
Lattice sequence L_PCWB2={ CWBg1,CWBg2,CWBg3,…,CWBgi,…,CWBgq(wherein L_PCWB2 is contained in L_PCWB, gi∈
[1,u]).Extract CWBgiSerial number O in Expressway Road raster series L_HBgi, then the corresponding sequence number sequence of L_PCWB2
For L_O={ Og1,Og2,Og3,…,Ogi,…,Ogq}.Calculate the average value Agi=(O of two neighboring element in L_Ogi+Og(i+1))/
2, then the sequence after L_O smoothing processing is L_OA={ Ag1,Ag2,Ag3,…,Agi,…,Ag(q-1)}.If there are bright for L_OA sequence
Aobvious growth trend then assert that the moving direction of this user is matched with highway and it is highway user.Become if increasing
Gesture is positive, then it is assumed that this user's moving direction is highway positive direction, and its track sets in monitoring time section is added
It is added in S_PE.If growth trend is negative, then it is assumed that this user's moving direction is highway negative direction, and by it in monitoring
Between track sets in section be added in S_NE.
Step 4.5: updating highway user set S_HU.It obtains T moment highway user collection and is combined into S_HU, if this
User is highway user and this user is not present in S_HU, then adds this user into S_HU.If this user is judged as non-height
There are this users in fast highway user and S_HU, then this user is deleted from S_HU.
Step 5: extracting each track sets in S_PE and S_NE, the positive negative sense that analytical calculation obtains each road grid is fast in real time
Degree.Wherein road grid real-time speed is the express highway section real-time speed that road grid includes.Wherein positive negative direction road
Road grid real-time speed calculation method is consistent.Below by taking the calculating of positive direction road grid real-time speed as an example, specific steps are such as
Under:
Step 5.1: corresponding speed container being set for road grid each in L_HB, for storing mobile phone in each track sets
Movement speed of the user between two tracing points.And speed is divided by Z grade according to the path length between two tracing points again,
It then include Z speed container in the speed container of each road grid, each speed container stores corresponding speed data.
Step 5.2: movement speed of the user between two tracing points in each track sets in calculating S_PE, according to two tracks
Distance between point carries out grade classification to the speed being calculated, and puts it into the corresponding speed of corresponding road grid and hold
In device.
Specific step is as follows:
Step 5.2.1: a certain track sets L_EP={ R in S_PE is obtained1,R2,R3,…,Ri,…,Rk(wherein RiIt indicates
I-th of tracing point in this track sets).Each tracing point in order traversal L_EP, calculate mobile phone user itself and subsequent each point it
Between movement speed, and the speed being calculated is included into the speed container of corresponding road grid.
Step 5.2.2: tracing point R is seti-1And RiCorresponding road grid is respectively B [e(i-1),f(i-1)] and B [ei,fi],
Obtain B [e(i-1),f(i-1)] and B [ei,fi] serial number in L_HB, if the difference of two serial numbers and the moving direction of this user are positive together
Or with being negative, then jump to 5.2.3.
Step 5.2.3: B [e is calculated(i-1),f(i-1)] and B [ei,fi] between distance Dis(i-1)i.Obtain B in L_HB
[e(i-1),f(i-1)] and B [ei,fi] between subsequence L_EHB={ B [s1,l1],B[s2,l2],…,B[si,li],…,B[sv,
lv], then B [e(i-1),f(i-1)] and B [ei,fi] between distance Dis(i-1)iIt can be calculated with following formula.
Step 5.2.4: tracing point R is seti-1And RiCorresponding time stamp T imeStamp(i-1)And TimeStampi, then user
In Ri-1And RiBetween movement speed Speed(i-1)i=Dis(i-1)i/(TimeStampi-TimeStamp(i-1))。
Step 5.2.5: the distance Dis obtained according to step 5.2.3(i-1)iBy Speed(i-1)iIt is put into B [e(i-1),f(i-1)]
With B [ei,fi] between subsequence L_EHB in the corresponding speed container of each road grid in.If Speed(i-1)iIt is corresponding
Grade is Level_Speed(i-1)i, rate sequence is L_Level={ Level1,Level2,Level3,…,Levelj,…,
Levelz(wherein LeveljIndicate j-th of speed class).Then Level_Speed can be determined by following formula(i-1)i.If being
Dis(i-1)iMore than or equal to DisNumjAnd it is less than DisNumj+1, then Level_Speed(i-1)iFor Levelj.Wherein DisNumj+1Greatly
In DisNumjIf Dis(i-1)iNot in any section, then this speed data is given up.
Step 5.3: the speed data in each speed container of road grid each in L_HB being gathered using clustering algorithm
Alanysis obtains the corresponding speed of each road grid difference distance grade.If road grid B [hi,wi] corresponding speed sets of containers
It is combined into S_B [hi,wi] _ Con={ B [hi,wi]_C1,B[hi,wi]_C2,B[hi,wi]_C3,…,B[hi,wi]_Cj,…,B[hi,
wi]_Cz, B [h is calculated using focusing solutions analysisi,wi]_CjIt is Con_ that middle speed data, which obtains the corresponding speed of grade j,
Speedj, then road grid B [hi,wi] different brackets sets of speeds are as follows:
S_Con_Speed={ Con_Speed1,Con_Speed2,Con_Speed3,…,Con_Speedj,…,Con_
Speedz}.Speed data sample number in each speed container is obtained, then road grid B [hi,wi] different brackets speed sample manifold
It is combined into S_Con_NS={ NS1,NS2,NS3,…,NSj,…NSz}.Then road grid B [hi,wi] speed total sample number B [hi,wi]_
SpeedSwatch_Sum=NS1+NS2+NS3+…+NSj+…+NSz.Road grid B [h is calculated using following formulai,wi] in real time
Speed B [hi,wi]_RTSpeed。
Step 5.4: the traffic behavior of express highway section therein is determined according to the real time execution speed of road grid.
Finally, it is stated that preferred embodiment above is only used to illustrate the technical scheme of the present invention and not to limit it, although logical
It crosses above preferred embodiment the present invention is described in detail, however, those skilled in the art should understand that, can be
Various changes are made to it in form and in details, without departing from claims of the present invention limited range.
Claims (5)
1. a kind of method for obtaining highway operating status in real time based on mobile phone signaling data, it is characterised in that: this method packet
Include following steps:
Step 1: according to highway actual path, delimiting rectangle and monitor region；It is the square grid of L to prison using length and width
It surveys region and carries out rasterizing processing；To monitor the region lower left corner as starting point, using geographical horizontal line as horizontal axis, vertical line is longitudinal axis division
Grid；Grid is expressed as B [x, y], wherein x line number where grid thus, y columns where grid thus；If highway is a certain
Driving direction, if the grid containing express highway section is road grid, thus obtains highway forward direction road grid to be positive
Lattice sequence L_HB={ B [h1,w1],B[h2,w2],B[h3,w3],…,B[hi,wi],…,B[hn,wn], wherein B [hi,wi] it is one
A road grid, hiLine number where grid thus, wiColumns where grid thus；
Step 2: obtaining the base station data in monitoring region, and the grating map determined according to step 1, establish base station information table；
Base station information table includes base station location area number LAC, base station cell number CellID, latitude and longitude of base station coordinate LOC, affiliated grid
CB, highway IOW, corresponding road grid CWB whether are covered；
Step 3: setting T moment user information collection and be combined into S_U, wherein user information includes user's unique identification UID and this user
Track sets before the T moment；It obtains and monitors all signaling datas that region generates in T to T+t monitoring time section, wherein
Signaling data includes user's unique identification UID, signaling data own base station position area number LAC, base station cell number CellID
With time stamp T imeStamp field；Then signaling data is classified by user and the data of each user is temporally stabbed
TimeStamp is ranked up, and filters pingpang handoff, and user's data co-located in continuous time period are closed
And to it construct track sets of the user in monitoring time section；User information set S_U is then updated according to UID, if S_U
Middle there are this users, then update the track sets of this user；If this user is not present in S_U, this user letter is added into S_U
Thus breath obtains the user information set S_U1 that track data updates in monitoring time section；
Step 4: being located at the user travelled on expressway is highway user, obtains T moment highway user set S_HU,
Wherein S_HU is contained in S_U；Effectively trip section segmentation is carried out to the track sets of user each in S_U1, and analyzes its final stage
The effectively track sets of trip section, judge whether it travels on expressway, and updating maintenance highway user set accordingly
Positive track sets set S_PE and negative sense track needed for S_HU and generation calculate each road grid real-time speed of highway
Arrangement set S_NE；
Step 5: extracting each track sets in S_PE and S_NE, analytical calculation obtains the positive negative sense real-time speed of each road grid；
Wherein road grid real-time speed is the express highway section real-time speed that road grid includes, wherein positive negative direction road grid
Lattice real-time speed calculation method is consistent.
2. a kind of method that highway operating status is obtained based on mobile phone signaling data in real time according to claim 1,
It is characterized by: in step 2, determine the affiliated grid CB in base station according to base station LOC first, then according to base station range,
Determine the grid set of base station covering；If containing road grid in the grid set of base station covering, by this base station IOW attribute value
It is set as iw, this IOW attribute value is otherwise set as ow；Wherein iw indicates that this base station range contains highway, and ow indicates this
Base station cannot cover highway；If base station IOW attribute value is iw, the road grid and grid where it of calculation base station covering
The distance of CB, taking apart from the smallest road grid is the grid CWB in base station maps to road；If base station IOW attribute value is ow,
Then its CWB value is null；If grid where certain base station is B [a, b], the road grid set of covering is S_CHB={ B
[c1,d1],B[c2,d2],B[c3,d3],…,B[ci,di],…,B[cm,dm]}；Calculate each grid and base station place grid in S_CHB
Distance S_D={ the D of lattice B [a, b]c1d1,Dc2d2,Dc3d3,…,Dcidi,…,Dcmdm, wherein DcidiFor B [ci,di] with B [a, b] it
Between distance, Dcidi=((ci-a)2+(di-b)2)(1/2), choosing the smallest grid of S_D intermediate value is grid of the base station maps to road
Lattice.
3. a kind of method that highway operating status is obtained based on mobile phone signaling data in real time according to claim 2,
It is characterized by: in step 4, specifically including:
Step 4.1: extracting the track sets of each user in S_U1, and carry out effectively trip section segmentation；Extract S_U1 in user most
Latter section of section of effectively going on a journey effectively jumps to step 4.2 if trip section signaling number is greater than certain threshold value if this, otherwise judges S_HU
In whether there is this user, this user if it exists, then according to user's driving direction by its track sets in monitoring time section
It is added in corresponding track sets set；
Step 4.2: obtaining effectively trip section track sets L_P={ P1,P2,P3,…,Pi,…,Pu, wherein PiIndicate track sequence
I-th signaling data in column, in chronological sequence sequence arranges between each signaling, obtains Pi base according to the base station information table of step 2
Stand corresponding road grid CWB, thus obtains the corresponding road raster series L_PCWB={ CWB of track sets1,CWB2,
CWB3,…,CWBi,…,CWBu, wherein CWBiFor PiCorresponding CWB, CWBiIt may be null；
Step 4.3: calculating the accounting of the CWB of non-null in L_PCWB sequence；Obtain the CWB that value in L_PCWB sequence is null
Number N1 and the non-null of value CWB number N2；Calculate the accounting R1=N2/ (N1+N2) of the CWB of non-null；If R1
Step 4.4 is then gone to greater than certain threshold value, otherwise assert this user non-freeway user, gos to step 4.5；
Step 4.4: extracting value in L_PCWB is that the CWB of non-null forms the road grid sequence of the corresponding non-empty of track sets P
Arrange L_PCWB2={ CWBg1,CWBg2,CWBg3,…,CWBgi,…,CWBgq, wherein L_PCWB2 is contained in L_PCWB, gi∈[1,
u]；Extract CWBgiSerial number O in Expressway Road raster series L_HBgi, then the corresponding sequence number sequence of L_PCWB2 is L_O
={ Og1,Og2,Og3,…,Ogi,…,Ogq}；Calculate the average value Agi=(O of two neighboring element in L_Ogi+Og(i+1))/2, then
Sequence after L_O smoothing processing is L_OA={ Ag1,Ag2,Ag3,…,Agi,…,Ag(q-1)}；If L_OA sequence exists apparent
Growth trend then assert that the moving direction of this user is matched with highway and it is highway user；If growth trend is
Just, then it is assumed that this user's moving direction is highway positive direction, and its track sets in monitoring time section is added to
In S_PE；If growth trend is negative, then it is assumed that this user's moving direction is highway negative direction, and by it in monitoring time section
Interior track sets are added in S_NE；
Step 4.5: updating highway user set S_HU；It obtains T moment highway user collection and is combined into S_HU, if this user
For highway user and this user is not present in S_HU, then adds this user into S_HU；If this user is judged as non-high-speed public affairs
There are this users in road user and S_HU, then this user is deleted from S_HU.
4. a kind of method that highway operating status is obtained based on mobile phone signaling data in real time according to claim 3,
It is characterized by: in steps of 5, by taking positive direction speed is sought as an example, step includes:
Step 5.1: corresponding speed container being set for road grid each in L_HB, for storing mobile phone user in each track sets
Movement speed between two tracing points；Speed is divided into Z grade according to the path length between two tracing points, then each road
It include Z speed container in the speed container of road grid, each speed container stores corresponding speed data；
Step 5.2: calculate in S_PE movement speed of the user between two tracing points in each track sets, according to two tracing points it
Between distance grade classification is carried out to the speed that is calculated, and put it into the corresponding speed container of corresponding road grid
In；
Step 5.3: cluster point being carried out to the speed data in each speed container of road grid each in L_HB using clustering algorithm
Analysis, obtains the corresponding speed of each road grid difference distance grade；If road grid B [hi,wi] corresponding speed container set is
S_B[hi,wi] _ Con={ B [hi,wi]_C1,B[hi,wi]_C2,B[hi,wi]_C3,…,B[hi,wi]_Cj,…,B[hi,wi]_
Cz, B [h is calculated using focusing solutions analysisi,wi]_CjIt is Con_Speed that middle speed data, which obtains the corresponding speed of grade j,j, then
Road grid B [hi,wi] different brackets sets of speeds be S_Con_Speed={ Con_Speed1,Con_Speed2,Con_
Speed3,…,Con_Speedj,…,Con_Speedz}；Speed data sample number in each speed container is obtained, then road grid B
[hi,wi] different brackets speed sample number collection is combined into S_Con_NS={ NS1,NS2,NS3,…,NSj,…NSz}；Then road grid B
[hi,wi] speed total sample number B [hi,wi] _ SpeedSwatch_Sum=NS1+NS2+NS3+…+NSj+…+NSz；Using following
Formula calculates road grid B [hi,wi] real-time speed B [hi,wi] _ RTSpeed:
Step 5.4: the traffic behavior of express highway section therein is determined according to the real time execution speed of road grid.
5. a kind of method that highway operating status is obtained based on mobile phone signaling data in real time according to claim 4,
It is characterized by: in the step 5.2, the specific steps are as follows:
Step 5.2.1: a certain track sets L_EP={ R in S_PE is obtained1,R2,R3,…,Ri,…,Rk, wherein RiIndicate this rail
I-th of tracing point in mark sequence；Each tracing point in order traversal L_EP calculates mobile phone user in its shifting between subsequent each point
Dynamic speed, and the speed being calculated is included into the speed container of corresponding road grid；
Step 5.2.2: tracing point R is seti-1And RiCorresponding road grid is respectively B [e(i-1),f(i-1)] and B [ei,fi], obtain B
[e(i-1),f(i-1)] and B [ei,fi] serial number in L_HB, if the moving direction of the difference of two serial numbers and this user are the same as being positive or together
It is negative, then jumps to 5.2.3；
Step 5.2.3: B [e is calculated(i-1),f(i-1)] and B [ei,fi] between distance Dis(i-1)i；Obtain B [e in L_HB(i-1),
f(i-1)] and B [ei,fi] between subsequence L_EHB={ B [s1,l1],B[s2,l2],…,B[si,li],…,B[sv,lv], then
B[e(i-1),f(i-1)] and B [ei,fi] between distance Dis(i-1)iIt can be calculated with following formula；
Step 5.2.4: tracing point R is seti-1And RiCorresponding time stamp T imeStamp(i-1)And TimeStampi, then user is in Ri-1
And RiBetween movement speed Speed(i-1)i=Dis(i-1)i/(TimeStampi-TimeStamp(i-1))；
Step 5.2.5: the distance Dis obtained according to step 5.2.3(i-1)iBy Speed(i-1)iIt is put into B [e(i-1),f(i-1)] and B
[ei,fi] between subsequence L_EHB in the corresponding speed container of each road grid in；If Speed(i-1)iCorresponding grade
For Level_Speed(i-1)i, rate sequence is L_Level={ Level1,Level2,Level3,…,Levelj,…,
Levelz, wherein LeveljIndicate j-th of speed class；Then Level_Speed can be determined by following formula(i-1)i；If being
Dis(i-1)iMore than or equal to DisNumjAnd it is less than DisNumj+1, then Level_Speed(i-1)iFor Levelj；Wherein DisNumj+1Greatly
In DisNumjIf Dis(i-1)iNot in any section, then this speed data is given up；
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611205475.XA CN106781479B (en) | 2016-12-23 | 2016-12-23 | A method of highway operating status is obtained based on mobile phone signaling data in real time |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611205475.XA CN106781479B (en) | 2016-12-23 | 2016-12-23 | A method of highway operating status is obtained based on mobile phone signaling data in real time |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106781479A CN106781479A (en) | 2017-05-31 |
| CN106781479B true CN106781479B (en) | 2019-03-22 |
Family
ID=58919182
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611205475.XA Active CN106781479B (en) | 2016-12-23 | 2016-12-23 | A method of highway operating status is obtained based on mobile phone signaling data in real time |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106781479B (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107463940B (en) * | 2017-06-29 | 2020-02-21 | 清华大学 | Vehicle type identification method and device based on mobile phone data |
| CN107644390B (en) * | 2017-09-26 | 2021-08-13 | 深圳大学 | Method and device for obtaining passenger flow data of track platform |
| CN107993452B (en) * | 2017-12-20 | 2021-06-08 | 杭州远眺科技有限公司 | Speed measurement method for detecting road passing speed on highway based on WIFI probe |
| CN108171992B (en) * | 2017-12-28 | 2020-12-01 | 重庆邮电大学 | Parallel highway vehicle speed calculation method based on mobile phone signaling big data |
| CN109996299B (en) * | 2017-12-30 | 2021-08-24 | 中国移动通信集团河北有限公司 | High-speed rail user identification method, device, equipment and medium |
| CN108322891B (en) * | 2018-01-11 | 2020-08-28 | 江苏欣网视讯软件技术有限公司 | Traffic area congestion identification method based on user mobile phone signaling |
| CN108492564B (en) * | 2018-04-18 | 2020-08-07 | 山东省交通规划设计院 | Method and system for measuring running speed of highway vehicle based on road network matching |
| CN109064745B (en) * | 2018-08-17 | 2020-10-27 | 中通服咨询设计研究院有限公司 | Grid-based highway vehicle screening method |
| CN109243175A (en) * | 2018-09-13 | 2019-01-18 | 广州丰石科技有限公司 | High speed algorithm based on mobile big data |
| CN110958558A (en) * | 2018-09-26 | 2020-04-03 | 北京融信数联科技有限公司 | Mobile big data-based mobile phone user space-time trajectory depicting method |
| CN111666360B (en) * | 2019-03-08 | 2023-08-01 | 北京四维图新科技股份有限公司 | Road construction state confirmation method and device |
| CN110708664B (en) * | 2019-10-11 | 2020-11-06 | 同帅科技（天津）有限公司 | Traffic flow sensing method and device, computer storage medium and electronic equipment |
| CN114333323B (en) * | 2022-01-05 | 2023-04-28 | 北京航空航天大学合肥创新研究院(北京航空航天大学合肥研究生院) | Expressway travel speed prediction method based on pressure characteristics |
| CN114881543B (en) * | 2022-07-06 | 2022-10-28 | 北京百度网讯科技有限公司 | Road grade determination method and device, electronic equipment and medium |
| CN117596551B (en) * | 2024-01-19 | 2024-04-09 | 浙江大学建筑设计研究院有限公司 | Green road network user behavior restoration method and device based on mobile phone signaling data |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008104673A1 (en) * | 2007-02-09 | 2008-09-04 | France Telecom | Method for estimating and signalling the density of mobile nodes in a road network |
| CN101510357A (en) * | 2009-03-26 | 2009-08-19 | 美慧信息科技（上海）有限公司 | Method for detecting traffic state based on mobile phone signal data |
| CN102332210A (en) * | 2011-08-04 | 2012-01-25 | 东南大学 | Method for extracting real-time urban road traffic flow data based on mobile phone positioning data |
| CN102708689A (en) * | 2012-06-19 | 2012-10-03 | 张家港市鸿嘉数字科技有限公司 | Real-time traffic monitoring system |
| CN104376718A (en) * | 2012-06-19 | 2015-02-25 | 张家港市鸿嘉数字科技有限公司 | Remote intelligent monitoring method for real-time traffic status |
| CN104778836A (en) * | 2015-04-03 | 2015-07-15 | 江苏省交通规划设计院股份有限公司 | Highway traffic state recognition method based on cellular signaling data quality perception |
| CN104978858A (en) * | 2015-06-15 | 2015-10-14 | 北京拓明科技有限公司 | Highway real-time traffic monitoring method based on mobile phone signaling |
| CN105070057A (en) * | 2015-07-24 | 2015-11-18 | 江苏省公用信息有限公司 | Real-time road traffic monitoring method |
| US9208682B2 (en) * | 2014-03-13 | 2015-12-08 | Here Global B.V. | Lane level congestion splitting |
| CN105303831A (en) * | 2015-10-20 | 2016-02-03 | 四川公用信息产业有限责任公司 | Method for determining congestion state of highway based on communication data |
| CN105608505A (en) * | 2015-12-22 | 2016-05-25 | 重庆邮电大学 | Cellular signaling data based track traffic travel mode identification method for resident |
| US9430942B2 (en) * | 2013-09-26 | 2016-08-30 | International Business Machines Corporation | Method and system for optimizing road traffic control in the presence of incidents |
| US9467806B2 (en) * | 2012-12-14 | 2016-10-11 | Omarco Network Solutions Limited | Improvements related to tracking moving objects |
-
2016
- 2016-12-23 CN CN201611205475.XA patent/CN106781479B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008104673A1 (en) * | 2007-02-09 | 2008-09-04 | France Telecom | Method for estimating and signalling the density of mobile nodes in a road network |
| CN101510357A (en) * | 2009-03-26 | 2009-08-19 | 美慧信息科技（上海）有限公司 | Method for detecting traffic state based on mobile phone signal data |
| CN102332210A (en) * | 2011-08-04 | 2012-01-25 | 东南大学 | Method for extracting real-time urban road traffic flow data based on mobile phone positioning data |
| CN102708689A (en) * | 2012-06-19 | 2012-10-03 | 张家港市鸿嘉数字科技有限公司 | Real-time traffic monitoring system |
| CN104376718A (en) * | 2012-06-19 | 2015-02-25 | 张家港市鸿嘉数字科技有限公司 | Remote intelligent monitoring method for real-time traffic status |
| US9467806B2 (en) * | 2012-12-14 | 2016-10-11 | Omarco Network Solutions Limited | Improvements related to tracking moving objects |
| US9430942B2 (en) * | 2013-09-26 | 2016-08-30 | International Business Machines Corporation | Method and system for optimizing road traffic control in the presence of incidents |
| US9208682B2 (en) * | 2014-03-13 | 2015-12-08 | Here Global B.V. | Lane level congestion splitting |
| CN104778836A (en) * | 2015-04-03 | 2015-07-15 | 江苏省交通规划设计院股份有限公司 | Highway traffic state recognition method based on cellular signaling data quality perception |
| CN104978858A (en) * | 2015-06-15 | 2015-10-14 | 北京拓明科技有限公司 | Highway real-time traffic monitoring method based on mobile phone signaling |
| CN105070057A (en) * | 2015-07-24 | 2015-11-18 | 江苏省公用信息有限公司 | Real-time road traffic monitoring method |
| CN105303831A (en) * | 2015-10-20 | 2016-02-03 | 四川公用信息产业有限责任公司 | Method for determining congestion state of highway based on communication data |
| CN105608505A (en) * | 2015-12-22 | 2016-05-25 | 重庆邮电大学 | Cellular signaling data based track traffic travel mode identification method for resident |
Non-Patent Citations (2)
| Title |
|---|
| 基于手机信令数据分析的路网运行监控系统;唐嘉立;《中国交通信息化》;20151231;正文第115-117页 |
| 基于手机信令的高速公路拥堵信息自动识别;沈潋，姜新农，周一;《中国交通信息化》;20151231;正文第119-122页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106781479A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106781479B (en) | A method of highway operating status is obtained based on mobile phone signaling data in real time | |
| CN110298500B (en) | Urban traffic track data set generation method based on taxi data and urban road network | |
| CN109410586A (en) | A kind of Traffic State Detection Method based on multivariate data fusion | |
| CN106197458A (en) | A kind of cellphone subscriber's trip mode recognition methods based on mobile phone signaling data and navigation route data | |
| WO2019047905A1 (en) | Road traffic analysis system, method and apparatus | |
| CN102521973B (en) | A kind of mobile phone switches the road matching method of location | |
| CN107766808A (en) | The method and system that Vehicle Object motion track clusters in road network space | |
| US20170091350A1 (en) | Near real-time modeling of pollution dispersion | |
| CN104851286B (en) | A kind of method of utilization public transit vehicle gps data dynamic evaluation road traffic state | |
| CN109561391B (en) | Highway service area pedestrian flow analysis method based on cellular network and Wi-Fi data | |
| CN105806351B (en) | A kind of road information estimation method and device | |
| CN106530716A (en) | Method for calculating highway section average speed based on mobile phone signaling data | |
| CN105355049B (en) | A kind of highway evaluation of running status method based on macroscopical parent map | |
| CN105788263B (en) | A kind of method by cellphone information predicted link congestion | |
| CN108961758B (en) | Road junction widening lane detection method based on gradient lifting decision tree | |
| CN104217593B (en) | A kind of method for obtaining road condition information in real time towards mobile phone travelling speed | |
| CN107784084A (en) | Road network generation method and system based on positioning data of vehicles | |
| CN106710208A (en) | Traffic state acquisition method and device | |
| CN107195180A (en) | A kind of traffic trip track extraction method and device based on the alert data of electricity | |
| CN106448173B (en) | A kind of long range trip traffic modal splitting method based on data in mobile phone | |
| CN108171993A (en) | A kind of vehicle on highway speed calculation method based on mobile phone signaling big data | |
| CN109859495A (en) | A method of overall travel speed is obtained based on RFID data | |
| CN108307316B (en) | User's high-precision method of discrimination is fast moved based on S1-MME data | |
| CN106651027A (en) | Internet regular bus route optimization method based on social network | |
| CN102779410A (en) | Parallel implementation method of multi-source heterogeneous traffic data fusion |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |