Evaluation method for suitability of land space development in mountainous areas

Abstract

The present invention discloses a method for evaluating the suitability of land and space development in mountainous areas, including: data acquisition and preprocessing: collecting and organizing basic data of the target area; calculating the importance of ecological protection: evaluating the importance of ecological protection from the perspectives of ecosystem service function importance and ecological vulnerability, and obtaining an ecological protection importance evaluation map; analyzing and evaluating the suitability of agricultural production and urban construction: constructing a preliminary dual-suitability scheme for agricultural production suitability zoning and urban construction suitability zoning; refining the dual-suitability: establishing a multi-threshold spatial morphological optimization algorithm to identify centrifugal distance, land use compactness, and area thresholds, adjusting the final boundary according to a function-dominant strategy, correcting fragmented small faces of mountainous land parcels, modifying narrow patches, and improving the degree of contiguousness; and outputting the evaluation results of the suitability of land and space development in mountainous areas. By optimizing the algorithm and function-dominant strategy, the present invention solves the problem of distinguishing "suitability" in complex geomorphic areas and improves the accuracy and practicality of the evaluation results.

Description

Technical Field

[0001] This invention relates to the field of land and space planning technology, and in particular to a method for evaluating the suitability of land and space development in mountainous areas. Background Technology

[0002] Existing methods for assessing the suitability of land development are mostly applicable to areas with relatively flat terrain. When applied to mountainous areas with complex topography, they may have the following shortcomings: (1) Single-function-oriented assessment methods are difficult to reconcile the contradiction between the multi-functional suitability of land and the uniqueness of suitability assessment results, resulting in low accuracy in the division of "dual-suitable areas"; (2) Traditional suitability assessment methods do not consider the characteristics of fragmented land parcels, vertical mountain terrain, and high incidence of disasters in mountainous areas, resulting in a large deviation between the assessment results and actual application needs; (3) In the process of balancing the dual goals of agricultural production and urban construction, it is difficult to accurately distinguish "dual-suitable areas", which can easily lead to land use conflicts and affect the rational development and efficient utilization of land. Existing methods have failed to effectively classify functional suitability when applied in mountainous areas and have not solved the problem of low contiguousness of land parcels. Therefore, it is urgent to construct a land space assessment model for mountainous areas with complex topography to solve the problem of difficulty in quantifying and classifying the suitability of land development in mountainous areas, solve the "dual suitability" problem, and improve the accuracy of suitability differentiation and computational efficiency. Summary of the Invention

[0003] To address the problems existing in the prior art, the purpose of this invention is to provide a method for evaluating the suitability of land development in mountainous areas. This invention solves the problem of distinguishing "suitability" in complex geomorphological areas by optimizing algorithms and functional-driven strategies, thereby improving the accuracy and practicality of the evaluation results.

[0004] To achieve the above objectives, the technical solution adopted by this invention is: a method for evaluating the suitability of land development in mountainous areas, comprising the following steps:

[0005] Step 1, Data Acquisition and Preprocessing: Collect and organize basic data for the target area;

[0006] Step 2: Calculate the importance of ecological protection: Evaluate the importance of ecological protection from two aspects: the importance of ecosystem service functions and ecological vulnerability, and obtain an ecological protection importance evaluation map;

[0007] Step 3: Analyze and evaluate the suitability of agricultural production and urban construction: Construct a preliminary plan for dual suitability based on agricultural production suitability zones and urban construction suitability zones;

[0008] Step 4, Dual Suitability Refinement: Establish a multi-threshold spatial morphology optimization algorithm to identify centrifugal distance, land use compactness and area threshold, adjust the final boundary according to the function-driven strategy, correct the fragmented small facets of mountainous plots, modify narrow patches, and improve the degree of contiguousness;

[0009] Step 5: Output the results of the suitability evaluation for land development in mountainous areas.

[0010] As a further improvement of the present invention, in step 1, the basic data includes: basic geography, land resources, water resources, disasters, and climate and meteorological data.

[0011] As a further improvement of the present invention, in step 2, the importance of the ecosystem service functions includes: the importance of water conservation function, the importance of soil and water conservation function, the importance of biodiversity maintenance function, and the importance of windbreak and sand fixation function; the ecological vulnerability includes vulnerability to soil erosion, vulnerability to desertification, and vulnerability to rocky desertification.

[0012] As a further improvement of the present invention, step 2 is specifically as follows:

[0013] Based on the principle of the resource and environmental carrying capacity short-plate effect, the evaluation results of the ecological importance of the target area are obtained by superimposing and analyzing the evaluation results of single factors:

[0014] ;

[0015] In the formula, Indicates the level of importance of ecological protection. These respectively indicate the importance levels of water conservation, soil and water conservation, biodiversity maintenance, and windbreak and sand fixation functions. These represent the vulnerability levels of soil erosion, desertification, and rocky desertification, respectively.

[0016] As a further improvement of the present invention, step 3 specifically includes the following steps:

[0017] Step 3.1: Conduct individual evaluations of agricultural production conditions from the aspects of land resources, water resources, light and heat conditions, environmental pollution, and climate disasters; overlay and discriminate the evaluation results of land resources, water resources, and light and heat conditions to obtain preliminary results of agricultural production condition levels; establish an agricultural production condition correction matrix, correct the preliminary evaluation results based on the evaluation results of environmental pollution and climate disasters, determine the key influencing factors of agricultural production, and obtain the results of agricultural production suitability zoning.

[0018] Step 3.2: The suitability assessment for urban construction is conducted comprehensively from the aspects of land resources, water resources, climate comfort, environmental capacity, geological hazards, and locational advantages. The assessment results of land resources, water resources, and climate comfort are superimposed and discriminated to obtain preliminary results of the urban construction condition level. A dual coupling factor of geological hazards and topographic relief is introduced to establish an urban construction condition correction matrix. The preliminary results are corrected based on the assessment results of environmental capacity, geological hazard risk, and locational advantages to determine the key influencing factors for urban construction and obtain the urban construction suitability zoning results.

[0019] Step 3.3: Use spatial masking technology to remove the distribution areas of extremely important level in the ecological protection importance assessment map, and construct preliminary plans for agricultural production suitability zoning and urban construction suitability zoning.

[0020] As a further improvement of the present invention, step 4 specifically includes the following steps:

[0021] Step 4.1: Calculate the plot area, land use compactness, and centrifugal distance;

[0022] Step 4.2: Based on the mountainous terrain data, use the sliding window analysis method to identify the centrifugal distance. and The patches, identifying the area of ​​aggregated land parcels and The patches; among which:

[0023] ;

[0024] In the formula, E is the total length of the patch boundary, and TA is the patch area;

[0025] Step 4.3: Correlation of Topographic Relief. Based on the centrifugal distance, the suitability level is optimized, and the preliminary plans for agricultural production suitability zones and urban construction suitability zones are downgraded. Here, S=1 indicates a one-level downgrade in suitability, and S=0 indicates no change in suitability level.

[0026] ;

[0027] In the formula, S represents the suitability level after centrifugal distance optimization;

[0028] Step 4.4: Based on suitability level S, according to the plot area and land compactness Further optimize the suitability level:

[0029] ;

[0030] In the formula, This indicates the suitability level after optimization of land parcel area and land use compactness. Indicates an unsuitable area. Indicates the suitable area. Indicates the aggregated land parcel area, Indicates land parcel Land use compactness, Indicates land parcel perimeter.

[0031] As a further improvement of the present invention, in step 4, the functional dominant strategy includes: (1) combining the priority protection of high-quality farmland to divide suitable areas for agricultural production; (2) combining historical and cultural protection lines, mineral resource distribution and urban planning to determine priority development areas for towns and divide suitable areas for urban construction.

[0032] As a further improvement of the present invention, step 5 is as follows: Generate an ecological protection importance distribution map, an agricultural production suitability distribution map, and an urban construction suitability distribution map in ArcGIS software; divide the dual-suitable zones according to the functional-driven strategy; output the final agricultural production suitability zoning map and urban construction suitability zoning map; and generate a comparative analysis statistical report of the dual-suitable zones before and after refined correction.

[0033] This invention addresses the challenge of accurately distinguishing the suitability of agricultural production and urban development in complex topographical areas. Driven by geographic information big data, it establishes a multi-threshold spatial morphology optimization algorithm (MTSMOA) and constructs a multi-dimensional evaluation method for the suitability of mountainous land development through a function-driven strategy. The method includes: ecological protection importance assessment, suitability classification for agricultural production and urban development, dynamic correction based on centrifugal distance and land use compactness, and refined division of "dual-suitable zones." This invention solves the difficulties in quantifying and classifying the suitability of mountainous land development and addresses the "dual-suitability" issue, significantly improving the accuracy of suitability evaluation and computational efficiency. It provides scientific decision-making support for land spatial planning and achieves synergistic optimization of resource development and ecological protection.

