Street Centralities and Land Use Intensities Based on Points of Interest ( POI ) in Shenzhen , China

Urban land use and transportation are closely associated. Previous studies have investigated the spatial interrelationship between street centralities and land use intensities using land cover data, thus neglecting the social functions of urban land. Taking the city of Shenzhen, China, as a case study, we used reclassified points of interest (POI) data to represent commercial, public service, and residential land, and then investigated the varying interrelationships between the street centralities and different types of urban land use intensities. We calculated three global centralities (“closeness”, “betweenness”, and “straightness”) as well as local centralities (1-km, 2-km, 3-km, and 5-km searching radiuses), which were transformed into raster frameworks using kernel density estimation (KDE) for correlation analysis. Global closeness and straightness are high in the urban core area, and roads with high global betweenness outline the skeleton of the street network. The spatial patterns of the local centralities are distinguished from the global centralities, reflecting local location advantages. High intensities of commercial and public service land are concentrated in the urban core, while residential land is relatively scattered. The bivariate correlation analysis implies that commercial and public service land are more dependent on centralities than residential land. Closeness and straightness have stronger abilities in measuring the location advantages than betweenness. The centralities and intensities are more positively correlated on a larger scale (census block). These findings of the spatial patterns and interrelationships of the centralities and intensities have major implications for urban land use and transportation planning.


