Next Article in Journal
Biologically Inspired Hierarchical Contour Detection with Surround Modulation and Neural Connection
Previous Article in Journal
Multiparametric Monitoring in Equatorian Tomato Greenhouses (III): Environmental Measurement Dynamics
Previous Article in Special Issue
Recent Surface Water Extent of Lake Chad from Multispectral Sensors and GRACE
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Methods of Population Spatialization Based on the Classification Information of Buildings from China’s First National Geoinformation Survey in Urban Area: A Case Study of Wuchang District, Wuhan City, China

1
School of Remote Sensing and Information Engineering, Wuhan University, Wuhan 430072, China
2
State Key Laboratory of Information Engineering in Surveying, Mapping & Remote Sensing, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(8), 2558; https://doi.org/10.3390/s18082558
Submission received: 25 June 2018 / Revised: 25 July 2018 / Accepted: 3 August 2018 / Published: 4 August 2018
(This article belongs to the Special Issue Spatial Analysis and Remote Sensing)

Abstract

:
Most of the currently mature methods that are used globally for population spatialization are researched on a single level, and are dependent on the spatial relationship between population and land covers (city, road, water area, etc.), resulting in difficulties in data acquisition and an inability to identify precise features on the different levels. This paper proposes a multi-level population spatialization method on the different administrative levels with the support of China’s first national geoinformation survey, and then considers several approaches to verify the results of the multi-level method. This paper aims to establish a multi-level population spatialization method that is suitable for the administrative division of districts and streets. It is assumed that the same residential house has the same population density on the district level. Based on this assumption, the least squares regression model is used to obtain the optimized prediction model and accurate population space prediction results by dynamically segmenting and aggregating house categories.In addition, it is assumed that the distribution of the population is relatively regular in communities that are spatially close to each other, and that the population densities on the street level are similar, so the average population density is assessed by optimizing the community and surrounding residential houses on the street level. Finally, the scientificalness and rationality of the proposed method is proved by spatial autocorrelation analysis, overlay analysis, cross-validation analysis and accuracy assessment methods.

1. Introduction

Population data are one of the most direct indicators of human activity [1]. With the development of China’s urbanization process from 1949 to 2015, the proportion of the urban population in China increased from approximately 10% to 57.35% [2]. The spatial distribution of the population, population flow, and population structure are becoming increasingly important for the development of cities. The spatial distribution of the population influences not only the urbanization process and living environment [3,4], but also the development plan of the regional public education system, medical facilities, and other services, which are related to people’s vital interests [5,6,7].
The spatial distribution of the population is affected by many factors, such as geographic location, land cover, convenience of road networks, water areas, and economic development [8,9]. Therefore, traditional research methods mainly fit spatial population distributions by studying the coupling relationship between regional population density and its influence factors. Liao Shunbao et al. [10] examined the correlation between the population density and land use in Tibet and Qinghai Province and proposed a spatial model of population through multi-source data fusion method. Du Guoming et al. [11] used the data from the fifth census of Shenyang City and residential areas data in order to simulate population distribution through the spatial interpolation method. Given the shortcomings of the current spatial methods for urban populations, Kang Tingjun et al. [12] developed a multi-agent-based urban population distribution method. Using North Korea’s district-level census data, Shi Tingting et al. [13] analyzed the relationship between North Korea’s population density and spatial factors, and then performed multiple regression analysis to spatial status of North Korea’s population density. Dong Chun et al. [14] combined population statistics data with geographical data and economic data to establish apopulation spatialization method, which examines the coupling relationship between population distribution and related factors in a certain region.
Remote sensing imagery provides a new idea for population spatial modeling [15,16]. Many scholars use the advantages of remote sensing imagery, including its multilevel nature and high degree of timeliness, combined with geographic information system (GIS) technology to buildpopulation spatialization model at different spatial levels [17,18,19,20]. Chen Qing et al. [21] studied the correlation between night-time remote sensing images and geographic factors, performing a population spatialization experiment in the highly efficient eco-economic region in the Yellow River Delta. Lo C P et al. [22] studied the relationship between the gray value of thematic mapper (TM) images in different bands and urban population density. Li Shujuan et al. [23] used high-resolution remote sensing image to extract building information for different functions and calculated the population accommodation coefficients of different buildings in order to establish the spatial distribution map of the urban population. Wang Shixin et al. [24] used three-dimensional (3D) reconstruction technology to identify and extract urban residents’ housing areas and elevation in ZY-3 images.
Typically, in population census data, the statistical unit is the administrative unit, thus, the statistical level is relatively coarse and the types of data are limited. Therefore, the selection and improvement of mathematical methods are crucial for obtaining high-precision population spatialization results. Commonly used methods include geostatistics methods [25], spatial regression models [26,27], spatial interpolation methods [28,29], and machine learning methods [30,31]. Holt et al. [32] used the improved population weight method to interpolate census data spatially, and this method can better explain the spatial distribution of the population within the census administrative division. Wang Keijing et al. [33] studied population spatialization by using multivariate statistical regression and geo-weighted regression (GWR) models. Cao Li-qin et al. [34] predicted the population of 76 districts or counties in Hubei Province in 2002 by using the neural network model to establish a relationship between the brightness of nighttime light data and the urban population.
The study of population spatialization has become more comprehensive given the integration of more data sources [35,36] and technological methods [37,38]. At present, there are a number of mature data sets of population spatialization achievements covering the world, countries, or regions, such as Landscan [39], Worldpop [40], and GHS-POP [41]. These data sets provide detailed and accurate results of population maps of dynamic population flow [42,43,44], age structure change [45,46,47], urbanization development [48,49,50], building or settlement characteristic information [51,52,53], and greatly promote the cross-study of population spatialization. By combining with other related fields, important data and method support are provided to guide the urban planning [42,54], to assess the risk of demographic risk [55,56], and to improve the population quality of life [57,58].
Facing thefact that the fusion of more and more data sources, the variety method of population spatialization and the difference perspective of population research, it is a very important direction of the future research to use suitable data and establish the population spatialization method to meet the needs of different administrative units. At present, many scholars have carried out a series of researches on data process and methods for the improvement of the data source precision [54], the cross validation of population spatialization method [45,47], and the evaluation of the experimental results [31]. Few people pay attention to the demand and difference of the population spatialization method under different administrative units. Based on this, this paper makes use of data sources to establish population spatialization methods under the perspective of different administrative units, and tries to establish a reasonable method system of population accuracy evaluation to verify the rationality of the experimental results.
The China’s first national geoinformation survey started in January 2013 and lasted three years. Its purpose was to systematically obtain the authoritative, objective, and accurate information on the geographic conditions of the country in order to provide an important data foundation for promoting ecological environmental protection and building a resource conserving and environmentally friendly society. Through synthetically using the global navigation satellite system (GNSS), aerospace remote sensing (RS) technology, GIS technology, and other modern surveying and mapping technology, the survey can dynamically and quantitatively recognize land surface morphology, land covers, build-up zones, and monitor the spatial distribution and development of resources, the environment, ecology, and economic factors. This data set mainly contains three types of data: land topography data (DEM, Slope data), land cover classification data (“LCA”, which contain 10 major categories, such as farmland, garden land, woodland, and more than 100 smaller categories), and social geographical units, including point of interests (“POIs ”, educational facilities, hospitals, and so on), administrative unit categories and vector boundaries and other urban integrated functional units (“BUCA”, “BUCP”). The greatest advantage of this data set is that it is highly accurate in building space information, such as building location, shape, and other characteristics, including building types and height. Such detailed information on building classification provides useful data for the study of population spatializationbased on the housing construction area.
The main innovations of this paper are as follows: (1) The spatial and attribute information of buildings in China’s first national geoinformation survey is fully mined. Through the combination of different administrative divisions and thresholds for the proportion of housing construction areas, this paper gives functional attributes to all buildings and screens out residential houses. (2) A multi-level population spatialization method that is applicable to different administrative unit levels is established. (3) Various methods are used to qualitatively and quantitatively study differences in the experimental results on thedifferent levels. Thecommon and differences are well analyzed and explained.

