A method and system for evaluating arable land reserves

Abstract

The present application relates to the technical field of cultivated land protection examination and evaluation, in particular to a cultivated land holding amount examination and evaluation method and system. The examination and evaluation method comprises the following steps: S1: obtaining a large database based on cultivated land protection; S2: screening out disaster destroyed cultivated land, cultivated land converted from non-cultivated land and three-region-three-line outflow cultivated land data in the large database respectively, and calculating disaster reduction number, cultivated land conversion reduction number and three-region-three-line outflow reduction number respectively; S3: standardizing the cultivated land protection data, calculating the degree of human controllability according to the standardized cultivated land protection data, and then calculating the disaster weight, cultivated land conversion weight and three-region-three-line outflow weight respectively, and obtaining the cultivated land protection examination and evaluation score through weighted calculation. The present application improves the data processing efficiency by splitting the data, makes the data more clear, improves the verification accuracy, realizes efficient and accurate quality evaluation of the operation results of multiple databases, and improves the utilization rate of data processing.

Description

Technical Field

[0001] This invention relates to the field of farmland protection assessment and evaluation technology, and in particular to a method and system for assessing and evaluating farmland reserves. Background Technology

[0002] In recent years, the state has conducted comprehensive assessments of the implementation of the responsibility system for farmland protection and food security in many provinces. With the completion of the third national land survey and the new round of territorial spatial planning, as well as the official implementation of the "three zones and three lines" delineation results, more precise new requirements have been put forward for the responsibility system for farmland protection and food security.

[0003] When assessing the farmland protection situation in various provinces, the state issued a document that scores five categories of indicators: farmland area; protection of permanent basic farmland; balance between farmland occupation and replenishment, balance between farmland entry and exit, and farmland restoration tasks; illegal occupation of farmland and permanent basic farmland; and bonus items. However, the document only provides a simple explanation of the indicators without establishing specific scoring methods and procedures. The rules are unclear and the procedures are ambiguous, making it vague for each province on how to operate the scoring indicators of the farmland protection responsibility system. Summary of the Invention

[0004] Therefore, it is necessary to propose a method and system for assessing and evaluating arable land area, addressing the problems of unclear rules and ambiguous indicators in the existing assessment and evaluation of arable land protection responsibilities.

[0005] This invention is achieved through the following technical solution: A method for assessing and evaluating arable land area includes the following steps:

[0006] S1: Obtain a large database based on farmland protection.

[0007] S2: Filter out data on disaster-damaged farmland, farmland converted from farmland, and farmland flowing out of the "Three Zones and Three Lines" in the large database, and calculate the reduction figures for disaster-damaged farmland, farmland converted from farmland, and farmland flowing out of the "Three Zones and Three Lines" respectively.

[0008] S3: The farmland protection assessment score is calculated using a weighted calculation method. The specific method is as follows:

[0009] S31: Standardize the processing of farmland protection data.

[0010] S32: Calculate the degree of human controllability based on standardized farmland protection data.

[0011] S33: Calculate the disaster damage weight, farmland conversion weight, and outflow weight for the three zones and three lines based on the degree of human control.

[0012] S34: The reduction figures for disaster damage, reduction figures for land conversion from farmland to forest, and reduction figures for land outflow from the three regions and three lines are weighted and calculated to obtain the farmland protection assessment score.

[0013] The aforementioned assessment method comprehensively considers farmland protection efforts resulting from natural disasters, ecological restoration, and the outflow of farmland due to the "Three Zones and Three Lines" policy. By intersecting and erasing multiple databases, it obtains data on farmland protected by natural disasters, ecological restoration, and the outflow of farmland due to the "Three Zones and Three Lines" policy. Compared to directly statistically analyzing relevant farmland data, this embodiment demonstrates more interconnected database processing, improves data utilization, and achieves higher statistical accuracy, thereby enhancing the objectivity of the assessment. By splitting the data, not only is data processing efficiency improved, but the data also becomes clearer, less prone to confusion, and increases verification accuracy. Furthermore, with the assistance of data algorithms, it is possible to accelerate the processing of large amounts of data, employing more optimized algorithms and models to process the data, thus achieving efficient and accurate quality assessment of the results from multiple databases.

[0014] Furthermore, in step S2, the method for calculating the adjustment reduction number is as follows:

[0015] S21: Preprocess the large database to obtain multiple land use patches with consistent proportions, namely, patches of the same caliber, urban and rural patches (without open mapping), non-protected arable land patches, ecological conversion of arable land patches, urban and rural land use patches, arable land protection patches for the "Three Zones and Three Lines" (referring to the three zones and three lines of protected arable land use), and permanent arable land patches for the "Three Zones and Three Lines." Perform an intersection operation on the patches of the same caliber and the urban and rural patches (without open mapping) to obtain the basic land use patches.

[0016] S22: Filter disaster-damaged farmland data, calculate the standard deviation of disaster-damaged farmland and the average value of disaster-damaged protected farmland, and calculate the disaster reduction amount based on the set coefficients. Then, the disaster reduction amount P... ZH Expressed as:

[0017]

[0018] In the formula, k1 and k2 are manually set coefficients, and σ(S ZH ) represents the standard deviation of the area of ​​farmland destroyed by disaster, and Avg() represents the average value. S is the target control area for disaster recovery. ZHTK The actual controlled area for disaster recovery.

