A forest area forest land identification method, device, equipment, medium and product

Abstract

The application discloses a forest area forest land identification method, device, equipment, medium and product, relates to the field of territorial space planning, and comprises the following steps: extracting forest land data in territorial change survey data; sequentially performing spatial outward expansion, fusion and spatial inward shrinkage on the forest land data to generate forest land fusion data; respectively performing partition statistics on the forest land data and the forest land fusion data based on a digital elevation model to obtain a forest land data elevation statistical table and a forest land fusion data elevation statistical table; connecting the forest land fusion data elevation statistical table to the forest land fusion data, and performing screening to obtain forest area data; performing spatial intersection analysis on the forest area data and the forest land data to obtain forest area forest land intersection data; connecting the forest land data elevation statistical table to the forest area forest land intersection data, and performing screening to obtain forest area forest land data. The application can improve the accuracy of forest area forest land identification.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of territorial space planning, and in particular to a forest area and forest land identification method, device, equipment, medium and product. BACKGROUND

[0002] Forest resources are important basic data for ecological protection, forestry management and territorial space planning. Accurate identification of forest area and forest land range is of great significance for forest resource monitoring, city health assessment, territorial space use control and carbon sink assessment.

[0003] As the basic data of national natural resource management, the territorial change survey data is authoritative, up-to-date and standardized, and its land class classification system provides reliable data support for forest area and forest land identification. However, the forest land data in the territorial change survey data only reflects the land class attribute, and is disturbed by city street trees and other polygons, and cannot directly identify the forest land in the forest area. Some places take the forest land larger than 100 hectares and concentrated in a piece as the principle of identifying the range of forest area, but this method does not take into account the fragmentation characteristics of forest land polygons in the forest area, nor does it take into account the topographic features, and the identification accuracy is not high. Therefore, how to identify the forest area and forest land based on the territorial change survey data has become a key technical requirement for city health assessment, carbon sink assessment and other fields. SUMMARY

[0004] The purpose of the present application is to provide a forest area and forest land identification method, device, equipment, medium and product, which can improve the accuracy of forest area and forest land identification.

[0005] To achieve the above purpose, the present application provides the following solutions.

[0006] In a first aspect, the present application provides a forest area and forest land identification method, comprising: extracting forest land data in territorial change survey data; sequentially performing spatial outward expansion, fusion and spatial inward contraction on the forest land data to generate forest land fusion data; based on a digital elevation model, respectively performing partition statistics on the forest land data and the forest land fusion data to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table; connecting the forest land fusion data elevation statistics table to the forest land fusion data and performing screening to obtain forest area data; performing spatial intersection analysis on the forest area data and the forest land data to obtain forest area and forest land intersection data; connecting the forest land data elevation statistics table to the forest area and forest land intersection data and performing screening to obtain forest area and forest land data.

[0007] In an embodiment, the forest land data in the territorial change survey data is extracted, specifically comprising: obtaining the territorial change survey data; extracting the territorial change survey data based on the land class name or land class code to obtain the forest land data.

[0008] In an embodiment, the woodland data is sequentially subjected to spatial expansion, fusion and spatial contraction to generate woodland fusion data, specifically comprising: determining a set buffer radius according to the woodland data; performing expansion on the woodland data based on the set buffer radius to generate woodland expansion data; performing spatial fusion on the woodland expansion data to generate woodland expansion fusion data; performing spatial contraction on the woodland expansion fusion data based on the set buffer radius to generate woodland fusion data.

[0009] In an embodiment, the woodland data and the woodland fusion data are respectively subjected to zonal statistics based on a digital elevation model to obtain a woodland data elevation statistics table and a woodland fusion data elevation statistics table, specifically comprising: performing projection conversion on a coordinate system of the digital elevation model to generate a projected version of the digital elevation model consistent with a projection coordinate system of the land change survey data; performing mathematical rounding operation on the projected version of the digital elevation model to generate an integer digital elevation model; performing zonal statistics on the woodland data based on the integer digital elevation model to respectively perform zonal statistics on the woodland data and the woodland fusion data to obtain a woodland data elevation statistics table and a woodland fusion data elevation statistics table.

