Agricultural land concentration calculation method and system considering peripheral connected elements

Abstract

The application discloses a kind of cultivated land concentrated piece degree calculation method and system considering peripheral connected elements, belong to cultivated land layout planning field.The method specifically includes: obtaining preprocessed cultivated land map spot vector data;Filter ditch, rural road, river and so on adjacent to cultivated land as connected element data;With cultivated land map spot as center, generate buffer zone, filter connected element and merge to form piece map spot vector data set;After fusing adjacent piece map spot, calculate the net area of cultivated land in each spot, and obtain piece degree result by hierarchical statistics.The method considers the correlation between cultivated land and peripheral land class, innovatively introduces peripheral connected element, solves the problem that traditional method is difficult to accurately analyze piece pattern under multi-scale, avoids grid data error through vector data modeling, and parameters can be adjusted according to regional characteristics, to provide scientific basis for cultivated land layout optimization.

Description

Technical Field

[0001] This invention belongs to the field of farmland layout planning, and in particular relates to a method and system for calculating the degree of contiguous farmland concentration that takes into account surrounding connectivity factors. Background Technology

[0002] The degree of contiguous farmland is an important and commonly used key indicator for measuring whether the quality and spatial distribution pattern of farmland are conducive to large-scale, mechanized farming and efficient management. It reflects the current state of farmland layout and the effectiveness of contiguous farmland management from a macro perspective, and is of great significance for ensuring national food security, improving agricultural production efficiency, reducing farming costs, and implementing effective farmland protection policies. Scientifically and effectively calculating the degree of contiguous farmland is a crucial foundation for optimizing farmland layout, promoting farmland protection, ensuring food security, and optimizing the spatial pattern of the national territory.

[0003] Existing research on methods for calculating the contiguousness of cultivated land both domestically and internationally mainly includes spatial connectivity calculation methods, morphological delineation methods, vector buffer methods, and landscape index methods. Early on, Zhou Shangyi et al. first proposed using spatial connectivity calculation methods to analyze farmland contiguousness. Guo Zihan and Yang Yongxia, based on this, designed a method combining vector and raster methods for calculating cultivated land contiguousness. Li Mengyang et al. proposed a multi-threshold spatial connectivity analysis method, combining landscape index analysis to identify the contiguous pattern of basic farmland. Wang Xiaoyan et al. used a method based on GIS and mathematical morphology for contiguousness analysis of concentrated cultivated land areas. Lu Xuejun et al. proposed a morphological closure algorithm, using multidisciplinary methods to explore calculation paths. Duan Gang proposed a vector-based buffer determination method for calculating contiguousness, which has been widely used in the contiguousness analysis of cultivated land vector plots.

[0004] However, although the aforementioned methods have been applied in many places, most of them are based on raster data or convert vector data into raster data before calculation, and are mostly used for studies at the county level and below. At larger scales, they are prone to data errors and redundancy, making them difficult to adapt to the precise and accurate analysis of contiguous farmland patterns across multiple scales, and unable to meet the needs of modern agriculture for refined management. Li Bo et al. designed a non-uniform grid contiguous clustering method that takes into account the differences in farmland area by optimizing the grid construction method, which addresses the problems of uneven data distribution and data redundancy, but data distortion still exists. Currently, the methods for calculating the concentration and contiguousness of farmland focus more on the farmland type itself, using simple buffer analysis and thresholds to determine contiguous areas, without fully considering the connectivity, attachment, and isolation relationships between farmland and surrounding related elements in agricultural scenarios, such as the irrigation role of ditches and water networks as agricultural ancillary facilities, and the obstruction effect of roads on large-scale contiguous farming by agricultural machinery. Furthermore, they cannot effectively distinguish and quantify neighboring elements of different types and different barriers, making it difficult to scientifically measure and reflect the spatial characteristics of farmland according to actual conditions.

[0005] Therefore, the problems with the current methods for calculating the contiguousness of arable land can be summarized in two aspects:

[0006] First, current common methods are mainly designed for raster data, which makes it difficult to adapt to multi-scale, especially large-scale, precise and effective analysis of concentrated and contiguous farmland patterns.

[0007] Second, ignoring the relationship between arable land and other surrounding elements, failing to view the agricultural system as a whole, focusing only on the single land type of arable land, and failing to consider in detail and quantitatively the connectivity and barrier effects of non-arable land elements adjacent to arable land patches on the spatial contiguousness of arable land may lead to measurement results that are inconsistent with the actual contiguous status and the needs of improvement.

[0008] In summary, there is a need to design a method for calculating the degree of contiguousness of cultivated land that can be adapted to various scales and takes into account the connectivity elements around cultivated land, so as to improve the accuracy of the analysis of the contiguous pattern of cultivated land. This is a technical problem that urgently needs to be solved. Summary of the Invention

[0009] The technical problem this invention aims to solve is that existing analyses of concentrated and contiguous farmland mainly rely on raster data, making it difficult to accurately analyze multi-scale, especially large-scale, farmland contiguous patterns. The invention addresses how to achieve accurate calculation of farmland concentration and contiguousness at multiple scales to reflect the overall farmland layout. Furthermore, existing technologies do not consider surrounding land types from a holistic perspective, resulting in a relatively singular analysis object that fails to objectively reflect the spatial distribution characteristics of farmland under actual farming scenarios and management needs. The invention seeks to effectively quantify the connectivity and segmentation of farmland space by various non-farmland elements surrounding farmland patches, thereby more realistically reflecting the degree of farmland concentration and contiguousness in actual management. The invention provides a method and system for calculating farmland concentration and contiguousness that considers surrounding connectivity elements.

[0010] To achieve the above-mentioned objectives, the present invention specifically adopts the following technical solution:

[0011] In a first aspect, the present invention provides a method for calculating the contiguousness of cultivated land considering surrounding connectivity factors, which includes the following steps:

[0012] S1. Obtain the crop patch vector dataset within the target area and preprocess it. The preprocessed crop patch vector dataset consists of multiple crop patches.

[0013] S2. For each type of adjacent land use, the adjacent land use types that have a connectivity effect on contiguous cultivation of arable land are taken as the data type of connectivity element, and a data type set of connectivity element data around arable land is formed; based on the data type of connectivity element, the vector data of arable land related land use patches in the target area is obtained and preprocessed.

[0014] S3. On the plane space where a cultivated land patch is located, take the cultivated land patch as the center and extend it outward by a certain distance to form a cultivated land vector buffer patch, and satisfy that the distance from all points in the cultivated land vector buffer patch to the cultivated land patch does not exceed the preset threshold for filtering distance of cultivated land surrounding elements; traverse all cultivated land patches, and form a cultivated land vector buffer zone from all the obtained cultivated land vector buffer patches.

