Spatiotemporal Evolution of Multiscale Urbanization Level In The Beijing-Tianjin-Hebei Region Using The Integration Of DMSP/OLS And NPP/VIIRS Night Light Datasets

: The level of urbanization is a key factor in urban development. In this study, to better characterize the level of urbanization, the panel entropy weight method is used to weight the factors of population, industry, and area to construct a composite indicator of urbanization. A panel regression between this composite index and the average night light values after fusion shows a strong correlation. An accuracy test indicates that the estimated value of fused average light as calculated by the urbanization level estimation model that adequately represents the urbanization level. On this basis, night light data is corrected for zero error on the pixel scale, and spatiotemporal evolution analyses are performed on the city and county scales. The standard deviation ellipse method is used to find that the spatial distribution pattern of the Beijing-Tianjin-Hebei urbanization level from 1995 to 2018 radiates and spreads to the northeast, with Beijing-Tianjin as the center. The spatial pattern shows a contracting trend that is strengthening year by year. Slope analyses show that areas with rapid urbanization growth are mainly concentrated in Beijing and Tianjin. The urbanization development speed of most counties in the Hebei Province is at a low level.


Introduction
Urbanization refers to the spatial and demographic process whereby towns and cities gain importance with the concentration of population in a given economy and society [1]. Throughout history, urbanization has always been a key factor in the development process [2]. According to China's National New Urbanization Plan (2014-2020), released in 2014, the urbanization rate of permanent residents in China will reach about 60% by 2020 [3]. It is expected that, under the Thirteenth Five-Year Plan, China will continue to promote a new type of people-centered urbanization [4], which has hitherto yielded fruitful results. According to the National Bureau of Statistics, at the end of 2019 there were 848 million permanent urban residents in China, 60.6 percent of the total population, 4.5 percentage points higher than at the end of 2015, reaching the planned target one year ahead of schedule. Currently, China's economic development has entered a new normal of structural adjustment, in which urbanization is considered an important engine of China's economic growth [5] The Beijing-Tianjin-Hebei (BTH) region includes Beijing, Tianjin, and 11 prefecturelevel cities in the Hebei Province. The total area is about 216,000 square kilometers, accounting for 2.3% of China's land area. In 2018, the permanent population was 118 million, accounting for 8% of China's total population. The BTH has the most dynamic economy, the highest degree of openness, and the largest number of foreign-born population [6]. At the same time, in this region, urbanization and regional development are extremely incompatible-including the capital Beijing, the first-tier city Tianjin, and the relatively backward cities in the Hebei Province. In order to achieve coordinated economic and social development in the BTH region and to achieve a balanced increase in the urbanization level, it is necessary to formulate more targeted and refined policies. This indicates that it is necessary to understand the temporal and spatial evolution characteristics of the urbanization level [3] of the BTH region on a fine scale (at least at the county level), analyze its temporal and spatial differences, and formulate relevant policies according to local conditions to achieve further urbanization.
Lin et al. [17] and Zhang et al. [18] have previously analyzed land urbanization, and Chen et al. [4]and Ju et al. [19] have conducted policy research on the direction of China's urbanization in terms of development rate.
Regarding the level of urbanization, Dong et al. [20] established a comprehensive evaluation model using GIS technology to study the scale of urban development in the BTH region. Wang et al. [21] reconciled the differences between official census data and other survey data to recalculate China's urban population and employment data. These results indicate that official statistics underestimate China's rate of urbanization. To predict the future state of urbanization in China, Gu et al. [22] proposed a system dynamics (SD) model for China's urbanization and conducted stock-flow tests and sensitivity analyses using real data from 1998 to 2013 to establish the validity of the model and simulation effects. Jin et al. [23] applied stochastic frontier analysis to evaluate the urbanization efficiency of the Yangtze River Economic Belt from 2004 to 2015 and analyzed the spatial correlation characteristics in these data.
The above investigations have produced no consensus on the authenticity and reliability of urbanization as measured by different indicators in China, which may affect our understanding of the process of urbanization and its impact on the environment [24]. Night light (NTL) data acquired by satellite has been used as an alternative measure to a series of socio-economic variables to investigate human activities and the development of urbanization [25]. A large number of studies have used DMSP/OLS (Defense Meteorological Satellite Program's Operational Linescan System) or NPP/VIIRS (National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite) data time series to monitor population density [6,26,27], urban growth [28], and carbon emissions [29,30].
Continuous satellite observation of NTL provides consistent and effective alternative measures for population and socio-economic dynamics in urbanization [31]. In applied study on urbanization level, Zhang et al. [32] used NTL data to draw dynamic maps of urbanization on the regional and global scales, employing time characteristics to determine differences in the urbanization trajectory. Yi et al. [33] used the DMSP/OLS NTL data from 1992 to 2010 to construct an urban light index to analyze and evaluate urbanization and urban expansion. To provide a reliable and accurate measurement of China's comprehensive urbanization level,Gao et al. [24] and Fu et al. [25] used DMSP/OLS data to establish a composite luminous index, analyzing dynamic changes in China's urbanization level from 1992 to 2012. Yang et al. [3] extracted time series NTL trajectories for each 1 km × 1 km area from 1992 to 2013 to characterize urbanization dynamics, establishing and comparing the distribution and spatial patterns of NTL trajectories for different administrative divisions in China.
Previous studies have found that DMSP/OLS data are subject to pixel saturation in urban centers, lack comparability among pixels, and expire in 2013 [30], which makes it impossible to measure the differences among urban centers in different provinces or changes in urbanization after the integration of the BTH region in 2014. NPP/VIIRS NTL data have the advantages of DMSP/OLS data [34] and evade the problems listed for DMSP/OLS data, enabling an expansion of research direction and application fields for NTL data [35]. Many scholars have used DMSP/OLS [36] or NPP/VIIRS [37] data for research, but relatively few studies have used integrated NTL datasets to determine urbanization.
Most studies of composite indicators representing urbanization level have given equal weight to each urbanization factor [24,25], but this is not adequate to the importance for each factor in the construction of composite indicators. In this paper, using panel data, the entropy method is used for weighting to increase the differential importance of each factor in the constructed composite urbanization index.
In previous research, an integration of the DMSP/OLS and NPP/VIIRS NTL datasets [29,34] from 1995 to 2018 has been utilized; they both ensured light efficiency and timeliness. Using studies of the urbanization level [24] carbon emissions [38], and other data, this paper used light data and a composite urbanization index to establish an estimation model for urbanization level with a fixed effect model, correcting light data at the pixel scale. To verify the model accuracy, a precision test of the estimation model was carried out on the city scale.
The spatio-temporal evolution of the urbanization level of the BTH region was analyzed at the pixel scale and the city scale. At the city scale, the standard deviation ellipse was utilized to describe the spatial distribution of urbanization.
With reference to the work of [38], corrected lighting data were retrieved at the county-level scale using the urbanization-level estimation model. Slope analysis was used to establish the spatiotemporal evolution and the growth rate characteristics of the urbanization level at the county level for the BTH region from 1995 to 2018 to provide a reference to assess the BTH urbanization policy.

Study Area
The BTH region was selected as the study area. Figure 1 shows the geographical location of Beijing, Tianjin, and Hebei in China.

Data Sources
The data sources include NTL data: DMSP/OLS and NPP/VIIRS. At the pixel scale, the spatial resolution was 1km × 1km. Lambert Azimuthal Equal Area Projection was employed in ArcGIS. Statistical data: population, industry and area, in terms of city scale (including Beijing, Tianjin, Shijiazhuang, and other 13 cities). Descriptive information on the data sources are given in Table 1.  (2) Image reprojection, resampling, and cropping were performed on the two kinds of NTL data.
(3) In order to solve the problem of pixel saturation of DMSP/OLS data and lack of comparability of data pixels, the DMSP/OLS NTL data was mutually corrected and fused.
The invariant target area was selected as Jixi City in Heilongjiang Province for mutual image correction, which showed relatively stable economic and social development from 1992 to 2013 [39].
An annual fusion of NTL data obtained from different satellites during the same year was carried out to resolve the problem of varying NTL images produced by different satellites.
(4) Compared with the DMSP/OLS data, the NPP/VIIRS sensor can detect weaker lights with higher temporal and spatial resolution, but it does not remove sporadic weak lights such as flares, which will cause higher image noise. In order to reduce the influence of noise, this paper takes the stable bright element area of the 2016 annual synthetic data as the invariable area and makes stability corrections to the 2012-2016 data.
(5) Time series correction of two NTL data. The basic assumption of the time series correction was consistent with the reality of China's rapid economic development, that is, the pixel DN value of the NTL data from the previous year should not be greater than that in the following year.
2) Data fusion processing Because the DMSP/OLS data covered 1992-2013, and the NPP/VIIRS data covered time since 2012, it was necessary to merge the two.
The fusion of the DMSP/OLS and NPP/VIIRS data was done following the literature [38], and county-scale DMSP/OLS data and NPP/VIIRS data for 2012 and 2013 were selected for regression fitting. A quadratic regression model was selected, and the fitting formula was given as shown in Equation (1). The R 2 of 0.8354 indicated a strong correlation between the two datasets.  In the ArcGIS software environment, the image was cut with a vector map mask on the BTH region to obtain NTL image data for the corresponding year. To avoid interference from the regional area, the average NTL values were selected for subsequent estimation of the urbanization level.

