A Study of the Time–space Evolution Characteristics of Urban–rural Integration Development in a Mountainous Area Based on Esda-gis: the Case of the Qinling-daba Mountains in China

The multi-index comprehensive evaluation method is used in this paper to estimate the urban–rural integration development level of the Qinling-Daba Mountains in China and build an evaluation index system that is composed by five subsystems and 18 basic indexes. The integration of the indexes is conducted through linear weight sum method and the weight of each index is determined through analytic network process to get the urban–rural integration development coordination indexes of each county in the Qinling-Daba Mountains. Meanwhile, the time–space evolution characteristic of urban–rural integration development in the Qinling-Daba Mountains in the past 10 years is studied through exploratory spatial data analysis and GIS technology. The results indicate that urban–rural integration displays a spatial imbalance, with a situation of polarization, i.e., high–level counties gathering with each other; the spatial aggregation rule of west–high, east–low is coupled with land conditions, industrial activity, and the transportation situation to affect the urban–rural development. Finally, we propose that development priorities should be circular industry, eco-agriculture, and tourism to accelerate urban–rural development and work towards a comprehensive modern transportation system and other infrastructure in the Qinling-Daba Mountains.


Introduction
Urban and rural areas are two entities that interact with and mutually influence each other in regional development.The urban-rural relationship is the most basic economic and social relationship of human beings.At present there is great disparity between urban and rural areas in terms of the level of development in China.As rapid industrialization and urbanization continue, urban-rural transformation and sustainable rural development have become important issues at the forefront of geographical research on the urban-rural relationship [1].
With the rapid development of new-type urbanization processes in China [2], achieving regional sustainable and balanced development has become an important part of overall harmonious urban-rural development [3], and a strategic mission to solve the urban-rural dual structure and reduce differences between urban and rural areas.Urban-rural integration development was put forward as a major strategy for sustainable rural development in the 18th Central Committee of the Communist Party of China, in 2012.In response to these policy propositions and requirements, the geographic study of coordinated urban-rural development and the construction of the modern countryside should be strengthened [4].It will construct a kind of effective mechanism between urban and rural areas,

Natural Resource and Environmental Basis
The Qinling-Daba Mountains are located in the south of Shaanxi province.A major ridge of the Qinling Mountains borders the Guanzhong Plain basin area and Micang Mountain, a major ridge of the Daba Mountains, borders Sichuan province.The Qinling-Daba Mountains stretch in the east-west direction with huge mountains, deep valleys, and numerous river valleys.The area of basin and flatland is less than 10%.The horizontal and vertical zonation of climate in the Qinling-Daba Mountains is noticeable.The Qinling-Daba Mountains have abundant forest specialties, mineral resources, and ecological diversity across their two climate zones, subtropical and warm temperate.In the southern Shaanxi Province, the rainfall is abundant, the topography is greatly changed, and rivers provide abundant water resources.Except for the broad valley basin along the Hanjiang River and a partial small graben basin concentrated with cultivated lands and densely populated, most areas are sparsely populated with good forest cover; a large area of old-growth forest is still reserved in the remote mountains.

Overview of Social and Economic Development
The Qinling-Daba Mountains are in three prefecture-level administrative districts (from west to east): Hanzhong, Ankang, and Shangluo.The national territorial area is 7.01 × 10 4 km 2 , accounting for 34.1% of the whole province, and the total population is 8,390,000, accounting for 22.5% of the whole province, with a population density of 120 persons/km 2 .Compared with the other two areas of Shaanxi Province (Guanzhong and Northern Shaanxi), the overall economy of the Qinling-Daba Mountains is still weak, and economic development is lagging.Since 2000, the total output of southern Shaanxi has been less than 1/5 that of the whole province, with the proportion reducing year after year.GDP per capita has also been at a minimum among the three areas, equal to 1/2 of the average level of the whole province.In the Qinling-Daba Mountains land supply capacity is insufficient, road conditions are poor, and the lack of economic development leaves citizens living in poor conditions and experiencing strict ecological constraints.The lack of high-level talent, few job opportunities, and undeveloped industrial system mean that the region faces competition from Guanzhong and northern Shaanxi.In addition, Southern Shaanxi is also the source of the national Middle Route Project of the South-to-North Water Diversion and Hanjiang-to-Weihe River Diversion Project of Shaanxi Province, which subjects the region to high environmental protection pressures.Mineral resources, the secondary industry that development depends on, are mainly distributed within ecological preservation areas and water conservation areas, so the contradiction between resource development and ecological protection is serious.

Data Sources and Pre-Processing
In order to grasp the urban-rural integration development characteristics of the Qinling-Daba Mountains, it is necessary to make identifications for a certain historical period.In this paper, the essential data of 2003-2013 are adopted to analyze the urban-rural integration development level of the territory year after year.The original data used in this research derive from the Shaanxi Statistical Yearbook, the Socioeconomic Statistical Yearbook of County (city) in China, the Shaanxi Regional Statistical Yearbook, the Hanzhong Statistical Yearbook, the Ankang Statistical Yearbook, and the Shangluo Statistical Yearbook.Any index that is not included in the individual statistical yearbooks is derived from the related government department.
Original data will be given a data standardization treatment first, including data chemotaxis processing and dimensionless method.The data chemotaxis processing is conducted through taking the reciprocal for index of inverse correlation, and the reciprocal would be the mapped interval of [0, 1] by the Min-Max method.

Multi-Index Comprehensive Evaluation Method
The connotations of urban-rural integration development are abundant, not only in the economic field but also in terms of society, culture, and the environment [15][16][17].In this highly complex situation, judgment depends on the technical support of a multi-index comprehensive evaluation method.The main elements of the urban-rural development evaluation index system were suggested by domestic scholars to reflect the contents of the human development index and the new millennium goals of the United Nations.For example, 17 specific indexes were selected and the regional urban-rural integration development evaluation index system was built from three dimensions [18]; three first class indicators and 23 second class indicators were adopted in the urban-rural integration development evaluation index system was designed to conduct an empirical analysis of the urban-rural integration development level of Nanjing [19].The indicator system and methodology for evaluating urban-rural integration development degree was established [20].On the operational level, the multi-index comprehensive evaluation turns a descriptive index with different evaluation objects and dimensions into a dimensionless, relative evaluation value to reach a general evaluation of this object.This approach is meant to solve the relative problems of different dimension indexes and reflect the general trends and proficiency of objects being evaluated, so as to make a fair comparison and solve sequencing problems, etc.
In this paper, the linear weight sum method is used to conduct an integration of the indexes Equation (1); the function of weight in the linear weight sum method is more obvious than in other synthesis methods.Linear compensation can be conducted between indexes; changes of index weight might reflect that there are some big and small functions each index plays in comprehensive assessment.The formula is as follows: where Z is the integrated assessment value gained by the object being evaluated, w is the weight of the index being evaluated, x is the evaluation value of a single indicator, and n is the number of the evaluation index.

Network Analysis Method
The evaluation of regional urban-rural integration development requires a multi-factor and multilevel system assessment.The influences on urban-rural integration development of each index in the indicator system must be determined, namely the weight.There are many factors that influence urban-rural integration development, involving all aspects of society and the economy, with mutual effects and dependence relationships among the influence factors.
The predecessor of the Analytic Network Process (ANP) is the Analytic Hierarchy Process (AHP) proposed by Professor Thomas L. Saaty of the University of Pittsburgh.In the AHP, internal factors in the same hierarchy are relatively independent.There is no mutual influence or dependent relationship between the factors [21,22].Nevertheless, various factors in urban-rural integration have a mutual impacted, making it hard to divide various factors into an internally independent, progressive, hierarchical structure.The interaction between factors mainly presents a network-like structure.The ANP exactly adapts to such a demand and is developed into a systematic decision-making method extended from the AHP.For an internally independent network structure, ANP may be more suitable.
The first step when conducting a network analysis is to construct the ANP structure, confirm a control layer, analyze the relationship between factors in accordance with control standards, and confirm the relationships between elements.
If there are control standards P 1 , P 2 , • • • , P m in a control layer and there are factor groups in the control layer are regarded as standards and factors u jk k = 1, 2, • • • , n j in the factor groups U j (j = 1, 2, • • • , N) as secondary standards.Concentrated factors U i will conduct pairwise comparison in line with the degree of influences to u jk Equation (2).The measurement of importance in pairwise comparison can also apply the rule of nine, and a judging matrix under standard P s is constructed.
Next, an ordering vector w is obtained by a group eigenvalue method.Moreover, the abovementioned characteristic vectors are subjected to a consistency check.Similarly, the ordering vector relative to other factors can be obtained.Moreover, the local weight vector matrix can be obtained Equation (3): The column vector in W ij is the importance ordering vector of factors If factors in U j will not be impacted by U i , W ij = 0, the factor set in the network layer will require a pairwise combination to repeat the abovementioned steps.Finally, a hypermatrix W under the control standard P s can be obtained Equation (4).Thus s = 1, 2, • • • , m, the amount of such a hypermatrix will be m: In the hypermatrix, every element is a sub-matrix.The sum of every column is 1, but W is not a normalized matrix.In order to make subsequent calculations more convenient, we should regard P s as a standard when comparing the importance of factors.The judging matrix of factor groups Equation ( 5) under the standard P s is: If factor groups are independent of U j , the importance vector is 0 and the weighting matrix is Equation ( 6 After weighting elements of hypermatrix W, the weighting hypermatrix W is shown as follows Equation (7).The sum of columns in W is 1, and In order to reflect the dependence relationships between factors in the network, it is necessary to perform a stabilizing treatment on the weighting hypermatrix to calculate the relative limit ordering vector of every weight hypermatrix Equation ( 8): If the limit is convergent and unique, the value in the corresponding line of the original matrix is the stable weight of the factors, namely, column W k is the relative limit ordering vector of factors in the network layer on the factor j under the standard P s .ANP is adopted to determine the weight of each index, and factors in the system are divided into two parts, a control layer and a network layer.The control layer also includes two parts, the target and the decision criterion.The decision criterion is independent, only controlled by target factors; there can be no decision criterion in the control layer, but a simple target is of the essence.The weight of each decision criterion can be determined through a traditional analytic hierarchy process.The network layer is composed of all factors governed by the control layer, and an interactional network structure is formed among the factors [23,24].The calculation steps of ANP are complex and the core work is solving hypermatrix [25,26].

Exploratory Spatial Data Analysis
Exploratory Spatial Data Analysis (ESDA) is a spatial mechanism of interaction that, through calculating the spatial autocorrelation of research objects, would be used for the description and visual analysis of a geographic phenomenon spatial pattern to discover the spatial agglomeration and spatial anomalies as well as reveal the object of study [27].ESDA includes two types of analytical perspectives, global spatial autocorrelation analysis and local spatial autocorrelation analysis [28][29][30].
Global spatial autocorrelation analysis is used for determining spatial features of geographic phenomena in the whole study area, whether or not spatial agglomeration exists.Moran's I index is a measurement index of global spatial autocorrelation [31].The value of Moran's I is within the scope of [−1, 1].Moran's I > 0 means positive spatial correlation; Moran's I < 0 means negative spatial correlation; Moran's I = 0 means uncorrelated.
Local indicators of spatial association (LISA) are used for measuring the spatial difference and significance between a certain place and the surrounding areas.Global spatial autocorrelation analysis only reveals the spatial dependency of geographical phenomena as a whole and cannot provide local spatial difference comparison, so it is necessary to conduct local spatial autocorrelation analysis [32].Local Moran's I index is generally used as a local spatial autocorrelation measurement index to study local spatial distribution patterns by combining the Moran Scatter Plot and LISA Aggregation Plot; the aim is to resolve the global Moran's I to each object of study [33].
Under the given significance level, spatial agglomeration will occur for similar values of geographic attributes if you want to show positive local spatial autocorrelation; spatial agglomeration will occur for dissimilar values of geographic attributes if you want to show negative local spatial autocorrelation.The "hot spot" and "cold spot" of local spatial agglomeration can be identified and the singular spatial value can be revealed by a LISA aggregation map, formed through a combination of Scatter Plot with LISA significance level [12].

Evaluation Object
The 25 county territories of the Qinling-Daba Mountains act as elementary units (Figure 1) when evaluating the urban-rural integration development level in view of the country's economy.

Evaluation Object
The 25 county territories of the Qinling-Daba Mountains act as elementary units (Figure 1) when evaluating the urban-rural integration development level in view of the country's economy.

Establishment of Index System
Urban-rural integration development aims at changing the dual model of "urban industry, rural and agricultural" by combining urban and rural developments with unified coordination, following recent economic and social development thought on the integration of workers and peasants, and solving existing urban and rural questions by developing foresight.
There is no uniform indicator weighting system for assessing the integrated development of urban and rural areas in China.This has been one of the hotspots studied by Chinese scholars in recent years.Existing studies on urban-rural integration mainly apply an aggregative indicator system to select an indicator.In view of the development connotations of urban-rural integration, generally speaking, indicator selection should reflect the differences in population, space, economy, society, and the environment between urban and rural areas [5].However, for undeveloped mountainous areas, besides a reduction of the urban-rural gap, it is extremely important to improve the overall development level [34][35][36][37][38].In this paper, urban-rural development shows a dual structure, with different population, space, economy, society, and environment.Different subsystems of human society can be formed based on the five aspects mentioned above, namely population, space, economy, society, and the environment.In order to achieve the goal of urban-rural integration, benign interaction between urban and rural areas must be realized for these five subsystems that coordinate with each other.Construction of a multi-index evaluation index system should emphasize a systematic, scientific, feasible, and comparable approach.When measuring the urban-rural integration development level of a population subsystem, the indexes that need to be selected are as follows: three basic indexes of urbanization rate (x 1 ), occupied population proportion of non-agricultural industries (x 2 ), and non-agricultural industries proportion in country working population (x 3 ); two basic indexes of urban density (x 4 ) and penetration rate of telephone (x 5 ) for space subsystem; seven basic indexes of Per Capita GDP (x 6 ), economic density (x 7 ), proportion of non-agricultural industries output value (x 8 ), agricultural mechanization level (x 9 ), fixed-asset investment per capita (x 10 ), total retailing of social consuming goods per capita (x 11 ), and urban and rural residents per capita income ratio (x 12 ) for economy subsystem; four basic indexes of teacher/pupil ratio in common elementary and middle schools (x 13 ), number of health agency beds per 10,000 people (x 14 ), number of health workers per 10,000 people (x 15 ), and number of social welfare adopting unit beds per 10,000 people (x 16 ) for society subsystem; and two basic indexes of green coverage ratio in built-up areas (x 17 ) and per capita public green areas in built-up areas (x 18 ) for the environmental subsystem.Therefore, when judging the urban-rural integration development level, we can start with the five subsystems that represent stages in the process of urban-rural integration development.

Influence Factor Weight
The urban-rural integration development evaluation index system of the Qinling-Daba Mountains has been built as mentioned earlier.There are 18 basic indexes in the system, which can be regarded as factors of the network layer.The network layer should include the target layer and criterion layer; in this paper, there is only one target layer in the control layer and no evaluation criteria layer.Factors of the network layer can be directly compared through a comprehensive assessment of urban-rural integration development level as the judgment standard (Figure 2).The ANP evaluation model is established through Super Decisions 2.0.8 software to find a matrix solution and obtain the weight of each influencing factor (Table 1).
The control layer points to the factor group box in the network layer, indicating that the control layer is impacted by the Factor Group in Figure 2. In the box Factor Group in the network layer, it refers to the factor set of describing the same property.The arrow points to Factor Group B from Factor Group A, indicating that Factor Group A is impacted by Factor Group B. If both ends of the connecting line between Factor Group A and Factor Group B have an arrow, this indicates an interaction between Factor Group A and Factor Group B. A circular arrow indicates the interaction of factors in the Factor Group.Corresponding values can be obtained by bringing each essential data point and influence factor weight into the formula of the linear weight sum method (Table 2).

Historical Evolution of the Urban-Rural Integration Development Level of the Qingling-Daba Mountains
The urban-rural integration development coordination index of the Qingling-Daba Mountains is used in this paper as a kind of spatial attribute, and ArcGIS 10.0 is used to analyze and simulate its spatial trends to get a three-dimensional figure (Figure 3).The arrow direction of the X-axis in the figure is east, the arrow direction of the Y-axis is north, and the Z-axis is the coordination index of urban-rural integration development.Points on the XZ plane are the projection of the urban-rural integration development coordination index of each county in longitude direction, while curves on the XZ plane are the spatial trend of the urban-rural integration development level in the longitude direction.Points on the YZ plane are the projection of the urban-rural integration development coordination index of each county in the latitude direction, while curves on the YZ plane are the spatial trend of the urban-rural integration development level in the latitude direction.
The result of spatial trends analysis indicates that the urban-rural integration development level of the Qinling-Daba Mountains shows a high-west, low-east trend in the longitude direction and a trend of high in the middle area, relatively low in the north, and lowest in the south in the latitude direction.Specifically speaking, in the longitude direction, counties in the west are higher than the counties in the middle; the east part on the whole is not as good as the western part, although the eastern section of the trend line starts to rise slightly from 2010.The spatial rule of the urban-rural integration also has a close connection with the industrial basis and transportation conditions of the evaluated unit.As a significant zone of the Third-line Construction, Hanzhong has a preferable industrial location; industrial activity gets weaker the further east you go in the southern part of Shaanxi Province.The western part is higher than the eastern part and the urban-rural integration development level is identical to the spatial rule for industrial development.Hanzhong was the earliest zone to have a railway and highway.Summing up the above, the high-west, low-east spatial clustering rule of urban-rural integration development in the Qinling-Daba Mountains has a strong correlation with the spatial characteristics of the land, industrial activity, and the transportation network.

Exploratory Spatial Data Analysis of the Urban-rural Integration Development Level in the Qinling-Daba Mountains
GeoDa 1.4.6 software is used to calculate the overall and partial spatial relativity evaluation index of urban-rural integration development coordination index for each county in the Qinling-Daba Mountains year by year [39,40].
The Queen neighborhood space matrix technique is used to determine spatial weight Equation ( 9) [41], that is: The Hanjiang River traverses the Qinling-Daba Mountains from the west to the east, then flows into Hubei Province.The main stream of Hanjiang River flows through six counties (Ningqiang, Mian, Nanzheng, Chenggu, Yang, and Xixiang) in the Hanzhong and four counties (Shiquan, Hanyin, Xunyang, and Baihe) in Ankang.A planting region of the leading grain crop of rice of the Qinling-Daba Mountains is also located along the main stream of the Hanjiang River.It is also an intensive belt of industries and cities and towns, and the main channel for communication, transportation, exchange of information, and personnel mobility.The intensive degree, average scale, and population density degree of the cities and towns along the main stream of Hanjiang River are all higher than those of the southern and northern mountains, providing more space for the social and economic development of the region.Urban and rural development cannot happen without an exchange of materials and information, a main railway, and a motorway.The national and provincial main roads located along the main stream of Hanjiang River provide obvious advantages in terms of communication and transportation.Under the joint action of the above advantages, the urban-rural integration development level of the Qinling-Daba Mountains reflects a spatial clustering characteristic: the zone along the main stream of Hanjiang River is higher than the southern and northern mountains.
The spatial rule of the urban-rural integration also has a close connection with the industrial basis and transportation conditions of the evaluated unit.As a significant zone of the Third-line Construction, Hanzhong has a preferable industrial location; industrial activity gets weaker the further east you go in the southern part of Shaanxi Province.The western part is higher than the eastern part and the urban-rural integration development level is identical to the spatial rule for industrial development.Hanzhong was the earliest zone to have a railway and highway.Summing up the above, the high-west, low-east spatial clustering rule of urban-rural integration development in the Qinling-Daba Mountains has a strong correlation with the spatial characteristics of the land, industrial activity, and the transportation network.

Exploratory Spatial Data Analysis of the Urban-rural Integration Development Level in the Qinling-Daba Mountains
GeoDa 1.4.6 software is used to calculate the overall and partial spatial relativity evaluation index of urban-rural integration development coordination index for each county in the Qinling-Daba Mountains year by year [39,40].
The Queen neighborhood space matrix technique is used to determine spatial weight Equation ( 9) [41], that is: Area j and i have common boundary or common verte 0 i = j or area J and I without common boundary or common vertex .
In this formula, w ij means the spatial weight of area j to area i.

Global Spatial Autocorrelation Analysis
At a significance level of α = 0.05, Z-value ≥ 1.96 can be tested.The annual p-value and Z-value of global Moran's I in counties in the Qinling-Daba Mountains can pass the test (Table 3).The global Moran's I value of the urban-rural integration development coordination index in the Qinling-Daba Mountains is all positive, indicating a positive correlation in space, namely high value zones are close to high value zones and the low value zones are close to low value zones.It can be seen that spatial clustering existed in terms of the urban-rural integration development of counties in the Qinling-Daba Mountains, but the clustering intensity was not strong (0.105 < Moran's I < 0.297).Among them, the year with the highest clustering intensity was 2011 and the year with the lowest clustering intensity was 2007.According to the changing Moran's I values over the years (Figure 4), the urban-rural integration development level in the Qinling-Daba Mountains can be divided into two phases.The first phase is from 2003 to 2007, when the spatial clustering extent of the urban-rural integration development decreased year after year.The polarization in space decreased, and there was a similar trend for regional development.The industrial development in the Qinling-Daba Mountains established a certain baseline during the Third-line Construction period in the late 1960s and 1970s, but the industry activity has reduced gradually due to inconvenient transportation, miscommunication of information, and slow update of techniques, etc.Thus, the spatial clustering of regional urban-rural integration development reduced, and the entire region was in a low-level, balanced development state.
According to the changing Moran's I values over the years (Figure 4), the urban-rural integration development level in the Qinling-Daba Mountains can be divided into two phases.The first phase is from 2003 to 2007, when the spatial clustering extent of the urban-rural integration development decreased year after year.The polarization in space decreased, and there was a similar trend for regional development.The industrial development in the Qinling-Daba Mountains established a certain baseline during the Third-line Construction period in the late 1960s and 1970s, but the industry activity has reduced gradually due to inconvenient transportation, miscommunication of information, and slow update of techniques, etc.Thus, the spatial clustering of regional urban-rural integration development reduced, and the entire region was in a low-level, balanced development state.The second phase is from 2008 to 2013, when the spatial characteristic of urban-rural integration development is aggregation and the spatial polarization starts to recover.The Shaanxi provincial party committee and provincial government put forward a regional development strategy of "breakthrough development in Southern Shaanxi", and confirmed that green industry should be the main direction in the Qinling-Daba Mountains.In the meantime, environmentally sustainable travel The second phase is from 2008 to 2013, when the spatial characteristic of urban-rural integration development is aggregation and the spatial polarization starts to recover.The Shaanxi provincial party committee and provincial government put forward a regional development strategy of "breakthrough development in Southern Shaanxi", and confirmed that green industry should be the main direction in the Qinling-Daba Mountains.In the meantime, environmentally sustainable travel has been another economic growth point for some counties.Moreover, the transportation situation has been tremendously improved.Thus, the overall economic and social development in the Qinling-Daba Mountains during this period has been rapid; some counties have green agriculture and ecological tourism resources in place, and have broken through the external traffic bottleneck.There has been a spatial clustering of transportation facilities and tourism resources.

Local Spatial Autocorrelation Analysis
The spatial relationship between a certain territory unit and its adjacent units can be analyzed through the drawing of a Moran scatter diagram.The Moran scatter diagram for each year in the Qinling-Daba Mountains (Figure 5) indicates that the urban-rural integration development level of 25 counties in the Qinling-Daba Mountains was clustered around two terminals.The Moran scatter diagram can be divided into four quadrants, corresponding to four types of local spatial forms between the county and its neighbors.The first quadrant (H-H) shows that regional units of high observed value are enclosed by a high-value area.The third quadrant (L-L) shows that regional units of low observed value are enclosed by a low-value area.The counties in the H-H and L-L quadrants are different every year (Table 4).It can be observed from Figures 3, 5, and 6, and Table 4 that the urban-rural integration development level of each county in the Qinling-Daba Mountains is clustered around two terminals-that is, the high-level counties cluster with the high-level counties and the low-level counties cluster with the low-level counties.Counties in the western part, located along the main stream of the Hanjiang River and the southern Qinling Mountains, are more likely to be in the first quadrant, while the opposite is true of the third quadrant.This shows that the urban-rural integration development of the Qinling-Daba Mountains lacks spatial balance.The development in this region is not balanced due to differences in resources endowment, land conditions, and transportation facilities in the 25 counties within the region; these are conditions that the urban-rural economies and social development rely on.The Qinling-Daba Mountains lack counties with strong leading effects, the spatial polarization is weak, and spatial clustering is still the overall trend for the urban-rural integration development during this period.

Conclusions and Discussion
Studies show that the urban-rural integration development level of the Qinling-Daba Mountains displays a spatial imbalance, with a "Higher in the West and Lower in the East" spatial variation and a situation of aggregating to the ends.Namely, high-level counties and low-level counties gather with each other, and the spatial characteristics couple with the land conditions, industrial activity,

Conclusions and Discussion
Studies show that the urban-rural integration development level of the Qinling-Daba Mountains displays a spatial imbalance, with a "Higher in the West and Lower in the East" spatial variation and a situation of aggregating to the ends.Namely, high-level counties and low-level counties gather with each other, and the spatial characteristics couple with the land conditions, industrial activity, and transportation situation to affect the urban-rural development.The socioeconomic development in the study area is restricted by the land conditions, ecological protection requirements, etc. [42].A new model for improving regional development and promoting a balance of urban and rural structures is necessary to solve a series of problems, such as urban and rural structural imbalance, rural development lag, etc. [43,44].
The spatial strategy for the integrated development of urban and rural areas in the Qinling-Daba Mountains of Southern Shaanxi should be concentrated on the Hanzhong Basin, the Yuehe River, and Shangdan Valley.These areas also have the highest integrated development of urban and rural areas.It is necessary to continue to develop the driving role of urban areas on rural development, form a core regional development area with stronger absorption and radiation, and encourage overall development.
The northern area of the Qinling-Daba Mountains in Southern Shaanxi is the ecological protection zone of the Qinling Mountains and has a relatively lower integrated development of urban and rural areas.The area is a water conservation area of Chinese South-to-North Water Diversion Project.Urban-rural development is strictly limited by the ecological protection policies.On the one hand, it is necessary to advocate the concept of a circular economy and an ecological economic development mode, to drive the optimization and promotion of traditional agriculture, develop a modern service industry with a focus on tourism, increase the income of branch agriculture industry, and help populations in mountainous areas overcome poverty.On the other hand, we must encourage populations in mountainous areas to move to flat areas with a good economic outlook along the Hanjiang River, increase urban capacity, and further reduce the gap between urban and rural development.
The improvement of regional traffic facilities is an important driver to promote the overall improvement of urban-rural development in the Qinling-Daba Mountains of Southern Shaanxi and reduce the gap between urban and rural development.Therefore, it is necessary to further speed up traffic system construction, especially the construction of modern traffic facilities such as airports, high-speed railways, and highways.

Figure 1 .
Figure 1.Division diagram of evaluation unit.

Figure 1 .
Figure 1.Division diagram of evaluation unit.

Figure 2 .
Figure 2. The structure of ANP comprehensive evaluation level of urban-rural integration development in the Qinling-Daba Mountains.

Figure 2 .
Figure 2. The structure of ANP comprehensive evaluation level of urban-rural integration development in the Qinling-Daba Mountains.
region.Urban and rural development cannot happen without an exchange of materials and information, a main railway, and a motorway.The national and provincial main roads located along the main stream of Hanjiang River provide obvious advantages in terms of communication and transportation.Under the joint action of the above advantages, the urban-rural integration development level of the Qinling-Daba Mountains reflects a spatial clustering characteristic: the zone along the main stream of Hanjiang River is higher than the southern and northern mountains.

Figure 3 .
Figure 3. Spatial trend analysis of urban-rural integration development level in the Qinling-Daba Mountains.

Figure 3 .
Figure 3. Spatial trend analysis of urban-rural integration development level in the Qinling-Daba Mountains.

Figure 4 .
Figure 4. Global Moran's I changing values of urban-rural integration development coordination in the Qinling-Daba Mountains.

Figure 4 .
Figure 4. Global Moran's I changing values of urban-rural integration development coordination in the Qinling-Daba Mountains.
spatial clustering is still the overall trend for the urban-rural integration development during this period.

Figure 5 .
Figure 5. Moran scatter plot of urban-rural integration development coordination index for each year in the Qinling-Daba Mountains.

Figure 6 .
Figure 6.LISA maps of urban-rural integration development in the Qinling-Daba Mountains.

Table 1 .
Global weighted table of comprehensive assessment factors in urban-rural integration development, Qinling-Daba Mountains.

Table 1 .
Global weighted table of comprehensive assessment factors in urban-rural integration development, Qinling-Daba Mountains.

Table 2 .
Urban-rural integration development coordination index table of each county in theQinling-Daba Mountains.

Table 3 .
Global Moran's I value table for the urban-rural integration development coordination index of the Qinling-Daba Mountains.