[0015] S4. Construct spatial filtering rules. Using the cultivated land vector buffer as a reference, filter the connected feature data that intersects with the cultivated land vector buffer in the preprocessed cultivated land associated land type vector data, and construct the cultivated land surrounding connected feature dataset from the connected feature data.

[0016] S5. Spatially merge the preprocessed cultivated land patch vector dataset in S1 with the cultivated land surrounding connectivity element dataset in S4 to obtain a contiguous patch vector dataset containing complete cultivated land patches and cultivated land surrounding connectivity elements. Merge the spatially adjacent contiguous patches in the contiguous patch vector dataset and segment the merged contiguous patches according to the target area boundary to obtain a contiguous patch fusion vector dataset.

[0017] S6. For a contiguous patch in the contiguous patch fusion vector dataset, spatially connect the contiguous patch and the cultivated patch that intersects with it according to the spatial intersection relationship, calculate the total net area of ​​the cultivated patch covered within the range of the contiguous patch, and take the summation result as the net area of ​​contiguous cultivated land within the contiguous patch. Traverse the contiguous patch fusion vector dataset to obtain the net area of ​​contiguous cultivated land corresponding to all contiguous patches.

[0018] S7. Preset multiple area classification intervals, count the number of contiguous patches in each area classification interval and the sum of the net areas of all contiguous cultivated land in that area classification interval, so as to obtain the result of the concentration and contiguousness of cultivated land in the target area.

[0019] Based on the above scheme, each step can be implemented in the following preferred manner.

[0020] As a preferred embodiment of the first aspect mentioned above, in step S1, the cultivated land parcel vector dataset is obtained by filtering the land use status data in the national land survey data. The filtering categories include three secondary categories: paddy fields, irrigated land, and dry land. The filtering rule is the land use code value of the land use status data. When preprocessing the cultivated land parcel vector dataset, it includes geometric repair of topological errors.

[0021] As a preferred embodiment of the first aspect mentioned above, in step S2, the connected element data types include three types: ditch land type, rural road land type, and river land type. Using the land type code value as a filtering rule, three types of cultivated land related land type vector data are selected from the land use status data: ditch vector data belonging to the ditch land type, rural road vector data belonging to the rural road land type, and river water surface vector data belonging to the river land type. The cultivated land related land type vector data is then preprocessed, and topological errors are geometrically corrected. The preprocessed ditch vector data, preprocessed rural road vector data, and preprocessed river water surface vector data are used to form ditch vector dataset, rural road vector dataset, and river water surface vector dataset, respectively.

[0022] As a preferred embodiment of the first aspect mentioned above, the formation process of the cultivated land vector buffer zone in step S3 is as follows:

[0023]

[0024] in, A represents a vector buffer patch of cultivated land; A represents any point in the planar space where the cultivated land patch is located. This represents the distance from the cultivated patch to point A. This represents a cultivated patch in the cultivated patch vector dataset; This indicates the threshold distance for filtering elements around cultivated land.

[0025] As a preferred embodiment of the first aspect above, the specific process of obtaining the dataset of connected elements surrounding cultivated land by spatial filtering rules in step S4 is as follows:

[0026] S41. Perform a spatial overlay and intersection operation on the cultivated land vector buffer and the ditch patch vector dataset selected in S2 to obtain the set of ditch connectivity elements associated with the cultivated land vector buffer around the cultivated land; wherein, the set of ditch connectivity elements consists of ditch patches in the ditch patch vector dataset that are located within the cultivated land vector buffer.

[0027] S42. Perform a spatial overlay and intersection operation on the cultivated land vector buffer and the rural road patch vector dataset selected in S2 to obtain the set of rural road connectivity elements associated with the cultivated land vector buffer around the cultivated land; wherein, the rural road connectivity element set consists of rural road patches in the rural road patch vector dataset that are located within the cultivated land vector buffer.

[0028] S43. Perform a spatial overlay and intersection operation between the cultivated land vector buffer and the river surface vector dataset selected in S2 to obtain the set of river surface connectivity elements associated with the cultivated land vector buffer around the cultivated land; wherein, the set of river surface connectivity elements consists of river surface patches in the river surface vector dataset that are located within the cultivated land vector buffer.

[0029] S44. Take the ditch patches in the ditch connectivity feature set, the rural road patches in the rural road connectivity feature set, and the river water surface patches in the river water surface connectivity feature set as connectivity feature data, and make up the farmland surrounding connectivity feature dataset from all connectivity feature data.

[0030] As a preferred embodiment of the first aspect mentioned above, in step S5, contiguous patches are merged based on the principle of spatial adjacency. The specific fusion process is as follows: when the boundaries between two contiguous patches intersect, the two contiguous patches with intersecting boundaries are merged into a new contiguous patch. Each contiguous patch in the contiguous patch vector dataset is traversed until each contiguous patch in the contiguous patch vector dataset is spatially separated from other contiguous patches. The final contiguous patch is used as the fused contiguous patch.

[0031] As a preferred embodiment of the first aspect, the specific process of obtaining the net area of ​​contiguous cultivated land in step S6 is as follows: In the preprocessed cultivated land vector dataset, cultivated land patches that intersect with the contiguous cultivated land patch in space are taken as candidate cultivated land patches, the net area of ​​each candidate cultivated land patch corresponding to the contiguous cultivated land patch is obtained, and the net areas of all candidate cultivated land patches are added together to obtain the net area of ​​contiguous cultivated land within the contiguous cultivated land patch.

[0032] Furthermore, in step S6, for the first... Patchwork pattern The corresponding net area of ​​contiguous arable land Calculate using the following formula:

[0033]

[0034] in, The i-th contiguous patch intersects with the i-th patch space. The net area of ​​a cultivated land patch is the net area of ​​a candidate cultivated land patch; n is the number of cultivated land patches that intersect with the i-th contiguous patch space.

[0035] As a preferred embodiment of the first aspect mentioned above, in step S7, multiple consecutive area grading intervals are set, the number of contiguous patches in each area grading interval is counted sequentially, and within each area grading interval, the net area of ​​contiguous cultivated land corresponding to all contiguous patches in that interval is added together, and the sum of the number of contiguous patches and the net area of ​​contiguous cultivated land contained in each area grading interval is used to obtain the result of the concentration and contiguousness of cultivated land in the target area.

