Land resource checking method and system

Abstract

The invention discloses a land resource checking method and system, and the method comprises the steps: firstly, determining a checking range of land resource checking from a land class pattern spot element layer in the land utilization data of an annual land change survey result, and according to the checking range, carrying out the checking of the land resources from the quality grade evaluation result data related to a land class to be checked; and a grading layer corresponding to the checking range is extracted. Then, performing spatial overlapping on the land class pattern spots of the land class to be checked and the grading pattern layers; and finally, judging whether a same land class grading unit overlapped with the land class pattern spot exists in the grading layer or not. And if not, comparing each feature element data of the land type pattern spot with the feature element data of each grading unit in the grading layer, selecting the grading unit of which the feature element is closest to the land type pattern spot, and endowing the land type pattern spot with utilization and other information corresponding to the grading unit. Therefore, batch processing of land resource checking is achieved, repeated operation is reduced, and the speed and efficiency of land resource checking are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of special types of information retrieval, and specifically relates to a method and system for land resource inventory. Background Art

[0002] Natural resources owned by the whole people refer to natural resources that are owned by the state, have scarcity, usefulness (including economic benefits, social benefits, and ecological benefits), and clear property rights. In 2024, the Ministry of Natural Resources issued the "Technical Guidelines for the Inventory of Natural Resources Owned by the Whole People", which stipulates the technical processes for the inventory of natural resources owned by the whole people, including natural resources assets of agricultural land owned by the whole people, and involves the inventory of physical quantities of assets, the inventory of use rights, the construction of price systems, the accounting of economic values, the construction of databases, etc. According to the requirements of the "Guidelines", the current inventory of land resources is mainly based on arcgis software, and data analysis is manually completed in a specific operation sequence. However, the data analysis process is fixed and has many steps, and manual operation is relatively cumbersome, resulting in slow inventory speed and low efficiency of land resources. Summary of the Invention

[0003] In view of the deficiencies of the prior art, the present invention proposes a method and system for land resource inventory, which can improve the speed and efficiency of land resource inventory. The specific technical solutions are as follows:

[0004] In the first aspect, a method for land resource inventory is provided. In the first implementable manner of the first aspect, it includes:

[0005] Obtain the land type patches and classification layers of the land types to be inventoried within the inventory scope, and perform spatial overlay on the land type patches and classification layers;

[0006] Determine whether there are same-land-type classification units overlapping with the land type patches in the classification layer;

[0007] In response to the absence of same-land-type classification units overlapping with the land type patches, obtain various characteristic element data corresponding to each same-land-type classification unit in the land type patches and classification layers;

[0008] By comparing the various characteristic element data corresponding to the land type patches and each classification unit, determine the classification unit with the closest characteristic elements to the land type patches, and assign the utilization grade information of this classification unit to the land type patches.

[0009] Combined with the first implementable manner of the first aspect, in the second implementable manner of the first aspect, it further includes:

[0010] In response to the existence of same-land-type classification units overlapping with the land type patches, assign the utilization grade information of the classification unit with the largest overlapping area with the land type patches to the land type patches.

[0011] Combined with the first or second implementation manner of the first aspect, in the third implementation manner of the first aspect, determining the grading unit with the closest characteristic elements to the land type patches includes:

[0012] According to the corresponding characteristic element data, respectively calculate the comprehensive grading scores corresponding to each of the grading units by using a grading algorithm;

[0013] Compare the comprehensive grading scores corresponding to each of the grading units, and select the grading unit corresponding to the maximum comprehensive grading score as the grading unit with the closest characteristic elements to the land type patches.

[0014] Combined with the third implementation manner of the first aspect, in the fourth implementation manner of the first aspect, calculating the comprehensive grading score corresponding to the grading unit includes:

[0015] According to the corresponding characteristic element data, respectively assign grading values corresponding to the characteristic elements of the grading unit by using corresponding grading assignment methods;

[0016] Combined with the grading values corresponding to each of the characteristic elements and the weights, calculate the comprehensive grading score corresponding to the grading unit by using a weighted algorithm.

[0017] Combined with the fourth implementation manner of the first aspect, in the fifth implementation manner of the first aspect, when the land type to be inventoried is cultivated land, the grading values corresponding to the characteristic elements of the grading unit are assigned as follows:

[0018] According to the characteristic element data corresponding to altitude, slope aspect, and / or slope, respectively calculate the data differences of altitude, slope aspect, and / or slope between the land type patches and the grading units;

[0019] Respectively match the data differences of altitude, slope aspect, and / or slope with the corresponding grading intervals;

[0020] According to the assignment rules of the grading intervals to which the data differences belong, respectively assign the grading values corresponding to altitude, slope aspect, and / or slope of the grading unit.

[0021] Combined with the fifth implementation manner of the first aspect, in the sixth implementation manner of the first aspect, the assignment rules include:

[0022] When the data difference corresponding to altitude belongs to (0, H g ), assign the grading value corresponding to altitude of the grading unit according to the distance attenuation algorithm;

[0023] And / or, when the data difference corresponding to slope belongs to (0, S g ), assign the grading value corresponding to slope of the grading unit according to the distance attenuation algorithm;

[0024] Among them, H g and S g are respectively the maximum differences in altitude and slope corresponding between the land use patches and the gradation units.

[0025] Combined with the fourth implementation manner of the first aspect, in the seventh implementation manner of the first aspect, when the land use type to be inventoried is arable land, the grading assignment corresponding to the characteristic elements of the gradation units includes:

[0026] Determine the maximum straight-line distance D x1max and the minimum straight-line distance D x1min between the land use patch and all gradation units according to the characteristic element data corresponding to the shortest distance;

[0027] When the shortest distance D x1y1 between the land use patch and the gradation unit belongs to (D x1min , D x1max ), assign the grading assignment corresponding to the shortest distance of the gradation unit according to the following calculation formula:

[0028]

[0029] Combined with the fourth implementation manner of the first aspect, in the eighth implementation manner of the first aspect, when the land use type to be inventoried is garden land, the grading assignment corresponding to the characteristic elements of the gradation units includes:

[0030] Calculate the data differences in road accessibility and / or labor distance between the land use patch and the gradation units respectively according to the characteristic element data corresponding to road accessibility and / or labor distance;

[0031] Match the data differences in road accessibility and / or labor distance with the corresponding grading intervals respectively;

[0032] Assign the grading assignments corresponding to the road accessibility and / or labor distance of the gradation units respectively according to the assignment rules of the grading intervals to which the data differences belong.

[0033] Combined with the eighth implementation manner of the first aspect, in the ninth implementation manner of the first aspect, the assignment rules include:

[0034] When the data difference corresponding to road accessibility belongs to (0, R g ), assign the grading assignment corresponding to road accessibility of the gradation unit according to the distance decay algorithm;

[0035] And / or, when the data difference corresponding to labor distance belongs to (0, W g ), assign the grading assignment corresponding to labor distance of the gradation unit according to the distance decay algorithm;

[0036] Among them, R g and W g are respectively the maximum differences in road accessibility and labor distance between the land use type patches and the grading units.

[0037] In a second aspect, a land resource inventory system is provided, including:

[0038] An overlay module configured to obtain the land use type patches and the grading layers of the land to be inventoried within the inventory scope, and perform spatial overlay on the land use type patches and the grading layers;

[0039] An analysis module configured to determine whether there are grading units of the same land use type that overlap with the land use type patches in the grading layer;

[0040] In response to the absence of grading units of the same land use type that overlap with the land use type patches, obtain various characteristic element data corresponding to each grading unit of the same land use type in the land use type patches and the grading layers;

[0041] By comparing the various characteristic element data corresponding to the land use type patches and each grading unit, determine the grading unit with the closest characteristic elements to the land use type patches, and assign the utilization grade information of the grading unit to the land use type patches.

