Construction method of high-density urban area ecological-recreation function composite corridor network

Abstract

The invention belongs to the field of urban ecological corridors, and discloses a construction method of a high-density urban area ecological-recreation function composite corridor network. The method comprises the following steps: inputting first phenological data and a first environmental factor into a maximum entropy model, and outputting grids of a suitable-to-grow area source ground of a target bird; adopting a natural breakpoint method to re-classify the grids of the source land of the suitable-for-growth area, and determining the grids as target ecological source land; determining a second environment factor greater than a preset correlation coefficient threshold and a contribution rate of the second environment factor by adopting a cutting method, and determining a landscape resistance surface of the target ecological source land based on the second environment factor and the contribution rate; calling a map software API interface, and screening out a target recreation space; constructing a cost grid surface of the target recreation space by adopting a spatial principal component analysis method and a cost factor; and determining an ecological corridor and a recreation type corridor by adopting a minimum path cost method, a landscape resistance surface and a cost grid surface. The urban ecological corridor with harmonious symbiosis of human and nature can be constructed.

Description

Technical Field

[0001] The present invention relates to the technical field of urban ecological corridors, and in particular to a method for constructing a high-density urban ecological-recreation function composite corridor network. Background Art

[0002] In recent years, the level of urbanization has achieved a leapfrog growth, which has promoted economic development. However, in the process of urban expansion, it has inevitably led to habitat fragmentation, biodiversity reduction, and reduced ecological service functions, which are also the main ecological problems brought about by the urbanization process. At the same time, the fragmentation of green space has made urban space more indifferent, leading to the indifference of social relationships and the isolation of social classes. Modern urban residents have a strong desire to find closeness to nature and a quiet pastoral life in the jungle of reinforced concrete.

[0003] In order to solve the above problems, urban ecological corridors can be constructed and the spatial layout of urban green spaces can be optimized. Compared with the planning of ecological corridors in natural areas, the construction of urban ecological corridors faces challenges such as complex land use, high intensity of human activity interference and diversified use needs. It shows strong construction limitations, functional complexity and diversity of behavior patterns. It is difficult to fully use the traditional ecological corridor planning theory and methods, especially in high-density cities with scarce land resources. Its construction is gradually transformed into a multi-objective and multi-functional model. Summary of the invention

[0004] Based on this, it is necessary to provide a method for constructing a high-density urban ecological-recreational function composite corridor network in response to the above technical problems.

[0005] In a first aspect, a method for constructing a high-density urban ecological-recreational function composite corridor network is provided, the method comprising:

[0006] Inputting a plurality of first phenological data and a plurality of first environmental factors into a pre-trained maximum entropy model, and outputting a grid of a source of a suitable habitat for target birds;

[0007] The grid of the source of the suitable habitat is reclassified by using the natural breakpoint method to obtain the candidate ecological source suitable for the target birds, and the candidate ecological source with an area greater than a preset threshold is determined as the target ecological source suitable for the target birds;

[0008] The jackknife method is used to determine a plurality of second environmental factors corresponding to the target birds and greater than a preset correlation coefficient threshold and a contribution rate of the plurality of second environmental factors, and based on the plurality of second environmental factors and the corresponding contribution rates, the landscape resistance surface of the target ecological source is determined;

[0009] Call the map software API application program interface to capture preset recreation keywords and screen out target recreation spaces that meet the preset recreation function requirements;

[0010] The spatial principal component analysis method is used to construct the cost grid surface of the target recreation space based on the cost factors of the target recreation space; the cost factors include the node service range, spatial construction potential, residents' preferred path and resource attractiveness;

[0011] The minimum path cost method is adopted to determine the ecological corridor of the target birds and the recreational corridor of humans based on the landscape resistance surface of the target ecological source and the cost grid surface of the target recreational space.

[0012] As an optional implementation, before inputting the plurality of first phenological data and the plurality of first environmental factors into a pre-trained maximum entropy model and outputting the grid of the source of the suitable habitat of the target bird, the method further includes:

[0013] A correlation analysis method is used to determine the correlation coefficient values ​​between multiple phenological data and multiple environmental factors, and multiple first phenological data and multiple first environmental factors that are less than the preset correlation coefficient threshold are obtained.

[0014] As an optional implementation manner, determining the landscape resistance surface of the target ecological source area according to the plurality of the second environmental factors and the corresponding contribution rates includes:

[0015] The knife-cut method is used to determine the second environmental factor and its corresponding contribution rate, and multiple second environmental factors whose contribution rates are greater than a preset contribution rate threshold are determined as multiple key factors constituting the landscape resistance surface;

[0016] A linear function obtained by adding the products of the plurality of key factors and the corresponding contribution rates is determined as the landscape resistance surface of the target ecological source area.

[0017] As an optional implementation, the spatial construction potential includes road space conversion potential, water system space conversion potential and linear space utilization potential, the residents' preferred paths include walking preference paths and cycling preference paths, and the spatial principal component analysis method is used to construct the cost grid surface of the target recreation space based on the cost factor of the target recreation space, including:

[0018] According to a plurality of preset distance service thresholds, assigning a node service range of each recreation node of the target recreation space, and determining a weight of the node service range based on the spatial principal component analysis method;

[0019] Assigning values ​​to the road space conversion potential of the target recreational space according to multiple road grades, and determining weights of the road space conversion potential based on the spatial principal component analysis method;

[0020] According to the distance threshold between each pixel and the water system space in each recreation node, the water system space conversion potential of the target recreation space is assigned a value, and based on the spatial principal component analysis method, the weight of the water system space conversion potential is determined;

[0021] Based on the color satellite image, the strip green space and the remaining other land in the target recreational space are identified, and based on the strip green space and the remaining other land, the linear space utilization potential of the target recreational space is assigned a value, and based on the spatial principal component analysis method, the weight of the linear space utilization potential is determined;

[0022] Assigning values ​​to the walking preference paths of the target recreational space according to the walking frequency of the residents' walking preference paths, and determining weights of the walking preference paths based on the spatial principal component analysis method;

[0023] Assigning values ​​to the cycling preference paths of the target recreational space according to the cycling frequency of the residents' cycling preference paths, and determining weights of the cycling preference paths based on the spatial principal component analysis method;

[0024] Assigning a value to the resource attraction of the target recreational space according to a predetermined level of resource attraction, and determining a weight of the resource attraction based on the spatial principal component analysis method;