2. Data Acquisition and Preprocessing

2.1. China’s Administrative Classification

The current administrative divisions in China are as follows: The first level is the provincial administrative units, which mainly include provinces, autonomous regions, municipalities, and special administrative regions. The second level is the prefectural administrative regions, which are divided into prefecture-level cities and regions. The third level is the county-level administrative units, which mainly include municipal districts, counties, and county-level cities. The fourth level is the township-level administrative units, which mainly include streets and townships. The fifth level is the village-level administrative units, which mainly include communities and administrative villages. The latter is a group-level administrative village that is divided into natural village groups and community residential groups [59,60]. The research objects in this paper are mainly the street and community units under the municipal district, corresponding to the fourth and fifth levels of the administrative division. Residential quarters, however, are not administrative units; rather, contains residential houses and are equipped with commercial outlets, culture and education, entertainment, and other public facilities [61].

2.2. Data Acquisition

The specific data needed for the multi-level population spatialization method are shown in Table 1: (1) Spatial vector data, including street administrative units, community administrative units, and housing construction features, are mainly collected from China’s first national geoinformation survey. (2) Demographic data, including resident population in sub-districts and resident population in sub-communities, are mainly collected from the Wuhan Statistical Yearbook in 2015, The Sixth National Population Census, and the Wuhan Community Demographic Census in 2013.

2.3. Data Preprocessing

Due to differences in the administrative vector boundaries, as collected by China’s first national geoinformation survey and the Community Demographic Census, the administrative boundaries corresponding to house buildings are not completely consistent. At the same time, as some community administrative boundaries are adjusted or merged with other communities, a small part of the “community vacancy” area may exist in the data results, that is, there are no community attributes for the house buildings in this area. In addition, because of the error in the boundaries of the residential quarters, some house buildings do not have attributes of residential communities. Therefore, data preprocessing mainly corrects the boundaries of administrative units, and defines the three types of attributes, including street, community, and residential quarters, clearly corresponding to the house buildings. The schematic diagram of data preprocessing is showed in Figure 1.
(1)
Take the vector boundary of non-residential quarters from the BUCA layer and the BUCP layer as the standard. Then, review all of the house buildings in the LCA layer and delete those belonging to non-residential quarters by using spatial location query and attribute query functions.
(2)
Take the vector boundary of residential quarters from the BUCA and BUCP layer as the standard to merge the house buildings in the LCA layer. Then, add the corresponding residential quarters attribute information to these house buildings through the spatial location query and attribute query functions.
(3)
Review the remaining house buildings in the LCA layer. Then, take the community vector boundary as the standard to merge house buildings belonging to the same community and add the corresponding community attribute information through the spatial location query and attribute query functions.
(4)
Use the street boundary to determine the street information for all communities through the spatial location query and attribute query functions. Then, add the corresponding street attribute information that is not available at the community level and residential quarters.

3. Multi-Level Method and Experimental Verification

Following data preprocessing, the multi-level population spatialization method is established for large-level regions, such as the district level and small-level areas, such as the street level. Then, space autocorrelation analysis, overlay analysis, and cross-validation analysis are used to verify the rationality of these methods. The workflow of multi-level population spatialization method and verification is showed in Figure 2.

3.1. The Method on the District Level

The main basis for the population spatialization method on the district level is as follows: sub-streets in China are mostly located in the urban centers or urban development areas; the layout of this region’s residential houses is more uniform than that in other areas. According to the standard of China’s first national geoinformation survey, residential quarters that are of the same type have high degree of similarity in terms of floor height and construction area. Therefore, the population density of different types of residential houses should be estimated by least squares regression through the classification that is based on the assumption that residential houses of the same type have the same population density [62].
As showed in Figure 3, the first step in processing the data is to calculate the areas of residential houses with attributes of the residential quarter, and then count the total areas of residential houses on each street and record them in vector R1. The second step is to calculate the house building areas that do not have community attributes but have street attributes, and then count the total areas of house buildings on each street and record them in vector R2. The next step is to count all of the house building areas with street attributes, and then count the total areas of house buildings on each street and record them in vector R3. The residential population of each street is evaluated by using the Wuhan Statistical Yearbook and The Sixth National Population Census.
Field surveys and expert verification indicate that residential houses are not extracted fully if the total house building area with residential quarter properties is less than 20% of the total street house building area and if the house building area with non-community properties exceeds 20% of the total street house building area. Thus, the ratios R1/R3 and R2/R3 must be calculated. If the results of the ratios meet the above conditions, then all of the residential houses have been extracted and residential house areas have been stored in R1. Otherwise, the house buildings with community attributes and non-residential quarter attributes should be regarded as residential houses and addedinR1. Then, the residential houses should be classified into multi-floor buildings, multi-floor independent buildings, low-floor buildings, and low-floor independent buildings (the detailed description of house buildings is showed in Table 2), according to the attributes and enter the final residential house areas into Equation (1).
y i = a 1 S 1 i + a 2 S 2 i + a 3 S 3 i + a 4 S 4 i + Δ i ( i = 1 n )
where y i represents the resident population count for each street, a 1 a 4 represent different housing density coefficients, and S 1 i S 4 i represent the areas of the different type of residential houses. Δ i represents random error term under the assumption of normal distribution and n represents the number of streets.
First, the coefficients and goodness of fit of the results must be evaluated. If some of the coefficients are negative or the overall goodness of fit is less than 0.5, multi-floor independent houses and multi-floor houses should be combined into a new type of house building, and low-floor independent houses and low-rise houses should be merged into another new type of house building. Then, the residential houses should be reclassified. The improved model is shown in Equation (2):
y i = a 1 ( S 1 i + S 2 i ) + a 2 ( S 3 i + S 4 i ) + Δ i ( i = 1 n )
The fitted population can be calculated through the coefficient results. However, because there is a certain deviation between the fitted population and the actual population count, the coefficients must be corrected [14,38] by using Equation (3).
a i = ( y i / y i ) 1 × a i ( i = 1 , 2 )  
Finally, it is appropriate to establish 250-m spatial grid cell, and the corrected coefficients are used to estimate the population count in each grid cell based on the areas of the different residential houses.

3.2. The Method on the Street Level