[0010] In an embodiment, the woodland fusion data elevation statistics table is connected to the woodland fusion data and screened to obtain forest region data, specifically comprising: connecting the woodland fusion data elevation statistics table to the woodland fusion data based on a unique identification code to obtain connected woodland fusion data; calculating the woodland area and the first elevation range difference of each plot in the connected woodland fusion data; the first elevation range difference is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected woodland fusion data; screening the connected woodland fusion data based on the woodland area, the first elevation range difference and the maximum elevation value to determine the forest region data.

[0011] In an embodiment, the woodland data elevation statistics table is connected to the forest region woodland intersection data and screened to obtain forest region woodland data, specifically comprising: connecting the woodland data elevation statistics table to the forest region woodland intersection data based on a unique identification code to obtain connected forest region woodland intersection data; calculating the second elevation range difference of each plot in the connected forest region woodland intersection data; the second elevation range difference is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected forest region woodland intersection data; screening the connected forest region woodland intersection data based on the second elevation range difference to obtain forest region woodland data.

[0012] In a second aspect, the present application provides a forest area forest land identification device, comprising: a forest land extraction module configured to extract forest land data from national land change survey data; a forest land processing module configured to sequentially perform spatial expansion, fusion and spatial contraction on the forest land data to generate forest land fusion data; an elevation statistics module configured to perform partition statistics on the forest land data and the forest land fusion data based on a digital elevation model to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table; a first screening module configured to connect the forest land fusion data elevation statistics table to the forest land fusion data and perform screening to obtain forest area data; a spatial intersection analysis module configured to perform spatial intersection analysis on the forest area data and the forest land data to obtain forest area forest land intersection data; and a second screening module configured to connect the forest land data elevation statistics table to the forest area forest land intersection data and perform screening to obtain forest area forest land data.

[0013] In a third aspect, the present application provides a computer device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the forest area forest land identification method.

[0014] In a fourth aspect, the present application provides a computer readable storage medium having a computer program stored thereon, wherein the computer program is executable by a processor to implement the forest area forest land identification method.

[0015] In a fifth aspect, the present application provides a computer program product comprising a computer program executable by a processor to implement the forest area forest land identification method.

[0016] According to the specific embodiments provided by the present application, the following technical effects are disclosed:

[0017] The present application provides a forest area forest land identification method, device, equipment, medium and product, which extracts forest land data from national land change data, optimizes the forest land data through spatial expansion, fusion and spatial contraction, and performs partition statistics in combination with a digital elevation model, connects the forest land fusion data elevation statistics table to the forest land fusion data, and performs screening to obtain forest area data; performs spatial intersection analysis on the forest area data and the forest land data to obtain forest area forest land intersection data; connects the forest land data elevation statistics table to the forest area forest land intersection data, and performs screening to obtain forest area forest land data, thereby improving the accuracy of forest area forest land identification. BRIEF DESCRIPTION OF DRAWINGS

[0018] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed in the embodiments. Obviously, the drawings described below only illustrate some of the embodiments of the present application, and all other drawings obtained by those of ordinary skill in the art without creative effort based on these drawings belong to the protection scope of the present application.

[0019] Figure 1 A schematic diagram of a forest land recognition method.

[0020] Figure 2 A schematic diagram of two map patches of forest land data.

[0021] Figure 3 A schematic diagram of two map patches of forest land expansion data generated by performing an outer expansion operation on the two map patches of forest land data through buffer analysis.

[0022] Figure 4 A schematic diagram of a map patch of forest land expansion fusion data generated by performing a spatial fusion operation on the two map patches of forest land expansion data.

[0023] Figure 5 A schematic diagram of a map patch of forest land fusion data generated by performing an inner shrinking operation on the map patch of forest land expansion fusion data in an application scenario.

[0024] Figure 6 A flowchart of a forest land recognition method.

[0025] Figure 7 A schematic diagram of a function module of a forest land recognition device.

[0026] Figure 8 A structural schematic diagram of a computer device provided by an embodiment of the present application. DETAILED DESCRIPTION

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