[0025] The sum of the product of the value assigned to the node service range of the target recreational space and the weight of the node service range, the product of the value assigned to the road space conversion potential and the weight of the road space conversion potential, the product of the value assigned to the water system space conversion potential and the weight of the water system space conversion potential, the product of the value assigned to the linear space utilization potential and the weight of the linear space utilization potential, the product of the value assigned to the walking preference path and the weight of the walking preference path, the product of the value assigned to the cycling preference path and the weight of the cycling preference path, and the product of the value assigned to the resource attraction and the weight of the resource attraction is determined as the cost grid surface of the target recreational space.

[0026] As an optional implementation, the formula for determining the resource attractiveness of the target recreational space is:

[0027]

[0028] Where f(x) represents the resource attraction at the spatial location x, r represents the path distance attenuation threshold, n represents the number of facilities whose path distance to location x is less than or equal to r, and d ix represents the shortest distance from point i to facility point x, and k represents the spatial weight function.

[0029] As an optional implementation, the minimum path cost method is used to determine the ecological corridor of the target birds and the recreational corridor of humans based on the landscape resistance surface of the target ecological source and the cost grid surface of the target recreational space, including:

[0030] The minimum path cost method is used to determine the minimum cost path between the landscape resistance surfaces of any two ecological nodes in the target ecological source area, and the minimum cost path is determined as the ecological corridor of the target bird;

[0031] The minimum path cost method is adopted to determine the minimum cost path between the cost grid surfaces of any two recreation nodes in each target recreation space, and the minimum cost path is determined as the recreational corridor for the human beings.

[0032] As an optional implementation, the method further includes:

[0033] Obtaining the bird diversity index of each ecological node in the target ecological source area;

[0034] Based on the bird diversity index, each of the ecological nodes is graded to obtain multiple levels of ecological nodes of each of the ecological nodes;

[0035] Based on the preset evaluation system and random forest classifier, the recreation nodes of the target recreation space are classified to obtain recreation nodes of multiple levels.

[0036] As an optional implementation, the method further includes:

[0037] Perform spatial overlay analysis on multiple ecological nodes of the above mentioned level and multiple recreational nodes of the above mentioned level to obtain multiple composite nodes and functional use levels of multiple composite nodes;

[0038] A spatial overlay analysis is performed on the ecological corridors and the recreational corridors to extract urban ecological corridors with complex functions. The urban ecological corridors with complex functions are based on the ecological corridors. The urban ecological corridors with complex functions are ecological corridors and recreational corridors connecting the same two nodes, ecological corridors and recreational corridors with adjacent spatial nodes, and ecological corridors and recreational corridors that run through the same residential area.

[0039] As an optional implementation, the method further includes:

[0040] Based on circuit theory, the ecological nodes of the target ecological source and the landscape resistance surface are input into the Pinchpoint Mapper and Barrier Mapper modules respectively, and the pinch points / obstacle spaces in the complex-functional urban ecological corridor and the pinch points / obstacle spaces in the ecological corridor are extracted respectively. The pinch points / obstacle spaces in the urban ecological corridor are used as the primary point protection / restoration space, and the pinch points / obstacle spaces in the ecological corridor are used as the secondary point protection / restoration space.

[0041] As an optional implementation manner, the suitable zone source includes a high suitable zone source, a medium suitable zone source, a low suitable zone source and a non-suitable zone source.

[0042] The present application provides a method for constructing a composite corridor network of ecological and recreational functions in a high-density urban area. The technical solution provided by the embodiments of the present application brings at least the following beneficial effects: by determining the target ecological source of target birds and the target recreational space of humans, constructing a landscape resistance surface and a cost grid surface, and based on the landscape resistance surface of the target ecological source and the cost grid surface of the target recreational space, determining the ecological corridor of the target birds and the recreational corridor of humans. Performing spatial superposition and overlap analysis on the ecological corridor and the recreational corridor, an urban ecological corridor is obtained. In this way, the living habits of organisms and the recreational needs of humans are integrated into the source selection process of the ecological and recreational corridors, which improves the accuracy and practicality of the source selection and meets the main function of species protection. It is of great practical significance to develop a composite functional urban ecological corridor construction method that takes into account both biodiversity protection and residents' recreational needs, and to build an urban ecological corridor where people and nature coexist in harmony.

[0043] It should be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0045] Figure 1 A flow chart of a method for constructing a high-density urban ecological-recreational function composite corridor network provided in an embodiment of the present application;

[0046] Figure 2 A schematic diagram of the distribution of a suitable habitat space and a target ecological source provided in an embodiment of the present application;

[0047] Figure 3 Aj is a schematic diagram of a grid diagram of an environmental factor provided in an embodiment of the present application;

[0048] Figure 3 numeral k is a schematic diagram of a comprehensive resistance surface provided in an embodiment of the present application;

[0049] Figure 4 A schematic diagram of the distribution of recreational nodes provided in an embodiment of the present application;

[0050] Figure 5 AG is a schematic diagram of a grid diagram of a cost factor provided in an embodiment of the present application;

[0051] Figure 5 h is a schematic diagram of a cost grid surface provided in an embodiment of the present application;

[0052] Figure 6 A schematic diagram of the distribution of an ecological corridor provided in an embodiment of the present application;

[0053] Figure 7 A schematic diagram of the distribution of a recreational corridor provided in an embodiment of the present application;

[0054] Figure 8 A schematic diagram of a multi-level ecological node provided in an embodiment of the present application;

[0055] Fig. 9 A schematic diagram of the classification of recreational nodes provided in an embodiment of the present application;

[0056] Fig.10 A schematic diagram of the classification of a composite source provided in an embodiment of the present application;

[0057] Fig.11 Figure a is a schematic diagram of the distribution of a composite functional urban ecological corridor provided in an embodiment of the present application;

[0058] Fig.11 bd in the figure is a schematic diagram of a composite type of a composite functional urban ecological corridor provided in an embodiment of the present application;

[0059] Fig.12 A schematic diagram of the distribution of a pinch point space provided in an embodiment of the present application;

[0060] Fig.13 A schematic diagram of the distribution of an obstacle space provided in an embodiment of the present application;

[0061] Fig.14 A schematic diagram of a method for constructing a high-density urban ecological-recreational function composite corridor network provided in an embodiment of the present application. DETAILED DESCRIPTION

[0062] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0063] The following will describe in detail a method for constructing a high-density urban ecological-recreation function composite corridor network provided by an embodiment of the present application in combination with a specific implementation method. Figure 1 A flowchart of a method for constructing a high-density urban ecological-recreation function composite corridor network provided in an embodiment of the present application, such as Figure 1 As shown, the specific steps are as follows:

[0064] Step 101, input a plurality of first phenological data and a plurality of first environmental factors into a pre-trained maximum entropy model, and output a grid of the source of the suitable habitat of the target bird.

[0065] In implementation, birds and humans can be selected as the target objects for the construction of urban ecological corridors. Before constructing urban ecological corridors, it is necessary to first determine the target ecological source of the target birds. Therefore, multiple first phenological data and multiple first environmental factors can be input into the pre-trained maximum entropy model to output the grid of the source of the target bird's suitable habitat. In the source selection process, birds and leisure people are used as target species, and the living habits of birds and human recreational needs are integrated into the source selection process of ecological and recreational corridors, aiming to improve the accuracy and practicality of source selection.

[0066] Before executing step 101, it is necessary to first determine multiple first phenological data and multiple first environmental factors. The specific steps are: use the correlation analysis method to determine the correlation coefficient values ​​between multiple phenological data and multiple environmental factors, and obtain multiple first phenological data and multiple first environmental factors that are less than a preset correlation coefficient threshold.

[0067] In implementation, 19 phenological data (bio_1-19), normalized difference vegetation index (NDVI), distance to road (DR), distance to water system (DW), vegetation canopy height (CH), building proportion (AP), night light index (NL) can be calculated through Spearman correlation analysis. When the absolute value of the correlation coefficient is between 0.1-0.3, 0.3-0.5 and greater than 0.5, it is considered that there is a weak correlation, moderate correlation and strong correlation between variables such as phenological data and environmental factors. Phenological data and environmental factors with a correlation coefficient lower than 0.5, that is, no significant correlation with other variables, are screened out, and the grid of the source of the suitable habitat of the target birds is determined through the maximum entropy MaxEnt model.

[0068] Step 102, using the natural breakpoint method, reclassify the grid of the source of the suitable habitat to obtain the alternative ecological source suitable for the target birds, and determine the alternative ecological source with an area greater than a preset area threshold as the target ecological source suitable for the target birds.

[0069] In the implementation, the natural breakpoint method is used to reclassify the grid of the source area of ​​the suitable habitat, and the source areas of the suitable habitat include high-suitability source areas, medium-suitability source areas, low-suitability source areas, and non-suitability source areas. High-suitability source areas and medium-suitability source areas are selected as alternative ecological source areas for the survival of target birds. For example, the natural breakpoint method of ArcGIS10.2 can be used for reclassification. Multiple alternative ecological source areas are compared with the preset area thresholds respectively, and the alternative ecological source areas with an area greater than the preset area threshold are determined as target ecological source areas for the survival of target birds. Figure 2 As shown, Figure 2 Including the study area boundaries and target ecological source areas of urban ecological corridors.

[0070] Step 103, using the jackknife method to determine multiple second environmental factors corresponding to the target birds that are greater than a preset correlation coefficient threshold and the contribution rates of the multiple second environmental factors, and based on the multiple second environmental factors and the corresponding contribution rates, determine the landscape resistance surface of the target ecological source.

[0071] In the implementation, the jackknife method is used to evaluate the environmental variable factors affecting the target birds, their contribution rate and replacement importance, and the key factors constituting the landscape resistance surface are screened out according to the importance of the environmental factors, and then a linear landscape resistance surface is constructed. The environmental factors affecting the target birds can be used as the second environmental factor. The correlation of the second environmental factor is greater than the preset correlation coefficient threshold. Figure 3 As shown, Figure 3 It includes the grid map of environmental factors and the comprehensive resistance surface, among which the comprehensive resistance surface is the landscape resistance surface.

[0072] Specifically, the specific steps of determining the landscape resistance surface of the target ecological source area based on multiple second environmental factors and corresponding contribution rates are as follows:

[0073] Step 1: Use the knife-cut method to determine the second environmental factor and its corresponding contribution rate, and determine multiple second environmental factors with contribution rates greater than a preset contribution rate threshold as multiple key factors constituting the landscape resistance surface.

[0074] In implementation, the jackknife method is used to determine the second environmental factor and its corresponding contribution rate. The contribution rate is compared with a preset contribution rate threshold, and multiple second environmental factors with contribution rates greater than the preset contribution rate threshold are determined as multiple key factors constituting the landscape resistance surface. The determination of the second environmental factor and its corresponding contribution rate using the jackknife method can be set and determined in the maximum entropy model.

[0075] Step 2: The linear function of the sum of the products of multiple key factors and corresponding contribution rates is determined as the landscape resistance surface of the target ecological source area.

[0076] In the implementation, the linear function of the product of multiple key factors and the corresponding contribution rate is determined as the landscape resistance surface of the target ecological source. The contribution degree is the weight, and the second environmental factor multiplied by the weight forms the final landscape resistance surface. In the process of resistance surface construction, the MaxEnt model and spatial principal component analysis technology are used to assign values ​​to resistance factors and construct landscape resistance surfaces to improve the objectivity of weight assignment.

[0077] Step 104, calling the map software API application program interface, grabbing preset recreation keywords, and screening out target recreation spaces that meet the preset recreation function requirements.

[0078] During implementation, the map software API (Application Program Interface) interface can be called to capture preset recreational keywords and screen out urban green spaces with recreational service functions. Based on the preset recreational function requirements, the target recreational space is screened out from the urban green space. Among them, urban green spaces are divided into park green spaces, production green spaces, protective green spaces, ancillary green spaces and other green spaces. Functional facility green spaces in park green spaces and ancillary green spaces, as well as scenic spots, nature reserves, scenic forests and the like in other green spaces are all important types of urban green spaces that can meet the leisure needs of urban residents, and are also the core objects for the selection of target recreational spaces in this application. For example, the API interface of Amap can be called through Python to capture data using "tourist attractions" and "leisure and entertainment" as preset recreational keywords, and to screen out urban green spaces with recreational service functions, such as Figure 4 As shown in the figure, it includes multiple recreation nodes. By referring to authoritative rankings such as Ctrip's reputation list and Dianping's must-go list, the preset recreational function requirements are that the number of comments is more than 300, the types are diverse, the local characteristics are highly popular, and the urban green space is used as the target recreation space. You can also refer to authoritative rankings such as Ctrip's reputation list and Dianping's must-go list to screen and consider urban green spaces based on the number of comments, characteristics, etc. to obtain the target recreation space.