Preprocessing of Statistical Data
For population, industry, and area statistics in Table 1, the following calculations were performed: The three indicators were standardized by min-max: where i x is the original statistical data, which i factor is the standardized data of the three urbanization factors,

Panel Entropy Weight Method
The entropy-weighting method is an objective method that has come into widespread use in recent years [40,41]. It is a method to determine the weight of indicators by using the amount of information provided by the entropy of various indicators. Entropy weighting method avoids the interference of human factors in the weight of each rating index and make the evaluation result more realistic. It overcomes the problem that the index weighting process is greatly affected by human factors in the current evaluation method. By calculating the entropy value of each index, it measures the size of the index information, to ensure that the established index can reflect most of the original information.
The smaller the entropy value of the factor, the greater the degree of variation in the factor value, the greater the amount of information provided, and the greater the weighting of the factor.
Assuming that the element number in the study area is p ,the number of factors used is q , then the factor matrix is p q  The entropy of the ith factor is defined as: where, i S is the entropy of the ith factor , and p is the number of elements, ij y is defined as: where ij M refers to the standardized factor value of the ith factor and the jth evaluation object. The entropy weight of the ith factor is defined as: where, i  is the entropy weight of the ith factor , 0 q is the number of evaluation factors.

Construction of Composite Urbanization Index
The urbanization composite index is calculated according to the following formula: x is the population, industry, and area statis-