[0036] Secondly, the present invention provides a system for calculating the contiguousness of cultivated land considering surrounding connectivity factors, comprising:

[0037] The cultivated land patch acquisition module is used to acquire the cultivated land patch vector dataset within the target area and preprocess it. The preprocessed cultivated land patch vector dataset consists of multiple cultivated land patches.

[0038] The module for obtaining farmland-related land use parcels is used to, for each type of adjacent land use, take the adjacent land use types that have a connectivity effect on contiguous farmland cultivation as the data type of connectivity element, and form a set of data type of connectivity element around farmland; based on the data type of connectivity element, obtain the vector data of farmland-related land use parcels in the target area, and preprocess them.

[0039] The cultivated land vector buffer acquisition module is used to expand outward a certain distance from the cultivated land patch in the planar space where the cultivated land patch is located, and to satisfy that the distance from all points in the cultivated land vector buffer patch to the cultivated land patch does not exceed the preset threshold for filtering distance of cultivated land surrounding elements; it traverses all cultivated land patches and forms a cultivated land vector buffer from all the obtained cultivated land vector buffer patches.

[0040] The connectivity feature acquisition module is used to construct spatial filtering rules. Taking the cultivated land vector buffer as a reference, it filters the connectivity feature data that intersects with the cultivated land vector buffer in the preprocessed cultivated land associated land type vector data, and constructs the cultivated land surrounding connectivity feature dataset from the connectivity feature data.

[0041] The contiguous patch fusion module is used to spatially merge the preprocessed cultivated land patch vector dataset with the cultivated land surrounding connectivity element dataset to obtain a contiguous patch vector dataset containing complete cultivated land patches and cultivated land surrounding connectivity elements. It merges spatially adjacent contiguous patches in the contiguous patch vector dataset and segments the merged contiguous patches according to the target area boundary to obtain a contiguous patch fusion vector dataset.

[0042] The module for calculating the net area of ​​contiguous cultivated land is used to spatially connect a contiguous patch in the contiguous patch fusion vector dataset with the cultivated land patches that intersect with it based on spatial intersection relationships, calculate the total net area of ​​cultivated land patches covered within the contiguous patch, and use the summation result as the net area of ​​contiguous cultivated land within the contiguous patch. It then iterates through the contiguous patch fusion vector dataset to obtain the net area of ​​contiguous cultivated land corresponding to all contiguous patches.

[0043] The result acquisition module is used to preset multiple area classification intervals, count the number of contiguous patches in each area classification interval and the sum of the net areas of all contiguous cultivated land in that area classification interval, thereby obtaining the result of the concentration and contiguousness of cultivated land in the target area.

[0044] Thirdly, the present invention provides a computer program product, including a computer program / instruction, which, when executed by a processor, can implement the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements as described in any of the solutions of the first aspect above.

[0045] Fourthly, the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements as described in any of the solutions of the first aspect above.

[0046] Compared with the prior art, the present invention has the following advantages:

[0047] This invention discloses a method for calculating the contiguousness of cultivated land considering surrounding connectivity elements. Addressing the shortcomings of traditional cultivated land contiguousness analysis, such as its difficulty in accurately analyzing contiguous cultivated land patterns at multiple scales, especially large scales, and its failure to fully consider the connectivity, attachment, and barrier relationships between cultivated land and surrounding elements in actual agricultural scenarios, this invention designs a method for calculating the contiguousness of cultivated land based on vector data. This method can be applied at multiple scales, including regions, counties, cities, and provinces, avoiding an imbalance between data errors and redundancy. Furthermore, considering the topography and cultivation characteristics of contiguous cultivated land, this invention takes a systemic approach, not limiting itself to a single cultivated land type, but incorporating surrounding land types into the contiguousness calculation process. Based on the different interactions between cultivated land and surrounding land types, it selects attachment and connecting land types such as ditches, rural roads, and rivers that are adapted to modern agricultural production and management methods, while removing barrier land types such as villages, highways, and railways that hinder large-scale cultivation. This results in a new cultivated land contiguousness calculation rule that better reflects actual concentrated cultivation scenarios and reduces errors caused by land type ambiguity in the original raster format. This invention provides a novel solution for calculating the contiguousness of arable land, enabling more scientific and accurate analysis of contiguous arable land patterns that better aligns with actual conditions and remediation needs. This provides clear targets and scientific basis for precisely implementing arable land layout optimization. Furthermore, model parameters such as distance thresholds can be scientifically adjusted based on the geographical characteristics, farming habits, and management needs of different regions, demonstrating good adaptability and potential for wider application. Attached Figure Description

[0048] Figure 1 This is a flowchart of the method of the present invention;

[0049] Figure 2This is a flowchart illustrating step S4 of the present invention;

[0050] Figure 3 This is a schematic diagram of step S5 of the present invention;

[0051] Figure 4 This is a schematic diagram of spatial intersection according to the present invention;

[0052] Figure 5 This is a schematic diagram of spatial merging according to the present invention;

[0053] Figure 6 This is a schematic diagram of spatial fusion according to the present invention;

[0054] Figure 7 This is a schematic diagram of a portion of the crop patch vector data layer layerGD in an embodiment of the present invention;

[0055] Figure 8 This is a schematic diagram of a portion of the cultivated land vector buffer layer (layerBufferGD) according to an embodiment of the present invention.

[0056] Figure 9 This is a schematic diagram of the connectivity elements within a certain range around the cultivated land in an embodiment of the present invention;

[0057] Figure 10 This is a schematic diagram of the fused continuous patch pattern according to an embodiment of the present invention. Detailed Implementation

[0058] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Many specific details are set forth in the following description to provide a thorough understanding of the present invention. However, the present invention can be practiced in many other ways different from those described herein, and those skilled in the art can make similar modifications without departing from the spirit of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below. Technical features in the various embodiments of the present invention can be combined accordingly without mutual conflict.

[0059] In the description of this invention, it should be understood that the terms "first" and "second" are used only for descriptive purposes and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Therefore, a feature defined with "first" and "second" may explicitly or implicitly include at least one of those features.

[0060] like Figure 1 As shown, in a preferred embodiment of the present invention, the method for calculating the contiguousness of cultivated land considering surrounding connectivity factors includes the following steps S1 to S7. The specific implementation process of each step will be described in detail below.

[0061] S1. Obtain the crop patch vector dataset (G) within the target area and preprocess it. The preprocessed crop patch vector dataset consists of multiple crop patches.

[0062] It should be noted that in step S1 of this invention, the cultivated land parcel vector dataset is obtained by filtering the land use status data in the national land survey data. The filtering categories include three secondary categories: paddy fields, irrigated land, and dry land. The filtering rule is the land use code value of the land use status data. When preprocessing the cultivated land parcel vector dataset, it includes geometric repair of topological errors.

[0063] In step S1 of this embodiment, the cultivated land plot vector dataset is filtered based on the land use status data in the national land survey data. The filtering categories include three secondary categories: paddy fields, irrigated land, and dry land. The filtering rule is that the land type code is 0101, 0102, or 0103.

[0064] S2. For each type of adjacent land use, the adjacent land use types that have a connectivity effect on contiguous cultivation of arable land are taken as the data type of connectivity element, and a data type set of connectivity element data around arable land is formed; based on the data type of connectivity element, the vector data of arable land related land use patches in the target area is obtained and preprocessed.

[0065] It should be noted that in step S2 of this invention, the connected element data types include three types: ditch land type, rural road land type, and river land type. Then, using the land type code value as the filtering rule, three types of cultivated land related land type vector data are selected from the land use status data: ditch vector data belonging to the ditch land type, rural road vector data belonging to the rural road land type, and river water surface vector data belonging to the river land type. The cultivated land related land type vector data is then preprocessed, and topological errors are geometrically corrected. The preprocessed ditch vector data, preprocessed rural road vector data, and preprocessed river water surface vector data are used to form ditch vector dataset, rural road vector dataset, and river water surface vector dataset, respectively.

[0066] In step S2 of this embodiment, for each type of adjacent land use, its potential connectivity and separation / barrier effects on contiguous farmland cultivation are analyzed. For example, contiguous farmland has characteristics such as mechanized operations, good drainage, and convenient irrigation, while construction land such as villages, roads, and railways have a barrier effect on farmland. Therefore, surrounding factors are not considered for this type of land use. After the above analysis, three types of farmland-related land use categories—ditch land, rural road land, and river land—are selected from all adjacent land use categories to construct a data set of farmland surrounding connectivity elements.

[0067] After forming a data set of connected elements around cultivated land, the land use status data is then filtered based on the connected element data type. Data belonging to the ditch land category is selected to form a ditch patch vector dataset, data belonging to the rural road land category is selected to form a rural road patch vector dataset, and data belonging to the river land category is selected to form a river water surface patch vector dataset.

[0068] Specifically, in step S2 of this embodiment, the above screening process is as follows:

[0069] S21. From the land use status data, filter out the ditch patch vector data according to the land use code value 1107 or 1107A to form the ditch patch vector dataset (Q).

[0070] S22. From the land use status data, select rural road patch vector data according to the land category code value 1006, thus forming a rural road patch vector dataset (N).

[0071] S23. From the land use status data, select the river water surface vector data according to the land category code value 1101 to form the river water surface vector dataset (H).

[0072] It should also be noted that the core of this invention is to consider the calculation of the concentration and contiguousness of cultivated land based on the connectivity elements surrounding cultivated land. Due to differences in farming traditions and topography, the connectivity elements surrounding cultivated land have different functions in different regions. For example, in the northern plains, rivers have relatively wide water surfaces and mainly serve the functions of navigation and transportation, which can have a blocking effect on contiguous cultivation. In the southern regions, the water network is dense and mainly serves the function of irrigation, which can assist agricultural production. Therefore, the applicability of rivers as connectivity elements can be appropriately considered in different regions, and the current screening rules for river elements around cultivated land can be replaced by setting different buffer zones.

[0073] S3. On the plane space where a cultivated land patch is located, take the cultivated land patch as the center and extend it outward by a certain distance to form a cultivated land vector buffer patch, and satisfy that the distance from all points in the cultivated land vector buffer patch to the cultivated land patch does not exceed the preset threshold for filtering distance of cultivated land surrounding elements; traverse all cultivated land patches, and form a cultivated land vector buffer zone from all the obtained cultivated land vector buffer patches.

[0074] It should be noted that in step S3 of this invention, a screening distance threshold for surrounding elements of cultivated land is set to quantify the actual effect of surrounding land types on cultivated land. If the distance is too long, the connectivity effect of the surrounding element is considered weaker than its blocking effect. A buffer analysis is then performed on the cultivated land patches based on this distance threshold to obtain the cultivated land vector buffer (BG). The cultivated land vector buffer consists of all cultivated land vector buffer patches (bg) that meet the conditions. For any point A within a cultivated land vector buffer patch, its distance d(g, A) to a cultivated land patch (g) must be less than or equal to the screening distance threshold for surrounding elements of cultivated land. If the distance exceeds the screening distance threshold, the point is considered not to be within the cultivated land vector buffer patch, and the corresponding surrounding element has a weak connectivity effect on the cultivated land, or may even act as a barrier.

[0075] In this embodiment, the formation process of the above-mentioned cultivated land vector buffer zone is as follows:

[0076]

[0077] in, A represents a vector buffer patch of cultivated land; A represents any point in the planar space where the cultivated land patch is located. This represents the distance from the cultivated patch to point A. This represents a cultivated patch in the cultivated patch vector dataset; This represents the threshold for filtering elements around cultivated land; in this embodiment, it is taken as... .

[0078] S4. Construct spatial filtering rules. Using the cultivated land vector buffer as a reference, filter the connected feature data that intersects with the cultivated land vector buffer in the preprocessed cultivated land associated land parcel vector data, and construct the cultivated land surrounding connected feature dataset from the connected feature data.

[0079] It should be noted that in step S4 of the present invention, spatial filtering rules are constructed, and connected element data within a certain range around the cultivated land that intersects with the cultivated land vector buffer are filtered out. For example, ditch patches, rural road patches and river surface patches that intersect with the cultivated land vector buffer are filtered out, thereby forming a set of ditch connected elements, a set of rural road connected elements and a set of river surface connected elements.

[0080] In step S4 of this embodiment, as Figure 2 As shown, the specific process of obtaining the connected feature dataset around cultivated land using spatial filtering rules is as follows:

[0081] S41. Filtering ditch patches around cultivated land: Perform spatial overlay and intersection operation on the cultivated land vector buffer and the ditch patch vector dataset filtered in S2 to obtain the ditch connectivity feature set (BQ) associated with the cultivated land vector buffer; wherein, the ditch connectivity feature set consists of ditch patches in the ditch patch vector dataset that are located within the cultivated land vector buffer.

[0082] In step S41 of this embodiment, the calculation process for obtaining the ditch connectivity element set is as follows:

[0083]

[0084] In the formula: This indicates finding the intersection.