[0019] S23: Filter the data on converted farmland, calculate the standard deviation of converted farmland and the average value of the controlled conversion of farmland for protection, and calculate the reduction in converted farmland based on the set coefficients. Then, the reduction in converted farmland, P... TG Expressed as:

[0020]

[0021] In the formula, σ(S) TG ) represents the standard deviation of the area of ​​farmland converted from cultivated land. S is the target area for ecological restoration through farmland conversion. ZHTKThe actual controlled area for ecological restoration of farmland.

[0022] S24: Filter the cultivated land data for the "Three Zones and Three Lines" (referring to the three zones and three lines of cultivated land outflow), calculate the standard deviation of cultivated land outflow from the "Three Zones and Three Lines" and the average value of the controlled cultivated land outflow from the "Three Zones and Three Lines" protected cultivated land, and calculate the reduction amount of outflow from the "Three Zones and Three Lines" based on the set coefficients. Then, the reduction amount P of outflow from the "Three Zones and Three Lines" is... GTK Expressed as:

[0023]

[0024] In the formula, σ(S) SSL The standard deviation of the area of ​​cultivated land flowing out of the three regions and three lines is denoted as . The target control area for the outflow of arable land in the three zones and three lines, S SSLTK This refers to the actual area of ​​cultivated land flowing out of the "three zones and three lines" (referring to specific administrative regions and economic zones).

[0025] Furthermore, in step S22, the method for filtering disaster-damaged farmland data is as follows:

[0026] S221: Select disaster-damaged patches from patches of the same caliber, and select unused farmland flow patches from basic land use patches. Perform an intersection operation on the disaster-damaged patches and the unused farmland flow patches to obtain disaster-damaged unused patches.

[0027] S222: Erase unprotected cultivated land parcels from the disaster-damaged unused parcels to obtain the disaster-damaged parcels. Read the fields of the disaster-damaged parcels and calculate the disaster-damaged cultivated land area. Then the disaster-damaged cultivated land area S... ZH This can be expressed as:

[0028] S ZH =S'[(LW∩ZH)-FBH]*β -2

[0029] LW=CH LW (DL), ZH=CH ZH (TKJ)

[0030] In the formula, S'[] is the area of ​​the map patch, LW is the farmland flow to unused map patch, ZH is the disaster-damaged map patch, FBH is the unprotected farmland map patch, β is the map patch scale, DL is the basic land type map patch, and CH() is the filtering operation.

[0031] Furthermore, in step S23, the method for filtering the data on converted farmland is as follows:

[0032] S231: Select farmland-to-forest-grass-wetland patches from the land use patches. Perform an intersection operation on the farmland-to-forest-grass-wetland patches and the ecological conversion of farmland to forest-grass-wetland patches to obtain the forest-grass-wetland conversion of farmland patches.

[0033] S232: Erase unprotected cultivated land patches in the forest-grassland wetland conversion map to obtain the converted cultivated land patches. Read the fields of the converted cultivated land patches and calculate the converted cultivated land area. The converted cultivated land area S is then calculated. TG Expressed as:

[0034] S TG =S'[(TG∩LCS)-FBH]*β -2

[0035] LCS=CH LCS (DL)

[0036] In the formula, TG represents the ecological conversion of farmland to forest and grassland wetland plots, and LCS represents the conversion of farmland to forest and grassland wetland plots.

[0037] Furthermore, in step S24, the calculation method for filtering the cultivated land data of the three zones and three lines is as follows:

[0038] 241: Erase urban and rural map patches from the basic land use map patches to obtain erased urban and rural map patches. Select non-arable urban and rural map patches from the erased urban and rural map patches.

[0039] 242: Erase the permanent cultivated land patches of the Three Zones and Three Lines in the cultivated land protection patch to obtain the general cultivated land patches of the Three Zones and Three Lines.

[0040] 243: Perform an intersection operation on non-cultivated land parcels in urban and rural areas and general cultivated land parcels along the "Three Zones and Three Lines" to obtain cultivated land parcels flowing out of the "Three Zones and Three Lines". Read the fields of the cultivated land parcels flowing out of the "Three Zones and Three Lines" and calculate the area of ​​cultivated land flowing out of the "Three Zones and Three Lines". Then the area of ​​cultivated land flowing out of the "Three Zones and Three Lines" is S. SSL Expressed as:

[0041] S SSL =S'[FG∩(SSB-SSY)]*β -2

[0042] FG = CH FG (DL-CZC)

[0043] In the formula, FG represents non-arable land parcels in urban and rural areas, SSB represents protected arable land parcels in the three zones and three lines, SSY represents permanent arable land parcels in the three zones and three lines, and CZC represents arable land parcels used in urban and rural areas.

[0044] Furthermore, in step S32, the formula for calculating the degree of human controllability is:

[0045]

[0046] In the formula, R ZH The `avg()` function calculates the average value of farmland protected from disasters, representing the degree of human control over its preservation. To protect the target area of ​​disaster-damaged farmland under a unified standard, SDZH To standardize the actual protected arable land area of ​​disaster-damaged farmland, R TG The degree of human control over the conversion of farmland to protected farmland. To protect the target area of ​​cultivated land under a unified standard for returning farmland to forest, SD TG To determine the actual protected arable land area under a unified standard for returning farmland to forest, R SSL The degree of human control over the outflow of arable land in the "three zones and three lines" (referring to areas with high population density and rural development zones). To protect the target area of ​​cultivated land outflowing under unified standards, SD SSL The actual protected arable land area under a unified standard for arable land outflow.