Urbanization Level Estimation Model
The larger the DN value of the integrated NTL data, the larger the urbanization level. That is, the DN value has a close linear correlation to urbanization level at the pixel level. This study assumes that the correlation between DN and the municipal composite urbanization index is consistent at the pixel level.
Based on the integrated NTL datasets, the estimation model for urbanization level is established as follows: where p U represents composite urbanization index at p element, p DN represents DN value at p element, a and b were the parameters of the linear regression equation. Relying on the accuracy of urbanization level estimation model, the zero-error method of municipal composite urbanization index was adopted to correct the pixel-scale data.

Standard Deviational Ellipse
The standard deviation elliptical (SDE) method is a classical method of analyzing the directional characteristics of spatial distribution. The SDE has been used to analyze the spatial distribution of various phenomena in urban studies [42]. The existing research shows that the operation of SDE is simple and easy to be implemented in ArcGIS. More importantly, this technique is effective in exploring the spatial distribution of geographic events [43]. It is a good choice to employ SDE to study the spatial pattern of urbanization level of BTH in different periods.
The standard deviation ellipse can be utilized to explore the centrality, expansion, orientation, and spatial morphology of the spatial distribution of urbanization from the global and spatial perspectives. The specific calculation is as follows.
The form of the standard deviation ellipse is as follows: where i x and i y are the coordinates of elements, X and Y are the average central coordinates of the elements, and n is the total number of elements.
The calculation method of rotation angle is as follows: where i x and i y are the difference between the mean center and coordinates of XY .
The standard deviations of the X and Y axes are as follows: The SDE method quantitatively represents the spatial distribution characteristics of the urbanization level across the spatial distribution range of the ellipse and the basic parameters of the center, long axis, short axis, and azimuth. The distribution range of the standard deviation ellipse gives the main range of the spatial distribution of the urbanization level.

Slope Analysis
Slope analysis [38] is to analyze the linear tendency of each pixel by using unitary linear regression model, which is widely used to analyze the year-over-year change of a certain variable. This study calculated the change in slope of the urbanization level from 1995 to 2018 by establishing an estimation model of the composite urbanization index. The changing trends of the urbanization level are analyzed on the county scale. The slope value was estimated using the least squares method, following the below formula: where t represents the total number of years,

