CN112948769B - City circle range determining method and system based on commuting big data - Google Patents

City circle range determining method and system based on commuting big data Download PDF

Info

Publication number
CN112948769B
CN112948769B CN202110260787.5A CN202110260787A CN112948769B CN 112948769 B CN112948769 B CN 112948769B CN 202110260787 A CN202110260787 A CN 202110260787A CN 112948769 B CN112948769 B CN 112948769B
Authority
CN
China
Prior art keywords
administrative
region
area
urban
order
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
CN202110260787.5A
Other languages
Chinese (zh)
Other versions
CN112948769A (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.)
China Academy Of Urban Planning & Design
Original Assignee
China Academy Of Urban Planning & Design
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 China Academy Of Urban Planning & Design filed Critical China Academy Of Urban Planning & Design
Priority to CN202110260787.5A priority Critical patent/CN112948769B/en
Publication of CN112948769A publication Critical patent/CN112948769A/en
Application granted granted Critical
Publication of CN112948769B publication Critical patent/CN112948769B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Business, Economics & Management (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Tourism & Hospitality (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Educational Administration (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Biology (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Remote Sensing (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method and a system for determining urban circle range based on commuting big data, which relate to the technical field of urban and regional space range quantitative analysis. The method can avoid the defect that the planning range of the urban circle is too generalized or the daily commute circle is emphasized on one side, and further accurately define the urban circle.

Description

City circle range determining method and system based on commuting big data
Technical Field
The invention relates to the technical field of quantitative analysis of urban and regional spatial ranges, in particular to a method and a system for determining urban circle range based on commuting big data.
Background
The cultivation of the urban ring which forms the commuting efficient and integrated development is the key for improving the gathering power, the radiation power, the innovation power and the competitiveness of the city, is a strong engine for promoting the regional coordinated development, and is a duplicator for stably increasing, expanding the internal demand and building up the national competitive advantage.
In the process of Chinese urbanization, due to the formation of the peripheral secondary center of the central city, the networked spatial pattern metropolitan area space definition is necessarily promoted to consider the evolution trend from the first-order abdominal land (central city commuting area) to the second-order abdominal land (peripheral radiation expansion area).
The tight commute relationship is a key element in the spatial extent definition of the metro ring. But the centripetal commuting rate commonly used in foreign countries and the 1-hour commuting circle of the mainstream of China define standards, and cannot be suitable for the characteristics of Chinese towns. Through big data demonstration research on Beijing, shanghai, guangzhou, chengdu, wuhan, xian, zheng Zhou and 3.2 hundred million commuting crowds, the data discovery: (1) The peripheral group of the Chinese city has relatively perfect city functions, the internal commuting proportion is more than 80-90%, the centripetal commuting rate is less than 10%, and large-range fluctuation exists, so that a definition standard is difficult to set; (2) The definition of the commuting circle of one hour is too extensive, 1+ X cities are often met, and the range of hundreds of thousands of square kilometers is greatly different from the concept of the urban circle; at present, in the main cities in China, the space scale of 1-hour commuting does not exceed 25 kilometers, and the urban circle of the city cannot be accurately defined by the existing urban circle range determining method.
Administrative management, resource management and government governance in China do not have a management level matched with a metropolitan area, so that the development of a region which develops township to a certain extent needs to be planned overall by taking a central city as a core. The calculation method of the urban area can guide the objective definition of the space range, which is particularly important for whether the resultant force and the collaborative development can be formed between cities in the urban area.
Disclosure of Invention
The invention aims to provide a method and a system for determining the range of a metro ring based on commuting big data, which can measure the internal relation to a peripheral grouping of the metro ring and the two-stage gradient of the range of the metro ring in the abdominal region, thereby more accurately determining the range of the metro ring.
In order to achieve the purpose, the invention provides the following scheme:
a metropolitan area range determination method based on commute big data comprises the following steps:
determining commuting OD data of a research city; the commuting OD data comprise resident residence information and resident employment place information of which at least one end is positioned in a study city;
determining a research area according to the commuting OD data, and dividing the research area according to map information of the research area to obtain an administrative area set contained in a metro ring potential space range; the administrative region set comprises a plurality of administrative regions; the research area is an administrative area covered by a resident residence or a resident employment place in the commuting OD data;
determining the central urban area of the study city as the central urban area range of the urban circle;
calculating a network connectivity index between any two administrative areas in the research area according to the commuting OD data;
calculating a first-order comprehensive index of the network connectivity of the abdominal region and the urban region of each first set administrative region, and determining a first-order abdominal region range of the urban ring according to the first-order comprehensive index of the network connectivity of the abdominal region and the urban region; the first set administrative region is each administrative region located in the research region range; the first-order abdominal region and urban region network connectivity comprehensive index is the sum of network connectivity indexes of each administrative region contained in the first set administrative region and the central urban region range;
determining closely-connected administrative regions of each second set administrative region; the tight connection administrative areas are administrative areas which are associated with the second set administrative area and have network connectivity indexes larger than a tight connection set threshold value in the research area; the second set administrative region is an administrative region which is within the research region and is not located within the central urban region range and the first-order abdominal region range;
calculating a second-order abdominal region and urban area network connectivity comprehensive index corresponding to each second set administrative region, and judging whether each second set administrative region belongs to a second-order abdominal standard or not according to the second-order abdominal region and urban area network connectivity comprehensive index so as to obtain a second-order abdominal range of the urban area; the second-order comprehensive index of the network connectivity of the abdominal region is the sum of the second set administrative region and the first-order comprehensive index of the network connectivity of the abdominal region of each administrative region closely connected with the second set administrative region;
and converging the central urban area range, the first-order abdominal region range and the second-order abdominal region range to obtain an urban circle space range.
Optionally, the research area is divided according to the map information of the research area to obtain an administrative area set included in the metro area potential space range, and the method specifically includes:
and dividing the research area according to the district-county administrative boundary in the map information of the research area to obtain an administrative area set contained in the urban circle potential space range.
Optionally, the determining the central urban area of the research city as the central urban area range of the urban circle specifically includes:
determining the central area of the research city as a preliminary central urban area of the urban circle according to the total planning range of the current city and the planning and regulated range of the territorial space;
determining commuting demographic distribution information for the research city;
and locally adjusting the preliminary central urban area according to the district administrative boundary in the map information of the research city and the commuting population distribution information so as to obtain the final central urban area range of the urban area.
Optionally, the calculating, according to the commute OD data, a network connectivity index between any two administrative areas in the research area specifically includes:
according to the commute OD data and formula
Figure GDA0004015433540000031
Calculating a network connectivity index between any two administrative areas in the research area;
wherein R is ij The network connectivity of the ith administrative area and the jth administrative interval;
Figure GDA0004015433540000032
Figure GDA0004015433540000033
in order to standardize the commuting connection strength index of the employment from the ith administrative district to the jth administrative district,
Figure GDA0004015433540000034
the commuting contact strength index of employment from the jth administrative district to the ith administrative district after standardization;
Figure GDA0004015433540000035
V ij the number of commutes from the i th administrative district to the j th administrative district,
Figure GDA0004015433540000036
the commute population is the sum of employment populations from the ith administrative area to the administrative areas, and n represents the total number of the administrative areas;
Figure GDA0004015433540000041
V ji the number of commutes from the j administrative district to the i administrative district,
Figure GDA0004015433540000042
the commute population is the sum of employment populations from the jth administrative area to each administrative area; max (R) ij ) The maximum value of the network connectivity in the commuting network of all administrative areas is obtained;
Figure GDA0004015433540000043
the network connectivity indexes of the ith administrative area and the jth administrative interval are normalized.
Optionally, the determining the first-order abdominal area range of the urban circle according to the first-order abdominal area urban area network connectivity comprehensive index specifically includes:
sequentially judging whether the first-order abdominal region urban area network connectivity comprehensive index is larger than a first-order abdominal region judgment threshold value or not;
and if so, determining a first-order abdominal region by a first set administrative region corresponding to the first-order comprehensive index of network connectivity of the urban area of the abdominal region, and further determining a first-order abdominal region range of the urban area.
Optionally, the second-order abdominal region standard is determined according to the second-order abdominal region urban area network connectivity comprehensive index, and then a second-order abdominal region range of the urban circle is obtained, specifically including:
sequentially judging whether the second-order abdominal region urban network connectivity comprehensive index is greater than a second-order abdominal region judgment threshold value;
and if so, determining a second set administrative district corresponding to the second-order abdominal region urban network connectivity comprehensive index as a second-order abdominal region, and further determining a second-order abdominal region range of the urban district.
A commute big data based metro area determination system comprising:
the data acquisition module is used for determining commuting OD data of a research city; the commuting OD data comprise resident residence information and resident employment place information of which at least one end is positioned in a research city;
the administrative region set determining module is used for determining a research region according to the commuting OD data, and dividing the research region according to map information of the research region to obtain an administrative region set contained in a metro ring potential space range; the administrative region set comprises a plurality of administrative regions; the research area is an administrative area covered by a resident residence or a resident employment place in the commuting OD data;
the central urban area range determining module is used for determining the central urban area of the research city as the central urban area range of the urban circle;
the network connectivity index calculation module is used for calculating a network connectivity index between any two administrative areas in the research area according to the commuting OD data;
the first-order abdominal region range determining module is used for calculating a first-order abdominal region and urban area network connectivity comprehensive index of each first set administrative region and determining a first-order abdominal region range of the urban area according to the first-order abdominal region and urban area network connectivity comprehensive index; the first set administrative region is each administrative region located in the research region range; the first-order abdominal region and urban region network connectivity comprehensive index is the sum of the network connectivity indexes of the first set administrative region and the administrative regions contained in the central urban region range;
the close connection administrative region determining module is used for determining the close connection administrative region of each second set administrative region; the closely-connected administrative regions are administrative regions which are associated with the second set administrative region and have network connectivity indexes larger than a closely-connected set threshold value in the research region; the second set administrative region is an administrative region which is in the research region and is not positioned in the central urban region range and the first-order abdominal region range;
the second-order abdominal region determining module is used for calculating a second-order abdominal region and urban area network connectivity comprehensive index corresponding to each second set administrative region, and judging whether each second set administrative region belongs to a second-order abdominal region standard or not according to the second-order abdominal region and urban area network connectivity comprehensive index so as to obtain a second-order abdominal region range of the urban area; the second-order comprehensive index of the network connectivity of the abdominal region and the urban area is the sum of the second set administrative area and the first-order comprehensive index of the network connectivity of the abdominal region and the urban area of each administrative area tightly connected with the second set administrative area;
and the urban circle space range determining module is used for converging the central urban area range, the first-order abdominal region range and the second-order abdominal region range to obtain an urban circle space range.
Optionally, the administrative region set determining module specifically includes:
and the administrative region set determining unit is used for dividing the research region according to the district-county administrative boundary in the map information of the research region to obtain an administrative region set contained in the urban area potential space range.
Optionally, the central urban area range determining module specifically includes:
the preliminary central urban area determining unit is used for determining the central area of the research city as the preliminary central urban area of the urban area according to the total planning range of the current city and the planning and regulation range of the territorial space;
a commute population distribution information determining unit, configured to determine commute population distribution information of the research city;
and the central urban area range determining unit is used for locally adjusting the initial central urban area according to the district and county administrative boundary in the map information of the research city and the commuting population distribution information so as to obtain the final central urban area range of the urban area.
Optionally, the network connectivity index calculating module specifically includes:
a network connectivity index calculation unit for calculating the network connectivity index according to the commuting OD data and the formula
Figure GDA0004015433540000061
Calculating a network connectivity index between any two administrative areas in the research area;
wherein R is ij The network connectivity of the ith administrative area and the jth administrative interval;
Figure GDA0004015433540000062
Figure GDA0004015433540000063
in order to standardize the commuting contact strength index of employment from the ith administrative district to the jth administrative district,
Figure GDA0004015433540000064
the commuting contact strength index of employment from the jth administrative district to the ith administrative district after standardization;
Figure GDA0004015433540000065
V ij the number of commuters who live in the jth administrative district from the ith administrative district,
Figure GDA0004015433540000066
the commute population is the sum of employment populations from the ith administrative area to the administrative areas, and n represents the total number of the administrative areas;
Figure GDA0004015433540000067
V ji the number of commutes from the j administrative district to the i administrative district,
Figure GDA0004015433540000068
the commute population sum from the j administrative district to the employment of each administrative district; max (R) ij ) The maximum value of the network connectivity in the commuting network of all administrative areas is obtained;
Figure GDA0004015433540000069
the network connectivity indexes of the ith administrative area and the jth administrative interval are standardized.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention determines the central urban area range, the first-order abdominal area range and the second-order abdominal area range based on the commuting OD data and the network connectivity index, avoids the defect that the urban circle planning range is too generalized or the daily commuting circle is emphasized on one side, and can accurately define the urban circle space range.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments will be briefly described below. The drawings in the following description are only some embodiments of the invention, and it is obvious to a person skilled in the art that other drawings can be obtained on the basis of these drawings without inventive step.
FIG. 1 is a flow chart of a metropolitan area range determination method based on commute big data according to the present invention;
FIG. 2 is a block diagram of a metropolitan area range determination system based on commute big data according to the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a method and a system for determining the urban circle range based on commuting big data, which can accurately define the urban circle.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
The urban ring is a urbanized space form which takes one or more urban central urban areas with certain regional influence and radiation force as a core, comprises peripheral secondary central cities, new cities, a batch of medium and small cities and small towns and spans one or more administrative bodies. And each level of towns in the urban area have close population, economy, traffic and information connection.
Commuter OD data: in the 1-6 months in 2019, information of resident residence places and information of resident employment places obtained by Baidu map location service and mobile phone data of mobile communication operators are converged by district, county and administrative districts to form commuting OD data, and two mutually independent data sources are subjected to cross verification and fusion analysis.
The urban circle emphasizes the regional scope of the complete urban function, often characterized by tight commute connections and internal job-hold balances. Research finds that the urban circle space of the main city of China is composed of three circle layers of a central city area, a central city commuting area and a peripheral radiation expansion area, and has the following characteristics: (1) The central city commuting area is in a space range closely connected with the central city area, the maximum space scale is 40 kilometers at present, but 20% -30% of the commuting area is not a complete city functional area, and the other end of the commuting area is outside the central city commuting area; (2) The expansion area which is radiated outwards by the central city commuting area has a spatial dimension of 60-80 kilometers, and more than 95% of internal work and live balance can be realized in a spatial range formed by the central city area, the central city commuting area and the peripheral radiation expansion area; (3) The 90 minutes elapsed is the upper limit of commute duty and is also a boundary for the metro ring space to extend.
In view of this, the present embodiment provides a city circle range determination method based on commute big data, as shown in fig. 1, including the following steps.
Step 101: determining commuting OD data of a research city; the commuting OD data takes commuting population with at least one end positioned in the administrative boundary range of the urban region of the central city as a research object, and the data content consists of resident residence information and resident employment information; the research city is a central city which has a certain population scale and economic total amount in an urban circle, has a certain area radiation influence and bears the area comprehensive service function in a certain range; the research city includes resident residences and/or resident employment sites.
Step 102: determining a research area according to the commuting OD data, and dividing the research area according to map information of the research area to obtain an administrative area set contained in a metro ring potential space range; the administrative region set comprises a plurality of administrative regions; the study area is an administrative area covered by the resident residence or the resident employment place in the commuting OD data.
Regarding the determination of the statistical units, there are generally two types of statistical units. One is that the county is used as a basic unit, such as countries in Europe, america and the like; the other is a village and town as a basic unit, such as taiwan.
If the county-level administrative districts are used as basic units, county-level administrative boundaries do not need to be broken, and as long as the overall level of the county meets the standards of the peripheral counties of the metropolitan circle, the whole county is brought into the range of the metropolitan circle. The method has the advantages of taking care of the integrity of the administrative district and causing no management confusion. The method has the defect that for some counties with irregular shapes, particularly for counties with complicated shapes and large difference of internal economic development, a large area of rural areas which do not meet the requirements of the metropolis circle can be classified into the metropolis circle.
In view of this, the dividing process of the present embodiment is: and dividing the research area according to the district-county administrative boundary in the map information of the research area to obtain an administrative area set contained in the urban circle potential space range.
Step 103: determining the central urban area of the research city as the central urban area range of the urban circle, and specifically comprising the following steps:
firstly, determining a central area of the research city as a preliminary central urban area of an urban area according to a current urban overall planning range and a national space planning regulation range; secondly, determining commuting population distribution information of the research city; finally, according to the district-county administrative boundaries in the map information of the research city and the commuting population distribution information, the preliminary central urban area is locally adjusted to obtain the final central urban area range of the urban circle; if the proportion of residential or employment population of the district and county within the central urban area range of the urban circle exceeds 30%, the administrative administration of the district and county is brought into the central urban area range of the urban circle; and if the commuting population distribution of the administrative scope of the district and the county is more than a set value, bringing the administrative scope of the district and the county into the central urban area scope of the urban circle.
The central urban area is an urban core built-up area which has a certain population scale and economic total amount in an urban circle, has a certain area radiation influence and bears the comprehensive service function of an area in a certain range.
Step 104: and calculating the network connectivity index between any two administrative areas in the research area according to the commuting OD data.
And the standardized network connectivity index is used as the measurement label of the commuting contact degree, so that the influence caused by the scale of administrative districts and the difference of internal commuting proportion is overcome.
(1) Administrative interval commuting contact:
and constructing a commuting connection index of the administrative region, and overcoming the influence of population scale difference of different administrative regions on measuring the closeness degree of the commuting connection through standardized processing of the commuting population number of the administrative region.
(2) Network connectivity index:
the network connectivity reflects the closeness of commuting connections in the administrative district network, and a network connectivity index is obtained by adopting maximum value standardization.
Therefore, step 104 specifically includes: according to the commute OD data and formula
Figure GDA0004015433540000091
Calculating the network connectivity between any two administrative regions in the research regionAnd (4) counting.
Wherein R is ij The network connectivity of the ith administrative area and the jth administrative interval;
Figure GDA0004015433540000092
Figure GDA0004015433540000101
in order to standardize the commuting connection strength index of the employment from the ith administrative district to the jth administrative district,
Figure GDA0004015433540000102
the commuting contact strength index of employment from the jth administrative district to the ith administrative district after standardization;
Figure GDA0004015433540000103
V ij the number of commutes from the i th administrative district to the j th administrative district,
Figure GDA0004015433540000104
the sum of commuting population from the i th administrative district to employment of each administrative district, wherein n represents the total number of the administrative districts;
Figure GDA0004015433540000105
V ji the number of commuters who live from the jth administrative district to the ith administrative district,
Figure GDA0004015433540000106
the commute population is the sum of employment populations from the jth administrative area to each administrative area; max (R) ij ) The maximum value of the network connectivity in the commuting network of all administrative areas is obtained;
Figure GDA0004015433540000107
the network connectivity indexes of the ith administrative area and the jth administrative interval are normalized.
Taking the example of the research city being beijing, the residential area being the Tongzhou area, the employment area being the sunny area, the eastern area, the western area, the hai lake area, the Fengtai area and the Shijing mountain area, the network connectivity index of the Tongzhou area and other administrative areas as shown in Table 1 is calculated.
TABLE 1 information relationship table between Tongzhou district and other administrative districts
Figure GDA0004015433540000108
Figure GDA0004015433540000111
Step 105: calculating a first-order comprehensive index of the network connectivity of the abdominal region and the urban region of each first set administrative region, and determining a first-order abdominal region range of the urban ring according to the first-order comprehensive index of the network connectivity of the abdominal region and the urban region; the first set administrative region is each administrative region located in the research region range; the first-order abdominal region and urban region network connectivity comprehensive index is the sum of the network connectivity indexes of the first set administrative region and the administrative regions contained in the central urban region range.
Taking the study city as beijing as an example, the first set administrative region is a transit region, the administrative regions to which the central urban region belongs include a sunny region, a east city region, a west city region, a hai lake region, a Fengtai region and a rocky mountain region, and the first-order comprehensive network connectivity index of the abdominal region corresponding to the transit region is a total number 2.51 referred to in table 1.
Then, sequentially judging whether the first-order abdominal region and urban area network connectivity comprehensive index corresponding to each first set administrative area is greater than a first-order abdominal region judgment threshold value; and if so, determining a first-order abdominal region by a first set administrative region corresponding to the first-order comprehensive index of network connectivity of the urban area of the abdominal region, and further determining a first-order abdominal region range of the urban area. Preferably, the first-order ventral discrimination threshold provided by the present embodiment is 1. And when the first set administrative region is within the range of the central urban area and the range of the first-order abdominal land, the first set administrative region is classified into the range of the central urban area.
Step 106: determining a closely-connected administrative region of each second set administrative region; the tight connection administrative areas are administrative areas which are associated with the second set administrative area and have network connectivity indexes larger than a tight connection set threshold value in the research area; the second set administrative area is an administrative area within the study area and not located within the central urban area and the first-order abdominal area.
And constructing a measurement method of the urban circle commuting relation from first order to second order by means of second-order connection of a complex network theory. It is more important to have "popular" Friends than just a large number of Friends, not just to count the number of Friends, but to count the number of more Friends that each friend can bring, tracking the Friends of the friend-called "Second-degree Friends" (Second-degree Friends).
Taking the research city as Beijing as an example, the second set administrative region includes administrative regions such as Sanhe city of Gallery city, gu' an county, xiang he county, hui nationality autonomous region of Dai factory, guangyang district, wuqing district of Tianjin city, juding city \28095, huai city, laishui county, and Zhang Koukou county, 28095, deer county, etc.
Taking the xianghe county as an example, the administrative regions associated with the xianghe county and having the network connectivity index greater than the tightly connected set threshold include the sunny region, the Tongzhou region, the Daxing region and the east city region, which are calculated according to the above formula and the commuting OD data.
Step 107: calculating a second-order abdominal region and urban area network connectivity comprehensive index corresponding to each second set administrative area, and judging whether each second set administrative area belongs to a second-order abdominal standard or not according to the second-order abdominal region and urban area network connectivity comprehensive index so as to obtain a second-order abdominal range of the urban circle; the second-order comprehensive index of the network connectivity of the abdominal region and the urban area is the sum of the second set administrative area and the first-order comprehensive index of the network connectivity of the abdominal region and the urban area of each administrative area tightly connected with the second set administrative area.
Taking the xianghe county as an example, the administrative regions tightly connected with the xianghe county are a township region, a Tongzhou region, a Daxing region and an east city region, and then the second-order comprehensive index of the abdominal region network connectivity is equal to the sum of the first-order comprehensive indexes of the township region, the Tongzhou region, the Daxing region and the east city region.
Then, judging whether the second-order abdominal region urban network connectivity comprehensive index is larger than a set second-order abdominal region judgment threshold value; and if so, determining a second set administrative district corresponding to the second-order abdominal region urban area network connectivity comprehensive index as the second-order abdominal region of the urban area. Preferably, the second-order abdominal discrimination threshold provided by the present embodiment is 1.
Step 108: converging the central urban area range, the first-order abdominal region range and the second-order abdominal region range to obtain an urban circle space range; the method specifically comprises the following steps:
aggregating the spatial ranges of administrative areas contained in the central urban area range, the first-order abdominal area range and the second-order abdominal area range, and determining the spatial range of the urban circle of the research city.
Wherein, the first-order abdominal region (urban circle commuting region) refers to a county-level administrative unit closely linked with the daily commute of the central city. The first-order abdominal land defines the range of the central city directly influencing the abdominal land, and reflects the centripetal gathering capacity and the urban spatial structure of the urban circle at the initial development stage.
Second-order abdominal areas (metro-ring radiating areas) refer to county-level administrative units closely related to the industrial functions of a central city, with spatial dimensions of a radius that may extend to 60-80 km and more, typically involving one or more neighboring districts. The second-order abdominal land defines the associated abdominal land of the urban circle, reflecting the regional radiation of the central city and the shifting of the industrial economy' reach and cross-regional interoperability.
Taking the Beijing city as an example, the method provided by the embodiment is combined to determine the Beijing city circle.
Commuting OD data: in 2019, in 1-6 months, information of resident residence places and information of employment places of residents are obtained by using Baidu map location services and mobile phone data of mobile communication operators, a commuting OD is formed by gathering district, county and administrative districts, and two mutually independent data sources are subjected to cross verification and fusion analysis.
Table 2: administrative region commuting contact basic data sample table
Figure GDA0004015433540000131
Figure GDA0004015433540000141
The determination result is: the Beijing City circle has a range of 2.7 ten thousand square kilometers and a space radius of 93 kilometers, and comprises 16 administrative districts of Beijing City, 5 administrative districts of Sanhe City, guan county, xianghe county, hui autonomous district of Dachang, guangyang county, wujin City Wuqing district, baoding City \28095, huaiyai county, 28095and Lu county of Zhongshui county, laishui county and Zhang Jiakou.
Wherein, the area of the central urban area is 1373 square kilometers, the space radius is 21 kilometers, and 53.2 percent of internal working and living balance in the area is realized; the first-order abdominal area is 7739 square kilometers, the space radius is 54 kilometers, and the first-order abdominal area comprises a Tongzhou area, a Daxing area, a mountain area, a gate head ditch area, a Chang Ping area and a cisternal area in Beijing city, and 80.1 percent of job-stop balance is realized in the range of a central urban area and a commuting area; the second-order abdominal area is 18148 square kilometers, and 99.8% of the balance of the position is achieved in the workplace.
Table 3: beijing city circle center city area space composition and commuting characteristic summary table
Figure GDA0004015433540000142
Table 4: summary table of first-order abdominal space composition and commuting characteristics of Beijing city circle
Figure GDA0004015433540000143
Figure GDA0004015433540000151
Table 5: second-order abdominal space composition and commuting characteristic summary table of Beijing city circle
Figure GDA0004015433540000152
Figure GDA0004015433540000161
To achieve the above object, the present invention further provides a system for determining a metro area based on commute big data, as shown in fig. 2, including:
a data acquisition module 201, configured to determine commuting OD data of a research city; the commuting OD data includes resident residence information and resident employment information of which at least one end is located in the research city.
An administrative region set determining module 202, configured to determine a research region according to the commuting OD data, and divide the research region according to map information of the research region to obtain an administrative region set included in a metro ring potential spatial range; the administrative region set comprises a plurality of administrative regions; the study area is an administrative area covered by the resident residence or the resident employment place in the commuting OD data.
And the central urban area range determining module 203 is used for determining the central urban area of the research city as the central urban area range of the urban circle.
And a network connectivity index calculation module 204, configured to calculate, according to the commuting OD data, a network connectivity index between any two administrative areas in the research area.
A first-order abdominal region range determining module 205, configured to calculate a first-order abdominal region and urban region network connectivity comprehensive index of each first set administrative area, and determine a first-order abdominal region range of the urban district according to the first-order abdominal region and urban region network connectivity comprehensive index; the first set administrative region is each administrative region located in the research region range; the first-order abdominal region and urban region network connectivity comprehensive index is the sum of the network connectivity indexes of the first set administrative region and the administrative regions contained in the central urban region range.
A tightly-connected administrative region determination module 206, configured to determine a tightly-connected administrative region of each second set administrative region; the tight connection administrative areas are administrative areas which are associated with the second set administrative area and have network connectivity indexes larger than a tight connection set threshold value in the research area; the second set administrative area is an administrative area within the study area and not located within the central urban area and the first-order abdominal area.
The second-order abdominal region range determining module 207 is configured to calculate a second-order abdominal region and urban area network connectivity comprehensive index corresponding to each second set administrative region, and determine whether each second set administrative region belongs to a second-order abdominal region standard according to the second-order abdominal region and urban area network connectivity comprehensive index, so as to obtain a second-order abdominal region range of the urban area; the second-order comprehensive index of the network connectivity of the abdominal region is the sum of the second set administrative region and the first-order comprehensive index of the network connectivity of the abdominal region of each administrative region closely connected with the second set administrative region.
A metropolitan area spatial range determination module 208, configured to converge the central urban area range, the first-order abdominal range, and the second-order abdominal range to obtain a metropolitan area spatial range.
The administrative region set determining module 202 specifically includes:
and the administrative region set determining unit is used for dividing the research region according to the district-county administrative boundaries in the map information of the research region to obtain an administrative region set contained in the urban circle potential space range.
The central urban area range determining module 203 specifically includes:
the preliminary central urban area determining unit is used for determining the central area of the research city as the preliminary central urban area of the urban area according to the total planning range of the current city and the planning and regulation range of the territorial space; a commute population distribution information determining unit, configured to determine commute population distribution information of the research city; and the central urban area range determining unit is used for locally adjusting the initial central urban area according to the district and county administrative boundary in the map information of the research city and the commuting population distribution information so as to obtain the final central urban area range of the urban area.
The network connectivity index calculating module 204 specifically includes:
a network connectivity index calculation unit for calculating the network connectivity index according to the commuting OD data and the formula
Figure GDA0004015433540000171
Calculating a network connectivity index between any two administrative areas in the research area; wherein R is ij The network connectivity of the ith administrative area and the jth administrative interval;
Figure GDA0004015433540000172
Figure GDA0004015433540000173
in order to standardize the commuting connection strength index of the employment from the ith administrative district to the jth administrative district,
Figure GDA0004015433540000174
the normalized commute contact strength index of employment from the j administrative district to the i administrative district;
Figure GDA0004015433540000181
V ij the number of commuters who live in the jth administrative district from the ith administrative district,
Figure GDA0004015433540000182
the sum of commuting population from the i th administrative district to employment of each administrative district, wherein n represents the total number of the administrative districts;
Figure GDA0004015433540000183
V ji the number of commutes from the j administrative district to the i administrative district,
Figure GDA0004015433540000184
the commute population is the sum of employment populations from the jth administrative area to each administrative area; max (R) ij ) The maximum value of the network connectivity in the commuting network of all administrative areas is obtained;
Figure GDA0004015433540000185
the network connectivity indexes of the ith administrative area and the jth administrative interval are standardized.
From the perspective of administrative management, resource management and government governance in China, there is no management level matched with an urban circle, so that the development of regions with the central city as a core for overall urbanization development to a certain degree is required. The objective definition of the spatial range of the urban area is particularly important for whether the resultant force and the collaborative development can be formed between cities in the urban area.
The evolution trend of urban circle multi-center and the formation of peripheral urban circle secondary center inevitably promote the networked spatial pattern, and urban circle space definition needs to consider the commuting connection from the first-order abdominal area (central urban commuting area) to the second-order abdominal area (peripheral radiation expansion area). The network connectivity index is more accurate than the centripetal commuting rate index in describing the closeness of the commuting relation among urban areas, and the first-order comprehensive index of the network connectivity of the abdominal region and the second-order comprehensive index of the network connectivity of the abdominal region are both more than 1, so that the network connectivity index can be used as the standard for defining the urban circle space of the Chinese city. In addition, with the urbanization development of China, the urban area range also has dynamic change characteristics, the urban area space range is divided and brought into the homeland space planning and urban physical examination work, and continuous tracking and dynamic monitoring have important significance
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (8)

1. A metropolitan area range determination method based on commute big data is characterized by comprising the following steps:
determining commuting OD data of a research city; the commuting OD data comprise resident residence information and resident employment place information of which at least one end is positioned in a research city;
determining a research area according to the commuting OD data, and dividing the research area according to map information of the research area to obtain an administrative area set contained in a metro ring potential space range; the administrative region set comprises a plurality of administrative regions; the research area is an administrative area covered by a resident residence or a resident employment place in the commuting OD data;
determining the central urban area of the research city as the central urban area range of the urban circle;
calculating a network connectivity index between any two administrative areas in the research area according to the commuting OD data;
calculating a first-order comprehensive index of the network connectivity of the abdominal region and the urban region of each first set administrative region, and determining a first-order abdominal region range of the urban circle according to the first-order comprehensive index of the network connectivity of the abdominal region and the urban region; the first set administrative region is each administrative region located in the research region range; the first-order abdominal region and urban region network connectivity comprehensive index is the sum of network connectivity indexes of each administrative region contained in the first set administrative region and the central urban region range;
determining closely-connected administrative regions of each second set administrative region; the tight connection administrative areas are administrative areas which are associated with the second set administrative area and have network connectivity indexes larger than a tight connection set threshold value in the research area; the second set administrative region is an administrative region which is in the research region and is not positioned in the central urban region range and the first-order abdominal region range;
calculating a second-order abdominal region and urban area network connectivity comprehensive index corresponding to each second set administrative region, and judging whether each second set administrative region belongs to a second-order abdominal standard or not according to the second-order abdominal region and urban area network connectivity comprehensive index so as to obtain a second-order abdominal range of the urban area; the second-order comprehensive index of the network connectivity of the abdominal region is the sum of the second set administrative region and the first-order comprehensive index of the network connectivity of the abdominal region of each administrative region closely connected with the second set administrative region;
converging the central urban area range, the first-order abdominal region range and the second-order abdominal region range to obtain an urban circle space range;
the calculating, according to the commuting OD data, a network connectivity index between any two administrative areas in the research area specifically includes:
according to the commute OD data and formula
Figure FDA0004015433530000021
Calculating a network connectivity index between any two administrative areas in the research area;
wherein R is ij The network connectivity of the ith administrative area and the jth administrative interval;
Figure FDA0004015433530000022
Figure FDA0004015433530000023
in order to standardize the commuting contact strength index of employment from the ith administrative district to the jth administrative district,
Figure FDA0004015433530000024
the normalized commute contact strength index of employment from the j administrative district to the i administrative district;
Figure FDA0004015433530000025
V ij the number of commutes from the i th administrative district to the j th administrative district,
Figure FDA0004015433530000026
the sum of commuting population from the i th administrative district to employment of each administrative district, wherein n represents the total number of the administrative districts;
Figure FDA0004015433530000031
V ji the number of commuters who live from the jth administrative district to the ith administrative district,
Figure FDA0004015433530000032
the commute population sum from the j administrative district to the employment of each administrative district; max (R) ij ) The maximum value of the network connectivity in the commuting network of all administrative areas is obtained;
Figure FDA0004015433530000033
the network connectivity indexes of the ith administrative area and the jth administrative interval are normalized.
2. The method according to claim 1, wherein the research area is divided according to the map information of the research area to obtain an administrative area set included in a metropolitan area potential space range, and the method specifically comprises:
and dividing the research area according to the district-county administrative boundary in the map information of the research area to obtain an administrative area set contained in the urban circle potential space range.
3. The method of claim 1, wherein the determining the central urban area of the study city as the central urban area range of the metropolitan circle comprises:
determining the central area of the research city as a preliminary central urban area of the urban circle according to the total planning range of the current city and the planning and regulated range of the territorial space;
determining commuting demographic distribution information for the research city;
and according to the district and county administrative boundaries in the map information of the research city and the commuting population distribution information, locally adjusting the initial central urban area to obtain the final central urban area range of the urban circle.
4. The method as claimed in claim 1, wherein the determining the first-order abdominal range of the metropolitan area based on the commute big data comprises:
sequentially judging whether the comprehensive index of the network connectivity of the first-order abdominal region urban area is greater than a first-order abdominal region judgment threshold value or not;
and if so, determining the first-order abdominal land as the first-order abdominal land according to the first-order abdominal land urban area network connectivity comprehensive index, and further determining the first-order abdominal land range of the urban circle.
5. The method according to claim 1, wherein the determining whether each second set administrative area belongs to a second-order abdominal standard is performed according to the second-order comprehensive index of network connectivity in the abdominal region and the urban area, so as to obtain a second-order abdominal range of the urban circle, specifically comprises:
sequentially judging whether the second-order abdominal region urban network connectivity comprehensive index is greater than a second-order abdominal region judgment threshold value;
and if so, determining a second set administrative district corresponding to the second-order abdominal region urban network connectivity comprehensive index as a second-order abdominal region, and further determining a second-order abdominal region range of the urban district.
6. An urban circle range determination system based on commute big data, comprising:
the data acquisition module is used for determining commuting OD data of a research city; the commuting OD data comprise resident residence information and resident employment place information of which at least one end is positioned in a research city;
the administrative region set determining module is used for determining a research region according to the commuting OD data, dividing the research region according to map information of the research region and obtaining an administrative region set contained in a potential spatial range of the urban circle; the administrative region set comprises a plurality of administrative regions; the study area is an administrative area covered by a resident residence or a resident employment place in the commuting OD data;
the central urban area range determining module is used for determining the central urban area of the research city as the central urban area range of the urban circle;
the network connectivity index calculation module is used for calculating the network connectivity index between any two administrative areas in the research area according to the commuting OD data;
the first-order abdominal region range determining module is used for calculating a first-order abdominal region and urban area network connectivity comprehensive index of each first set administrative region and determining a first-order abdominal region range of the urban area according to the first-order abdominal region and urban area network connectivity comprehensive index; the first set administrative region is each administrative region located in the research region range; the first-order abdominal region and urban region network connectivity comprehensive index is the sum of the network connectivity indexes of the first set administrative region and the administrative regions contained in the central urban region range;
the close connection administrative region determining module is used for determining the close connection administrative region of each second set administrative region; the tight connection administrative areas are administrative areas which are associated with the second set administrative area and have network connectivity indexes larger than a tight connection set threshold value in the research area; the second set administrative region is an administrative region which is in the research region and is not positioned in the central urban region range and the first-order abdominal region range;
the second-order abdominal region determining module is used for calculating a second-order abdominal region and urban area network connectivity comprehensive index corresponding to each second set administrative region, and judging whether each second set administrative region belongs to a second-order abdominal region standard or not according to the second-order abdominal region and urban area network connectivity comprehensive index so as to obtain a second-order abdominal region range of the urban area; the second-order comprehensive index of the network connectivity of the abdominal region and the urban area is the sum of the second set administrative area and the first-order comprehensive index of the network connectivity of the abdominal region and the urban area of each administrative area tightly connected with the second set administrative area;
the urban circle space range determining module is used for converging the central urban area range, the first-order abdominal region range and the second-order abdominal region range to obtain an urban circle space range;
the network connectivity index calculation module specifically includes:
a network connectivity index calculation unit for calculating the network connectivity index according to the commuting OD data and the formula
Figure FDA0004015433530000061
Calculating a network connectivity index between any two administrative areas in the research area;
wherein R is ij The network connectivity of the ith administrative area and the jth administrative interval;
Figure FDA0004015433530000062
Figure FDA0004015433530000063
in order to standardize the commuting contact strength index of employment from the ith administrative district to the jth administrative district,
Figure FDA0004015433530000064
inhabit the jth administrative district to the jth administrative district after standardizationThe commuting connection strength index of employment of the i administrative areas;
Figure FDA0004015433530000065
V ij the number of commutes from the i th administrative district to the j th administrative district,
Figure FDA0004015433530000066
the commute population is the sum of employment populations from the ith administrative area to the administrative areas, and n represents the total number of the administrative areas;
Figure FDA0004015433530000067
V ji the number of commutes from the j administrative district to the i administrative district,
Figure FDA0004015433530000068
the commute population sum from the j administrative district to the employment of each administrative district; max (R) ij ) The maximum value of the network connectivity in the commuting network of all administrative areas is obtained;
Figure FDA0004015433530000069
the network connectivity indexes of the ith administrative area and the jth administrative interval are standardized.
7. The commute big data-based metro area range determination system according to claim 6, wherein said administrative region set determination module specifically comprises:
and the administrative region set determining unit is used for dividing the research region according to the district-county administrative boundary in the map information of the research region to obtain an administrative region set contained in the urban area potential space range.
8. The metropolitan area coverage determination system based on big commute data of claim 6, wherein the central metropolitan area coverage determination module specifically comprises:
a preliminary central urban area determining unit, configured to determine a central area of the study city as a preliminary central urban area of an urban circle according to a current urban overall planning range and a national space planning regulation range;
a commute population distribution information determining unit, configured to determine commute population distribution information of the research city;
and the central urban area range determining unit is used for locally adjusting the initial central urban area according to the district-county administrative boundary in the map information of the research city and the commuting population distribution information so as to obtain the final central urban area range of the urban area.
CN202110260787.5A 2021-03-10 2021-03-10 City circle range determining method and system based on commuting big data Active CN112948769B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110260787.5A CN112948769B (en) 2021-03-10 2021-03-10 City circle range determining method and system based on commuting big data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110260787.5A CN112948769B (en) 2021-03-10 2021-03-10 City circle range determining method and system based on commuting big data