[0047] Furthermore, in step S33, the disaster weight ω ZK Weight of land conversion ω TG Three-zone, three-line outflow weight ω SSL The calculation formulas are as follows:

[0048] ω ZH =R ZH / (R ZH +R TG +R SSL )

[0049] ω TG =R TG / (R ZH +R TG +R SSL )

[0050] ω SSL =R SSL / (R ZH +R TG +R SSL ).

[0051] Furthermore, in step S34, the farmland protection assessment score P S The calculation formula is:

[0052] P S =ω ZH P ZH +ω TG P TG +ω SSL P SSL

[0053] In the formula, P ZH P TG P SSL These are the reduction figures for disaster damage, reduction figures for land conversion from farmland to forest, and reduction figures for outflow from the "three regions and three lines" (referring to areas with high population density and rural development).

[0054] This invention also provides a system for assessing and evaluating arable land area, comprising a data acquisition module, a preprocessing module, a map patch processing module, and a calculation module.

[0055] The data acquisition module is used to collect data from a large database based on farmland protection.

[0056] The preprocessing module is used for: 1. Converting the large database into multiple land use patches with consistent proportions. 2. Standardizing farmland protection data.

[0057] The patch processing module is used to perform intersection, erasure, and filtering operations on different patches in order to retain the required patch data.

[0058] The calculation module is used for: 1. Calculating the reduction amount of disaster-damaged farmland based on the selected disaster-damaged farmland data. 2. Calculating the reduction amount of farmland converted from cultivation based on the selected farmland conversion data. 3. Calculating the reduction amount of farmland outflow from the "Three Zones and Three Lines" based on the selected farmland data from the "Three Zones and Three Lines" (referring to the three zones and three lines of farmland protection). 4. Calculating the degree of human controllability based on the standardized farmland protection data. 5. Calculating the disaster damage weight, farmland conversion weight, and outflow weight from the "Three Zones and Three Lines" based on the degree of human controllability. 6. Performing a weighted calculation on the disaster damage reduction amount, farmland conversion reduction amount, and outflow weight from the "Three Zones and Three Lines" based on the disaster damage weight, farmland conversion weight, and outflow weight to obtain the farmland protection assessment score.

[0059] Compared with the prior art, the present invention has the following beneficial effects:

[0060] The evaluation method of this invention comprehensively considers farmland protection caused by natural disasters, ecological restoration, and farmland outflow due to the "Three Zones and Three Lines" policy. By performing intersection and erasure operations on multiple databases, it obtains data on farmland protected by natural disasters, ecological restoration, and farmland outflow due to the "Three Zones and Three Lines" policy. Compared to directly statistically analyzing relevant farmland data, this embodiment demonstrates more interconnected database processing, improves data utilization, and achieves higher statistical accuracy, thereby enhancing the objectivity of the evaluation score. By splitting the data, it not only improves processing efficiency but also makes the data clearer, less prone to data confusion, and increases verification accuracy. Furthermore, with the assistance of data algorithms, it can accelerate the processing of large amounts of data, employing more optimized algorithms and models to process the data, thus achieving efficient and accurate quality assessment of the results from multiple databases.

[0061] This invention, based on data from disaster-induced reductions, land conversion reductions, and outflow reductions from the "Three Zones and Three Lines" (referring to areas with concentrated agricultural land use and forestry), further employs a data aggregation and calculation method. Each reduction data point is assigned a weight and scored separately. Finally, a weighted calculation method is used to obtain the final assessment score for arable land preservation. This assessment method not only refines the scoring rules for different arable land protection categories and fully considers actual changes in arable land protection, but also sets the weights for each protection category based on the degree of human controllability. This allows for a more objective and accurate evaluation of arable land protection status and provides important guidance for setting assessment standards for arable land protection. Attached Figure Description

[0062] Figure 1 This is a flowchart illustrating the steps of the method for assessing and evaluating arable land area in Embodiment 1 of the present invention.

[0063] Figure 2 A flowchart for calculating disaster reduction figures;

[0064] Figure 3 A flowchart for calculating the reduction of farmland due to land conversion;

[0065] Figure 4 A flowchart for calculating the outflow reduction figures for the three zones and three lines;

[0066] Figure 5 A schematic diagram of the framework structure of the assessment and evaluation system for arable land retention. Detailed Implementation

[0067] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0068] It should be noted that when a component is said to be "installed on" another component, it can be directly on the other component or it may be in a component that is centered on it. When a component is said to be "set on" another component, it can be directly set on the other component or it may also be in a component that is centered on it. When a component is said to be "fixed to" another component, it can be directly fixed to the other component or it may also be in a component that is centered on it.

[0069] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. The terminology used herein in the specification of this invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The term "or / and" as used herein includes any and all combinations of one or more of the associated listed items.

[0070] Example 1: Please refer to Figure 1 This embodiment provides a method for assessing and evaluating arable land reserves, including the following steps:

[0071] S1: Obtain a large database based on farmland protection. The large database includes a database with the same scope, a database of urban and rural areas that are not open, a non-protection database, a database of ecological farmland conversion, a database of urban and rural land use, a database of farmland protection in the "Three Zones and Three Lines" area, and a database of permanent farmland in the "Three Zones and Three Lines" area.