[0085] S42. Filtering rural road patches around cultivated land: Perform spatial overlay and intersection operation on the cultivated land vector buffer and the rural road patch vector dataset filtered in S2 to obtain the rural road connectivity feature set (BN) associated with the cultivated land vector buffer; wherein, the rural road connectivity feature set consists of rural road patches in the rural road patch vector dataset that are located within the cultivated land vector buffer.

[0086] In step S42 of this embodiment, the calculation process for obtaining the rural road connectivity element set is as follows:

[0087]

[0088] S43. Filtering river surface patches around cultivated land: Perform spatial overlay and intersection operation on the cultivated land vector buffer and the river surface patch vector dataset filtered in S2 to obtain the river surface connectivity feature set (BH) associated with the cultivated land vector buffer; wherein, the river surface connectivity feature set consists of river surface patches in the river surface patch vector dataset that are located within the cultivated land vector buffer.

[0089] In step S43 of this embodiment, the calculation process for obtaining the river surface connectivity element set is as follows:

[0090]

[0091] S44. Take the ditch patches in the ditch connectivity feature set, the rural road patches in the rural road connectivity feature set, and the river surface patches in the river surface connectivity feature set as connectivity feature data, and form a farmland periphery connectivity feature dataset (BU) from all connectivity feature data.

[0092] In step S44 of this embodiment, the ditch connectivity feature set, the rural road connectivity feature set, and the river surface connectivity feature set actually represent the ditch, rural road, and river surface features connected to the cultivated land. These three types of features are also connectivity feature data. Therefore, this embodiment spatially merges the ditch connectivity feature set, the rural road connectivity feature set, and the river surface connectivity feature set, that is, takes the union of the three, thereby forming the cultivated land surrounding connectivity feature dataset. The process of taking the intersection in S41~S43 is as follows: Figure 4 As shown, the process of taking the union of S44 is as follows: Figure 5 As shown.

[0093] S5. Spatially merge the preprocessed cultivated land parcel vector dataset in S1 with the cultivated land surrounding connectivity element dataset (BU) in S4 to obtain a contiguous parcel vector dataset (U) containing complete cultivated land parcels and cultivated land surrounding connectivity elements. Merge the spatially adjacent contiguous parcels in the contiguous parcel vector dataset, and segment the merged contiguous parcels according to the target area boundary to obtain a contiguous parcel fusion vector dataset (D).

[0094] It should be noted that in step S5 of the present invention, as Figure 3 As shown, based on the connectivity features of cultivated land patches themselves and their connectivity features with surrounding connected elements, contiguous patches are aggregated to form a contiguous patch vector dataset (U). Then, the contiguous patches in this dataset are fused to obtain a fused contiguous patch vector dataset (D). The calculation process for obtaining the contiguous patch vector dataset is as follows:

[0095]

[0096] It should be noted that in step S5 of the present invention, the contiguous patches are fused based on the principle of spatial adjacency. The specific fusion process is as follows: when the boundaries between two contiguous patches intersect, the two contiguous patches with intersecting boundaries are merged into a new contiguous patch. Every contiguous patch in the contiguous patch vector dataset is traversed until every contiguous patch in the contiguous patch vector dataset is spatially separated from other contiguous patches. The final contiguous patch is the fused contiguous patch.

[0097] In this embodiment, for a contiguous patch in the contiguous patch vector dataset, if there are other contiguous patches in the dataset that intersect with this contiguous patch at their boundaries, the contiguous patch and the intersecting contiguous patches are merged into a new contiguous patch. This new contiguous patch replaces the two contiguous patches before merging and is added to the contiguous patch vector dataset. The dataset continues to determine if there are any contiguous patches intersecting with its boundaries. If so, the intersecting contiguous patches are merged with the new contiguous patch. If not, the next contiguous patch in the contiguous patch vector dataset is used as the basis for the search, and the above process is repeated. This continues until every contiguous patch in the contiguous patch vector dataset no longer has a spatial adjacency relationship with any other contiguous patch. The final contiguous patch is the fused contiguous patch.

[0098] S6. For a contiguous patch in the contiguous patch fusion vector dataset, spatially connect the contiguous patch and the cultivated land patches that intersect with it according to the spatial intersection relationship, calculate the total net area of ​​the cultivated land patches covered within the contiguous patch, and take the summation result as the net area of ​​contiguous cultivated land within the contiguous patch. Traverse the contiguous patch fusion vector dataset to obtain the net area of ​​contiguous cultivated land corresponding to all contiguous patches, thereby obtaining the scale of contiguous cultivated land for each contiguous patch.

[0099] It should be noted that, in step S6 of the present invention, the specific process of obtaining the net area of ​​contiguous cultivated land in a contiguous patch is as follows: In the preprocessed cultivated patch vector dataset, cultivated patches that intersect with the contiguous patch space are taken as candidate cultivated patches, the net area of ​​each candidate cultivated patch corresponding to the contiguous patch is obtained, and the net areas of all candidate cultivated patches are added together to obtain the net area of ​​contiguous cultivated land in the contiguous patch, which is taken as the scale of the contiguous cultivated land in the contiguous patch.

[0100] In step S6 of this embodiment, based on the contiguous patch fusion vector dataset, it is spatially connected with cultivated patches according to the spatial intersection relationship, and the net area of ​​cultivated patches covered within the range of each contiguous patch is calculated to accurately obtain the cultivated land area, which is used as the cultivated land scale of the contiguous patch.

[0101] For the Patchwork pattern The corresponding net area of ​​contiguous arable land Calculate using the following formula:

[0102]

[0103] in, The i-th contiguous patch intersects with the i-th patch space. The net area of ​​a cultivated land patch is the net area of ​​a candidate cultivated land patch; n is the number of cultivated land patches that intersect with the i-th contiguous patch space.

[0104] S7. Preset multiple area classification intervals, count the number of contiguous patches in each area classification interval and the sum of the net areas of all contiguous cultivated land in that area classification interval, so as to obtain the result of the concentration and contiguousness of cultivated land in the target area.

[0105] It should be noted that in step S7 of the present invention, multiple consecutive area grading intervals are set, the number of contiguous patches in each area grading interval is counted sequentially, and within each area grading interval, the net area of ​​contiguous cultivated land corresponding to all contiguous patches in that interval is added together, and the sum of the number of contiguous patches and the net area of ​​contiguous cultivated land contained in each area grading interval is used as the result of the concentration and contiguousness of cultivated land in the target area.