As the second smallest administrative unit in China, the community covers between a few to more than ten residential quarters, and the community areas are much smaller than the street areas. Therefore, the residential house areas of the community are far smaller than the street areas. The least squares method is suitable when the sample size is large and the variation in the independent variable is small, the resulting error of this method may be too large in order to fit the estimated community population. The population distribution is relatively regular in a community with close spatial distances, and the type of residential houses is uniform; thus, residential houses that are in close distances to each other have similar population densities. In this paper, it is reasonable to calculate the average population density of residential houses on the community level [37,63].
Using the first step of the Figure 3, the total area of residential houses in each residential quarter should be counted and recorded in vector C1. Then, the total area of house buildings in each community should be counted and recorded in vector C2. The residential population of each community from the Wuhan Community Demographic Census is recorded in vector P. The workflow of the population spatialization method on the street level is showed in Figure 4.
Field surveys and expert verification indicate that residential houses are not fully extracted if the total house building area with residential quarter properties is less than 20% of the house building area with non-residential quarter properties, and if the average population density of residential houses is above one person per square meter. Thus, the ratios R4 = C1/C2 and R5 = P/C1 must be calculated. If the results of the ratios meet the above conditions, then all of the residential houses have been extracted and residential house areas have been stored in C1. Otherwise, the house buildings that have community attributes but do not have residential quarter attributes should be regarded as residential houses, and the area should be added the area into C1.Finally, the final results of the residential house areas should be entered into Equation (4).
y i = P i / C 1 i ( i = 1 , n )
Finally, it is appropriate to establish 50-m spatial grid cell, count the areas of residential houses in each grid cell and calculate the total population.

3.3. Experimental Result Verification

The results of different administrative units are obtained by using the multi-level population space model, but whether the results of the experiments are reasonable requires further verification. Therefore, in this paper, three qualitative and quantitative analysis methods are utilized. Detailed procedures are described in Section 3.3.1, Section 3.3.2 and Section 3.3.3.

3.3.1. Spatial Autocorrelation Analysis

Spatial autocorrelation refers to the potential interdependence of some variables within observation data in the same distribution area. The local Moran index (Lisa index) [64,65] is a classical algorithm that is used to detect local spatial autocorrelation. It can reflect the spatial aggregation condition well by calculating the index of the correlation between the spatial unit and the neighborhood. At the same time, this method can take the variability of the local state into account when calculating the global index. Therefore, in this paper, this method is used to explore the spatial autocorrelation of the population.

3.3.2. Spatial Overlay Analysis

Overlay analysis creates new feature layers by stacking two-layer or multilayer map elements, not only generating new spatial relationships between multiple features, but also linking their attributes [14,26]. Geographic locations (governments), educational resources (primary and secondary schools), medical and health resources (hospitals on the different levels, community service centers), and road networks (traffic trunks) are selected in this paper. By calculating the average shortest distance and coverage degree, this paper analyzes the dominant resource advantage and influence degree in the population agglomeration areas, demonstrating the rationality of the population spatialization results.
The average shortest distance refers to the hierarchical statistics of population spatialization results, and it is calculated as the average distance between the different population levels to the nearest government, educational resources, medical and health resources, and road network. Quantifying the linear relationship between different population levels and impact factors can reflect the rationality of the method, at the same time, can provide a deeper understanding of the impact of different factors on population distribution.
Regarding these features as the center, a multi-layered buffer zone was established (the multilayer buffer schematic is showed in Figure 5). The service capabilities of the impact factors were explored through the statistics of the population covered by different buffer zones. For the point-like features, such as governments, educational resources, and medical and health resources, the domain expansion search is carried out by the topological relation by using the grid cell of the point-like feature as the center, and the total population in the buffer zone is counted.
For the road network feature, buffers of 50 m, 100 m, 150 m, and 200 m are established on the district level, and the total population in the buffer zone is counted. To solve the problem of grid size mismatch between the two levels, this paper spatially aggregates the population results on the street level. The size of the merged grid is consistent with the population grid on the district level.

3.3.3. Cross-Validation Analysis

Cross-validation is a model validation technique that is used to assess how the results of a statistical analysis will generalize to an independent data set [47,66]. In this paper, the two levels of the population spatialization results are used to analyze the rationality of the multi-level method and to explain the difference. The cross-validation formula is shown in Equation (5):
R e l a t i v e   e r r o r = 2 × ( C S i S S i ) / ( C S i + S S i ) ( i = 1 , 2 . n )
where C S i is the population of each grid cell on the district level. S S i is the population of each grid cell on the street level. n is the number of grid cells. To solve the problem of grid size mismatch between the two results, this paper spatially aggregates the population grid cells on the street level, such that it is consistent with the population grid cell on the district level.

4. Results and Discussion

Wuchang District is one of the areas in downtown of Wuhan City and it is adjacent to the Yangtze River and the Han River. This district is the political center of Hubei Province and is also the place where universities and talents converge. The total area is 107.76 km2, and the center is 30°33’56”in the north latitude and 114°18’90” in the east. The urbanization rate of Wuchang District reached 96.2% in 2010, and the district has consistently maintained an urbanization rate of 100% in recent years. The population in Wuchang District has maintained a growth rate of approximately 1% since 2010, implying that the population has steadily increased. In 2015, there are 14 street-level governments, 196 educational resources, 105 medical and health facilities, and a road network length of 338.84 km.
According to the administrative divisions, Wuchang District consists of the following 14 streets Baishazhou, Huanghelou, Jiyuqiao, Liangdao, Luojiashan, Nanhu, Shouyi Road, Shuiguohu, Xujiapeng, Yangyuan, Zhonghua Road, Zhongnan Road, Ziyang, and Shidong. The geographic location of Wuchang District is shown in Figure 6.
Wuchang District is selected as an experimental area for the following reasons:
(1)
Wuchang District is located in the central urban area of Wuhan City, where the buildings are more concentrated and the types of buildings are more complicated. Therefore, the method will be scientific and universal if it has highly accurate results.
(2)
Wuchang District contains 14 streets and 195 communities, and the data and information on house buildings are adequate.
(3)
Wuhan City has conducted the Community Demographic Census since 2013. The granularity of statistical units is small, the population data sources are adequate, and the recency of the data is sufficient.
(4)
Wuchang District has made many efforts to rebuild house buildings in recent years. If we can extract the spatial distribution of residential houses accurately, and remove abandoned buildings and other types of buildings to calculate the population spatialization results accurately, the study can provide important reference values for other cities with rapid urbanization.

4.1. The Results of the Population Spatialization Method on the District Level