[0072] Among them, the "consistent caliber database" refers to data collected within the same scope. Taking a specific city as an example, it statistically analyzes the land use status of all land within that city, with each statistical analysis maintaining a completely consistent scope and type. The "town and village database" refers to a database that only includes statistics on the current state of land use, without including specific town and village names. The "non-protection database" refers to a database that does not include protected arable land in land use statistics. The "ecological conversion of arable land" database includes databases of arable land converted to other uses as permitted by law. The "town and village land use database" is a database of arable land used in urban and village construction. In the "three zones and three lines" framework, the "three zones" refer to the three functional zones of agriculture, ecology, and urban areas, and the "three lines" refer to permanent basic farmland, ecological protection red lines, and urban development boundaries.

[0073] S2: Filter data on disaster-damaged farmland, farmland converted from farmland, and farmland flowing out of the "Three Zones and Three Lines" (referring to designated areas for agricultural development and rural revitalization) from the large database, and calculate the corresponding reduction amounts for each. The specific method is as follows:

[0074] S21: Preprocess the data to convert various types of data into land use patches with consistent scale. Specifically, convert the above database into patches of the same caliber (TKJ), urban and rural patches (BDK, unopened), non-protected cultivated land patches (FBH), ecological conversion of cultivated land patches (TG), urban and rural land use patches (CZC), cultivated land protection patches for the "Three Zones and Three Lines" (SSB), and permanent cultivated land patches for the "Three Zones and Three Lines" (SSY).

[0075] An intersection operation is performed on the same-caliber map patch TKJ and the unopened urban / village map patch BDK to obtain the basic land use map patch DL. The intersection operation of map patches is similar to the intersection of sets, that is, retaining the common parts of the two. The same-caliber database and the unopened urban / village database are databases statistically derived using different methods for the same area and the same land use category. Retaining the common parts of the two effectively ensures the accuracy of the data statistics.

[0076] S22: Filter the data on disaster-damaged farmland and calculate the disaster reduction amount. Please combine with... Figure 2 The specific calculation method for disaster damage reduction is as follows:

[0077] S221: Select disaster-damaged patches ZH from the same-caliber patch TKJ. Select unused farmland migration patches LW from the basic land use patch DL. The selection process is for a specific category of data. For example, mark farmland data with disaster damage markers in the same-caliber patch TKJ, and then select the corresponding disaster-damaged patch ZH based on the range of the disaster damage markers.

[0078] An intersection operation is performed on the disaster-damaged patch ZH and the unused cultivated land patch LW to obtain the disaster-damaged unused patch ZH-LW.

[0079] S222: Erase the non-protected cultivated land parcel FBH in the disaster-damaged unused parcel ZH-LW to obtain the disaster-damaged parcel ZH-KC. The erasure operation is similar to subtraction of sets, that is, deleting the part of the former that also belongs to the latter. By performing the erasure operation on the disaster-damaged unused parcel, removing the part that does not belong to protected cultivated land, and retaining the disaster-damaged cultivated land within the protected cultivated land area, the protection status of disaster-damaged cultivated land can be observed more intuitively and effectively.

[0080] The field ZH-KC of the disaster-damaged land parcel is read to calculate the area of ​​the damaged farmland. Calculating the farmland area from the parcel can be done by first using computer software for graphic recognition to obtain the area of ​​the parcel itself, and then enlarging it according to the parcel's scale to obtain the actual farmland area. Alternatively, only the parcel area can be calculated, reducing the computational workload without affecting the final score. In this embodiment, the calculated actual farmland area can be retained and recorded for further data analysis.

[0081] The area of ​​farmland destroyed by the disaster is S ZH This can be expressed as:

[0082] S ZH =S'[(LW∩ZH)-FBH]*β -2

[0083] LW=CH LW (DL), ZH=CH ZH (TKJ)

[0084] In the formula, S'[] is the area of ​​the map patch, LW is the farmland flowing to the unused map patch, ZH is the disaster-damaged map patch, FBH is the unprotected farmland map patch, β is the map patch scale, and CH() is the filtering operation.

[0085] S223: Summarize the data on the area of ​​farmland damaged by disaster and calculate the disaster reduction figure. The disaster reduction figure P is then calculated. ZH This can be expressed as:

[0086]

[0087] In the formula, k1 and k2 are manually set coefficients, and σ(S ZH ) represents the standard deviation of the area of ​​farmland destroyed by disaster, and Avg() represents the average value. S is the target control area for disaster recovery. ZHTK The actual controlled area for disaster recovery.

[0088] in, N represents the statistical years, i = 1, 2, ..., N. This represents the average area of ​​farmland destroyed by disaster.

[0089] S23: Filter the data on farmland converted from cultivation and calculate the reduction in cultivated land area. Please combine this with... Figure 3 The specific calculation method for the reduction of farmland due to land conversion is as follows:

[0090] S231: Select LCS (Local Land Flow to Forest / Grassland Wetland) patches from the land use maps. Perform an intersection operation on the LCS of cultivated land flowing to forest / grassland wetland patches and the TG (Ecological Return to Cultivated Land) patches to obtain the LCS-TG of forest / grassland wetland return to cultivated land patches.

[0091] S232: In the forest-grassland wetland conversion map patch LCS-TG, erase the non-protected cultivated land patch FBH to obtain the converted cultivated land patch TG-KC. Read the fields of the converted cultivated land patch TG-KC and calculate the converted cultivated land area. Then the converted cultivated land area S TG This can be expressed as:

[0092] S TG =S'[(TG∩LCS)-FBH]*β -2

[0093] LCS=CH LCS (DL)

[0094] In the formula, TG represents the ecological conversion of farmland to forest and grassland wetland plots, and LCS represents the conversion of farmland to forest and grassland wetland plots.