[0034] The beneficial effects of this invention are:

[0035] When conducting suitability assessments for land development in mountainous areas, this invention improves the accuracy of delineating "dual-suitable zones" in mountainous areas to over 70%, reducing the area after correction by 42.8%; the accuracy of identifying important ecological protection zones is improved, accounting for 31.68% of the total area; and the method for assessing the suitability of land development in mountainous areas yields assessment results for the importance of ecological protection, suitability for agricultural production, and suitability for urban construction, supporting the scientific delineation of "three zones and three lines," reducing land use conflicts, and coordinating the contradictions between resource development and ecological protection. Attached Figure Description

[0036] Figure 1 This is an overall architecture diagram of an embodiment of the present invention. Detailed Implementation

[0037] The embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

[0038] Example

[0039] like Figure 1 As shown, a method for evaluating the suitability of land development in mountainous areas includes the following steps:

[0040] Step 1: Data Acquisition and Preprocessing. Collect and organize basic data for the target area, including basic geography, land resources, water resources, disasters, climate and meteorology data.

[0041] Step 2: Calculate the importance of ecological protection. The importance of ecological protection is evaluated from two aspects: the importance of ecosystem service functions and ecological vulnerability. Based on the principle of the "weakest link effect" of resource and environmental carrying capacity, the evaluation results of individual factors are overlaid to obtain the evaluation result of the importance of ecological protection.

[0042] Step 3: Analyze and evaluate the suitability of agricultural production and urban construction. Individual evaluations of land resources, climate, environment, disasters, and location are conducted sequentially. Discrimination and correction matrices are used for grade correction. An integrated evaluation is then performed. Spatial masking technology is used to remove areas of extremely important grade from the ecological protection importance evaluation map. Preliminary plans for agricultural production suitability zoning and urban construction suitability zoning are then constructed.

[0043] Step 4, “Dual Suitability” Refinement. Using the Multi-Threshold Spatial Morphology Optimization (MTSMOA) algorithm, land use compactness and area thresholds are used to adjust the final boundaries according to the function-driven strategy, correcting fragmented facets in mountainous areas and modifying narrow patches to improve contiguousness;

[0044] Step 5: Output the suitability assessment results for land development in mountainous areas. Generate an ecological protection importance zoning map based on the assessment results from Step 2. Using the revised results from Step 4, divide the area into "dual-suitable zones" according to the functional-driven strategy, and output the final agricultural production suitability zoning map and urban construction suitability zoning map. A comparative analysis report is obtained by comparing the results before and after the refined revision of the "dual-suitable zones" in Step 4.

[0045] In step 1, the basic data of the target area includes basic geography, land resources, water resources, environment, ecology, disasters, climate and meteorology data.

[0046] In step 2, the importance of the ecosystem service functions includes the importance of water conservation, soil and water conservation, biodiversity maintenance, and windbreak and sand fixation. Ecological vulnerability mainly includes vulnerability to soil erosion, desertification, and rocky desertification.

[0047] Based on the "weakest link effect" principle of resource carrying capacity, the evaluation results of each individual element are superimposed and analyzed to obtain the evaluation results of the ecological protection importance of the target area:

[0048] ;

[0049] In the formula, Indicates the level of importance of ecological protection. These respectively indicate the importance levels of water conservation, soil and water conservation, biodiversity maintenance, and windbreak and sand fixation functions. These represent the vulnerability levels of soil erosion, desertification, and rocky desertification, respectively.

[0050] Step 3 specifically includes the following steps:

[0051] Step 3.1: Conduct a single-item evaluation of agricultural production conditions from the aspects of land resources, water resources, light and heat conditions, environmental pollution, and climate disasters.

[0052] (1) Land resource evaluation: The land grade is classified according to the slope. When the soil silt content is between 60% and 80%, it has a significant impact on crops, while when it exceeds 80%, it is completely unsuitable for planting crops. Therefore, based on the slope grading results, if the soil silt content grade is poor, it is classified as poor; if the silt content grade is relatively poor, the slope grade is downgraded by 1 level.

[0053] (2) Water resources assessment: Based on long-term precipitation observation data, a multi-year average precipitation distribution layer is obtained through spatial interpolation. For areas where the proportion of transit water sources is large and it is difficult to fully reflect agricultural water supply conditions by relying solely on the total local water resources, the modulus of the total water consumption control index of county-level administrative regions is used for calculation.

[0054] (3) Evaluation of light and heat conditions: Based on long-term temperature observation data, spatial interpolation is performed and corrected by altitude to obtain the active accumulated temperature.

