Method for evaluating ecological suitability of agricultural land

Abstract

The present invention relates to the technical field of ecological suitability evaluation of agricultural land, and disclosed is a method for evaluating the ecological suitability of agricultural land. The present invention comprises: determining values of indicators of agricultural land, wherein the values of the indicators of the agricultural land comprise slope, elevation, slope length, profile curvature, soil thickness, plan curvature, the minimum number of days of continuous flooding, the maximum number of days of continuous flooding, the cumulative number of days of flooding, the average number of days of continuous flooding, minimum flooding depth, maximum flooding depth, cumulative flooding depth, and average flooding depth. In the present invention, the natural ecology of agricultural land, i.e., whether regional water, soil and climate conditions meet the ecological requirements of the agricultural land, is evaluated, so that it is ensured that the agricultural land can guarantee the sustainability of an ecosystem and environment; the ecological status thereof can be quickly and accurately identified, weights of indicators in an element layer can be determined, and finally, a system for evaluating the ecological suitability of agricultural land is constructed; the method can reflect the ecological suitability of agricultural land in a timely, macroscopic, and dynamic manner.

Description

A method for evaluating the ecological suitability of agricultural land Technical Field

[0001] The present invention belongs to the field of agricultural land ecological suitability evaluation, and in particular relates to an agricultural land ecological suitability evaluation method. Background Art

[0002] Ecological suitability refers to the degree to which, within a given habitat, the space provided by climate, ecology, topography, hydrology, and other factors, as well as the material cycles and energy flows, are suitable for the survival and positive evolution of biological communities. It can also be understood as the maximum intensity of human development permitted by a habitat while ensuring the structural and functional integrity of the ecosystem. Existing research on the ecological suitability of agricultural land has rarely considered indicators such as climate and land conditions, resulting in inaccurate assessments of agricultural land ecological suitability.

[0003] In view of this, the present invention is proposed.

[0004] Summary of the Invention

[0005] The technical problem to be solved by the present invention is to overcome the deficiencies of the prior art and provide a method for evaluating the ecological suitability of agricultural land.

[0006] In order to solve the above technical problems, the basic concept of the technical solution adopted by the present invention is:

[0007] A method for evaluating the ecological suitability of agricultural land comprises the following steps:

[0008] Step 1, determine the index values ​​of agricultural land, including slope, elevation, slope length, profile curvature, soil thickness, plane curvature, shortest continuous flooding days, longest continuous flooding days, cumulative flooding days, average continuous flooding days, minimum flooding depth, maximum flooding depth, cumulative flooding depth, and average flooding depth;

[0009] Step 2: Determine the weight of each indicator and calculate the value of each factor based on the value and weight of each indicator. The factors include topographic factors and hydrological factors. Topographic factor indicators include slope, elevation, slope length, profile curvature, soil thickness, and plane curvature. Hydrological factor indicators include the shortest continuous flooding days, the longest continuous flooding days, the cumulative flooding days, the average continuous flooding days, the minimum flooding depth, the maximum flooding depth, the cumulative flooding depth, and the average flooding depth.

[0010] Step 3: Determine the weight of each factor and calculate the comprehensive index of ecological suitability evaluation based on the value and weight of each factor;

[0011] Step 4: For different agricultural lands, select the required indicators from the indicator system, use the information entropy-corrected hierarchical analysis method and the Delphi method to determine the weight value of each indicator, and construct a regional agricultural land ecological suitability measurement indicator evaluation system based on the characteristics of agricultural land. Determine the assignment rules for different indicators, collect data on regional related indicators, and spatialize and rasterize them based on geographic information system technology (GIS) technology to construct a raster dataset with consistent coordinate system and spatial range. According to the evaluation criteria of different indicators, normalize and assign values ​​respectively;

[0012] Step 5: By performing weighted summation of the indicator weights and the assigned values ​​pixel by pixel, a comprehensive index map of the ecological suitability evaluation of different agricultural lands is generated. On this basis, the spatial data reclassification technology and the comprehensive index of ecological suitability evaluation are used to divide the comprehensive index map of the ecological suitability evaluation of agricultural land into five levels: very suitable, relatively suitable, basically suitable, unsuitable, and very unsuitable.