Publications (2)

Publication Number Publication Date
CN112948769A CN112948769A (en) 2021-06-11
CN112948769B true CN112948769B (en) 2023-01-31

Family

ID=76229185

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110260787.5A Active CN112948769B (en) 2021-03-10 2021-03-10 City circle range determining method and system based on commuting big data

Country Status (1)

Country Link
CN (1) CN112948769B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116233823B (en) * 2023-05-10 2023-08-01 深圳市城市交通规划设计研究中心股份有限公司 Identification method of cross-city commute ring, electronic equipment and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017133110A1 (en) * 2016-02-02 2017-08-10 东南大学 Urban dynamic spatial structure circle layer definition method
CN110533038A (en) * 2019-09-04 2019-12-03 广州市交通规划研究院 A method of urban vitality area and inner city Boundary Recognition based on information data

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017133110A1 (en) * 2016-02-02 2017-08-10 东南大学 Urban dynamic spatial structure circle layer definition method
CN110533038A (en) * 2019-09-04 2019-12-03 广州市交通规划研究院 A method of urban vitality area and inner city Boundary Recognition based on information data

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于手机信令数据的城市群地区都市圈空间范围多维识别——以京津冀为例;赵鹏军等;《城市发展研究》;20190926;第26卷(第9期);全文 *
基于通勤大数据的城市职住分离度研究――《2020年全国主要城市通勤监测报告》核心指标分析;冉江宇等;《城市交通》;20200925;第18卷(第05期);全文 *
都市圈多维界定及其空间匹配关系研究――以京津冀地区为例;姚永玲等;《城市发展研究》;20200726;第27卷(第07期);全文 *