[0095] S233: Summarize the data on the area of ​​farmland converted from cultivation and calculate the reduction in farmland area. The reduction in farmland area is P. TG This can be expressed as:

[0096]

[0097] In the formula, σ(S) TG ) represents the standard deviation of the area of ​​farmland converted from cultivated land. S is the target area for ecological restoration through farmland conversion.ZHTK The actual controlled area for ecological restoration of farmland.

[0098] S24: Filter the cultivated land data for the three zones and three lines, and calculate the outflow reduction amount for the three zones and three lines. Please combine this with... Figure 4 The specific calculation method for the outflow reduction in the three zones and three lines is as follows:

[0099] S241: Erase the town / village map patch CZC from the basic land cover patch DL to obtain the town / village erased map patch CZC-CC. Select the town / village non-cultivated map patch FG from the town / village erased map patch CZC-CC.

[0100] S242: Erase the permanent cultivated land patch SSY of the Three Zones and Three Lines in the cultivated land protection patch SSB to obtain the general cultivated land patch SS-YG of the Three Zones and Three Lines.

[0101] S243: Perform an intersection operation on the non-cultivated land parcels FG in urban and rural areas and the general cultivated land parcels SS-YG along the "Three Zones and Three Lines" to obtain the cultivated land parcels SSL flowing out of the "Three Zones and Three Lines". Read the fields of the cultivated land parcels SSL flowing out of the "Three Zones and Three Lines" and calculate the area of ​​cultivated land flowing out of the "Three Zones and Three Lines". Then the area of ​​cultivated land flowing out of the "Three Zones and Three Lines" is S. SSL This can be expressed as:

[0102] S SSL =S'[FG∩(SSB-SSY)]*β -2

[0103] FG = CH FG (DL-CZC)

[0104] In the formula, FG represents non-arable land parcels in urban and rural areas.

[0105] S244: Summarize the data on the area of ​​cultivated land flowing out of the "Three Zones and Three Lines" and calculate the reduction amount of cultivated land outflow from the "Three Zones and Three Lines". Then, the reduction amount P of cultivated land outflow from the "Three Zones and Three Lines" is calculated. GTK This can be expressed as:

[0106]

[0107] In the formula, σ(S) SSL The standard deviation of the area of ​​cultivated land flowing out of the three regions and three lines is denoted as . The target control area for the outflow of arable land in the three zones and three lines, S SSLTK This refers to the actual area of ​​cultivated land flowing out of the "three zones and three lines" (referring to specific administrative regions and economic zones).

[0108] S3: Obtain data on changes in farmland protection and stable farmland use, combine them with comprehensive control data, set comprehensive control weights, farmland protection weights and stable farmland use weights according to the degree of human controllability, and calculate farmland protection assessment scores using a weighted calculation method.

[0109] The specific calculation method for the farmland protection assessment score is as follows:

[0110] S31: Standardize farmland protection data. For example, set a target area of ​​S for farmland protection for disaster-damaged farmland. ZHA After investing cost C ZH After artificial adjustments, the actual area of ​​disaster-damaged farmland protected was S. ZHB Based on the unified cost C0, the target protected arable land area for disaster-damaged arable land under the unified standard is calculated as follows:

[0111]

[0112] The actual protected arable land area under the unified standard for disaster-damaged arable land is:

[0113] SD ZH =S ZHB *C0 / C ZH

[0114] Correspondingly, the target protected area of ​​cultivated land under the unified standard for returning cultivated land to forest is:

[0115]

[0116] The actual protected arable land area under the unified standard for returning farmland to forest is:

[0117] SD TG =S TGB *C0 / C TG

[0118] The target protected arable land area under unified standards for outflowed arable land is:

[0119]

[0120] The actual protected arable land area under the unified standard for outflowed arable land is:

[0121] SD SSL =S SSLB *C0 / C SSL

[0122] In the formula, S TGA S TGB These are the target protected arable land area and the actual protected arable land area, respectively, S. SSLA S SSLB C represents the target protected arable land area and the actual protected arable land area, respectively, for arable land outflowing from the area. TG C SSL These are the costs invested in the conversion of farmland to non-farmland and the transfer of farmland from other areas.

[0123] S32: Calculate the degree of human controllability under a uniform standard. The formula for calculating the degree of human controllability is:

[0124]

[0125] In the formula, R ZH The `avg()` function calculates the average value of farmland protected from disasters, representing the degree of human control over its preservation. To protect the target area of ​​disaster-damaged farmland under a unified standard, SD ZH To standardize the actual protected arable land area of ​​disaster-damaged farmland, R TG The degree of human control over the conversion of farmland to protected farmland. To protect the target area of ​​cultivated land under a unified standard for returning farmland to forest, SD TG To determine the actual protected arable land area under a unified standard for returning farmland to forest, R SSL The degree of human control over the outflow of arable land in the "three zones and three lines" (referring to areas with high population density and rural development zones). To protect the target area of ​​cultivated land outflowing under unified standards, SD SSL The actual protected arable land area under a unified standard for arable land outflow.

[0126] S33: Calculate the disaster weight ω based on the degree of human controllability. ZK Weight ω of land conversion TG Three-zone, three-line outflow weight ω SSL The formula is expressed as:

[0127] ω ZH =R ZH / (R ZH +R TG +R SSL )

[0128] ω TG =R TG / (R ZH +R TG +R SSL )