[0028] In order to make the above-mentioned purposes, features and advantages of the present application more obvious and easy to understand, the present application will be further described in detail below with reference to the drawings and specific embodiments.

[0029] In an exemplary embodiment, as Figure 6As shown, a forest area forest land identification method is provided, which is executed by a computer device, specifically, can be executed by a terminal or a server, or can be executed by a terminal and a server together, in the embodiments of the present application, the method is applied to the server as an example, and includes the following steps.

[0030] Step 601: Extract forest land data in the land change survey data.

[0031] Step 602: Perform spatial expansion, fusion and spatial shrinkage on the forest land data in sequence to generate forest land fusion data.

[0032] Step 603: Based on the digital elevation model, respectively perform partition statistics on the forest land data and the forest land fusion data to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table.

[0033] Step 604: Connect the forest land fusion data elevation statistics table to the forest land fusion data and perform screening to obtain forest area data.

[0034] Step 605: Perform spatial intersection analysis on the forest area data and the forest land data to obtain forest area forest land intersection data.

[0035] Step 606: Connect the forest land data elevation statistics table to the forest area forest land intersection data and perform screening to obtain forest area forest land data.

[0036] The forest land data is extracted by using the specification of the land change data, the forest land data is optimized by spatial expansion, fusion and spatial shrinkage, and then the digital elevation model is used for partition statistics, the forest land fusion data elevation statistics table is connected to the forest land fusion data, and screening is performed to obtain forest area data, the forest area data and the forest land data are subjected to spatial intersection analysis to obtain forest area forest land intersection data, and the forest land data elevation statistics table is connected to the forest area forest land intersection data, and screening is performed to obtain forest area forest land data, thereby improving the accuracy of forest area forest land identification.

[0037] In one exemplary embodiment, the forest land data in the land change survey data is extracted, specifically including: acquiring land change survey data; extracting the land change survey data based on the land class name or land class code to obtain forest land data.

[0038] In an exemplary embodiment, the woodland data is sequentially subjected to spatial expansion, fusion and spatial contraction to generate woodland fusion data, specifically comprising: determining a set buffer radius according to the woodland data; expanding the woodland data based on the set buffer radius to generate woodland expansion data; performing spatial fusion on the woodland expansion data to generate woodland expansion fusion data; and performing spatial contraction on the woodland expansion fusion data based on the set buffer radius to generate the woodland fusion data.

[0039] In an exemplary embodiment, the woodland data and the woodland fusion data are respectively subjected to zonal statistics based on a digital elevation model to obtain a woodland data elevation statistics table and a woodland fusion data elevation statistics table, specifically comprising: performing projection conversion on a coordinate system of the digital elevation model to generate a projected version of the digital elevation model consistent with a projection coordinate system of the land change survey data; performing mathematical rounding operation on the projected version of the digital elevation model to generate an integer digital elevation model; and performing zonal statistics on the woodland data based on the integer digital elevation model to obtain the woodland data elevation statistics table and the woodland fusion data elevation statistics table.

[0040] In an exemplary embodiment, the woodland fusion data elevation statistics table is connected to the woodland fusion data, and screening is performed to obtain woodland data, specifically comprising: connecting the woodland fusion data elevation statistics table to the woodland fusion data based on a unique identification code to obtain connected woodland fusion data; calculating the woodland area and the first elevation range difference of each plot in the connected woodland fusion data; the first elevation range difference is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected woodland fusion data; and performing screening based on the woodland area, the first elevation range difference and the maximum elevation value according to the connected woodland fusion data to determine the woodland data.

[0041] In an exemplary embodiment, the woodland data elevation statistics table is connected to the woodland intersection data, and screening is performed to obtain woodland data, specifically comprising: connecting the woodland data elevation statistics table to the woodland intersection data based on a unique identification code to obtain connected woodland intersection data; calculating the second elevation range difference of each plot in the connected woodland intersection data; the second elevation range difference is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected woodland intersection data; and performing screening based on the second elevation range difference according to the connected woodland intersection data to obtain the woodland data.

[0042] The application utilizes land change survey data and a digital elevation model, adopts technical means such as 'expansion-fusion-constriction' optimization, elevation range statistics and multi-threshold dynamic screening, and can accurately identify forest areas and forest land from land change survey data, thereby providing reliable data support for city health assessment, carbon sink assessment and the like.

[0043] As shown in Figure 1 As another exemplary embodiment, a specific process of the forest area and forest land identification method in actual application is also provided, including the following steps.

[0044] Forest land data is extracted from land change survey data according to land class names or land class codes.

[0045] Geographic Information System (GIS) spatial expansion, fusion and constriction operations are sequentially performed on the forest land data to generate forest land fusion data.

[0046] Based on a digital elevation model, partition statistics in GIS spatial analysis are respectively performed on the forest land data and the forest land fusion data, and relevant elevation information is counted, so as to finally obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table.

[0047] The forest land fusion data elevation statistics table is connected to the forest land fusion data, and forest area data is screened out according to elevation threshold parameters, area threshold parameters and the like.

[0048] GIS spatial intersection analysis is performed on the forest area data and the forest land data to obtain forest area and forest land intersection data.

[0049] The forest land fusion data elevation statistics table is connected to the forest area and forest land intersection data, and forest area and forest land data is screened out through an elevation threshold parameter.

[0050] Further, forest land data is extracted from land change survey data according to land class names or land class codes, specifically including: extracting map patches with land class names of 'tree forest land' or'shrub forest land' or 'bamboo forest land' or 'other forest land' from land change survey data to form forest land data; or extracting map patches with land class codes of '0301' or '0302' or '0305' or '0307' from land change survey data to form forest land data.

[0051] Further, the forest land data is sequentially subjected to GIS spatial expansion, fusion and contraction operations to generate forest land fusion data, specifically including: setting a buffer radius according to the forest land data; performing an expansion operation on the forest land data through GIS buffer analysis to generate forest land expansion data, wherein the buffer radius is the set buffer radius; performing a spatial fusion operation on the forest land expansion data to generate forest land expansion fusion data; performing a contraction operation on the forest land expansion fusion data through GIS buffer analysis to generate forest land fusion data, wherein the buffer radius is the negative of the set buffer radius; and assigning a unique identification code to the forest land fusion data.

[0052] The method for generating forest land fusion data is further described below through specific application scenarios.

[0053] Figure 2 For two polygons of forest land data, Figure 3 For two polygons of forest land expansion data generated by performing a buffer analysis operation on the two polygons of forest land data, Figure 4 For a polygon of forest land expansion fusion data generated by performing a spatial fusion operation on the two polygons of forest land expansion data, Figure 5 For a polygon of forest land fusion data generated by performing a contraction operation on the polygon of forest land expansion fusion data.

[0054] Further, based on a digital elevation model, the forest land data and the forest land fusion data are subjected to partition statistics in GIS spatial analysis, and relevant elevation information is counted to ultimately obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table, specifically including: performing a projection conversion operation on the coordinate system of the digital elevation model to generate a projected version of the digital elevation model that is consistent with the projection coordinate system of the national land change survey data; performing a mathematical rounding operation on the projected version of the digital elevation model to generate an integer digital elevation model; combining the integer digital elevation model, performing partition statistics in GIS spatial analysis on the forest land data, and counting relevant elevation information to obtain the forest land data elevation statistics table, wherein the relevant elevation information includes the maximum elevation value within a polygon and the minimum elevation value within a polygon, and the forest land data elevation statistics table should also retain the unique identification code information in the forest land data; combining the integer digital elevation model, performing partition statistics in GIS spatial analysis on the forest land fusion data, and counting relevant elevation information to obtain the forest land fusion data elevation statistics table, wherein the relevant elevation information includes the maximum elevation value within a polygon and the minimum elevation value within a polygon, and the forest land fusion data elevation statistics table should also retain the unique identification code information in the forest land fusion data. The forest land data elevation statistics table and the forest land fusion data elevation statistics table both include relevant elevation information.

