CN106909692B - Method for calculating urban public facility coverage radiation index - Google Patents

Method for calculating urban public facility coverage radiation index Download PDF

Info

Publication number
CN106909692B
CN106909692B CN201710138930.7A CN201710138930A CN106909692B CN 106909692 B CN106909692 B CN 106909692B CN 201710138930 A CN201710138930 A CN 201710138930A CN 106909692 B CN106909692 B CN 106909692B
Authority
CN
China
Prior art keywords
service
attribute
calculating
resident
field
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
Application number
CN201710138930.7A
Other languages
Chinese (zh)
Other versions
CN106909692A (en
Inventor
杨伯钢
张伟松
刘博文
高骞
王淼
王凯松
汤大顺
任海英
刘清丽
关丽
陈娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Surveying and Mapping
Original Assignee
Beijing Institute of Surveying and Mapping
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Surveying and Mapping filed Critical Beijing Institute of Surveying and Mapping
Priority to CN201710138930.7A priority Critical patent/CN106909692B/en
Publication of CN106909692A publication Critical patent/CN106909692A/en
Application granted granted Critical
Publication of CN106909692B publication Critical patent/CN106909692B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services

Abstract

The invention discloses a method for calculating the coverage radiation index of urban public facilities, which comprises the following steps: step S1, acquiring city basic data and importing the city basic data into a file geographic database; step S2, sparse population space point data is extracted; step S3, constructing a network data set; step S4, creating an OD cost matrix; step S5, splitting the field and the associated attribute; step S6, calculating the ratio of the service supply quantity to the resident demand quantity of the service point according to the service position; step S7, calculating a ratio of service supply quantity to resident demand quantity of the resident according to the location of the resident; and step S8, summarizing the urban public facility coverage radiation indexes and drawing a public facility coverage radiation index map. The invention reflects the balance degree of the urban public facilities by calculating the covering radiation indexes of the urban public facilities such as schools, hospitals, cultural facilities and the like in a certain radius range, and embodies the service capability of the urban public facilities in a certain space range.

Description

Method for calculating urban public facility coverage radiation index
Technical Field
The invention relates to a method for planning urban public facilities, in particular to a method for calculating a coverage radiation index of the urban public facilities, and belongs to the technical field of urban planning and geographic information systems.
Background
The urban public facilities are the supporting carriers of urban social service industry, and are various public and service facilities for providing public service products for citizens, and specifically comprise facilities for education, medical health, culture, sports, transportation, social welfare and guarantee, administrative management and community service, postal telecommunication, commercial financial service and the like. Public service facilities play a vital role in meeting the increasing demands of people for material and cultural life and comprehensively improving the quality of life. Therefore, fair configuration of public service facilities is an important goal for planning and layout thereof, and especially in the context of policies for basic public service equalization, this goal is more important.
The existing methods for evaluating the balance of public facilities are generally realized by calculating the accessibility of space, and mainly comprise the following steps: (1) a proportional method, which calculates the ratio of the total amount of service resources (facility number, service number, etc.) to the total amount of service population, the higher the ratio is, the better the accessibility is; (2) the closest distance method, which assumes that residents always select the nearest public facilities, the closer the distance, the better the accessibility; (3) the Huff model method.
However, the proportional method cannot reveal the accessibility change in a large research unit, and is only suitable for evaluating the distribution difference between areas (such as different provinces, cities and counties) of public facilities; the nearest distance method neglects the consideration of residents on other factors such as facility scale and the like, and along with the improvement of a traffic network and a travel mode, the consideration of the distance factors in the process of selecting facilities by the residents is gradually reduced; the Huff model method is relatively abstract and difficult to understand, and meanwhile, the travel friction coefficient expressing the distance attenuation relation between the facility and the demand point is not easy to determine.
Disclosure of Invention
Aiming at the defects of the existing method for evaluating the balance of public facilities, the invention provides a method for effectively calculating the coverage radiation index of the urban public facilities such as schools, hospitals, cultural entertainment facilities and the like in order to realize the purposes of fair configuration and balanced distribution of the urban public facilities, reasonably considers the problem of trans-regional function between people and the public facilities, is suitable for both areas and areas, is intuitive and easy to understand, and can accurately evaluate the service radiation range, the service capability and the balance of the public facilities in one city.
The technical scheme adopted by the invention for solving the technical problems is as follows: a method for calculating the covering radiation index of urban public facilities is characterized by comprising the following steps: the method comprises the following steps:
step S1, acquiring city basic data and importing the city basic data into a file geographic database;
step S2, sparse population space point data is extracted;
step S3, constructing a network data set;
step S4, creating an OD cost matrix;
step S5, splitting the field and the associated attribute;
step S6, calculating the ratio of the service supply quantity to the resident demand quantity of the service point according to the service position;
step S7, calculating a ratio of service supply quantity to resident demand quantity of the resident according to the location of the resident;
and step S8, summarizing the urban public facility coverage radiation indexes and drawing a public facility coverage radiation index map.
In step S1, the city base data includes the following spatial data: the system comprises urban public facilities, roads, residential area boundaries, population and administrative divisions, wherein the file geographic database is a file geographic database of ArcGIS software.
In step S2, the population meeting the public facility action condition is screened out, and is divided into blocks according to the minimum residential area boundary region, and the residential area boundary region is divided into a plurality of grids, and all the population points of each grid are gathered to the central point of the grid.
The step S2 specifically includes the following steps:
(1) screening population data with the age within a certain interval;
(2) dividing a population element data set according to the minimum residential area boundary region;
(3) dividing each residential area boundary area into a plurality of quadrangles;
(4) finding all human mouth points contained in each quadrilateral region R;
(5) for the population point included in each quadrilateral region R, the central coordinate position P (X) of the quadrilateral region is calculatedR，YR) And a total population number N, wherein the total population attribute is added to the central point P and is stored.
The step S3 specifically includes the following steps:
(1) creating an empty set of element data;
(2) importing road data in the area into the created element data set, wherein the spatial reference of the road must be consistent with the spatial reference of the element data set;
(3) setting related parameters and constructing a network data set; the relevant parameters of the network data set are set as follows: setting a network data set name; establishing connectivity at the superposition end point of the line element without constructing a turning type, using an elevation field and establishing a driving direction; the attribute "Length" is added, with the unit set to meters.
The step S4 specifically includes the following steps:
(1) acquiring a network data set class;
(2) setting the name of an OD cost matrix layer;
(3) setting the impedance cost attribute as Length and the search range as 3000 meters;
(4) setting a starting point as a residential area element class and setting a destination as a service point element class;
(5) solving a network data set;
the form of a Name field attribute value in a line element type attribute in the OD cost matrix layer is r-s, wherein r is the Name of a residential area, and s is the Name of a service point; the attributes of the residential area element class and the service point element class include a Name field, and the service point element class attribute table includes a required supply amount attribute.
The supply quantity attribute in the service point element class comprises the number of jobs required by the operation of the urban public facility, the number of the urban public facilities and the area of the urban public facility, which index is specifically used is determined according to the type of the urban public facility, and if the urban public facility is a school or a hospital, the supply quantity attribute in the service point element class can be one or more of the number of jobs (the number of service persons) required by the operation of the urban public facility and the area of the urban public facility; in case the urban public facility is a cultural entertainment facility, the supply amount attribute in the service point element class may be one or more of the number of persons served for the operation of the urban public facility, the number of urban public facilities, and the area of the urban public facility.
The method for calculating the supply quantity in the service point element class comprises the following steps:
if the supply quantity attribute adopts an index, the supply quantity
Figure BDA0001242231550000041
Where C is the supply, X is the value of the index, and X is the sum of the values of the index in the city.
If the supply quantity attribute adopts multiple indexes, the supply quantity
Figure BDA0001242231550000042
In the formula, C is supply quantity, n is the number of supply quantity attribute indexes, xi is the value of the ith index, and X is the sum of the values of the ith index in the city.
In step S5, the Name field in the line element class attribute of the OD cost matrix is divided into two parts, and the attribute values of the two parts of fields are associated with the line element class; the step S5 specifically includes the following steps:
(1) saving line element class L in the OD cost matrix layer;
(2) adding a field R and a field S to the line element class L;
(3) the attribute value of the Name field of the split line element class L is r and s;
(4) newly adding field attribute assignment, assigning a residential area name R as R, and assigning a service point name S as S;
(5) the line element class L and the residential area element class are associated according to the attribute value of the field R, and the line element class L and the service point element class are associated according to the attribute value of the field S. The purpose of splitting the fields and associated attributes is to extract useful information from the OD cost matrix map layers, the residential site element classes, and the service point element classes and establish the connections between them.
In step S6, the resident demand in each service point threshold range is counted separately, the resident demand is connected to the line element class L attribute table, and the ratio between the service supply amount and the resident demand is calculated by the utility location; the step S6 specifically includes the following steps:
(1) adding field resident demand and field supply demand ratio SR in the service point element class attribute table;
(2) for each service point element F in the service point element classiReading service Point element FiService point name NiAnd supply volume C of urban public facilitiesi；
(3) For the field S of the line element class L of the OD cost matrix, inquiring the attribute value in the field S and the N of the service point element according to the attributeiAll element sets Fs with the same value;
(4) for each element Fs in FsjAcquiring FsjNumber of people served P within a threshold range of servicej；
(5) Calculating and storing the resident demand, wherein the calculation formula of the resident demand is as follows:
Figure BDA0001242231550000051
in the formula, SiThe resident demand of the ith service point, n is the total number of elements in Fs, PjThe number of persons served by the jth element in Fs;
(6) calculating and storing a ratio of service supply quantity to resident demand quantity within a service point threshold range, wherein the calculation formula of the ratio of service supply quantity to resident demand quantity within the service point threshold range is as follows:
Figure BDA0001242231550000052
in the formula, SRiSupply-to-demand ratio for the ith service Point, CiIs the supply volume of the ith service Point, SiThe resident demand amount of the ith residential point.
In step S7, summarizing ratio data of service supply and resident demand according to the resident position, and obtaining a result of radiation coverage index of the public facility to all residents; the step S7 specifically includes the following steps:
(1) adding a field coverage radiation index A in a residential area element class attribute table;
(2) for each residential area element F of the residential area element classiReading FiName of residential area Ni；
(3) For the field R of the line element class L of the OD cost matrix, the attribute value in the field R and the N of the residential area element are inquired according to the attributeiAll element sets Fs with the same value;
(4) for each element Fs of FsjAcquiring FsjSupply demand ratio SR ofj；
(5) Calculating FiThe covering radiation index is calculated by the formula:
Figure BDA0001242231550000061
in the formula, AiIs the covering radiation index of the ith element in Fs, n is the total number of elements in Fs, SRjIs the supply demand ratio of the jth element in Fs.
Preferably, in step S3, the network data set is created by Arc Engine redevelopment or created by tools of ArcGIS software;
in step S4, the OD cost matrix is constructed by Arc Engine redevelopment or created by tools of ArcGIS software;
in step S8, the utility coverage radiance index thematic map is drawn using ArcMap software of the ArcGIS platform.
The invention has the beneficial effects that: the invention reflects the balance degree of the urban public facilities by calculating the covering radiation indexes of the urban public facilities such as schools, hospitals, cultural facilities and the like in a certain radius range, and embodies the service capability of the urban public facilities in a certain space range.
The invention provides a method for rarefying population space point-like data, wherein service object groups of different age groups are screened out according to specific public facility types in the rarefying method, invalid redundant data are eliminated, and a calculation result is more accurate. The human mouth is partitioned according to the minimum administrative area boundary, and the memory occupation during data processing is reduced. All the population points in each quadrangle are gathered to one point, so that the problem of internal population space distribution is solved, the data volume is further reduced under the condition of ensuring the precision, and the calculation efficiency is improved.
The invention provides a method for calculating the coverage radiation index of urban public facilities, which comprises the steps of firstly calculating the ratio of the service supply quantity and the resident demand quantity of a service point by utilizing the relationship between the service point element class and an OD cost matrix line image layer, then finding out all service points corresponding to each resident on the basis, summarizing the ratio data of the service supply quantity and the resident demand quantity according to the resident positions, and then obtaining the coverage radiation index result of the public facilities to all residents.
Compared with the existing method for evaluating the balance of public facilities, the method makes full use of the existing data of the sixth census and the first geographical national condition census (other data sources can be used as long as the data requirements are met), considers the complex interaction between supply and demand in different areas, and is suitable for both areas and units in the areas. Moreover, the search of public facilities by the method crosses administrative boundaries, and the problem of cross-regional interaction between people and urban public facilities is more reasonably considered.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a flow chart of a method for sparse spatial punctate data of a population according to the present invention;
FIG. 3 is a flow chart of a method of calculating a ratio of service supply to resident demand according to the present invention;
FIG. 4 is a flow chart of a method for thinning out spatial point-like data of a population according to the present invention.
Detailed Description
In order to clearly explain the technical features of the present invention, the following detailed description of the present invention is provided with reference to the accompanying drawings. The following disclosure provides many different embodiments, or examples, for implementing different features of the invention. To simplify the disclosure of the present invention, the components and arrangements of specific examples are described below. Furthermore, the present invention may repeat reference numerals and/or letters in the various examples. This repetition is for the purpose of simplicity and clarity and does not in itself dictate a relationship between the various embodiments and/or configurations discussed. It should be noted that the components illustrated in the figures are not necessarily drawn to scale. Descriptions of well-known components and processing techniques and procedures are omitted so as to not unnecessarily limit the invention.
Aiming at the defects of the existing method for evaluating the balance of public facilities, in order to realize the purposes of fair configuration and balanced distribution of urban public facilities, the invention reasonably considers the problem of trans-regional action between people and public facilities, is suitable for both areas and areas, is intuitive and easy to understand, and can accurately evaluate the service radiation range, the service capacity and the balance of the public facilities in one city.
As shown in fig. 1, a method for calculating an urban public service coverage radiation index according to the present invention is mainly used for coverage measurement to determine urban public service deficient areas so as to help management departments to better define public service deficient areas, make proper planning and effective public policy, and comprises the following steps:
step S1, acquiring city basic data and importing the city basic data into a file geographic database;
step S2, sparse population space point data is extracted;
step S3, constructing a network data set;
step S4, creating an OD cost matrix;
step S5, splitting the field and the associated attribute;
step S6, calculating the ratio of the service supply quantity to the resident demand quantity of the service point according to the service position;
step S7, calculating a ratio of service supply quantity to resident demand quantity of the resident according to the location of the resident;
and step S8, summarizing the urban public facility coverage radiation indexes and drawing a public facility coverage radiation index map.
The invention includes schools, hospitals, cultural entertainment facilities and the like, and the method for calculating the coverage radiation index of the urban public facilities is described in detail by taking primary schools as an example. As shown in fig. 1 to 4, a method for calculating an elementary school covering radiation index according to the present invention includes the steps of:
and S1, acquiring spatial data such as elementary schools, roads, residential area boundaries, population and administrative districts, and importing the data into a file geographic database of ArcGIS software.
S2, thinning population punctate data;
because the amount of population punctiform spatial data is huge, the memory consumption is high, and the calculation time is long, the population punctiform spatial data needs to be thinned, so that the program efficiency is improved. The invention uses a regional thinning method, has good data processing efficiency and can meet the precision requirement, and the thinning process is as follows as shown in figure 2:
s21, screening population data with the age within a certain range (such as 6-12 years selected by primary schools);
s22, dividing the population element data set according to the minimum boundary area (such as streets, cells and the like);
s23, dividing a certain boundary area into a plurality of quadrangles, and traversing all the boundary areas;
s24, finding all human mouth points contained in a certain quadrangle R, and traversing all quadrangles;
and S25, calculating the central coordinate position P (XR, YR) and the total population number N of the population points contained in a certain quadrangle R, storing and traversing all quadrangles.
S3, constructing a network data set;
arc Engine is a simple, application independent, Arc Objects programming environment for building a complete class library of embedded GIS components for custom applications. The invention constructs a network data set based on the interface provided by Arc Engine.
S31, creating an element data set;
and S32, importing Road data. Importing road data in the area into the created element data set, wherein the spatial reference of the road must be consistent with the spatial reference of the element data set;
and S33, setting related parameters and constructing a network data set.
The network data set is constructed according to the elements which serve as network sources and have the connectivity strategies and the attributes associated with the network sources, and the specific construction process is as follows:
s331, acquiring a working space of road data;
s332, newly building a network data set class in the working space, and setting a network data set name;
s333, adding road element classes participating in network data set;
s334, establishing connectivity at the superposition end point of the online element without constructing a turning type, using an elevation field and establishing a driving direction in the network data set;
and S335, adding an attribute 'Length' in the network data set, wherein the unit is set to be meter.
S4, creating an OD cost matrix;
an origin-destination (OD) cost matrix is used to find and measure the minimum cost path in the network from multiple origination points to multiple destination points, and the present invention uses length cost, which is created based on the interface provided by Arc Engine. The creation process is as follows:
s41, acquiring a network data set class;
s42, setting the name of the OD cost matrix layer;
s43, setting the impedance cost attribute as Length and the search range as 3000 meters;
s44, setting the starting point as a residential area element class, setting the destination as a service point element class, wherein the attributes of the residential area element class and the service point element class both contain a Name field;
and S45, solving the network data set.
S5, splitting fields and associated attributes;
the Name field attribute value in the line element class attribute of the OD cost matrix created in step S4 is in the form of r-S, where r is the Name of the residential area and S is the Name of the service point, whereupon the Name field is split into two parts. The split field and associated attributes are intended to extract useful information from the OD cost matrix map layer, the residential site element class and the service point element class and establish a connection between them. The method comprises the following steps:
s51, saving line element class L in the OD cost matrix layer;
s52, adding a field R and a field S to the line element class L;
s53, the attribute value of the Name field of the split line element class L is r and S;
s54, assigning the attribute of the newly added field, wherein the R is assigned as R, and the S is assigned as S;
s55, the line element class L and the residential area element class are associated according to the attribute value of the field R (i.e., the residential area name), and L and the service point element class are associated according to the attribute value of the field S (i.e., the service point name).
S6, calculating the ratio of service supply quantity to resident demand quantity according to the service position;
s61, adding a field supply demand ratio SR in the service point element class attribute table;
s62, for a certain service point element F in the service point element classiReading FiService point name NiAnd supply amount Ci；
S63, for the field S of the line element class L of the OD cost matrix, inquiring the sum N thereof according to the attributeiAll elements Fs of the same value;
s64, for some element Fs in FsjAcquiring FsjNumber of people served P within a threshold range of servicej；
S65, repeating the step S64, and traversing all elements in the Fs;
s66, calculating the demand of residents
Figure BDA0001242231550000111
Where n is the total number of elements in Fs, PjThe number of persons served by the jth element in Fs;
s67, calculating the ratio of the service supply quantity to the resident demand quantity within the service point threshold value range
Figure BDA0001242231550000112
And S68, repeating the steps S62-S67, and traversing all the service point elements in the service point element class.
S7, summarizing the ratio of the service supply quantity to the resident demand quantity according to the resident positions;
s71, adding a field coverage radiation index A in the attribute table of the element class of the residential area;
s72, for a certain residential area element F of the residential area element classiReading FiName of residential area Ni；
S73, for the field R of the line element class L of the OD cost matrix, inquiring the sum N thereof according to the attributeiAll elements Fs of the same value;
s74, a specific element Fs of FsjAcquiring FsjSupply demand ratio SR ofj(by the step S4, the line element class L can acquire the service supply quantity-to-resident demand quantity ratio SR calculated at S6);
s75, repeating the step S73, and traversing all elements in the Fs;
s76, calculating FiCovering radiation index of
Figure BDA0001242231550000113
Where n is the total number of elements in Fs, SRjA supply demand ratio for the jth element in Fs;
s77, repeating the steps S72-S76, and traversing all the elements in the residential area element class.
S8, drawing a special map of the radiation index covered by the public facilities, and drawing the special map of the radiation index covered by the public facilities by a hierarchical color method for the ratio data of service supply quantity to resident demand quantity by utilizing the symbolic system of the layer attribute of the ArcMAP software of the ArcGIS platform
Moreover, the scope of the present application is not intended to be limited to the particular embodiments of the process, machine, manufacture, composition of matter, means, methods and steps described in the specification. As one of ordinary skill in the art will readily appreciate from the disclosure of the present invention, processes, machines, manufacture, compositions of matter, means, methods, or steps, presently existing or later to be developed, that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized according to the present invention. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufacture, compositions of matter, means, methods, or steps.