[0106] To better demonstrate the specific implementation and technical effects of the present invention, the method for calculating the concentration and contiguousness of cultivated land considering surrounding connectivity elements, as shown in steps S1 to S7 of the above preferred implementation, is used to conduct a calculation and analysis of the concentration and contiguousness of cultivated land in a certain area.

[0107] Example

[0108] The specific implementation process of the method for calculating the concentration and contiguousness of cultivated land considering surrounding connectivity elements used in this embodiment is as described above and will not be repeated here. In this embodiment, the method proposed in this invention is used below, in conjunction with the accompanying drawings, to calculate and analyze the concentration and contiguousness of cultivated land in a certain area.

[0109] First, this embodiment uses a land use map layer based on land use status data within a certain province. Cultivated land parcel vector data is selected according to the rule that the attribute field "Land Use Code" value is "0101", "0102", or "0103". This includes three subcategories: paddy fields, irrigated land, and dry land, which is layerGD. Figure 7 As shown. A routine topology check is performed on the layerGD layer, including topology issues such as self-intersection of polygons and overlapping features. Geometric repairs are then performed on topology errors to ensure the topological correctness of the polygon vectors.

[0110] Next, the vector data of cultivated land associated with different land use categories are acquired, including land use categories such as ditches, rural roads, and river surfaces. This vector data is then preprocessed. Specifically, based on the land use category layer, ditch vector data is selected according to the rule that the "Land Use Code" value is "1107" or "1107A", and named layerGQ; rural road vector data is selected according to the rule that the "Land Use Code" value is "1006", and named layerNCDL; and river surface vector data is selected according to the rule that the "Land Use Code" value is "1101", and named layerHL. A routine topology check is performed on these layers, including checking for topological issues such as self-intersection of land use categories and overlapping features. Geometric corrections are then performed to ensure the topological correctness of the land use vectors.

[0111] Next, buffer analysis is performed on the layerGD layer with a certain distance as the radius to generate a farmland vector buffer. In this embodiment, a buffer radius of 10 meters is selected. The farmland vector buffer patches are then merged according to the rule of spatial adjacency; that is, spatially adjacent patches are merged into one large patch to avoid overlap between buffer patches. The newly generated merged layer is named layerBufferGD, as shown below. Figure 8 As shown.

[0112] Then, related features within a certain range around the cultivated land are obtained. Specifically, spatial intersection analysis is performed on layerBufferGD and layerGQ to generate ditch vector data within the cultivated land vector buffer range, named layerInGQ; spatial intersection analysis is performed on layerBufferGD and layerNCDL to generate rural road vector data within the cultivated land vector buffer range, named layerInNCDL; and spatial intersection analysis is performed on layerBufferGD and layerHL to generate river vector data within the cultivated land vector buffer range, named layerInHL. The related feature patches filtered through the cultivated land vector buffer are shown below. Figure 9 As shown. Since this method focuses on the connectivity between cultivated lands, the related features around the cultivated land in this embodiment are processed as a whole. Spatial merging is performed on layerInGQ, layerInNCDL, and layerInHL, concentrating them into one layer to obtain the related feature patches within a 10-meter radius around the cultivated land, named layerGLYS.

[0113] By spatially merging layerGD and layerGLYS, we obtain the vector data of cultivated land patches and related feature patches within a 10-meter radius around them, which is named layerBGD.

[0114] The image patches in layerBGD are merged according to the rule of spatial adjacency. For details of the merging process, please refer to [reference needed]. Figure 6 In other words, if a patch and its surrounding patches have intersecting boundaries, they are merged into a larger patch. This larger patch is then used as a base to search for other patches with intersecting boundaries, and this process is repeated until the patch and its surrounding patches are spatially separated. The resulting spatially merged patches are non-overlapping, non-overlapping, and form independent blocks. Finally, the merged patches undergo topology checks and repairs to avoid topology problems such as gaps and self-intersections.

[0115] Considering that actual farmland management projects are typically managed at the county level or below, this embodiment trims the spatially fused vector data according to the county-level administrative boundaries within the region, ultimately generating contiguous patch fused vector data, named layerLP, and adding a new field "Net Farmland Area" to record the net farmland area. The fused contiguous patches are shown below. Figure 10 As shown.

[0116] Spatially connect layerLP and layerGD according to the rules of spatial intersection. For each contiguous patch in layerLP, sum the net area of ​​the cultivated patches in layerGD that are connected to it, and record the total net area in the "Net Cultivated Land Area" field to ensure that the actual cultivated land area in the contiguous cultivated land patch is consistent with the original cultivated land patch area.

[0117] Based on different area classifications, information on contiguous cultivated land parcels of different sizes is summarized and statistically analyzed. This embodiment employs a four-level classification and aggregation, statistically analyzing the area and number of contiguous parcels within the following intervals: net cultivated land area greater than or equal to 500 mu and less than 1000 mu; greater than or equal to 1000 mu and less than 5000 mu; greater than or equal to 5000 mu and less than 10000 mu; and greater than or equal to 10000 mu. This data serves as the result for calculating the concentration and contiguousness of cultivated land in that region. Furthermore, this embodiment stipulates that parcels with a net cultivated land area of ​​less than 500 mu are not included in the statistical scope due to their small contiguous size.

[0118] This embodiment employs the technical methods proposed in this invention to scientifically and effectively measure and analyze the concentration and contiguousness of cultivated land within a sample area over three consecutive years. It effectively identifies areas within the sample area with a relatively good pattern of concentrated and contiguous cultivated land and areas with relatively fragmented cultivated land distribution, making the analysis results more consistent with the actual situation. Furthermore, through the annual change analysis of cultivated land concentration and contiguousness, it reflects the evolution of cultivated land patterns and the effectiveness of cultivated land protection and comprehensive land consolidation efforts in recent years. It also identifies shortcomings in cultivated land layout under the current management status, providing important data results and technical methodological support for cultivated land protection and consolidation and food security.

[0119] It should also be noted that the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements in the above embodiments can essentially be executed by a computer program or module. Therefore, similarly, based on the same inventive concept, another preferred embodiment of the present invention also provides a system for calculating the contiguousness of cultivated land considering surrounding connectivity elements, corresponding to the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements provided in the above embodiments, comprising:

[0120] The cultivated land patch acquisition module is used to acquire the cultivated land patch vector dataset within the target area and preprocess it. The preprocessed cultivated land patch vector dataset consists of multiple cultivated land patches.

[0121] The module for obtaining farmland-related land use parcels is used to, for each type of adjacent land use, take the adjacent land use types that have a connectivity effect on contiguous farmland cultivation as the data type of connectivity element, and form a set of data type of connectivity element around farmland; based on the data type of connectivity element, obtain the vector data of farmland-related land use parcels in the target area, and preprocess them.

[0122] The cultivated land vector buffer acquisition module is used to expand outward a certain distance from the cultivated land patch in the planar space where the cultivated land patch is located, and to satisfy that the distance from all points in the cultivated land vector buffer patch to the cultivated land patch does not exceed the preset threshold for filtering distance of cultivated land surrounding elements; it traverses all cultivated land patches and forms a cultivated land vector buffer from all the obtained cultivated land vector buffer patches.

[0123] The connectivity feature acquisition module is used to construct spatial filtering rules. Taking the cultivated land vector buffer as a reference, it filters the connectivity feature data that intersects with the cultivated land vector buffer in the preprocessed cultivated land associated land type vector data, and constructs the cultivated land surrounding connectivity feature dataset from the connectivity feature data.

[0124] The contiguous patch fusion module is used to spatially merge the preprocessed cultivated land patch vector dataset with the cultivated land surrounding connectivity element dataset to obtain a contiguous patch vector dataset containing complete cultivated land patches and cultivated land surrounding connectivity elements. It merges spatially adjacent contiguous patches in the contiguous patch vector dataset and segments the merged contiguous patches according to the target area boundary to obtain a contiguous patch fusion vector dataset.

[0125] The module for calculating the net area of ​​contiguous cultivated land is used to spatially connect a contiguous patch in the contiguous patch fusion vector dataset with the cultivated land patches that intersect with it based on spatial intersection relationships, calculate the total net area of ​​cultivated land patches covered within the contiguous patch, and use the summation result as the net area of ​​contiguous cultivated land within the contiguous patch. It then iterates through the contiguous patch fusion vector dataset to obtain the net area of ​​contiguous cultivated land corresponding to all contiguous patches.

[0126] The result acquisition module is used to preset multiple area classification intervals, count the number of contiguous patches in each area classification interval and the sum of the net areas of all contiguous cultivated land in that area classification interval, thereby obtaining the result of the concentration and contiguousness of cultivated land in the target area.

[0127] It is understood that the method for calculating the contiguousness of cultivated land considering surrounding connectivity factors described in S1-S7 above can essentially be implemented by a computer program. Therefore, based on the same inventive concept, another preferred embodiment of the present invention also provides a computer program product corresponding to the method for calculating the contiguousness of cultivated land considering surrounding connectivity factors provided in the above embodiments. This product includes a computer program / instruction that, when executed by a processor, can implement the method for calculating the contiguousness of cultivated land considering surrounding connectivity factors as described in the above embodiments.

[0128] Therefore, based on the same inventive concept, another preferred embodiment of the present invention also provides a computer-readable storage medium corresponding to the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements provided in the above embodiments. The storage medium stores a computer program, which, when executed by a processor, can realize the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements in the above embodiments.

[0129] It is understood that the aforementioned storage media may include random access memory (RAM) or non-volatile memory (NVM), such as at least one disk storage device. Furthermore, the storage media may also be various media capable of storing program code, such as USB flash drives, external hard drives, magnetic disks, or optical discs.

[0130] It should also be noted that those skilled in the art will understand that, for the sake of convenience and brevity, the specific working process of the system described above can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here. In the embodiments provided in this application, the division of steps or modules in the system and method is merely a logical functional division, and there may be other division methods in actual implementation. For example, multiple modules or steps may be combined or integrated together, and a module or step may also be split.

[0131] The embodiments described above are merely preferred embodiments of the present invention and are not intended to limit the invention. Those skilled in the art can make various changes and modifications without departing from the spirit and scope of the invention. Therefore, all technical solutions obtained through equivalent substitution or transformation fall within the protection scope of the present invention.

Claims

1. A method for calculating the contiguousness of cultivated land considering surrounding connectivity factors, characterized in that, Includes the following steps: S1. Obtain the crop patch vector dataset within the target area and preprocess it. The preprocessed crop patch vector dataset consists of multiple crop patches. S2. For each type of adjacent land use, the adjacent land use types that have a connectivity effect on contiguous cultivation of arable land are taken as the data type of connectivity element, and a data type set of connectivity element data around arable land is formed; based on the data type of connectivity element, the vector data of arable land related land use patches in the target area is obtained and preprocessed. S3. On the plane space where a cultivated land patch is located, take the cultivated land patch as the center and extend it outward by a certain distance to form a cultivated land vector buffer patch, and satisfy that the distance from all points in the cultivated land vector buffer patch to the cultivated land patch does not exceed the preset threshold for filtering distance of cultivated land surrounding elements; traverse all cultivated land patches, and form a cultivated land vector buffer zone from all the obtained cultivated land vector buffer patches. S4. Construct spatial filtering rules. Using the cultivated land vector buffer as a reference, filter the connected feature data that intersects with the cultivated land vector buffer in the preprocessed cultivated land associated land type vector data, and construct the cultivated land surrounding connected feature dataset from the connected feature data. S5. Spatially merge the preprocessed cultivated land patch vector dataset in S1 with the cultivated land surrounding connectivity element dataset in S4 to obtain a contiguous patch vector dataset containing complete cultivated land patches and cultivated land surrounding connectivity elements. Merge the spatially adjacent contiguous patches in the contiguous patch vector dataset and segment the merged contiguous patches according to the target area boundary to obtain a contiguous patch fusion vector dataset. S6. For a contiguous patch in the contiguous patch fusion vector dataset, spatially connect the contiguous patch and the cultivated patch that intersects with it according to the spatial intersection relationship, calculate the total net area of ​​the cultivated patch covered within the range of the contiguous patch, and take the summation result as the net area of ​​contiguous cultivated land within the contiguous patch. Traverse the contiguous patch fusion vector dataset to obtain the net area of ​​contiguous cultivated land corresponding to all contiguous patches. S7. Preset multiple area classification intervals, count the number of contiguous patches in each area classification interval and the sum of the net areas of all contiguous cultivated land in that area classification interval, so as to obtain the result of the concentration and contiguousness of cultivated land in the target area.

2. The method for calculating the contiguousness of cultivated land considering surrounding connectivity factors as described in claim 1, characterized in that, In step S1, the cultivated land parcel vector dataset is obtained by filtering the land use status data in the national land survey data. The filtering categories include three secondary categories: paddy fields, irrigated land, and dry land. The filtering rule is the land use code value of the land use status data. When preprocessing the cultivated land parcel vector dataset, geometric repair of topological errors is performed.