Introduction
Urban land use and transportation are two crucial subsystems within urban systems that mutually interact and influence each other [1][2][3][4].Urban roads facilitate socioeconomic activities taking place on urban land, and urban land use, in turn, influences the travel behavior of citizens [5][6][7].Urban socioeconomic activities and land use are never evenly distributed in self-organized cities nor planned cities, and urban roads are the same.Therefore, the spatial correlation between urban land use and urban roads has long been a major research interest of urban researchers and planners [8][9][10].
Street networks, composed of urban roads, are called the skeleton of a city [11].A street network is a typical complex network with two different representations, duel representation and primal representation, which can be characterized and quantified by metrics and topological measurements [12,13].Centrality is a fundamental concept in graph theory and network science, which identifies the most important vertices within a graph, for instance, the key infrastructure nodes in urban street networks [14,15].Street centrality captures location advantages in a city and plays a crucial role in shaping the intraurban variation of urban structures and land uses [16].Previous studies have examined the interrelationship between different centrality indicators ("betweenness", "closeness", and "straightness") and economic activities or general land uses intensities within a city [17][18][19].For example, Wang et al. [16] analyzed the relationship between street centralities and land use intensities, which are represented by population densities and employment densities.We seek to further understand how street centralities specifically influence different types of urban land use, such as commercial land, residential land, and industrial land.
However, to the best of our knowledge, there are few studies focusing on the influence of centrality indicators on different types of urban land use.Rui and Ban (2014) found variant correlations between different centralities and land use types in Stockholm using land-cover data [20].Chaudhuri and Clarke [21] used land-cover data (urban areas, natural areas, agricultural land use, and water bodies) to investigate the spatiotemporal dynamics of the coupling between land use (change) and street networks.However, land cover data can only quantify the spatial extent of land use and cannot characterize the intensity of land use.Urban land use can be differentiated either by their physical properties or social functions [22,23].Remote sensing has long been used to acquire the physical properties of land use and to classify land into land cover types [24], which neglects knowing the social functions of urban land use [22].Emerging types of big data, such as mobile phone data, points of interest (POI), trajectories, and social media data, have been used to acquire fine urban land use classifications [22,[25][26][27][28][29][30][31], providing solutions for acquiring different types of urban land use.
A point of interest (POI) is a specific point location, referring to all geographical entities that can be abstracted as points.Previous studies have suggested that POI data is capable of describing land use at a disaggregated level and has a finer grain than a conventional land use map [32].Using POI data can not only identify the types of urban land use, but also their intensities.In this study, taking the city of Shenzhen, China, as a case study, we use points of interest (POI) to represent different urban land use types and their intensities, namely, commercial land, public service land, and residential land.We focus on the distinguished interrelationships between street centralities and different urban land use intensities.

Study Area
Shenzhen, located in the south of China, is one of the youngest megacities in China (Figure 1).Thanks to China's Open Policy reform, Shenzhen has experienced rapid growth over the past few decades.In 2016, the permanent resident population in Shenzhen was 12 million, with a gross domestic product (GDP) of approximately 300 billion USD [33].There are ten districts in Shenzhen, namely: Luohu, Futian, Nanshan, Yantian, Baoan, Longgang, Longhua, Guangming, Pingshan, and Dapeng (Figure 1), among which Luohu, Futian, Nanshan, and Yantian are the four initial special economic zones in Shenzhen [34].Presently, Futian and Luohu are the downtown areas and Nanshan is the high-technology zone in Shenzhen City [35].
Shenzhen is not built on China's vast plains, but instead around hills along the coastline (Figure 1).Restricted by its special terrain, Shenzhen implemented a multicluster urban master plan at the beginning of its establishment [32,35].In each urban cluster, planners have assigned and designed employment and residential areas to minimize commuting across clusters.Under such circumstances, the roads carrying commuting across clusters are particularly important in the street network.

Street Network
Street centralities' characterization heavily relies on the high-quality data of urban roads.OpenStreetMap, which uses edits and updates of roads across the world done by volunteers, can be a vital approach for acquiring urban roads, and has been successfully used in other studies [36][37][38][39][40]. Nevertheless, there are two obstacles (typology and file format) in building a street network using OpenStreetMap road data.Recently, Boeing has developed a useful toolbox (OSMnx) based on Python programming, which can easily download a street network from OpenStreetMap, check and correct the network topology, and analyze the street networks [38].Specifically, this toolbox produces nodes only at locations where roads intersect, thus excluding overpasses and tunnels.This process results in a nonplanar graph, while most of the previous studies simplify the street networks as planar networks [38].
We used OSMnx to download the drivable urban roads in Shenzhen City on 8 February 2018 [38].There are 28,064 nodes (intersections) and 41,789 edges (links between two intersections) in the street network.The road length (length of edges) and kernel density of the nodes (using kernel density estimation (KDE) in a default searching radius) are presented in Figure 2a,b, respectively.The road length varies from 0.2 m to 13.6 km, with a mean of 223.1 m.The total length of the roads is 9323 km, and the average road density in Shenzhen is 4.67 km/km 2 (the land area in Shenzhen City is 1997 km 2 ).The length of the edges is relatively short in the core area, while long edges are roads across districts (Figure 2a).A high density of nodes are concentrated in the core area (Futian, Nanshan, and Luohu districts) (Figure 2b).In the following section, we analyzed the correlation between the street centralities and land use intensities in Shenzhen, excluding the Dapeng district, because of the sparse roads there (Figure 2a).

Street Network
Street centralities' characterization heavily relies on the high-quality data of urban roads.OpenStreetMap, which uses edits and updates of roads across the world done by volunteers, can be a vital approach for acquiring urban roads, and has been successfully used in other studies [36][37][38][39][40]. Nevertheless, there are two obstacles (typology and file format) in building a street network using OpenStreetMap road data.Recently, Boeing has developed a useful toolbox (OSMnx) based on Python programming, which can easily download a street network from OpenStreetMap, check and correct the network topology, and analyze the street networks [38].Specifically, this toolbox produces nodes only at locations where roads intersect, thus excluding overpasses and tunnels.This process results in a nonplanar graph, while most of the previous studies simplify the street networks as planar networks [38].
We used OSMnx to download the drivable urban roads in Shenzhen City on 8 February 2018 [38].There are 28,064 nodes (intersections) and 41,789 edges (links between two intersections) in the street network.The road length (length of edges) and kernel density of the nodes (using kernel density estimation (KDE) in a default searching radius) are presented in Figure 2a,b, respectively.The road length varies from 0.2 m to 13.6 km, with a mean of 223.1 m.The total length of the roads is 9323 km, and the average road density in Shenzhen is 4.67 km/km 2 (the land area in Shenzhen City is 1997 km 2 ).The length of the edges is relatively short in the core area, while long edges are roads across districts (Figure 2a).A high density of nodes are concentrated in the core area (Futian, Nanshan, and Luohu districts) (Figure 2b).In the following section, we analyzed the correlation between the street centralities and land use intensities in Shenzhen, excluding the Dapeng district, because of the sparse roads there (Figure 2a).

Points of Interest (POI)
We collected the POI data in Shenzhen City from Amap (https://www.amap.com/) on 28 February 2018, with a total of 611,122 records.These POIs are divided into 15 primary categories and 81 subcategories.The numbers and proportions of POIs in the 15 primary categories are presented in Figure 3.The core function of a city is to provide urban residents with housing and employment, as well as public services.On the other hand, from the perspective of urban land functions, the most important land types are commercial land, residential land, and industrial land.Unfortunately, the POI data we obtained cannot accurately identify the industrial land in Shenzhen City (Figure 3).Therefore, we summarized the POI data into three categories, namely, commercial and business facilities, residential sites, and administration and public services (see detailed information in Figure 3).We used the three categories of POIs to represent commercial land, residential land, and public service land, respectively (Figure 3).

Points of Interest (POI)
We collected the POI data in Shenzhen City from Amap (https://www.amap.com/) on 28 February 2018, with a total of 611,122 records.These POIs are divided into 15 primary categories and 81 subcategories.The numbers and proportions of POIs in the 15 primary categories are presented in Figure 3.The core function of a city is to provide urban residents with housing and employment, as well as public services.On the other hand, from the perspective of urban land functions, the most important land types are commercial land, residential land, and industrial land.Unfortunately, the POI data we obtained cannot accurately identify the industrial land in Shenzhen City (Figure 3).Therefore, we summarized the POI data into three categories, namely, commercial and business facilities, residential sites, and administration and public services (see detailed information in Figure 3).We used the three categories of POIs to represent commercial land, residential land, and public service land, respectively (Figure 3).

Points of Interest (POI)
We collected the POI data in Shenzhen City from Amap (https://www.amap.com/) on 28 February 2018, with a total of 611,122 records.These POIs are divided into 15 primary categories and 81 subcategories.The numbers and proportions of POIs in the 15 primary categories are presented in Figure 3.The core function of a city is to provide urban residents with housing and employment, as well as public services.On the other hand, from the perspective of urban land functions, the most important land types are commercial land, residential land, and industrial land.Unfortunately, the POI data we obtained cannot accurately identify the industrial land in Shenzhen City (Figure 3).Therefore, we summarized the POI data into three categories, namely, commercial and business facilities, residential sites, and administration and public services (see detailed information in Figure 3).We used the three categories of POIs to represent commercial land, residential land, and public service land, respectively (Figure 3).

Street Centrality Measures
Centrality plays a crucial role in understanding the structural properties of the complex relational network [15].The multiple centrality assessment (MCA) method has been broadly used to analyze urban street networks based on primal graphs, in which the road intersections are turned into nodes, and the connections between the nodes have been turned into edges [16,18,19,41].There are three centrality indicators, namely, the closeness centrality, betweenness centrality, and straightness centrality, which have been widely used in related studies [16,[18][19][20].
Closeness centrality (C C ) quantifies to what extent a node (i) is close to all of the other nodes along the shortest paths of the network, which is defined as follows: where N is the total number of nodes in the network and d ij is the shortest path length between node i and j.
Betweenness centrality (C B ) quantifies the centrality of a node (i) by counting how many times the node (i) is traversed by the shortest paths of all of the pairs of nodes in the network, which is defined as follows: where n jk is the number of shortest paths between nodes j and k, and n jk (i) is the number of these shortest paths that traverse node i.
Straightness centrality (C S ) quantifies the extent of a deviation of the distance from the shortest path and the straight line between node i and j, which is defined as follows: where d Eucl ij is the Euclidean distance between node i and j.A higher straightness centrality (C s ) of node (i) means a better directness to the node, resulting a higher commuting efficiency.
We calculated these centrality indicators at each node using urban network analysis tools [41], which can be executed as a toolbox in a GIS software.The centrality indicator of a node is largely affected by the searching radius.The default searching radius is infinite, and when calculating a node's centrality indicator, the arithmetic will traverse all of the nodes in the road network, resulting in global centralities.Local centralities with specific searching radiuses (1 km, 2 km, 3 km, and 5 km) are also calculated, which will be presented in the discussion section.

Kernel Density Estimation (KDE)
As the calculation of street centralities is based on nodes (points), and POI data is also the point data, which are unevenly distributed, we need centralities and land use intensities across the whole study area in order to perform a correlation analysis.Among the numerous methods of spatial smoothing and spatial interpolation, kernel density estimation (KDE) has been commonly used to transform the results of street centralities and land use intensities into continuous raster frameworks [17,18].KDE is a nonparametric way to estimate the probability density function of a random variable.In GIS, KDE uses the density within a range (window) of each observation in order to calculate the value at the center of the window by weighting nearby observations more than distant observations based on a kernel function [17].The definition of the equation and a detailed explanation of KDE can be found in a previous relevant paper [17].
We used KDE to acquire the raster frameworks of the centrality indicators and POI-based urban land use intensities.We first used the default searching radius of the KDE method in ArcGIS10.2,which is the shortest of the width or height of the output extent in the output spatial reference divided by 30.In this study, the searching radius in KDE ranges from 1145 m (for commercial and business facilities) to 1400 m (for residential sites).We calculated the urban land use intensities under different searching radiuses in KDE (1 km, 2 km, 3 km, and 5 km).The centrality indicators and land use intensities under different searching radiuses in KDE are highly correlated.Thus, we only compare the correlations between the centrality indicators and land use intensities using the raster frameworks under the default searching radius in KDE.The use of the default searching radius in KDE is to avoid any confusion with the searching radiuses in the calculation of the local centralities.The cell size for the output raster dataset is defined as the default value in ArcGIS10.2,which is the shorter of the width or height of the output extent in the output spatial reference divided by 250.In this study, the height of the output extent is shorter, which is approximately 45 km, resulting in a cell size of approximately 180 m.
After acquiring raster frameworks, we used the 1 km grid and subdistricts (census blocks) in Shenzhen to summarize the centrality indicators and land use intensities.Then, we used the Pearson's correlation analysis to investigate the correlation between the street centralities and land use intensities.We first presented the results of correlation analysis using a 1 km grid (N = 1621), and we compared that with the correlation results using census blocks (N = 58) in order to check the influences of the units in the correlation analysis.The overall flowchart of this study is shown in Figure 4. We used KDE to acquire the raster frameworks of the centrality indicators and POI-based urban land use intensities.We first used the default searching radius of the KDE method in ArcGIS10.2,which is the shortest of the width or height of the output extent in the output spatial reference divided by 30.In this study, the searching radius in KDE ranges from 1145 m (for commercial and business facilities) to 1400 m (for residential sites).We calculated the urban land use intensities under different searching radiuses in KDE (1 km, 2 km, 3 km, and 5 km).The centrality indicators and land use intensities under different searching radiuses in KDE are highly correlated.Thus, we only compare the correlations between the centrality indicators and land use intensities using the raster frameworks under the default searching radius in KDE.The use of the default searching radius in KDE is to avoid any confusion with the searching radiuses in the calculation of the local centralities.The cell size for the output raster dataset is defined as the default value in ArcGIS10.2,which is the shorter of the width or height of the output extent in the output spatial reference divided by 250.In this study, the height of the output extent is shorter, which is approximately 45 km, resulting in a cell size of approximately 180 m.
After acquiring raster frameworks, we used the 1 km grid and subdistricts (census blocks) in Shenzhen to summarize the centrality indicators and land use intensities.Then, we used the Pearson's correlation analysis to investigate the correlation between the street centralities and land use intensities.We first presented the results of correlation analysis using a 1 km grid (N = 1621), and we compared that with the correlation results using census blocks (N = 58) in order to check the influences of the units in the correlation analysis.The overall flowchart of this study is shown in Figure 4.

Spatial Patterns of Street Centralities
We first calculated three global centrality indicators at each node using urban network analysis tools [41].Then, the centralities of the edges are summarized using the average of the centralities of the two nodes linking to this edge [16].The spatial distributions of the three centrality indicators using KDE in the default searching radiuses are presented in Figure 5.The frequency distributions of the three centralities indicators at the node level are presented in Figure 6.

Spatial Patterns of Street Centralities
We first calculated three global centrality indicators at each node using urban network analysis tools [41].Then, the centralities of the edges are summarized using the average of the centralities of the two nodes linking to this edge [16].The spatial distributions of the three centrality indicators using KDE in the default searching radiuses are presented in Figure 5.The frequency distributions of the three centralities indicators at the node level are presented in Figure 6.Global closeness centrality (Glob C C ) presents a concentric ring form, which decreases from the center to the outside (Figure 5a).The highest value of Glob C C in the core area reflects the location advantage within the whole network (Figure 5a).The results of KDE show that high values of Glob C C are concentrated in the urban core area, namely, the Futian, Luohu, and Nanshan Districts (Figure 5b).In Longhua District, which is relatively far from the urban core area, there is also a high-value area of Glob C C .This is the Foxconn Technology Park, where urban roads are dense and have a high Glob C c (Figure 5b).
The spatial patterns of global betweenness centralities (Glob C B ) are obviously different from those of Glob C C .Most of the roads show a very low Glob C B (Figures 5c and 6b), and only important roads have a high Glob C B , which can be more clearly seen in the KDE results (Figure 5d).Shennan Avenue and North Ring Avenue are the east-west arteries, playing an important role in connection within the core urban area; thus, the two roads present a high Glob C B .Other highways or expressways connecting different urban clusters also have a high Glob C B .Those roads with a low   Global closeness centrality (Glob C C ) presents a concentric ring form, which decreases from the center to the outside (Figure 5a).The highest value of Glob C C in the core area reflects the location advantage within the whole network (Figure 5a).The results of KDE show that high values of Glob C C are concentrated in the urban core area, namely, the Futian, Luohu, and Nanshan Districts (Figure 5b).In Longhua District, which is relatively far from the urban core area, there is also a high-value area of Glob C C .This is the Foxconn Technology Park, where urban roads are dense and have a high Glob C c (Figure 5b).
The spatial patterns of global betweenness centralities (Glob C B ) are obviously different from those of Glob C C .Most of the roads show a very low Glob C B (Figures 5c and 6b), and only important roads have a high Glob C B , which can be more clearly seen in the KDE results (Figure 5d).Shennan Avenue and North Ring Avenue are the east-west arteries, playing an important role in connection within the core urban area; thus, the two roads present a high Glob C B .Other highways or expressways connecting different urban clusters also have a high Glob C B .Those roads with a low Global closeness centrality (Glob C C ) presents a concentric ring form, which decreases from the center to the outside (Figure 5a).The highest value of Glob C C in the core area reflects the location advantage within the whole network (Figure 5a).The results of KDE show that high values of Glob C C are concentrated in the urban core area, namely, the Futian, Luohu, and Nanshan Districts (Figure 5b).In Longhua District, which is relatively far from the urban core area, there is also a high-value area of Glob C C .This is the Foxconn Technology Park, where urban roads are dense and have a high Glob C c (Figure 5b).
The spatial patterns of global betweenness centralities (Glob C B ) are obviously different from those of Glob C C .Most of the roads show a very low Glob C B (Figures 5c and 6b), and only important roads have a high Glob C B , which can be more clearly seen in the KDE results (Figure 5d).Shennan Avenue and North Ring Avenue are the east-west arteries, playing an important role in connection within the core urban area; thus, the two roads present a high Glob C B .Other highways or expressways connecting different urban clusters also have a high Glob C B .Those roads with a low Glob C B are not predominant in the whole network, but they could have important connections in the local areas.This dynamic will be reflected in the local centralities.
Global straightness centralities (Glob C S ) are high in the core urban area and in the northeast and west of the city.The typical road length in the urban core areas is short, partly resulting in the relatively high Glob C S .In contrast to the frequency distribution of Glob C B , more nodes have a higher Glob C S (Figure 6b,c).The KDE result shows that a high Glob C S is concentrated in the core urban area, similar to the spatial pattern of Glob C C .

Spatial Distributions of Urban Land Use Intensities
We used the reclassified POI data to represent the following three types of urban land use: commercial land, public service land, and residential land.We used KDE to generate the land use intensity raster frameworks shown in Figure 7.The frequency distributions of land use intensities in a 1 km grid are shown in Figure 8.
Compared to the spatial patterns of the centrality indicators, high intensities of urban land use are more concentrated in smaller areas (Figure 7).Specific the three types of POI-based urban land use intensities, there are also apparent disparities in their spatial patterns.The intensities of commercial land and public services land are extremely concentrated in the following two urban areas: the Luohu-Futian area and Nanshan-Baoan area (Figure 7a,b).In contrast, the agglomeration of residential land is not as obvious as the other two, and its spatial distribution is relatively scattered.In addition to the centralized distribution of residential land in the core urban area, Longhua and Longgang, which are relatively far from the core urban area, also have a large amount of residential land (Figure 7c).Glob C B are not predominant in the whole network, but they could have important connections in the local areas.This dynamic will be reflected in the local centralities.Global straightness centralities (Glob C S ) are high in the core urban area and in the northeast and west of the city.The typical road length in the urban core areas is short, partly resulting in the relatively high Glob C S .In contrast to the frequency distribution of Glob C B , more nodes have a higher Glob C S (Figure 6b,c).The KDE result shows that a high Glob C S is concentrated in the core urban area, similar to the spatial pattern of Glob C C .

Spatial Distributions of Urban Land Use Intensities
We used the reclassified POI data to represent the following three types of urban land use: commercial land, public service land, and residential land.We used KDE to generate the land use intensity raster frameworks shown in Figure 7.The frequency distributions of land use intensities in a 1 km grid are shown in Figure 8.
Compared to the spatial patterns of the centrality indicators, high intensities of urban land use are more concentrated in smaller areas (Figure 7).Specific to the three types of POI-based urban land use intensities, there are also apparent disparities in their spatial patterns.The intensities of commercial land and public services land are extremely concentrated in the following two urban areas: the Luohu-Futian area and Nanshan-Baoan area (Figure 7a,b).In contrast, the agglomeration of residential land is not as obvious as the other two, and its spatial distribution is relatively scattered.In addition to the centralized distribution of residential land in the core urban area, Longhua and Longgang, which are relatively far from the core urban area, also have a large amount of residential land (Figure 7c).

Correlations between Street Centrality and Land Use Intensity
We summarized the street centralities and land use intensities in each 1 km grid.Scatter plots of street centralities and land use intensities excluding the grids with zero values of centralities or intensities are presented in Figure 9.The straight line is the linear regression line between the street centralities and land use intensities.Overall, the street centralities and land use intensities are positively correlated.We used a bivariate correlation analysis (Pearson's r) to test the degree of correlation between the street centralities and land use intensities using SPSS 21.We also calculated the Pearson's r of the logarithmic form of centralities or land use intensities.The results of the bivariate correlation analysis are presented in Table 1.

Correlations between Street Centrality and Land Use Intensity
We summarized the street centralities and land use intensities in each 1 km grid.Scatter plots of street centralities and land use intensities excluding the grids with zero values of centralities or intensities are presented in Figure 9.The straight line is the linear regression line between the street centralities and land use intensities.Overall, the street centralities and land use intensities are positively correlated.We used a bivariate correlation analysis (Pearson's r) to test the degree of correlation between the street centralities and land use intensities using SPSS 21.We also calculated the Pearson's r of the logarithmic form of centralities or land use intensities.The results of the bivariate correlation analysis are presented in Table 1.

Correlations between Street Centrality and Land Use Intensity
We summarized the street centralities and land use intensities in each 1 km grid.Scatter plots of street centralities and land use intensities excluding the grids with zero values of centralities or intensities are presented in Figure 9.The straight line is the linear regression line between the street centralities and land use intensities.Overall, the street centralities and land use intensities are positively correlated.We used a bivariate correlation analysis (Pearson's r) to test the degree of correlation between the street centralities and land use intensities using SPSS 21.We also calculated the Pearson's r of the logarithmic form of centralities or land use intensities.The results of the bivariate correlation analysis are presented in Table 1.Generally, the street centralities and land use intensities are significantly correlated in all four forms (x~y, Lnx~y, x~Lny, and Lnx~Lny), and their logarithmic forms (Lnx~Lny) show the highest correlation (Table 1).The street centralities have varying impacts on the different types of urban land use.From the perspective of centralities, Glob C C and Glob C S have a higher Pearson's r with land use intensities than that of Glob C B , which implies that the intensities of urban land use are more dependent on the Glob C C and Glob C S of the street network.As previously mentioned, Glob C B measures the nodes that play key roles in connecting the entire street network, and the nodes with an important connectivity may not be active places for socioeconomic activities, which leads to the relatively low correlation between the Glob C B and urban land use intensities.
In terms of different urban land use, commercial land and public service land have higher Pearson's r values with centralities than that of residential land.Compared with the centralized distribution of commercial land and public service land, the spatial distribution of residential land is scattered, resulting in a lower correlation with street centralities.

Local Street Centralities
To test the correlation between the local street centralities and the land use intensities, we calculated three centrality indicators for each code under different searching radiuses of 1 km, 2 km, 3 km, and 5 km.The spatial distributions of the KDE of the local centralities are presented in Figure 10.Generally, the spatial patterns of the local centralities are different from those of the global centralities.The local centralities show the location advantage of the nodes within a local area.In particular, the local C B measures the connectivity of a road within a local area, whose spatial patterns are more scattered compared to global C B .The high values of local C C and local C S are mainly concentrated in the core area, but there are also obvious locally high values far from the core area (Figure 10).With the increase in the searching radius in the centralities calculations, the spatial pattern of the centralities is closer to the global centralities pattern, thus showing a transformation of local characteristics to the overall features.
We also summarized the local street centralities in the 1 km grid and analyzed their correlations between the local centralities under different searching radiuses and urban land use intensities.Due to the similar degree of correlation in the different forms (Table 1), we only tested the linear correlation between the local centralities and intensities in this study (Table 2).The linear correlation coefficients between the global centralities and land use intensities are also included in Table 2 for comparison.The correlations between local C C (local C B ) and the land use intensities show a decreasing trend as the searching radius increases from 1 km to 5 km, while the correlations between local C S and the land use intensities tend to increase as the searching radius increases (Table 2).
In comparisons of the global and local centralities, local C C 1km is more correlated with the intensities of the commercial land and public service land than global C C , but it is opposite for the residential land.Global C B has a higher correlation with land use intensities, indicating the higher impact on the urban land use intensities compared to local C B .However, local C S has a higher correlation with urban land use intensities.In summary, the optimal scales of correlation between the commercial land and the three centralities are local C C 1km , global C B , and local C S 5km , which are shared by public service land.The optimal scales of correlation between the residential land and the three centralities are global C C , global C B , and local C S 5km .These Pearson's r values under the optimal scales are in bold in Table 2.
commercial land and the three centralities are local C C 1km, global C B , and local C S 5km, which are shared by public service land.The optimal scales of correlation between the residential land and the three centralities are global C C , global C B , and local C S 5km.These Pearson's r values under the optimal scales are in bold in Table 2.

of Statistical Units
To test the impact of the different statistical units on the correlation between the street centralities and land use intensities, we summarized the street centralities and land use intensities using the Shenzhen subdistrict administrative units (census blocks).The spatial distributions of the street centralities and land use intensities at the census block level are shown in Figure 11.We only calculated the linear correlation between the street centralities and urban land use intensities under the optimal scales, which are presented in Table 3. Dapeng District is excluded from the correlation analysis for consistency, resulting in 58 census blocks in total.
The street centralities in the census blocks show an apparent concentration of high centralities in the core area (blue ellipses) in Figure 11a,c,e.The high intensities of commercial land and public service land are also concentrated in the core urban area (Figure 11b,d).However, for residential land, the Longhua and Minzhi blocks have high intensities, which are away from the core area (Figure 11f).Urban land intensities in Nanshan Disctrict (red circles), the high-technology zone, are relatively low compared to the Futian and Luohu districts, the initial special economic zones.The results of the correlation analysis show that the street centralities and urban land use intensities at the subdistrict level have a higher Pearson's r compared with that in the 1 km grid (Tables 2 and 3).The subdistricts that have a larger scale than the 1 km grid and can eliminate local noise, resulting in a higher degree of correlation.

Effect of Statistical Units
To test the impact of the different statistical units on the correlation between the street centralities and land use intensities, we summarized the street centralities and land use intensities using the Shenzhen subdistrict administrative units (census blocks).The spatial distributions of the street centralities and land use intensities at the census block level are shown in Figure 11.We only calculated the linear correlation between the street centralities and urban land use intensities under the optimal scales, which are presented in Table 3. Dapeng District is excluded from the correlation analysis for consistency, resulting in 58 census blocks in total.
The street centralities in the census blocks show an apparent concentration of high centralities in the core area (blue ellipses) in Figure 11a,c,e.The high intensities of commercial land and public service land are also concentrated in the core urban area (Figure 11b,d).However, for residential land, the Longhua and Minzhi blocks have high intensities, which are away from the core urban area (Figure 11f).Urban land intensities in Nanshan Disctrict (red circles), the high-technology zone, are relatively low compared to the Futian and Luohu districts, the initial special economic zones.The results of the correlation analysis show that the street centralities and urban land use intensities at the subdistrict level have a higher Pearson's r compared with that in the 1 km grid (Tables 2 and 3).The subdistricts that have a larger scale than the 1 km grid and can eliminate local noise, resulting in a higher degree of correlation.

6.
Urban land use structure and transport systems mutually interact.Previous studies on the interrelationships between street centralities and land use intensities relied on land cover data, neglecting the social functions of urban land use.Taking Shenzhen City as a case study, this study contributed to distinguishing the different types of urban land use using reclassified points of interest (POI) data and investigating their varying correlations with three street centralities (closeness, betweenness, and straightness).The high values of the global closeness and straightness centralities are concentrated in the urban core areas, while the high values of the global betweenness centralities highlight the backbone of the road network.As the searching radius increases, the spatial patterns of the local centralities are more similar to the global centralities, showing a transformation of local characteristics to overall features.The results of our correlation analysis show a general strong correlation between urban land intensities and street centralities.Specifically, commercial land and public service land are more correlated with street centralities than residential land, because the spatial distribution of residential land is relatively scattered compared to the other two types of urban land use.We also found that the closeness and straightness centralities have stronger abilities in measuring locational advantages than the betweenness centrality.The statistical units can affect the correlation between the street centralities and the land use intensities, and they are more correlated in a larger scale unit because of their reduced noise and ignored extraneous details.The analysis of the spatial patterns of street centralities and urban land use intensities, as well as their interdependence can help the implementation of urban land use planning and transportation planning.
This study uses POI data to characterize the different types of urban land use.However, the results are not verified by the determined data, which may reduce the feasibility of the result.Secondly, Shenzhen City is adjacent to Dongguan City and Huizhou City, and the residents of the three cities commute between each other frequently, especially at the junction of Shenzhen and Dongguan.This study only obtained the street network within the administrative boundaries of Shenzhen, resulting in a bias for the street centralities calculation at the edge of the city boundary, also known as the "edge effect" [42].Finally, we only tested correlations between the street centralities and land use intensities and their variations.Their spatial distributions and interrelationships between them are complex, which can be influenced by regulation, history, as well as many other factors [21].Further studies could help explain the mechanism of interdependence between urban land use and transportation.

Figure 1 .
Figure 1.Location of Shenzhen City in China and its spatial extent.

Figure 1 .
Figure 1.Location of Shenzhen City in China and its spatial extent.

Figure 2 .
Figure 2. Spatial distribution of the road net and road lengths in Shenzhen City (a) and density of nodes in Shenzhen City (b).The nodes density used the kernel density estimation (KDE) in a default searching radius.

Figure 3 .
Figure 3. Numbers and proportions of 15 primary categories of points of interest (POIs) in Shenzhen City.These POIs are divided into commercial and business facilities (61.6%), residential sites (6.9%), administration and public services (8.1%), and others (23.4%).

Figure 2 .
Figure 2. Spatial distribution of the road net and road lengths in Shenzhen City (a) and density of nodes in Shenzhen City (b).The nodes density used the kernel density estimation (KDE) in a default searching radius.

Figure 2 .
Figure 2. Spatial distribution of the road net and road lengths in Shenzhen City (a) and density of nodes in Shenzhen City (b).The nodes density used the kernel density estimation (KDE) in a default searching radius.

Figure 3 .
Figure 3. Numbers and proportions of 15 primary categories of points of interest (POIs) in Shenzhen City.These POIs are divided into commercial and business facilities (61.6%), residential sites (6.9%), administration and public services (8.1%), and others (23.4%).

Figure 3 .
Figure 3. Numbers and proportions of 15 primary categories of points of interest (POIs) in Shenzhen City.These POIs are divided into commercial and business facilities (61.6%), residential sites (6.9%), administration and public services (8.1%), and others (23.4%).

Figure 4 .
Figure 4. Flowchart of analysis of interrelationship between street centralities and POI-based urban land use intensities.

Figure 4 .
Figure 4. Flowchart of analysis of interrelationship between street centralities and POI-based urban land use intensities.

Figure 5 .
Figure 5. Spatial distributions of global centrality indicators of street networks in Shenzhen.The centrality indicators of each edge are averaged from two nodes linking to this edge, and they are divided into five levels using natural breaks in ArcGIS10.2(a,c,e); kernel density estimation of centrality indicators at nodes (b,d,f).

Figure 6 .
Figure 6.Frequency distribution of global centrality indicators at each node, global closeness (a), global betweenness (b), and global straightness (c).

Figure 5 .
Figure 5. Spatial distributions of global centrality indicators of street networks in Shenzhen.The centrality indicators of each edge are averaged from two nodes linking to this edge, and they are divided into five levels using natural breaks in ArcGIS10.2(a,c,e); kernel density estimation of centrality indicators at nodes (b,d,f).

Figure 5 .
Figure 5. Spatial distributions of global centrality indicators of street networks in Shenzhen.The centrality indicators of each edge are averaged from two nodes linking to this edge, and they are divided into five levels using natural breaks in ArcGIS10.2(a,c,e); kernel density estimation of centrality indicators at nodes (b,d,f).

Figure 6 .
Figure 6.Frequency distribution of global centrality indicators at each node, global closeness (a), global betweenness (b), and global straightness (c).

Figure 6 .
Figure 6.Frequency distribution of global centrality indicators at each node, global closeness (a), global betweenness (b), and global straightness (c).

Figure 7 .
Figure 7. Spatial distributions of intensities of commercial and business facilities (a), public and administrative services (b), and residential sites (c).These represent commercial land use intensities, public service land use intensities, and residential land use intensities, respectively.

Figure 7 .
Figure 7. Spatial distributions of intensities of commercial and business facilities (a), public and administrative services (b), and residential sites (c).These represent commercial land use intensities, public service land use intensities, and residential land use intensities, respectively.

Figure 8 .
Figure 8. Frequency distributions of land use intensities using kernel density estimation (KDE) in a 1 km grid.(a) Commercial land use intensity, (b) public service land use intensity, (c) residential land use intensity.

Figure 9 .
Figure 9. Scatter plots of street centralities and land use intensities.The same vertical column shares the same centrality indicator: Closeness, Betweenness, and Straightness, from left to right, respectively.The same horizontal row shares the same land use intensity: commercial land, public service land, and residential land, from top to bottom, respectively.The straight lines are linear regression lines.

Figure 8 .
Figure 8. Frequency distributions of land use intensities using kernel density estimation (KDE) in a 1 km grid.(a) Commercial land use intensity, (b) public service land use intensity, (c) residential land use intensity.

Figure 8 .
Figure 8. Frequency distributions of land use intensities using kernel density estimation (KDE) in a 1 km grid.(a) Commercial land use intensity, (b) public service land use intensity, (c) residential land use intensity.

Figure 9 .
Figure 9. Scatter plots of street centralities and land use intensities.The same vertical column shares the same centrality indicator: Closeness, Betweenness, and Straightness, from left to right, respectively.The same horizontal row shares the same land use intensity: commercial land, public service land, and residential land, from top to bottom, respectively.The straight lines are linear regression lines.

Figure 9 .
Figure 9. Scatter plots of street centralities and land use intensities.The same vertical column shares the same centrality indicator: Closeness, Betweenness, and Straightness, from left to right, respectively.The same horizontal row shares the same land use intensity: commercial land, public service land, and residential land, from top to bottom, respectively.The straight lines are linear regression lines.

Figure 10 .
Figure 10.Spatial patterns of KDE of the local street centralities using searching radiuses of 1 km, 2 km, 3 km, and 5 km.The same vertical column shares the same centrality indicator: Closeness, Betweenness, and Straightness, from left to right, respectively.The same horizontal row shares the same searching radius: 1 km, 2 km, 3 km, and 5 km, from top to bottom, respectively.

Figure 10 .
Figure 10.Spatial patterns of KDE of the local street centralities using searching radiuses of 1 km, 2 km, 3 km, and 5 km.The same vertical column shares the same centrality indicator: Closeness, Betweenness, and Straightness, from left to right, respectively.The same horizontal row shares the same searching radius: 1 km, 2 km, 3 km, and 5 km, from top to bottom, respectively.

Figure 11 .
Figure 11.Spatial distributions of street centralities and land use intensities at the subdistrict level (census block) in Shenzhen City.(a) Closeness at 1 km searching radius, (b) Commercial land use intensity, (c) Global Betweenness, (d) Public service land use intensity, (e) Straightness at 2 km searching radius, (f) Residential land use intensity.

Figure 11 .
Figure 11.Spatial distributions of street centralities and land use intensities at the subdistrict level (census block) in Shenzhen City.(a) Closeness at 1 km searching radius, (b) Commercial land use intensity, (c) Global Betweenness, (d) Public service land use intensity, (e) Straightness at 2 km searching radius, (f) Residential land use intensity.

Table 1 .
Pearson's r between global street centralities and land use intensities in a 1 km grid.Publ, and Resi mean commercial land, public service land, and residential land, respectively.** p < 0.01.

Table 2 .
Pearson's r between local street centralities and land use intensities in the 1 km grid.Comm, Publ, and Resi mean commercial land, public service land, and residential land, respectively.## Global centralities are also included for the convenience of comparison.** p < 0.01. #

Table 2 .
Pearson's r between local street centralities and land use intensities in the 1 km grid.
# Comm, Publ, and Resi mean commercial land, public service land, and residential land, respectively.## Global centralities are also included for the convenience of comparison.** p < 0.01.

Table 3 .
Pearson's r between street centralities and land use intensities at the subdistrict level (census block) (N = 58).

Table 3 .
Pearson's r between street centralities and land use intensities at the subdistrict level (census block) (N = 58).