[0055] (4) Environmental pollution assessment: Based on the results of a detailed survey of soil pollution in the region, the content of major pollutants at each location is analyzed and spatial interpolation is performed.

[0056] (5) Climate Disaster Assessment: Collect and organize historical data on various climate elements and meteorological disasters, and statistically analyze the frequency of occurrence of individual meteorological disasters each year, including drought, flooding, high-temperature heat damage, low-temperature cold damage, and strong wind disasters, and classify their risks. The agricultural production suitability evaluation index system is shown in Table 1:

[0057] Table 1

[0058] ;

[0059] The evaluation results of land resources, water resources, and light and heat conditions are overlaid and discriminated to obtain preliminary results of the agricultural production condition level; the resource-based discrimination matrix is ​​shown in Table 2:

[0060] Table 2

[0061] ;

[0062] An agricultural production condition correction matrix was established. The preliminary evaluation results were corrected based on the environmental pollution and climate disaster assessment results to determine the key influencing factors for agricultural production. Factors rated as "good" or "relatively good" were classified as suitable, while those rated as "average," "relatively poor," or "poor" were classified as unsuitable. The agricultural production condition correction matrix is ​​shown in Table 3.

[0063] Table 3

[0064] ;

[0065] Step 3.2: The suitability assessment for urban construction is conducted by comprehensively evaluating aspects such as land resources, water resources, climate comfort, environmental capacity, geological hazards, and locational advantages.

[0066] (1) Land resource assessment: The land is divided according to the slope. The slope classification result is the basis. If the elevation grade is poor, the land resource grade is poor. If the elevation grade is relatively poor, the slope grade is reduced by 1 level. For areas with poor topographic relief, the assessment result is reduced by 2 levels. If the topographic relief is relatively poor, the assessment result is reduced by 1 level.

[0067] (2) Water resources assessment: Classification is based on the total water resources modulus of county-level administrative regions. For areas where the total local water resources are difficult to fully reflect the urban water supply conditions, the total water consumption control index modulus of county-level administrative regions is used for classification.

[0068] (3) Climate comfort evaluation: The temperature and humidity index (THI) is used to characterize climate comfort. Combined with the climate characteristics of Sichuan Province, the climate comfort is classified.

[0069] (4) Environmental capacity assessment: including atmospheric environmental capacity and water environmental capacity assessment. Spatial interpolation of calm days and average wind speed is used to assess atmospheric environmental capacity; water environmental capacity is calculated for major pollutants such as chemical oxygen demand and ammonia nitrogen.

[0070] (5) Geological hazard assessment: including seismic hazard, landslide, debris flow susceptibility, ground subsidence susceptibility and ground collapse susceptibility assessment. The highest level among active faults, landslides, debris flows, ground subsidence and ground collapse is taken as the geological hazard hazard level.

[0071] (6) Location Advantage Evaluation: An optimized location advantage algorithm was independently developed, which calculates travel time to the central city via "point-to-point" calculations, improving computational speed by 90% and accuracy by 30%. The urban construction suitability evaluation index system is shown in Table 4:

[0072] Table 4

[0073] ;

[0074] The results of the evaluations of land resources, water resources, and climate comfort are superimposed and judged to obtain preliminary results of the urban construction condition level.

[0075] A dual-coupling factor of geological hazards and topographic relief was introduced to establish a correction matrix for urban construction conditions. The preliminary results were corrected based on the evaluation results of environmental capacity, geological hazard risk, and locational advantage. Key influencing factors for urban construction were determined, and urban construction conditions rated as good or relatively good were classified as suitable, while those rated as average, relatively poor, or poor were classified as unsuitable. The correction matrix for urban construction conditions is shown in Table 5.

[0076] Table 5

[0077] ;

[0078] Step 3.3: Use spatial masking technology to remove the extremely important areas in the ecological protection importance assessment map, and construct preliminary plans for agricultural production suitability zones and urban construction suitability zones.

[0079] In step 4, the Multi-Threshold Spatial Morphology Optimization Algorithm (MTSMOA) specifically includes the following steps:

[0080] Step 4.1: Calculate the plot area, land compactness, and centrifugal distance.

[0081] Step 4.2: Based on the mountainous terrain data, use the sliding window analysis method to identify the centrifugal distance. and The patches, identifying the area of ​​aggregated land parcels and plaques.

[0082] Centrifugal distance This indicates the degree of deviation between the shape of the patch and a circle of the same area. (Using a circle as a reference)