Also Published As

Publication number Publication date
CN112948769A (en) 2021-06-11

Similar Documents

Publication Publication Date Title
CN109299438B (en) Public transport facility supply level evaluation method based on network appointment data
WO2017133110A1 (en) Urban dynamic spatial structure circle layer definition method
CN111581315B (en) Public service facility reachability calculation method and device
Zheng et al. Exploring both home-based and work-based jobs-housing balance by distance decay effect
CN102117383A (en) Method for diagnosing river pollution in real time
CN108495254B (en) Traffic cell population characteristic estimation method based on signaling data
CN102033932A (en) Integrated design-oriented urban rail transit station sorting method
Bertsiou et al. Multi‐criteria analysis and Geographic Information Systems methods for wind turbine siting in a North Aegean Island
CN112348344A (en) Public transport reachable index calculation method
Lv et al. Polycentric urban development and its determinants in China: A geospatial big data perspective
CN112365391A (en) Land diversity measurement method based on 'homeland survey' data
CN112948769B (en) City circle range determining method and system based on commuting big data
CN107194541A (en) A kind of power distribution network power supply zone method based on adaptive weighting Voronoi diagram
CN113222327A (en) Method for analyzing block accessibility based on service capability
CN116796904A (en) Method, system, electronic equipment and medium for predicting new line passenger flow of rail transit
CN115049215A (en) Method for constructing urban traffic flow spatial distribution model
CN112533140B (en) Shared bicycle distribution condition evaluation method based on index
CN114037239A (en) Potential model employment reachability analysis method based on multi-source big data
Xing et al. Sustainable development evaluation of urban traffic system
CN104700616B (en) Urban traffic pressure Forecasting Methodology and system
CN104794164A (en) Method for recognizing settlement parking spaces meeting social parking requirement on basis of open source data
CN113850430A (en) Bus trunk network identification method based on bus card swiping data
CN113744525A (en) Traffic distribution prediction method based on feature extraction and deep learning
CN113032693A (en) 15-minute life circle dividing method considering medical facility service capacity
Sun et al. A Novel Model for Accessibility of Bus Station and Its Application to Bottlenecks Identification: A Case Study in Harbin

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