[0129] ω SSL =R SSL / (R ZH +R TG +R SSL )

[0130] The higher the degree of human control, the greater the influence on the overall regulation of arable land, and therefore the higher its weight in the assessment score.

[0131] S34: Calculate the farmland protection assessment score P based on the set weights and the control parameters. S The formula is expressed as:

[0132] PS =ω ZH P ZH +ω TG P TG +ω SSL P SSL .

[0133] The evaluation method in this embodiment comprehensively considers farmland protection caused by natural disasters, ecological restoration, and farmland outflow due to the "Three Zones and Three Lines" policy. By performing intersection and erasure operations on multiple databases, it obtains data on farmland protected by natural disasters, ecological restoration, and farmland outflow due to the "Three Zones and Three Lines" policy. Compared to directly statistically analyzing relevant farmland data, this embodiment demonstrates more interconnected database processing, improves data utilization, and achieves higher statistical accuracy, thereby enhancing the objectivity of the evaluation score. By splitting the data, it not only improves processing efficiency but also makes the data clearer, less prone to data confusion, and increases verification accuracy. Furthermore, with the assistance of data algorithms, it can accelerate the processing of large amounts of data, employing more optimized algorithms and models to process the data, thus achieving efficient and accurate quality assessment of the results from multiple databases.

[0134] Furthermore, based on data from disaster-related reductions, land conversion reductions, and outflow reductions from the "three zones and three lines" (referring to areas with concentrated agricultural land use and forestry), a data aggregation and calculation method is employed. Weights are assigned to each reduction data point, and individual assessments are conducted. Finally, a weighted calculation method is used to obtain the final assessment score for cultivated land preservation. This embodiment's assessment method not only refines the scoring rules for cultivated land protection categories and fully considers actual changes in cultivated land protection, but also sets weights for each protection category based on the degree of human controllability. This allows for a more objective and accurate evaluation of cultivated land protection status and provides important guidance for setting assessment standards for cultivated land protection.

[0135] Example 2: Please refer to Figure 5 This embodiment provides a system for assessing and evaluating arable land reserves, which can be controlled using the method for assessing and evaluating arable land reserves described in Embodiment 1. The system includes a data acquisition module, a preprocessing module, a map processing module, and a calculation module.

[0136] The data acquisition module is used to collect data from a large database based on farmland protection. This large database includes databases with consistent coverage, closed databases for urban and rural areas, non-protected databases, ecological farmland conversion databases, urban and rural land use databases, farmland protection databases for the "Three Zones and Three Lines" (referring to designated areas and three protection zones), and permanent farmland databases for the "Three Zones and Three Lines."

[0137] The preprocessing module is used for: 1. Converting the large database into multiple land use patches with consistent scale. These are respectively denoted as patches of the same caliber, urban and rural patches (not opened), non-protected cultivated land patches, ecological conversion of cultivated land patches, urban and rural land use patches, cultivated land protection patches for the "Three Zones and Three Lines" and permanent cultivated land patches for the "Three Zones and Three Lines".

[0138] Second, standardize the processing of farmland protection data.

[0139] The patch processing module is used to perform intersection, erasure, and filtering operations on different patches in order to retain the required patch data.

[0140] In this embodiment, the patch processing module is used to filter data on disaster-damaged farmland, farmland converted from cultivation to forest, and farmland in the "Three Zones and Three Lines" region. The specific operation process is as follows:

[0141] Filtering data on disaster-damaged farmland: Disaster-damaged patches are selected from patches of the same caliber, and unused farmland due to land flow are selected from basic land use patches. Intersection is performed on the disaster-damaged patches and the unused farmland due to land flow patches to obtain disaster-damaged unused patches. Non-protected farmland patches are erased from the disaster-damaged unused patches to obtain patches after disaster deduction. Fields of the disaster-damaged patches are read to calculate the area of ​​disaster-damaged farmland.

[0142] Filtering data on converted farmland: Filter farmland parcels from land use maps to identify those that flow from farmland to forest / grassland wetlands. Perform an intersection operation on these farmland-to-forest / grassland wetland parcels and the ecological conversion farmland parcels to obtain forest / grassland wetland conversion farmland parcels. Erase non-protected farmland parcels from these forest / grassland wetland conversion farmland parcels to obtain the parcels after farmland conversion deduction. Read the fields of the parcels after farmland conversion deduction and calculate the area of ​​converted farmland.

[0143] Filtering cultivated land data under the "Three Zones and Three Lines" framework: Erase urban and rural land use patches from the basic land category map to obtain erased urban and rural patches. Filter out non-cultivated urban and rural patches from these erased patches. Erase permanent cultivated land patches under the "Three Zones and Three Lines" cultivated land protection map to obtain general cultivated land patches under the "Three Zones and Three Lines" framework. Perform an intersection operation between the non-cultivated urban and rural patches and the general cultivated land patches under the "Three Zones and Three Lines" framework to obtain cultivated land outflow patches under the "Three Zones and Three Lines" framework. Read the fields of the outflow cultivated land patches under the "Three Zones and Three Lines" framework and calculate the outflow area of ​​cultivated land under the "Three Zones and Three Lines" framework.