[0013] Optionally, step 4, and using the entropy weight method to calculate the weight of the agricultural land represented by each evaluation indicator includes: using the maximum value normalization method to process the evaluation matrix composed of the selected multiple evaluation indicators to obtain a standard matrix; according to the obtained standard matrix, calculating the information entropy value of the evaluation indicator.

[0014] Optionally, the slope calculation formula in step 1 is:

[0015] S is the slope; fx is the rate of change of elevation in the x-direction; fy is the rate of change of elevation in the y-direction;

[0016] The slope length calculation formula is:

[0017] L is the slope length; DEM is the DEM vector data of the water-fluctuation zone; S is the slope;

[0018] The calculation formulas for profile curvature and plane curvature are:

[0019] Where Cp is the profile curvature; Cc is the plane curvature; Z = f(x1, y1) is the function defined in the DEM data coordinate system of the drawdown zone, fx1 is the first-order partial derivative of x1, fy1 is the first-order partial derivative of y1, fx1x1 is the second-order partial derivative of x1, fy1y1 is the second-order partial derivative of y1, and fx1y1 is the first-order partial derivative of x1 and then the first-order partial derivative of y1.

[0020] Optionally, the method for determining the values ​​of the longest continuous flooding days, the shortest continuous flooding days, the average continuous flooding days, the cumulative flooding days, the maximum flooding depth, the minimum flooding depth, the average flooding depth, and the cumulative flooding depth in step 1 is as follows: daily water level data from the agricultural land hydrological station, calculate the longest continuous flooding days in the drawdown zone (the longest number of continuous flooding days in a year), the shortest continuous flooding days (the shortest number of continuous flooding days in a year), the average continuous flooding days, the cumulative flooding days (the total number of days flooded in a year), the maximum flooding depth (the maximum depth of flooding in a year), the minimum flooding depth (the minimum depth of flooding in a year), the average flooding depth, and the cumulative flooding depth (the total flooding depth in a year) hydrological factors.

[0021] Optionally, the formula for calculating the average duration of flooding in step 1 is: Where Dav is the average number of days of continuous flooding; dt is the total number of days of continuous flooding in a year; f is the number of continuous flooding in a year;

[0022] The formula for calculating the average flooding depth is: Where Haf is the average flooding depth; hf is the sum of the flooding depths for each day of the year; and df is the number of flooding days in a year.

[0023] Optionally, step 3 also includes testing the effectiveness of each indicator. The testing steps are as follows: the indicators are screened using the principal component analysis method, and the principal components with a cumulative variance contribution rate greater than 95% are selected as model construction indicators. Due to the different magnitudes of the indicators, the range method is used to standardize the data. The formula is as follows: When the indicator is positive, that is, the larger the indicator data, the better the ecological suitability, including elevation, slope length, soil thickness, etc., the following formula is used for standardization:

[0024] When the indicator is negative, that is, the larger the indicator data, the worse the ecological suitability, including slope, profile curvature, plane curvature, longest continuous flooding days, shortest continuous flooding days, average continuous flooding days, cumulative flooding days, maximum flooding depth, minimum flooding depth, average flooding depth, cumulative flooding depth, etc., the following formula is used for standardization:

[0025] Where Pij is the standardized result of each indicator, Xij is the actual value, and Xmax and Xmin represent the maximum and minimum values, respectively. Testing the effectiveness of indicators can reduce the correlation between model indicators and avoid reducing the accuracy of model results due to duplication of indicators.

[0026] Optionally, the indicator weights in step 2 are determined by spatial principal components, and the formula is as follows:

[0027] Where Hj is the common factor variance of the indicator; Wj is the indicator weight; λ is the transformed main factor eigenvector; j is the number of original indices; k is the number of components; and m is the total number of main components.

[0028] Optionally, the calculation of the comprehensive index of ecological suitability evaluation in step 3 adopts a multi-objective integration method to calculate the comprehensive index. The calculation formulas of the criterion layer and the target layer are:

[0029] Where I is the comprehensive index of ecological suitability evaluation; Ai is the weight of the evaluation index; and Ui is the standardized value of the evaluation index.

[0030] Optionally, the indicators in step 4 covering natural ecological agricultural land can simultaneously meet the goals of economic development and ecological protection. For different agricultural lands, based on expert opinions, the required indicators are selected and the indicator levels are divided.