Claims (10)

1. A method for calculating the covering radiation index of urban public facilities is characterized by comprising the following steps: the method comprises the following steps:
step S1, acquiring city basic data and importing the city basic data into a file geographic database;
step S2, sparse population space point data is extracted;
step S3, constructing a network data set;
step S4, creating an OD cost matrix;
step S5, splitting the field and the associated attribute;
step S6, calculating the ratio of the service supply quantity to the resident demand quantity of the service point according to the service position;
step S7, calculating a ratio of service supply quantity to resident demand quantity of the resident according to the location of the resident;
step S8, summarizing urban public facility coverage radiation indexes and drawing a public facility coverage radiation index map;
the step S3 specifically includes the following steps:
(1) creating an empty set of element data;
(2) importing road data in the area into the created element data set, wherein the spatial reference of the road must be consistent with the spatial reference of the element data set;
(3) setting related parameters and constructing a network data set;
the step S4 specifically includes the following steps:
(1) acquiring a network data set class;
(2) setting the name of an OD cost matrix layer;
(3) setting the impedance cost attribute as Length and the search range as 3000 meters;
(4) setting a starting point as a residential area element class and setting a destination as a service point element class;
(5) solving a network data set;
in step S5, the Name field in the line element class attribute of the OD cost matrix is divided into two parts, and the attribute values of the two parts of fields are associated with the line element class;
in step S6, the resident demand in each service point threshold range is counted separately, the resident demand is connected to the line element class L attribute table, and the ratio between the service supply amount and the resident demand is calculated by the utility location;
in step S7, the ratio data of the service supply amount and the resident demand amount is summarized by the resident position, and the result of the radiation coverage index of the public facility to all the residents is obtained.
2. The method of claim 1, wherein in step S1, the city basic data comprises the following spatial data: the system comprises urban public facilities, roads, residential area boundaries, population and administrative divisions, wherein the file geographic database is a file geographic database of ArcGIS software.
3. The method as claimed in claim 2, wherein in step S2, the population meeting the requirement of a public facility is screened out, and is divided into blocks according to the minimum residential area boundary region, and the residential area boundary region is divided into a plurality of grids, and all population points of each grid are gathered to the central point of the grid.
4. The method for calculating the radiation index of urban public facility coverage according to claim 3, wherein the step S2 comprises the following steps:
(1) screening population data with the age within a certain interval;
(2) dividing a population element data set according to the minimum residential area boundary region;
(3) dividing each residential area boundary area into a plurality of quadrangles;
(4) finding all human mouth points contained in each quadrilateral region R;
(5) for the population point included in each quadrilateral region R, the central coordinate position P (X) of the quadrilateral region is calculatedR，YR) And a total population number N, wherein the total population attribute is added to the central point P and is stored.
5. The method for calculating the radiation index of urban public facility coverage according to claim 4, wherein the relevant parameters of the network data set are set as follows: setting a network data set name; establishing connectivity at the superposition end point of the line element without constructing a turning type, using an elevation field and establishing a driving direction; the attribute "Length" is added, with the unit set to meters.
6. The method for calculating urban public facility coverage radiation index according to claim 5,
the form of a Name field attribute value in a line element type attribute in the OD cost matrix layer is r-s, wherein r is the Name of a residential area, and s is the Name of a service point; the attributes of the residential area element class and the service point element class include a Name field, and the service point element class attribute table includes a required supply amount attribute.
7. The method for calculating the radiation index of urban public facility coverage according to claim 6, wherein the step S5 comprises the following steps:
(1) saving line element class L in the OD cost matrix layer;
(2) adding a field R and a field S to the line element class L;
(3) the attribute value of the Name field of the split line element class L is r and s;
(4) newly adding field attribute assignment, assigning a residential area name R as R, and assigning a service point name S as S;
(5) the line element class L and the residential area element class are associated according to the attribute value of the field R, and the line element class L and the service point element class are associated according to the attribute value of the field S.
8. The method for calculating the radiation index of urban public facility coverage according to claim 7, wherein the step S6 comprises the following steps:
(1) adding field resident demand and field supply demand ratio SR in the service point element class attribute table;
(2) for each service point element F in the service point element classiReading service Point element FiService point name NiAnd supply volume C of urban public facilitiesi；
(3) For the field S of the line element class L of the OD cost matrix, inquiring the attribute value in the field S and the N of the service point element according to the attributeiAll element sets Fs with the same value;
(4) for each element Fs in FsjAcquiring FsjNumber of people served P within a threshold range of servicej；
(5) Calculating and storing the resident demand, wherein the calculation formula of the resident demand is as follows:
Figure FDA0002251618820000041
in the formula, SiThe resident demand of the ith service point, n is the total number of elements in Fs, PjThe number of persons served by the jth element in Fs;
(6) calculating and storing a ratio of service supply quantity to resident demand quantity within a service point threshold range, wherein the calculation formula of the ratio of service supply quantity to resident demand quantity within the service point threshold range is as follows:
Figure FDA0002251618820000042
in the formula, SRiSupply-to-demand ratio for the ith service Point, CiIs the supply volume of the ith service Point, SiThe resident demand amount of the ith residential point.
9. The method for calculating the radiation index of urban public facility coverage according to claim 8, wherein the step S7 comprises the following steps:
(1) adding a field coverage radiation index A in a residential area element class attribute table;
(2) for each residential area element F of the residential area element classiReading FiName of residential area Ni；
(3) For the field R of the line element class L of the OD cost matrix, the attribute value in the field R and the N of the residential area element are inquired according to the attributeiAll element sets Fs with the same value;
(4) for each element Fs of FsjAcquiring FsjSupply demand ratio SR ofj；
(5) Calculating FiThe covering radiation index is calculated by the formula:
Figure FDA0002251618820000051
in the formula, AiIs the covering radiation index of the ith element in Fs, n is the total number of elements in Fs, SRjIs the supply demand ratio of the jth element in Fs.
10. The method for calculating the radiation index of urban public facility coverage according to any one of claims 1 to 9,
in step S3, the network data set is created by Arc Engine redevelopment or created by tools of ArcGIS software;
in step S4, the OD cost matrix is constructed by Arc Engine redevelopment or created by tools of ArcGIS software;
in step S8, the utility coverage radiance index thematic map is drawn using ArcMap software of the ArcGIS platform.
CN201710138930.7A 2017-03-09 2017-03-09 Method for calculating urban public facility coverage radiation index Active CN106909692B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710138930.7A CN106909692B (en) 2017-03-09 2017-03-09 Method for calculating urban public facility coverage radiation index

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710138930.7A CN106909692B (en) 2017-03-09 2017-03-09 Method for calculating urban public facility coverage radiation index