[0144] The calculation module is used for: 1. Calculating the reduction amount of disaster-damaged farmland based on the selected disaster-damaged farmland data. 2. Calculating the reduction amount of farmland converted from cultivation based on the selected farmland conversion data. 3. Calculating the reduction amount of farmland outflow from the "Three Zones and Three Lines" based on the selected farmland data from the "Three Zones and Three Lines" (referring to the three zones and three lines of farmland protection). 4. Calculating the degree of human controllability based on the standardized farmland protection data. 5. Calculating the disaster damage weight, farmland conversion weight, and outflow weight from the "Three Zones and Three Lines" based on the degree of human controllability. 6. Performing a weighted calculation on the disaster damage reduction amount, farmland conversion reduction amount, and outflow weight from the "Three Zones and Three Lines" based on the disaster damage weight, farmland conversion weight, and outflow weight to obtain the farmland protection assessment score. The calculation process of each step of the calculation module is shown in Example 1 and will not be elaborated here.

[0145] The technical features of the above embodiments can be combined in any way. For the sake of brevity, 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.

[0146] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the invention patent. 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 all fall within the protection scope of the present invention. Therefore, the protection scope of this invention patent should be determined by the appended claims.

Claims

1. A method for assessing and evaluating arable land area, characterized in that, Includes the following steps: S1: Obtain a large database based on farmland protection; S2: Filter out data on disaster-damaged farmland, farmland converted from farmland, and farmland outflowing from the "Three Zones and Three Lines" in the large database, and calculate the reduction figures for disaster-damaged farmland, farmland converted from farmland, and farmland outflowing from the "Three Zones and Three Lines" respectively; S3: The farmland protection assessment score is calculated using a weighted calculation method; the specific method is as follows: S31: Standardize the processing of farmland protection data; S32: Calculate the degree of human controllability based on standardized farmland protection data; S33: Calculate the disaster damage weight, farmland conversion weight, and outflow weight in the three zones and three lines according to the degree of human control. S34: The reduction figures for disaster damage, the reduction figures for land conversion from farmland to forest, and the reduction figures for land outflow from the three zones and three lines are weighted according to the weights to obtain the farmland protection assessment score.

2. The method for assessing and evaluating arable land area according to claim 1, characterized in that, In step S2, the method for calculating the adjustment reduction number is as follows: S21: Preprocess the large database to obtain multiple land use patches with consistent proportions, namely, patches of the same diameter, urban and rural patches without opening, non-protected arable land patches, ecological conversion of arable land patches, urban and rural land use patches, arable land protection patches for the three zones and three lines, and permanent arable land patches for the three zones and three lines; perform an intersection operation on the patches of the same diameter and the urban and rural patches without opening to obtain basic land use patches; S22: Filter disaster-damaged farmland data, calculate the standard deviation of disaster-damaged farmland and the controlled average value of disaster-damaged protected farmland, and calculate the disaster reduction amount based on the set coefficients; then the disaster reduction amount P... ZH Expressed as: ● P ZH =k1σ(S ZH )+k2Avg(|S ZHTK -S ZHTK |) In the formula, k1 and k2 are manually set coefficients, and σ(S ZH ) represents the standard deviation of the area of ​​farmland destroyed by disaster, and Avg() represents the average value. ● S ZHTK S is the target control area for disaster recovery. ZHTK The actual controlled area for disaster recovery; S23: Filter the data on farmland converted from cultivation, calculate the standard deviation of the converted farmland and the average value of the controlled conversion of farmland for protection, and calculate the reduction amount of farmland converted from cultivation based on the set coefficients; then the reduction amount P is calculated. TG Expressed as: ● P TG =k1σ(S TG )+k2Avg(|S TGTK -S TGTK |) In the formula, σ(S) TG ) represents the standard deviation of the area of ​​farmland converted from cultivated land. ● S ZHTK S is the target area for ecological restoration through farmland conversion. ZHTK The actual controlled area for ecological restoration of farmland converted from cultivated land; S24: Filter the cultivated land data of the "Three Zones and Three Lines" (referring to the three zones and three lines of cultivated land outflow), calculate the standard deviation of cultivated land outflow from the "Three Zones and Three Lines" and the average value of cultivated land protection outflow from the "Three Zones and Three Lines," and calculate the reduction amount of cultivated land outflow from the "Three Zones and Three Lines" based on the set coefficients; then calculate the reduction amount P of cultivated land outflow from the "Three Zones and Three Lines." GTK Expressed as: ● P GTK =k1σ(S SSL )+k2Avg(|S SSLTK -S SSLTK |) In the formula, σ(S) SSL The standard deviation of the area of ​​cultivated land flowing out of the three regions and three lines is denoted as . ● S SSLTK The target control area for the outflow of arable land in the three zones and three lines, S SSLTK This refers to the actual area of ​​cultivated land flowing out of the "three zones and three lines" (referring to specific administrative regions and economic zones).

3. The method for assessing and evaluating arable land area according to claim 2, characterized in that, In step S22, the method for filtering disaster-damaged farmland data is as follows: S221: Select disaster-damaged patches from patches of the same caliber, and select unused farmland flow patches from basic land use patches; perform an intersection operation on disaster-damaged patches and unused farmland flow patches to obtain disaster-damaged unused patches. S222: Erase non-protected cultivated land parcels from the disaster-damaged unused parcels to obtain the disaster-damaged parcels; read the fields of the disaster-damaged parcels and calculate the disaster-damaged cultivated land area; then the disaster-damaged cultivated land area S ZH This can be expressed as: S ZH =S'[(LW∩ZH)-FBH]*β -2 LW=CH LW (DL),ZH=CH ZH (TKJ) In the formula, S'[] is the area of ​​the map patch, LW is the farmland flow to unused map patch, ZH is the disaster-damaged map patch, FBH is the unprotected farmland map patch, β is the map patch scale, DL is the basic land type map patch, and CH() is the filtering operation.