Results of Estimated Urbanization Level Model
The panel data for the three selected urbanization level factors were established using the entropy weight method, and the entropy weights were 0.1772, 0.3119, and 0.5109, respectively. According to Equation (9) in 3.4.2, the composite index of the urbanization level was obtained: where, i  is the population, industry, and area statistics.
The weighted composite urbanization index is shown in Table 2. The composite urbanization index and the mean values of the fused NTL data were used to establish the fixed-effects model between the indexes, as shown in Table 3.  Table 4.  Table 5 shows the relative errors and the average relative errors of 13 cities in the estimation model across four years. Among the relative errors, only the estimate for Beijing in 1995 exceeded 25%, indicating an overestimated urbanization level for Beijing in 1995. The average relative errors for the 13 cities were all less than 25%, indicating a good estimation accuracy for the model. The overall city precision percentage is shown in Table 6. There were 12, 13, 13, and 13 cities with a relative error of less than 25%, accounting for 92.3%, 100%, 100%, and 100% of the total precision in 1995, 2003, 2010, and 2016 respectively. The average estimation accuracy of more than 98% was relatively high, and the precision of the urbanization level estimation model on the municipal level met the requirements. Figure 3 plots that, from 1995 to 2018, the overall urbanization level of BTH was strengthened. Areas with high levels of urbanization were obviously concentrated in the city centers of BTH. Areas with low levels of urbanization were distributed across in Baoding, Chengde, and Zhangjiakou and were concentrated in the northwestern part of the BTH region in an inverted "L" shape. After this fusion, the characteristics of NTL data distribution in the BTH region were more consistent with the differences in the urbanization process.

Spatiotemporal Dynamics on a Municipal Scale
As shows in Figure 4, from 1995 to 2003, the urbanization level of Langfang City rose from a higher to the highest level, and that of Tangshan City and Handan City rose from a medium to a higher level. The changes in Langfang and Tangshan reflected the way that Beijing and Tianjin drove the development of neighboring provinces and cities. From 2003 to 2010, overall urbanization did not significantly change. From 2010 to 2018, the urbanization of Xingtai City rose from a lower to a medium level. Overall, from 1995 to 2018, the level of urbanization in BTH grew. High-level urbanization areas throughout this period were always concentrated near Beijing and Tianjin (Langfang City and Tangshan City were added in 2003), while low-level areas were always concentrated in Chengde City and Zhangjiakou City in the Hebei Province. Figure 5 shows the spatiotemporal evolution of the standard deviation ellipse from 1995 to 2018.  Table 7 presents changes in various parameters in the standard deviation ellipse from 1995 to 2018. Among them, Center X and Center Y represent the center points of the ellipse. XstdDist and YstdDist represent the length of the X axis and the length of the Y axis. Rotation represents the direction angle of the ellipse. In reference to spatial rotation changes, the rotation value rose from 27.7723 in 1995 to 27. The magnitude of the decline indicates that areas of high-level urbanization are increasingly concentrated in the main directions around Beijing and Tianjin.
From 1995 to 2018, the ellipse-flattening rate increased year-over-year, such that the spatial distribution pattern of the urbanization level of BTH still radiated and diffused its influence from Beijing-Tianjin to the periphery, especially toward the northeast. However, the spread of this influence shows tightening trend, and it is strengthening year-over-year. From 1995 to 2018, the areas with the lowest level of urbanization were concentrated in Zhangjiakou City, Chengde City, and Qinglong Manchu Autonomous Region in Qinhuangdao City, as well as counties in the northwestern part of Baoding City. They were distributed in an L shape on the northwestern edge of BTH.

