Construction of Ecological and Recreation Patterns in Rural Landscape Space: A Case Study of the Dujiangyan Irrigation District in Chengdu, China

: The rapid expansion of urbanization has promoted the prosperity of the economy and society but has also caused rural ecological problems. This study takes the Dujiangyan Irrigation District as an example, to construct ecological and recreation patterns, and it breaks through the pattern construction mode of using a single ecological factor that has been used in previous studies. We analyzed the impact of landscape connectivity and area thresholds at different scales on the selection of source areas, and integrated urban construction and human disturbance factors into resistance surface construction. Finally, a comprehensive landscape network combining “ecology and recreation” was determined through the minimum cumulative resistance model. Multiscale landscape connectivity analysis and area threshold setting greatly promoted the results of source identiﬁcation. After optimization, we identiﬁed four ecological corridors and twenty-seven recreational corridors, and the regional landscape security pattern was signiﬁcantly improved compared with the previous upper planning content. Therefore, this study provides a reference for regional long-term planning and has reference signiﬁcance for the spatial protection and utilization of rapidly urbanizing areas.


Introduction
Since the 1980s, China's rural land has expanded unprecedentedly with the development of society and the economy. This expansion trend has replaced the traditional farming and living space in rural areas with more modern and highly intensive production and lifestyles and has led to a series of ecological environment problems, such as the intensification of landscape fragmentation and the weakening of ecosystem services [1,2]. In this context, China has proposed the rural revitalization strategy, taking "ecological revitalization" as an important part of rural revitalization, bringing opportunities for the development of rural areas [3,4]. Therefore, from the perspective of landscape patterns, it is of great significance to explore the risks faced by landscape patterns, balance the relationship between rural ecological environmental protection and development and solve the contradiction in development from the perspective of environmental friendliness to optimize resource allocation, and promote the coordinated development of space.
At present, landscape pattern research is mature in ecology, and the application scenarios are also more extensive. By analyzing the characteristics and driving factors of spatio-temporal changes in landscape patterns, some researchers have discussed and analyzed the contribution degree of various elements in landscape patterns to space and the correlation, among various factors [5,6]. By analyzing the coupling degree between landscape patterns and ecosystem service function and value, some scholars have proposed key elements and specific measures to improve ecosystem service function [7]. However, most existing landscape pattern research relies on landscape index analysis, and landscape Land 2022, 11, 383 3 of 21 study relies on the landscape evaluation results, selects the landscape index to analyze the landscape pattern connectivity and area threshold, and determines the scope of ecological and recreational source areas, respectively, with objective data, hoping to improve the accuracy and scientific nature of source area identification [31][32][33]. Secondly, based on the data of land-use types, roads, water systems, landform and geomorphic conditions, the ecological and recreational resistance surfaces are constructed, respectively, and the resistance surfaces integrating multiple elements can make the expression of ecological process closer to the real space [34]. Finally, the ecological and recreational corridors are identified based on the MCR model, and the comprehensive expression of spatial development, protection and utilization is realized by the compound landscape corridor [35] (Figure 1).

Research Design
Relying on the basic model of "source identification; resistance surface construction; corridor extraction", this study expects to build a comprehensive pattern integrating ecology and recreation, so as to realize the coordinated development of ecological protection and recreation development in this area. Firstly, for the identification of source areas, this study relies on the landscape evaluation results, selects the landscape index to analyze the landscape pattern connectivity and area threshold, and determines the scope of ecological and recreational source areas, respectively, with objective data, hoping to improve the accuracy and scientific nature of source area identification [31][32][33]. Secondly, based on the data of land-use types, roads, water systems, landform and geomorphic conditions, the ecological and recreational resistance surfaces are constructed, respectively, and the resistance surfaces integrating multiple elements can make the expression of ecological process closer to the real space [34]. Finally, the ecological and recreational corridors are identified based on the MCR model, and the comprehensive expression of spatial development, protection and utilization is realized by the compound landscape corridor [35] (Figure 1).

Study Area
Dujiangyan city is located west of Chengdu, the capital of China's Sichuan Province. It received its name from the world-famous Dujiangyan Irrigation System. Dujiangyan Irrigation System is one of the few world heritage sites and has three main categories: World Natural Heritage, World Cultural Heritage and World Irrigation Engineering Heritage. Its water conservancy culture enjoys a strong reputation worldwide. The irrigation area affected by the Dujiangyan Irrigation System exceeds 7100 km 2 . Approxi-

Study Area
Dujiangyan city is located west of Chengdu, the capital of China's Sichuan Province. It received its name from the world-famous Dujiangyan Irrigation System. Dujiangyan Irrigation System is one of the few world heritage sites and has three main categories: World Natural Heritage, World Cultural Heritage and World Irrigation Engineering Heritage. Its water conservancy culture enjoys a strong reputation worldwide. The irrigation area affected by the Dujiangyan Irrigation System exceeds 7100 km 2 . Approximately 254 km 2 of land is within the eastern plain area of Dujiangyan city, which in turn is located within the core area of the Dujiangyan Irrigation System (30 • 44 -31 • 02 N, 103 • 0 -103 • 47 E). The area has the geomorphological characteristics of a typical plain irrigation area and was selected as the scope of the study in this paper. The region covers three towns (approximately 80 villages) with a total population of 242,100. This area is an important implementation area of the rural revitalization strategy of Chengdu. By choosing this area as the research focus of this paper, we have a clear understanding of how landscape space should be protected and utilized under the background of the rural revitalization strategy (Figure 2). 103°0′-103°47′ E). The area has the geomorphological characteristics of a typical plain irrigation area and was selected as the scope of the study in this paper. The region covers three towns (approximately 80 villages) with a total population of 242,100. This area is an important implementation area of the rural revitalization strategy of Chengdu. By choosing this area as the research focus of this paper, we have a clear understanding of how landscape space should be protected and utilized under the background of the rural revitalization strategy (Figure 2).