Publications (2)

Publication Number Publication Date
CN106909692A CN106909692A (en) 2017-06-30
CN106909692B true CN106909692B (en) 2020-04-14

Family

ID=59187797

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710138930.7A Active CN106909692B (en) 2017-03-09 2017-03-09 Method for calculating urban public facility coverage radiation index

Country Status (1)

Country Link
CN (1) CN106909692B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609748A (en) * 2017-08-18 2018-01-19 深圳航天智慧城市系统技术研究院有限公司 A kind of method and system that Development Assessment is carried out to urban facilities
CN107563803B (en) * 2017-08-24 2020-12-08 北京工商大学 Quotient circle dividing method based on distance consumption grid
CN108428056A (en) * 2018-03-13 2018-08-21 重庆市地理信息中心 The evaluation method and device of community's life's circle
CN110852547B (en) * 2019-08-19 2023-04-07 上海元卓信息科技有限公司 Public service facility grading method based on position data and clustering algorithm
CN111553566A (en) * 2020-04-10 2020-08-18 东莞市城建规划设计院 Method for defining service range of urban public service facility
CN111709642A (en) * 2020-06-16 2020-09-25 北京市测绘设计研究院 Method for evaluating layout balance of barrier-free facilities in city
CN112036640B (en) * 2020-08-31 2022-10-28 河南大学 Division method and device for public facility service area
CN112396289A (en) * 2020-10-12 2021-02-23 中国兵器科学研究院 Equipment system logistics support capability assessment method, device, equipment and medium
CN113011768B (en) * 2021-03-31 2021-12-03 中国城市建设研究院有限公司 Public facility data processing method, system, electronic device and medium
CN113435769A (en) * 2021-07-07 2021-09-24 北京市测绘设计研究院 Method for calculating barrier-free facility index of city
CN116467788B (en) * 2023-05-10 2023-11-21 悉地（苏州）勘察设计顾问有限公司 Barrier-free environment construction management method and system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002157209A (en) * 2000-11-17 2002-05-31 Dorikomu:Kk Retrieval system using three-dimensional virtual space
JP2003044864A (en) * 2001-07-27 2003-02-14 Techno Com:Kk Device for processing regional information and method for displaying regional information
CN102750411A (en) * 2012-06-19 2012-10-24 中国地质大学（武汉） Urban dynamic micro-simulation method based on multi-agent discrete choice model
CN102750427A (en) * 2012-07-23 2012-10-24 中国地质大学（武汉） Urban traffic planning simulation method and system based on multi-agent activity model
CN102779453A (en) * 2012-06-15 2012-11-14 清华大学 Demonstration device and demonstration method for urban public security multi-disaster secondary derivation coupling simulation
CN103455856A (en) * 2013-09-12 2013-12-18 中国科学院地理科学与资源研究所 Technical method for urban space function identification

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002157209A (en) * 2000-11-17 2002-05-31 Dorikomu:Kk Retrieval system using three-dimensional virtual space
JP2003044864A (en) * 2001-07-27 2003-02-14 Techno Com:Kk Device for processing regional information and method for displaying regional information
CN102779453A (en) * 2012-06-15 2012-11-14 清华大学 Demonstration device and demonstration method for urban public security multi-disaster secondary derivation coupling simulation
CN102750411A (en) * 2012-06-19 2012-10-24 中国地质大学（武汉） Urban dynamic micro-simulation method based on multi-agent discrete choice model
CN102750427A (en) * 2012-07-23 2012-10-24 中国地质大学（武汉） Urban traffic planning simulation method and system based on multi-agent activity model
CN103455856A (en) * 2013-09-12 2013-12-18 中国科学院地理科学与资源研究所 Technical method for urban space function identification

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
城市公共服务设施布局的均衡性研究-以北京市城六区医疗设施为例;刘静等;《城市发展研究》;20160526;正文第6-11页 *
面向智慧城市建设的居民公共服务需求研究-以河北省石家庄市为例;赵勇等;《地理科学进展》;20150429;正文第473-481页 *