Spatiotemporal Dynamics from County Scale
The spatiotemporal development of the urbanization in each city presents the characteristics of point-piece-flat evolution. Although the level of urbanization in most regions remains below the middle level, compared with 1995, the scope has greatly expanded, covering almost the entirety of the central and southern regions.
(2) Number of counties above the medium level As shown in Figure 7, from 1995 to 2018, the number of highest-level urbanization counties has been increasing in BTH. The reason for the instability of the number of higher-level and medium-level counties is that the urbanization level of some counties at this level has increased (rather than decreased).  As shown in Figure 8, from 1995 to 2018, the cities with a significant increase in the number of counties above the medium level were Shijiazhuang City, Tangshan City, Handan City, and Langfang City. In Baoding City and Qinhuangdao City, the number of counties above the medium level rose from one in 1995 to three in 2018. The number of counties at this level in Chengde City and Zhangjiakou City rose from zero in 1995 to one in 2018. The number of counties in Xingtai City and Cangzhou City remained one, while there were zero in Hengshui City.
It is worth noting that although the number of counties above the medium level in Zhangjiakou City and Chengde City featured an increasing trend from 1995 to 2018, few counties were above the medium level. Similarly, in Xingtai City, Cangzhou City, and Hengshui City, the overall level of urbanization was low.    In Figure 10, it can be seen that from 1995 to 2018, zero counties had the lowest-level urbanization in Langfang City. In Tangshan City, Handan City, Xingtai City, Cangzhou City, and Hengshui City the number of counties with the lowest level of urbanization fell to zero between 1995 and 2018. In Shijiazhuang City and Qinhuangdao City, the number fell to one. In addition, Baoding City, Zhangjiakou City, and Chengde City had a large number of counties with the lowest urbanization level in 2018, especially Zhangjiakou City and Chengde City.  To further explore the growth rate in urbanization for each county in the BTH region from 1995 to 2018, a slope analysis of the estimated urbanization level for 24 consecutive years from 1995 to 2018 was conducted. The slope results were divided into five groups, using the natural fracture method in ArcGIS. The specific groupings are given in Table 8.  Following the grouping in Table 8, the slopes are shown in Figure 11. As shown in Figure 11, The rapid-growth regions are concentrated in the city centers of Beijing, Tianjin City, and certain prefecture-level cities in the Hebei Province (Yuhua District in Shijiazhuang City, Lunan District in Tangshan City, Beidaihe District in Qinhuangdao City, Congtai District in Handan City, and Nancheng District in Baoding City). In general, the urbanization level of most counties in the Hebei Province grew only slowly. Figure 12 shows the proportion of counties corresponding to each growth level.

Medium-growth
Relatively-slow-growth Slow-growth

Medium-growth
Relatively-slow-growth Slow-growth (c) Hebei Province  Figure 12 shows that more than half of districts in Beijing and Tianjin City had rapid growth and relatively rapid growth.
By contrast, the counties with rapid growth and relatively rapid growth only accounted for 13% of the Hebei province, and 77% of the counties in the urbanization growth rate were below medium growth.
In Figure 13, the slow growth and relatively slow growth counties were concentrated in Xingtai City, Baoding City, Zhangjiakou City, Chengde City, Cangzhou City, and Hengshui City, especially in Zhangjiakou City and Chengde City.

Medium-growth
Relatively-slow-growth Slow-growth For Shijiazhuang City, Qinhuangdao City, and Handan City, different growth rates correspond to different numbers of counties, indicating that the urbanization development of these three cities was relatively unbalanced.
For Tangshan City and Langfang City, the number of counties was concentrated in relatively rapid growth, medium growth, and relatively slow growth, and the quantity gap among these was small, indicating that their urbanization development was relatively balanced.

Analysis and Discussion
Unlike the equal weighting in existing research [24,25], the composite urbanization index constructed by the panel entropy weighting can reflect the importance of each factor. On this basis, the urbanization level estimation model constructed had good results and high accuracy (Table 6). Limited by the availability of statistical indicators, the urbanization level estimation model can only perform accuracy testing on a city scale. However, as China's statistical information improves, data at the county scale or smaller will continue to develop, allowing better scale model testing.
With this model, zero error correction on the pixel scale can be carried out on the fused NTL data. The corrected light map can represent the spatiotemporal evolution analysis of the urbanization level on the city or county scale. Fused NTL data can overcome the shortcomings of a single source such as short timelines and through zero-error correction at the pixel scale, so NTL data can be flexibly applied to the estimation of urbanization levels.
From 1995 to 2018, areas with high levels of urbanization in BTH were concentrated in Beijing, Tianjin, and the city centers of prefecture-level cities in the Hebei Province (Figures 3 and 4). Urban areas with low levels of urbanization were concentrated in the northwestern part of the Hebei Province, distributed in an "L" shape on the northwestern edge of the BTH region, and these could be areas for subsequent urbanization work. Using standard deviation ellipse analysis ( Figure 5, Table 7), it was found that the influence of Beijing and Tianjin on their surrounding areas is concentrated in the northeast, such as in Langfang and Tangshan, and the scope of their influence was shrinking year-over-year. This indicates that the scope of influence of Beijing and Tianjin's ability to promote the development of the surrounding cities has been shrinking year-over-year. The number of counties with the highest level of urbanization in the three cities of BTH was increasing year-over-year (Figure 7), and the number of counties at the lowest level was decreasing

Slow-growth
Hengshui year-over-year (Figure 9). In general, urbanization was constantly improving in BTH.
Combined with the results of slope analysis (Figure 12), Beijing and Tianjin were generally at higher levels of urbanization development, and more than 50% of the districts and counties were fast and relatively fast. The number of counties in the Hebei Province is decreasing with the increase in level, showing a pyramidal distribution. High-level counties in the Hebei Province were mainly concentrated in Shijiazhuang City, Tangshan, Handan City, and Langfang City (Figure 8). Combined with a slope analysis (Figure 13a,b,d,j), the urbanization of Shijiazhuang and Handan was relatively unbalanced, while that of Tangshan and Langfang was relatively balanced.
Xingtai City, Cangzhou City, and Hengshui City featured a small number of counties above the medium level of development (Figure 8), and the number of counties with the lowest level of development was also small (Figure 10), indicating that the county level data were mostly concentrated in the lower and medium levels. Combined with slope analysis (Figure 13e,i,k), the urbanization development rate of these three cities was slow, reflecting a lack of momentum for urbanization development in these three cities, and targeted industries and policies are urgently needed to drive the development of urbanization.
Low-level counties in the Hebei Province were concentrated to the northwest of Baoding City, Zhangjiakou City, and Chengde City ( Figure 10). Combined with slope analysis (Figure 13f,g,h), Baoding City, Chengde City and Zhangjiakou City showed a low level of urbanization and low growth rate overall, which can be a target for future development.

Conclusions and Suggestions
(1) The fixed-effects model showed that the composite urbanization index thus produced had a strong correlation with the fused NTL data. The urbanization level estimation model constructed had good results and high accuracy. The empirical results showed that after correction, it is appropriate to use night light data to indicate the level of urbanization. An effective attempt had been made for subsequent urban development research.
(2) All the results indicated that the urbanization levels of Beijing City and Tianjin City are both the highest and the fastest growing. The spatiotemporal results at the county level showed that the urbanization levels of Shijiazhuang City and Handan City are the highest in the Hebei Province. For these areas with high levels, the focus can be on how the quality of urbanization can be improved and how modern city construction can be carried out. It is necessary to give full play to the role of high-level urbanization in the surrounding areas and to promote the development of neighboring areas. However, slope analysis showed that the urbanization development of Shijiazhuang City and Handan City is not balanced. In the follow-up development, it is necessary to pay attention to their own balanced development and enable the fast-growing districts (counties) to drive the neighboring districts (counties) to achieve the coordinated development of the entire city.
(3) Both the results at city and county scales indicated that the urbanization level of Chengde City and Zhangjiakou City is at the lowest level. At the county scale a large number of counties where Baoding City was at the lowest level of urbanization. These areas are largely marginal areas of the BTH region. In the future, it is necessary to establish the dominant position of local cities, improve the level of infrastructure, optimize the industrial layout, enhance the attractiveness of the population, and improve the overall improvement of regional urbanization. (4) Although the spatiotemporal results did not show that Xingtai City, Cangzhou City and Hengshui City are at low urbanization level, the slope results showed that the urbanization growth of these three cities is slow. The first task is to solve how to improve rapid urbanization, how to stimulate the vitality of urban development and narrow the gap with other cities. For example, technological innovation and industrial upgrading should be used to drive economic growth, stimulate domestic demand and then increase exports, increase labor rates, and thereby increase the level of urbanization. (5) The spatiotemporal results at the county scale showed that the number of counties with the lowest levels of urbanization in Langfang City and Tangshan City are both zero. Slope analysis showed that the development of the two cities is relatively balanced. For Langfang City and Tangshan City, the urban construction plan should be reasonably planned to ensure the continuous and stable development of the city's level of urbanization, and the quality of the urbanization level can be improved.