In the experiment, it is not satisfactory to divide residential houses into four types because the results do not meet a conditional judgment. Thus, it is helpful to reclassify the residential houses. The correlation between the street population and estimated population when using the least squares regression model is shown in Figure 7.
The Figure 7 shows a clear linear correlation between the resident population and estimated population. The fitting coefficient reaches 0.936 and the goodness of fit is 0.725, which not only satisfies the conditional judgment, but also verifies that the experimental results have good accuracy.
According to the relative proportion between the estimated population and street population, it is reasonable to use Equation (3) to correct the population count of each type of house building. The corrected coefficient results are shown in Table 3.
As shown in Table 3, the relative proportion of almost all streets is basically less than 30%, except for Nanhu. The field investigation found that the main reason for this result is that the construction of Nanhu has been developing rapidly in recent years; the areas of house buildings have been increasing significantly, while the resident occupancy rate remains relatively low. The average fitting error in the Wuchang District is only 13.03%, indicating that the use of this method on the district level is reasonable.
According to the calculations of the population spatialization method on the district level, the number of population grid cells in Wuchang District is 1300 and the total population is approximately 1.21 million. The overall accuracy of the experiment reaches 99.98%, as the actual resident population is 1.182 million. The 250-m spatial population result for Wuchang District is shown in Figure 8.
As shown in Figure 8a, the areas with a large population are mainly distributed in the northern and middle-central areas of Wuchang District, including Yangyuan, Zhongnan Road, Zhonghua Road, Huanghelou, and Liangdao. The local correlation analysis result of the population spatialization is shown in Figure 8b. Most areas of the Wuchang District do not have an obvious spatial correlation of the population, especially Xujiapeng, Jiyuqiao, Luojiashan, Baishazhou, and Nanhu, where the population is small and the residents are scattered. In regions where the correlation is obvious, the vast majority of regions satisfy the “High-High” condition, and they are concentrated in Yangyuan, Shuiguohu, Zhonghua Road, Huanghelou, Liangdao, and Zhongnan Road. The residential houses in these areas are relatively compact, and large-level residential quarters have a spatial structure of “adjacent”. Therefore, these areas mainly provide living space that meets the needs of residents in terms of regional planning.
As shown in Figure 8c and Figure 9 and Table 4, as the population increases, these factors obviously reduce the service distance, and the experimental results have better fit accuracy. Medical and health resources are most sensitive to population concentration, and, when compared with other features, these services are associated with a higher percentage of the population in the nearest buffers. With the expansion of urbanization, the areas surrounding the government are mainly replaced by commercial land, and the population coverage in different buffer areas shows that the population tends to move outward from the center of the street, but the location of governments remains highly attractive. The reasonable distribution and perfect construction of educational resources and road networks have a lesser influence on the population level. As shown by the overlay, the central region with a larger population has obvious advantages in terms of location and the strong construction of public service facilities. In the southern parts of Zhongnan Road and Shuiguohu, areas that are located on the main road of the Wuchang District, the road network is well developed and traffic is convenient. Although these areas are far from the government, they are also highly attractive to the population, leading to the development of education and medical facilities.

4.2. The Results of the Population Spatialization Method on the Street Level

According to the calculations of the population spatialization method on the street level, the number of population grid cells in Wuchang District is 28,599 and the total population is approximately 1.22 million. The overall accuracy of the experiment reached 99.97%, as the actual resident population is 1.182 million. The experimental results are shown in Figure 10.
As shown in Figure 10a, the areas with large populations in Wuchang District are mainly distributed in the central region, including parts of Luojiashan, Zhongnan Road, Huanghelou, Shouyi Road, and Shuiguohu. In contrast to the results that are presented in Section 3.1, these results show that the development of road networks in some parts of Luojiashan and Shuiguohu is relatively common, but the population is also relatively large. This result is mainly due to the proximity of colleges and universities to these areas, as educational and medical resources are abundant. Areas with a small population are also more concentrated, and are mainly distributed in Xujiapeng, Baishazhou, and Nanhu. Educational and medical and health resources in these areas are relatively scarce.
Similarly, the local correlation analysis on the results of the population spatialization is shown in Figure 10b. The figure shows no obvious aggregation in most areas of Wuchang District. In addition, the “High-High” condition and “Low-Low” condition regions reflect the strong correlation of the population distribution in some areas of Wuchang District. The differentiation between high-aggregation regions and low-aggregation regions is also obvious. When comparing the results of the spatial autocorrelation on the district and street levels, except for some streets in Luojiashan and Xujiapeng, Shouyi Road, and Nanhu, the distribution of “High-High” population aggregation areas obtained by the two methods is very similar. The “Low-Low” population clustering characteristics on the street scale are more obvious than those on the district scale.
As shown in Table 5, more than half of the regions in Wuchang District have no population. Most grid cells have less than 25 people, while the grid cells with more than 300 people account for only approximately 3% of the total grid cells. The number of grid cells shows a significant decrease as the number of people increases, indicating that living space in Wuchang District is relatively decentralized and that the number of population-concentrated areas is relatively small.
When compared with the use of the method on the district level, which shows a strong correlation between population results and these features, the use of the method on the street level can better reflect the geospatial uncertainty. Of the examined features, medical and health resources are the most sensitive to population concentration. Furthermore, as compared with the other features, these services are associated with a higher percentage of the population in the nearest buffers. The government has no apparent sensitivity to the extent of population aggregation. Distance increases as the number of people increasing, and the trend is rapidly decreasing in areas with a population of more than 1200. Educational resources and road networks have a relatively small impact on the population and cover nearly 80% of the population in the first buffer area. As shown by the overlay, most of the areas with large populations are close to governmental locations with medical and health and educational resources and well-developed road networks. Although road networks are not well developed in some parts of Luojiashan and Shuiguohu, there are universities and abundant educational and medical care resources nearby, so the population is large too. The Baishazhou and Nanhu streets are far from the urban center. These areas have poor road networks, and the population is small.

4.3. The Results of Cross-Validation Analysis

The average deviations of the population that were covered by government, educational resources, medical and health resources, and road networks in buffer zones were 7.98%, 0.91%, 3.68%, and 7.56%, respectively, and the correlation coefficient of the results that were obtained by the two methods was 0.59. This analysis shows that the population within the coverage of the impact factor is relatively small, and that the results of the two methods are highly consistent. A thematic map in the population in Wuchang District on the different levels is shown in Figure 11.
Data on differences in the population results are collected, and the results are shown in Table 6.
Table 6 shows that the population results obtained from the two levels are not significantly different.61.84% of the population difference value is between −0.4 and 0.4, and only approximately 9% of the results are less than −1 or greater than 1. Large population differences are concentrated in the marginal areas of Wuchang District, such as Xujiapeng, Zhongnan Road, Baishazhou, Shidong, and Luojiashan. When combined with the areas of residential houses, the type of residential houses, and the residents of the community, these communities that are far from the center of Wuchang District have a small number of people and large areas of residential houses. Therefore, the estimated population on the district level is larger than that on the street level. Regarding Luojiashan, the distributions of the population and the residential houses are much more concentrated than those of the other streets because of the large terrain undulations. Additionally, several communities have many residents and small areas of houses, so the estimated population on the street level is greater than that on the district level.
The following conclusions can be drawn that are based on the above analysis: the population spatialization method on the district level can better highlight the spatial distribution of the population from the macro perspective. This method focuses on the impact of different types of residential houses on population density. This method is suitable when house areas are sufficiently large and the distribution of population and types of houses are diversified. Meanwhile, the population spatialization method on the street level can better display the present situation of the spatial distribution of the population from the micro perspective. This method focuses on calculating the impact of the residential houses within a small region on the population of the community. It is suitable for areas where the type of residential houses is simple and the house areas are small.

4.4. The Evaluation of Population Fit Accuracy