[0042] Beneficial effects: By using the land resource inventory method and system of the present invention, the land use type patches of the land to be inventoried within the inventory scope can be batch spatially overlapped with the grading layers in the quality grade evaluation results of the same land use type. Through the spatial overlap results, the grading unit closest to each land use type patch to be inventoried can be determined, and the national utilization grade information or orchard grade information of the closest grading unit can be automatically assigned to each land use type patch to be inventoried to determine the quality grade corresponding to each land use type patch, thereby realizing the batch processing of land resource inventory, reducing repeated operations, and improving the speed and efficiency of land resource inventory. Description of the Drawings

[0043] In order to more clearly illustrate the specific embodiments of the present invention, the drawings required for the specific embodiments will be briefly introduced below. In all the drawings, the components or parts are not necessarily drawn according to the actual ratio.

[0044] Figure 1 It is a flowchart of the land resource inventory method provided by an embodiment of the present invention;

[0045] Figure 2 It is a system block diagram of the land resource inventory system provided by an embodiment of the present invention. Detailed Embodiments

[0046] The embodiments of the technical solution of the present invention will be described in detail below with reference to the accompanying drawings. The following embodiments are only used to illustrate the technical solution of the present invention more clearly, so they are only examples and cannot be used to limit the protection scope of the present invention.

[0047] It should be understood that a "land type patch" refers to a single land type plot divided by administrative boundaries, ownership boundaries, and single linear land features.

[0048] It should also be understood that "national utilization grade, etc." is a grade division of land resources based on various factors such as the natural conditions, economic demands, policy orientations, and utilization status of the land. At present, according to the investigation and evaluation work of the grade of cultivated land quality, the cultivated land in the country is divided into 15 quality grades, among which the cultivated land of grade 1 has the best quality and that of grade 15 is the worst. The cultivated land in the country is usually divided into high-quality land, high-grade land, medium-grade land, and low-grade land according to grades 1-4, 5-8, 9-12, and 13-15. The grade of cultivated land directly affects the economic benefits of agricultural production and crop yields.

[0049] The "grading layer", also known as the "grading unit layer", refers to a layer that divides land into different grades or types of grading units according to specific methods and standards in land resource management or evaluation. This layer usually also includes relevant index information for grade evaluation. The grading layer of cultivated land comes from the annual updated evaluation results of the grade of cultivated land quality, including information such as the land type, ownership, location, economic coefficient, national utilization grade, and grading index of cultivated land. The grading layer of orchards comes from the orchard grading results, including information such as the land type, ownership, location, national orchard grade, provincial orchard grade, and grading index of orchards.

[0050] Embodiment 1

[0051] As Figure 1 shown in the flowchart of the land resource inventory method, the inventory method includes:

[0052] Step 1: Obtain the land type patches and grading layers of the land types to be inventoried within the inventory scope, and perform spatial superposition on the land type patches and grading layers;

[0053] Step 2: Determine whether there are grading units of the same land type overlapping with the land type patches in the grading layer;

[0054] In response to the absence of grading units of the same land type overlapping with the land type patches, obtain various characteristic element data corresponding to each grading unit of the same land type in the land type patches and grading layers;

[0055] By comparing the land type patches and various characteristic element data corresponding to each grading unit, determine the grading unit with the closest characteristic elements to the land type patches, and assign the utilization grade information of this grading unit to the land type patches.

[0056] Specifically, first, the inventory scope of this land resource inventory can be determined from the land use data in the annual land change survey results, specifically from the land type polygon element layer. And according to the inventory scope, the grading layer corresponding to the inventory scope can be extracted from the quality grading evaluation result data related to the land types to be inventoried. For example, if the land type to be inventoried is arable land, the grading layer corresponding to the inventory scope can be extracted from the annual updated evaluation result data of the arable land quality grading. If the land type to be inventoried is orchards, the grading layer corresponding to the inventory scope can be extracted from the orchard grading result data. Then, the land type polygons of the land types to be inventoried within the inventory scope can be spatially overlapped with the grading layer. Finally, it can be judged whether there are same-land-type grading units overlapping with the land type polygons in the grading layer according to the spatial overlap result. If not, the data of each characteristic element corresponding to the land type polygon can be compared with the data of the characteristic elements corresponding to each grading unit in the grading layer respectively, and the grading unit with the characteristic elements closest to the land type polygon can be selected, and the utilization grade information corresponding to this grading unit can be assigned to this land type polygon. Thus, the batch processing of land resource inventory is realized, reducing repeated operations and improving the speed and efficiency of land resource inventory.