[0031] Optionally, in step 4, the weight of each indicator is determined, and it is decomposed into multiple levels from top to bottom according to different attributes, and the indicators in the same level are subordinate to the indicators of the upper level. The design of the present invention includes four levels, the highest level is the target level, which is the "evaluation standard of the ecological suitability measurement index of a certain agricultural land" here, the second level is the sub-target level, which is the "natural ecology" and "social economy" goals respectively, the third level is the criterion level, and the fourth level is the indicator level. Starting from the second level, a pair comparison matrix is ​​constructed for the same-level indicators that are subordinate to the same upper-level indicators.

[0032] After adopting the above technical solution, the present invention has the following beneficial effects compared with the prior art. Of course, any product implementing the present invention does not necessarily need to achieve all of the advantages described below at the same time:

[0033] By analyzing the natural ecology of agricultural land, that is, whether the regional water, soil, and climate conditions meet the ecological needs of agricultural land, it is ensured that agricultural land can guarantee the sustainability of the ecosystem and environment, and its ecological status can be identified quickly and accurately to determine whether it is suitable for vegetation development. Topographic factors and hydrological elements are used as element layers, and the weights of each element layer are clarified. Parameters such as elevation, slope, length, soil layer thickness, longest continuous flooding days, shortest continuous flooding days, average continuous flooding days, cumulative flooding days, and maximum flooding depth are used as indicator layers to determine the weights of each indicator in the element layer. Finally, an agricultural land ecological suitability evaluation system is constructed. This method can reflect the ecological suitability of agricultural land in a timely, macroscopic and dynamic manner.

[0034] The specific embodiments of the present invention will be described in further detail below with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] The drawings described below are only some embodiments. A person skilled in the art can derive other drawings based on these drawings without inventive effort. In the drawings:

[0036] FIG1 is a flow chart of a method according to an embodiment of the present invention.

[0037] It should be noted that these drawings and textual descriptions are not intended to limit the conceptual scope of the present invention in any way, but rather to illustrate the concept of the present invention for those skilled in the art by referring to specific embodiments. DETAILED DESCRIPTION

[0038] The present invention will now be described in further detail with reference to the accompanying drawings.

[0039] Referring to FIG1 , this embodiment provides a method for evaluating the ecological suitability of agricultural land, including the following steps:

[0040] Step 1, determine the index values ​​of agricultural land, including slope, elevation, slope length, profile curvature, soil thickness, plane curvature, shortest continuous flooding days, longest continuous flooding days, cumulative flooding days, average continuous flooding days, minimum flooding depth, maximum flooding depth, cumulative flooding depth, and average flooding depth;

[0041] Step 2: Determine the weight of each indicator and calculate the value of each factor based on the value and weight of each indicator. The factors include topographic factors and hydrological factors. The indicators of topographic factors include slope, elevation, slope length, profile curvature, soil thickness, and plane curvature. The indicators of hydrological factors include the shortest continuous flooding days, the longest continuous flooding days, the cumulative flooding days, the average continuous flooding days, the minimum flooding depth, the maximum flooding depth, the cumulative flooding depth, and the average flooding depth.

[0042] Step 3: Determine the weight of each factor and calculate the comprehensive index of ecological suitability evaluation based on the value and weight of each factor;

[0043] Step 4: For different agricultural lands, select the required indicators from the indicator system, use the information entropy-corrected hierarchical analysis method and the Delphi method to determine the weight value of each indicator, and construct a regional agricultural land ecological suitability measurement indicator evaluation system based on the characteristics of agricultural land. Determine the assignment rules for different indicators, collect data on regional related indicators, and spatialize and rasterize them based on geographic information system technology (GIS) technology to construct a raster dataset with consistent coordinate system and spatial range. According to the evaluation criteria of different indicators, normalize and assign values ​​respectively;

[0044] Step 5: By performing weighted summation of the indicator weights and the assigned values ​​pixel by pixel, a comprehensive index map of the ecological suitability evaluation of different agricultural lands is generated. On this basis, the spatial data reclassification technology and the comprehensive index of ecological suitability evaluation are used to divide the comprehensive index map of the ecological suitability evaluation of agricultural land into five levels: very suitable, relatively suitable, basically suitable, unsuitable, and very unsuitable.