4. The method for assessing and evaluating arable land area according to claim 2, characterized in that, In step S23, the method for filtering the data on converted farmland is as follows: S231: Select farmland-to-forest-grass-wetland patches from land use patches; perform an intersection operation on farmland-to-forest-grass-wetland patches and ecological conversion of farmland to obtain forest-grass-wetland conversion of farmland patches; S232: Erase the non-protected cultivated land patches in the forest-grassland wetland conversion map to obtain the converted cultivated land patches; read the fields of the converted cultivated land patches and calculate the converted cultivated land area; then the converted cultivated land area S TG Expressed as: S TG =S'[(TG∩LCS)-FBH]*β -2 LCS=CH LCS (DL) In the formula, TG represents the ecological conversion of farmland to forest and grassland wetland plots, and LCS represents the conversion of farmland to forest and grassland wetland plots.

5. The method for assessing and evaluating arable land area according to claim 2, characterized in that, In step S24, the calculation method for filtering the cultivated land data of the three zones and three lines is as follows: 241: Erase urban and rural map patches from the basic land use map patches to obtain urban and rural map patches after erasure; select non-arable urban and rural map patches from the urban and rural map patches after erasure; 242: Erase the permanent cultivated land patches of the Three Zones and Three Lines in the cultivated land protection patch to obtain the general cultivated land patches of the Three Zones and Three Lines; 243: Perform an intersection operation on non-cultivated map patches in urban and rural areas and general cultivated map patches in the "Three Zones and Three Lines" to obtain cultivated map patches flowing out of the "Three Zones and Three Lines". Read the fields of cultivated land outflow from the "Three Zones and Three Lines" data, and calculate the area of ​​cultivated land outflowing from the "Three Zones and Three Lines" data; then the cultivated land outflow area S from the "Three Zones and Three Lines" data is calculated. SSL Expressed as: S SSL =S’[FG∩(SSB-SSY)]*β -2 FG=CH FG (DL-CZC) In the formula, FG represents non-arable land parcels in urban and rural areas, SSB represents protected arable land parcels in the three zones and three lines, SSY represents permanent arable land parcels in the three zones and three lines, and CZC represents arable land parcels used in urban and rural areas.

6. The method for assessing and evaluating arable land area according to claim 1, characterized in that, In step S32, the formula for calculating the degree of human controllability is: In the formula, R ZH The `avg()` function calculates the average value of farmland protected from disasters, representing the degree of human control over its preservation. To protect the target area of ​​disaster-damaged farmland under a unified standard, SD ZH To standardize the actual protected arable land area of ​​disaster-damaged farmland, R TG To ensure the degree of human control over the conversion of farmland to protected arable land, S ● D TG To protect the target area of ​​cultivated land under a unified standard for returning farmland to forest, SD TG To determine the actual protected arable land area under a unified standard for returning farmland to forest, R SSL To assess the degree of human control over the outflow of arable land in the "Three Zones and Three Lines" area, S ● D SSL To protect the target area of ​​cultivated land outflowing under unified standards, SD SSL The actual protected arable land area under a unified standard for arable land outflow.

7. The method for assessing and evaluating arable land area according to claim 6, characterized in that, In step S33, the disaster weight ω ZK Weight of land conversion ω TG Three-zone, three-line outflow weight ω SSL The calculation formulas are as follows: ω ZH =R ZH / (R ZH +R TG +R SSL ) ω TG =R TG / (R ZH +R TG +R SSL ) ω SSL =R SSL / (R ZH +R TG +R SSL )。 8. The method for assessing and evaluating arable land area according to claim 7, characterized in that, In step S34, the farmland protection assessment score P S The calculation formula is: P S =ω ZH P ZH +oh TG P TG +oh SSL P SSL In the formula, P ZH P TG P SSL These are the reduction figures for disaster damage, reduction figures for land conversion from farmland to forest, and reduction figures for outflow from the "three regions and three lines" (referring to areas with high population density and rural development).

9. A system for assessing and evaluating arable land area, comprising the method for assessing and evaluating arable land area as described in any one of claims 1 to 8, characterized in that, The performance evaluation system includes: The data acquisition module is used to collect data from a large database based on farmland protection. The preprocessing module is used for:

1. converting the large database into multiple land use patches with consistent proportions; 2. standardizing the farmland protection data; The patch processing module is used to perform intersection, erasure, and filtering operations on different patches in order to retain the required patch data; The calculation module is used for:

1. Calculating the amount of disaster-damaged farmland to be reduced based on the selected disaster-damaged farmland data; 2. Calculating the amount of farmland converted from cultivation to other uses based on the selected farmland conversion data; 3. Calculating the amount of farmland outflow from the "Three Zones and Three Lines" based on the selected farmland data from the "Three Zones and Three Lines"; 4. Calculating the degree of human controllability based on the standardized farmland protection data; 5. Calculating the disaster damage weight, farmland conversion weight, and farmland outflow weight from the "Three Zones and Three Lines" based on the degree of human controllability; 6. Performing a weighted calculation on the amount of disaster-damaged farmland to be reduced, farmland conversion to other uses, and farmland outflow from the "Three Zones and Three Lines" based on the disaster damage weight, farmland conversion weight, and farmland outflow weight to obtain the farmland protection assessment score.

Images