[0057] In this embodiment, optionally, it further includes:

[0058] In response to the existence of a same-land-type grading unit overlapping with the land type polygon, the utilization grade information of the grading unit with the largest overlapping area with the land type polygon is assigned to the land type polygon.

[0059] Specifically, if there are multiple grading units in the grading layer that overlap with the land type polygon and have the same land category, the overlapping area between each grading unit and the land type polygon can be calculated, and the grading unit with the largest overlapping area can be selected, and the utilization grade information corresponding to this grading unit can be assigned to the land type polygon.

[0060] In this embodiment, optionally, in step 2, to determine the grading unit with the characteristic elements closest to the land type polygon, it includes:

[0061] According to the corresponding characteristic element data, the comprehensive grading scores corresponding to each grading unit are calculated respectively by using a grading algorithm;

[0062] The comprehensive grading scores corresponding to each grading unit are compared, and the grading unit corresponding to the largest comprehensive grading score is selected as the grading unit with the characteristic elements closest to the land type polygon.

[0063] Specifically, when selecting the classification unit that is closest to the characteristic elements of the land use patch, based on the data of multiple characteristic elements corresponding to the land use patch and each classification unit, such as slope, elevation, and shortest distance, a classification algorithm can be used to calculate the comprehensive classification scores corresponding to each classification unit respectively. Then, compare the comprehensive classification scores corresponding to each classification unit, and select the classification unit with the highest comprehensive classification score as the classification unit with the closest characteristic elements to the land use patch.

[0064] Since when determining the comprehensive classification score of the classification unit, multiple significant characteristic elements of the land are comprehensively considered, it is ensured that the characteristics such as the topography and land use status type of the selected classification unit are closest to the land use patch, improving the accuracy of the land resource inventory, and thus improving the speed and efficiency of the land resource inventory.

[0065] In this embodiment, optionally, calculating the comprehensive classification score corresponding to the classification unit includes:

[0066] According to the corresponding characteristic element data, respectively assign the classification values corresponding to the characteristic elements of the classification unit by using the corresponding classification and assignment methods;

[0067] Combined with the classification values and weights corresponding to each characteristic element, use a weighted algorithm to calculate the comprehensive classification score corresponding to the classification unit.

[0068] Specifically, when calculating the comprehensive classification score of the classification unit, based on the data of the characteristic elements of the land use patch and the classification unit, a classification and assignment method matching the characteristic elements can be used to respectively determine the classification value of the classification unit for each characteristic element. For example, based on the slope data of the land use patch and the classification unit, a classification and assignment method matching the slope can be used to determine the classification value of the classification unit for the slope. Then, based on the classification values and weights of the classification unit for each characteristic element, a weighted algorithm can be used to calculate the comprehensive classification score of the classification unit.

[0069] In this embodiment, optionally, when the land type to be inventoried is arable land, assigning the classification values corresponding to the characteristic elements of the classification unit includes:

[0070] According to the characteristic element data corresponding to elevation, slope aspect, and / or slope, respectively calculate the data differences of elevation, slope aspect, and / or slope between the land use patch and the classification unit;

[0071] Respectively match the data differences of elevation, slope aspect, and / or slope with the corresponding classification intervals;

[0072] According to the assignment rules of the classification interval to which the data difference belongs, respectively assign the classification values corresponding to elevation, slope aspect, and / or slope of the classification unit.

