Method for realizing cultivated land subsidy issuing, electronic equipment and storage medium
Through the combination of multi-scale segmentation model and target classification model, subsidized plots in cultivated land are accurately identified and divided, which solves the problem of unclear plot edges in the existing technology, and improves the accuracy and effectiveness of farmland subsidies.
Patent Information
- Application Number
- CN202411923164.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-12-25
- Publication Date
- 2025-04-18
AI Technical Summary
In the monitoring of cultivated land, it is difficult for the prior art to accurately identify and divide land that is doped with small areas of other crops or land types, resulting in the effective issuance of cultivated land subsidies.
The remote sensing image is divided using a multi-scale segmentation model, the plot attribute features are extracted, the target feature vector is generated, and the target feature vector is input into the pre-acquisitioned target classification model, the geographical area type is determined, and the area of the subsidized plot is further divided into area grids, and the area of the subsidized plot is accurately identified.
The estimation accuracy of the subsidized land area has been improved, the effective distribution of cultivated land subsidies has been ensured, and the sustainable development of cultivated land has been promoted.
Smart Images
Figure CN119851153A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the field of image processing technology, and in particular to a method, electronic equipment and storage medium for implementing the distribution of cultivated land subsidies. Background Art
[0002] With the increasing emphasis on food security, the monitoring of non-agricultural and non-grain land is increasingly valued. Carrying out non-agricultural and non-grain land monitoring is crucial to the effective issuance of farmland subsidies, and can also promote and ensure the sustainable development of farmland. The subsidized land is mainly grain fields, which refer to relatively concentrated and contiguous plots based on permanent basic farmland, with good water and soil resources, relatively complete infrastructure, and crops such as wheat, corn, and rice.
[0003] At present, remote sensing technology is usually used to monitor cultivated land areas. High-resolution satellites are launched into space to obtain qualified images, and cultivated land plot information is identified through remote sensing. When segmenting image data, most plots can be segmented more accurately because the crops in subsidized plots are relatively concentrated and continuous plots, and different crops are significantly different. However, some plots are mixed with small areas of other crops or other plot types, resulting in unclear plot edges. Therefore, when segmenting plots, they cannot be accurately identified and segmented, affecting the effective issuance of cultivated land subsidies. Summary of the invention
[0004] In response to the above technical problems, the present invention provides a method, electronic device and storage medium for implementing the issuance of cultivated land subsidies, which can improve the accuracy of estimating the area of subsidized land, thereby realizing the effective issuance of cultivated land subsidies.
[0005] According to a first aspect of the present invention, there is provided a method for implementing the payment of cultivated land subsidies, comprising the following steps:
[0006] Based on several remote sensing images of the target geographical area acquired within a preset time period, a multi-scale segmentation model is used to divide each remote sensing image of the target geographical area into several sub-geographical area images to obtain several sub-geographical area images corresponding to each sub-geographical area.
[0007] For a plurality of sub-geographical area images corresponding to each sub-geographical area, plot attribute features are extracted from each sub-geographical area image to obtain a target feature vector corresponding to each sub-geographical area.
[0008] For any sub-geographical area, the target feature vector of the sub-geographical area is input into a pre-acquired target classification model, and the geographical area type of the sub-geographical area is determined according to the output result of the target classification model; the geographical area types include subsidizable plots, non-subsidizable plots and pending plots.
[0009] According to the geographic area type of each sub-geographic area, several sub-geographic areas corresponding to the subsidized plots are screened out, and each screened sub-geographic area is divided into several regional grids.
[0010] The target feature vector corresponding to each regional grid is obtained, and the target feature vector corresponding to each regional grid is input into the target classification model to determine the geographical area type of each regional grid, so as to realize the distribution of cultivated land subsidies according to the geographical area type of each regional grid.
[0011] According to a second aspect of the present invention, a non-transitory computer-readable storage medium is provided, in which at least one instruction or at least one program is stored, and the at least one instruction or the at least one program is loaded and executed by a processor to implement the above-mentioned method for implementing the payment of cultivated land subsidies.
[0012] According to a third aspect of the present invention, there is provided an electronic device comprising a processor and the above-mentioned non-transitory computer-readable storage medium.
[0013] The present invention has at least the following beneficial effects:
[0014] The present invention provides a method for implementing the issuance of cultivated land subsidies. First, each remote sensing image of a target geographical area is divided into a plurality of sub-geographical area images. For each of the plurality of sub-geographical area images corresponding to each sub-geographical area, plot attribute features are extracted for each sub-geographical area image to obtain a target feature vector corresponding to each sub-geographical area. The target feature vector includes a plurality of feature dimensions that can be used to distinguish different geographical area types, so that the identified land object type is more comprehensive and accurate. Then, the target feature vector is input into a target classification model to obtain the geographical area type of the sub-geographical area, and it is determined whether it is a subsidized plot. Each sub-geographical area corresponding to the subsidized plot is divided into a plurality of regional grids, and the geographical area type of each regional grid is determined according to the target feature vector of each regional grid. The issuance of cultivated land subsidies is implemented according to the total area of the regional grids belonging to the subsidized plot. The present invention can improve the estimation accuracy of the area of the subsidized plot by performing the above-mentioned refinement processing on the sub-geographical area belonging to the subsidized plot, thereby realizing the effective issuance of cultivated land subsidies. BRIEF DESCRIPTION OF THE DRAWINGS
[0015] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.
[0016] Figure 1 A flow chart of a method for implementing the issuance of cultivated land subsidies provided in an embodiment of the present invention. DETAILED DESCRIPTION
[0017] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.
[0018] The embodiment of the present invention provides a method for implementing the distribution of cultivated land subsidies, such as Figure 1 As shown, the method comprises the following steps:
[0019] S100, based on a plurality of remote sensing images of the target geographical area acquired within a preset time period, a multi-scale segmentation model is used to divide each remote sensing image of the target geographical area into a plurality of sub-geographical area images, so as to acquire a plurality of sub-geographical area images corresponding to each sub-geographical area.
[0020] Preferably, the preset time period is a time period during which several pre-selected crops are in the growth stage. For example, wheat is planted in November and harvested in May-June, and its intermediate growth stage can be used as the preset time period. Rice is planted in summer and harvested in autumn, and the intermediate growth stage of rice can be used as another preset time period. It should be noted that according to the subsidy issuance principle, subsidies can be issued as long as there is farming within the year. Therefore, images can be taken during the growing seasons of different crops, and the plots that need subsidies can be screened out based on the shooting results of each preset time period.
[0021] Furthermore, the edge detection algorithm in the multi-scale segmentation model is used to divide each remote sensing image into several sub-geographical area images; it can be understood that: according to the different types of crops or landforms in the target geographical area, the target geographical area is divided into several plots.
[0022] Specifically, a plurality of sub-geographical area images corresponding to each sub-geographical area are obtained through the following steps:
[0023] S101, based on the number of sub-geographical area images corresponding to each remote sensing image, the number with the highest frequency among the several numbers is used as the target number corresponding to the sub-geographical area images, so as to screen out several remote sensing images corresponding to the target number. For example, since several remote sensing images of the target geographical area are obtained, some of them are not clear, resulting in different numbers and areas of the segmented sub-geographical area images, and therefore screening is required.
[0024] S102, for the screened remote sensing images, obtain the image center point coordinates of each sub-geographical area image in each remote sensing image, and cluster all the image center point coordinates according to the target number to obtain the target number of image center point clusters; it can be understood as: clustering several sub-geographical area images belonging to the same sub-geographical area into one class.
[0025] Furthermore, the k-means clustering model is used to cluster the coordinates of all image center points, where the value of k is the value of the target number.
[0026] S103, for any image center point cluster, delete the sub-geographic area images corresponding to the image center point coordinates whose distance from the cluster center of the image center point cluster is greater than a first preset distance threshold, and use the sub-geographic area images corresponding to the remaining several image center point coordinates in the image center point cluster as several sub-geographic area images of the same sub-geographic area.
[0027] As mentioned above, considering the problem of remote sensing image clarity, the remote sensing images were screened twice. First, the remote sensing images were screened for the first time according to the number of divided sub-geographical areas, and the remote sensing images obtained were clearer. Then, according to the center points of several sub-geographical area images of each sub-geographical area, the sub-geographical area images with large center point deviations were deleted, so that more accurately divided sub-geographical area images could be obtained, thereby making the extracted features of the sub-geographical areas closer to the actual situation.
[0028] S200 , for a plurality of sub-geographical area images corresponding to each sub-geographical area, extract land parcel attribute features from each sub-geographical area image to obtain a target feature vector corresponding to each sub-geographical area.
[0029] Furthermore, the characteristic dimensions corresponding to the target characteristic vector include normalized vegetation index, green light normalized difference vegetation index, normalized water index, soil adjusted vegetation index, band mean and brightness. The characteristic value corresponding to each characteristic dimension is the average value of the corresponding features extracted from the images of several sub-geographical regions.
[0030] Specifically, the normalized difference vegetation index NDVI = (Nir-Red) / (Nir+Red), wherein Nir is the reflectivity of the near infrared band, and Red is the reflectivity of the red light band.
[0031] Green Normalized Difference Vegetation Index GNDVI = (Nir-Green) / (Nir+Green), where Green is the reflectivity of the green light band.
[0032] Normalized Difference Water Index NDWI = (Green - Nir) / (Green + Nir).
[0033] Soil Adjusted Vegetation Index SAVI = (1 + L) × (Nir - Red) / (Nir + Red + L), where L is the soil adjustment coefficient, ranging from 0 to 1, and is adjusted according to the surface vegetation density.
[0034] As mentioned above, since different geographical areas may have crops, pond water bodies, wastelands, etc., it is necessary to obtain the vegetation, water and soil conditions on the surface of the geographical area. By obtaining the characteristic values corresponding to the above-mentioned multiple characteristic dimensions, the distinction between land feature types can be more comprehensive and accurate, thereby improving the accuracy of geographical area type division.
[0035] S300, for any sub-geographical area, input the target feature vector of the sub-geographical area into the pre-acquired target classification model, and determine the geographical area type of the sub-geographical area according to the output result of the target classification model; the geographical area types include subsidized plots, non-subsidized plots and pending plots. In a specific implementation, subsidized plots include plots for growing grain and oil crops, soybeans, potatoes, and vegetables, non-subsidized plots are mainly plots that have been changed to pond water bodies, construction land, and plots for growing perennial cash crops, such as planting flowers, seedlings, etc., and pending plots are mainly wasteland, simple Dapeng, and other plots. When obtaining geographical area samples, a number of geographical area samples are obtained for each type of plot.
[0036] In a specific embodiment, the target classification model is obtained through the following steps:
[0037] S301, based on a number of geographic area samples obtained within a preset time period, clustering is performed on the feature vector samples corresponding to each geographic area sample, each obtained geographic area sample cluster is used as a sub-sample set to train a number of preset classification models, and a key classification model corresponding to each geographic area sample cluster is determined according to the geographic area type output result of each preset classification model.
[0038] Specifically, the geographic area samples refer to image samples obtained by segmenting remote sensing image samples according to a multi-scale segmentation model.
[0039] Furthermore, the feature vector samples corresponding to the geographic area samples are obtained in the same manner as the target feature vectors corresponding to the sub-geographical areas, which will not be described in detail herein.
[0040] Specifically, the feature vector samples corresponding to each geographic region sample may be clustered using an unsupervised clustering method, for example, hierarchical clustering, K-means clustering or fuzzy clustering, so as to cluster samples with similar feature vectors into the same cluster.
[0041] In a specific embodiment, the key classification model corresponding to each geographical area sample cluster is determined by the following steps:
[0042] S3011, divide the several feature vector samples in each geographic area sample cluster into sub-training sets and sub-testing sets. Since the vegetation index of the same type of ground objects is roughly the same at the same shooting time, and the number of ground object types is small, the proportion of the training set can be appropriately reduced, and the proportion of the sub-training set and the sub-testing set can be set to 1:1.
[0043] S3012, inputting a plurality of feature vector samples in a plurality of sub-training sets into a plurality of preset classification models, so that each preset classification model is learned according to the input plurality of feature vector samples and the geographical area type corresponding to each feature vector sample.
[0044] Preferably, there are at least two types of preset classification models, for example, a neural network model, a random forest prediction model, a logistic regression model, etc.
[0045] S3013, inputting several feature vector samples in several sub-test sets into several learned preset classification models, and calculating the prediction accuracy score of each preset classification model for the geographic area sample cluster based on the difference between the predicted result of the geographic area type corresponding to each feature vector sample and the actual geographic area type corresponding to each feature vector sample.
[0046] Specifically, the prediction accuracy score of any preset classification model for the geographic region sample cluster meets the following conditions:
[0047] F=∑ m i＝1 (2T i / (2T i +N i +P i )) / m, where F represents the prediction accuracy score of the preset classification model for the geographical area sample cluster, T i N is the number of correct predictions among several feature vector samples predicted as the i-th geographic area type, i is the number of incorrect predictions in the feature vector samples predicted as the i-th geographic area type, P i is the number of feature vector samples that are not predicted to be the i-th geographical area type but belong to the i-th geographical area type, m is the number of geographical area types, and the value range of i is 1 to m. In this embodiment, m=3.
[0048] In a specific embodiment, the random forest prediction model is used as an example to illustrate as follows:
[0049] For all sub-training sets, several sub-training sets are randomly selected with replacement to form an input sample set, that is, the feature vector samples in the extracted several sub-training sets are used to form an input sample set, and the number of input sample sets is equal to the number of decision trees in the random forest prediction model. Each input sample set is separately input into a decision tree, so that the decision tree learns the input sample set, and the best feature is selected for splitting according to the information gain each time the split is performed until the training is completed; then for any sub-test set, several feature vector samples in the sub-test set are input into each decision tree model in the random forest prediction model, and the classification result of the geographical area type corresponding to each feature vector sample can be obtained according to any decision tree model, and voting is performed according to the output results of all decision tree models to obtain the final geographical area type corresponding to each feature vector sample, which is compared with the actual geographical area type corresponding to each feature vector sample, so as to obtain the above parameters, so as to calculate the prediction accuracy score of the random forest prediction model for the geographical area sample cluster.
[0050] S3014, using the preset classification model corresponding to the maximum prediction accuracy score as the key classification model corresponding to the geographic area sample cluster.
[0051] As mentioned above, since different preset classification models have different classification effects on different geographic area types, it is necessary to screen out key classification models with good classification effects for each geographic area type. In the above process, the gap between the test results of each preset classification model and the actual geographic area type is taken into account, so that the screened key classification model is more suitable for the corresponding geographic area sample cluster, thereby improving the prediction accuracy in subsequent geographic area type predictions.
[0052] S302: Determine a target sample cluster corresponding to each sub-geographical area from a plurality of geographic area sample clusters according to a target feature vector corresponding to each sub-geographical area.
[0053] Specifically, step S302 includes the following steps:
[0054] S3021, obtaining a central feature vector corresponding to the cluster center of each geographical area sample cluster.
[0055] S3022: For any sub-geographical region, calculate the distance value between the target feature vector corresponding to the sub-geographical region and the central feature vector corresponding to the cluster center of each geographic region sample cluster.
[0056] Furthermore, a direct distance method is used to calculate the distance between the target feature vector and the central feature vector. Those skilled in the art are aware of the specific calculation process of calculating the sample distance in the cluster, which will not be described in detail here.
[0057] S3023: Determine the geographic area sample cluster corresponding to the minimum distance value as the target sample cluster corresponding to the sub-geographic area.
[0058] S303: Acquire a target classification model corresponding to each sub-geographical area according to the target sample cluster corresponding to each sub-geographical area and the key classification model corresponding to each target sample cluster.
[0059] As mentioned above, since the feature vector samples in each sub-geographical area sample cluster are relatively similar and are suitable for a certain preset classification model, by calculating the distance between the target feature vector and the central feature vector, the geographic area sample cluster that is most similar to the target feature vector can be found, thereby finding the preset classification model with the highest accuracy for it, thereby improving the accuracy of judging the geographic area type.
[0060] S400: According to the geographic region type of each sub-geographic region, a plurality of sub-geographic regions corresponding to the subsidized land parcels are screened out, and each of the screened out sub-geographic regions is divided into a plurality of regional grids. For example, the side length of the regional grid is 3 m.
[0061] S500, obtaining the target feature vector corresponding to each regional grid, and inputting the target feature vector corresponding to each regional grid into the target classification model, determining the geographical area type of each regional grid, so as to realize the distribution of cultivated land subsidies according to the geographical area type of each regional grid.
[0062] In a specific embodiment, the method of distributing farmland subsidies according to the geographical area type of each regional grid includes the following steps:
[0063] S501, screening out a number of regional grids belonging to subsidized plots of land according to the geographical area type of each regional grid.
[0064] S502, calculating the total area of the selected regional grids, and distributing farmland subsidies according to the calculation results.
[0065] As mentioned above, since the crops in the subsidized plots are relatively concentrated and continuous plots, most of the plots can be well segmented, but some plots are mixed with small areas of other varieties, especially with unclear edges, and cannot be accurately identified and segmented when the multi-scale segmentation model is used. Therefore, it is necessary to further process the sub-geographical areas to achieve a more accurate estimate of the subsidized area.
[0066] In summary, the present invention provides a method for implementing the issuance of cultivated land subsidies. First, each remote sensing image of the target geographical area is divided into a number of sub-geographical area images. For each of the sub-geographical area images corresponding to each sub-geographical area, land attribute features are extracted for each sub-geographical area image to obtain a target feature vector corresponding to each sub-geographical area. The target feature vector includes a plurality of feature dimensions that can be used to distinguish different geographical area types, so that the identified land object type is more comprehensive and accurate. Then, the target feature vector is input into a target classification model to obtain the geographical area type of the sub-geographical area, and it is determined whether it is a subsidized plot. Each sub-geographical area corresponding to the subsidized plot is divided into a number of regional grids, and the geographical area type of each regional grid is determined according to the target feature vector of each regional grid. The issuance of cultivated land subsidies is implemented according to the total area of the regional grids belonging to the subsidized plot. The present invention can improve the estimation accuracy of the area of the subsidized plot by performing the above-mentioned refinement processing on the sub-geographical area belonging to the subsidized plot, thereby realizing the effective issuance of cultivated land subsidies.
[0067] An embodiment of the present invention also provides a non-transitory computer-readable storage medium, which can be set in an electronic device to store at least one instruction or at least one program related to implementing a method in a method embodiment. The at least one instruction or the at least one program is loaded and executed by the processor to implement the method for implementing the distribution of cultivated land subsidies provided in the above embodiment.
[0068] An embodiment of the present invention further provides an electronic device, comprising a processor and the aforementioned non-transitory computer-readable storage medium.
[0069] Although some specific embodiments of the present invention have been described in detail by way of example, it will be appreciated by those skilled in the art that the above examples are for illustration only and are not intended to limit the scope of the present invention. It will also be appreciated by those skilled in the art that various modifications may be made to the embodiments without departing from the scope and spirit of the present invention. The scope of the present invention is defined by the appended claims.
Claims
1. A method for distributing cultivated land subsidies, characterized in that: The method comprises the following steps: Based on a plurality of remote sensing images of the target geographical area acquired within a preset time period, a multi-scale segmentation model is used to divide each remote sensing image of the target geographical area into a plurality of sub-geographical area images, so as to obtain a plurality of sub-geographical area images corresponding to each sub-geographical area; For a plurality of sub-geographical area images corresponding to each sub-geographical area, land attribute features are extracted from each sub-geographical area image to obtain a target feature vector corresponding to each sub-geographical area; For any sub-geographical area, input the target feature vector of the sub-geographical area into a pre-acquired target classification model, and determine the geographical area type of the sub-geographical area according to the output result of the target classification model; the geographical area type includes subsidized plots, non-subsidized plots and pending plots; According to the geographic area type of each sub-geographical area, a number of sub-geographical areas corresponding to the subsidized plots are screened out, and each of the screened sub-geographical areas is divided into a number of regional grids; The target feature vector corresponding to each regional grid is obtained, and the target feature vector corresponding to each regional grid is input into the target classification model to determine the geographical area type of each regional grid, so as to realize the distribution of cultivated land subsidies according to the geographical area type of each regional grid.
2. The method for implementing the payment of cultivated land subsidies according to claim 1, characterized in that: The following steps are used to obtain several sub-geographical area images corresponding to each sub-geographical area: Based on the number of sub-geographical area images corresponding to each remote sensing image, the number with the highest frequency of occurrence among several numbers is used as the target number corresponding to the sub-geographical area images, so as to screen out several remote sensing images corresponding to the target number; For the selected remote sensing images, the image center point coordinates of each sub-geographical area image in each remote sensing image are obtained, and all the image center point coordinates are clustered according to the target number to obtain the target number of image center point clusters; For any image center point cluster, the sub-geographic area images corresponding to the image center point coordinates whose distance from the cluster center of the image center point cluster is greater than a first preset distance threshold are deleted, and the sub-geographic area images corresponding to the remaining several image center point coordinates in the image center point cluster are taken as several sub-geographic area images of the same sub-geographic area.
3. The method for implementing the payment of cultivated land subsidies according to claim 1, characterized in that: The characteristic dimensions corresponding to the target characteristic vector include normalized vegetation index, green light normalized difference vegetation index, normalized water index, soil adjusted vegetation index, band mean and brightness.
4. The method for implementing the payment of cultivated land subsidies according to claim 1, characterized in that: Obtain the target classification model through the following steps: Based on a number of geographic area samples obtained within a preset time period, clustering is performed on the feature vector samples corresponding to each geographic area sample, using each acquired geographic area sample cluster as a subsample set to train a number of preset classification models, and determining a key classification model corresponding to each geographic area sample cluster according to the geographic area type output result of each preset classification model; According to the target feature vector corresponding to each sub-geographical area, a target sample cluster corresponding to each sub-geographical area is determined from a plurality of geographic area sample clusters; According to the target sample cluster corresponding to each sub-geographical area and the key classification model corresponding to each target sample cluster, the target classification model corresponding to each sub-geographical area is obtained.
5. The method for implementing the payment of cultivated land subsidies according to claim 4, characterized in that: The step of determining a target sample cluster corresponding to each sub-geographical area from a plurality of geographic area sample clusters according to the target feature vector corresponding to each sub-geographical area comprises the following steps: Get the central feature vector corresponding to the cluster center of each geographic area sample cluster: For any sub-geographical area, calculate the distance value between the target feature vector corresponding to the sub-geographical area and the central feature vector corresponding to the cluster center of each geographical area sample cluster; The geographic area sample cluster corresponding to the minimum distance value is determined as the target sample cluster corresponding to the sub-geographic area.
6. The method for implementing the payment of cultivated land subsidies according to claim 1, characterized in that: The method of distributing farmland subsidies according to the geographical area type of each regional grid includes the following steps: According to the geographical area type of each regional grid, a number of regional grids belonging to subsidized plots are screened out; Calculate the total area of several selected regional grids and distribute cultivated land subsidies based on the calculation results.
7. A non-transitory computer-readable storage medium, wherein at least one instruction or at least one program is stored in the storage medium, characterized in that: The at least one instruction or the at least one program is loaded and executed by the processor to implement the method for issuing cultivated land subsidies as described in any one of claims 1-6.
8. An electronic device, characterized in that: The invention comprises a processor and the non-transitory computer-readable storage medium as claimed in claim 7.