[0045] Step 4 of this embodiment, and using the entropy weight method to calculate the weight of the agricultural land represented by each evaluation indicator, includes: using the maximum value normalization method to process the evaluation matrix composed of the selected multiple evaluation indicators to obtain a standard matrix; and calculating the information entropy value of the evaluation indicator based on the obtained standard matrix.

[0046] The slope calculation formula in step 1 of this embodiment is:

[0047] S is the slope; fx is the rate of change of elevation in the x direction; fy is the rate of change of elevation in the y direction;

[0048] The slope length calculation formula is:

[0049] L is the slope length; DEM is the DEM vector data of the water-fluctuation zone; S is the slope;

[0050] The calculation formulas for profile curvature and plane curvature are:

[0051] Where Cp is the profile curvature; Cc is the plane curvature; Z = f(x1, y1) is the function defined in the DEM data coordinate system of the drawdown zone, fx1 is the first-order partial derivative of x1, fy1 is the first-order partial derivative of y1, fx1x1 is the second-order partial derivative of x1, fy1y1 is the second-order partial derivative of y1, and fx1y1 is the first-order partial derivative of x1 and then the first-order partial derivative of y1.

[0052] The method for determining the values ​​of the longest continuous flooding days, the shortest continuous flooding days, the average continuous flooding days, the cumulative flooding days, the maximum flooding depth, the minimum flooding depth, the average flooding depth, and the cumulative flooding depth in step 1 of this embodiment is as follows: the daily water level data of the agricultural land hydrological station is used to calculate the hydrological factors of the longest continuous flooding days (the longest number of continuous flooding days in a year), the shortest continuous flooding days (the shortest number of continuous flooding days in a year), the average continuous flooding days, the cumulative flooding days (the total number of days flooded in a year), the maximum flooding depth (the maximum depth of flooding in a year), the minimum flooding depth (the minimum depth of flooding in a year), the average flooding depth, and the cumulative flooding depth (the total flooding depth in a year) in the drawdown zone.

[0053] The calculation formula for the average continuous flooding days in step 1 of this embodiment is:

[0054] Where Dav is the average number of days of continuous flooding; dt is the total number of days of continuous flooding in a year; f is the number of continuous flooding in a year;

[0055] The formula for calculating the average flooding depth is: Where Haf is the average flooding depth; hf is the sum of the flooding depths for each day of the year; and df is the number of flooding days in a year.

[0056] Step 3 of this embodiment also includes testing the effectiveness of each indicator. The testing steps are as follows: the indicators are screened using the principal component analysis method, and the principal components with a cumulative variance contribution rate greater than 95% are selected as model construction indicators. Due to the different magnitudes of the indicators, the range method is used to standardize the data. The formula is as follows: When the indicator is positive, that is, the larger the indicator data, the better the ecological suitability, including elevation, slope length, soil thickness, etc., the following formula is used for standardization:

[0057] When the indicator is negative, that is, the larger the indicator data, the worse the ecological suitability, including slope, profile curvature, plane curvature, longest continuous flooding days, shortest continuous flooding days, average continuous flooding days, cumulative flooding days, maximum flooding depth, minimum flooding depth, average flooding depth, cumulative flooding depth, etc., the following formula is used for standardization:

[0058] Where Pij is the standardized result of each indicator, Xij is the actual value, and Xmax and Xmin represent the maximum and minimum values, respectively. Testing the effectiveness of indicators can reduce the correlation between model indicators and avoid reducing the accuracy of model results due to duplication of indicators.

[0059] In step 2 of this embodiment, the index weight is determined by using the spatial principal component, and the formula is as follows:

[0060] Where Hj is the common factor variance of the indicator; Wj is the indicator weight; λ is the transformed main factor eigenvector; j is the number of original indices; k is the number of components; and m is the total number of main components.

[0061] In step 3 of this embodiment, the calculation of the comprehensive index of ecological suitability evaluation adopts a multi-objective integration method to calculate the comprehensive index. The calculation formulas of the criterion layer and the target layer are:

[0062] Where I is the comprehensive index of ecological suitability evaluation; Ai is the weight of the evaluation index; and Ui is the standardized value of the evaluation index.

[0063] In step 4 of this embodiment, the indicators covering natural ecological agricultural land can simultaneously meet the goals of economic development and ecological protection. For different agricultural lands, the required indicators are selected based on expert opinions and divided into indicator levels.