3. The method for calculating the contiguousness of cultivated land considering surrounding connectivity factors as described in claim 1, characterized in that, In step S2, the connected feature data types include three types: ditch land type, rural road land type, and river land type. Using the land type code value as the filtering rule, three types of cultivated land related land type vector data are selected from the land use status data: ditch vector data belonging to the ditch land type, rural road vector data belonging to the rural road land type, and river water surface vector data belonging to the river land type. The cultivated land related land type vector data is then preprocessed, and topological errors are geometrically corrected. The preprocessed ditch vector data, preprocessed rural road vector data, and preprocessed river water surface vector data are used to form ditch vector dataset, rural road vector dataset, and river water surface vector dataset, respectively.

4. The method for calculating the contiguousness of cultivated land considering surrounding connectivity factors as described in claim 3, characterized in that, In step S4, the specific process of obtaining the connected feature dataset around cultivated land using spatial filtering rules is as follows: S41. Perform a spatial overlay and intersection operation on the cultivated land vector buffer and the ditch patch vector dataset selected in S2 to obtain the set of ditch connectivity elements associated with the cultivated land vector buffer around the cultivated land; wherein, the set of ditch connectivity elements consists of ditch patches in the ditch patch vector dataset that are located within the cultivated land vector buffer. S42. Perform a spatial overlay and intersection operation on the cultivated land vector buffer and the rural road patch vector dataset selected in S2 to obtain the set of rural road connectivity elements associated with the cultivated land vector buffer around the cultivated land; wherein, the rural road connectivity element set consists of rural road patches in the rural road patch vector dataset that are located within the cultivated land vector buffer. S43. Perform a spatial overlay and intersection operation between the cultivated land vector buffer and the river surface vector dataset selected in S2 to obtain the set of river surface connectivity elements associated with the cultivated land vector buffer around the cultivated land; wherein, the set of river surface connectivity elements consists of river surface patches in the river surface vector dataset that are located within the cultivated land vector buffer. S44. Take the ditch patches in the ditch connectivity feature set, the rural road patches in the rural road connectivity feature set, and the river water surface patches in the river water surface connectivity feature set as connectivity feature data, and make up the farmland surrounding connectivity feature dataset from all connectivity feature data.

5. The method for calculating the contiguousness of cultivated land considering surrounding connectivity factors as described in claim 1, characterized in that, In step S5, contiguous patches are merged based on the principle of spatial adjacency. The specific fusion process is as follows: when the boundaries of two contiguous patches intersect, the two contiguous patches with intersecting boundaries are merged into a new contiguous patch. Every contiguous patch in the contiguous patch vector dataset is traversed until every contiguous patch in the contiguous patch vector dataset is spatially separated from other contiguous patches. The final contiguous patch is the fused contiguous patch.

6. The method for calculating the contiguousness of cultivated land considering surrounding connectivity factors as described in claim 1, characterized in that, In step S6, the specific process of obtaining the net area of ​​contiguous cultivated land in a contiguous patch is as follows: In the preprocessed cultivated patch vector dataset, cultivated patches that intersect with the contiguous patch space are taken as candidate cultivated patches, the net area of ​​each candidate cultivated patch corresponding to the contiguous patch is obtained, and the net areas of all candidate cultivated patches are added together as the net area of ​​contiguous cultivated land in the contiguous patch.

7. The method for calculating the contiguousness of cultivated land considering surrounding connectivity factors as described in claim 1, characterized in that, In step S7, multiple consecutive area grading intervals are set, the number of contiguous patches in each area grading interval is counted sequentially, and the net area of ​​contiguous cultivated land corresponding to all contiguous patches in each area grading interval is added together. The sum of the number of contiguous patches and the net area of ​​contiguous cultivated land in each area grading interval is used to obtain the result of the concentration and contiguousness of cultivated land in the target area.

8. A system for calculating the contiguousness of cultivated land considering surrounding connectivity factors, characterized in that, include: The cultivated land patch acquisition module is used to acquire the cultivated land patch vector dataset within the target area and preprocess it. The preprocessed cultivated land patch vector dataset consists of multiple cultivated land patches. The module for obtaining farmland-related land use parcels is used to, for each type of adjacent land use, take the adjacent land use types that have a connectivity effect on contiguous farmland cultivation as the data type of connectivity element, and form a set of data type of connectivity element around farmland; based on the data type of connectivity element, obtain the vector data of farmland-related land use parcels in the target area, and preprocess them. The cultivated land vector buffer acquisition module is used to expand outward a certain distance from the cultivated land patch in the planar space where the cultivated land patch is located, and to satisfy that the distance from all points in the cultivated land vector buffer patch to the cultivated land patch does not exceed the preset threshold for filtering distance of cultivated land surrounding elements; it traverses all cultivated land patches and forms a cultivated land vector buffer from all the obtained cultivated land vector buffer patches. The connectivity feature acquisition module is used to construct spatial filtering rules. Taking the cultivated land vector buffer as a reference, it filters the connectivity feature data that intersects with the cultivated land vector buffer in the preprocessed cultivated land associated land type vector data, and constructs the cultivated land surrounding connectivity feature dataset from the connectivity feature data. The contiguous patch fusion module is used to spatially merge the preprocessed cultivated land patch vector dataset with the cultivated land surrounding connectivity element dataset to obtain a contiguous patch vector dataset containing complete cultivated land patches and cultivated land surrounding connectivity elements. It merges spatially adjacent contiguous patches in the contiguous patch vector dataset and segments the merged contiguous patches according to the target area boundary to obtain a contiguous patch fusion vector dataset. The module for calculating the net area of ​​contiguous cultivated land is used to spatially connect a contiguous patch in the contiguous patch fusion vector dataset with the cultivated land patches that intersect with it based on spatial intersection relationships, calculate the total net area of ​​cultivated land patches covered within the contiguous patch, and use the summation result as the net area of ​​contiguous cultivated land within the contiguous patch. It then iterates through the contiguous patch fusion vector dataset to obtain the net area of ​​contiguous cultivated land corresponding to all contiguous patches. The result acquisition module is used to preset multiple area classification intervals, count the number of contiguous patches in each area classification interval and the sum of the net areas of all contiguous cultivated land in that area classification interval, thereby obtaining the result of the concentration and contiguousness of cultivated land in the target area.

9. A computer program product comprising a computer program / instructions, characterized in that, When the computer program / instruction is executed by the processor, it can implement the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The storage medium stores a computer program, which, when executed by a processor, implements the method for calculating the contiguousness of cultivated land considering surrounding connectivity elements as described in any one of claims 1 to 7.

Images