[0055] Further, the forest land data elevation statistics table is connected to the forest land fusion data, and forest area data is screened according to an elevation threshold parameter, an area threshold parameter, etc., specifically including: connecting the forest land data elevation statistics table to the forest land fusion data through a unique identification code shared by the two; performing an area calculation operation on the connected forest land fusion data to obtain forest area; performing an elevation range range difference calculation on the connected forest land fusion data, the elevation range range difference being a difference between a maximum elevation value and a minimum elevation value in each polygon; determining a maximum elevation threshold, an elevation range range difference threshold, and a forest area threshold according to an actual situation of a research area; screening, from the connected forest land fusion data, polygons whose maximum elevation value exceeds the maximum elevation threshold to generate forest land data meeting the maximum elevation threshold; screening, from the connected forest land fusion data, polygons whose elevation range range difference exceeds the elevation range range difference threshold to generate forest land data meeting the elevation range range difference threshold; screening, from the connected forest land fusion data, polygons whose forest area exceeds the forest area threshold to generate forest land data meeting the area threshold; and taking an intersection of the forest land data meeting the maximum elevation threshold, the forest land data meeting the elevation range range difference threshold, and the forest land data meeting the area threshold to generate forest area data.

[0056] Further, the forest land data elevation statistics table is connected to the forest land fusion data, and forest area data is screened according to an elevation threshold parameter, an area threshold parameter, etc., specifically including: connecting the forest land data elevation statistics table to the forest land fusion data through a unique identification code shared by the two; performing an area calculation operation on the connected forest land fusion data to obtain forest area; performing an elevation range range difference calculation on the connected forest land fusion data, the elevation range range difference being a difference between a maximum elevation value and a minimum elevation value in each polygon; determining a maximum elevation threshold, an elevation range range difference threshold, and a forest area threshold according to an actual situation of a research area; screening, from the connected forest land fusion data, polygons whose maximum elevation value exceeds the maximum elevation threshold to generate forest land data meeting the maximum elevation threshold; screening, from the connected forest land fusion data, polygons whose elevation range range difference exceeds the elevation range range difference threshold to generate forest land data meeting the elevation range range difference threshold; screening, from the connected forest land fusion data, polygons whose forest area exceeds the forest area threshold to generate forest land data meeting the area threshold; and taking an intersection of the forest land data meeting the maximum elevation threshold, the forest land data meeting the elevation range range difference threshold, and the forest land data meeting the area threshold to generate forest area data.

[0057] The application extracts forest land data from land change survey data according to land class names or land class codes; performs GIS spatial expansion, fusion and shrinkage operations on the forest land data in sequence to generate forest land fusion data; performs partition statistics in GIS spatial analysis on the forest land data and the forest land fusion data respectively based on a digital elevation model, and counts relevant elevation information, to finally obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table; connects the forest land fusion data elevation statistics table to the forest land fusion data, and filters out forest area data according to an elevation threshold parameter, an area threshold parameter and the like; performs GIS spatial intersection analysis on the forest area data and the forest land data to obtain forest area forest land intersection data; connects the forest land data elevation statistics table to the forest area forest land intersection data, and filters out forest area forest land data through an elevation threshold parameter; and uses land change survey data and a digital elevation model to accurately identify forest area forest land from land change survey data, thereby providing reliable data support for city health assessment, carbon sink assessment and the like.

[0058] Based on the same inventive concept, the embodiments of the application also provide a forest area forest land identification device for implementing the forest area forest land identification method described above. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme described in the above method, so the specific limitations in one or more forest area forest land identification device embodiments provided below can refer to the limitations of the forest area forest land identification method described above, which will not be repeated here.

[0059] In one exemplary embodiment, as shown in Figure 7 A forest area forest land identification device is provided, including: a forest land extraction module configured to extract forest land data from land change survey data; a forest land processing module configured to perform spatial expansion, fusion and spatial shrinkage on the forest land data in sequence to generate forest land fusion data; an elevation statistics module configured to perform partition statistics on the forest land data and the forest land fusion data respectively based on a digital elevation model to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table; a first filtering module configured to connect the forest land fusion data elevation statistics table to the forest land fusion data and perform filtering to obtain forest area data; a spatial intersection analysis module configured to perform spatial intersection analysis on the forest area data and the forest land data to obtain forest area forest land intersection data; and a second filtering module configured to connect the forest land data elevation statistics table to the forest area forest land intersection data and perform filtering to obtain forest area forest land data.

[0060] In an exemplary embodiment, a computer device is provided, which can be a server or a terminal, and an internal structure diagram thereof can be as shown in Figure 8As shown in the figure. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. Among them, the processor, the memory and the input / output interface are connected through the system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capability. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store forest area forest land identification data. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with the terminal outside through network connection. The computer program is executed by the processor to realize a forest area forest land identification method.

[0061] Those skilled in the art can understand that, Figure 8 The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement. In one exemplary embodiment, a computer device is provided, including a memory and a processor, the memory storing a computer program, and the processor executing the computer program to realize the above-mentioned method embodiments.

[0062] In one exemplary embodiment, a computer readable storage medium is provided, storing a computer program, which is executed by a processor to realize the above-mentioned method embodiments.

[0063] In one exemplary embodiment, a computer program product is provided, including a computer program, which is executed by a processor to realize the above-mentioned method embodiments.

[0064] It should be noted that the user information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties, and the collection, use and processing of related data need to comply with relevant regulations.

[0065] In the present application, all actions of obtaining signals, information or data are carried out in compliance with the corresponding data protection regulations and policies of the country where the device is located, and with the authorization given by the corresponding device owner.

[0066] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of the methods can be included. Any reference to a memory, a database or other medium used in the embodiments provided in the present application can include at least one of a non-volatile and a volatile memory. The non-volatile memory can include a read-only memory (ROM), a magnetic tape, a floppy disk, a flash memory, an optical storage, a high-density embedded non-volatile memory, a resistive random access memory (ReRAM), a magnetoresistive random access memory (MRAM), a ferroelectric random access memory (FRAM), a phase change memory (PCM), a graphene memory, etc. The volatile memory can include a random access memory (RAM) or an external cache memory, etc. As an illustration but not limitation, the RAM can be in various forms, such as a static random access memory (SRAM) or a dynamic random access memory (DRAM), etc.

[0067] The database involved in the embodiments provided in the present application can include at least one of a relational database and a non-relational database. The non-relational database can include a distributed database based on a blockchain, etc., without being limited thereto. The processor involved in the embodiments provided in the present application can be a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., without being limited thereto.

[0068] The technical features of the above embodiments can be combined arbitrarily. In order to make the description concise, all possible combinations of the technical features in the above embodiments are not described, however, as long as the combinations of the technical features do not exist contradictory, they should be considered as the scope of the present application.

[0069] The principles and implementation modes of the present application are described by applying specific examples in the present application. The above-mentioned embodiments are only used to help understand the method and its core idea of the present application; meanwhile, for those skilled in the art, according to the idea of the present application, the specific implementation mode and application range can be changed. In conclusion, the content of the present application should not be understood as a limitation.

Claims

1. A method for identifying forest land in a forest area, characterized by, The forest region forest land identification method comprises the following steps: extracting forest land data in the land change survey data; sequentially performing spatial outward expansion, fusion and spatial inward shrinkage on the forest land data to generate forest land fusion data; based on a digital elevation model, respectively performing partition statistics on the forest land data and the forest land fusion data to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table; connecting the forest land fusion data elevation statistics table to the forest land fusion data and performing screening to obtain forest region data; connecting the forest land fusion data elevation statistics table to the forest land fusion data and performing screening to obtain forest region data, specifically comprising: connecting the forest land fusion data elevation statistics table to the forest land fusion data based on a unique identification code to obtain connected forest land fusion data; calculating the forest area and the first elevation range range of each plot in the connected forest land fusion data; the first elevation range range is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected forest land fusion data; screening based on the forest area, the first elevation range range and the maximum elevation value according to the connected forest land fusion data to determine the forest region data; performing spatial intersection analysis on the forest region data and the forest land data to obtain forest region forest land intersection data; connecting the forest land data elevation statistics table to the forest region forest land intersection data and performing screening to obtain forest region forest land data; connecting the forest land data elevation statistics table to the forest region forest land intersection data and performing screening to obtain forest region forest land data, specifically comprising: connecting the forest land data elevation statistics table to the forest region forest land intersection data based on a unique identification code to obtain connected forest region forest land intersection data; calculating the second elevation range range of each plot in the connected forest region forest land intersection data; the second elevation range range is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected forest region forest land intersection data; screening based on the second elevation range range according to the connected forest region forest land intersection data to obtain forest region forest land data.

