CN105117816A - City impedance calculation method based on points of interest - Google Patents
City impedance calculation method based on points of interest Download PDFInfo
- Publication number
- CN105117816A CN105117816A CN201510430023.0A CN201510430023A CN105117816A CN 105117816 A CN105117816 A CN 105117816A CN 201510430023 A CN201510430023 A CN 201510430023A CN 105117816 A CN105117816 A CN 105117816A
- Authority
- CN
- China
- Prior art keywords
- interest
- plot
- center
- individual
- point
- 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
Landscapes
- Traffic Control Systems (AREA)
Abstract
The invention relates to a city impedance calculation method based on points of interest, comprising the following steps: (1) acquiring data of points of interest in a city from an Internet medium, including the name, type, latitude and longitude positions and operation time data of each point of interest; (2) sorting the operation time data of the points of interest, and outputting a standard format; (3) using city road network data to mark out plot units and determine the center of each plot with a polygon completely closed by a road system in the city as the basis, and mapping the points of interest to the corresponding plot units based on the latitude and longitude positions; (4) planning the routes from the center of each plot to all the points of interest within the reachable range, and obtaining route parameters; and (5) calculating the impedance value corresponding to each plot unit in the city according to the route planning result by considering the route parameters and population parameters in the corresponding plot unit. The method is beneficial to quantifying the impedance values of different plots of a city in different periods of time, thus providing reference for traffic management and urban land use planning.
Description
Technical field
The present invention relates to planning technology field, smart city, particularly a kind of city impedance computation method based on point of interest, for urban traffic control and plot planning provide reference.
Background technology
For the calculating of Urban Traffic impedance, recent domestic scholar has carried out corresponding research for this, and proposes some computing method.Wherein based on the utility models of Discrete Choice Model, the utility value that one is parameter with population or density is arranged to each terminal in region, the point that individual cognition selects effectiveness maximum is as choice for traveling, the size of impedance is differentiated with this, but have ignored complicated urban road situation, make the accuracy of method not high.And traditional BPR speed flow models of being obtained by regretional analysis and improving one's methods, although add the link travel time function of the factors such as distance, running time, expense, but all need a large amount of data as support in realization, comprise the demographic data etc. in the geographic distribution positional information of urban public utilities, the road Traffic Volume of different time sections, research range, this is very unobtainable in the many countries comprising China.
And along with the arrival of large data age, millions of user provides a large amount of points of interest and (is called for short POI, PointofInterest) data and data of registering, these data are that smart city service item provides data source, equally also for the research of city impedance computation method is provided convenience.Utilize the data of these magnanimity, the data input of alternative classic method, and overcome the inadequate problem of real-time, set up the impedance computation method of the uniqueness of a set of data Network Based.
Summary of the invention
The object of the present invention is to provide a kind of city impedance computation method based on point of interest, the method is conducive to the resistance value quantizing different plot, city different time sections, thus provides reference for traffic flow management and urban land_use plan.
For achieving the above object, technical scheme of the present invention is: a kind of city impedance computation method based on point of interest, comprises the following steps:
Steps A: gather incity, city interest point data from internet medium, comprising the title of point of interest, type, longitude and latitude position, business hours data;
Step B: data preparation is carried out to the business hours of point of interest diversified in specifications, outputting standard form;
Step C: utilize city road network data, based on the polygon that roadnet encloses completely in city, marks off Land unit, determines the center, plot of each Land unit; Described point of interest is mapped to correspondingly in module unit by longitude and latitude position one by one;
Step D: path planning is carried out to all points of interest in center to coverage pattern, each plot, and way to acquire parameter, comprising road network distance, number of turns, angle of turn, crossing parameter through path;
Step e: according to route programming result, considers the path parameter of all points of interest in center, plot to coverage pattern and the demographic parameter correspondingly in module unit, calculates the resistance value that in city, each Land unit is corresponding.
Further, in step e, the computing method of the resistance value that each Land unit is corresponding are as follows:
Step e 1: the road network Distance geometry calculating all points of interest in center, plot to coverage pattern
d i , computing formula is as follows:
Wherein
d i represent the
ithe road network distance sum of center, plot all points of interest in coverage pattern of individual Land unit,
d i,j be
ithe center, plot of individual Land unit is in coverage pattern
jindividual point of interest road network distance through path,
den i for
ithe sum of all points of interest in center, the plot coverage pattern of individual Land unit;
Step e 2: the angle of turn parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
a i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest angle of turn parameter through path,
lfor on path
lsecondary turning,
θ l be
lthe angle of turn of secondary turning,
d l be
lsecondary turning and
lroad section length between turning for-1 time,
for regulatory factor;
Step e 3: the number of turns parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
t i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest number of turns parameter through path;
Step e 4: the crossing parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
i i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest crossing parameter through path,
pfor on path
pindividual crossing,
turn p be
pthe mode of advancing of individual crossing, comprises right-hand rotation
turnRight, keep straight on
straight, turn left
turnLefthree kinds of modes, corresponding turning cost coefficient is respectively 1,2,3,
d p be
pindividual crossing and
ppath between-1 crossing;
Step e 5: calculate the demographic parameter of Land unit in setting-up time section, computing formula is as follows:
Wherein
p i,t represent the
iindividual Land unit exists
tdemographic parameter in time period,
kfor in center, the plot coverage pattern of this Land unit and
tin time period business the
kindividual point of interest,
sign k,t represent the
kindividual point of interest exists
tsum of registering in time period,
s i be
ithe area of individual Land unit;
Step e 6: according to above-mentioned parameter, calculate the resistance value of Land unit, computing formula is as follows:
Wherein
m i,t represent the
iindividual Land unit exists
tresistance value in time period,
den i for
ithe sum of all points of interest in center, the plot coverage pattern of individual Land unit,
den i,t for
iin center, the plot coverage pattern of individual Land unit and
tthe sum of all points of interest of business in time period.
The invention has the beneficial effects as follows by by network point of interest and register data and city road network data, route transportation planning method combines, substitute road route information and the population trip data parameter of Conventional impedance method, construct a city impedance computation method based on social networks point of interest, the method is applicable to the plot impedance computation of all cities different time sections, quantize the trip cost of different plot, city different time, for citizens' activities provide guidance, the real-time traffic administration in different plot is carried out and city Global motion planning provides reference for Urban Planner, there is very strong practicality and wide application prospect.
Accompanying drawing explanation
Fig. 1 is the realization flow figure of the embodiment of the present invention.
Embodiment
Below in conjunction with drawings and the specific embodiments, the present invention is described in further detail.
The present invention is by introducing the city interest point data of internet mass and data of registering, and combine with city road network data, the plot divided with city is for base unit, planning plot is to the path of point of interest, and obtained the parameters such as road network distance on outbound path, number of turns, angle of turn, crossing and population by rational computing formula, thus calculate city impedance.
The present invention is a kind of city impedance computation method based on point of interest, as shown in Figure 1, comprises the following steps:
Steps A: comment on from masses the internet media such as net and gather incity, city interest point data, comprising the data such as title, type, longitude and latitude position, business hours of point of interest.
Step B: carry out data preparation to the business hours of point of interest diversified in specifications, exports shape as the standard format of " 8:00-20:00 ".
Step C: utilize city road network data, based on the polygon that roadnet encloses completely in city, marks off Land unit, determines the center, plot of each Land unit.Mapped to one by one correspondingly in module unit by longitude and latitude position by described point of interest, impedance computation will be discussed in units of Land unit.
Step D: path planning is carried out to all points of interest in center to coverage pattern, each plot, and way to acquire parameter, comprising the parameter such as road network distance, number of turns, angle of turn, crossing through path.Coverage pattern refers to center, plot for the center of circle, and setting Radius parameter, the region in radius is the coverage pattern at this center, plot, and the point of interest in radius is the point of interest in this center, plot coverage pattern.
Step e: according to route programming result, considers the path parameter of all points of interest in center, plot to coverage pattern and the demographic parameter correspondingly in module unit, calculates the resistance value that in city, each Land unit is corresponding.Concrete, the computing method of the resistance value that each Land unit is corresponding are as follows:
Step e 1: the road network Distance geometry calculating all points of interest in center, plot to coverage pattern
d i , computing formula is as follows:
Wherein
d i represent the
ithe road network distance sum of center, plot all points of interest in coverage pattern of individual Land unit,
d i,j be
ithe center, plot of individual Land unit is in coverage pattern
jindividual point of interest road network distance through path,
den i for
ithe sum of all points of interest in center, the plot coverage pattern of individual Land unit.
Step e 2: the angle of turn parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
a i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest angle of turn parameter through path,
lfor on path
lsecondary turning,
θ l be
lthe angle of turn of secondary turning,
d l be
lsecondary turning and
lroad section length between turning for-1 time,
for regulatory factor.
Step e 3: the number of turns parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
t i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest number of turns parameter through path.
Step e 4: the crossing parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
i i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest crossing parameter through path,
pfor on path
pindividual crossing,
turn p be
pthe mode of advancing of individual crossing, comprises right-hand rotation
turnRight, keep straight on
straight, turn left
turnLefthree kinds of modes, corresponding turning cost coefficient is respectively 1,2,3,
d p be
pindividual crossing and
ppath between-1 crossing.
Step e 5: calculate the demographic parameter of Land unit in setting-up time section, computing formula is as follows:
Wherein
p i,t represent the
iindividual Land unit exists
tdemographic parameter in time period,
kfor in center, the plot coverage pattern of this Land unit and
tin time period business the
kindividual point of interest,
sign k,t represent the
kindividual point of interest exists
tsum of registering in time period,
s i be
ithe area of individual Land unit.
Step e 6: according to above-mentioned parameter, calculate the resistance value of Land unit, computing formula is as follows:
Wherein
m i,t represent the
iindividual Land unit exists
tresistance value in time period,
den i for
ithe sum of all points of interest in center, the plot coverage pattern of individual Land unit,
den i,t for
iin center, the plot coverage pattern of individual Land unit and
tthe sum of all points of interest of business in time period.
Be more than preferred embodiment of the present invention, all changes done according to technical solution of the present invention, when the function produced does not exceed the scope of technical solution of the present invention, all belong to protection scope of the present invention.
Claims (2)
1., based on a city impedance computation method for point of interest, it is characterized in that, comprise the following steps:
Steps A: gather incity, city interest point data from internet medium, comprising the title of point of interest, type, longitude and latitude position, business hours data;
Step B: data preparation is carried out to the business hours of point of interest diversified in specifications, outputting standard form;
Step C: utilize city road network data, based on the polygon that roadnet encloses completely in city, marks off Land unit, determines the center, plot of each Land unit; Described point of interest is mapped to correspondingly in module unit by longitude and latitude position one by one;
Step D: path planning is carried out to all points of interest in center to coverage pattern, each plot, and way to acquire parameter, comprising road network distance, number of turns, angle of turn, crossing parameter through path; Coverage pattern refers to center, plot for the center of circle, setting Radius parameter, and the region in radius is the coverage pattern at this center, plot, and the point of interest in radius is the point of interest in this center, plot coverage pattern;
Step e: according to route programming result, considers the path parameter of all points of interest in center, plot to coverage pattern and the demographic parameter correspondingly in module unit, calculates the resistance value that in city, each Land unit is corresponding.
2. a kind of city impedance computation method based on point of interest according to claim 1, it is characterized in that, in step e, the computing method of the resistance value that each Land unit is corresponding are as follows:
Step e 1: the road network Distance geometry calculating all points of interest in center, plot to coverage pattern
d i , computing formula is as follows:
Wherein
d i represent the
ithe road network distance sum of center, plot all points of interest in coverage pattern of individual Land unit,
d i,j be
ithe center, plot of individual Land unit is in coverage pattern
jindividual point of interest road network distance through path,
den i for
ithe sum of all points of interest in center, the plot coverage pattern of individual Land unit;
Step e 2: the angle of turn parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
a i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest angle of turn parameter through path,
lfor on path
lsecondary turning,
θ l be
lthe angle of turn of secondary turning,
d l be
lsecondary turning and
lroad section length between turning for-1 time,
for regulatory factor;
Step e 3: the number of turns parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
t i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest number of turns parameter through path;
Step e 4: the crossing parameter calculating each point of interest in center, plot to coverage pattern, computing formula is as follows:
，
Wherein
i i,j represent the
iin center to coverage pattern, individual plot
jindividual point of interest crossing parameter through path,
pfor on path
pindividual crossing,
turn p be
pthe mode of advancing of individual crossing, comprises right-hand rotation
turnRight, keep straight on
straight, turn left
turnLefthree kinds of modes, corresponding turning cost coefficient is respectively 1,2,3,
d p be
pindividual crossing and
ppath between-1 crossing;
Step e 5: calculate the demographic parameter of Land unit in setting-up time section, computing formula is as follows:
Wherein
p i,t represent the
iindividual Land unit exists
tdemographic parameter in time period,
kfor in center, the plot coverage pattern of this Land unit and
tin time period business the
kindividual point of interest,
sign k,t represent the
kindividual point of interest exists
tsum of registering in time period,
s i be
ithe area of individual Land unit;
Step e 6: according to above-mentioned parameter, calculate the resistance value of Land unit, computing formula is as follows:
Wherein
m i,t represent the
iindividual Land unit exists
tresistance value in time period,
den i for
ithe sum of all points of interest in center, the plot coverage pattern of individual Land unit,
den i,t for
iin center, the plot coverage pattern of individual Land unit and the sum of all points of interest of doing business within the t time period.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510430023.0A CN105117816B (en) | 2015-07-22 | 2015-07-22 | A kind of city impedance computation method based on point of interest |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510430023.0A CN105117816B (en) | 2015-07-22 | 2015-07-22 | A kind of city impedance computation method based on point of interest |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105117816A true CN105117816A (en) | 2015-12-02 |
| CN105117816B CN105117816B (en) | 2018-11-27 |
Family
ID=54665798
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510430023.0A Expired - Fee Related CN105117816B (en) | 2015-07-22 | 2015-07-22 | A kind of city impedance computation method based on point of interest |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105117816B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106844642A (en) * | 2017-01-20 | 2017-06-13 | 杭州电子科技大学 | A kind of method that the density of population in road network grid is calculated based on GIS |
| CN107346313A (en) * | 2016-05-05 | 2017-11-14 | 清华大学 | The method and apparatus that virtual face is excavated |
| CN108287856A (en) * | 2017-01-10 | 2018-07-17 | 北京三快在线科技有限公司 | Service range determines method, apparatus and electronic equipment |
| CN109308549A (en) * | 2017-11-20 | 2019-02-05 | 章文洁 | A kind of City complex activity analysis method and system |
| CN113935525A (en) * | 2021-10-11 | 2022-01-14 | 上海交通大学 | Dynamic traffic flow distribution method based on POI influence |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103531024A (en) * | 2013-10-28 | 2014-01-22 | 武汉旭云科技有限公司 | Dynamic traffic network urban road feature model and modeling method thereof |
-
2015
- 2015-07-22 CN CN201510430023.0A patent/CN105117816B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103531024A (en) * | 2013-10-28 | 2014-01-22 | 武汉旭云科技有限公司 | Dynamic traffic network urban road feature model and modeling method thereof |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107346313A (en) * | 2016-05-05 | 2017-11-14 | 清华大学 | The method and apparatus that virtual face is excavated |
| CN107346313B (en) * | 2016-05-05 | 2020-11-27 | 清华大学 | Method and device for mining virtual surface |
| CN108287856A (en) * | 2017-01-10 | 2018-07-17 | 北京三快在线科技有限公司 | Service range determines method, apparatus and electronic equipment |
| CN108287856B (en) * | 2017-01-10 | 2021-05-25 | 北京三快在线科技有限公司 | Service range determining method and device and electronic equipment |
| CN106844642A (en) * | 2017-01-20 | 2017-06-13 | 杭州电子科技大学 | A kind of method that the density of population in road network grid is calculated based on GIS |
| CN106844642B (en) * | 2017-01-20 | 2020-02-18 | 杭州电子科技大学 | Method for calculating population density in road network grid based on GIS |
| CN109308549A (en) * | 2017-11-20 | 2019-02-05 | 章文洁 | A kind of City complex activity analysis method and system |
| CN113935525A (en) * | 2021-10-11 | 2022-01-14 | 上海交通大学 | Dynamic traffic flow distribution method based on POI influence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105117816B (en) | 2018-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cooper et al. | Using multiple hybrid spatial design network analysis to predict longitudinal effect of a major city centre redevelopment on pedestrian flows | |
| Zhai et al. | Using mobile signaling data to exam urban park service radius in Shanghai: methods and limitations | |
| CN105117816A (en) | City impedance calculation method based on points of interest | |
| Wang et al. | Estimating dynamic origin-destination data and travel demand using cell phone network data | |
| CN108564226A (en) | A kind of public bus network optimization method based on taxi GPS and mobile phone signaling data | |
| CN109688532A (en) | A kind of method and device dividing city function region | |
| CN102700576B (en) | Urban Rail Transit passenger flow monitoring method | |
| Li et al. | Gravity model in dockless bike-sharing systems within cities | |
| HK1067844A1 (en) | System and method for providing traffic information using operational data of a wireless network | |
| CN103063222B (en) | With the navigation system and its operating method for turning to restriction scheme | |
| CN104303013B (en) | Location-based numerical data method and product | |
| CN109520499B (en) | Method for realizing regional real-time isochrones based on vehicle GPS track data | |
| CN103823951B (en) | Method for quantifying characteristics of road network and effect of characteristics of road network on land utilization | |
| CN110414795B (en) | Newly-increased high-speed rail junction accessibility influence method based on improved two-step mobile search method | |
| Huang | Data integration for urban transport planning | |
| Zhang et al. | Trip route planning for bicycle-sharing systems | |
| Xiao et al. | Does the culture service supply of green spaces match the demand of residents in a new district? A perspective from China. | |
| CN108710996A (en) | Gather region hotel addressing appraisal procedure in hotel based on tourism trip time and space usage | |
| Xie et al. | The accessibility of Nanjing urban park based on GIS | |
| Kahraman et al. | In the effects of accessibility factors on land values in the CBD of Izmir | |
| Li | Application of Spatial‐Temporal Behavioral Trajectory Analysis in the Space Design of Digital Villages | |
| Dunn et al. | Evaluating walkability in the age of open data: OpenStreetMap and community-level transportation analysis | |
| Haijuan et al. | Evaluation of accessibility to urban public sports facilities: A GIS approach based on network analysis model | |
| Kovács et al. | Measuring visitor flows using mobile positioning data in three Hungarian second-tier cities | |
| CN104063830B (en) | A kind of exchange method of transport information |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| 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: 20181127 Termination date: 20210722 |