[0079] Step 105, using the spatial principal component analysis method, based on the cost factors of the target recreation space, construct a cost grid surface of the target recreation space; the cost factors include node service range, spatial construction potential, residents' preferred path and resource attractiveness.

[0080] In the implementation, the spatial principal component analysis method SPCA (spatial principal component analysis) is used to objectively assign values ​​to cost factors and determine the weights of cost factors, so that the cost grid surface can be constructed. The cost factors include four dimensions: node service range, spatial construction potential, residents' preferred path, and resource attractiveness. In the process of building the resistance surface, the MaxEnt model and spatial principal component analysis technology are used to assign values ​​to resistance factors and construct the cost grid surface to improve the objectivity of weight assignment.

[0081] Specifically, the spatial construction potential includes road space conversion potential, water system space conversion potential and linear space utilization potential, and the residents' preferred paths include walking preferred paths and cycling preferred paths. The specific steps for executing step 106 are as follows:

[0082] Step three, assign a value to the node service range of each recreation node in the target recreation space according to multiple preset distance service thresholds, and determine the weight of the node service range based on the spatial principal component analysis method.

[0083] In implementation, the node service range of each recreational node in the target recreational space is assigned according to multiple preset distance service thresholds. Then, based on the spatial principal component analysis method, the weight of the node service range is determined. For example, the node service range selects 500, 1350, 3000 and 5000m as the preset distance service thresholds, which represent the distances of 5min walk, 15min walk or 5min bicycle, 10min bicycle respectively. The Euclidean distance tool of ArcGIS software can also be used to analyze the service range of recreational nodes. For example, the node service range of less than 500m is assigned a value of 10, the node service range of 500m-1350m is assigned a value of 20, the node service range of 1350m-3000m is assigned a value of 30, the node service range of 3000m-5000m is assigned a value of 40, and the node service range of more than 5000m is assigned a value of 50. The weight of the node service range is 0.056. Figure 5 As shown, a grid diagram is obtained in which the cost factor is the service range of the node.

[0084] Step 4: Assign values ​​to the road space conversion potential of the target recreational space according to multiple road grades, and determine the weight of the road space conversion potential based on the spatial principal component analysis method.

[0085] In implementation, spatial construction potential includes road space conversion potential, water system space conversion potential and linear space utilization potential. First, determine the assignment and weight of road space conversion potential. According to multiple road grades, the road space conversion potential of the target recreational space is graded and assigned, and the weight of the road space conversion potential is determined based on the spatial principal component analysis method. Among them, the multiple road grades are non-road areas, expressways, main roads, secondary roads and branch roads. The value assigned to non-road areas is 1000, the value assigned to expressways is 40, the value assigned to main roads is 30, the value assigned to secondary roads is 20, and the value assigned to branch roads is 10. The weight of the road space conversion potential is determined to be 0.216. Figure 5 As shown, the cost factor is a grid map of the road space conversion potential.

[0086] Step five: According to the distance threshold from each pixel to the water system space in each recreational node, the water system space transformation potential of the target recreational space is assigned, and the weight of the water system space transformation potential is determined based on the spatial principal component analysis method.

[0087] In implementation, spatial construction potential includes road space transformation potential, water system space transformation potential and linear space utilization potential. Then determine the assignment and weight of the water system space transformation potential. According to the distance threshold from each pixel to the water system space in each recreational node, the water system space transformation potential of the target recreational space is graded and assigned, and then the weight of the water system space transformation potential is determined based on the spatial principal component analysis method. For example, the distance from the pixel to the water system space can be 0-50m, 50m-100m, 100m-300m, 300m-600m, and greater than 600m. Assign values ​​according to different distances, 0-50m is assigned a value of 10, 50m-100m is assigned a value of 20, 100m-300m is assigned a value of 30, 300m-600m is assigned a value of 40, and greater than 600m is assigned a value of 50. The weight of the water system space transformation potential is 0.278. Figure 5 As shown, the cost factor is a grid map of the spatial transformation potential of the water system.

[0088] Step six: Identify the strip green space and other remaining land in the target recreational space based on color satellite images, assign values ​​to the linear space utilization potential of the target recreational space based on the strip green space and other remaining land, and determine the weight of the linear space utilization potential based on the spatial principal component analysis method.

[0089] In implementation, spatial construction potential includes road space conversion potential, water system space conversion potential and linear space utilization potential. Finally, the assignment and weight of linear space utilization potential are determined. Based on color satellite images, the strip green space and remaining other land in the target recreational space are identified, and based on the strip green space and remaining other land, the linear space utilization potential of the target recreational space is assigned. Based on the spatial principal component analysis method, the weight of linear space utilization potential is determined. For example, ENVI software is used to visually interpret the high-resolution true color satellite images of the study area to identify the strip green space and remaining other land in the city. The assignment of the strip green space is 10, and the assignment of the remaining other land can be 500. The weight of the linear space utilization potential is determined to be 0.103. Figure 5 As shown, the resulting grid map shows that the cost factor is a linear space utilization potential.

[0090] Step seven: assign values ​​to the preferred walking paths of the target recreational space according to the walking frequency of residents' preferred walking paths, and determine the weights of the preferred walking paths based on the spatial principal component analysis method.

[0091] In implementation, residents' preferred paths include walking preference paths and cycling preference paths. First, determine the assignment and weight of walking preference paths. Collect the walking frequency of users' walking preference paths, and also classify and assign values ​​to the walking preference paths of the target recreational space through the walking routes recommended in the Keep application. Then, based on the spatial principal component analysis method, determine the weight of the walking preference path. Among them, the point-to-line function of ArcGIS software can also be used and vectorized. Further, spatial rasterization can be performed, and the natural breakpoint classification method can be used to classify it and assign resistance and weight. For example, walking preference paths can be divided into high-frequency use sections, higher-frequency use sections, medium-frequency use sections, low-frequency use sections and other sections according to walking frequency. The assignment value of high-frequency use sections is 10, the assignment value of higher-frequency use sections is 20, the assignment value of medium-frequency use sections is 30, the assignment value of low-frequency use sections is 40, and the assignment value of other sections is 1000. The weight of the walking preference path is determined to be 0.127. For example Figure 5 As shown, the grid map of the cost factor for the walking preference path is obtained.

[0092] Step 8: Assign values ​​to the target recreational space’s cycling preference paths according to the cycling frequency of residents’ cycling preference paths, and determine the weights of the cycling preference paths based on the spatial principal component analysis method.

[0093] In implementation, residents' preferred paths include walking preference paths and cycling preference paths. Then, determine the assignment and weight of the cycling preference paths. Collect the riding frequency of users' cycling preference paths, and also classify and assign values ​​to the cycling preference paths of the target recreational space through the cycling routes recommended in the Keep application. Then, based on the spatial principal component analysis method, determine the weight of the cycling preference path. Among them, the point-to-line function of ArcGIS software can also be used and vectorized, and further, spatial rasterization can be performed, and the natural breakpoint classification method can be used to classify it and assign resistance and weight. For example, the cycling preference path can be divided into high-frequency use sections, higher-frequency use sections, medium-frequency use sections, low-frequency use sections and other sections according to the frequency of riding. The assignment value of high-frequency use sections is 10, the assignment value of higher-frequency use sections is 20, the assignment value of medium-frequency use sections is 30, the assignment value of low-frequency use sections is 40, and the assignment value of other sections is 1000. Determine the weight of the cycling preference path to be 0.195. For example Figure 5 As shown, the cost factor is a grid map of the cycling preference path.

[0094] Step nine: assign values ​​to the resource attractiveness of the target recreational space according to the predetermined level of resource attractiveness, and determine the weight of the resource attractiveness based on the spatial principal component analysis method.

[0095] In implementation, resource attractiveness can be expressed through functional density. The denser the distribution of POI (Point of Interest), the higher the overall attraction value of the plot. The resource attractiveness is determined by the kernel density analysis method. The resource attractiveness of the target recreational space can be graded and assigned according to the predetermined level of resource attractiveness, and the weight of resource attractiveness can be determined based on the spatial principal component analysis method. For example, the level of resource attractiveness can be high, relatively high, medium, relatively low and low. Among them, the high value is assigned to 10, the relatively high value is assigned to 20, the medium value is assigned to 30, the relatively low value is assigned to 40 and the low value is assigned to 50. The weight of resource attractiveness is determined to be 0.025. For example Figure 5 As shown, a grid diagram is obtained in which the cost factor is the resource attractiveness.

[0096] Furthermore, the formula for determining the resource attractiveness of the target recreational space is:

[0097]

[0098] Where f(x) represents the resource attraction at the spatial location x, r represents the path distance attenuation threshold, n represents the number of facilities whose path distance to location x is less than or equal to r, and d ix represents the shortest distance from point i to facility point x, and k represents the spatial weight function.

[0099] Step ten, determine the sum of the product of the assigned value of the node service range of the target recreational space and the weight of the node service range, the product of the assigned value of the road space conversion potential and the weight of the road space conversion potential, the product of the assigned value of the water system space conversion potential and the weight of the water system space conversion potential, the product of the assigned value of the linear space utilization potential and the weight of the linear space utilization potential, the product of the assigned value of the walking preference path and the weight of the walking preference path, the product of the assigned value of the cycling preference path and the weight of the cycling preference path, and the product of the assigned value of the resource attraction and the weight of the resource attraction as the cost grid surface of the target recreational space.

[0100] In implementation, the values ​​of various cost factors are weighted and superimposed to form a cost grid map. That is to say, the product of the value of the node service range of the target recreational space and the weight of the node service range, the value of the road space conversion potential and the weight of the road space conversion potential, the value of the water system space conversion potential and the weight of the water system space conversion potential, the value of the linear space utilization potential and the weight of the linear space utilization potential, the value of the walking preference path and the weight of the walking preference path, the value of the cycling preference path and the weight of the cycling preference path, and the value of the resource attraction and the weight of the resource attraction are summed to determine the cost grid surface of the target recreational space. Among them, the cost factors can also be weighted and superimposed based on ArcGIS software to form a cost grid map, such as Figure 5 The cost raster surface shown.

[0101] As an optional implementation, the construction cost of a single grid can also be determined, and then the cost grid surface can be determined by cumulative addition. In this way, the formula for determining the cost grid surface of the target recreation space is:

[0102]

[0103] Where, E represents the construction cost index of the i-th evaluation grid; P i,j represents the value of the jth cost factor of the ith unit, W j Represents the weight of each cost factor.

[0104] Step 106, using the minimum path cost method, based on the landscape resistance surface of the target ecological source and the cost grid surface of the target recreational space, determine the ecological corridor for the target birds and the recreational corridor for humans.

[0105] In the implementation, the landscape resistance surface of the target ecological source and the cost grid surface of the target recreational space were selected, and the ecological corridors of potential target birds and the recreational corridors of humans were simulated based on the minimum cost path method. Figure 6 The ecological corridors and Figure 7 The recreational corridor shown.

[0106] Specifically, the specific steps of executing step 106 are as follows:

[0107] Step A, using the minimum path cost method, determine the minimum cost path between the landscape resistance surfaces of any two ecological nodes in the target ecological source area, and determine the minimum cost path as the ecological corridor of the target bird.

[0108] In the implementation, the minimum path cost method is used to determine the minimum cost path between the landscape resistance surfaces of any two ecological nodes in the target ecological source, and the minimum cost path is determined as the ecological corridor of the target bird. After selecting the target ecological source and landscape resistance surface, the Linkage PathwaysTool tool in the Linkage Mapper 3.1.0 plug-in can also be used to simulate the minimum cost path. For each ecological node, the ecological node is regarded as a circuit node, and then the CWD (Cost Weighted Distance) from all pixels in the landscape resistance surface to the ecological node is calculated to form a cost-weighted distance surface. The minimum cost weighted distance calculated is used to determine the minimum cost path from one ecological node to another. Among them, the maximum number of adjacent connections for each ecological node is limited to 4, and the cutoff cost weighted distance threshold is set to 200,000. It can be set according to actual conditions and is not limited here.

[0109] Step B, using the minimum path cost method, determines the minimum cost path between the cost grid surfaces of any two recreation nodes in each target recreation space, and determines the minimum cost path as the recreation corridor for humans.

[0110] In the implementation, the minimum path cost method is used to determine the minimum cost path between the cost grid surfaces of any two recreation nodes in each target recreation space, and the minimum cost path is determined as the recreational corridor for humans. After selecting the target recreation space and the cost grid surface, the Linkage PathwaysTool tool in the Linkage Mapper 3.1.0 plug-in can also be used to simulate the minimum cost path. For each recreation node, the recreation node is regarded as a circuit node, and then the CWD from all pixels in the cost grid surface to the recreation node is calculated to form a cost-weighted distance surface. The minimum cost weighted distance calculated is used to determine the minimum cost path from one recreation node to another. Among them, the maximum number of adjacent connections for each ecological node is limited to 4, and the cutoff cost weighted distance threshold is set to 200,000. It can be set according to the actual situation and is not restricted here.

[0111] Furthermore, ecological nodes and recreational nodes can be classified into protection / development categories, and surface source areas can be compounded and classified according to different compound types. The specific steps are as follows:

[0112] Step C, obtaining the bird diversity index of each ecological node in the target ecological source area.

[0113] During implementation, the bird diversity index of each ecological node in the target ecological source area is obtained, and the ecological nodes can be graded later according to the bird diversity index.

[0114] Step D: Based on the bird diversity index, each ecological node is graded to obtain multiple levels of ecological nodes of each ecological node.

[0115] In the implementation, based on the bird diversity index, each ecological node is graded to obtain multiple levels of ecological nodes of each ecological node. Among them, the multiple levels of ecological nodes can be the first-level protected ecological source area (A1), the second-level protected ecological source area (A2) and the third-level protected ecological source area (A3), such as Figure 8 shown.

[0116] Step E: Based on the preset evaluation system and random forest classifier, the recreation nodes of the target recreation space are classified to obtain recreation nodes of multiple levels.

[0117] In the implementation, based on the preset evaluation system, the recreation nodes of the target recreation space are classified using random forest to obtain multiple levels of recreation nodes. The multiple levels of recreation nodes can be "excellent recreation nodes" (B1), "good recreation nodes" (B2), "general recreation nodes (B3)", such as Fig. 9 As shown. The evaluation indicators of the evaluation system include node cognitive perception, user emotional attitude, node direct popularity, node potential popularity, node fun, cost-effectiveness and aesthetics. Among them, TF-IDF is term frequency-inverse document frequency, which is often used in the field of information retrieval and data mining to evaluate the importance of words in a certain document in the corpus. The collected comments on the recreational nodes are used for statistics, and the top 20 keywords with TF-IDF values ​​are selected as core subject words. Then the frequency of the core subject words appearing in each node is calculated as the node cognitive perception score of this embodiment. TF-IDF i，j =TF I，J ×IDF i，j, where TF-IDF means word frequency-inverse document frequency, TF means word frequency, and IDF means inverse document frequency. Secondly, the TF-IDF value of each recreational node is sorted and calculated, and the adjectives in the top 100 keywords are screened, and the words without emotional color are eliminated. Then the screened words are classified into positive, neutral and negative emotional attitudes, and the emotional values ​​are calculated. The specific rule is to score 1 point for positive emotional words, 0 points for neutral emotional words, and -1 points for negative emotional words. The total score is calculated as the emotional attitude score of the user in this embodiment. Then, the node direct popularity selects the indicator that can directly reflect the popularity of recreational resources, that is, the number of node comments on popular travel websites such as Ctrip and Dianping as the direct popularity score of the node. The potential popularity of the node is the degree of attention paid to the recreational resources, that is, the number of people who may visit in the future. The browsing popularity of the node on the Weibo website is obtained as the potential popularity score of the node. Finally, the node interest / cost-effectiveness / aesthetics respectively reflect the degree to which the recreational resources are interesting and playable, economically reasonable, and beautiful in landscape. The Houyi collector is used to collect the average values ​​of the interest / cost-effectiveness / scenery scores of the target recreational nodes on Ctrip.com as the node interest / cost-effectiveness / aesthetics scores.

[0118] Furthermore, by spatially overlaying the nodes according to different composite types, we can obtain nine composite types and usage instructions as shown in Table 1, which are used to clarify the priority of composite source function use. Fig.10 As shown, the first-level composite source includes the first-level ecological source-first-level recreation node, and the first-level ecological source-second-level recreation node; the second-level composite source includes the first-level ecological source-third-level recreation node, the second-level ecological source-first-level recreation node, and the second-level ecological source-second-level recreation node; the third-level composite source includes the second-level ecological source-third-level recreation node, the third-level ecological source-first-level recreation node, the third-level ecological source-second-level recreation node, and the third-level ecological source-third-level recreation node.

[0119] Table 1

[0120]

[0121] Furthermore, according to the potential corridor simulation paths of ecological corridors and recreational corridors, the composite type of single-function urban ecological corridors can be set to obtain the spatial location distribution of the route selection of composite-function urban ecological corridors. The composite-function urban ecological corridor path is based on the ecological corridor path and can meet the main function of species protection. The specific steps are as follows:

[0122] Step F, spatially overlaying and analyzing multiple levels of ecological nodes and multiple levels of recreational nodes to obtain multiple composite nodes and functional use levels of multiple composite nodes.

[0123] During implementation, multiple levels of ecological nodes and multiple levels of recreational nodes are spatially overlaid and analyzed to obtain multiple composite nodes and functional usage levels of multiple composite nodes.

[0124] Step G, perform spatial overlay analysis on ecological corridors and recreational corridors to extract urban ecological corridors with composite functions, which are based on ecological corridors. Among them, urban ecological corridors with composite functions are ecological corridors and recreational corridors connecting the same two nodes, ecological corridors and recreational corridors adjacent to spatial nodes, and ecological corridors and recreational corridors running through the same residential area.

[0125] In the implementation, the spatial superposition analysis of ecological corridors and recreational corridors was carried out to extract urban ecological corridors with complex functions. Fig.11 As shown, Figure a is a composite-functional urban ecological corridor. The composite-functional urban ecological corridor is defined as the linear area where the ecological corridor and the recreational corridor overlap or are adjacent. Among them, ArcGIS software can be used to perform spatial overlay analysis on the potential paths of the ecological / recreational corridors to extract the potential paths of the composite-functional urban ecological corridor. The composite-functional urban ecological corridor should focus on species protection and have slow-moving service functions. Combined with the actual urban construction goals and the feasibility of corridor construction, the composite-functional urban ecological corridor is based on the ecological corridor. Among them, the composite-functional urban ecological corridors are ecological corridors and recreational corridors connecting the same two nodes, ecological corridors and recreational corridors adjacent to spatial nodes, and ecological corridors and recreational corridors that run through the same residential area. As shown Fig.11 As shown in Figures BD, there are three complex types of urban ecological corridors with complex functions.

[0126] Furthermore, based on circuit theory, the ecological nodes and landscape resistance surfaces of the target ecological source areas are input into the Pinchpoint Mapper and Barrier Mapper modules respectively, and the pinch points / barrier spaces in the urban ecological corridors with complex functions and the pinch points / barrier spaces in the ecological corridors are extracted respectively. The pinch points / barrier spaces in the urban ecological corridors are used as the primary point protection / restoration spaces, and the pinch points / barrier spaces in the ecological corridors are used as the secondary point protection / restoration spaces.

[0127] In the implementation, based on the circuit theory, the ecological nodes and landscape resistance surfaces of the target ecological source are input into the Pinchpoint Mapper and Barrier Mapper modules respectively, and the pinch points / barrier spaces in the urban ecological corridor with complex functions and the pinch points / barrier spaces in the ecological corridor are extracted respectively. Among them, based on the Pointpoint Mapper tool, the "Pairwise" pair mode function is used to obtain the spatial distribution map of the ecological pinch points. Considering the connectivity of the overall landscape, the weighted cost distance can be set to 2000m. Using the Barrier Mapper tool, the percentage improvement score of the minimum cost path is used as the identification method of the barrier space. The core parameters of the module Minimum Detection Radius can be set to 150, Maximum Detection Radius can be set to 300, the detection step can be set to Wie50, and the combination method is Sum. Since the ecological pinch point is an area with high current density in the urban ecological corridor, it usually means that the possibility of passage in the area is high, or there is no other alternative path available. At the same time, these areas are very likely to become weak links that cut off the connection between ecological sources, and need to be given priority attention and protection. Ecological barrier points are areas where organisms are hindered from moving between sources. They are determined by calculating the size of the current recovery after removing the barrier points. Removing this area can enhance the connectivity between sources. Therefore, the pinch point space and barrier space are divided into primary pinch point / barrier space and secondary pinch point / barrier space, respectively. Among them, the pinch point space / barrier space in the ecological corridor of the complex functional city is regarded as the primary pinch point space / barrier space, and the pinch point / barrier space in the ecological corridor of the ecological city is regarded as the secondary pinch point space / barrier space. Fig.12 and Fig.13 The schematic diagram of the distribution of pinch point space and the schematic diagram of the distribution of obstacle space are shown. In this way, the spatial location of ecological pinch points and obstacle points can be clearly identified to provide guidance for the construction of key protection / restoration areas of the corridor.

[0128] Fig.14 A schematic diagram of a method for constructing a high-density urban ecological-recreational function composite corridor network provided in an embodiment of the present application. Fig.14As shown in the figure, ecological corridors can be constructed based on species protection to clearly indicate species (target birds). Then, correlation analysis was performed on phenological data and environmental factors to screen out environmental factors such as night light index, distance from roads, vegetation canopy height, distance from water systems, building proportion and normalized vegetation index. The MaxEnt model was used to extract ecological sources and construct comprehensive resistance surfaces (landscape resistance surfaces). The LCP minimum cost path model was used to simulate and extract potential ecological corridors to obtain ecological corridors. At the same time, recreational corridors were constructed based on slow travel, and leisure humans were selected as the target objects. The API interface of Gaode Map was used to obtain recreational resources, and recreational node screening was completed based on recreational resources. SPCA was used to perform principal component analysis on the component factors of recreational nodes to construct cost grid surfaces. Among them, the component factors include node service scope, spatial construction potential, residents' preferred path and resource attractiveness. The LCP minimum cost path model was used to simulate and extract potential recreational corridors to obtain recreational corridors. Furthermore, based on the surface ecological base for protection / development of ecological corridors and recreational corridors, the ecological source classification and recreational node classification are spatially superimposed and analyzed to obtain nine composite types. Linear ecological corridors are obtained based on ecological corridors and recreational corridors. Linear ecological corridors can be composite corridors that connect the same pair of source areas with target source areas, and are adjacent to spatial areas and pass through the same residential area. Point repair can also be performed to obtain point repair space, and the pinch point area and obstacle area are obtained based on circuit theory. In this way, a composite functional urban ecological corridor is completely obtained. In this way, the ecological corridor of the ecological city and the ecological corridor of the recreational city are spatially superimposed and analyzed, and the types are compounded from the perspectives of the surface ecological base, linear ecological corridor, and point repair space to complete the construction of a composite functional urban ecological corridor network. The construction of composite functional urban ecological corridors through the present invention is of great significance for alleviating land use conflicts, protecting urban species, and improving the quality of life of residents.

[0129] The embodiment of the present application provides a method for constructing an urban ecological corridor, and the target service object is selected according to the land use characteristics of high-density cities. Bird point data and urban green spaces with recreational functions are obtained as ecological nodes and recreational nodes respectively; environmental factors and cost factors are selected according to biological living habits and urban construction potential to construct landscape resistance surfaces and cost grid surfaces respectively. Potential ecological corridors and recreational corridors are extracted based on the minimum cost path model; the ecological corridors and recreational corridors are spatially superimposed and analyzed, and the types are compounded from the perspectives of planar ecological base, linear ecological corridors, and point-shaped restoration space to complete the construction of a composite-functional urban ecological corridor network. The construction of composite-functional urban ecological corridors through the present invention is of great significance for alleviating land use conflicts, protecting urban species, and improving the quality of life of residents. The present invention adopts a dual perspective of integrating ecological and recreational types. Through this method, it is intended to develop a composite-functional urban ecological corridor construction method that takes into account both biodiversity protection and residents' recreational needs, and to build an urban ecological corridor where people and nature live in harmony, which has important practical significance.

[0130] It should be understood that although Figure 1 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 1 At least part of the steps may include multiple steps or multiple stages. These steps or stages are not necessarily performed at the same time, but can be performed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed in turn or alternately with other steps or at least part of the steps or stages in other steps.

[0131] It can be understood that the same / similar parts between the various embodiments of the above method in this specification can refer to each other, and each embodiment focuses on the differences from other embodiments. For related points, please refer to the description of other method embodiments.

[0132] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

[0133] It should also be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for display, data for analysis, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties.

[0134] Each embodiment in this specification is described in a related manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0135] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0136] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims

1. A method for constructing a high-density urban ecological-recreational function composite corridor network, characterized in that: The method comprises: Inputting a plurality of first phenological data and a plurality of first environmental factors into a pre-trained maximum entropy model, and outputting a grid of a source of a suitable habitat for target birds; The grid of the source of the suitable habitat is reclassified by using the natural breakpoint method to obtain the candidate ecological source suitable for the target birds, and the candidate ecological source with an area greater than a preset threshold is determined as the target ecological source suitable for the target birds; The jackknife method is used to determine a plurality of second environmental factors corresponding to the target birds and greater than a preset correlation coefficient threshold and a contribution rate of the plurality of second environmental factors, and based on the plurality of second environmental factors and the corresponding contribution rates, the landscape resistance surface of the target ecological source is determined; Call the map software API application program interface to capture preset recreation keywords and screen out target recreation spaces that meet the preset recreation function requirements; The spatial principal component analysis method is used to construct the cost grid surface of the target recreation space based on the cost factors of the target recreation space; the cost factors include the node service range, spatial construction potential, residents' preferred path and resource attractiveness; The minimum path cost method is adopted to determine the ecological corridor of the target birds and the recreational corridor of humans based on the landscape resistance surface of the target ecological source and the cost grid surface of the target recreational space.

2. The method according to claim 1, characterized in that: Before inputting the plurality of first phenological data and the plurality of first environmental factors into the pre-trained maximum entropy model and outputting the grid of the source of the suitable habitat of the target bird species, the method further comprises: A correlation analysis method is used to determine the correlation coefficient values ​​between multiple phenological data and multiple environmental factors, and multiple first phenological data and multiple first environmental factors that are less than the preset correlation coefficient threshold are obtained.

3. The method according to claim 1, characterized in that Determining the landscape resistance surface of the target ecological source area according to the plurality of the second environmental factors and the corresponding contribution rates includes: The knife-cut method is used to determine the second environmental factor and its corresponding contribution rate, and multiple second environmental factors whose contribution rates are greater than a preset contribution rate threshold are determined as multiple key factors constituting the landscape resistance surface; A linear function obtained by adding the products of the plurality of key factors and the corresponding contribution rates is determined as the landscape resistance surface of the target ecological source area.

4. The method according to claim 1, characterized in that: The spatial construction potential includes road space conversion potential, water system space conversion potential and linear space utilization potential; the residents' preferred paths include walking preference paths and cycling preference paths; the spatial principal component analysis method is used to construct the cost grid surface of the target recreation space based on the cost factor of the target recreation space, including: According to a plurality of preset distance service thresholds, assigning a node service range of each recreation node of the target recreation space, and determining a weight of the node service range based on the spatial principal component analysis method; Assigning values ​​to the road space conversion potential of the target recreational space according to multiple road grades, and determining weights of the road space conversion potential based on the spatial principal component analysis method; According to the distance threshold between each pixel and the water system space in each recreation node, the water system space conversion potential of the target recreation space is assigned a value, and based on the spatial principal component analysis method, the weight of the water system space conversion potential is determined; Based on the color satellite image, the strip green space and the remaining other land in the target recreational space are identified, and based on the strip green space and the remaining other land, the linear space utilization potential of the target recreational space is assigned a value, and based on the spatial principal component analysis method, the weight of the linear space utilization potential is determined; Assigning values ​​to the walking preference paths of the target recreational space according to the walking frequency of the residents' walking preference paths, and determining weights of the walking preference paths based on the spatial principal component analysis method; Assigning values ​​to the cycling preference paths of the target recreational space according to the cycling frequency of the residents' cycling preference paths, and determining weights of the cycling preference paths based on the spatial principal component analysis method; Assigning a value to the resource attraction of the target recreational space according to a predetermined level of resource attraction, and determining a weight of the resource attraction based on the spatial principal component analysis method; The sum of the product of the value assigned to the node service range of the target recreational space and the weight of the node service range, the product of the value assigned to the road space conversion potential and the weight of the road space conversion potential, the product of the value assigned to the water system space conversion potential and the weight of the water system space conversion potential, the product of the value assigned to the linear space utilization potential and the weight of the linear space utilization potential, the product of the value assigned to the walking preference path and the weight of the walking preference path, the product of the value assigned to the cycling preference path and the weight of the cycling preference path, and the product of the value assigned to the resource attraction and the weight of the resource attraction is determined as the cost grid surface of the target recreational space.

5. The method according to claim 4, characterized in that The formula for determining the resource attractiveness of the target recreational space is: Where f(x) represents the resource attraction at the spatial location x, r represents the path distance attenuation threshold, n represents the number of facilities whose path distance to location x is less than or equal to r, and d ix represents the shortest distance from point i to facility point x, and k represents the spatial weight function.

6. The method according to claim 1, characterized in that The method of using the minimum path cost to determine the ecological corridor of the target birds and the recreational corridor of humans based on the landscape resistance surface of the target ecological source and the cost grid surface of the target recreational space includes: The minimum path cost method is used to determine the minimum cost path between the landscape resistance surfaces of any two ecological nodes in the target ecological source area, and the minimum cost path is determined as the ecological corridor of the target bird; The minimum path cost method is adopted to determine the minimum cost path between the cost grid surfaces of any two recreation nodes in each target recreation space, and the minimum cost path is determined as the recreational corridor for the human beings.

7. The method according to claim 1, characterized in that The method further comprises: Obtaining the bird diversity index of each ecological node in the target ecological source area; Based on the bird diversity index, each of the ecological nodes is graded to obtain multiple levels of ecological nodes of each of the ecological nodes; Based on the preset evaluation system and random forest classifier, the recreation nodes of the target recreation space are classified to obtain recreation nodes of multiple levels.

8. The method according to claim 7, characterized in that The method further comprises: Perform spatial overlay analysis on multiple ecological nodes of the above mentioned level and multiple recreational nodes of the above mentioned level to obtain multiple composite nodes and functional use levels of multiple composite nodes; A spatial overlay analysis is performed on the ecological corridors and the recreational corridors to extract urban ecological corridors with complex functions. The urban ecological corridors with complex functions are based on the ecological corridors. The urban ecological corridors with complex functions are ecological corridors and recreational corridors connecting the same two nodes, ecological corridors and recreational corridors with adjacent spatial nodes, and ecological corridors and recreational corridors that run through the same residential area.

9. The method according to claim 8, characterized in that The method further comprises: Based on circuit theory, the ecological nodes of the target ecological source and the landscape resistance surface are input into the Pinchpoint Mapper and Barrier Mapper modules respectively, and the pinch points / obstacle spaces in the complex-functional urban ecological corridor and the pinch points / obstacle spaces in the ecological corridor are extracted respectively. The pinch points / obstacle spaces in the urban ecological corridor are used as the primary point protection / restoration space, and the pinch points / obstacle spaces in the ecological corridor are used as the secondary point protection / restoration space.

10. The method according to claim 1, characterized in that The suitable growth zone sources include high suitable growth zone sources, medium suitable growth zone sources, low suitable growth zone sources and non-suitable growth zone sources.

Images