[0073] Specifically, in this embodiment, when the land type to be inventoried is arable land, altitude, slope aspect, slope gradient, and shortest distance can be selected as characteristic elements to calculate the comprehensive grading scores of grading units. Among them, when determining the grading assignments of grading units regarding altitude, slope aspect, and slope gradient, the altitude data difference, slope aspect data difference, and slope gradient data difference between the grading unit and the land type patch can be calculated respectively based on the altitude data, slope aspect data, and slope gradient data of the grading unit and the land type patch. Then, these three data differences are respectively matched with the grading intervals corresponding to altitude, slope aspect, and slope gradient, and the grading assignments of the grading unit regarding altitude, slope aspect, and slope gradient are respectively given according to the assignment rules corresponding to the matched grading intervals.

[0074] In this embodiment, optionally, the assignment rules include:

[0075] When the data difference corresponding to altitude belongs to (0, H g ), the grading assignment corresponding to altitude of the grading unit is given according to the distance decay algorithm;

[0076] And / or, when the data difference corresponding to slope gradient belongs to (0, S g ), the grading assignment corresponding to slope gradient of the grading unit is given according to the distance decay algorithm;

[0077] Among them, H g , S g are respectively the maximum differences corresponding to altitude and slope gradient between the land type patch and the grading unit.

[0078] Specifically, the assignment rules corresponding to altitude include:

[0079] When the altitude data difference between the grading unit and the land type patch satisfies H y -H x | = 0, H xys = 100;

[0080] When the altitude data difference between the grading unit and the land type patch satisfies H y -H x | = H g ), H xys = 0;

[0081] When the altitude data difference between the grading unit and the land type patch satisfies 0 < |H y -H x | < H g , the grading assignment corresponding to altitude of the grading unit is given by using the distance decay algorithm. The specific calculation formula is as follows:

[0082]

[0083] The assignment rules corresponding to slope gradient include:

[0084] When the slope data difference between the grading unit and the land use patch satisfies |S y - S x | = 0, S xys = 100;

[0085] When the slope data difference between the grading unit and the land use patch satisfies |S y - S x | = S g , S xys = 0;

[0086] When the slope data difference between the grading unit and the land use patch satisfies 0 < |S y - S x | < S g , the distance attenuation algorithm is used to assign the grading value corresponding to the slope to the grading unit. The specific calculation formula is as follows:

[0087]

[0088] The assignment rules corresponding to the slope aspect include:

[0089] It is set that the shady slope, semi-shady slope, flat land, semi-sunny slope, and sunny slope are 1, 2, 3, 4, and 5 respectively.

[0090] When the slope aspect data difference between the grading unit and the land use patch satisfies |O y - O x | = 0, O xys = 100;

[0091] When the slope aspect data difference between the grading unit and the land use patch satisfies |O y - O x | = 4, O xys = 0;

[0092] When the slope aspect data difference between the grading unit and the land use patch satisfies |O y - O x | = 1, O xys = 75;

[0093] When the slope aspect data difference between the grading unit and the land use patch satisfies |O y - O x | = 2, O xys = 50;

[0094] When the slope aspect data difference between the grading unit and the land use patch satisfies |O y - O x | = 3, O xys = 25.

[0095] Among them, H y, H x are the elevation data corresponding to the grading units and land use patches respectively, S y , S x are the slope data corresponding to the grading units and land use patches respectively, O y , O x are the aspect data corresponding to the grading units and land use patches respectively.

[0096] In this embodiment, the grading unit assignment rule for the shortest distance includes:

[0097] Determine the maximum straight-line distance D xmax and the minimum straight-line distance D xmin between the land use patch and all grading units according to the characteristic element data corresponding to the shortest distance;

[0098] When the shortest distance D xy between the land use patch and the grading unit belongs to (D x1min , D x1max ), assign the grading value corresponding to the shortest distance to the grading unit according to the following calculation formula:

[0099]