[0083] ;

[0084] In the formula, E is the total length of the patch boundary, and TA is the patch area.

[0085] Step 4.3: Correlation of Topographic Relief Based on centrifugal distance to optimize suitability level, the preliminary plans for agricultural production suitability zones and urban construction suitability zones are downgraded, where S=1 indicates a one-level reduction in suitability and S=0 indicates no change in suitability level.

[0086] ;

[0087] In the formula, S represents the suitability level after centrifugal distance optimization.

[0088] Step 4.4: Based on suitability level S, according to the plot area and land compactness The suitability level was optimized again.

[0089] ;

[0090] In the formula, This indicates the suitability level after optimization of land parcel area and land use compactness. Indicates an unsuitable area. Indicates the suitable area. Indicates the aggregated land parcel Area, in km² 2 , Indicates land parcel Land use compactness, Indicates land parcel Perimeter, measured in km.

[0091] The functional-dominant strategies include: (1) dividing suitable agricultural production areas by prioritizing the protection of high-quality farmland; and (2) determining priority development areas for towns and dividing suitable urban construction areas by combining historical and cultural protection lines, mineral resource distribution and urban planning.

[0092] In step 5, based on the evaluation results obtained in steps 2 and 4, an ecological protection importance distribution map, an agricultural production suitability distribution map, and an urban construction suitability distribution map are generated in ArcGIS software, and a comparative statistical report on the "dual-suitable area" before and after refined correction is formed.

[0093] The following uses the method of this embodiment to conduct a suitability assessment for land development in mountainous areas:

[0094] Regarding the importance of ecological protection, the area of ​​the ecologically extremely important zone in a certain province is 154,000 km². 2 This accounts for 31.68% of the region's area. The area of ​​the ecological protection important zone is 214,000 km². 2 This accounts for 44.02% of the region's area. Regarding agricultural suitability, the area of ​​suitable agricultural production zones is 63,500 km². 2 This accounts for 13.07% of the total area. Areas unsuitable for agricultural production are widely distributed, covering 268,600 km². 2 This accounts for 55.25% of the total area. Regarding suitability for urban development, the area suitable for urban development is 82,000 km². 2 This accounts for 16.86% of the region's area. Meanwhile, areas unsuitable for urban development cover 250,100 km². 2 It accounts for 51.46% of the region's area.

[0095] The "dual suitability" issue was refined, and the area of ​​the "dual suitability zone" before the revision was 59,800 km². 2 This accounts for 12.30% of the total area. The revised area of ​​the "dual-suitable zone" has decreased to 34,200 km². 2 This accounts for 7.04% of the total area. Except for a certain plain economic zone, the accuracy of other economic zones all exceeded 70%. Among them, the area of ​​the dual-suitable zone in a certain economic zone ranged from 49.90 km². 2 Reduced to 3.30 km 2 The accuracy of the distinction reached 93.38%. The area of ​​suitable agricultural production zones increased by 62,800 km² after the correction. 2 The revised suitable urban construction area has decreased by 37,200 km² compared to the previous version. 2 The "dual-suitable zone" amendment effectively mitigated potential land use conflicts between agricultural production and urban development.

[0096] The embodiments described above are merely illustrative of specific implementations of the present invention, and while the descriptions are detailed, they should not be construed as limiting the scope of the present invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention.

Claims

1. A method for evaluating the suitability of land development in mountainous areas, characterized in that, Includes the following steps: Step 1, Data Acquisition and Preprocessing: Collect and organize basic data for the target area; Step 2: Calculate the importance of ecological protection: Evaluate the importance of ecological protection from two aspects: the importance of ecosystem service functions and ecological vulnerability, and obtain an ecological protection importance evaluation map; Step 3: Analyze and evaluate the suitability of agricultural production and urban construction: Construct a preliminary plan for dual suitability based on agricultural production suitability zones and urban construction suitability zones; Step 4, Dual Suitability Refinement: Establish a multi-threshold spatial morphology optimization algorithm to identify centrifugal distance, land use compactness and area threshold, adjust the final boundary according to the function-driven strategy, correct the fragmented small facets of mountainous plots, modify narrow patches, and improve the degree of contiguousness; Step 4 specifically includes the following steps: Step 4.1: Calculate the plot area, land use compactness, and centrifugal distance; Step 4.2: Based on the mountainous terrain data, use the sliding window analysis method to identify the centrifugal distance D. i >0.4 and D i For patches ≤0.4, identify aggregated plot area A. i ≤0.5 and A i Patches >0.5; among which: In the formula, E is the total length of the patch boundary, and TA is the patch area; Step 4.3: Correlation of Topographic Relief. Based on the centrifugal distance, the suitability level is optimized, and the preliminary plans for agricultural production suitability zones and urban construction suitability zones are downgraded. Here, S=1 indicates a one-level downgrade in suitability, and S=0 indicates no change in suitability level. In the formula, S represents the suitability level after centrifugal distance optimization; Step 4.4: Based on the suitability level S, and according to the plot area A i and land compactness C i Further optimize the suitability level: In the formula, MTSMOA i This indicates the suitability level after optimization of land parcel area and land use compactness, with S0 representing an unsuitable area. 1-R Indicates the suitable area, A i C represents the area of ​​plot i after aggregation. i F represents the land use compactness of plot i. i Indicates the perimeter of plot i; Step 5: Output the results of the suitability evaluation for land development in mountainous areas.

2. The method for evaluating the suitability of land development in mountainous areas according to claim 1, characterized in that, In step 1, the basic data includes: basic geography, land resources, water resources, disasters, and climate and meteorological data.

3. The method for evaluating the suitability of land development in mountainous areas according to claim 1, characterized in that, In step 2, the importance of the ecosystem service functions includes: the importance of water conservation function, the importance of soil and water conservation function, the importance of biodiversity maintenance function, and the importance of windbreak and sand fixation function; the ecological vulnerability includes vulnerability to soil erosion, vulnerability to desertification, and vulnerability to rocky desertification.

4. The method for evaluating the suitability of land development in mountainous areas according to claim 3, characterized in that, Step 2 is described in detail below: Based on the principle of the resource and environmental carrying capacity short-plate effect, the evaluation results of the ecological importance of the target area are obtained by superimposing and analyzing the evaluation results of single factors: IE i =MAX{TQ i ,TA i ,TB i ,TG i ,TS i ,TD i ,TR i } In the formula, EI i Indicates the level of importance of ecological protection, TQ i TA i TB i TG i These respectively represent the importance levels of water conservation, soil and water conservation, biodiversity maintenance, and windbreak and sand fixation functions. TS i TD i TR i These represent the vulnerability levels of soil erosion, desertification, and rocky desertification, respectively.

5. The method for evaluating the suitability of land development in mountainous areas according to claim 1, characterized in that, Step 3 specifically includes the following steps: Step 3.1: Conduct individual evaluations of agricultural production conditions from the aspects of land resources, water resources, light and heat conditions, environmental pollution, and climate disasters; overlay and discriminate the evaluation results of land resources, water resources, and light and heat conditions to obtain preliminary results of agricultural production condition levels; establish an agricultural production condition correction matrix, correct the preliminary evaluation results based on the evaluation results of environmental pollution and climate disasters, determine the key influencing factors of agricultural production, and obtain the results of agricultural production suitability zoning. Step 3.2: The suitability assessment for urban construction is conducted comprehensively from the aspects of land resources, water resources, climate comfort, environmental capacity, geological hazards, and locational advantages. The assessment results of land resources, water resources, and climate comfort are superimposed and discriminated to obtain preliminary results of the urban construction condition level. A dual coupling factor of geological hazards and topographic relief is introduced to establish an urban construction condition correction matrix. The preliminary results are corrected based on the assessment results of environmental capacity, geological hazard risk, and locational advantages to determine the key influencing factors for urban construction and obtain the urban construction suitability zoning results. Step 3.3: Use spatial masking technology to remove the distribution areas of extremely important level in the ecological protection importance assessment map, and construct preliminary plans for agricultural production suitability zoning and urban construction suitability zoning.

6. The method for evaluating the suitability of land development in mountainous areas according to claim 1, characterized in that, In step 4, the functional-dominant strategy includes: (1) dividing suitable agricultural production areas by prioritizing the protection of high-quality farmland; and (2) determining priority development areas for towns and dividing suitable urban construction areas by combining historical and cultural protection lines, mineral resource distribution and urban planning.

7. The method for evaluating the suitability of land development in mountainous areas according to claim 6, characterized in that, Step 5 is as follows: Generate an ecological protection importance distribution map, an agricultural production suitability distribution map, and an urban construction suitability distribution map in ArcGIS software. Divide the dual-suitability zones according to the functional-driven strategy, output the final agricultural production suitability zoning map and urban construction suitability zoning map, and generate a comparative analysis statistical report of the dual-suitability zones before and after refined correction.

Images