This paper uses the 1-km population grid data set of China from the National Earth System Science Data Sharing Infrastructure (http://www.geodata.cn), which established multivariate statistical models for populations in 1-km pixels in 2010 in China based on the correlations between the population and land use types. Urban population density, traffic conditions, DEM, and so on were used for model correction and forty counties with township population data from eastern, western, and central of China were chosen for precision verification.
To solve the problem of grid size mismatch between the two results, this paper spatially aggregates the population results on the street level and the district level, and the size of the merged grid is consistent with the 1-km population grid data set. The number of effective population grids has been adjusted to 66 by excluding the population grid around Wuchang District and avoiding a large population error due to the lack of data on housing construction in other urban areas. The fit analysis was performed on the estimated population on the district level and on the street level. The obtained results are shown in Figure 12.
As shown in Figure 12, the goodness of fit on the street level is slightly lower than that on the district level, and the fitting coefficient on the street level is closer to 1 than that on the district level. The results on the district level have higher accuracy, but the results on the street level have a lower coefficient sensitivity. Thus, the results of the two levels are highly accurate and they have their own advantages.

5. Conclusions

This paper fully mines the Geographical National Monitoring Data to establish a multi-level population spatialization method for the different administrative levels. It also uses spatial analysis methods to explore correlations and influencing features, and finally verifies the accuracy of the fit results.
(1)
The average fitting error in Wuchang District is only 13.03%, the fitting coefficient reaches 0.936, and the overall population estimation accuracy reaches 99.98% after building reclassification. The overall population estimation accuracy is 99.97% on the street level. The results truly reflect the spatial distribution of the population on the different levels.
(2)
The spatial correlation in most areas of Wuchang District is not obvious through the spatial autocorrelation analysis results. In regions where the correlation is obvious, there is a large proportion of areas in the “High-High” condition on both levels. The distributions of the “High-High” population aggregation areas that were obtained by the two methods are highly similar, and the “Low-Low” areas on the street level are more obvious.
(3)
In most areas, geographical location and road network are the dominant features that promote population aggregation. In other areas, the availability of public service can attract population aggregation, despite less convenience in terms of traffic.
(4)
The population results that were obtained from the two levels are not significantly different; more than 60% of the population difference value is between −0.4 and 0.4. The average deviations of the population covered by different factors in the buffer zones were 7.98%, 0.91%, 3.68%, and 7.56%, and the correlation coefficient of the results obtained by the two methods was 0.59.
(5)
When comparing the accuracy of experimental results against the 1-km population grid data, the fitting coefficient is 1.324, and the goodness of fit is 0.422 on the district level, while the fitting coefficient is 1.236 and the goodness of fit is 0.300 on the street level. The results of the two levels are highly accurate and have their own advantages.
However, the experiment found some deviations from relying solely on China’s first national geoinformation survey. The extraction of building construction may consider non-residential buildings or non-demolished houses, which may impact the calculation of population density. Future research will consider taxi trajectory data to solve this problem. The population density on the street or community level can be divided through the mining of taxi trajectories and analysis of residents’ activity patterns and activity areas. Then, differences in the spatial distribute of residents can be explored and the impact of different types of residential houses and regional locations on population activities can be analyzed. Finally, a model is established to estimate the population distribution status in residential areas at a smaller scale, and to improve the accuracy of population spatialization.

Author Contributions

L.L. provided research ideas, designed experimental methods and completed the writing of the first draft and the visualization of the results. J.L. provided the theoretical basis, conducted paper revision and directed the experiment. Z.J. and L.Z. were responsible for data collection, field research and experimental results rationality assessment. P.Z. contributed analysis tools.

Funding

This research was supported by the National Natural Science Foundation of China (NSFC), Nos. 41471109 and 41501139; Fundamental Research Funds for the Central Universities of China, No. 2042016kf0066.

Acknowledgments

The China’s first national geoinformation survey and Wuhan Community Demographic Census in 2013 were supported by Wuhan Institute of Surveying and Mapping and the 1-km population grid dataset was supported by National Earth System Science Data Sharing Infrastructure, National Science & Technology Infrastructure of China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, X.; Huang, Y.; Dong, P.; Jiang, D.; Liu, H. An updating system for the gridded population database of china based on remote sensing, GIS and spatial database technologies. Sensors 2009, 9, 1128–1140. [Google Scholar] [CrossRef] [PubMed]
  2. Chen, Y.-G. The corresponding relationships between the demographic transition, urbanization process, and industrial development of China. Geogr. Res. 2010, 29, 2109–2120. [Google Scholar]
  3. Chen, X. The effect of population agglomeration on heterogeneity of regional labor productivity. Popul. Res. 2015, 39, 85–95. [Google Scholar]
  4. Kopachevsky, I.; Kostyuchenko, Y.V.; Stoyka, O. Land use drivers of population dynamics in tasks of security management and risk assessment. Int. J. Math. Eng. Manag. Sci. 2016, 1, 18–24. [Google Scholar]
  5. Liang, H.; Liu, Y. Study on spatial-temporal change and simulation of population in Beijing based on census data. Acta Geogr. Sin. 2014, 69, 1487–1495. [Google Scholar]
  6. Wang, Q.; Xiu, C.-L.; Wei, Y.; Zhang, L.-J. Evaluation of job accessibility of shenyang using gaussian based 2-step floating catchment areas method. J. Hum. Geogr. 2015, 30, 78–82. [Google Scholar]
  7. Qi, L.-L.; Zhou, S.-H.; Yan, X.-P. Endpoint attractive factors of medical facilities’ accessibility: Based on GPS floating car data in Guangzhou. Acta Geogr. Sin. 2014, 34, 580–586. [Google Scholar]
  8. Xiao, H.; Tian, H.-Y.; Zhu, P.-J.; Yu, H.-K. The dynamic simulation and forecast of urban population distribution based on the multi—Agent system. Prog. Geogr. 2010, 29, 347–354. [Google Scholar]
  9. Wu, W.; Gao, X. Population density functions of chinese cities: A review. Prog. Geogr. 2010, 29, 968–974. [Google Scholar]
  10. Liao, S.; Sun, J.L. GIS based spatialization of population census data in Qinghai-Tibet plateau. Acta Geogr. Sin. 2003, 58, 25–33. [Google Scholar]
  11. Du, G.; Yu, F.; Zhang, S.-W. Urban population distribution modeling and the pattern identifying—Taking Shenyang city as an example. Geo Inf. Sci. 2011, 12, 34–39. [Google Scholar]
  12. Kang, T.-J.; Zhang, X.-C.; Zhao, Y.; Wang, H.-Y.; Zhang, W. Agent-based urban population distribution model. Acta Geogr. Sin. 2012, 32, 790–797. [Google Scholar]
  13. Shi, T.-T.; Yang, X.-M.; Lan, R.-Q. Spatial simulation and influence factors detection of population density in North Korea. J. Geomat. Sci. Technol. 2017, 34, 79–84. [Google Scholar]
  14. Dong, C.; Liu, J.-P.; Zhao, R.; Wang, G.-X. An discussion on correlation of geographical parameter with spatial population distribution. Remote Sens. Inf. 2002, 4, 61–64. [Google Scholar]
  15. Kraus, S.P.; Senger, L.W.; Ryerson, J.M. Estimating population from photographically determined residential land use types. Remote Sens. Environ. 1974, 3, 35–42. [Google Scholar] [CrossRef]
  16. Iisaka, J.; Hegedus, E. Population estimation from Landsat imagery. Remote Sens. Environ. 1982, 12, 259–272. [Google Scholar] [CrossRef]
  17. Schneider, A.; Mertes, C.M.; Tatem, A.J.; Tan, B.; Sulla-Menashe, D.; Graves, S.J.; Patel, N.N.; Horton, J.A.; Gaughan, A.E.; Rollo, J.T.; et al. A new urban landscape in East–Southeast Asia, 2000–2010. Environ. Res. Lett. 2015, 10. [Google Scholar] [CrossRef]
  18. Lwin, K.; Murayama, Y. A GIS approach to estimation of building population for micro-spatial analysis. Trans. GIS 2009, 13, 401–414. [Google Scholar] [CrossRef]
  19. Briggs, D.J.; Gulliver, J.; Fecht, D.; Vienneau, D.M. Dasymetric modelling of small-area population distribution using land cover and light emissions data. Remote Sens. Environ. 2007, 108, 451–466. [Google Scholar] [CrossRef]
  20. Ural, S.; Hussain, E.; Shan, J. Building population mapping with aerial imagery and GIS data. Int. J. Appl. Earth Obs. Geoinform. 2011, 13, 841–852. [Google Scholar] [CrossRef]
  21. Chen, Q.; Hou, X.-Y.; Wu, L. Compared of population spatialization models based on land use data and DMSP/OLS data respectively: A case study in the efficient ecological economic of the yellow river delta. J. Hum. Geogr. 2014, 29, 94–100. [Google Scholar]
  22. Lo, C.P. Application of landsat TM data for quality of life assessment in an urban environment. Comput. Environ. Urban Syst. 1997, 21, 259–276. [Google Scholar] [CrossRef]
  23. Li, S.; Wang, L.; Dong, N. Simulation of urban small-areas population space-time distribution based on building extraction: Taking Beijing Donghuamen Subdistrict as an example. J. Geo-Inf. Sci. 2013, 15, 19–28. [Google Scholar] [CrossRef]
  24. Wang, S.; Tian, Y.; Zhou, Y.; Liu, W.; Lin, C. Fine-scale population estimation by 3D reconstruction of urban residential buildings. Sensors 2016, 16. [Google Scholar] [CrossRef] [PubMed]
  25. Flowerdew, R.; Green, M. Statistical methods for inference between incompatible zonal systems. Accuracy Spatial Databases 1989, 239–247. [Google Scholar]
  26. Novack, T.; Kux, H.; Freitas, C. Estimation of population density of census sectors using remote sensing data and spatial regression. In Geocomputation, Sustainability and Environmental Planning; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  27. Zhang, Z.-B.; Yang, Y.; Ju, C.-P.; Deng, C.-Y. The spatial evolution of hui population and its social response in Lanzhou. Acta Geogr. Sin. 2014, 34, 921–929. [Google Scholar]
  28. Wu, C.; Murray, A.T. A cokriging method for estimating population density in urban areas. Comput. Environ. Urban Syst. 2005, 29, 558–579. [Google Scholar] [CrossRef]
  29. Goodchild, M.F.; Anselin, L.; Deichmann, U. A framework for the areal interpolation of socioeconomic data. Environ. Plan. A 1993, 25, 383–397. [Google Scholar] [CrossRef]
  30. Chen, F.-Q. Application of Support Vector Machine in the Analysis of Population Data; Liaoning Normal University: Dalian, China, 2016. [Google Scholar]
  31. Stevens, F.R.; Gaughan, A.E.; Linard, C.; Tatem, A.J. Disaggregating census data for population mapping using random forests with remotely-sensed and ancillary data. PLoS ONE 2015, 10. [Google Scholar] [CrossRef] [PubMed]
  32. Holt, J.B.; Lo, C.P.; Hodler, T.W. Dasymetric estimation of population density and areal interpolation of census data. Cartogr. Geogr. Inf. Sci. 2004, 31, 103–121. [Google Scholar] [CrossRef]
  33. Wang, K.; Cai, H.; Yang, X. Multiple level spatialization of demographic data with multi-factor linear regression and geographically weighted regression models. Porg. Geogr. 2016, 35, 1494–1505. [Google Scholar]
  34. Cao, L.-Q.; Li, P.-X.; Zhang, L.-P. Urban population estimation based on the DMSP/OLS night-time satellite data—A case of Hubei province. Remote Sens. Inf. 2009, 1, 83–87. [Google Scholar]
  35. Ruktanonchai, N.W.; DeLeenheer, P.; Tatem, A.J.; Alegana, V.A.; Caughlin, T.T.; Erbach-Schoenberg, E.; Lourenço, C.; Ruktanonchai, C.W.; Smith, D.L. Identifying malaria transmission foci for elimination using human mobility data. PLoS Comput. Biol. 2016, 12. [Google Scholar] [CrossRef] [PubMed]
  36. Thakur, G.S.; Kuruganti, T.; Bobrek, M.; Killough, S.; Nutaro, J.; Liu, C.; Lu, W. Real-time urban population monitoring using pervasive sensor network. In Proceedings of the 24th ACM SIGSPATIAL International Conference on Advances in Geographic Information Systems, Burlingame, CA, USA, 31 October–3 November 2016. [Google Scholar]
  37. Bharti, N.; Lu, X.; Bengtsson, L.; Wetter, E.; Tatemc, A.J. Remotely measuring populations during a crisis by overlaying two data sources. Int. Health 2015, 7, 90–98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Ye, J.; Yang, X.H.; Jiang, D.; Urban, M.L. The grid scale effect analysis on town leveled population statistical data spatialization. J. Geo-Inf. Sci. 2010, 12, 40–46. [Google Scholar]
  39. Bhaduri, B.; Bright, E.; Coleman, P.; Urban, M.L. LandScan USA: A high-resolution geospatial and temporal modeling approach for population distribution and dynamics. Geo J. 2007, 69, 103–117. [Google Scholar] [CrossRef]
  40. Tatem, A.J. WorldPop, open data for spatial demography. Sci. Data 2017, 4. [Google Scholar] [CrossRef] [PubMed]
  41. Thomson, D.R.; Stevens, F.R.; Ruktanonchai, N.W.; Tatem, A.J.; Castro, M.C. GridSample: An R package to generate household survey primary sampling units (PSUs) from gridded population data. Int. J. Health Geogr. 2017, 16. [Google Scholar] [CrossRef] [PubMed]
  42. Thakur, G.; Sims, K.; Mao, H.; Piburn, J.; Sparks, K.; Urban, M.; Stewart, R.; Weber, E.; Bhaduri, B. Utilizing geo-located sensors and social media for studying population dynamics and land classification. In Human Dynamics Research in Smart and Connected Communities; Springer: Cham, Switzerland, 2018. [Google Scholar]
  43. Hughes, C.; Zagheni, E.; Abel, G.; Wiśniowski, A.; Sorichetta, A.; Weber, I.; Tatem, A.J. Inferring Migrations: Traditional Methods and New Approaches Based on Mobile Phone, Social Media, and Other Big Data; European Commission: Brussels, Belgium, 2016. [Google Scholar]
  44. Buckee, C.O.; Tatem, A.J.; Metcalf, C.J.E. Seasonal population movements and the surveillance and control of infectious diseases. Trends Parasitol. 2017, 33, 10–20. [Google Scholar] [CrossRef] [PubMed]
  45. Cai, Q.; Rushton, G.; Bhaduri, B.; Bright, E.; Coleman, P. Estimating small-area populations by age and sex using spatial interpolation and statistical inference methods. Trans. GIS 2006, 10, 577–598. [Google Scholar] [CrossRef]
  46. Utazi, C.E.; Thorley, J.; Alegana, V.A.; Ferrari, M.J.; Takahashi, S.; Metcalf, C.J.E.; Lessler, J.; Tatem, A.J. High resolution age-structured mapping of childhood vaccination coverage in low and middle income countries. Vaccine 2018, 36, 1583–1591. [Google Scholar] [CrossRef] [PubMed]
  47. Bosco, C.; Alegana, V.; Bird, T.C.; Pezzulo, L.; Bengtsson, A.; Sorichetta, J.; Steele, G.; Hornby, C.; Ruktanonchai, N.; Ruktanonchai, E.; et al. Exploring the high-resolution mapping of gender-disaggregated development indicators. J. R. Soc. Interface 2017, 14. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, X.; Wang, M. How polycentric is urban China and why? A case study of 318 cities. Landsc. Urban Plan. 2016, 151, 10–20. [Google Scholar] [CrossRef]
  49. Patel, N.N.; Stevens, F.R.; Huang, Z.; Gaughan, A.E.; Elyazar, I.; Tatem, A.J. Improving large area population mapping using geotweet densities. Trans. GIS 2017, 21, 317–331. [Google Scholar] [CrossRef] [PubMed]
  50. Gaughan, A.E.; Stevens, F.R.; Huang, Z.; Nieves, J.J.; Sorichetta, A.; Lai, S.; Ye, X.; Linard, C.; Hornby, G.M.; Hay, S.I.; et al. Spatiotemporal patterns of population in mainland China, 1990 to 2010. Sci. Data 2016, 3. [Google Scholar] [CrossRef] [PubMed]
  51. Yuan, J. Automatic building extraction in aerial scenes using convolutional networks. arXiv, 2016; arXiv:1602. [Google Scholar]
  52. Vijayaraj, V.; Bright, E.A.; Bhaduri, B.L. High resolution urban feature extraction for global population mapping using high performance computing. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Barcelona, Spain, 23–28 July 2007. [Google Scholar]
  53. Patel, N.N.; Angiuli, E.; Gamba, P.; Gaughan, A.; Lisini, G.; Stevens, F.R.; Tatem, A.J.; Trianni, G. Multitemporal settlement and population mapping from Landsat using Google Earth Engine. Int. J. Appl. EarthObs. Geoinform. 2015, 35, 199–208. [Google Scholar] [CrossRef] [Green Version]
  54. Mao, H.; Thakur, G.; Bhaduri, B. Exploiting mobile phone data for multi-category land use classification in Africa. In Proceedings of the 2nd ACM SIGSPATIAL Workshop on Smart Cities and Urban Analytics, Burlingame, CA, USA, 31 October–3 November 2016. [Google Scholar]
  55. Linard, C.; Alegana, V.A.; Noor, A.M.; Snow, R.W.; Tatem, A.J. A high resolution spatial population database of Somalia for disease risk mapping. Int. J. Health Geogr. 2010, 9. [Google Scholar] [CrossRef] [PubMed]
  56. Dobson, J.E.; Bright, E.A.; Coleman, P.R.; Worley, B.A. LandScan: A global population database for estimating populations at risk. Photogramm. Eng. Remote Sens. 2000, 66, 849–857. [Google Scholar]
  57. Nieves, J.J.; Stevens, F.R.; Gaughan, A.E.; Linard, C.; Sorichetta, A.; Hornby, G.; Patel, N.N.; Tatem, A.J. Examining the correlates and drivers of human population distributions across low-and middle-income countries. J. R. Soc. Interface 2017, 14. [Google Scholar] [CrossRef] [PubMed]
  58. Alegana, V.A.; Wright, J.; Pezzulo, C.; Tatem, A.J.; Atkinson, P.M. Treatment-seeking behaviour in low-and middle-income countries estimated using a Bayesian model. BMC Med. Res. Methodol. 2017, 17. [Google Scholar] [CrossRef] [PubMed]
  59. Wei, H.K. Administrative hierarchy and growth of city scale in China. Urban Environ. Stud. 2014, 1, 4–17. [Google Scholar]
  60. Tian, S.S. An Introduction to the Administrative Division of China; Peking University Press: Beijing, China, 2005. [Google Scholar]
  61. Zhou, Y.-P.; Shi, F.; Xu, J.-G.; Xu, X.-Y. Study on the planning and design of residential area. Modern UrbanRes. 2016, 3, 2–8. [Google Scholar]
  62. Carè, A.; Garatti, S.; Campi, M.C. A coverage theory for least squares. J. R. Stat. Soc. 2017, 79, 1367–1389. [Google Scholar] [CrossRef]
  63. Zhang, Y.; Liu, L. Understanding temporal pattern of human activities using Temporal Areas of Interest. Appl. Geogr. 2018, 94, 95–106. [Google Scholar] [CrossRef]
  64. Engen, S.; Lee, A.M.; Sæther, B.E. Spatial distribution and optimal harvesting of an age-structured population in a fluctuating environment. Math. Biosci. 2018, 296, 36–44. [Google Scholar] [CrossRef] [PubMed]
  65. Anselin, L. Local indicators of spatial association—LISA. Geogr. Anal. 1995, 27, 93–115. [Google Scholar] [CrossRef]
  66. Xu, J.; Li, A.; Li, D.; Liu, Y.; Du, Y.; Pei, T.; Ma, T.; Zhou, C. Difference of urban development in China from the perspective of passenger transport around spring festival. Appl. Geogr. 2017, 87, 85–96. [Google Scholar] [CrossRef]
Figure 1. Schematic Diagram of Data Preprocessing.
Figure 1. Schematic Diagram of Data Preprocessing.
Sensors 18 02558 g001
Figure 2. Workflow of Multi-level Population Spatialization Method and Verification.
Figure 2. Workflow of Multi-level Population Spatialization Method and Verification.
Sensors 18 02558 g002
Figure 3. Workflow of the Population Spatialization Method on the District Level.
Figure 3. Workflow of the Population Spatialization Method on the District Level.
Sensors 18 02558 g003
Figure 4. Workflow of the Population Spatialization Method on the Street Level.
Figure 4. Workflow of the Population Spatialization Method on the Street Level.
Sensors 18 02558 g004
Figure 5. Multilayer Buffer Schematic.
Figure 5. Multilayer Buffer Schematic.
Sensors 18 02558 g005
Figure 6. The Geographic Location of Wuchang District in Wuhan, China.
Figure 6. The Geographic Location of Wuchang District in Wuhan, China.
Sensors 18 02558 g006
Figure 7. The Fit of the Street Population and Estimated Population.
Figure 7. The Fit of the Street Population and Estimated Population.
Sensors 18 02558 g007
Figure 8. The Experimental Results and Analysis on the District Level ((a) Wuchang 250-m Population Grid; (b) The Result of Local Moran’s I; and, (c) Spatial Overlay Analysis).
Figure 8. The Experimental Results and Analysis on the District Level ((a) Wuchang 250-m Population Grid; (b) The Result of Local Moran’s I; and, (c) Spatial Overlay Analysis).
Sensors 18 02558 g008
Figure 9. The Fit of Population and the Average Nearest Distance to Features ((a) result on the district level; and, (b) result on the street level).
Figure 9. The Fit of Population and the Average Nearest Distance to Features ((a) result on the district level; and, (b) result on the street level).
Sensors 18 02558 g009aSensors 18 02558 g009b
Figure 10. The Experimental Results and Analysis on the Street Level ((a) Wuchang 50-m Population Grid; (b) The Result of Local Moran’s I; and, (c) Spatial Overlay Analysis).
Figure 10. The Experimental Results and Analysis on the Street Level ((a) Wuchang 50-m Population Grid; (b) The Result of Local Moran’s I; and, (c) Spatial Overlay Analysis).
Sensors 18 02558 g010aSensors 18 02558 g010b
Figure 11. Population Difference According to Cross-validation Analysis.
Figure 11. Population Difference According to Cross-validation Analysis.
Sensors 18 02558 g011
Figure 12. The Correlation between the Grid Population Data set and the Estimated Population.
Figure 12. The Correlation between the Grid Population Data set and the Estimated Population.
Sensors 18 02558 g012
Table 1. Detailed data usage list.
Table 1. Detailed data usage list.
Data NameData SourcesData InterpretationData Format
House buildingThe LCA layer in the Land Cover Classification DataHouse buildings refer to the urban and rural areas of residential areas of housing construction, according to the attributes can be divided into 5 typesShapefile
Urban integrated functional Units (point/area)The BUCP and BUCA layers in Social Geographical UnitsThe Space unit divided by function and ownership within the urban residential areas. Including residential quarters and non-residential quarters (industrial and mining enterprises, institutions and companies). The difference between BUCA layer and BUCP layer is that BUCA layer contains surface vector features while BUCP contains point vector elementsShapefile
District and street level administrative divisionsThe BOUA5 and BOUA6 layers in Social Geographical UnitsVector data with region and district and street boundaryShapefile
Community resident population and Community level administrative divisionWuhan Community Demographic CensusCommunity vector boundary with community resident populationShapefile
Street resident populationThe Sixth National Population CensusAll streets resident population data in WuhanExcel
District resident populationWuhan Statistical YearbookAll counties resident population in WuhanExcel
Table 2. Detailed Description of House Buildings.
Table 2. Detailed Description of House Buildings.
multi-floor buildingsMore than 10 m in height or over four floors, construction area is more than 1600 m2. Mostly in densely populated areas in the central city
multi-floor independent buildingsMore than 10 m in height or over four floors, construction area is more than 200 m2. Most buildings are scattered with low population density
low-floor buildingsLess than 10 m in height or four floors below, construction area is more than 1600 m2. Mostly high-grade residential quarters or planning township gathering area
low-floor independent buildingsLess than 10 m in height or four floors below, construction area is more than 200 m2. Most buildings are in rural areas where the economy is lagging behind and there are no plans for housing construction.
abandoned house buildingAbandoned buildings after the migration
Table 3. Statistics of Corrected Coefficient Divided by Street.
Table 3. Statistics of Corrected Coefficient Divided by Street.
NameStreet Population (P)Multi-Floor + Multi-Floor IndependentLow-Floor + Low-Floor IndependentEstimates Population (P)Fitting Error (%)Coefficient 1Coefficient 2
Baishazhou78,676535,235792,363.499,206.82−26.100.0875270.040169
Huanglelou52,713282,799307,857.446,805.2411.210.1242990.057044
Jiyuqiao61,329495,492.2142,078.461,882.9−0.900.109380.050198
Liangdao64,008393,565.1458,598.366,665.45−4.150.1059680.048632
Luojiashan62,574564,579.4228,263.973,873.29−18.050.0934870.042904
Nanhu53,159839,957.991,927.397,360.69−83.140.0602610.027655
Shidong461811,322.3974,066.455001.169−8.300.1019120.04677
Shouyi Road69,872511,765225,943.767,926.762.780.1135290.052102
Shuiguohu172,0071,397,723263,553.6167,613.12.550.1132610.051979
Xujiapeng122,1291,171,513557,920.3157,556.8−29.010.0855510.039262
Yangyuan109,245840,912.3357,030.3110,893.7−1.510.1087270.049898
Zhonghua Road44,693276,456.8370,415.949,273.92−10.250.1001070.045942
Zhongnan Road234,4791,216,299550,660.3162,13230.850.1596170.073253
Ziyang52,770253,579.2242,408.340,265.2523.700.1446440.066381
Table 4. Population Coverage by Multiple Buffers.
Table 4. Population Coverage by Multiple Buffers.
Coverage Degree (%)District LevelStreet Level
GovernmentEducational ResourcesMedical and Health ResourcesRoad NetworkGovernmentEducational ResourcesMedical and Health ResourcesRoad Network
Seed2.1720.6613.28--2.0720.1111.68--
First (50 m) buffer16.1982.4063.9992.6715.5082.9662.8276.84
Second (100 m) buffer34.2196.8187.7095.3232.0095.9585.3590.66
Third (150 m) buffer51.6999.2295.1296.9746.1599.2193.6892.63
Forth (200 m) buffer67.1599.7097.7297.6857.8510097.4293.96
Table 5. Statistics of the 50-m Grid Population is obtained through counting the population spatialization results on the street level.
Table 5. Statistics of the 50-m Grid Population is obtained through counting the population spatialization results on the street level.
PopulationCount (N)Proportion (%)
015,71154.94
(0,25]353012.34
(25,50]23428.19
(50,100]319211.16
(100,150]16375.72
(150,200]7442.60
(200,300]7452.60
(300,450]4411.54
(450,700]1400.49
(700,1200]1040.36
>1200 130.06
Total28,599100
Table 6. Statistics of Population Differences on the Different Levels.
Table 6. Statistics of Population Differences on the Different Levels.
Difference Count (P)Proportion (%)
<−1.0624.77
[−1.0–−0.8)332.54
[−0.8–−0.4)735.61
[−0.4–0)13310.23
041832.15
(0–0.4]25319.46
(0.4–0.8]20916.08
(0.8–1.0]685.23
>1.0513.93
Total13001

Share and Cite

MDPI and ACS Style

Li, L.; Li, J.; Jiang, Z.; Zhao, L.; Zhao, P. Methods of Population Spatialization Based on the Classification Information of Buildings from China’s First National Geoinformation Survey in Urban Area: A Case Study of Wuchang District, Wuhan City, China. Sensors 2018, 18, 2558. https://doi.org/10.3390/s18082558

AMA Style

Li L, Li J, Jiang Z, Zhao L, Zhao P. Methods of Population Spatialization Based on the Classification Information of Buildings from China’s First National Geoinformation Survey in Urban Area: A Case Study of Wuchang District, Wuhan City, China. Sensors. 2018; 18(8):2558. https://doi.org/10.3390/s18082558

Chicago/Turabian Style

Li, Linze, Jiansong Li, Zilong Jiang, Lingli Zhao, and Pengcheng Zhao. 2018. "Methods of Population Spatialization Based on the Classification Information of Buildings from China’s First National Geoinformation Survey in Urban Area: A Case Study of Wuchang District, Wuhan City, China" Sensors 18, no. 8: 2558. https://doi.org/10.3390/s18082558

APA Style

Li, L., Li, J., Jiang, Z., Zhao, L., & Zhao, P. (2018). Methods of Population Spatialization Based on the Classification Information of Buildings from China’s First National Geoinformation Survey in Urban Area: A Case Study of Wuchang District, Wuhan City, China. Sensors, 18(8), 2558. https://doi.org/10.3390/s18082558

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop