CN117095042B - Method and device for determining earth type transfer space difference of earth surface coverage - Google Patents
Method and device for determining earth type transfer space difference of earth surface coverage Download PDFInfo
- Publication number
- CN117095042B CN117095042B CN202310967151.3A CN202310967151A CN117095042B CN 117095042 B CN117095042 B CN 117095042B CN 202310967151 A CN202310967151 A CN 202310967151A CN 117095042 B CN117095042 B CN 117095042B
- Authority
- CN
- China
- Prior art keywords
- data
- index data
- class transfer
- ground
- earth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012546 transfer Methods 0.000 title claims abstract description 187
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000008859 change Effects 0.000 claims abstract description 51
- 238000012216 screening Methods 0.000 claims abstract description 29
- 238000004458 analytical method Methods 0.000 claims description 40
- 238000007621 cluster analysis Methods 0.000 claims description 30
- 230000007704 transition Effects 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000006870 function Effects 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 7
- 238000007619 statistical method Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 238000006116 polymerization reaction Methods 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 21
- 238000009826 distribution Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 238000007620 mathematical function Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012732 spatial analysis Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/762—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using clustering, e.g. of similar faces in social networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30184—Infrastructure
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The application provides a method and a device for determining earth type transfer spatial diversity of earth surface coverage, and relates to the field of geographic measurement. The method for determining the earth type transfer spatial diversity of the earth surface coverage comprises the following steps: obtaining index data of the earth class transfer according to historical earth surface coverage data of the target area; obtaining a characteristic curve based on the index data; obtaining target index data according to the index data corresponding to the part of the characteristic curve, wherein the slope change rate of the part does not exceed the screening threshold value; determining the spatial diversity of the ground class transfer of the target area according to the target index data; the space diversity determining method provided by the embodiment of the application can determine the space diversity situation more accurately, and is beneficial to the management of land resources.
Description
Technical Field
The application relates to the field of geographic measurement, in particular to a method and a device for determining earth type transfer spatial diversity of earth surface coverage.
Background
Land Cover (Land Cover) refers to the type of different features covered by different areas of the earth's surface, such as forests, grasslands, farmlands, urban structures, etc. The distribution and variation of surface coverage reflects the progress of economic activities of the human society. Land Use Change refers to a transition or Change process between surface coverage types that describes the transition of a region from one surface coverage type to another. Spatial diversity (Spatial Differentiation) refers to the difference in surface coverage or earth transfer between different regions on the earth's surface.
The determination of the space difference of the ground transfer not only can embody the changes of nature, humanity, society, economy, ecology and the like of each region, but also can provide guidance and suggestion for the data production in the future; at present, in the process of judging the space diversity situation of the ground class transfer, the obtained ground class transfer matrix is abnormally sparse, the distribution situation of the ground class transfer cannot be determined from the all-ground class change layer, and the effectiveness and the global property of the space diversity of the ground class transfer are to be improved.
Disclosure of Invention
The embodiment of the application aims to provide a method and a device for determining the space diversity of class transfer of a surface coverage land, which are used for extracting main class transfer from the existing change condition of the class transfer and performing cluster analysis, outlier analysis and the like based on the main class transfer, so that the accuracy of determining the space diversity of the class transfer is effectively improved.
In a first aspect, an embodiment of the present application provides a method for determining a ground class transfer spatial diversity of a surface coverage, where the method includes: obtaining index data of the earth class transfer according to historical earth surface coverage data of the target area; obtaining a characteristic curve based on the index data; obtaining target index data according to the index data corresponding to the part of the characteristic curve, wherein the slope change rate of the part does not exceed the screening threshold value; and determining the spatial diversity of the ground class transition of the target area according to the target index data.
In the implementation process, the embodiment of the application obtains the index data of the land class transfer through the data of the earth surface coverage for a plurality of years; the main ground type transfer types are accurately screened out by calculating the slope change rate of the characteristic curve fitted based on the index data, the ground type transfer space difference situation of the target area is accurately revealed, and useful information and support are provided for land utilization planning and resource management.
Optionally, in an embodiment of the present application, obtaining the index data of the earth class transfer according to the historical earth surface coverage data of the target area includes: according to the earth surface coverage data, calculating a ground class transfer matrix of the target area; and respectively obtaining target area index data based on the ground class transfer matrix.
In the implementation process, the embodiment of the application can further analyze and obtain index data reflecting the spatial diversity of the regional class transfer of the target area based on the method for calculating the regional class transfer matrix between years based on the historical earth surface coverage data of many years. Based on the historical data, a land class transfer matrix is used for future land utilization planning and prediction.
Optionally, in an embodiment of the present application, obtaining the characteristic curve based on the index data includes: arranging index data in a reverse order, and integrating each index data to obtain a plurality of integrated data; performing curve fitting on the plurality of integral data to obtain a characteristic curve; wherein the characteristic curve comprises a logarithmic function curve or an inverse proportion function curve.
In the implementation process, a plurality of integral data can be obtained by integrating the index data; the method can reflect the cumulative effect in the surface coverage type transfer process, so that the characteristics of the surface coverage type transfer can be extracted and characterized better. And curve fitting is carried out on the plurality of integral data, so that a proper mathematical function can be found to describe the characteristic curve of the surface coverage type transfer, and the dynamic change rule of the surface coverage type transfer is better expressed.
Optionally, in an embodiment of the present application, the index data includes ground class transfer number data and ground class transfer area data; the characteristic curves comprise a ground class transfer quantity characteristic curve and a ground class transfer area characteristic curve; obtaining target index data according to the index data corresponding to the part of the characteristic curve with the slope change rate not exceeding the screening threshold value, wherein the target index data comprises: and acquiring the intersection of the ground class transfer quantity data corresponding to the part of the ground class transfer quantity characteristic curve, the change rate of which does not exceed the screening threshold value, and the ground class transfer area data corresponding to the part of the ground class transfer area characteristic curve, the change rate of which does not exceed the screening threshold value, so as to obtain target index data.
In the implementation process, the method for determining the space diversity of the earth surface coverage earth type transfer provided by the embodiment of the application can be used for more accurately determining the space diversity of the earth type transfer by more finely screening important earth type transfer types by considering two aspects of the earth type transfer quantity and the earth type transfer area.
Optionally, in an embodiment of the present application, analyzing, according to the target index data, a spatial variance of a ground class transition of the target area includes: performing cluster analysis on the target index data to obtain a cluster analysis result; and determining the spatial diversity of the ground class transfer of the target area according to the clustering analysis result.
Optionally, in an embodiment of the present application, analyzing, according to a result of the cluster analysis, a spatial variance of a ground class transition of the target area includes: and respectively carrying out heat analysis on the ground class transfer quantity and the ground class transfer area according to the clustering analysis result.
In the implementation process, the embodiment of the application can reveal the spatial distribution condition of different ground class transfer types in the target area through cluster analysis and ground class average area statistics, and know the characteristics and differences of the ground class transfer in different areas. The method is not only suitable for the analysis of the earth class transfer of the target area, but also can split the nationwide data into corresponding subsets for the spatial analysis of the nationwide earth class transfer, and provides a wider research view angle.
Optionally, in an embodiment of the present application, before performing cluster analysis on the target index data to obtain a cluster analysis result, the method further includes: statistical analysis and polymerizability analysis are performed on the target index data to evaluate the data quality of the target index data.
In the implementation process, the earth surface coverage earth type transfer space differentiation method provided by the embodiment of the application evaluates the quality and reliability of target index data through statistical analysis and polymerizability analysis, and avoids inaccurate or unreliable clustering results caused by poor data quality.
In a second aspect, an embodiment of the present application provides a device for determining a difference in a ground class transfer space covered by a ground surface, where the device includes: the system comprises a characteristic curve generation module, a target data acquisition module and an analysis result generation module; the characteristic curve generation module is used for obtaining index data of the ground class transfer according to historical ground surface coverage data of the target area; obtaining a characteristic curve based on the index data; the target data acquisition module is used for acquiring target index data according to the index data corresponding to the part of the characteristic curve, the slope change rate of which does not exceed the screening threshold value; the analysis result generation module is used for determining the spatial difference of the ground class transition of the target area according to the target index data.
In a third aspect, an embodiment of the present application provides an electronic device, where the electronic device includes a memory and a processor, where the memory stores program instructions, and when the processor reads and executes the program instructions, the processor performs the steps in any of the foregoing implementation manners.
In a fourth aspect, embodiments of the present application further provide a computer readable storage medium having stored therein computer program instructions that, when read and executed by a processor, perform the steps of any of the above implementations.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments of the present application will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of a method for determining a ground class transition spatial diversity of a surface coverage according to an embodiment of the present application;
FIG. 2 is a flow chart of generating a characteristic curve according to an embodiment of the present application;
fig. 3 is a flowchart of spatial diversity acquisition provided in an embodiment of the present application;
fig. 4 is a schematic block diagram of a determining device for determining a ground type transfer spatial diversity of earth surface coverage according to an embodiment of the present application;
fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application. For example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions. In addition, functional modules in the embodiments of the present invention may be integrated together to form a single part, or each module may exist alone, or two or more modules may be integrated to form a single part.
The applicant finds that the earth surface coverage is carefully classified in the research process, and the earth surface coverage is as high as hundreds of types, if the annual land class transfer condition data is taken as a research sample, the quantity of the earth surface coverage is exponentially increased after the earth surface coverage is arranged and combined; the ground class transfer matrix obtained by combining the analysis unit and the ground class transfer change is abnormally sparse, the distribution condition cannot be analyzed from the all-ground class change layer, and the effective result cannot be obtained by cluster analysis.
Based on the above, the application provides a method for determining the space diversity of the earth surface covered earth transfer, which extracts main earth transfer from the existing earth transfer change condition and performs cluster analysis, outlier analysis and the like based on the main earth transfer, thereby effectively improving the accuracy of determining the space diversity of the earth transfer.
Before the present application, the basic concepts of earth surface coverage, earth class transfer, and spatial diversity will be briefly described.
Earth's surface coverage is a complex of various types of matter on the earth's surface and its natural attributes and characteristics, the distribution and variation of which reflect the progress of economic activities of the human society. The resource and ecological environment change trend can be obtained through the surface coverage, effective geographic information support is provided for constructing an ecological safety strategy pattern and enhancing ecological environment protection and treatment, and important basis is provided for regional planning, climate change, ecological system evaluation and other national important strategies. Therefore, a method for determining the distribution condition of the earth surface coverage earth type change is established, objective and practical earth type change can be intuitively embodied, and directions can be provided for data production.
The earth class transfer refers to the process of conversion or change between earth surface coverage types. It describes the transition of a region from one surface coverage type to another. For example, forests are cut down for use as farmlands, farmlands are replaced by city developments, and the like.
Spatial diversity refers to the difference in surface coverage or earth transfer between different regions on the earth's surface. The spatial diversity can reveal the development change existing in the objective world, and the development change conditions in different areas are different; the determination of the spatial diversity not only can embody the changes of nature, humane, society, economy, ecology and the like in each region, but also can provide guidance and suggestion for the data production in the future, for example, the ground surface coverage classification and the like can be finely adjusted according to the change condition of development and the important monitoring ground surface features so as to meet the actual region, for example, the ground surface coverage classification is finely adjusted according to the change condition of the abnormal ground surface features.
Referring to fig. 1, fig. 1 is a flowchart of a method for determining a difference in earth type transfer space of earth surface coverage according to an embodiment of the present application; the method for determining the earth type transfer spatial diversity of the earth surface coverage comprises the following steps:
step S100: and obtaining index data of the earth class transfer according to the historical earth surface coverage data of the target area.
In the above step S100, the index data of the earth class transition is obtained from the historic earth surface coverage data of the target area, wherein the index data may be understood as the feature data reflecting the spatial diversity of the earth class transition.
The historical earth surface coverage data studied in the present application refers to earth surface coverage data of a period of time, and may be earth surface coverage data of years, earth surface coverage data of months, earth surface coverage data of days, or the like; in the embodiment of the application, only the data covered by the earth surface for years is used as the data sample for analysis, and other sample data are analyzed by using the determination method for the earth surface covered by the earth surface for the transition space diversity, which is provided by the embodiment of the application, within the protection scope of the application.
Step S200: based on the index data, a characteristic curve is obtained.
In step S200, a characteristic curve is obtained based on index data reflecting the difference in the class transition space.
Step S300: and obtaining target index data according to the index data corresponding to the part of the characteristic curve, wherein the slope change rate of the part does not exceed the screening threshold value.
In the step S300, after the curve is fitted, the slope change rate of the characteristic curve is calculated; further, index data corresponding to a portion where the change rate of the screening slope does not exceed the screening threshold is screened as target index data.
It should be noted that, in the embodiment of the present application, the purpose of setting the screening threshold is to screen the required index data as an important analysis object; for example, it is intended in the present application to screen out major ground class transitions.
In the above procedure, the formula may be:
calculating a slope change rate, wherein Δk represents the slope change rate, Δy' i+1 ，Δy‘ i Indicating the slope change of points i+1 and i, y' i+1 、y‘ i 、y‘ i-1 The slopes of points i+1, i-1 are shown, and y is the ordinate of the characteristic curve.
Step S400: and determining the spatial diversity of the ground class transition of the target area according to the target index data.
In the step S400, the spatial diversity of the land class transfer of the target area is determined according to the screened target index data, for example, the land class transfer type is mainly distributed in a certain place, or the area of the land class transfer type in a certain place is determined.
As can be seen from fig. 1, in the embodiment of the present application, the data of the earth class transfer index is obtained by covering the data with earth surface over a plurality of years; the main ground type transfer types are accurately screened out by calculating the slope change rate of the characteristic curve fitted based on the index data, the ground type transfer space difference situation of the target area is accurately revealed, and useful information and support are provided for land utilization planning and resource management.
In an alternative embodiment, obtaining the index data of the earth class transfer based on the historical earth surface coverage data of the target area may be achieved by:
according to the earth surface coverage data, calculating a ground class transfer matrix of the target area; and respectively obtaining target area index data based on the ground class transfer matrix.
The earth class transfer matrix (Land Use Transition Matrix) is a matrix used to describe transfer relationships between earth surface coverage types. It records the probability or proportion of transitioning from one surface coverage type to another surface coverage type over a particular period of time.
For example, the target area is divided into m research units, the earth coverage data is taken as earth coverage data of many years, and the annual earth class transfer matrix of the target area can be expressed as follows:
wherein U is i Represents the ith unit, T j Representing the j-th transition place class.
Therefore, the embodiment of the application can further analyze and obtain index data reflecting the spatial diversity of the regional class transfer of the target area based on the method for calculating the regional class transfer matrix of the years based on the historical earth surface coverage data of the years. Based on historical data, the land use planning and prediction are performed in the future by using a land class transfer matrix, so that land resources are better managed and protected.
Referring to fig. 2, fig. 2 is a flow chart for generating a characteristic curve according to an embodiment of the present application; in an alternative implementation of the embodiments of the present application, step 200: based on the index data, obtaining the characteristic curve may include the steps of:
step 210: the index data are arranged in reverse order, and each index data is integrated to obtain a plurality of integrated data.
In the step S210, taking the data sample as the earth surface coverage data for many years as an example, the obtained index data may be index data between years. It should be understood that the integration operation in this application refers to the accumulation/summation operation; for example, for the n-th index data after the reverse arrangement, the integration operation performed on the n-th index data is to sum the 1 st to n-th index data, resulting in integration data about the n-th index data.
Step S220: and performing curve fitting on the plurality of integral data to obtain a characteristic curve.
In the above step S220, the obtained index data are arranged in reverse order, and then each index data is integrated to obtain a plurality of integrated data. Further, fitting a curve, such as a logarithmic function curve or an inverse proportion function curve, results in a characteristic curve.
As can be seen from fig. 2, by integrating the index data, a plurality of integrated data can be obtained; the method can reflect the cumulative effect in the surface coverage type transfer process, so that the characteristics of the surface coverage type transfer can be extracted and characterized better. And curve fitting is carried out on the plurality of integral data, so that a proper mathematical function can be found to describe the characteristic curve of the surface coverage type transfer, and the dynamic change rule of the surface coverage type transfer is better expressed.
Further, in an optional implementation manner of the embodiment of the present application, the index data includes land class transfer number data and land class transfer area data; the characteristic curves comprise a ground class transfer quantity characteristic curve and a ground class transfer area characteristic curve. Step S300: according to the index data corresponding to the part of the characteristic curve, the slope change rate of which does not exceed the screening threshold value, the target index data is obtained by the following modes:
and acquiring the intersection of the ground class transfer quantity data corresponding to the part of the ground class transfer quantity characteristic curve, the change rate of which does not exceed the screening threshold value, and the ground class transfer area data corresponding to the part of the ground class transfer area characteristic curve, the change rate of which does not exceed the screening threshold value, so as to obtain target index data.
If the index data includes the land class transfer number data and the land class transfer area data, then, taking a certain statistical unit as an example, the land class transfer number is as shown in table 1:
TABLE 1
In Table 1, U i Represents a specific statistical unit, T1 represents the last phase, and T2 representsNext time phase, C 1 To C n Representing N different land cover types, N ij C representing phase from T1 i C of phase to T2 j Is a number of spots.
The ground class transfer areas are shown in table 2:
TABLE 2
In Table 2, T1 represents the previous phase, T2 represents the next phase, C 1 To C n Represents n different land cover types, A ij C representing phase from T1 i C of phase to T2 j Is a pattern spot area of (c).
Thus, the land class transfer number data and the land class transfer area data are obtained, and the land class distribution arrangement of the above land class transfer number data and land class transfer area data is changed according to the operation unit and the geographical grid, respectively, as shown in table 3:
TABLE 3 Table 3
| Statistics unit | Type of ground class transfer | Quantity of | Area of |
| U i | C m C n | N mn | A mn |
| … | … | … | … |
First behavior statistics unit U in Table 3 i In land type C m Becomes land type C n Is of the number N mn Area A mn The method comprises the steps of carrying out a first treatment on the surface of the The type of the ground class transfer of the statistical unit and the corresponding number and area can be visually seen through the table 3, which is beneficial to the use of the data required later.
Further, after integrating the ground class transfer quantity data and the ground class transfer area data respectively, fitting a curve to obtain a ground class transfer quantity characteristic curve and a ground class transfer area characteristic curve, where the ground class transfer quantity characteristic curve and the ground class transfer area characteristic curve may be logarithmic functions such as y=a× lnx +b, such as inverse proportional functions of y=x/(a×x+b).
According to the above formulaCalculating the slope change rate of the curve; wherein Δk represents the slope change rate, Δy' i+1 ，Δy‘ i Indicating the slope change of points i+1 and i, y' i+1 、y‘ i 、y‘ i-1 The slopes of points i+1, i-1 are shown.
After the slope change rates of the ground class transfer quantity characteristic curve and the ground class transfer area characteristic curve are obtained, the ground class transfer type is screened by considering the two aspects of the ground class area and the ground class quantity and judging whether the slope change rate exceeds a screening threshold value.
Specifically, the ground class transfer type corresponding to the data which does not exceed the screening threshold value is screened out according to the ground class transfer quantity characteristic curve, and then the ground class transfer type corresponding to the data which does not exceed the screening threshold value is screened out according to the ground class transfer area characteristic curve; further, an intersection of two ground class transfer types is obtained, and the ground class transfer type corresponding to the target index data corresponding to the intersection is obtained.
Illustratively, the screening threshold may be set to 0.99, i.e., the selection curve changes slowly, approaching the data prior to the level.
For example, from the aspect of the land area, the type of land change obtained by screening is [ "0391_0120", "0311_0120", "0120_0131" ], wherein "0391_0120" represents the land with the code 0391 of the previous year, and the next year becomes the land with the code 0120; from the aspect of the number of the ground classes, the ground class change types obtained by screening are [ "0311_0120", "0391_0120", "0120_0521" ], and the ground class change types after intersection are [ "0391_0120", "0311_0120", "0120_0131", "0120_0521" ].
Therefore, according to the method for determining the space diversity of the earth surface coverage earth type transfer, which is provided by the embodiment of the application, the important earth type transfer types are more finely screened by considering two aspects of the earth type transfer quantity and the earth type transfer area, so that the space diversity condition of the earth type transfer can be more accurately determined.
Referring to fig. 3, fig. 3 is a flowchart of spatial diversity acquisition provided in an embodiment of the present application; in an alternative implementation of the embodiment of the present application, step S400: according to the target index data, determining the spatial diversity of the ground class transfer of the target area can be realized by the following steps:
step S410: and carrying out cluster analysis on the target index data to obtain a cluster analysis result.
In the step S410, cluster analysis is performed on the screened target index data to obtain a cluster analysis result; according to the screening result of the target index data, the national data are split into corresponding subsets, and the land area and the land quantity are clustered respectively.
Step S420: and determining the spatial diversity of the ground class transfer of the target area according to the clustering analysis result.
In the step S420, further, the data subset clustering result is subjected to the group average area statistics and the clustering heat analysis, and the representative group is selected for analysis, so as to determine the spatial diversity of the group transfer of the target area.
In an alternative embodiment, analyzing the spatial diversity of the class shifts of the target region based on the cluster analysis results, comprising: and respectively carrying out heat analysis on the ground class transfer quantity and the ground class transfer area according to the clustering analysis result.
For example, four subsets of data are obtained by splitting according to the variable clustering result, and if the subset 1 is one of the subsets, the clustering result of the subset 1 is based on the formula: ground class average area = ground class area/ground class number, the ground class average area statistics are calculated.
For example, the clustering heat analysis is performed on the clustering result of the subset 1, and it is known from the area average statistics and heat map of the transfer area of each region that the area value of 0391_0120 is highest in the 3 region, and is mainly distributed in the northeast, north and northwest regions of China, 0120_0131 is mainly distributed in the 4 and 5 regions, 0120_0311 is mainly distributed in the 1 and 5 regions, 0120_0391 is mainly distributed in the 3 and 5 regions, X0311_0120 is mainly distributed in the 3 and 5 regions, X0311_0391 is mainly distributed in the 1 and 5 regions, and X0391_0120 is mainly distributed in the 3 and 5 regions.
Further, 0391_0120 is selected as a representative ground class in the cluster nodes of the subset 1 for analysis, and the areas and the numbers of the transferred ground classes are analyzed by adopting a natural breakpoint method and the like, so that the distribution of the areas and the numbers is consistent, and the areas and the numbers of the ground classes are more in areas such as northeast China, shandong China, inner Mongolia and Yunnan.
That is, based on the cluster heat map, the area average statistics of the transfer land class of each area and the heat map, the situation that the land class transfer types corresponding to the target index data are distributed nationwide can be obtained.
In some embodiments, the abnormal value of the transfer of the ground class of the target area can be obtained according to the cluster analysis result, and the number of times of the transfer of the ground class in the abnormal value result according to the operation unit statistics in each class can be determined.
Therefore, the embodiment of the application can reveal the spatial distribution condition of different types of the ground class transfer in the target area through cluster analysis and ground class average area statistics, and know the characteristics and differences of the ground class transfer in different areas. The method is not only suitable for the analysis of the earth class transfer of the target area, but also can split the nationwide data into corresponding subsets for the spatial analysis of the nationwide earth class transfer, and provides a wider research view angle.
In an optional implementation manner of the embodiment of the present application, before performing cluster analysis on the target index data to obtain a cluster analysis result, the method further includes: statistical analysis and polymerizability analysis are performed on the target index data to evaluate the data quality of the target index data.
Statistical analysis may include calculating maximum, minimum, median, mean, standard deviation, etc., and when the data distribution has large variance and mean, normalization, etc., is required.
The analysis of the polymerizability is mainly to perform clustering evaluation analysis on data. And calculating the hopkins index of the data, and carrying out clustering evaluation. When the value of the hopkins statistic is <0.5, it is shown that the data is highly aggregated; further, for index data with good statistics and polymerizability, a kmeans method may be used for clustering.
Note that the hopkins index is a statistical index for evaluating the quality of clustering. It is used to determine whether the data set has a clustered structure, i.e., whether the data points tend to cluster together to form clusters, or whether the data points are evenly distributed throughout the data space. kmeans is a common clustering algorithm that divides data points into K clusters, each cluster being represented by an average of its internal data points (cluster center).
That is, the earth surface coverage earth type transfer spatial diversity method provided by the embodiment of the application evaluates the quality and reliability of target index data through statistical analysis and polymerizability analysis, and avoids inaccurate or unreliable clustering results caused by poor data quality.
Referring to fig. 4, fig. 4 is a schematic block diagram of a determining device for determining a ground type transfer spatial diversity of a ground surface coverage according to an embodiment of the present application; the apparatus 100 for determining the earth type transition spatial diversity of the earth surface coverage comprises: a characteristic curve generation module 110, a target data acquisition module 120, and an analysis result generation module 130.
The characteristic curve generating module 110 is configured to obtain index data of the earth class transfer according to historical earth surface coverage data of the target area; based on the index data, a characteristic curve is obtained.
The target data obtaining module 120 is configured to obtain target index data according to the index data corresponding to the portion where the slope change rate of the characteristic curve does not exceed the screening threshold;
the analysis result generation module 130 is configured to determine a spatial variance of the ground class transition of the target area according to the target index data.
In an alternative embodiment, the characteristic curve generating module 110 is specifically configured to calculate a ground class transfer matrix of the target area according to the ground surface coverage data in a process of obtaining the index data of the ground class transfer according to the historical ground surface coverage data of the target area; and respectively obtaining target area index data based on the ground class transfer matrix.
In an alternative embodiment, the characteristic curve generating module 110 is specifically configured to, in the process of obtaining the characteristic curve based on the index data: arranging index data in a reverse order, and integrating each index data to obtain a plurality of integrated data; performing curve fitting on the plurality of integral data to obtain a characteristic curve; wherein the characteristic curve comprises a logarithmic function curve or an inverse proportion function curve.
In an alternative embodiment, the index data includes ground class transfer quantity data and ground class transfer area data; the characteristic curves comprise a ground class transfer quantity characteristic curve and a ground class transfer area characteristic curve; the target data obtaining module 120 is specifically configured to, in a process of obtaining target index data according to the index data corresponding to a portion where the slope change rate of the characteristic curve does not exceed the screening threshold value: and acquiring the intersection of the ground class transfer quantity data corresponding to the part of the ground class transfer quantity characteristic curve, the change rate of which does not exceed the screening threshold value, and the ground class transfer area data corresponding to the part of the ground class transfer area characteristic curve, the change rate of which does not exceed the screening threshold value, so as to obtain target index data.
In an alternative embodiment, the analysis result generation module 130 is specifically configured to, in a process of analyzing the spatial diversity of the ground class transition of the target area according to the target index data: performing cluster analysis on the target index data to obtain a cluster analysis result; and determining the spatial difference of class transfer of the target area according to the clustering analysis result.
In an alternative embodiment, the determining device 100 for determining the earth type transfer spatial diversity covered by the earth surface is further configured to, before performing cluster analysis on the target index data to obtain a cluster analysis result: statistical analysis and polymerizability analysis are performed on the target index data to evaluate the data quality of the target index data.
In an alternative embodiment, the analysis result generation module 130 is specifically configured to, in a process of analyzing the spatial diversity of class transitions of the target area according to the cluster analysis result: and respectively carrying out heat analysis on the ground class transfer quantity and the ground class transfer area according to the clustering analysis result.
Referring to fig. 5, fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present application. An electronic device 300 provided in an embodiment of the present application includes: a processor 301 and a memory 302, the memory 302 storing machine-readable instructions executable by the processor 301, which when executed by the processor 301 perform the method as described above.
Based on the same inventive concept, embodiments of the present application also provide a computer readable storage medium, where a computer program instruction is stored, and when the computer program instruction is read and executed by a processor, the steps in any of the above implementations are performed.
The computer readable storage medium may be any of various media capable of storing program codes, such as random access Memory (Random Access Memory, RAM), read Only Memory (ROM), programmable Read Only Memory (Programmable Read-Only Memory, PROM), erasable Read Only Memory (Erasable Programmable Read-Only Memory, EPROM), electrically erasable Read Only Memory (Electric Erasable Programmable Read-Only Memory, EEPROM), and the like. The storage medium is used for storing a program, the processor executes the program after receiving an execution instruction, and the method executed by the electronic terminal defined by the process disclosed in any embodiment of the present invention may be applied to the processor or implemented by the processor.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other manners. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be other manners of division in actual implementation, and for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some communication interface, device or unit indirect coupling or communication connection, which may be in electrical, mechanical or other form.
Further, the units described as separate units may or may not be physically separate, and units displayed as units may or may not be physical units, may be located in one place, or may be distributed over a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
Furthermore, functional modules in various embodiments of the present application may be integrated together to form a single portion, or each module may exist alone, or two or more modules may be integrated to form a single portion.
Alternatively, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed on a computer, produces a flow or function in accordance with embodiments of the present invention, in whole or in part.
The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in or transmitted from one computer-readable storage medium to another, for example, by wired (e.g., coaxial cable, optical fiber, digital Subscriber Line (DSL)), or wireless (e.g., infrared, wireless, microwave, etc.).
In this document, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising … …" does not exclude the presence of other like elements in a process, method, article or apparatus that comprises the element.
The foregoing is merely exemplary embodiments of the present application and is not intended to limit the scope of the present application, and various modifications and variations may be suggested to one skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.
Claims (7)
1. A method for determining a spatial diversity of earth transfer of earth surface coverage, the method comprising:
obtaining index data of the earth class transfer according to historical earth surface coverage data of the target area;
arranging the index data in a reverse order, and integrating each index data to obtain a plurality of integrated data; performing curve fitting on a plurality of integral data to obtain a characteristic curve; the characteristic curves comprise logarithmic function curves or inverse proportion function curves, and the characteristic curves comprise ground class transfer quantity characteristic curves and ground class transfer area characteristic curves; the index data comprises ground class transfer quantity data and ground class transfer area data;
acquiring the land class transfer quantity data corresponding to the part of the land class transfer quantity characteristic curve, the change rate of which does not exceed the screening threshold value, and the intersection of the land class transfer area data corresponding to the part of the land class transfer area characteristic curve, the change rate of which does not exceed the screening threshold value, so as to obtain target index data;
performing cluster analysis on the target index data to obtain a cluster analysis result; and determining the spatial difference of the ground class transfer of the target area according to the clustering analysis result.
2. The method of claim 1, wherein obtaining the index data for the earth class transfer based on historical earth surface coverage data for the target area comprises:
calculating a ground class transfer matrix of the target area according to the ground surface coverage data;
and respectively obtaining the index data of the target area based on the ground class transfer matrix.
3. The method according to claim 1, wherein before performing cluster analysis on the target index data to obtain a cluster analysis result, the method further comprises:
and carrying out statistical analysis and polymerization analysis on the target index data to evaluate the data quality of the target index data.
4. The method of claim 1, wherein analyzing the spatial diversity of the geoclass transfer of the target region based on the cluster analysis results comprises:
and respectively carrying out heat analysis on the ground class transfer quantity and the ground class transfer area according to the clustering analysis result.
5. A device for determining a spatially diverse ground-based transition of a surface coverage, the device comprising: the system comprises a characteristic curve generation module, a target data acquisition module and an analysis result generation module;
the characteristic curve generation module is used for obtaining index data of the earth class transfer according to historical earth surface coverage data of the target area; arranging the index data in a reverse order, and integrating each index data to obtain a plurality of integrated data; performing curve fitting on a plurality of integral data to obtain a characteristic curve; the characteristic curves comprise logarithmic function curves or inverse proportion function curves, and the characteristic curves comprise ground class transfer quantity characteristic curves and ground class transfer area characteristic curves; the index data comprises ground class transfer quantity data and ground class transfer area data;
the target data acquisition module is used for acquiring the land class transfer quantity data corresponding to the part of the land class transfer quantity characteristic curve, the change rate of which is not more than the screening threshold value, and the intersection of the land class transfer area data corresponding to the part of the land class transfer area characteristic curve, the change rate of which is not more than the screening threshold value, so as to obtain target index data;
the analysis result generation module is used for carrying out cluster analysis on the target index data to obtain a cluster analysis result; and determining the spatial difference of the ground class transfer of the target area according to the clustering analysis result.
6. An electronic device comprising a memory and a processor, the memory having stored therein program instructions which, when executed by the processor, perform the steps of the method of any of claims 1-4.
7. A computer readable storage medium, characterized in that the computer readable storage medium has stored therein computer program instructions which, when executed by a processor, perform the steps of the method according to any of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310967151.3A CN117095042B (en) | 2023-08-02 | 2023-08-02 | Method and device for determining earth type transfer space difference of earth surface coverage |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310967151.3A CN117095042B (en) | 2023-08-02 | 2023-08-02 | Method and device for determining earth type transfer space difference of earth surface coverage |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117095042A CN117095042A (en) | 2023-11-21 |
| CN117095042B true CN117095042B (en) | 2024-02-23 |
Family
ID=88776256
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310967151.3A Active CN117095042B (en) | 2023-08-02 | 2023-08-02 | Method and device for determining earth type transfer space difference of earth surface coverage |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117095042B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103049916A (en) * | 2013-01-18 | 2013-04-17 | 国家基础地理信息中心 | Method for detecting land cover changes based on spectral slope differences |
| CN108537441A (en) * | 2018-04-09 | 2018-09-14 | 中国地质大学（武汉） | Land Use Transition economic society Effect Evaluation and monitoring information integrated system |
| CN109739943A (en) * | 2018-12-14 | 2019-05-10 | 中国测绘科学研究院 | Change statistical processing methods towards natural resources vector ground mulching |
| CN110837875A (en) * | 2019-11-18 | 2020-02-25 | 国家基础地理信息中心 | Method and device for judging quality abnormity of earth surface coverage data |
| CN112750299A (en) * | 2019-10-29 | 2021-05-04 | 华为技术有限公司 | Traffic flow analysis method and device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102020122010B4 (en) * | 2020-08-24 | 2023-05-04 | Bareways GmbH | METHOD AND SYSTEM FOR DETERMINING A CONDITION OF A GEOGRAPHIC LINE |
-
2023
- 2023-08-02 CN CN202310967151.3A patent/CN117095042B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103049916A (en) * | 2013-01-18 | 2013-04-17 | 国家基础地理信息中心 | Method for detecting land cover changes based on spectral slope differences |
| CN108537441A (en) * | 2018-04-09 | 2018-09-14 | 中国地质大学（武汉） | Land Use Transition economic society Effect Evaluation and monitoring information integrated system |
| CN109739943A (en) * | 2018-12-14 | 2019-05-10 | 中国测绘科学研究院 | Change statistical processing methods towards natural resources vector ground mulching |
| CN112750299A (en) * | 2019-10-29 | 2021-05-04 | 华为技术有限公司 | Traffic flow analysis method and device |
| CN110837875A (en) * | 2019-11-18 | 2020-02-25 | 国家基础地理信息中心 | Method and device for judging quality abnormity of earth surface coverage data |
Non-Patent Citations (5)
| Title |
|---|
| Spatial Differentiation of Land Use and Landscape Pattern Changes in the Beijing-Tianjin-Hebei Area;Yafer Wang;《Aerospace Information Research Institute 》;1-15 * |
| Transforming the autocorrelation function of a time series to detect land cover change;B.P.Salmon等;《2016 IEEE International Geoscience and Remote Sensing Symposium》;5181-5184 * |
| 吉林省土地利用空间分异变化研究;于海霞等;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》(第2022年第4期期);C038-1474 * |
| 江西省土地利用变化及其对人类活动的响应;徐羽;钟业喜;冯兴华;徐丽婷;郑林;吴巍;;水土保持研究(第01期);181-186 * |
| 面向地理国情监测的地表覆盖变化信息提取方法;程滔;李广泳;陶舒;周旭;毕凯;;测绘地理信息(第03期);103-107 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117095042A (en) | 2023-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Woolf et al. | Statewide modeling of bobcat, Lynx rufus, habitat in Illinois, USA | |
| Betts et al. | Thresholds in songbird occurrence in relation to landscape structure | |
| Koch et al. | Revealing areas of high nature conservation importance in a seasonally dry tropical forest in Brazil: Combination of modelled plant diversity hot spots and threat patterns | |
| Potvin et al. | A Bayesian hierarchical modeling framework for correcting reporting bias in the US tornado database | |
| Bastin et al. | Separating grazing and rainfall effects at regional scale using remote sensing imagery: A dynamic reference-cover method | |
| Pressey et al. | Scheduling conservation action in production landscapes: priority areas in western New South Wales defined by irreplaceability and vulnerability to vegetation loss | |
| Shi et al. | Local spatial modeling of white-tailed deer distribution | |
| Mortuza et al. | Historical and future drought in Bangladesh using copula-based bivariate regional frequency analysis | |
| Vasconcelos et al. | Identifying priority areas for island endemics using genetic versus specific diversity–the case of terrestrial reptiles of the Cape Verde Islands | |
| Bashari et al. | Risk of fire occurrence in arid and semi-arid ecosystems of Iran: an investigation using Bayesian belief networks | |
| Hysa et al. | A GIS based method for indexing the broad-leaved forest surfaces by their wildfire ignition probability and wildfire spreading capacity | |
| Gonçalves et al. | Exploring the spatiotemporal dynamics of habitat suitability to improve conservation management of a vulnerable plant species | |
| CN112907113B (en) | Vegetation change cause identification method considering spatial correlation | |
| Zaniolo et al. | Automatic design of basin-specific drought indexes for highly regulated water systems | |
| Fernandes et al. | A novel downscaling approach to predict plant invasions and improve local conservation actions | |
| Masoudvaziri et al. | Impact of geophysical and anthropogenic factors on wildfire size: a spatiotemporal data-driven risk assessment approach using statistical learning | |
| Deepthi et al. | Performance assessment of general circulation models: application of compromise programming method and global performance indicator technique | |
| Yost | Probabilistic modeling and mapping of plant indicator species in a Northeast Oregon industrial forest, USA | |
| Griffin et al. | Ranking Mahalanobis distance models for predictions of occupancy from presence‐only data | |
| CN117095042B (en) | Method and device for determining earth type transfer space difference of earth surface coverage | |
| Lindley et al. | Inter-comparison of interpolated background nitrogen dioxide concentrations across Greater Manchester, UK | |
| Rosamond et al. | Representing a large region with few sites: The Quality Index approach for field studies | |
| CN115439287A (en) | Geological disaster risk evaluation method based on machine learning | |
| Croicu et al. | A fast spatial multiple imputation procedure for imprecise armed conflict events | |
| CN115345488A (en) | Wetland ecosystem service evaluation method, device, server and medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |