CN108920481A - A kind of road network method for reconstructing and system based on mobile phone location data - Google Patents
A kind of road network method for reconstructing and system based on mobile phone location data Download PDFInfo
- Publication number
- CN108920481A CN108920481A CN201810361223.9A CN201810361223A CN108920481A CN 108920481 A CN108920481 A CN 108920481A CN 201810361223 A CN201810361223 A CN 201810361223A CN 108920481 A CN108920481 A CN 108920481A
- Authority
- CN
- China
- Prior art keywords
- road
- road network
- mobile phone
- point
- triangle
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012549 training Methods 0.000 claims abstract description 47
- 235000007164 Oryza sativa Nutrition 0.000 claims abstract description 13
- 235000009566 rice Nutrition 0.000 claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000010801 machine learning Methods 0.000 claims abstract description 8
- 238000010276 construction Methods 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims description 17
- 238000010586 diagram Methods 0.000 claims description 14
- 241000209094 Oryza Species 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 6
- 238000007637 random forest analysis Methods 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000004220 aggregation Methods 0.000 claims description 2
- 230000002776 aggregation Effects 0.000 claims description 2
- 230000006399 behavior Effects 0.000 claims description 2
- 238000012706 support-vector machine Methods 0.000 claims description 2
- 230000002123 temporal effect Effects 0.000 claims description 2
- 238000005520 cutting process Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 3
- 240000007594 Oryza sativa Species 0.000 abstract 1
- 238000000605 extraction Methods 0.000 description 4
- 239000012141 concentrate Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/029—Location-based management or tracking services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Educational Administration (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Traffic Control Systems (AREA)
Abstract
A kind of road network method for reconstructing and system based on mobile phone location data, including：Raw data set is obtained, raw data set includes training dataset and predictive data set；Raw data set is to obtain tracing point according to mobile phone location data；To training dataset and predictive data set be utilized respectively two dimensions of time and density carry out feature constructions with combine, and generation training set and forecast set therewith；Using the training set of generation, sampled point disaggregated model is obtained based on the algorithm in machine learning；The disaggregated model learnt is applied to forecast set, tracing point is divided into waypoint and non-rice habitats point；The building of trade road is clicked through to the road identified based on constraint Delaunay triangulation network；According to topological relation, road axis is modified, to complete the reconstruction of road network.The invention proposes the new paragons of fining the road network production and update of a kind of high quality, will provide technical support for seamless navigation technology especially fine-grade road network aspect, and expand the coverage area of path navigation service.
Description
Technical field
The present invention relates to space-time trajectory data digging technology, especially a kind of road network based on mobile phone location data is rebuild
Method and system belong to GIS-Geographic Information System and intelligent transportation research field.
Background technique
Road is the important infrastructure in city, during rapid urban, how to be replenished in time and more new town
Road data is a problem urgently to be resolved to keep its accuracy, timeliness and degree of perfection.Residential communities are advised in city
Mould constantly expands, and internal passageway net is increasingly fine and complicated, but existing cell-level road information fine degree is not enough even
Missing causes outdoor thin portion navigation in residential block to there is " blind area ", causes very big obstacle to the development of outdoor seamless navigation technology.
At the same time, road network structure information is the core content of electronic map basic data, and even more electronic map realizes navigation feature, intelligence
The data basis of energy traffic system, therefore, accurately and in time more new town fine-grade road network is an important problem (city
City's fining road network refers to the cell-level branch that Intra-cell is respectively closed in city).
The convenience obtained with crowd-sourced track data and its real-time updated are greatly improved, academia and work
Industry can be readily available a large amount of track datas comprising space time information, so that the road network based on space-time trajectory data mentions
Method is taken to be possibly realized.Currently, there are mainly two types of Automatic Road Extraction methods：One is mention from vehicle GPS track data
Take road information；Another kind is to extract road network information based on high-resolution remote sensing image.But existing road network extracting method master
It concentrates in the research of each major trunk roads in city, and the fine road network of cell-level is usually ignored in city.Wherein, remote sensing image because
Remote for image-forming range, coverage area is big, while portraying the scarce capacity of different spaces scale fine information, based on such data
Processing method is only capable of preferably solving the information extraction of general major trunk roads, can not finally solve fining road network building and ask
Topic；In addition, the acquisition of GPS track data is limited inside the building range such as general school, cell to a certain extent, and pedestrian
Mobile trajectory data randomness is strong, and sample frequency is low, comprising various types of ground objects such as road, building and squares, so that being based on
Mobile trajectory data extracts road network, and there are numerous difficulties.
Summary of the invention
The technical problems to be solved by the invention are for road fine degree in existing cell-level road network not enough and part
The problem of missing, provides a kind of fining Road network extraction method, can be convenient the building work for accurately completing fine-grade road network
Make.This method will be converted into two classes of waypoint Yu non-rice habitats point the problem of carrying out road network reconstruction based on extensive pedestrian's data
Classification problem carries out the building of feature set using the various combination of two dimensions of time and density first, generates training set, then
Learn sampled point disaggregated model using random forests algorithm, application class model divides data set later, obtains road
Point set finally borrows constraint Delaunay triangulation network to having identified that road point set carries out road building and mention with road axis
It takes, to complete the building for refining road network in closed area.It includes the following steps：
Step 1) obtains raw data set, and raw data set includes training dataset and predictive data set；Wherein original number
It is that tracing point is obtained according to mobile phone location data according to collection；
Step 2) is utilized respectively two dimensions of time and density to training dataset and predictive data set and carries out feature construction
With combine, and generate training set and forecast set therewith；
Step 3) learns to obtain sampled point disaggregated model based on the algorithm in machine learning using the training set generated；
The disaggregated model learnt is applied to forecast set by step 4), and tracing point is divided into waypoint and non-rice habitats point；
Step 5) clicks through the building of trade road to the road identified based on constraint Delaunay triangulation network；
Step 6) modifies road axis, to complete the reconstruction of road network according to topological relation.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, training in the step 1
Data set is obtained according to following manner：The distribution pattern for obtaining multiple closing Intra-cell tracing points, selects tracing point distribution
The different multiple nonoverlapping small-scale regions of mode form training dataset.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, density in the step 2
Characteristic value includes：Binding time dimension carries out Delaunay Triangulation, the institute that each tracing point is connected to training dataset
There is the average length of adjacent triangle edges.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, density in the step 2
Characteristic value includes：Binding time dimension constructs Voronoi diagram, Voronoi diagram corresponding to each tracing point to training dataset
The inverse of middle area of a polygon.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, the time in the step 2
Specific step is as follows with the feature construction of two latitudes of density and combination：
By analyzing the behavior pattern of people, based on the space-time Density Distribution difference structure for closing small internal road and non-rice habitats
Feature space is built, that is, passes through track points amount, the adjoining triangle edges length based on constraint Delaunay triangulation network in buffer area
It is poor with the aggregation of three kinds of geometrical characteristic joint time dimension expression tracing point spatial and temporal distributions of area of polygon in Voronoi diagram
It is different；The identification of road waypoint is carried out based on the feature set building classifier that above-mentioned three kinds different geometry dimension control conditions are established.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, learn in the step 3
Algorithm includes random forest, support vector machines.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, for training dataset
In track point data generated based on road net data incomplete in OSM according to certain radius using corresponding OSM data
Line buffer area falls in the tracing point in buffer area labeled as 1, falls in the tracing point except buffer area labeled as 0, is based on this, will
The tracing point that training data is concentrated is divided into two class of waypoint and non-rice habitats point, label when as subsequent training machine learning model
Data.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, the step 5 is specifically wrapped
It includes：
The road point set building constraint Delaunay triangulation network formed using the road point extracted in step 4, root
According to adjacent triangle number, the triangle in the triangulation network is divided into 3 classes, respectively correspond the entrance of road, the trunk of road and
Intersection is respectively handled it, road axis is extracted, the specific implementation process is as follows described：
Building constraint Delaunay triangulation network, the triangle in the triangulation network is divided into the road point set identified
3 classes：Only having adjacent triangle on one side is A class, and it is B class that, which there is adjacent triangle on both sides, and there is adjacent triangle on three sides
For C class, A class triangle connects the vertex of the midpoint on unique proximity side corresponding thereto, and B class triangle connects in two neighbouring sides
Point, C class triangle connect the midpoint of center of gravity and three sides, in search connection procedure, first search for, terminate since C class triangle
In A class triangle or C class triangle, network a line is obtained, until having handled C class triangle；It uses and is opened from A class triangle again
Begin to search for, terminates at A class triangle, obtain all sides, to form road axis network.
Further, in the road network method for reconstructing of the invention based on mobile phone location data, the step 6 is specifically wrapped
It includes：
Optimization and perfect road center line, carry out the road axis near intersection processing is straightened, to complete to seal
Closed region refines the reconstruction of road network, and specific processing mode is as described below：
For the arc road network near intersection, since crosspoint, straight line is cropped with certain length, then compare sanction
The angle of rear remaining line is cut, when the angle of two lines meets a certain range, two lines of connection are into a line, when being crossroad,
Node of the intersection point as net is acquired, road network is ultimately formed, to complete the reconstruction of closed area road network.
The present invention is to solve its technical problem, additionally provides a kind of road network reconstructing system based on mobile phone location data,
Road network reconstruction is carried out using the road network method for reconstructing based on mobile phone location data of any of the above-described.
Implement the road network method for reconstructing and system of the invention based on mobile phone location data, the present invention is according to the row of people
Mining analysis is carried out to track data according to space-time density feature for mode, using included in a large amount of sparse mobile phone location data
In semantic information, propose it is a kind of constructed using mobile phone location data fining road network new method, biography can be improved
Unite road network update efficiency and reduce maintenance cost, substantially reduce road network extraction time, can be make up cell-level road
The not perfect and missing problem of information provides a kind of new mode.
Detailed description of the invention
Present invention will be further explained below with reference to the attached drawings and examples, in attached drawing：
Fig. 1 is flow chart of the method for the present invention；
Fig. 2 is the schematic diagram that the present invention calculates each area of a polygon in Voronoi diagram；
Fig. 3 is the schematic diagram for the average length that the present invention calculates the adjacent triangle edges of tracing point in Delaunay triangulation network.
Specific embodiment
For a clearer understanding of the technical characteristics, objects and effects of the present invention, now control attached drawing is described in detail
A specific embodiment of the invention.
The present invention relates to a kind of road network method for reconstructing based on mobile phone location data, carries out at fining road network
Reason.The present invention is analyzed by the distribution pattern to multiple closing Intra-cell tracing points first, then selects different community
The small range regions of interior multiple representative distribution patterns constructs training dataset, and remainder incorporates prediction data into
Collection for the track point data that training data is concentrated, carries out manual sort, is divided into waypoint and non-rice habitats point two at the same time
Class, the label value as training data；Using two dimensions of time and density carry out feature constructions with combine, wherein time dimension
It can be by being divided into 7 different periods (6 for one day:00-8:00,08:00-12:00,12:00-14:00,14:00-17:
00,17:00-20:00,20:00-24:00,24:00-06:00) it realizes, in terms of density, can be constructed from three angles
Feature：Points amount in track in buffer area, the area of each polygon based on Voronoi diagram and based on constraint Delaunay tri-
The average length of the adjacent triangle edges of each tracing point is netted at angle, combines all of above feature and corresponding with each tracing point
Label value forms training set；Using the training set of generation, based on the random forests algorithm study sampled point classification in sorting algorithm
Model；Later, the disaggregated model acquired is applied in predictive data set, for the road waypoint in the point data of identification prediction track
With non-rice habitats point, road point set is obtained；Finally in road point set, using constraint Delaunay triangulation network, road is extracted
Center line completes the building of fining road network in closing cell.
In order to make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with attached drawing and specific implementation
Example, the present invention will be described in further detail.It should be appreciated that specific embodiment described herein is only to explain this hair
It is bright, it is not intended to limit the present invention.
As shown in Figure 1, Fig. 1 is the side provided by the invention for carrying out fining road network based on mobile phone location data and rebuilding
Method flow chart, the described method comprises the following steps：
Step 1, raw data set is obtained, raw data set includes training dataset and predictive data set；Wherein original number
It is that tracing point is obtained according to mobile phone location data according to collection.Preferably, concentrate selection tracing point distribution pattern different from initial data
Multiple nonoverlapping small-scale regions form training datasets, remaining track point data is divided into predictive data set；Meanwhile it is right
In the track point data that training data is concentrated, needs to carry out manual sort, be divided into two class of waypoint and non-rice habitats point, as subsequent
Label data when training machine learning model.
Embodiment is implemented as follows：
For raw data set, tracing point distribution situation therein is analyzed first, it is multiple from different area of space selections
Different distribution patterns, such as：Center Assembled distribution, ribbon distribution, random discrete distribution etc., by a variety of distribution patterns
Track point data is merged into as training dataset, remaining track point data is divided into predictive data set；Meanwhile for training data
Corresponding OSM data can be used in the track point data of concentration, based on incomplete road net data in OSM according to certain half
Diameter generates line buffer area, falls in the tracing point in buffer area labeled as 1, falls in the tracing point except buffer area labeled as 0, is based on
This, is divided into two class of waypoint and non-rice habitats point for the tracing point that training data is concentrated, when as subsequent training machine learning model
Label data.
Wherein, the radius of line buffer area can be determined according to the mean radius of cell-level branch in city, Intra-cell road
The radius on road is generally within 8 meters.
Step 2, two dimensions of time and density are utilized respectively to training dataset and predictive data set to complete feature set
Building with combine, and generate training set and forecast set therewith.Wherein time dimension can be by being divided into 7 differences for one day
Period can carry out construction feature from three angles in terms of density：Points amount in track in buffer area, based on Voronoi diagram
The area of each polygon and the average length that triangle edges are abutted based on each tracing point of constraint Delaunay triangulation network, group
All of above feature and label value corresponding with each tracing point are closed, training set is formed.
Specific building process is as follows：
According to the trip of people rule, on time dimension, 7 periods was divided by one day, are 06 respectively:00-08:
00,8:00-12:00,12:00-14:00,14:00-17:00,17:00-20:00,20:00-24:00,24:00-06:00；?
In terms of density, buffer area radius takes 3 different values according to road conditions respectively：8m, 12m, 16m, by 3 buffer area radiuses and 7
A period carries out feature two-by-two and combines, and constructs 21 features, for indicating some tracing point in special time period, radius is
Density value in 8/12/16 meter, quantitative formula are as follows：
Wherein, { 6,8,12 } Radius=, Point={ P1,P2,....,Pn, Time=(6,8), (8,12), (12,
14), (14,17), (17,20), (20,24), (0,6) },It indicates in period TaUnder, tracing point PiIt is R in radiusbModel
The number of interior included tracing point is enclosed, the combination of different time sections and different radii feature is as shown in the table：
Later according to 7 different periods, training dataset is divided into 7 parts, construct respectively corresponding Voronoi diagram and
Delaunay triangulation network is constrained, for the Voronoi diagram that building obtains, the area of each polygon in figure is calculated, calculates therewith
The inverse of area, if Fig. 2 is a part of a Voronoi diagram, the present invention needs to calculate polygon R corresponding to tracing point P
Area, and the inverse and time dimension of area are combined, as 7 features of expression track vertex neighborhood density, are calculated public
Formula is as follows：
In above-mentioned formula,Indicate point PiThe area of corresponding polygon；Then for constraint Delaunay triangle
Net calculates the average length of the adjacent triangle edges of each tracing point, if Fig. 3 is a part of a Delaunay triangulation network, this
Invention needs to calculate the average length of all adjacent edges of tracing point P, and by it in conjunction with time dimension, adjacent as expression tracing point
Other 7 features of domain density, calculation formula are as follows：
In above-mentioned formula, Neighbor_P indicates point PiThe set of all consecutive points, L indicate the distance between consecutive points, such as
A in Fig. 3, b, c, d, e, f, g, h；Finally, by above-described 35 features and corresponding label described in step 1
Value merges, and generates training set.
Step 3, using the training set of generation, based in Python machine learning packet scikit-learn
RandomForestClassifier module, on the basis of partial parameters in adjustment algorithm, sampling that learning data is concentrated
Point disaggregated model.
Step 4, the sampled point disaggregated model learnt based on training set is applied in predictive data set, it is pre- for identifying
The road waypoint and non-rice habitats point in the point data of track are surveyed, road point set therein is extracted.
Step 5, using the road point set building constraint Delaunay triangulation network extracted in step 4, according to adjoining
Triangle in the triangulation network is divided into 3 classes by triangle number, is respectively corresponded the entrance of road, the trunk of road and road and is handed over
Prong is respectively handled it, road axis is extracted, the specific implementation process is as follows described：
Building constraint Delaunay triangulation network, the triangle in the triangulation network is divided into the road point set identified
3 classes：Only having adjacent triangle on one side is A class, and it is B class that, which there is adjacent triangle on both sides, and there is adjacent triangle on three sides
For C class, A class triangle connects the vertex of the midpoint on unique proximity side corresponding thereto, and B class triangle connects in two neighbouring sides
Point, C class triangle connect the midpoint of center of gravity and three sides, in search connection procedure, first search for, terminate since C class triangle
In A class triangle or C class triangle, network a line is obtained, until having handled C class triangle；It uses and is opened from A class triangle again
Begin to search for, terminates at A class triangle, obtain all sides, to form road axis network.
Step 6, optimization and perfect road center line, carry out the road axis near intersection processing is straightened, thus
The reconstruction of closed area fining road network is completed, specific processing mode is as described below：
For the arc road network near intersection, since crosspoint, straight line is cropped with certain length, then compare sanction
The angle of rear remaining line is cut, when the angle of two lines meets a certain range, two lines of connection are into a line, when being crossroad,
Node of the intersection point as net is acquired, road network is ultimately formed, to complete the reconstruction of closed area fining road network.
The embodiment of the present invention is described with above attached drawing, but the invention is not limited to above-mentioned specific
Embodiment, the above mentioned embodiment is only schematical, rather than restrictive, those skilled in the art
Under the inspiration of the present invention, without breaking away from the scope protected by the purposes and claims of the present invention, it can also make very much
Form, all of these belong to the protection of the present invention.
Claims (10)
1. a kind of road network method for reconstructing based on mobile phone location data, which is characterized in that include the following steps：
Step 1) obtains raw data set, and raw data set includes training dataset and predictive data set；Wherein raw data set
To obtain tracing point according to mobile phone location data；
Step 2) is utilized respectively two dimensions of time and density to training dataset and predictive data set and carries out feature construction and group
It closes, and generates training set and forecast set therewith；
Step 3) learns to obtain sampled point disaggregated model based on the algorithm in machine learning using the training set generated；
The disaggregated model learnt is applied to forecast set by step 4), and tracing point is divided into waypoint and non-rice habitats point；
Step 5) clicks through the building of trade road to the road identified based on constraint Delaunay triangulation network；
Step 6) modifies road axis, to complete the reconstruction of road network according to topological relation.
2. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that
Training dataset is obtained according to following manner in the step 1：Obtain the distribution of multiple closing Intra-cell tracing points
Mode, the multiple nonoverlapping small-scale regions for selecting tracing point distribution pattern different form training dataset.
3. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that the step
Density feature value includes in 2：Binding time dimension carries out Delaunay Triangulation, company of each tracing point institute to training dataset
The average length of all of its neighbor triangle edges connect.
4. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that the step
Density feature value includes in 2：Binding time dimension constructs Voronoi diagram to training dataset, corresponding to each tracing point
The inverse of area of a polygon in Voronoi diagram.
5. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that the step
Specific step is as follows for the feature construction of two latitudes of time and density and combination in 2：
By analyzing the behavior pattern of people, constructed based on the space-time Density Distribution difference for closing small internal road and non-rice habitats special
Levy space, i.e., by track points amount in buffer area, the adjoining triangle edges length based on constraint Delaunay triangulation network and
The aggregation of three kinds of geometrical characteristic joint time dimension expression tracing point spatial and temporal distributions of the area of polygon is poor in Voronoi diagram
It is different；The identification of road waypoint is carried out based on the feature set building classifier that above-mentioned three kinds different geometry dimension control conditions are established.
6. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that the step
Learning algorithm includes random forest, support vector machines in 3.
7. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that for training
Track point data in data set, using corresponding OSM data, based on incomplete road net data in OSM according to certain half
Diameter generates line buffer area, falls in the tracing point in buffer area labeled as 1, falls in the tracing point except buffer area labeled as 0, is based on
This, is divided into two class of waypoint and non-rice habitats point for the tracing point that training data is concentrated, when as subsequent training machine learning model
Label data.
8. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that the step
5 specifically include：
The road point set building constraint Delaunay triangulation network formed using the road point extracted in step 4, according to neighbour
Triangle number is connect, the triangle in the triangulation network is divided into 3 classes, respectively corresponds the entrance of road, the trunk of road and road
Intersection is respectively handled it, road axis is extracted, the specific implementation process is as follows described：
Building constraint Delaunay triangulation network, is divided into 3 classes for the triangle in the triangulation network in the road point set identified：
Only having adjacent triangle on one side is A class, and it is B class that, which there is adjacent triangle on both sides, and it is C that, which there is adjacent triangle on three sides,
Class, A class triangle connect the vertex of the midpoint on unique proximity side corresponding thereto, and B class triangle connects the midpoint on two neighbouring sides,
C class triangle connects the midpoint of center of gravity and three sides, in search connection procedure, first searches for since C class triangle, terminates at A
Class triangle or C class triangle, obtain network a line, until having handled C class triangle；Again using since A class triangle
Search, terminates at A class triangle, obtains all sides, to form road axis network.
9. the road network method for reconstructing according to claim 1 based on mobile phone location data, which is characterized in that the step
6 specifically include：
Optimization and perfect road center line, carry out the road axis near intersection processing is straightened, to complete enclosed area
Domain refines the reconstruction of road network, and specific processing mode is as described below：
For the arc road network near intersection, since crosspoint, straight line is cropped with certain length, then after comparing cutting
The angle of remaining line, when the angle of two lines meets a certain range, two lines of connection are into a line, when being crossroad, acquire
Node of the intersection point as net, ultimately forms road network, to complete the reconstruction of closed area road network.
10. a kind of road network reconstructing system based on mobile phone location data, which is characterized in that using any one of claim 1-9
The road network method for reconstructing based on mobile phone location data carries out road network reconstruction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810361223.9A CN108920481B (en) | 2018-04-20 | 2018-04-20 | Road network reconstruction method and system based on mobile phone positioning data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810361223.9A CN108920481B (en) | 2018-04-20 | 2018-04-20 | Road network reconstruction method and system based on mobile phone positioning data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108920481A true CN108920481A (en) | 2018-11-30 |
| CN108920481B CN108920481B (en) | 2020-11-27 |
Family
ID=64404350
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810361223.9A Expired - Fee Related CN108920481B (en) | 2018-04-20 | 2018-04-20 | Road network reconstruction method and system based on mobile phone positioning data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108920481B (en) |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109684384A (en) * | 2018-12-29 | 2019-04-26 | 太原科技大学 | A kind of track data space-time density Analysis System and its analysis method |
| CN109934865A (en) * | 2019-03-12 | 2019-06-25 | 重庆市勘测院 | Road network center line extraction method |
| CN110008872A (en) * | 2019-03-25 | 2019-07-12 | 浙江大学 | A kind of road network extracting method of combination track of vehicle and remote sensing images |
| CN110096535A (en) * | 2019-05-16 | 2019-08-06 | 北京航空航天大学 | A kind of traffic elasticity communication mode method for digging based on elastic triangle theory |
| CN110580388A (en) * | 2019-08-28 | 2019-12-17 | 北方工业大学 | Channel network extraction method based on crowd-sourced trajectory data |
| CN110633342A (en) * | 2019-07-29 | 2019-12-31 | 武汉光庭信息技术股份有限公司 | Lane topology network generation method |
| CN110688958A (en) * | 2019-09-28 | 2020-01-14 | 中国测绘科学研究院 | Intersection identification method based on GoogleLeNet neural network |
| CN110988880A (en) * | 2019-12-12 | 2020-04-10 | 南京莱斯电子设备有限公司 | Geographic information extraction and target tracking method based on SMR target track |
| CN111426329A (en) * | 2020-03-12 | 2020-07-17 | 深圳小库科技有限公司 | Road generation method and device, computer equipment and storage medium |
| CN111833630A (en) * | 2019-12-31 | 2020-10-27 | 北京嘀嘀无限科技发展有限公司 | Method and device for determining data release position, storage medium and electronic equipment |
| CN112013862A (en) * | 2020-07-31 | 2020-12-01 | 深圳大学 | Pedestrian network extraction and updating method based on crowdsourcing trajectory |
| CN112566025A (en) * | 2020-12-10 | 2021-03-26 | 南京市城市与交通规划设计研究院股份有限公司 | Bus passenger flow source destination identification method and device, electronic equipment and storage medium |
| CN112598724A (en) * | 2021-03-01 | 2021-04-02 | 武大吉奥信息技术有限公司 | Improved TIN-based planar element center line extraction method |
| CN113008251A (en) * | 2021-02-22 | 2021-06-22 | 湖南大学 | Digital map updating method for unstructured roads in closed area |
| CN113052084A (en) * | 2021-03-26 | 2021-06-29 | 中国地质大学（武汉） | Community-level vector road network extraction method based on mobile phone positioning data |
| WO2022031222A1 (en) * | 2020-08-01 | 2022-02-10 | Grabtaxi Holdings Pte. Ltd. | Processing appratus and method for generating route navigation data |
| CN117891892A (en) * | 2024-03-14 | 2024-04-16 | 南京师范大学 | Forest guard track optimization algorithm based on historical patrol data |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8396257B2 (en) * | 2008-07-25 | 2013-03-12 | Navteq B.V. | End user image open area maps |
| CN105069824A (en) * | 2015-08-11 | 2015-11-18 | 中南大学 | GPS data based automatic construction method and system for open-pit mine road network |
| CN106528740A (en) * | 2016-11-04 | 2017-03-22 | 中科宇图科技股份有限公司 | Delaunay triangular net-based road center line extraction method |
| CN106778605A (en) * | 2016-12-14 | 2017-05-31 | 武汉大学 | Remote sensing image road net extraction method under navigation data auxiliary |
-
2018
- 2018-04-20 CN CN201810361223.9A patent/CN108920481B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8396257B2 (en) * | 2008-07-25 | 2013-03-12 | Navteq B.V. | End user image open area maps |
| CN105069824A (en) * | 2015-08-11 | 2015-11-18 | 中南大学 | GPS data based automatic construction method and system for open-pit mine road network |
| CN106528740A (en) * | 2016-11-04 | 2017-03-22 | 中科宇图科技股份有限公司 | Delaunay triangular net-based road center line extraction method |
| CN106778605A (en) * | 2016-12-14 | 2017-05-31 | 武汉大学 | Remote sensing image road net extraction method under navigation data auxiliary |
Non-Patent Citations (2)
| Title |
|---|
| 杨伟 等: "基于众源轨迹数据的道路中心线提取", 《地理与地理信息科学》 * |
| 杨伟 等: "运用约束Delaunay三角网从众源轨迹线提取道路边界", 《测绘学报》 * |
Cited By (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109684384A (en) * | 2018-12-29 | 2019-04-26 | 太原科技大学 | A kind of track data space-time density Analysis System and its analysis method |
| CN109934865B (en) * | 2019-03-12 | 2020-11-27 | 重庆市勘测院 | Road network central line extraction method |
| CN109934865A (en) * | 2019-03-12 | 2019-06-25 | 重庆市勘测院 | Road network center line extraction method |
| CN110008872A (en) * | 2019-03-25 | 2019-07-12 | 浙江大学 | A kind of road network extracting method of combination track of vehicle and remote sensing images |
| CN110008872B (en) * | 2019-03-25 | 2021-04-06 | 浙江大学 | Road network extraction method combining vehicle track and remote sensing image |
| CN110096535A (en) * | 2019-05-16 | 2019-08-06 | 北京航空航天大学 | A kind of traffic elasticity communication mode method for digging based on elastic triangle theory |
| CN110096535B (en) * | 2019-05-16 | 2021-06-29 | 北京航空航天大学 | Elastic traffic propagation mode mining method based on elastic triangle theory |
| CN110633342A (en) * | 2019-07-29 | 2019-12-31 | 武汉光庭信息技术股份有限公司 | Lane topology network generation method |
| CN110580388A (en) * | 2019-08-28 | 2019-12-17 | 北方工业大学 | Channel network extraction method based on crowd-sourced trajectory data |
| CN110580388B (en) * | 2019-08-28 | 2023-04-07 | 北方工业大学 | Channel network extraction method based on crowd-sourced trajectory data |
| CN110688958A (en) * | 2019-09-28 | 2020-01-14 | 中国测绘科学研究院 | Intersection identification method based on GoogleLeNet neural network |
| CN110988880A (en) * | 2019-12-12 | 2020-04-10 | 南京莱斯电子设备有限公司 | Geographic information extraction and target tracking method based on SMR target track |
| CN111833630A (en) * | 2019-12-31 | 2020-10-27 | 北京嘀嘀无限科技发展有限公司 | Method and device for determining data release position, storage medium and electronic equipment |
| CN111426329A (en) * | 2020-03-12 | 2020-07-17 | 深圳小库科技有限公司 | Road generation method and device, computer equipment and storage medium |
| CN111426329B (en) * | 2020-03-12 | 2021-10-22 | 深圳小库科技有限公司 | Road generation method and device, computer equipment and storage medium |
| CN112013862A (en) * | 2020-07-31 | 2020-12-01 | 深圳大学 | Pedestrian network extraction and updating method based on crowdsourcing trajectory |
| CN112013862B (en) * | 2020-07-31 | 2022-06-17 | 深圳大学 | Pedestrian network extraction and updating method based on crowdsourcing trajectory |
| US11869348B2 (en) | 2020-08-01 | 2024-01-09 | Grabtaxi Holdings Pte. Ltd. | Processing apparatus and method for generating route navigation data |
| WO2022031222A1 (en) * | 2020-08-01 | 2022-02-10 | Grabtaxi Holdings Pte. Ltd. | Processing appratus and method for generating route navigation data |
| CN112566025A (en) * | 2020-12-10 | 2021-03-26 | 南京市城市与交通规划设计研究院股份有限公司 | Bus passenger flow source destination identification method and device, electronic equipment and storage medium |
| CN112566025B (en) * | 2020-12-10 | 2021-10-01 | 南京市城市与交通规划设计研究院股份有限公司 | Bus passenger flow source destination identification method and device, electronic equipment and storage medium |
| CN113008251B (en) * | 2021-02-22 | 2022-11-25 | 湖南大学 | Digital map updating method for unstructured roads in closed area |
| CN113008251A (en) * | 2021-02-22 | 2021-06-22 | 湖南大学 | Digital map updating method for unstructured roads in closed area |
| CN112598724A (en) * | 2021-03-01 | 2021-04-02 | 武大吉奥信息技术有限公司 | Improved TIN-based planar element center line extraction method |
| CN113052084A (en) * | 2021-03-26 | 2021-06-29 | 中国地质大学（武汉） | Community-level vector road network extraction method based on mobile phone positioning data |
| CN117891892A (en) * | 2024-03-14 | 2024-04-16 | 南京师范大学 | Forest guard track optimization algorithm based on historical patrol data |
| CN117891892B (en) * | 2024-03-14 | 2024-05-31 | 南京师范大学 | Forest guard track optimization method based on historical patrol data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108920481B (en) | 2020-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108920481A (en) | A kind of road network method for reconstructing and system based on mobile phone location data | |
| CN112133090A (en) | Multi-mode traffic distribution model construction method based on mobile phone signaling data | |
| CN112561146A (en) | Large-scale real-time traffic flow prediction method based on fuzzy logic and depth LSTM | |
| CN109635748B (en) | Method for extracting road characteristics in high-resolution image | |
| Wang et al. | Quality analysis of open street map data | |
| CN106897417A (en) | A kind of construction method of the city space holographic map based on the fusion of multi-source big data | |
| CN109189949B (en) | A kind of population distribution calculation method | |
| CN102521973A (en) | Road matching method for mobile phone switching positioning | |
| WO2021013190A1 (en) | Meteorological parameter-based high-speed train positioning method and system in navigation blind zone | |
| CN108921173B (en) | A kind of deep learning method of combination OSM and remote sensing image extraction overpass | |
| Tang et al. | Exploring urban travel patterns using density-based clustering with multi-attributes from large-scaled vehicle trajectories | |
| CN115512216A (en) | City functional area fine recognition method coupling block space-time characteristics and ensemble learning | |
| CN109447160A (en) | A kind of method of image and vector road junction Auto-matching | |
| CN107656987A (en) | A kind of subway station function method for digging based on LDA models | |
| CN108090624B (en) | Urban ecological safety simulation and prediction method for improving cellular automaton | |
| CN111539432B (en) | Method for extracting urban road by using multi-source data to assist remote sensing image | |
| CN113932821B (en) | Track map matching method based on continuous window average direction characteristics | |
| CN109657614A (en) | Traffic accident situ of taking photo by plane reconnoitres middle road and automatic vehicle identification method | |
| CN106971535A (en) | A kind of urban traffic blocking index calculating platform based on Floating Car GPS real time datas | |
| CN105844031B (en) | A kind of urban transportation gallery recognition methods based on mobile phone location data | |
| CN115098798A (en) | Multi-element traffic flow urban space organization pattern recognition method fusing map data | |
| CN110210667A (en) | O-D spatial and temporal distributions prediction technique based on mutative scale Geographical Weighted Regression Model | |
| CN113903173A (en) | Vehicle track feature extraction method based on directed graph structure and LSTM | |
| CN109800903A (en) | A kind of profit route planning method based on taxi track data | |
| CN106778547A (en) | The construction method and device of remote sensing image terrain classification pattern library |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201127 |