Also Published As

Publication number Publication date
CN106909692A (en) 2017-06-30

Similar Documents

Publication Publication Date Title
CN106909692B (en) Method for calculating urban public facility coverage radiation index
CN103971184B (en) Transmission line of electricity path generating method based on GIS Spatial GISs
Biba et al. A new method for determining the population with walking access to transit
CN103984997B (en) Power transmission engineering addressing selection method based on GIS spatial informations
CN105677804A (en) Determination of authority stations and building method and device of authority station database
WO2017101277A1 (en) City cognitive map generating method based on internet word frequency
WO2022198963A1 (en) Big data-based commercial space quality evaluation method and system, device, and medium
CN109840452A (en) A kind of grid covering scene automatic identifying method and device
CN111709642A (en) Method for evaluating layout balance of barrier-free facilities in city
CN107506433A (en) Urban development space general layout Scene Simulation system
CN109146360A (en) Grid method for building up and device and allocator and device
CN114881550B (en) Resource equalization allocation method and system based on spatial information
Qian et al. Incorporating landscape diversity into greenway alignment planning
CN110807287B (en) Construction method for dividing service grids in emergency repair point arrangement scene based on finite element grids
CN106951453A (en) A kind of geographical entity coding method of quick renewal and data sharing
Yuan et al. From design to digital model: A quantitative analysis approach to Garden Cities theory
CN114154880A (en) Method for evaluating walking accessibility of public service facility on street unit level
CN110413855A (en) A kind of region entrance Dynamic Extraction method based on taxi drop-off point
CN111881573B (en) Population space distribution simulation method and device based on urban inland inundation risk assessment
Zheng et al. Spatial Accessibility of Urban Green Space Based on Multiple Research Scales: A Case Study of Futian District, Shenzhen.
CN110458309A (en) A kind of net based on practical road network environment about share-car bus station position method
Srikanth et al. Complexity science-based spatial performance analyses of UNStudio/DP Architects’ SUTD Campus and WOHA's Kampung Admiralty
Dhamaniya Methodology of maintenance criteria of PMGSY roads in second phase of planning-a case study
CN114723316B (en) Reachability evaluation method and system for urban public facilities based on GIS
CN117079124B (en) Urban and rural landscape image quantification and promotion method based on community differentiation

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