2. The forest area forest land identification method according to claim 1, characterized by, extracting forest land data in the land change survey data, specifically comprising: obtaining land change survey data; extracting the land change survey data based on a land class name or a land class code to obtain forest land data.

3. The forest area forest land identification method according to claim 1, characterized by, sequentially performing spatial outward expansion, fusion and spatial inward shrinkage on the forest land data to generate forest land fusion data, specifically comprising: determining a set buffer radius according to the forest land data; performing outward expansion on the forest land data based on the set buffer radius to generate forest land outward expansion data; performing spatial fusion on the forest land outward expansion data to generate forest land outward expansion fusion data; performing spatial inward shrinkage on the forest land outward expansion fusion data based on the set buffer radius to generate forest land fusion data.

4. The forest area forest land identification method according to claim 1, characterized by, based on a digital elevation model, respectively performing partition statistics on the forest land data and the forest land fusion data to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table, specifically comprising: performing projection conversion on the coordinate system of the digital elevation model to generate a digital elevation model projection version consistent with the projection coordinate system of the land change survey data; Perform a mathematical rounding operation on the digital elevation model projection version to generate an integer digital elevation model; Based on the integer digital elevation model, the forest land data and the forest land fusion data are respectively subjected to zonal statistics to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table.

5. A forest land identification device characterized by comprising: The forest land recognition device comprises: A forest land extraction module is configured to extract forest land data from national land change survey data; A forest land processing module is configured to sequentially perform spatial outward expansion, fusion and spatial inward contraction on the forest land data to generate forest land fusion data; An elevation statistics module is configured to perform zonal statistics on the forest land data and the forest land fusion data based on a digital elevation model to obtain a forest land data elevation statistics table and a forest land fusion data elevation statistics table; A first screening module is configured to connect the forest land fusion data elevation statistics table to the forest land fusion data and perform screening to obtain forest region data, and to connect the forest land fusion data elevation statistics table to the forest land fusion data and perform screening to obtain forest region data, specifically comprising: connecting the forest land fusion data elevation statistics table to the forest land fusion data based on a unique identification code to obtain connected forest land fusion data; calculating the forest land area and a first elevation range range of each plot in the connected forest land fusion data; the first elevation range range is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected forest land fusion data; and performing screening based on the forest land area, the first elevation range range and the maximum elevation value according to the connected forest land fusion data to determine forest region data; A spatial intersection analysis module is configured to perform spatial intersection analysis on the forest region data and the forest land data to obtain forest region forest land intersection data; A second screening module is configured to connect the forest land data elevation statistics table to the forest region forest land intersection data and perform screening to obtain forest region forest land data, and to connect the forest land data elevation statistics table to the forest region forest land intersection data and perform screening to obtain forest region forest land data, specifically comprising: connecting the forest land data elevation statistics table to the forest region forest land intersection data based on a unique identification code to obtain connected forest region forest land intersection data; calculating a second elevation range range of each plot in the connected forest region forest land intersection data; the second elevation range range is the difference between the maximum elevation value and the minimum elevation value of each plot in the connected forest region forest land intersection data; and performing screening based on the second elevation range range according to the connected forest region forest land intersection data to obtain forest region forest land data.

6. A computer device comprising: A memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor executes the computer program to implement the forest region forest land recognition method of any one of claims 1-4.

7. A computer readable storage medium having stored thereon a computer program, characterized in that The computer program is executed by the processor to implement the forest region forest land recognition method of any one of claims 1-4.

8. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to implement the forest region forest land recognition method of any one of claims 1-4.

Images