[0100] Specifically, when determining the grading value of the grading unit for the shortest distance, the maximum value and the minimum value can be selected from the shortest distances between the land use patch and each grading unit as the maximum straight-line distance D xmax and the minimum straight-line distance D xmin respectively. Then, compare the shortest distances between each grading unit and the land use patch with the maximum straight-line distance D xmax and the minimum straight-line distance D xmin respectively.

[0101] When the shortest distance data between the grading unit and the land use patch satisfies D xy / D xmax = 1, D xys = 0;

[0102] When the shortest distance data between the grading unit and the land use patch satisfies D xy / D xmin = 1, D xys = 100;

[0103] When the shortest distance data between the grading unit and the land use patch satisfies D xmin < D xy < D xmax , the specific calculation formula for the grading value of the grading unit corresponding to the shortest distance is as follows:

[0104]

[0105] After determining the grading values of the grading units for slope, aspect, altitude, and shortest distance, the comprehensive grading score of the grading unit can be calculated based on the grading values and weights corresponding to these four characteristic elements. The specific calculation formula is as follows:

[0106] Z xy = H xys ×H w + O xys ×O w + S xys ×S w + D xys ×D w ;

[0107] Among them, H w , O w , S w , D w are the weights corresponding to altitude, aspect, slope, and shortest distance respectively.

[0108] After calculating the comprehensive grading scores of all grading units for the cultivated land type patches, the grading unit with the highest comprehensive grading score can be selected from all grading units, and the national utilization grade information corresponding to this grading unit is assigned to the land type patches to complete the inventory of the cultivated land quality grades within the inventory scope.

[0109] Example 2

[0110] Example 2 is generally the same as Example 1, and its main difference is that when the land type to be inventoried is garden land, the grading values corresponding to the characteristic elements assigned to the grading unit include:

[0111] According to the characteristic element data corresponding to the road accessibility and / or labor distance, calculate the data differences of the road accessibility and / or labor distance between the land type patches and the grading units respectively;

[0112] Match the data differences of the road accessibility and / or labor distance with the corresponding grading intervals respectively;

[0113] According to the assignment rules of the grading intervals to which the data differences belong, assign the grading values corresponding to the road accessibility and / or labor distance of the grading units respectively.

[0114] Specifically, when inventorying the quality grades of garden land, when there are no grading units in the grading layer that overlap with the garden land type patches and have the same land category. The comprehensive grading score of the grading unit can be calculated using the road accessibility and labor distance as characteristic elements.

[0115] Specifically, first, based on the road accessibility data and labor distance data of the land use map patches and the grading units, the differences in road accessibility and labor distance between the land use map patches and the grading units are calculated respectively. Then, the differences in road accessibility and labor distance are respectively matched with the grading intervals corresponding to the road accessibility and labor distance. After that, according to the assignment rules corresponding to the matched grading intervals, the grading assignments of the grading units regarding road accessibility and labor distance are respectively given.

[0116] Combining the grading assignments of the grading units regarding road accessibility and labor distance, a weighted algorithm is used to calculate the comprehensive grading score of the grading units. Based on the comprehensive grading score, the grading unit closest to the garden land use map patch can be selected, and the corresponding national garden land grade information of the grading unit is assigned to the garden land use map patch.

[0117] In this embodiment, optionally, the assignment rules include:

[0118] When the data difference corresponding to the road accessibility belongs to (0, R g ), the grading assignment corresponding to the road accessibility of the grading unit is given according to the distance decay algorithm;

[0119] And / or, when the data difference corresponding to the labor distance belongs to (0, W g ), the grading assignment corresponding to the labor distance of the grading unit is given according to the distance decay algorithm;

[0120] Wherein, R g , W g are respectively the maximum differences in road accessibility and labor distance between the land use map patch and the grading unit.

[0121] Specifically, the assignment rules regarding road accessibility include:

[0122] When the road accessibility data difference satisfies |R y -R x | = 0, R xys = 100;

[0123] When the road accessibility data difference satisfies |R y -R x | = R g ), R xys = 0;

[0124] When the road accessibility data difference satisfies 0 < |R y -R x | < R g ), the grading assignment corresponding to the road accessibility of the grading unit is given by using the distance decay algorithm. The specific calculation formula is as follows:

[0125]

[0126] The assignment rules for the labor distance include:

[0127] When the difference in labor distance data satisfies |W y -W x | = 0, then W xys = 100;

[0128] When the difference in labor distance data satisfies |W y -W x | = W g , then R xys = 0;

[0129] When the difference in labor distance data satisfies 0 < |W y -W x | < W g , the distance attenuation algorithm is used to assign the corresponding graded values of road accessibility to the grading units. The specific calculation formula is as follows:

[0130]

[0131] After calculating the graded values of the labor distance and road accessibility for the grading units, the comprehensive graded score of the grading units for the garden land type patches can be calculated according to the graded values and weights corresponding to the labor distance and road accessibility. The specific calculation formula is as follows:

[0132] Z x ′ y = W xys ×W w + R xys ×R w ;

[0133] Among them, W w , R w are the weights corresponding to the labor distance and road accessibility respectively.

[0134] As Figure 2 shown in the system block diagram of the land resource inventory system, this inventory system includes:

[0135] An overlay module, configured to obtain the land type patches and grading layers of the land types to be inventoried within the inventory scope, and perform spatial overlay on the land type patches and grading layers;

[0136] An analysis module, configured to determine whether there are grading units of the same land type that overlap with the land type patches in the grading layer;

[0137] In response to the absence of grading units of the same land type that overlap with the land type patches, obtain the data of various characteristic elements corresponding to each grading unit of the same land type in the land type patches and grading layers;

[0138] By comparing the data of various characteristic elements corresponding to the land use patches and each grading unit, determine the grading unit with the closest characteristic elements to the land use patches, and assign the utilization grade information of this grading unit to the land use patches.

[0139] Specifically, the inventory system includes a superposition module and an analysis module. Among them, the superposition module can determine the inventory scope of this land resource inventory from the land use patch element layer in the land use data of the annual land use change survey results. And according to the inventory scope, extract the grading layer corresponding to the inventory scope from the data of the annual update and evaluation results of the cultivated land quality grade. The analysis module can perform spatial overlap between the land use patches of the land types to be inventoried within the inventory scope and the grading layer, and judge whether there are grading units of the same land type overlapping with the land use patches according to the spatial overlap result. If not, the data of various characteristic elements corresponding to the land use patches can be respectively compared with the data of the characteristic elements corresponding to each grading unit in the grading layer, select the grading unit with the closest characteristic elements to the land use patches, and assign the utilization grade information corresponding to this grading unit to the land use patches. Thus, batch processing of land resource inventory is realized, reducing repeated operations and improving the speed and efficiency of land resource inventory.

[0140] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than limiting it; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be covered by the scope of the claims and the description of the present invention.

Claims

1. A method for land resource inventory, characterized in that, Including: Obtain the land type patches and grading layers of the land types to be inventoried within the inventory scope, and perform spatial overlay on the land type patches and grading layers; Determine whether there are same-land-type grading units overlapping with the land type patches in the grading layer; In response to the non-existence of same-land-type grading units overlapping with the land type patches, obtain various characteristic element data corresponding to each same-land-type grading unit in the land type patches and grading layers; By comparing various characteristic element data corresponding to the land type patches and each grading unit, determine the grading unit with the closest characteristic elements to the land type patches, and assign the utilization grade information of the grading unit to the land type patches.

2. The land resource inventory method according to claim 1, characterized in that, It also includes: In response to the existence of same-land-type grading units overlapping with the land type patches, assign the utilization grade information of the grading unit with the largest overlapping area with the land type patches to the land type patches.

3. The land resource inventory method according to claim 1 or 2, characterized in that, Determining the grading unit with the closest characteristic elements to the land type patches includes: According to the corresponding characteristic element data, respectively calculate the comprehensive grading scores corresponding to each grading unit by using a grading algorithm; Compare the comprehensive grading scores corresponding to each grading unit, and select the grading unit corresponding to the largest comprehensive grading score as the grading unit with the closest characteristic elements to the land type patches.

4. The land resource inventory method according to claim 3, wherein Calculating the comprehensive grading score corresponding to the grading unit includes: According to the corresponding characteristic element data, respectively assign grading values corresponding to each characteristic element of the grading unit by using corresponding grading assignment methods; Combined with the grading values and weights corresponding to each characteristic element, calculate the comprehensive grading score corresponding to the grading unit by using a weighted algorithm.

5. The land resource inventory method according to claim 4, wherein, When the land type to be inventoried is cultivated land, the grading values corresponding to the characteristic elements of the grading unit include: According to the characteristic element data corresponding to altitude, slope direction, and / or slope, respectively calculate the data differences of altitude, slope direction, and / or slope between the land type patches and the grading unit; Respectively match the data differences of altitude, slope direction, and / or slope with the corresponding grading intervals; According to the assignment rules of the grading intervals where the data differences belong, respectively assign the grading values corresponding to altitude, slope direction, and / or slope of the grading unit.

6. The land resource inventory method according to claim 5, characterized in that, The assignment rules include: When the data difference corresponding to the altitude belongs to (0, H g ), the grading value corresponding to the altitude of the grading unit is assigned according to the distance attenuation algorithm; and / or, when the data difference corresponding to the slope belongs to (0, S g ), assign the graded value corresponding to the slope of the grading unit according to the distance attenuation algorithm; Among them, H g and S g are respectively the maximum differences corresponding to elevation and slope between land type patches and grading units.

7. The land resource inventory method according to claim 4, characterized in that When the land type to be inventoried is cultivated land, the grading values corresponding to the characteristic elements of the grading unit include: Determine the maximum linear distance D and the minimum linear distance D between the land use patch and all the grading units according to the characteristic element data corresponding to the shortest distance x1max and x1min ; The shortest distance D between the local land type patches and the grading units x1y1 Belonging to (D x1min , D x1max ), the grading assignment corresponding to the shortest distance of the grading unit is given according to the following calculation formula:

8. The land resource inventory method according to claim 4, characterized in that, When the land type to be inventoried is garden land, the grading values corresponding to the characteristic elements of the grading unit include: According to the characteristic element data corresponding to road accessibility and / or labor distance, respectively calculate the data differences of road accessibility and / or labor distance between the land type patches and the grading unit; Respectively match the data differences of road accessibility and / or labor distance with the corresponding grading intervals; According to the assignment rules of the grading intervals where the data differences belong, respectively assign the grading values corresponding to road accessibility and / or labor distance of the grading unit.

9. The land resource inventory method according to claim 8, wherein, The assignment rules include: When the data difference corresponding to the road accessibility belongs to (0, R g ), the grading assignment corresponding to the road accessibility of the grading unit is given according to the distance decay algorithm; and / or, when the data difference corresponding to the labor distance belongs to (0, W g ), assign the grading value corresponding to the labor distance of the grading unit according to the distance attenuation algorithm; Among them, R g and W g are respectively the maximum differences in road accessibility and labor distance between land use patches and grading units.

10. A land resource inventory system, characterized in that, Including: An overlay module configured to obtain the land type patches and grading layers of the land types to be inventoried within the inventory scope, and perform spatial overlay on the land type patches and grading layers; An analysis module, configured to determine whether there are equalization units of the same land type that overlap with the land type patches in the equalization layer; In response to the absence of equalization units of the same land type that overlap with the land type patches, obtain various characteristic element data corresponding to the land type patches and each equalization unit of the same land type in the equalization layer; By comparing the various characteristic element data corresponding to the land type patches and each equalization unit, determine the equalization unit with the closest characteristic elements to the land type patches, and assign the utilization grade information of the equalization unit to the land type patches.

Images