[0064] In step 4 of this embodiment, the weight of each indicator is determined, and it is decomposed into multiple levels from top to bottom according to different attributes, and the indicators in the same level are subordinate to the indicators of the upper level. The design of the present invention includes four levels, the highest level is the target level, which is the "evaluation standard of the ecological suitability measurement index of a certain agricultural land" here, the second level is the sub-target level, which includes the "natural ecology" and "socio-economic" goals respectively, the third level is the criterion level, and the fourth level is the indicator level. Starting from the second level, a pair comparison matrix is ​​constructed for the indicators of the same level that are subordinate to the same upper level indicators.

[0065] The present invention ensures that agricultural land can guarantee the sustainability of the ecosystem and environment by analyzing the natural ecology of agricultural land, that is, whether the regional water, soil, and climate conditions meet the ecological needs of agricultural land. It can quickly and accurately identify its ecological status and determine whether it is suitable for vegetation development. It uses terrain factors and hydrological elements as element layers, clarifies the weight of each element layer, and uses parameters such as elevation, slope, length, soil layer thickness, longest continuous flooding days, shortest continuous flooding days, average continuous flooding days, cumulative flooding days, and maximum flooding depth as indicator layers to determine the weight of each indicator in the element layer, and finally constructs an agricultural land ecological suitability evaluation system. This method can reflect the ecological suitability of agricultural land in a timely, macroscopic and dynamic manner.

[0066] The contents not described in detail in this specification belong to the prior art known to those skilled in the art.

[0067] The present invention is not limited to the above-described embodiments. Any structural changes made under the guidance of the present invention, which have the same or similar technical solutions as the present invention, should be understood to fall within the scope of protection of the present invention. The technologies, shapes, and structural parts not described in detail in the present invention are all well-known technologies.

Claims

1. A method for evaluating the ecological suitability of agricultural land, characterized in that: The steps include: Step 1, determine the index values ​​of agricultural land, including slope, elevation, slope length, profile curvature, soil thickness, plane curvature, shortest continuous flooding days, longest continuous flooding days, cumulative flooding days, average continuous flooding days, minimum flooding depth, maximum flooding depth, cumulative flooding depth, and average flooding depth; Step 2, determine the weight of each indicator, and calculate the value of each element according to the value and weight of each indicator. The elements include terrain factors and hydrological elements. The indicators of terrain factors include slope, elevation, slope length, profile curvature, soil thickness, and plane curvature. The indicators of hydrological elements include the shortest continuous flooding days, the longest continuous flooding days, the cumulative flooding days, the average continuous flooding days, the minimum flooding depth, the maximum flooding depth, the cumulative flooding depth, and the average flooding depth. Step 3, determine the weight of each factor, and calculate the comprehensive index of ecological suitability evaluation according to the value and weight of each factor; Step 4: For different agricultural lands, select the required indicators in the indicator system, use the information entropy modified hierarchical analysis method and Delphi method to determine the weight value of each indicator, combine the characteristics of agricultural land, construct the regional agricultural land ecological suitability measurement indicator evaluation system, determine the assignment rules of different indicators, collect data on regional related indicators, and perform spatialization and rasterization based on geographic information system technology (GIS) to construct a raster data set with consistent matching of coordinate system and spatial range, and normalize and assign values ​​according to the evaluation standards of different indicators; Step 5: By performing weighted summation of the indicator weights and the assigned values ​​pixel by pixel, a comprehensive index map of ecological suitability evaluation of different agricultural lands is generated. On this basis, the spatial data reclassification technology and the comprehensive index of ecological suitability evaluation are used to divide the comprehensive index map of ecological suitability evaluation of agricultural land into five levels: very suitable, relatively suitable, basically suitable, unsuitable, and very unsuitable.

2. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: Step 4, using the entropy weight method to calculate the weight of the agricultural land represented by each evaluation index includes: using the maximum value standardization method to process the evaluation matrix composed of the selected multiple evaluation indicators to obtain a standard matrix; according to the obtained standard matrix, calculating the information entropy value of the evaluation index.

3. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: The slope calculation formula in step 1 is: S is the slope; fx is the rate of change of elevation in the x direction; fy is the rate of change of elevation in the y direction; The slope length calculation formula is: L is the slope length; DEM is the DEM vector data of the water-falling zone; S is the slope; The calculation formulas for profile curvature and plane curvature are: In the formula, Cp is the profile curvature; Cc is the plane curvature; Z=f(x1, y1) is the function defined in the coordinate system of the drawdown zone DEM data, fx1 is the first-order partial derivative of x1, fy1 is the first-order partial derivative of y1, fx1x1 is the second-order partial derivative of x1, fy1y1 is the second-order partial derivative of y1, and fx1y1 is the first-order partial derivative of x1 and then the first-order partial derivative of y1.

4. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: The method for determining the values ​​of the longest continuous flooding days, the shortest continuous flooding days, the average continuous flooding days, the cumulative flooding days, the maximum flooding depth, the minimum flooding depth, the average flooding depth, and the cumulative flooding depth in step 1 is as follows: the daily water level data of the agricultural land hydrological station is used to calculate the longest continuous flooding days in the drawdown zone (the longest number of continuous flooding days in a year), the shortest continuous flooding days (the shortest number of continuous flooding days in a year), the average continuous flooding days, the cumulative flooding days (the total number of days flooded in a year), the maximum flooding depth (the maximum depth of flooding in a year), the minimum flooding depth (the minimum depth of flooding in a year), the average flooding depth, and the cumulative flooding depth (the total number of flooding depths in a year) hydrological factors.

5. The method for evaluating the ecological suitability of agricultural land according to claim 4, characterized in that: The calculation formula for the average duration of flooding in step 1 is: In the formula, Dav is the average number of days of continuous flooding; dt is the total number of days of continuous flooding in a year; f is the number of continuous flooding in a year; The average flood depth calculation formula is: In the formula, Haf is the average hf is the average flooding depth; hf is the sum of the flooding depths for each day of the year; df is the number of flooding days in a year.

6. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: Step 3 also includes testing the effectiveness of each indicator. The testing steps are as follows: the indicators are screened using the principal component analysis method, and the principal components with a cumulative variance contribution rate greater than 95% are selected as model construction indicators. Due to the different magnitudes of the indicators, the range method is used to standardize the data. The formula is as follows: When the indicator is a positive indicator, that is, the larger the indicator data, the better the ecological suitability, including elevation, slope length, soil thickness, etc., the following formula is used for standardization: When the index is negative, that is, the larger the index data, the worse the ecological suitability, including slope, profile curvature, plane curvature, longest continuous flooding days, shortest continuous flooding days, average continuous flooding days, cumulative flooding days, maximum flooding depth, minimum flooding depth, average flooding depth, cumulative flooding depth, etc., the following formula is used for standardization: In the formula, Pij is the standardized result of each indicator, Xij is the actual value, and Xmax and Xmin represent the maximum and minimum values, respectively. Performing indicator validity testing can reduce the correlation between model indicators and avoid the reduction of the accuracy of model results due to indicator duplication.

7. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: In step 2, the indicator weights are determined by spatial principal components, and the formula is as follows: In the formula, Hj is the common factor variance of the indicator; Wj is the indicator weight; λ is the transformed main factor eigenvector; j is the number of original indexes; k is the number of components; and m is the total number of main components.

8. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: The calculation of the comprehensive index of ecological suitability evaluation in step 3 adopts a multi-objective integration method to calculate the comprehensive index. The calculation formulas of the criterion layer and the target layer are: In the formula, I is the comprehensive index of ecological suitability evaluation; Ai is the weight of the evaluation index, and Ui is the standardized value of the evaluation index.

9. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: The indicators in step 4 cover natural ecological agricultural land and can meet the goals of economic development and ecological protection at the same time. For different agricultural lands, expert opinions are integrated to select the required indicators and divide the indicator levels.

10. The method for evaluating the ecological suitability of agricultural land according to claim 1, characterized in that: In step 4, the weight of each indicator is determined, and it is decomposed into multiple levels from top to bottom according to different attributes, and the indicators in the same level are subordinate to the indicators in the upper level. The design of the present invention includes four levels, the highest level is the target level, which is "evaluation standard of ecological suitability measurement index of a certain agricultural land" here, the second level is the sub-target level, which is "natural ecology" and "social economy" goals respectively, the third level is the criterion level, and the fourth level is the indicator level. Starting from the second level, the indicators of the same level that are subordinate to the same upper level indicators are constructed into a pair comparison matrix.

Images