ਰ 峖 Ѐ Types of Data
The data types of this paper mainly include: (1) digital elevation model (DEM) data from the geospatial data cloud (http://www.gscloud.cn/ (accessed on 6 July 2021)), which are mainly used for the calculation of topographic factors; (2) land-use remote sensing images from the geospatial Data Cloud Platform of the Computer Network Information Center, Chinese Academy of Sciences (http://www.gscloud.cn/ (accessed on 8 July 2021)). The Landsat 8 OLI-TIRS image for 28 July 2020 (Scene No. 129-039), during a period of good plant growth and low cloud cover, was selected as the interpretation data source, which was mainly used for processing and analysis to obtain land-use types; (3) Open Street Map (OSM) data from the Open Street Map website (https://www.openstreetmap.org/ (accessed on 25 November 2020)), which were used for spatial network analysis and human activity disturbance analysis; (4) the Normalized Difference Vegetation Index (NDVI), which originated from the Chinese Academy of Sciences Resources and Environment Data Service Center (http://www.resdc.cn/Default.aspx (accessed on 8 July 2021)) and was mainly used for measuring regional vegetation coverage; (5) point of interest (POI) data that were derived from AmAP (https://lbs.amap.com/ (accessed on 27 May 2021)) and mainly used to analyze the distribution of various service facilities in the region; (6) the water system map of the Dujiangyan Irrigation District (including the irrigation system) and survey data of landscape resources, which were obtained from the Dujiangyan Irrigation District Management Committee and mainly used for water system analysis; and (7) the survey data

Types of Data
The data types of this paper mainly include: (1) digital elevation model (DEM) data from the geospatial data cloud (http://www.gscloud.cn/ (accessed on 6 July 2021)), which are mainly used for the calculation of topographic factors; (2) land-use remote sensing images from the geospatial Data Cloud Platform of the Computer Network Information Center, Chinese Academy of Sciences (http://www.gscloud.cn/ (accessed on 8 July 2021)). The Landsat 8 OLI-TIRS image for 28 July 2020 (Scene No. 129-039), during a period of good plant growth and low cloud cover, was selected as the interpretation data source, which was mainly used for processing and analysis to obtain land-use types; (3) Open Street Map (OSM) data from the Open Street Map website (https://www.openstreetmap.org/ (accessed on 25 November 2020)), which were used for spatial network analysis and human activity disturbance analysis; (4) the Normalized Difference Vegetation Index (NDVI), which originated from the Chinese Academy of Sciences Resources and Environment Data Service Center (http://www.resdc.cn/Default.aspx (accessed on 8 July 2021)) and was mainly used for measuring regional vegetation coverage; (5) point of interest (POI) data that were derived from AmAP (https://lbs.amap.com/ (accessed on 27 May 2021)) and mainly used to analyze the distribution of various service facilities in the region; (6) the water system map of the Dujiangyan Irrigation District (including the irrigation system) and survey data of landscape resources, which were obtained from the Dujiangyan Irrigation District Management Committee and mainly used for water system analysis; and (7) the survey data of landscape space feeling, which were mainly used for the evaluation of recreational space. The usage method and process of specific data sources are specifically introduced in the research methods.

Construction of Ecological and Recreational Space
Spatial construction is the cornerstone used to realize pattern optimization. In view of the current situation of natural and cultural landscapes in the research region, considering the accessibility of data and referring to relevant literature, the study started from the two aspects of the "natural environment and degree of human disturbance". The evaluation factors of ecological space included elevation, slope, slope aspect, water area, land use, vegetation coverage, disturbance of development and construction to the ecological environment (road distribution, residential area and population density distribution). For recreational space, this paper selected five recreational factors: distribution of tourism resources and facilities, distribution of cultural relics and historical sites and potential tourism resources, human demand for recreational activities and leisure space. Finally, the evaluation framework of ecological and recreational space was constructed based on the elements and references of relevant literature [36][37][38][39][40][41][42].
Then, for the weight determination of each index, the weighting method of subjective and objective combinations of the weighting method (analytic hierarchy process (AHP) and entropy weight method) was performed using ArcGIS software to achieve layer stacking to complete the construction of ecological and recreational spaces. Among them, the weight result calculated by the AHP method shows that the heterogeneity ratio of the criterion layer (B) combination CR = 0.0518 and that the consistency ratio of the subcriterion layer (C) combination CR = 0.0538, both of which are less than the threshold value of 0.1, indicating that the judgment matrix passes the consistency test and the weight calculated is meaningful. Finally, according to the result of superposition, the natural discontinuities' classification method (Jenks) on spatial analysis level 5 drawings was used for classification; the categories included the following values: 1 (not sensitive area); 3 (low sensitivity); 5 (moderately sensitive area); 7 (highly sensitive area); and 9 (extremely sensitive area), with lower grades indicating less sensitive landscape space and more benefits for the development of the space and construction (see details from Table A1 in Appendix A) [41,42].

Analysis of Multiscale Landscape Connectivity
Landscape connectivity refers to the degree to which landscape promotes or hinders material flow, energy flow and information flow. Maintaining good landscape connectivity has a dominant influence on regional ecological process development [36,41]. However, landscape connectivity has a strong scale dependence. When the scale changes, smaller patches in landscape space are constantly fused by adjacent or connected patches to form new spatial patches, and landscape connectivity changes [43][44][45]. Therefore, through repeated attempts at landscape spatial suitability at different scales, seven grid files with different sampling scales of 30, 60, 120, 240, 480, 960 and 1800 m were set. Based on previous research experience, six landscape indices, CONHESION, AI, CONTAG, SHDI, SHEI and DIVISION, were selected. The changes in landscape connectivity under different scale conditions were calculated in Fragstats [43,[46][47][48][49]. Notably, CONHESION, AI and CONTAG indicate that dominant patches in landscape space are well connected, while larger SHDI and DIVISION values indicate the landscape is more fragmented and less connected.

Principal Component Analysis
Principal component analysis (PCA) is a method used to simplify and assign weights to high-dimensional variables. PCA can assign the influence degree of relevant variables to corresponding principal component factors. This study used this method to objectively evaluate the overall connectivity of landscape space at different scales and calculated the comprehensive index within the research scope, according to the index of each constraint factor. The composite index is defined as a weighted sum of the principal component, and the right to reuse a variance contribution rate of each principal component is related; finally, this method obtains the ecological source to identify the best space dimension [36,49]. The formula is as follows: where ESI is the spatial index of unit i, P ij is the Jth index of unit i, w j is the weight of each index. However, because the sources and dimensions of different indicators are not unified, the indicators are not comparable; thus, they need to be standardized. Then, the landscape index can be substituted into the weight formula of the spatial PCA results to complete the calculation of the overall connectivity composite index. The min-max normalization method was used [50]. The formula is as follows: where X i is the measured value of index i, X max is the maximum value of index i, X min is the minimum value of index i, and K i is the standardized value of the indicator.

Construction of the Resistance Surface
The resistance surface refers to the carrier of the obstructed degree of species migration in the corridor, and its construction principle is constructed by calculating the cost of material flow, energy flow and information flow spreading among different landscape patches under the action of various resistance elements [29]. Relevant studies usually adopt land-use type data and set different resistance coefficients for different land-use types, such as forestland, water systems and other green space resistance coefficients, which take the minimum value, and construction land and other land-use types, which take the maximum value [34,35]. Therefore, the method of using land-use type as a dominant factor in the resistance surface, and using road classification, water and terrain factors as auxiliary correction, was completed by referencing the relevant literature to determine the assignment of elements. The resistance value was incorporated into the data to form the resistance surface (Table 1). After the construction of the resistance surface was completed, the Jenks method was adopted to divide the cumulative cost of the landscape into the three resistance grades of low, medium and high according to the overall resistance [48,49].

Identification of Corridors
A corridor is a channel connecting different sources and an important guarantee to meet the flow of regional matter, energy and information. The identification and optimization of corridors is conducive to enhancing the stability and integrity of regional functions [50,51]. The MCR model is usually used to analyze the identification of corridors. Its principle is to determine the minimum cost distance between different "sources" and connecting targets through resistance, which has wide adaptability. Therefore, this paper uses the MCR model to calculate ecological and recreational corridors [12][13][14]52]. The formula is as follows: where MCR is the minimum cumulative resistance value, D ij is the spatial distance of species from source j to landscape unit i, and R i is the resistance coefficient of landscape unit i to species movement. ƒ represents a positive correlation between minimum cumulative resistance and ecological processes.

Spatial Analysis of a Single Factor
The spatial distribution of single elements in the Dujiangyan Irrigation District revealed the influence and effect of different influence factors on the overall ecology and recreational space of the region ( Figure 3). Among the 13 data indicators, except C8 and C9, the data were all positively distributed (that is, the higher the value was, the higher the spatial sensitivity was). On the overall spatial distribution, the single factors of C1, C2 and C4 were natural factors, such as distribution in the space, and were relatively fragmented. C5, C7, C8 and other factors were greatly influenced by human activities, as evident in the town domain space and the surrounding space. The difference in the C5 and C7 of the two drawings, the central areas and the northern areas showed obvious low information. Additionally, C8 showed obvious aggregation in the high-value area around the town affected by human activities.
According to the eight single factor analysis drawings of ecological protection importance, it was found that the regional water network was densely covered, there were good water landscape resources and the vegetation coverage was high, showing a good aggregation of ecological resources. However, the influence range of disturbance factors has exceeded the limits of the towns, which is not conducive to the protection of ecological space.

Analysis of Ecological and Recreational Spatial Patterns
The spatial pattern of the evaluation results of the study area in accordance with the sensitivity was divided into five grades, 1, 3, 5, 7, and 9 (the higher the value was, the higher the sensitivity), to obtain the spatial pattern of the ecology and tourism of the Dujiangyan Irrigation District. The proportion of ecological space grades from low to high was 0.001%, 20.59%, 74.21%, 5.18%, and 0.02%, respectively ( Figure 4a). The proportion of recreational space grade from low to high was 2.21%, 22.59%, 50.61%, 24.49% and 0.10%, respectively ( Figure 4b).
Regarding the sensitive areas on the spatial distribution of ecological space, highly sensitive areas were mainly concentrated in the northern area, the eastern mountain area and the southern area along the river, and the ecological safety space was relatively low. The areas with low sensitivity were mainly located in the middle, north and southwest of the region and formed a more obvious spatial boundary with the town. The grades of recreational space were scattered in space, showing an obvious trend of spatial aggregation.

Analysis of Ecological and Recreational Spatial Patterns
The spatial pattern of the evaluation results of the study area in accordance with the sensitivity was divided into five grades, 1, 3, 5, 7, and 9 (the higher the value was, the Regarding the sensitive areas on the spatial distribution of ecological space, highly sensitive areas were mainly concentrated in the northern area, the eastern mountain area and the southern area along the river, and the ecological safety space was relatively low. The areas with low sensitivity were mainly located in the middle, north and southwest of the region and formed a more obvious spatial boundary with the town. The grades of recreational space were scattered in space, showing an obvious trend of spatial aggregation.

Changes in Landscape Connectivity and Scale
By analyzing the different scales of landscape ecology and the tourism connectivity index, the CONHESION, AI, CONTAG, SHDI, SHEI, and DIVISION landscape connectivity measure indices were selected, and different scales were obtained by computing the scatterplot [31][32][33]38]. The connectivity of ecological and recreational spaces showed obvious differences in scale. The AI and CONTAG indices of ecological space both first decreased and then increased, and the CONTAG index decreased sharply at the 120 m scale and then showed a fluctuating trend. The SHDI, DIVISION and SHEI all first increased and then decreased, with high values at the 960 m scale (Figure 5a). For recreational space, except for SHEI, the other indices showed a downward trend, and SHEI showed almost no change at 30-240 m but increased gradually after the 240 m and 1800 m peaks (Figure 5b).

Changes in Landscape Connectivity and Scale
By analyzing the different scales of landscape ecology and the tourism connectivity index, the CONHESION, AI, CONTAG, SHDI, SHEI, and DIVISION landscape connectivity measure indices were selected, and different scales were obtained by computing the scatterplot [31][32][33]38]. The connectivity of ecological and recreational spaces showed obvious differences in scale. The AI and CONTAG indices of ecological space both first decreased and then increased, and the CONTAG index decreased sharply at the 120 m scale and then showed a fluctuating trend. The SHDI, DIVISION and SHEI all first increased and then decreased, with high values at the 960 m scale (Figure 5a). For recreational space, except for SHEI, the other indices showed a downward trend, and SHEI showed almost no change at 30-240 m but increased gradually after the 240 m and 1800 m peaks (Figure 5b). Based on the index analysis results of the six indicators of ecological and recreational spaces corresponding to different scales, PCA was further conducted. By analyzing the statistical results of the variance contribution rate (see details from Appendix A, Table  A2), the principal component matrix of one ecological space and two recreational spaces was extracted (Table 2). Principal component 1 of ecological and recreational space had higher loading coefficients of AI, CONTAG and COHESION, which reflected the aggregation connectivity of the ecological security zone. The higher DIVISION of recreational space reflected the degree of separation and fragmentation of recreational space. This Based on the index analysis results of the six indicators of ecological and recreational spaces corresponding to different scales, PCA was further conducted. By analyzing the statistical results of the variance contribution rate (see details from Appendix A, Table A2), the principal component matrix of one ecological space and two recreational spaces was extracted ( Table 2). Principal component 1 of ecological and recreational space had higher loading coefficients of AI, CONTAG and COHESION, which reflected the aggregation connectivity of the ecological security zone. The higher DIVISION of recreational space reflected the degree of separation and fragmentation of recreational space. This was also consistent with the spatial distribution of recreational sensitivity mentioned above. In this study, the weight was calculated by the PCA method, and then the relationship between principal components and research items was calculated to obtain the "component score coefficient matrix". The expression order of the principal component function of ecological (4) and recreational (5)  Since recreation includes two principal component matrices, its comprehensive score should be calculated by multiplying the variance interpretation rate and component score [53]. The formula is as follows: 0.605 * Recreational principal component factor 1 + 0.395 * Recreational principal component factor 2 Then, standardized landscape index data of ecological space were substituted into Formula (4), and recreational space data were substituted into Formulas (5) and (6) to obtain the overall connectivity composite index at different scales [43,54]. Based on the comprehensive index of ecological and recreational space, this study selected the inflection point (120 m) of ecological space decline as the best scale for ecological source identification (Figure 6a). For recreational space, the overall declining trend was significant, so the 30 m scale was selected as the best scale for recreational space source identification (Figure 6b). obtain the overall connectivity composite index at different scales [43,54]. Based on the comprehensive index of ecological and recreational space, this study selected the inflection point (120 m) of ecological space decline as the best scale for ecological source identification (Figure 6a). For recreational space, the overall declining trend was significant, so the 30 m scale was selected as the best scale for recreational space source identification (Figure 6b).

Set the Minimum Area Threshold
The size of landscape patches will have an impact on the integrity of space, so it is necessary to determine the minimum patch area in the process of source identification [48]. In this study, by setting different area thresholds, the number of patches above the threshold was statistically analyzed to obtain a scatter diagram of the change in patch number and area ratio in ecological and recreational space. For ecological space, when the area threshold increased to 0.6 km 2 , the proportion of ecological patches and area decreased gently, indicating that although more patches were removed with the increase in area, the area of these patches was small and had little impact on the overall ecological pattern. Therefore, 0.6 km 2 was selected as the area threshold for ecological source identification (Figure 7a). For the recreational space, the decrease ratio between the number of spatial patches and area gradually decreased after 0.8 km 2 , so 0.8 km 2 was selected as the area threshold of the recreational space (Figure 7b).

Source Identification
After completing the screening of the scale and minimum area threshold, the spatial patches in this study were further screened by querying the evaluation and assignment attributes of patches, taking five points (moderately sensitive area) as the boundary, and The size of landscape patches will have an impact on the integrity of space, so it is necessary to determine the minimum patch area in the process of source identification [48]. In this study, by setting different area thresholds, the number of patches above the threshold was statistically analyzed to obtain a scatter diagram of the change in patch number and area ratio in ecological and recreational space. For ecological space, when the area threshold increased to 0.6 km 2 , the proportion of ecological patches and area decreased gently, indicating that although more patches were removed with the increase in area, the area of these patches was small and had little impact on the overall ecological pattern. Therefore, 0.6 km 2 was selected as the area threshold for ecological source identification (Figure 7a). For the recreational space, the decrease ratio between the number of spatial patches and area gradually decreased after 0.8 km 2 , so 0.8 km 2 was selected as the area threshold of the recreational space (Figure 7b). comprehensive index of ecological and recreational space, this study selected the inflection point (120 m) of ecological space decline as the best scale for ecological source identification (Figure 6a). For recreational space, the overall declining trend was significant, so the 30 m scale was selected as the best scale for recreational space source identification (Figure 6b).

Set the Minimum Area Threshold
The size of landscape patches will have an impact on the integrity of space, so it is necessary to determine the minimum patch area in the process of source identification [48]. In this study, by setting different area thresholds, the number of patches above the threshold was statistically analyzed to obtain a scatter diagram of the change in patch number and area ratio in ecological and recreational space. For ecological space, when the area threshold increased to 0.6 km 2 , the proportion of ecological patches and area decreased gently, indicating that although more patches were removed with the increase in area, the area of these patches was small and had little impact on the overall ecological pattern. Therefore, 0.6 km 2 was selected as the area threshold for ecological source identification (Figure 7a). For the recreational space, the decrease ratio between the number of spatial patches and area gradually decreased after 0.8 km 2 , so 0.8 km 2 was selected as the area threshold of the recreational space (Figure 7b).

Source Identification
After completing the screening of the scale and minimum area threshold, the spatial patches in this study were further screened by querying the evaluation and assignment attributes of patches, taking five points (moderately sensitive area) as the boundary, and

Source Identification
After completing the screening of the scale and minimum area threshold, the spatial patches in this study were further screened by querying the evaluation and assignment attributes of patches, taking five points (moderately sensitive area) as the boundary, and a total of 32 source patches were obtained, including 5 ecological source patches and 27 recreational source patches, which were divided into 2 levels of source patches through sensitivity classification. Among them, the ecological source area was 182.6 km 2 (accounting for 71.89% of the total area), and the recreational source area was 185.95 km 2 (accounting for 73.21% of the total area). Among the ecological sources, the secondary ecological sources covered a large area, while the primary sources were distributed on the slopes in the north of the region and a small amount in the areas with rich vegetation in the south (Figure 8a). The spatial distribution of recreational sources was basically consistent with the evaluation content, showing the characteristics of a clustered distribution in space (Figure 8b). Finally, this study identified the ecological/recreational source points at each level through the ArcGIS point transfer tool. The regional ecological/recreational source points were mainly concentrated in the east and south, and the source points were less distributed in the west due to their proximity to the Dujiangyan city area and more cities and towns.
(accounting for 73.21% of the total area). Among the ecological sources, the secondary ecological sources covered a large area, while the primary sources were distributed on the slopes in the north of the region and a small amount in the areas with rich vegetation in the south (Figure 8a). The spatial distribution of recreational sources was basically consistent with the evaluation content, showing the characteristics of a clustered distribution in space (Figure 8b). Finally, this study identified the ecological/recreational source points at each level through the ArcGIS point transfer tool. The regional ecological/recreational source points were mainly concentrated in the east and south, and the source points were less distributed in the west due to their proximity to the Dujiangyan city area and more cities and towns.

Identification of Ecological and Recreational Corridors
The resistance surface is an important factor that hinders the extension and diffusion of the corridor. In this study, the resistance surface was generated according to the mosaic superposition of resistance interference factors, and it was reclassified into high-, medium-and low-resistance areas. The areas of ecological high-, medium-and low-resistance areas were 58.07 (22.86%), 111.06 (43.72%) and 85.32 km 2 (33.59%), respectively. The areas of high-, medium-and low-resistance areas of recreation were 7.08 (2.79%), 64.91 (25.55%) and 182.51 km 2 (71.86%), respectively. In the ecological space, the higher resistance value was mainly distributed in the middle of the region. There are more expressways, urban expressways and construction land in the region, which greatly hinders the migration of species. Additionally, the towns in the south of the region contribute more high resistance space. In the recreational space, there are few high-resistance

Identification of Ecological and Recreational Corridors
The resistance surface is an important factor that hinders the extension and diffusion of the corridor. In this study, the resistance surface was generated according to the mosaic superposition of resistance interference factors, and it was reclassified into high-, mediumand low-resistance areas. The areas of ecological high-, medium-and low-resistance areas were 58.07 (22.86%), 111.06 (43.72%) and 85.32 km 2 (33.59%), respectively. The areas of high-, medium-and low-resistance areas of recreation were 7.08 (2.79%), 64.91 (25.55%) and 182.51 km 2 (71.86%), respectively. In the ecological space, the higher resistance value was mainly distributed in the middle of the region. There are more expressways, urban expressways and construction land in the region, which greatly hinders the migration of species. Additionally, the towns in the south of the region contribute more high resistance space. In the recreational space, there are few high-resistance areas, and only a few expressways are distributed in the middle of the region, which is due to the small number of regional expressway entrances and exits.
According to the MCR model, the minimum cost path of flow between the central points of ecological/recreational sources was calculated, that is, the corridor between the sources. There were four ecological corridors in the Dujiangyan Irrigation District, with a total length of 158.38 km. A total of 27 recreational corridors with a total length of 245.29 km were obtained, including 18 primary corridors and 9 secondary corridors. The final result integrated the main rivers in the region as a supplement to the ecological/recreational sources and the identification of corridors was completed. Overall, the regional ecological corridor avoided the main urban construction land and was mainly distributed along farmland and forestland. In terms of ecological corridor classification, one primary corridor ran through the north and south. A secondary ecological corridor was mainly located in the east of the region. In the south, because the ecological space is relatively fragile and there are few central points of ecological sources, there was no distribution of corridors (Figure 9a). The recreational corridor was not limited by the resistance of towns and roads and showed a more miscellaneous spatial distribution structure. Such distribution will increase the construction cost of the region and is not conducive to the sustainable development of space. Therefore, the corridor needs to be further optimized to make the corridor space adapt to local development (Figure 9b). total length of 158.38 km. A total of 27 recreational corridors with a total length of 245.29 km were obtained, including 18 primary corridors and 9 secondary corridors. The final result integrated the main rivers in the region as a supplement to the ecological/recreational sources and the identification of corridors was completed. Overall, the regional ecological corridor avoided the main urban construction land and was mainly distributed along farmland and forestland. In terms of ecological corridor classification, one primary corridor ran through the north and south. A secondary ecological corridor was mainly located in the east of the region. In the south, because the ecological space is relatively fragile and there are few central points of ecological sources, there was no distribution of corridors (Figure 9a). The recreational corridor was not limited by the resistance of towns and roads and showed a more miscellaneous spatial distribution structure. Such distribution will increase the construction cost of the region and is not conducive to the sustainable development of space. Therefore, the corridor needs to be further optimized to make the corridor space adapt to local development (Figure 9b).

Comprehensive Pattern of Ecology and Recreation in the Dujiangyan Irrigation District
Since the construction and development of recreational corridors mostly relies on urban and rural roads, greenways and other elements, the study adopted the proximity principle, combined the original miscellaneous recreational corridor network with the regional road network, and extended the relevant recreational corridors to form the final ecological and recreational space patterns ( Figure 10). After the comprehensive analysis

Comprehensive Pattern of Ecology and Recreation in the Dujiangyan Irrigation District
Since the construction and development of recreational corridors mostly relies on urban and rural roads, greenways and other elements, the study adopted the proximity principle, combined the original miscellaneous recreational corridor network with the regional road network, and extended the relevant recreational corridors to form the final ecological and recreational space patterns ( Figure 10). After the comprehensive analysis of the ecological and recreational patterns of the Dujiangyan Irrigation District, the recreational corridors formed a regular network distributed along the road, with a total length of 155.03 km. The ecological and recreational corridors also showed a certain degree of stacking area in space, especially in the first-class ecological corridor.
of the ecological and recreational patterns of the Dujiangyan Irrigation District, the recreational corridors formed a regular network distributed along the road, with a total length of 155.03 km. The ecological and recreational corridors also showed a certain degree of stacking area in space, especially in the first-class ecological corridor.

Analysis of Regional Landscape Patterns
Through the evaluation and construction of ESPs in the ecological and recreational spaces of the Dujiangyan Irrigation District, the evaluation results showed that the ecological level of the region was relatively low, and the average value was downstream of the five-level classification, which was 3.36 (positive distribution). The grade distribution of recreational space of 2.97 (negative distribution) was high. Therefore, in the existing space of rapid development, how to protect and repair the ecological environment and improve the overall ecological environment quality, especially in cities and towns and around roads and water system corridors, is an urgent problem to be solved in regional development. There were also a large number of expressways and urban truck roads running in the east-west direction of the region, which reduced spatial connectivity, seriously breaking up the ecological space on both sides of the road, and increasing the resistance to the ecological process [55][56][57].
Using the MCR model, this paper determined the ecological and recreational corridors, which are an important part of landscape patterns, and obtained the comprehensive spatial pattern of the combination of points, lines and surfaces. By sorting the ecological and recreational corridors and the relationship between nodes and elements in space, we proposed the ecological spatial pattern of "three axes, two belts and four cores" and the

Analysis of Regional Landscape Patterns
Through the evaluation and construction of ESPs in the ecological and recreational spaces of the Dujiangyan Irrigation District, the evaluation results showed that the ecological level of the region was relatively low, and the average value was downstream of the five-level classification, which was 3.36 (positive distribution). The grade distribution of recreational space of 2.97 (negative distribution) was high. Therefore, in the existing space of rapid development, how to protect and repair the ecological environment and improve the overall ecological environment quality, especially in cities and towns and around roads and water system corridors, is an urgent problem to be solved in regional development. There were also a large number of expressways and urban truck roads running in the east-west direction of the region, which reduced spatial connectivity, seriously breaking up the ecological space on both sides of the road, and increasing the resistance to the ecological process [55][56][57].
Using the MCR model, this paper determined the ecological and recreational corridors, which are an important part of landscape patterns, and obtained the comprehensive spatial pattern of the combination of points, lines and surfaces. By sorting the ecological and recreational corridors and the relationship between nodes and elements in space, we proposed the ecological spatial pattern of "three axes, two belts and four cores" and the recreational spatial pattern of "four axes and multiple points". Specifically, the primary ecological axis runs through the site from northeast to south, maintaining the ecological connectivity between cultivated land and forestland. Additionally, the secondary axis extends from the north to the middle and runs through the place with large ecological source coverage in the south, ensuring the multidirectional continuation of the regional ecological space. The recreational main axis connects the north along the road to the west Land 2022, 11, 383 15 of 21 of the region, forming an important route for regional tourism, while the recreational secondary axis is composed of two lines running from east to west in the north and one line running from north to south in the south, ensuring the integrity of the connection between regional rural tourism and scenic resources (Figure 11a). recreational spatial pattern of "four axes and multiple points". Specifically, the primary ecological axis runs through the site from northeast to south, maintaining the ecological connectivity between cultivated land and forestland. Additionally, the secondary axis extends from the north to the middle and runs through the place with large ecological source coverage in the south, ensuring the multidirectional continuation of the regional ecological space. The recreational main axis connects the north along the road to the west of the region, forming an important route for regional tourism, while the recreational secondary axis is composed of two lines running from east to west in the north and one line running from north to south in the south, ensuring the integrity of the connection between regional rural tourism and scenic resources (Figure 11a).

Optimization of Multiscale Landscape Connectivity
For source identification, in the past, the identification of ecological source sites was mostly used to directly select patches with large areas and good habitat quality as the source sites, or the traditional particle size analysis method was usually used for screening. There was little consideration of the scale problem and some subjectivity was involved [58,59]. This paper used the three-step source identification method of "connectivity index analysis; patch threshold determination; evaluation result classification", which greatly improved the objectivity and accuracy of source identification. This study used the landscape connectivity indices of CONHESION, AI, CONTAG, SHDI, SHEI, and DIVISION to obtain the best research scale. The scale of 120 m ecological space and 30 m recreational space selected in the text were used as the scales of source identification. The two-scale spaces had good landscape connectivity, and the comprehensive

Optimization of Multiscale Landscape Connectivity
For source identification, in the past, the identification of ecological source sites was mostly used to directly select patches with large areas and good habitat quality as the source sites, or the traditional particle size analysis method was usually used for screening. There was little consideration of the scale problem and some subjectivity was involved [58,59]. This paper used the three-step source identification method of "connectivity index analysis; patch threshold determination; evaluation result classification", which greatly improved the objectivity and accuracy of source identification. This study used the landscape connectivity indices of CONHESION, AI, CONTAG, SHDI, SHEI, and DIVISION to obtain the best research scale. The scale of 120 m ecological space and 30 m recreational space selected in the text were used as the scales of source identification. The two-scale spaces had good landscape connectivity, and the comprehensive connectivity value reached 0.8, indicating that the regional ecology and recreational sources had been best displayed on the spatial scale, which had a good relevance for guiding the construction of regional ecological space. Subsequent research can make the fitting degree of the final, overall, connectivity comprehensive index higher through the selection of more scales and indices to improve the accuracy of the final scale selection results.

Effect of Patch Size on Source Recognition
The determination of patch scale has a great impact on the number of patches and the functions carried by a single patch area, and the determination of its area should adapt to the local space [60]. Most of the relevant studies are aimed at cities and other such administrative spatial units. The determination of the spatial patch area was mostly in the range of 0.1-10 km 2 , but the total area of the study area in this paper was only 254 km 2 . After the area threshold analysis, the patch area scale of an ecological source area of 0.6 km 2 and a recreational source area of 0.8 km 2 could more accurately retain the source patch [61,62]. Under such a scale, follow-up urban and rural planning and construction can, according to such a scale, be used to carry out ecological restoration projects of ecological space and ecological sources, integrate tourism resources and reasonably develop rural tourism to realize the protection of regional ecological space and the development of recreational space to improve the level of regional ecological security. However, the final identification result of this paper is that there were only five patches in the ecological source, which was significant compared with the number of recreational patches. In a follow-up study, more screening conditions can be added to reduce the difference in the patch area identification of multiple indicators and improve the balance between the data.

Optimization of the Multifactor Resistance Surface Construction Method
For the identification of resistance surfaces, relevant studies often use a single type of land-use data to construct the resistance surfaces [51,52]. This leads to the separation of the resistance surface from the real space. Therefore, the geographical factors, water areas and road factors in this study strengthen the connection with the real space to a certain extent. However, there are still some problems in this study. The most important problem is that the road cannot fully reflect the impact of rich human activities in real space on the flow process of material and energy. Therefore, follow-up research should consider more social development factors, such as socio-economic differences and functional spatial differences in different regions, and should realize the complete expression of resistance factors in real space.

Implications to the Development of the Dujiangyan Irrigation District
By analyzing the ecological corridor planning mentioned in the relevant upper planning, the establishment of an ecological corridor in the upper planning was based on the river corridor. The similarity between the two landscape patterns lies in the addition of a water space corridor with the water system as the main body, which reflects the water network characteristics of the Dujiangyan Irrigation District. It also shows similarity in a small number of slope areas in the northern part of the area. The difference is that the corridor in the upper planning involves fewer elements and shows a relatively single ecological space, while our research results integrate factors such as land use, development and construction interference in the evaluation stage, reflecting the complexity of space under the influence of multiple factors. Therefore, our research results can be used as a reference and can supplement the construction of regional ecological and recreational space ( Figure 11).
In the follow-up development of the Dujiangyan Irrigation District, it is necessary to propose a new long-term optimization plan for landscape space, emphasizing forwardlooking protection and utilization. In terms of specific improvement content, the protection and restoration of the source area should be considered first, with cultivated land, forestland and the water area as matrix elements for protection and restoration, to improve the quality and coverage rate of the ecological source area. At the same time, the fragmented ecological land should be integrated to improve the stability of the ecosystem through scale change. Second, the connectivity between the source areas should be improved, green space construction along the corridor and the connectivity on both sides of the road should be strengthened, and the ecological environment quality should be improved. Finally, for the areas requiring recreational development, the construction content should be strictly controlled to reduce interference with the ecological environment and the coordinated development of ecological environment protection and development should be realized.

Conclusions
Taking Dujiangyan Irrigation District as an example, this study constructed regional comprehensive landscape patterns through ESPs from two dimensions of ecology and recreation. Based on the multiscale connectivity analysis and area threshold identification, the source area was constructed and the resistance surface was resolved to obtain the "eco-recreation" integrated landscape network based on the corridor.
The results show that (1) The development of regional recreation space has a good foundation, but the average value of the ecological space grade is only 3.36, which is low. In subsequent regional development, special attention should be paid to the construction of ecological landscapes in towns, roads, water corridors and other low-grade areas; (2) The study identified 120 and 30 m as the most suitable scales for ecological and recreational space and 0.6 and 0.8 km 2 as the minimum area threshold of ecological and recreational space, which had a great promotion effect on the accuracy of source identification. At the same time, the construction of resistance surfaces can improve the accuracy of resistance surfaces to a certain extent by integrating roads, water systems and topographies; and (3) By constructing ecological and recreational networks, this study identified four ecological corridors with a total length of 158.38 km in the Dujiangyan Irrigation District and 27 recreational corridors with a total length of 245.29 km, which can be used as an important space for subsequent regional development.
Based on the ecological network analysis, this study suggested that in the future ecological space planning of the region, priority should be given to ecological land consolidation, integration of scattered ecological sources and delineation of protected areas. Meanwhile, ecological corridors should be planned to improve the connectivity of ecological space, so as to ensure the ecological security pattern. In terms of recreation space, this study suggests that appropriate development should be carried out according to the recreation corridor, under the premise of maintaining ecological protection, so as to better utilize regional recreation resources and promote the coordinated development of space protection and utilization.
Appendix A Table A1. Grade assignment and spatial evaluation weight. The density of traffic and villages is low, with a population of 2000-3000 7

Criterion Layer (A) Sub-Criterion Layer (B) Weight Factor Layer (C) Conditions for Reclassification Value Weight
The density of traffic, villages and towns is moderate, 3000 < population ≤ 4000 5 The density of traffic, villages and towns is very high, and the population is more than 5000