CN116934561B - Quick traffic service population measuring and calculating method, equipment and medium - Google Patents
Quick traffic service population measuring and calculating method, equipment and medium Download PDFInfo
- Publication number
- CN116934561B CN116934561B CN202311197120.0A CN202311197120A CN116934561B CN 116934561 B CN116934561 B CN 116934561B CN 202311197120 A CN202311197120 A CN 202311197120A CN 116934561 B CN116934561 B CN 116934561B
- Authority
- CN
- China
- Prior art keywords
- population
- data
- rapid transit
- demographic
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000008878 coupling Effects 0.000 claims abstract description 5
- 238000010168 coupling process Methods 0.000 claims abstract description 5
- 238000005859 coupling reaction Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 30
- 238000012216 screening Methods 0.000 claims description 20
- 238000010586 diagram Methods 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 10
- 238000000691 measurement method Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 2
- 239000010410 layer Substances 0.000 description 30
- 238000004364 calculation method Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 3
- 229910000831 Steel Inorganic materials 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 229940054107 isochron Drugs 0.000 description 2
- MOYKHGMNXAOIAT-JGWLITMVSA-N isosorbide dinitrate Chemical compound [O-][N+](=O)O[C@H]1CO[C@@H]2[C@H](O[N+](=O)[O-])CO[C@@H]21 MOYKHGMNXAOIAT-JGWLITMVSA-N 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 239000002344 surface layer Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
- G08G1/0129—Traffic data processing for creating historical data or processing based on historical data
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a rapid traffic service population measuring and calculating method, equipment and medium, and belongs to the technical field of traffic service. The method comprises the steps of obtaining geographic longitude and latitude coordinates of different rapid transit stations; taking the geographic longitude and latitude coordinates of the different rapid transit stations as the center to obtain the corresponding isochrones of the different rapid transit stations; fusing the isochronous circles corresponding to the different rapid transit stations to obtain fused isochronous circles; acquiring population raster data; and superposing and coupling the fused isochrone and population raster data to obtain population data covered by the fused isochrone. By the method, the electronic equipment and the computer medium for running the method, the population served by the rapid transit can be accurately measured and calculated.
Description
Technical Field
The invention relates to the technical field of traffic service, in particular to a rapid traffic service population measuring and calculating method, device and medium.
Background
The traffic waiting time circle refers to a region which starts from a specific place and can be reached in a certain time by selecting a specified traffic mode. The traffic equal time circle is the geospatial expression of travel time consumption, and can be used for analyzing regional traffic analysis, urban job and occupancy separation, urban traffic convenience analysis, public service facility accessibility analysis and the like.
On the basis of the traffic equal time circle, the population number covered by the traffic equal time circle is calculated, and the method has important guiding significance for traffic planning and municipal planning. At present, no related algorithm is used for measuring and calculating the population number of the rapid transit service.
Therefore, how to measure the population data covered by the traffic circle, and to provide a rapid traffic service population measurement method are urgent to be solved by those skilled in the art.
Disclosure of Invention
In view of the above, the present invention provides a rapid traffic service population measuring and calculating method, apparatus and medium for solving at least some of the technical problems in the background art.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
the invention discloses a rapid traffic service population measuring and calculating method, which comprises the following steps:
obtaining geographic longitude and latitude coordinates of different rapid transit stations;
taking the geographic longitude and latitude coordinates of the different rapid transit stations as the center to obtain the corresponding isochrones of the different rapid transit stations;
fusing the isochronous circles corresponding to the different rapid transit stations to obtain fused isochronous circles;
acquiring population raster data;
and superposing and coupling the fused isochrone and population raster data to obtain population data covered by the fused isochrone.
Further, the method for acquiring the geographic longitude and latitude coordinates of different rapid transit stations specifically comprises the following steps:
acquiring POI (point of interest) point class data based on a digital line drawing map element vector data set, wherein the POI point class data comprises four attributes including a position name, a position address, a position coordinate and a position class;
and acquiring geographic longitude and latitude coordinates of different rapid transit stations based on the POI point data.
Further, the different rapid transit stops include civil airports, high-speed rail stations, and highway toll stations.
Further, the method for acquiring the isochronous circles corresponding to the different rapid transit stations specifically comprises the following steps:
taking geographic longitude and latitude coordinates of different rapid traffic stations as centers, and calculating an area reachable with the centers as starting points in a designated time;
and returning the reachable area as a polygonal outline which can be displayed on the map, wherein the polygonal outline is an isochronal circle corresponding to different rapid transit stations in the appointed time.
Further, the step of calculating the reachable area with the center as the starting point within the specified time specifically includes:
subtracting the set transit time from the designated time to serve as effective traffic time;
and calculating the reachable area in the effective traffic time as the practically reachable area taking the center as the starting point in the appointed time.
Further, fusing the equal time circles corresponding to the different rapid transit stations specifically comprises:
and obtaining vector image layers corresponding to the equal time circles of different rapid transit stations, and merging and fusing the vector image layers to obtain a fused equal time circle vector image.
Further, the demographic raster data is obtained by:
based on internet handset positioning information, obtaining a population thermodynamic diagram, and calculating the population quantity per 5 square kilometers on the population thermodynamic diagram as population raster data in combination with the population total number.
Further, the fused isochrone and population raster data are overlapped and coupled to obtain population data covered by the fused isochrone, which specifically comprises the following steps:
constructing an external rectangle of each polygonal contour based on the polygonal contour of the fused equal-time circle;
carrying out coarse screening statistics on population data by a longitude and latitude numerical comparison method to obtain a population statistics data point set meeting the condition of the circumscribed rectangular range;
using a space vector connection algorithm to carry out fine screening statistics on population data, and acquiring all population data point sets contained in the polygonal outline from all population data point sets meeting the condition of the circumscribed rectangular range;
traversing each demographic data point binary group in the demographic data fine screening condition point set, inquiring a demographic information table corresponding to demographic raster data, returning demographic data values corresponding to the binary groups in the demographic information table, and calculating demographic aggregate values covered by all polygonal outlines in the polygonal outline set, namely the fused demographic data covered by the isochrone.
In another aspect, the invention discloses an electronic device comprising a memory and a processor, the memory storing a computer program, the processor implementing the rapid transit service population measurement method according to any of the preceding claims when executing the computer program.
Furthermore, the invention also discloses a computer readable storage medium, on which a computer program is stored, which when being executed by a processor, implements the rapid traffic service population measuring method according to any one of the invention.
Compared with the prior art, the invention discloses a rapid traffic service population measuring and calculating method, device and medium, which have the following beneficial effects:
the invention provides a large-scale traffic equal-time population measuring and calculating method, which can be used for overlapping traffic equal-time-period population data to obtain the number of fast traffic service population covered by the traffic equal-time-period.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only embodiments of the present invention, and that other drawings can be obtained according to the provided drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic overall flow chart of a rapid transit service population measurement method according to an embodiment of the present invention.
Fig. 2 is a schematic diagram of demographic raster data according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of a geographic range of a rapid transit service according to an embodiment of the present invention.
Fig. 4 is a schematic diagram of rapid transit population density provided by an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Example 1
The embodiment of the invention discloses a rapid traffic service population measuring and calculating method, which is shown in figure 1 and comprises the following steps:
obtaining geographic longitude and latitude coordinates of different rapid transit stations;
taking the obtained geographic longitude and latitude coordinates of different rapid transit stations as the center, and obtaining the corresponding isochrones of the different rapid transit stations;
fusing the isochronous circles corresponding to the different rapid transit stations to obtain fused isochronous circles;
acquiring population raster data;
and superposing and coupling the fused isochrone and population raster data to obtain population data covered by the fused isochrone.
Specifically, the geographical longitude and latitude coordinates of different rapid transit stations can be obtained by the following steps:
firstly, acquiring POI (point of interest) class data based on a digital line drawing map element vector data set, wherein the POI class data generally comprises four attributes of a position name, a position address, a position coordinate and a position class; and then, obtaining the geographic longitude and latitude coordinates of different rapid transit stations by using POI point data.
In the embodiment of the invention, different rapid transit stations refer to rapid transit stations of different categories, including civil airports, high-speed rail stations and expressway toll stations.
In the embodiment of the invention, the method for acquiring the corresponding isochronous circles of the different rapid transit stations specifically comprises the following steps:
firstly, taking geographical longitude and latitude coordinates of different rapid transit stations as centers, calculating an area which can be reached in a specified time by taking the centers as starting points, such as an area which can be reached in one hour.
And then returning the reachable area in the appointed time as a polygonal outline which can be displayed on the map, wherein the obtained polygonal outline is an isochronal circle corresponding to different rapid traffic stations in the appointed time.
As a preferred embodiment, in the step of calculating the reachable area with the geographical longitude and latitude coordinates of the rapid transit station as the starting point within the specified time, an appropriate amount of transit time needs to be uniformly considered, for example, when measuring and calculating the reachable area within one hour, each six minutes of time for reaching the rapid transit station and leaving the rapid transit station can be considered to be subtracted as effective traffic time; that is, the actual reachable area within 48 minutes is calculated as the actual reachable area of the hour circle by subtracting the set transit time of 12 minutes from one hour and taking the actual 48 minutes as the effective traffic time of the hour circle of one hour.
The invention needs to fuse the equal time circles corresponding to different rapid traffic stations, and specifically comprises the following steps:
and obtaining vector image layers corresponding to the equal time circles of different rapid transit stations, and merging and fusing the vector image layers to obtain a fused equal time circle vector image.
In the embodiment of the invention, population raster data is utilized to represent population density space distribution, and the population raster data is obtained through the following steps:
based on the internet handset positioning information, a population thermodynamic diagram is obtained, and the population quantity per 5 square kilometers on the population thermodynamic diagram is calculated as population raster data (i.e. population density space distribution data) in combination with the population total number.
In the embodiment of the invention, finally, overlapping and coupling are needed to be carried out on the fused isochrone and population raster data so as to obtain the population data covered by the fused isochrone, which comprises the following steps:
1) Constructing an external rectangle of each polygonal contour based on the polygonal contour of the fused equal-time circle;
2) The population data is roughly screened and counted by a longitude and latitude numerical comparison method, and the specific process is as follows:
recall all demographic data point sets S meeting the circumscribed rectangular range condition, the demographic data point sets S being represented by the following expression:
S=S 1 ∪S 2 ∪…∪S i ∪…S |S| ，
S i ={(x 1 ,y 1 ),(x 2 ,y 2 ),…,(x j ,y j ),…,(x |Si| ,y |Si| ) (x) wherein j ,y j ) Coarse screening of statistical Point set S to satisfy population data i Demographic data point (longitude, latitude) doublet, |s i I is a point set S meeting population data coarse screening statistical conditions i The number of elements in S is a population statistics data point set meeting the condition of the circumscribed rectangle range (namely a population data coarse screening statistics condition point set S i And/s| represents the number of elements in all demographic data point sets S that meet the circumscribed rectangular range condition.
Wherein, the population data coarse screening statistical point set S is satisfied i Mid-demographic data point (longitude, latitude) doublet (x j ,y j ) The method meets the following conditions:wherein X is min The minimum longitude of the circumscribing rectangle for the polygonal outline; x is X max The maximum longitude of the circumscribing rectangle for the polygonal outline; y is Y min The minimum latitude of the rectangle is circumscribed by the polygonal outline; y is Y max The maximum latitude of the rectangle is circumscribed by the polygonal outline.
3) And (3) carrying out fine screening statistics on population data by using a space vector connection algorithm, wherein the specific process is as follows:
recall all demographic data point sets contained within the polygonal outline in all demographic data point sets S satisfying the circumscribed rectangular range conditionI.e.
；
Satisfy->Wherein, the method comprises the steps of, wherein,representing space vector join algorithm,/->Screening the statistical condition point set for satisfying the population data>Demographic data point (longitude, latitude) doublet of (i.e.)>Screening the statistical condition point set for satisfying the population data>The number of the elements in the steel sheet is equal to the number of the elements in the steel sheet,for all demographic data point sets contained inside the polygonal outline (satisfy demographic data undersize statistical condition point set +.>Union) of->Representation ofAll population data point sets contained inside the polygonal outline +.>The number of elements in the matrix.
4) Traversing each set of points meeting demographic data undersize statistical conditionsDemographic data point doublet in (a)Querying a demographic information table corresponding to the demographic raster data, and returning a demographic data value corresponding to the binary group in the demographic information table>And calculating a demographic summary value Z of polygon coverage in the polygon profile set, namely fused demographic data of the isochronal circle coverage. The specific expression is as follows: />。
In the embodiment of the invention, the polygon contour set is represented by a symbol P, and p= { P 1 ,p 2 ,…,p i ,…,p |p| And |p| is the number of elements in the polygon vector set P, i.e. the number of equal time circles (the number of polygon vectors at airports, toll stations, high-speed rail stations).
Because the vector space operation recalled by the points in the polygon has high complexity and calculated amount, the vector space operation process is split into two steps of coarse screening and fine screening, wherein the calculation of the coarse screening step adopts non-space operation, namely a numerical calculation mode, and the calculation complexity, the memory use amount and the calculated amount of about 99.9 percent can be greatly reduced; the fine screening step adopts vector space connection operation on the basis of the residual 0.1% calculated amount, so that the comprehensive complexity and calculated amount of the algorithm are in a moderate and controllable range, the operation time of the algorithm is short, the occupied memory is less, and a common personal computer can meet the operation requirement of an algorithm program and efficiently process large-scale data.
According to the basic nature of the set, there are no repeating elements within the set, and thus the setIn the construction process, repeated (longitude and latitude) binary groups generated in recall of the polygons with space coincidence can be removed, so that convenience is brought to calculation of the demographic summary value Z covered by the polygons in the polygon vector set P, and calculation accuracy is ensured.
At a specific operation level, the above process can also be expressed as: and taking the population raster data as a bottom layer, taking the fused isochron vector diagram as an upper layer, programming by using python language, overlapping the layers by an independently developed overlapping algorithm, and acquiring population data covered by the fused isochron according to the overlapped layers.
The calculation process of the superposition algorithm is described as follows:
reading population raster data and fused isochronous circle vector data;
defining the same space coordinate system for the population raster data and the fused isochrone vector data;
traversing the fused equal-time circles, and screening population grid points meeting the set longitude and latitude conditions for each equal-time circle;
extracting population grid points which meet the set longitude and latitude conditions, and detecting repeatability;
recording longitude and latitude and population data corresponding to the population grid points through repeated detection;
and accumulating the population data corresponding to all the extracted grid points to obtain the population data covered by the fused isochrone.
Example 2
The embodiment of the invention provides an electronic device, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor realizes the rapid traffic service population measuring and calculating method according to any one of the embodiment 1 of the invention when executing the computer program.
The processor referred to in this invention may be a central processing unit (Central Processing Unit, CPU), or other general purpose processor, digital signal processor (Digital Signal Processor, DSP), application specific integrated circuit (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory referred to in the present invention may mainly include a memory program area and a memory data area, wherein the memory program area may store an operating system, an application program required for at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, phonebook, etc.) created according to the use of the handset, etc. In addition, the memory may include high-speed random access memory, and may also include non-volatile memory, such as a hard disk, memory, plug-in hard disk, smart Media Card (SMC), secure Digital (SD) Card, flash Card (Flash Card), at least one disk storage device, flash memory device, or other volatile solid-state storage device.
Example 3
The invention also discloses a computer readable storage medium, on which a computer program is stored, which when executed by a processor implements the rapid traffic service population measuring method according to any one of the embodiments 1 of the invention.
The computer readable storage medium of the present invention may be any form of storage medium readable by a processor of a computer apparatus, including but not limited to, non-volatile memory, ferroelectric memory, etc., on which a computer program is stored.
The inventive principles of this invention are further explained below in connection with a 1 hour traffic hour fast population measuring method.
First, the rapid traffic connotation is decomposed, and the rapid traffic is decomposed into civil aviation service, high-speed rail service and expressway service. The population from the usual place, which can obtain the service within 1 hour using highway traffic, is defined as the population which enjoys the rapid traffic service within 1 hour, namely the population which can reach civil airports within 1 hour, can reach high-speed rail sites within 1 hour and can reach highway toll booths within 1 hour.
And secondly, obtaining the geographic longitude and latitude of the three types of rapid transit stations. The method comprises the steps that all the civil aviation airport, the high-speed railway station and the expressway toll gate adopt part of own data to check the POI data from the internet online map.
Wherein, POI (Point of Interest, commonly called point of interest) refers broadly to point class data in internet electronic map, generally comprising four attributes of name, address, coordinates, and category, originating from point class map element vector data set in basic mapping effort DLG (Digital Line Graphic, digital line drawing map) product. In the embodiment of the invention, POI information of national civil aviation airport, high-speed rail station/railway station and expressway toll station is obtained.
When the high-speed railway stations are acquired, the high-speed railway stations are matched through a high-speed railway line as a buffer zone, and the common-speed railway stations and the high-speed railway stations are distinguished through whether the high-speed railway stations are in the high-speed railway line buffer zone or not.
An isochrone centered around the three types of stations is then obtained.
And D, taking all the longitude and latitude coordinates of the stations obtained in the step two as centers, and solving an outward-expansion isochrone. Here, considering the traffic modes of aviation travel, high-speed rail travel and highway travel, which are not door-to-door, to the tail end, when the measurement and calculation are carried out for 1 hour, a proper amount of transit time is uniformly considered for 12 minutes, namely, transit time of 6 minutes is adopted at each of the two ends, and the time of 48 minutes, which can be reached, is actually obtained.
The working principle of acquiring a 48-minute equivalent circle centered on each station by using a mapbox online map API is as follows: an area reachable from a certain position within a specified time is calculated, and the reachable area is returned as a polygonal contour that can be displayed on a map. The polygonal outline is calculated using a grid, the outline returning format being GeoJSON Feature Collection. The usage restrictions of the Mapbox equal time circle API include: up to 300 requests per minute; only 1 coordinate is supported per request; up to 4 isochrones are supported per request; the longest time that can be specified for the isochrone profile is 60 minutes; the maximum distance that can be specified for the isochrone contour is 100000 meters.
Vector processing is then performed on the hours and administrative divisions.
And (3) loading the population raster data and the airport equal time circle, the high-speed railway equal time circle and the expressway toll gate equal time circle acquired from the step three by adopting arc GIS software. Meanwhile, in order to calculate provincial population data, a national administrative division shp file is loaded.
And combining three layers of the time circle of the high-speed rail service, the time circle of the airport service and the time circle of the expressway service by applying the function of combining vector layers, selecting a default value by a CRS parameter, and automatically adding the combined and output layers into the existing layers.
Because overlapping of different types of isochrones during demographics can result in repeated computation of underlying grid demographics, operating on overlapping isochrones merges into a complete surface layer. And opening [ fusion ] in the arc GIS software toolbox, inputting a layer, selecting a layer after merging, setting fusion field parameters to be full selection, and automatically adding the layer after merging into the existing layer. The fused isochrones do not have mutual overlapping parts, so repeated calculation of the data of the mouth in subsequent operation is avoided.
All administrative areas in the country are divided into a plurality of surface layers according to provincial administrative area boundaries. Opening in the toolbox (crossing), setting an input layer as a fused layer, setting an overlapped layer as a nationwide administrative division layer, keeping other parameter settings default, automatically adding the intersected layer in the existing layer, and renaming an equal-time circle coverage map.
And finally, acquiring the population number under the coverage range of the equal time circle.
Loading population distribution data of China. The population density space distribution data of the national academy of sciences is adopted, the data is a population thermodynamic diagram estimated based on the position information of Internet personnel, the population quantity per 5 square kilometers estimated by the national population quantity is combined, the space resolution is 5km, the data unit is ten thousands of people per 5 square kilometers, and the data is in a grid format. As shown in fig. 2.
The function of 'partition statistics' in the toolbox is selected, an input layer is set to be equal-circle coverage, a grid layer is selected to be population grid data, statistical information to be calculated is set to be sum, and the layer is automatically added in the existing layer after partition statistics. Right clicking on the layer icon selects the export-vector layer save as the format parameter selection csv file. After the population total number is exported, the population total number covered by the time circles of different provincial administrative areas and the like can be obtained.
The invention can be presented in two modes, namely, only the geographical range which can be served by the rapid transit is presented, as shown in fig. 3; in the second mode, population raster data is used as a bottom layer in arcgis, corresponding traffic equal-time circles to be presented are used as an upper layer, a graph is rendered, a final population density schematic diagram is shown in fig. 4, at this time, pure white parts in the graph are uncovered parts of the equal-time circles, and the color part is light and dark, so that population density is represented.
Through the method, the method for measuring and calculating the population ratio of the rapid transit service in 1 hour can be established for the national space range, and the method can measure and calculate the population ratio of the rapid transit service in 1 hour of the designated geographic range or region in the country and correspondingly perform visual presentation.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (8)
1. A rapid transit service population measuring and calculating method, comprising the steps of:
obtaining geographic longitude and latitude coordinates of different rapid transit stations;
taking the geographic longitude and latitude coordinates of the different rapid transit stations as the center to obtain the corresponding isochrones of the different rapid transit stations;
fusing the isochronous circles corresponding to the different rapid transit stations to obtain fused isochronous circles; fusing the equal time circles corresponding to the different rapid transit stations, which comprises the following steps: the method comprises the steps of obtaining vector image layers corresponding to equal time circles of different rapid transit stations, merging and fusing the vector image layers to obtain fused equal time circle vector images;
acquiring population raster data;
superposing and coupling the fused isochrone and population raster data to obtain population data covered by the fused isochrone; the method specifically comprises the following steps:
1) Constructing an external rectangle of each polygonal contour based on the polygonal contour of the fused equal-time circle;
2) The population data is roughly screened and counted by a longitude and latitude numerical comparison method, and the specific process is as follows:
recall all demographic data point sets S meeting the circumscribed rectangular range condition, the demographic data point sets S being represented by the following expression:
S＝S 1 ∪S 2 ∪…∪S i ∪…S |S| ，
S i ＝{(x 1 ,y 1 ),(x 2 ,y 2 ),…,(x j ,y j ),…,(x |Si| ,y |Si| ) (x) wherein j ,y j ) Coarse screening of statistical Point set S to satisfy population data i In demographic data point longitude and latitude binary group, |S i I is a point set S meeting population data coarse screening statistical conditions i S is a population statistics data point set which meets the condition of the circumscribed rectangle scope, and S is the number of elements in the population statistics data point set S which meets the condition of the circumscribed rectangle scope;
wherein, satisfy population data coarse screeningStatistics Point set S i Mid-demographic data point longitude and latitude doublets (x j ,y j ) The method meets the following conditions:
wherein X is min The minimum longitude of the circumscribing rectangle for the polygonal outline; x is X max The maximum longitude of the circumscribing rectangle for the polygonal outline; y is Y min The minimum latitude of the rectangle is circumscribed by the polygonal outline; y is Y max The maximum latitude of the rectangle is circumscribed by the polygonal outline;
3) And (3) carrying out fine screening statistics on population data by using a space vector connection algorithm, wherein the specific process is as follows:
recall all demographic data point sets contained within the polygonal outline in all demographic data point sets S satisfying the circumscribed rectangular range conditionI.e.
Satisfy->Wherein Sjoin (·, ·) represents the space vector join algorithm, ++>Screening the statistical condition point set for satisfying the population data>In demographic data point longitude and latitude binary group,/->Screening the statistical condition point set for satisfying the population data>The number of elements in->For all demographic data points contained inside the polygonal outline, +.>Representing the set of all demographic data points contained inside the polygonal outline +.>The number of elements in the matrix;
4) Traversing each set of points meeting demographic data undersize statistical conditionsDemographic data point doublet in (a)Querying a demographic information table corresponding to the demographic raster data, and returning a demographic data value corresponding to the binary group in the demographic information table +.>And calculating a demographic summary value Z of polygon coverage in the polygon profile set, namely fused population data of the isochronal circle coverage, wherein the specific expression is as follows: />
2. The rapid transit service population measurement method of claim 1, wherein obtaining geographic longitude and latitude coordinates of different rapid transit stations comprises:
acquiring POI (point of interest) point class data based on a digital line drawing map element vector data set, wherein the POI point class data comprises four attributes including a position name, a position address, a position coordinate and a position class;
and acquiring geographic longitude and latitude coordinates of different rapid transit stations based on the POI point data.
3. The rapid transit service population measurement method of claim 1, wherein the different rapid transit stops include civil airports, high-speed rail stops, and highway tollgates.
4. The rapid transit service population measurement method according to claim 1, wherein the obtaining of the isochronous circles corresponding to the different rapid transit stations comprises:
taking geographic longitude and latitude coordinates of different rapid traffic stations as centers, and calculating an area reachable with the centers as starting points in a designated time;
and returning the reachable area as a polygonal outline which can be displayed on the map, wherein the polygonal outline is an isochronal circle corresponding to different rapid transit stations in the appointed time.
5. The rapid transit service population measurement method according to claim 4, wherein the step of calculating the region reachable within the specified time with the center as the starting point comprises:
subtracting the set transit time from the designated time to serve as effective traffic time;
and calculating the reachable area in the effective traffic time as the practically reachable area taking the center as the starting point in the appointed time.
6. The rapid transit service population measurement method of claim 1, wherein the population raster data is obtained by:
based on internet handset positioning information, obtaining a population thermodynamic diagram, and calculating the population quantity per 5 square kilometers on the population thermodynamic diagram as population raster data in combination with the population total number.
7. An electronic device comprising a memory and a processor, the memory storing a computer program, the processor implementing the rapid transit service population measurement method of any one of claims 1-6 when executing the computer program.
8. A computer readable storage medium having stored thereon a computer program, which when executed by a processor implements the rapid transit service population measuring method of any one of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311197120.0A CN116934561B (en) | 2023-09-18 | 2023-09-18 | Quick traffic service population measuring and calculating method, equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311197120.0A CN116934561B (en) | 2023-09-18 | 2023-09-18 | Quick traffic service population measuring and calculating method, equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116934561A CN116934561A (en) | 2023-10-24 |
| CN116934561B true CN116934561B (en) | 2023-12-15 |
Family
ID=88382925
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311197120.0A Active CN116934561B (en) | 2023-09-18 | 2023-09-18 | Quick traffic service population measuring and calculating method, equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116934561B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110298778A (en) * | 2019-06-10 | 2019-10-01 | 东南大学 | One kind being based on raster data and the modified traffic zone demographics method of area ratio |
| CN114898013A (en) * | 2022-07-15 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | Traffic isochronal generation method, electronic device and storage medium |
| CN115630988A (en) * | 2022-12-22 | 2023-01-20 | 北京大学深圳研究生院 | Land road comprehensive traffic accessibility measuring and calculating method and device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2070006B1 (en) * | 2006-09-08 | 2015-07-01 | FortiusOne, Inc. | System and method for web enabled geo-analytics and image processing |
-
2023
- 2023-09-18 CN CN202311197120.0A patent/CN116934561B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110298778A (en) * | 2019-06-10 | 2019-10-01 | 东南大学 | One kind being based on raster data and the modified traffic zone demographics method of area ratio |
| CN114898013A (en) * | 2022-07-15 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | Traffic isochronal generation method, electronic device and storage medium |
| CN115630988A (en) * | 2022-12-22 | 2023-01-20 | 北京大学深圳研究生院 | Land road comprehensive traffic accessibility measuring and calculating method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116934561A (en) | 2023-10-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhu et al. | Analysis of the impact of the MRT system on accessibility in Singapore using an integrated GIS tool | |
| US9230167B2 (en) | Methods, apparatus, and articles of manufacture to measure geographical features using an image of a geographical location | |
| CN108665149B (en) | Method and system for measuring and calculating degree of mixing based on urban functions | |
| CN107679189A (en) | A kind of point of interest update method, device, server and medium | |
| CN110413855B (en) | Region entrance and exit dynamic extraction method based on taxi boarding point | |
| CN108053240B (en) | Method and device for generating vehicle-mounted advertisement delivery bus route scheme | |
| CN114548811B (en) | Airport reachability detection method and device, electronic equipment and storage medium | |
| CN112129310A (en) | Map display method and device | |
| CN112288311A (en) | Convenient and fast residential area supporting facility metering method based on POI data | |
| CN115640639A (en) | Urban updating scheme simulation and cost estimation method based on digital urban design | |
| JP2003208465A (en) | Method of providing required time information, execution system therefor, and processing program therefor | |
| Shi | Research on accessibility and equity of urban transport based on multisource big data | |
| CN110555432A (en) | Method, device, equipment and medium for processing interest points | |
| CN110362640A (en) | A kind of method for allocating tasks and device based on electronic map data | |
| CN115100012A (en) | Method for calculating walking accessibility of rail transit station | |
| CN116934561B (en) | Quick traffic service population measuring and calculating method, equipment and medium | |
| CN108038728B (en) | Bus body advertisement delivery line scheme generation method and device and electronic equipment | |
| CN108122083A (en) | Cargo distribution scheduling method and device | |
| CN117114210A (en) | Barrier-free public facility layout optimization method, device, equipment and storage medium | |
| CN117196197A (en) | Public transportation site layout optimization method | |
| WO2008072429A1 (en) | System for displaying image data associated with map information | |
| CN110132233A (en) | Current relief map drawing practice under a kind of CASS environment based on point cloud data | |
| CN113487341B (en) | Urban business strategy data processing method | |
| CN115311887A (en) | Method, equipment and storage medium for constructing high-quality public service resource accessibility model | |
| CN114677854A (en) | Public transport accessibility evaluation method and device |
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 |