Movatterモバイル変換


[0]ホーム

URL:


CN113487465A - City overlapping structure characteristic detection method and system based on label propagation algorithm - Google Patents

City overlapping structure characteristic detection method and system based on label propagation algorithm
Download PDF

Info

Publication number
CN113487465A
CN113487465ACN202110691291.3ACN202110691291ACN113487465ACN 113487465 ACN113487465 ACN 113487465ACN 202110691291 ACN202110691291 ACN 202110691291ACN 113487465 ACN113487465 ACN 113487465A
Authority
CN
China
Prior art keywords
urban
layer
data
overlapping
city
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110691291.3A
Other languages
Chinese (zh)
Other versions
CN113487465B (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 University of Geosciences Wuhan
Original Assignee
China University of Geosciences Wuhan
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 University of Geosciences WuhanfiledCriticalChina University of Geosciences Wuhan
Priority to CN202110691291.3ApriorityCriticalpatent/CN113487465B/en
Publication of CN113487465ApublicationCriticalpatent/CN113487465A/en
Application grantedgrantedCritical
Publication of CN113487465BpublicationCriticalpatent/CN113487465B/en
Expired - Fee Relatedlegal-statusCriticalCurrent
Anticipated expirationlegal-statusCritical

Links

Images

Classifications

Landscapes

Abstract

Translated fromChinese

本发明涉及城市规划领域,提供一种基于标签传播算法的城市重叠结构特征检测方法及系统,包括:获得城市网格和处理后的轨迹数据;将城市网格和处理后的轨迹数据进行加权匹配,获得四层有向加权网络;将四层有向加权网络输入图划分模型进行无监督训练,获得四层城市社区结构,提取四层城市社区结构中的重叠结构;通过兴趣点数据构建测量指标,通过测量指标对重叠结构进行识别,获得重叠结构的土地使用特点和空间交互模式。本发明将网络科学中图划分的方法引入城市规划,具有较好的效益,同时能够批量化、自动化的进行城市结构划分;并且能充分挖掘隐藏在城市居民活动中城市的空间交互信息,挖掘其土地利用特征和其空间交互关系。

Figure 202110691291

The invention relates to the field of urban planning, and provides a method and system for detecting overlapping structural features of cities based on a label propagation algorithm, including: obtaining urban grids and processed trajectory data; weighting and matching the urban grids and the processed trajectory data , obtain a four-layer directed weighted network; perform unsupervised training on the input graph partition model of the four-layer directed weighted network to obtain a four-layer urban community structure, and extract the overlapping structure in the four-layer urban community structure; construct measurement indicators through interest point data , the overlapping structures are identified through measurement indicators, and the land use characteristics and spatial interaction patterns of the overlapping structures are obtained. The invention introduces the method of graph division in network science into urban planning, which has good benefits, and at the same time, it can divide the urban structure in batches and automatically; Land use characteristics and their spatial interactions.

Figure 202110691291

Description

City overlapping structure characteristic detection method and system based on label propagation algorithm
Technical Field
The invention relates to the field of urban planning, in particular to a method and a system for detecting urban overlapping structure characteristics based on a label propagation algorithm.
Background
In the past thirty years, sufficient labor, good infrastructure and cheap land brought by the urbanization process lay a foundation for rapid development of economy. But can not avoid a plurality of problems in the urbanization process of China. The urban problem is particularly serious for some provincial cities or metropolis. The urban diseases are mainly manifested by traffic jam, housing shortage, water supply shortage, energy shortage, environmental deterioration and the like, which cause burden to cities and even restrict the development of the cities. The development structure of the city is closely related to the life and economy of urban residents, the urban spatial structure is determined by combining human activities through scientific means, an operable, scientific and reasonable spatial structure analysis method is provided, and the method becomes an important direction of digital urban research.
The urban structure is gradually complicated and diversified, obvious hierarchy and overlapping are achieved, the urban areas with different hierarchies have obvious hierarchical overlapping relations, the hierarchical overlapping relations are discussed from human activities, the regional change and the spatial distribution of the urban spatial structure can be gradually mastered from local parts to the whole, the interaction between the urban overlapping structure and other plots is far larger than the interaction between the urban overlapping structure and other plots, the urban overlapping structure can be understood as a junction area of urban spatial interaction, and the interaction can be calculated through human behaviors. Some experts have made relevant researches on city structure division methods, which can be divided into statistical survey-based methods and model-based methods, but currently, researches on city structure hierarchy and city overlap are still few. The method based on statistical survey combines survey statistics and expert evaluation modes for carrying out the demarcation, namely, in the process of demarcation of the urban structure, on the basis of on-site survey statistical results, a plurality of experts which have certain cognition on the city and have higher representativeness and authority are selected for carrying out the evaluation. This method is generally associated with great subjectivity and high time, labor and capital costs. The method based on the model is used for defining urban areas through scientific data analysis and big data mining methods under the support of public-source geographic big data, and provides an operable, scientific and reasonable space optimization model. The multi-source geographic data has the advantages of large data volume, strong current, rich sources, low cost and the like. Based on the characteristic of collecting the multi-source geographic data from bottom to top, researchers can easily obtain the spatial-temporal information which is wide in city range, abundant in mass and based on individuals, so that fine geographic analysis and modeling are achieved, and better service is provided for researching city structures.
The above is only for the purpose of assisting understanding of the technical aspects of the present invention, and does not represent an admission that the above is prior art.
Disclosure of Invention
The invention mainly aims to solve the technical problems of high subjectivity, high time, labor and capital cost in the prior art, and can detect the level overlapping structure in a city and identify the characteristics of the level overlapping structure.
In order to achieve the above object, the present invention provides a method for detecting urban overlapping structure features based on a label propagation algorithm, comprising the steps of:
s1: obtaining city map data, taxi track data and interest point data in a research area, and carrying out city unit division on the city map data to obtain city grids; preprocessing the taxi track data to obtain processed track data;
s2: carrying out weighted matching on the urban grids and the processed trajectory data to obtain a four-layer directed weighting network;
s3: carrying out unsupervised training on the four-layer directed weighting network input graph division model, obtaining four-layer urban community structures after the training is finished, and extracting an overlapping structure in the four-layer urban community structures;
s4: and constructing a measurement index through the point of interest data, and identifying the overlapped structure through the measurement index to obtain the land use characteristics and the space interaction mode of the overlapped structure.
Preferably, step S1 is specifically:
s11: processing the urban map data through GIS software, performing spatial fishing net analysis on urban areas of the urban map data, and dividing the urban areas into urban grids; the city grid comprises a plurality of 500mx500m city grid cells;
s12: removing the point data and invalid point data which are not in the urban area in the taxi track data to obtain removed track data;
s13: and extracting the data of the point of getting on or off the vehicle in each piece of the eliminated track data, wherein the processed track data is a set of the data of the point of getting on or off the vehicle.
Preferably, step S2 is specifically:
s21: matching each getting-on and getting-off point data in the processed track data with each city grid unit, simulating each city grid unit into a graph node, and simulating the number of interaction times among the city grid units into the weight of an edge;
s22: the directed weighting network is composed of a plurality of city grid units and interactive relations among the city grid units, and the interactive relations are related to the interactive times among the city grid units;
s23: and path layering is carried out on the processed track data, a first layer is used when the track path is less than 3km, a second layer is used when the track path is less than 5km, a third layer is used when the track path is less than 9km, all the track paths are fourth layers, corresponding directed weighting networks are respectively constructed on the first layer, the second layer, the third layer and the fourth layer, and the four-layer directed weighting networks are obtained.
Preferably, step S3 is specifically:
s31: initializing the memory of each node in the graph partitioning model by using the id of the corresponding node, and obtaining a corresponding unique label by each node;
s32: selecting a certain node as a listener node;
s33: all adjacent nodes of the listener node send own unique labels to the listener node, and the listener node selects the most popular label from all the received labels;
s34: repeating the steps S32-S33 for n times, traversing all the nodes and obtaining the most popular labels of all the nodes;
s35: post-processing all the labels of the nodes to obtain the four-layer urban community structure, and evaluating the division result of the four-layer urban community structure through an overlapping modularity function, wherein the overlapping modularity function specifically comprises the following steps:
Figure BDA0003126271390000031
wherein m is the weight sum of edges in the network, A is the weighted adjacent matrix of the network, and if an edge exists between the node v and the node w, A isvwThe weight of vw edge is, otherwise, 0 is adopted; k is a radical ofv,kwRespectively, the out-degree weight of the node v and the in-degree weight of the node w, Ov,OwThe number of communities to which the node v and the node w belong respectively;
s36: and extracting an overlapping structure in the four-layer city community structure.
Preferably, in step S4, the measurement index includes: richness, simpson index and entropy measure index;
the land use condition and the function type can be identified through the richness; the land mixing condition can be identified through the Simpson index and the entropy measurement index;
the formula of the richness is specifically as follows:
Figure BDA0003126271390000032
Fi,ldenotes the enrichment index, n, of POIs of class I in the ith plotl,iRepresenting the number of type I land use types in the ith plot, niIs the number of all POIs in the ith parcel. N is a radical oflN is the total number of POIs of class i, and N is the total number of POIs in the whole research area;
the simpson index and the entropy measurement index are expressed by a hill index, and the formula specifically comprises:
Figure BDA0003126271390000033
in the formula, D represents the value of the Hill index, piRepresenting the proportion of the i-th POI; when q is 1, it represents entropy, and higher values indicate that the distribution of POI species is more disordered, and lower values indicate that the distribution of POI species is more ordered; q is 2, the inverse of the happson index, which measures the probability that two POIs randomly selected from an urban area belong to the same category; therefore, the abundance of the POI and the relative abundance of different types of POI are considered, and the lower the value is, the higher the mixed utilization degree of the land is, and the higher the value is, the lower the mixed utilization degree of the land is.
Preferably, in step S4;
the land use characteristics are land use types and land mixing degrees, the functional areas are structures which show certain functional attributes in urban areas, such as residential areas and scenic areas, and the functional structures and the function mixing degrees of the overlapped structures can be measured through the richness; the POI mixing degree of the overlapped structure can be calculated through the Simpson index and the entropy measurement index, and the POI mixing degree can reflect the vitality of the land;
the space interaction mode represents the interaction condition of the overlapping area and the adjacent community, the interaction strength is expressed through the interaction times, an interaction network is built, the travel mode of people is analyzed through the trajectory flow from one functional area to another functional area, the travel mode has a fixed time law such as the interaction mode from a living area to a working area in the morning rush hour and the interaction mode from the living area to a leisure area in the holiday, the interaction hotspot area of the overlapping area is analyzed through the interaction strength, and the functional structure of the interaction area is further identified through the richness.
A city overlapping structure feature detection system based on a label propagation algorithm comprises the following modules:
the data acquisition module is used for acquiring city map data, taxi track data and interest point data in a research area, and performing city unit division on the city map data to obtain city grids; preprocessing the taxi track data to obtain processed track data;
the four-layer directed weighting network acquisition module is used for performing weighted matching on the urban grids and the processed trajectory data to acquire a four-layer directed weighting network;
the overlapping structure extraction module is used for carrying out unsupervised training on the four-layer directed weighting network input graph division model, obtaining four-layer urban community structures after the training is finished, and extracting an overlapping structure in the four-layer urban community structures;
and the identification module is used for constructing a measurement index through the point of interest data, identifying the overlapped structure through the measurement index and obtaining the land use characteristics and the space interaction mode of the overlapped structure.
The invention has the following beneficial effects:
1. the invention adopts an analogy reasoning method, introduces a method for dividing the network science into city planning, has better benefit, and can divide the city structure in batch and automatically;
2. the invention can fully mine the spatial interaction information of the city hidden in the urban resident activity, pay attention to the overlapping structure existing in the network, and mine the land utilization characteristics and the spatial interaction relationship thereof.
Drawings
FIG. 1 is a flow chart of a method according to an embodiment of the present invention;
FIG. 2 is a system block diagram according to an embodiment of the present invention;
the implementation, functional features and advantages of the objects of the present invention will be further explained with reference to the accompanying drawings.
Detailed Description
It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Referring to fig. 1, the invention provides a method for detecting characteristics of an urban overlapping structure based on a label propagation algorithm, which is an extension method based on a model method, and applies the idea of network science to urban division on the basis of the previous space division model research, and detects the overlapping structure existing in the city; the method is based on the idea of label propagation in graph division, adopts an analogy reasoning method, aims at the hierarchy and the overlapping property of cities in the urbanization process, discusses the hierarchy overlapping property from human activities, and is beneficial to gradually grasp the regional change and the spatial distribution of the city space structure from local to whole; the method effectively introduces the concept of overlapping nodes in the complex network into city division, combines taxi track data, performs community detection on long and short distance city community structures, fully excavates the space interaction information of cities hidden in city resident activities, pays attention to the overlapping structures in the city community structures, and excavates the land utilization type and the space interaction relationship thereof;
the method specifically comprises the following steps:
s1: obtaining city map data, taxi track data and interest point data in a research area, and carrying out city unit division on the city map data to obtain city grids; preprocessing the taxi track data to obtain processed track data;
s2: carrying out weighted matching on the urban grids and the processed trajectory data to obtain a four-layer directed weighting network;
s3: carrying out unsupervised training on the four-layer directed weighting network input graph division model, obtaining four-layer urban community structures after the training is finished, and extracting an overlapping structure in the four-layer urban community structures;
s4: and constructing a measurement index through the point of interest data, and identifying the overlapped structure through the measurement index to obtain the land use characteristics and the space interaction mode of the overlapped structure.
In this embodiment, step S1 specifically includes:
s11: processing the urban map data through GIS software, performing spatial fishing net analysis on urban areas of the urban map data, and dividing the urban areas into urban grids; the city grid comprises a plurality of 500mx500m city grid cells; in this embodiment, 4853 city grid cells are obtained in total;
s12: removing the point data and invalid point data which are not in the urban area in the taxi track data to obtain removed track data;
s13: extracting the boarding and alighting point data in each eliminated track data, wherein the processed track data is a set of the boarding and alighting point data; in this embodiment, 793253 pieces of boarding and alighting point data are obtained in total.
In this embodiment, step S2 specifically includes:
s21: matching each getting-on and getting-off point data in the processed track data with each city grid unit, simulating each city grid unit into a graph node, and simulating the number of interaction times among the city grid units into the weight of an edge;
s22: the directed weighting network is composed of a plurality of city grid units and interactive relations among the city grid units, and the interactive relations are related to the interactive times among the city grid units; each city grid unit interacts with a plurality of other city grid units;
s23: path layering is carried out on the processed track data, a first layer is used when the track path is less than 3km, a second layer is used when the track path is less than 5km, a third layer is used when the track path is less than 9km, all the track paths are fourth layers, corresponding directed weighting networks are respectively constructed on the first layer, the second layer, the third layer and the fourth layer, and the four layers of directed weighting networks are obtained;
in the specific implementation, the route layering threshold is determined according to the track proportion of each layer of route, and the tracks less than 3km, less than 5km and less than 9km respectively account for 29.3301%, 51.1679% and 76.0485% of the total track; meanwhile, the network structure of the four layers of directional weighting networks can be changed according to requirements.
In this embodiment, step S3 specifically includes:
s31: initializing the memory of each node in the graph partitioning model by using the id of the corresponding node, and obtaining a corresponding unique label by each node;
s32: selecting a certain node as a listener node;
s33: all adjacent nodes of the listener node send own unique labels to the listener node, and the listener node selects the most popular label from all the received labels;
s34: repeating the steps S32-S33 for n times, traversing all the nodes and obtaining the most popular labels of all the nodes;
in the specific implementation, when a four-layer directed weighting network input graph division model is subjected to unsupervised training, an SLPA model parameter needs to be set, in the graph division model, the iteration number of the model is set to be 100, the number of communities divided into the smallest layers is set to be 3, in order to keep the result of each division consistent, the random seed is set to be 5140727168289296997, meanwhile, through contrast detection, a random seed value is used, the modularity fluctuation result of the community division result is less than 0.1, and the normal modularity value is 0.3-0.7. In addition, r parameters for controlling the output of the overlapped communities are determined through comparison, wherein r belongs to {0.01,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45 and 0.5}, the modularity curve inflection point is determined when r is 0.1 through comparison experiments, and finally, r is selected to be 0.1 through training parameters, and the parameter selection can be changed according to requirements;
s35: post-processing all the labels of the nodes to obtain the four-layer urban community structure, and evaluating the division result of the four-layer urban community structure through an overlapping modularity function, wherein the overlapping modularity function specifically comprises the following steps:
Figure BDA0003126271390000071
wherein m is the weight sum of edges in the network, A is the weighted adjacent matrix of the network, and if an edge exists between the node v and the node w, A isvwThe weight of vw edge is, otherwise, 0 is adopted; k is a radical ofv,kwRespectively, the out-degree weight of the node v and the in-degree weight of the node w, Ov,OwThe number of communities to which the node v and the node w belong respectively;
in the specific implementation, the overlapping modularity function results of the four-layer urban community structure fluctuate up and down at 0.65, 0.6, 0.55 and 0.3 respectively, and considering the reasons of overlapping nodes and path lengths, it can be explained that the more the overlapping nodes are identified, the more chaotic the area is, the area possibly belongs to a plurality of areas, the result of community division is not good, the longer the path is, the larger the scale of community division is, and the interaction of short paths can influence the division result; for the overlapped nodes, the first community is the community with the maximum label probability, and so on;
s36: and extracting an overlapping structure in the four-layer city community structure.
In this embodiment, in step S4, the point-of-interest data is obtained through the gold API, and the data is reclassified; adopting a crawler method to acquire POI data of each category, wherein the data comprises 622206 data, and the POI data comprises attributes such as POI name, longitude and latitude, large category, medium category, small category, address and the like; meanwhile, referring to urban land use categories, reclassifying POIs into 16 categories;
the measurement index includes: richness, simpson index and entropy measure index;
the land use condition and the function type can be identified through the richness; the land mixing condition can be identified through the Simpson index and the entropy measurement index;
the formula of the richness is specifically as follows:
Figure BDA0003126271390000072
Fi,ldenotes the enrichment index, n, of POIs of class I in the ith plotl,iRepresenting the number of type I land use types in the ith plot, niIs the number of all POIs in the ith parcel. N is a radical oflN is the total number of POIs of class i, and N is the total number of POIs in the whole research area;
the simpson index and the entropy measurement index are expressed by a hill index, and the formula specifically comprises:
Figure BDA0003126271390000081
formula (II)Wherein D represents the value of the Hill index, piRepresenting the proportion of the i-th POI; when q is 1, it represents entropy, and higher values indicate that the distribution of POI species is more disordered, and lower values indicate that the distribution of POI species is more ordered; q is 2, the inverse of the happson index, which measures the probability that two POIs randomly selected from an urban area belong to the same category; therefore, the abundance of the POI and the relative abundance of different types of POI are considered, and the lower the value is, the higher the mixed utilization degree of the land is, and the higher the value is, the lower the mixed utilization degree of the land is.
In this embodiment, in step S4;
the land use characteristics are land use types and land mixing degrees, the functional areas are structures which show certain functional attributes in urban areas, such as residential areas and scenic areas, and the functional structures and the function mixing degrees of the overlapped structures can be measured through the richness; the POI mixing degree of the overlapped structure can be calculated through the Simpson index and the entropy measurement index, and the POI mixing degree can reflect the vitality of the land;
the space interaction mode represents the interaction condition of the overlapping area and the adjacent community, the interaction strength is expressed through the interaction times, an interaction network is built, the travel mode of people is analyzed through the trajectory flow from one functional area to another functional area, the travel mode has a fixed time law such as the interaction mode from a living area to a working area in the morning rush hour and the interaction mode from the living area to a leisure area in the holiday, the interaction hotspot area of the overlapping area is analyzed through the interaction strength, and the functional structure of the interaction area is further identified through the richness.
Referring to fig. 2, the present invention provides a system for detecting characteristics of an urban overlapping structure based on a label propagation algorithm, including the following modules:
thedata acquisition module 10 is configured to acquire city map data, taxi track data, and interest point data in a research area, and perform city unit division on the city map data to obtain city grids; preprocessing the taxi track data to obtain processed track data;
a four-layer directional weightingnetwork obtaining module 20, configured to perform weighted matching on the city grid and the processed trajectory data to obtain a four-layer directional weighting network;
the overlapstructure extraction module 30 is configured to perform unsupervised training on the four-layer directed weighting network input graph partitioning model, obtain a four-layer urban community structure after the training is completed, and extract an overlap structure in the four-layer urban community structure;
and theidentification module 40 is configured to construct a measurement index according to the point-of-interest data, identify the overlapping structure according to the measurement index, and obtain land use characteristics and a space interaction mode of the overlapping structure.
It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or system. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or system that comprises the element.
The above-mentioned serial numbers of the embodiments of the present invention are merely for description and do not represent the merits of the embodiments. In the unit claims enumerating several means, several of these means may be embodied by one and the same item of hardware. The use of the words first, second, third and the like do not denote any order, but rather the words first, second and the like may be interpreted as indicating any order.
The above description is only a preferred embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by using the contents of the present specification and the accompanying drawings, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.

Claims (7)

Translated fromChinese
1.一种基于标签传播算法的城市重叠结构特征检测方法,其特征在于,包括步骤:1. an urban overlapping structure feature detection method based on label propagation algorithm, is characterized in that, comprises the steps:S1:获取研究区域内的城市地图数据、出租车轨迹数据和兴趣点数据,对所述城市地图数据进行城市单元划分,获得城市网格;对所述出租车轨迹数据进行预处理,获得处理后的轨迹数据;S1: Acquire city map data, taxi trajectory data, and point-of-interest data in the study area, divide the city map data into city units, and obtain an urban grid; preprocess the taxi trajectory data, and obtain the processed trajectory data;S2:将所述城市网格和所述处理后的轨迹数据进行加权匹配,获得四层有向加权网络;S2: Perform weighted matching on the urban grid and the processed trajectory data to obtain a four-layer directed weighted network;S3:将所述四层有向加权网络输入图划分模型进行无监督训练,训练完成后获得四层城市社区结构,提取所述四层城市社区结构中的重叠结构;S3: Perform unsupervised training on the input graph partition model of the four-layer directed weighted network, obtain a four-layer urban community structure after the training is completed, and extract the overlapping structure in the four-layer urban community structure;S4:通过所述兴趣点数据构建测量指标,通过所述测量指标对所述重叠结构进行识别,获得所述重叠结构的土地使用特点和空间交互模式。S4: Constructing a measurement index based on the point of interest data, identifying the overlapping structure through the measurement index, and obtaining land use characteristics and a spatial interaction pattern of the overlapping structure.2.根据权利要求1所述的基于标签传播算法的城市重叠结构特征检测方法,其特征在于,步骤S1具体为:2. the urban overlapping structure feature detection method based on label propagation algorithm according to claim 1, is characterized in that, step S1 is specifically:S11:通过GIS软件处理所述城市地图数据,对所述城市地图数据的城市区域进行空间渔网分析,将所述城市区域划分为所述城市网格;所述城市网格包括多个500mx500m的城市网格单元;S11: Process the city map data through GIS software, perform spatial fishnet analysis on the city area of the city map data, and divide the city area into the city grid; the city grid includes a plurality of 500mx500m cities grid cells;S12:剔除所述出租车轨迹数据中不在所述城市区域的点数据和无效的点数据,获得剔除后的轨迹数据;S12: Eliminate the point data and invalid point data that are not in the urban area in the taxi trajectory data, and obtain the eliminated trajectory data;S13:提取每条所述剔除后的轨迹数据中的上下车点数据,所述处理后的轨迹数据为所述上下车点数据的集合。S13: Extract the pick-up and drop-off point data in each piece of the excluded trajectory data, and the processed trajectory data is a set of the pick-up point data.3.根据权利要求2所述的基于标签传播算法的城市重叠结构特征检测方法,其特征在于,步骤S2具体为:3. the urban overlapping structure feature detection method based on label propagation algorithm according to claim 2, is characterized in that, step S2 is specifically:S21:将所述处理后的轨迹数据中各所述上下车点数据与各所述城市网格单元进行匹配,将各所述城市网格单元类比为图节点,将各城市网格单元间的交互次数类比为边的权重;S21: Match each of the pick-up and drop-off point data in the processed trajectory data with each of the urban grid units, compare each of the urban grid units as a graph node, and compare the difference between the urban grid units The number of interactions is analogous to the weight of the edge;S22:有向加权网络由多个所述城市网格单元和各所述城市网格单元间的交互关系组成,所述交互关系与所述城市网格单元间的交互次数有关;S22: The directed weighted network is composed of a plurality of the urban grid units and the interaction relationship between each of the urban grid units, and the interaction relationship is related to the number of interactions between the urban grid units;S23:将所述处理后的轨迹数据进行路程分层,轨迹路程小于3km的为第一层,轨迹路程小于5km的为第二层,轨迹路程小于9km的第三层,全部轨迹路程为第四层,对所述第一层、所述第二层、所述第三层和所述第四层分别构建对应的有向加权网络,获得所述四层有向加权网络。S23: Perform route layering on the processed trajectory data, the first layer with the trajectory distance less than 3km, the second layer with the trajectory distance less than 5km, the third layer with the trajectory distance less than 9km, and the fourth layer with all the trajectory distances layer, respectively constructing corresponding directed weighted networks for the first layer, the second layer, the third layer and the fourth layer to obtain the four-layer directed weighted network.4.根据权利要求1所述的基于标签传播算法的城市重叠结构特征检测方法,其特征在于,步骤S3具体为:4. the urban overlapping structure feature detection method based on label propagation algorithm according to claim 1, is characterized in that, step S3 is specifically:S31:将所述图划分模型中各节点的内存均用对应节点的id初始化,各节点获得对应的唯一标签;S31: Initialize the memory of each node in the graph partitioning model with the id of the corresponding node, and each node obtains a corresponding unique label;S32:选择某一节点作为监听器节点;S32: select a node as the listener node;S33:所述监听器节点的所有相邻节点均向所述监听器节点发送自己的唯一标签,所述监听器节点在收到的所有标签中选择最流行标签;S33: All adjacent nodes of the listener node send their own unique tags to the listener node, and the listener node selects the most popular tag among all the tags received;S34:重复步骤S32-S33共n次,遍历所有节点,获得所有节点的最流行标签;S34: Repeat steps S32-S33 for a total of n times, traverse all nodes, and obtain the most popular labels of all nodes;S35:对各所述节点的所有标签进行后处理,获得所述四层城市社区结构,通过重叠模块度函数能够评价所述四层城市社区结构的划分结果,所述重叠模块度函数具体为:S35: Perform post-processing on all the labels of each of the nodes to obtain the four-layer urban community structure. The division result of the four-layer urban community structure can be evaluated through an overlapping modularity function. The overlapping modularity function is specifically:
Figure FDA0003126271380000021
Figure FDA0003126271380000021
其中,m为网络中边的权重和,A为网络的带权邻接矩阵,若节点v到节点w之间存在一条边,则Avw为vw边的权重,反之为0;kv,kw分别为节点v的出度权重和和节点w的入度权重和,Ov,Ow分别为节点v和节点w所属的社区数;Among them, m is the weight sum of the edges in the network, A is the weighted adjacency matrix of the network, if there is an edge between node v and node w, then Avw is the weight of the vw edge, otherwise it is 0; kv , kw are the out-degree weight sum of node v and the in-degree weight sum of node w, respectively, Ov , Ow are the number of communities to which node v and node w belong;S36:提取所述四层城市社区结构中的重叠结构。S36: Extract the overlapping structure in the four-layer urban community structure.5.根据权利要求1所述的基于标签传播算法的城市重叠结构特征检测方法,其特征在于,步骤S4中,所述测量指标包括:丰富度、辛普森指数和熵测量指标;5. The method for detecting urban overlapping structure features based on a label propagation algorithm according to claim 1, wherein in step S4, the measurement index comprises: richness, Simpson index and entropy measurement index;通过所述丰富度能够识别土地使用情况和功能类型;通过所述辛普森指数和所述熵测量指标可以识别土地混合状况;Land use and functional type can be identified by the abundance; land mixing can be identified by the Simpson index and the entropy measure;所述丰富度的公式具体为:The formula for the richness is specifically:
Figure FDA0003126271380000022
Figure FDA0003126271380000022
Fi,表示第i个地块中的第l类POI的富集指数,nl,表示第i个地块中的第l类土地利用类型的数量,ni是第i个地块中所有POI的数量。Nl为第l类POI的总数,N是整个研究区域内POI的总数;Fi, represents the enrichment index of the l-th type of POI in the i-th plot, nl, represents the number of the l-th type of land use type in the i-th plot, ni is all the Number of POIs. Nl is the total number of POIs of type l, and N is the total number of POIs in the entire study area;所述辛普森指数和所述熵测量指标通过希尔指数来表示,公式具体为:The Simpson index and the entropy measurement index are represented by the Hill index, and the formula is specifically:
Figure FDA0003126271380000031
Figure FDA0003126271380000031
公式中D代表希尔指数的值,pi代表第i类POI所占比例;当q=1时,它是代表熵,值越高说明POI种类分布越无序,越低代表POI种类分布越有序;q=2时是幸普森指数的逆值,辛普森指数衡量的是从一个城市区域中随机选择的两个POI属于同一类别的概率;因此,它既考虑了POI的丰富度,又考虑了不同类型POI的相对丰度,值越低说明土地的混合利用程度越高,值越高说明土地的混合利用程度越低。In the formula, D represents the value of the Hill index, and pi represents the proportion of the i-th type of POI; when q=1, it represents the entropy. Ordered; q=2 is the inverse of the Simpson index, which measures the probability that two POIs randomly selected from an urban area belong to the same class; thus, it takes into account both the abundance of POIs and the Considering the relative abundance of different types of POIs, the lower the value, the higher the degree of mixed use of land; the higher the value, the lower the degree of mixed use of land.
6.根据权利要求5所述的基于标签传播算法的城市重叠结构特征检测方法,其特征在于,步骤S4中;6. the urban overlapping structure feature detection method based on label propagation algorithm according to claim 5, is characterized in that, in step S4;所述土地使用特点为土地使用的类型及土地混合程度,功能区为城市区域中表现出一定功能属性的结构,如居住区、风景区,通过所述丰富度可以度量所述重叠结构的功能结构及功能混合程度;通过所述辛普森指数和熵测量指标能够计算所述重叠结构的POI混合程度,所述POI混合程度可以反映地块活力;The land use characteristics are the type of land use and the degree of land mixing, and functional areas are structures in urban areas that exhibit certain functional attributes, such as residential areas and scenic areas, and the functional structure of the overlapping structures can be measured through the richness. and functional mixing degree; the POI mixing degree of the overlapping structure can be calculated by the Simpson index and the entropy measurement index, and the POI mixing degree can reflect the vitality of the plot;所述空间交互模式代表所述重叠区域与邻近社区的交互情况,通过交互的次数来表现交互的强弱,构建交互网络,通过一个功能区到另一个功能区的轨迹流来分析人们的出行模式,所述出行模式如早高峰从居住区至工作区,节假日从居住区至休闲区此类具有固定时间规律的交互模式,通过交互强度分析所述重叠区域交互热点区域并进一步使用所述丰富度识别交互区域的功能结构。The spatial interaction pattern represents the interaction between the overlapping area and neighboring communities. The strength of interaction is represented by the number of interactions, an interaction network is constructed, and people’s travel patterns are analyzed through the trajectory flow from one functional area to another. , the travel patterns such as morning rush hour from residential area to work area, holidays from residential area to leisure area and other interactive modes with fixed time rules, analyze the overlapping area interaction hotspot area by interaction strength and further use the richness Identify the functional structure of the interaction area.7.一种基于标签传播算法的城市重叠结构特征检测系统,其特征在于,包括以下模块:7. An urban overlapping structure feature detection system based on a label propagation algorithm, characterized in that it comprises the following modules:数据获取模块,用于获取研究区域内的城市地图数据、出租车轨迹数据和兴趣点数据,对所述城市地图数据进行城市单元划分,获得城市网格;对所述出租车轨迹数据进行预处理,获得处理后的轨迹数据;The data acquisition module is used to acquire city map data, taxi trajectory data and point-of-interest data in the research area, divide the city map data into city units, and obtain city grids; preprocess the taxi trajectory data , obtain the processed trajectory data;四层有向加权网络获取模块,用于将所述城市网格和所述处理后的轨迹数据进行加权匹配,获得四层有向加权网络;A four-layer directed weighted network acquisition module is used to perform weighted matching between the urban grid and the processed trajectory data to obtain a four-layer directed weighted network;重叠结构提取模块,用于将所述四层有向加权网络输入图划分模型进行无监督训练,训练完成后获得四层城市社区结构,提取所述四层城市社区结构中的重叠结构;The overlapping structure extraction module is used to perform unsupervised training on the input graph partition model of the four-layer directed weighted network, obtain a four-layer urban community structure after the training is completed, and extract the overlapping structure in the four-layer urban community structure;识别模块,用于通过所述兴趣点数据构建测量指标,通过所述测量指标对所述重叠结构进行识别,获得所述重叠结构的土地使用特点和空间交互模式。The identification module is configured to construct a measurement index based on the interest point data, identify the overlapping structure through the measurement index, and obtain the land use characteristics and spatial interaction pattern of the overlapping structure.
CN202110691291.3A2021-06-222021-06-22City overlapping structure characteristic detection method and system based on label propagation algorithmExpired - Fee RelatedCN113487465B (en)

Priority Applications (1)

Application NumberPriority DateFiling DateTitle
CN202110691291.3ACN113487465B (en)2021-06-222021-06-22City overlapping structure characteristic detection method and system based on label propagation algorithm

Applications Claiming Priority (1)

Application NumberPriority DateFiling DateTitle
CN202110691291.3ACN113487465B (en)2021-06-222021-06-22City overlapping structure characteristic detection method and system based on label propagation algorithm

Publications (2)

Publication NumberPublication Date
CN113487465Atrue CN113487465A (en)2021-10-08
CN113487465B CN113487465B (en)2022-09-30

Family

ID=77935765

Family Applications (1)

Application NumberTitlePriority DateFiling Date
CN202110691291.3AExpired - Fee RelatedCN113487465B (en)2021-06-222021-06-22City overlapping structure characteristic detection method and system based on label propagation algorithm

Country Status (1)

CountryLink
CN (1)CN113487465B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN118627750A (en)*2024-07-042024-09-10苏州市中遥数字科技有限公司 A multi-dimensional image processing system based on high-resolution remote sensing data

Citations (9)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN103729475A (en)*2014-01-242014-04-16福州大学Multi-label propagation discovery method of overlapping communities in social network
US20140337356A1 (en)*2013-05-082014-11-13Yahoo! Inc.Identifying Communities Within A Social Network Based on Information Propagation Data
CN107784598A (en)*2017-11-212018-03-09山西大学A kind of network community discovery method
CN107818534A (en)*2017-10-312018-03-20武汉大学A kind of mankind's activity network area division methods with space constraint
CN108446862A (en)*2018-03-292018-08-24山东科技大学The three-stage policy algorithm of overlapping community detection in a kind of community network
CN109493119A (en)*2018-10-192019-03-19南京图申图信息科技有限公司A kind of city commercial center identification method and system based on POI data
CN109614458A (en)*2018-12-202019-04-12中国人民解放军战略支援部队信息工程大学 Method and Device for Mining Urban Community Structure Based on Navigation Data
CN110111575A (en)*2019-05-162019-08-09北京航空航天大学A kind of Forecast of Urban Traffic Flow network analysis method based on Complex Networks Theory
CN111698743A (en)*2020-06-092020-09-22嘉兴学院Complex network community identification method fusing node analysis and edge analysis

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
US20140337356A1 (en)*2013-05-082014-11-13Yahoo! Inc.Identifying Communities Within A Social Network Based on Information Propagation Data
CN103729475A (en)*2014-01-242014-04-16福州大学Multi-label propagation discovery method of overlapping communities in social network
CN107818534A (en)*2017-10-312018-03-20武汉大学A kind of mankind's activity network area division methods with space constraint
CN107784598A (en)*2017-11-212018-03-09山西大学A kind of network community discovery method
CN108446862A (en)*2018-03-292018-08-24山东科技大学The three-stage policy algorithm of overlapping community detection in a kind of community network
CN109493119A (en)*2018-10-192019-03-19南京图申图信息科技有限公司A kind of city commercial center identification method and system based on POI data
CN109614458A (en)*2018-12-202019-04-12中国人民解放军战略支援部队信息工程大学 Method and Device for Mining Urban Community Structure Based on Navigation Data
CN110111575A (en)*2019-05-162019-08-09北京航空航天大学A kind of Forecast of Urban Traffic Flow network analysis method based on Complex Networks Theory
CN111698743A (en)*2020-06-092020-09-22嘉兴学院Complex network community identification method fusing node analysis and edge analysis

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PALLA G ETAL.: "Uncovering the overlapping community structure of complex networks in nature and society", 《NATURE》*
吴建 等: "基于图遍历的局部社区发现算法", 《计算机应用研究》*
毕崇武 等: "基于标签语义关联的城市社群发现研究", 《现代情报》*

Cited By (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN118627750A (en)*2024-07-042024-09-10苏州市中遥数字科技有限公司 A multi-dimensional image processing system based on high-resolution remote sensing data

Also Published As

Publication numberPublication date
CN113487465B (en)2022-09-30

Similar Documents

PublicationPublication DateTitle
CN114969007B (en) A method for identifying urban functional areas based on functional mixing degree and ensemble learning
Xu et al.Quantitative analysis of spatial vitality and spatial characteristics of urban underground space (UUS) in metro area
BlanchardMathematical analysis of urban spatial networks
CN113806419B (en)Urban area function recognition model and recognition method based on space-time big data
CN108133302A (en)A kind of public bicycles potential demand Forecasting Methodology based on big data
Karimi et al.Urban expansion modeling using an enhanced decision tree algorithm
Yang et al.How to improve urban transportation planning in big data era? A practice in the study of traffic analysis zone delineation
CN110176141A (en)A kind of traffic zone division method and system based on POI and traffic characteristic
AhmedModelling spatio-temporal urban land cover growth dynamics using remote sensing and GIS techniques: A case study of Khulna City
CN113642757A (en)Internet of things charging pile construction planning method and system based on artificial intelligence
CN114861277A (en)Long-time-sequence national soil space function and structure simulation method
CN112954623B (en)Resident occupancy rate estimation method based on mobile phone signaling big data
CN115965171A (en)Micro-park site selection method based on ant colony optimization algorithm
Kang et al.Potential of urban land use by autonomous vehicles: Analyzing land use potential in Seoul capital area of Korea
CN112419711A (en) Prediction method of parking demand in closed parking lot based on improved GMDH algorithm
Natera Orozco et al.Quantifying life quality as walkability on urban networks: The case of Budapest
CN112949997A (en)System and method for community portrayal in urban planning design
Motieyan et al.A novel spatial index using spatial analyses and hierarchical fuzzy expert system for obtaining green TOD: a case study in Tehran city
CN113487465B (en)City overlapping structure characteristic detection method and system based on label propagation algorithm
Gao et al.Integrating multisource geographic big data to delineate urban growth boundary: a case study of Changsha
CN114139827B (en)Intelligent perception and optimization method for urban functional area function performance
CN113158084B (en)Method, device, computer equipment and storage medium for processing movement track data
CN109800903A (en)A kind of profit route planning method based on taxi track data
CN116415756B (en)Urban virtual scene experience management system based on VR technology
FengImplementation of Decision Support System for Ecological Environment Planning of Urban Green Space

Legal Events

DateCodeTitleDescription
PB01Publication
PB01Publication
SE01Entry into force of request for substantive examination
SE01Entry into force of request for substantive examination
GR01Patent grant
GR01Patent grant
CF01Termination of patent right due to non-payment of annual fee

Granted publication date:20220930

CF01Termination of patent right due to non-payment of annual fee

[8]ページ先頭

©2009-2025 Movatter.jp