Method and device for historical reconstruction of spatial distribution of land use change including conversion information

Abstract

The present application relates to a method and device for historical reconstruction of the spatial distribution of land use changes including conversion information, which belongs to the field of land use spatial data processing technology, and solves the problem of deviation between the existing land use conversion area based on remote sensing reconstruction and the statistical area of ​​ground survey. By integrating the long-term series of historical land use reconstruction data of the study area, the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area is determined, and the proportion of the area that does not undergo type conversion from each year to the next year is calculated. The statistical area of ​​land use in each year is ranked according to the proportion of the total area of ​​the study area, and the proportion of the area transferred out and transferred in for each land use type is iteratively calculated to obtain the reconstruction data of the spatial distribution of land use including conversion information. The present application achieves the effect of providing reliable land use conversion areas while ensuring statistical reliability consistent with survey statistical data.

Description

Technical Field

[0001] The present application relates to the technical field of land spatial distribution data processing, and in particular to a method and device for historical reconstruction of spatial distribution of land use changes including conversion information. Background Art

[0002] Long-term land use spatial data sets are of great value in revealing the historical trend of land use changes, predicting future changes, supporting land use planning and ecological protection. However, the long-term land use spatial distribution data at the global and regional levels are relatively limited, and these data often lack information on the conversion between different land use types. For example, the "Global Environmental History Dataset" (HYDE dataset) or the land use spatial distribution dataset based on statistical data only provide the change area of ​​land use by type, and do not include information such as the conversion area between different land use types; the "Land Use Assimilation Dataset" (LUH2 dataset) assumes the priority of conversion between different land use types and sets the minimum area conversion rule to simulate the conversion area between artificial land use types and natural land use types. However, there is still a significant difference between the above information and the actual land use conversion area.

[0003] Therefore, for the historical reconstruction of the spatial distribution of land use changes, it is crucial to solve the data sources and quantitative methods for describing the spatial distribution of land use type conversion. With the promotion and application of large-scale remote sensing image data and inversion technology, it can provide more real and reliable information on the distribution of objects, and conveniently become an important basis for the historical reconstruction of the spatial distribution of land use changes. Land use data based on remote sensing means can provide the area of ​​mutual conversion of land use in different spatial locations. However, the coverage period of remote sensing data is limited, which restricts the application scope of historical reconstruction of the spatial distribution of land use changes based solely on remote sensing data. At the same time, in the actual data processing process, the uncertainty of remote sensing inversion technology itself and the limited high-resolution remote sensing data have limited the reliability of land use remote sensing mapping, resulting in a significant deviation between the land use conversion area based on remote sensing reconstruction and the statistical area of ​​ground survey. Summary of the invention

[0004] In view of the above analysis, an embodiment of the present invention aims to provide a method and device for historical reconstruction of the spatial distribution of land use changes including conversion information, so as to solve the problems of missing conversion information between different land use types in the implementation results of the prior art and obvious deviation between the land use conversion area reconstructed by remote sensing and the statistical area of ​​ground survey.

[0005] The present application embodiment provides a method for historical reconstruction of spatial distribution of land use changes including conversion information, comprising the steps of:

[0006] Obtaining a remote sensing probability matrix of land use conversion; wherein the remote sensing probability matrix includes the average proportion of the area of ​​each land use type conversion to the total area of ​​the study area;

[0007] Based on the remote sensing probability matrix, the proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type is calculated, and a calculation rule for the proportion of the area without type conversion is established;

[0008] According to the remote sensing probability matrix and the calculation rules, the long-term series of land use reconstruction historical data of the study area are integrated to determine the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and the proportion of the area that does not undergo type conversion from each year to the next year is calculated based on this;

[0009] The land use statistical area of ​​each year is ranked according to its proportion to the total area of ​​the study area, and the area ratio of each land use type transferred out and in is iteratively calculated to obtain the reconstructed data of the spatial distribution of land use containing conversion information.

[0010] The historical reconstruction method of the spatial distribution of land use change containing conversion information in the embodiment of the present application calculates the proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, establishes a calculation rule for the proportion of the area without type conversion, integrates the long-term series of land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rule for the proportion of the area without type conversion, determines the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculates the proportion of the area without type conversion from each year to the next year, sorts the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, and iteratively calculates the proportion of the area transferred out and transferred in for each land use type, and obtains the reconstruction data of the spatial distribution of land use containing conversion information. Based on this, it not only provides a reliable land use conversion area, but also ensures the statistical reliability consistent with the survey statistical data.

[0011] As one of the optional embodiments, the process of obtaining the remote sensing probability matrix of land use includes the steps of:

[0012] Obtain the land use data covering the study area year by year and reclassify the land use types;

[0013] Superimpose the land use data of adjacent years, calculate the area of ​​each land use type that has not been converted and the area of ​​conversion between different land use types, and establish the land use change conversion probability matrix for each year;

[0014] The average value of the land use change transfer matrix for each year was calculated to obtain the remote sensing probability matrix of land use conversion.

[0015] As one of the optional embodiments, the process of calculating the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix and establishing a calculation rule for the area ratio without type conversion is as follows:

[0016] The proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type is as follows;

[0017] ;

[0018] in, ;

[0019] Among them, k, i and j represent the numbers of land use types. It refers to the proportion of the area where the type of land use k is converted from a certain initial year to the next year to the total area of ​​the study area. represents the proportion of land use area of ​​type k in the initial year to the total area of ​​the study area; and is the average proportion of the area of ​​each land use type conversion in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types;

[0020] The remote sensing probability matrix is:

[0021]

[0022] Among them, A is the remote sensing probability matrix of land use conversion. Aij in the remote sensing probability matrix represents the average proportion of the area converted from type i land use to type j land use to the total area of ​​the study area. If i=j, it is the proportion of the area that has not been converted.

[0023] As one of the optional embodiments, the long-term land use reconstruction historical data of the study area are integrated according to the remote sensing probability matrix and the calculation rules to determine the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and the proportion of the area that does not undergo type conversion from each year to the next year is calculated as follows:

[0024] For any year, the land use type conversion calculation of the study area is performed, the statistical area of ​​each type of land use in the study area in year y and year y+1 is obtained, and its proportion of the total area of ​​the study area is calculated, which are marked as and ; Where k represents the number of the land use type;

[0025] Calculate the proportion of area where type k land use does not change from year y to year y+1 ,as follows:

[0026] ;

[0027] in, It represents the ratio of the theoretically unchanged area of ​​each land use type from year y to year y+1 to the average area of ​​the corresponding land use type.

[0028] As one of the optional embodiments, the process of sorting the land use statistical area of ​​each year by the proportion of the total area of ​​the study area, iteratively calculating the area ratio of each land use type transferred out and transferred in, and obtaining the reconstructed data of the land use spatial distribution containing the conversion information includes the steps of:

[0029] The statistical area of ​​land use in each year is sorted from small to large according to the proportion of the total area of ​​the study area, and the area ratio of each land use type is iteratively calculated to obtain the statistical probability matrix of land use conversion from each year to the next year;

[0030] Based on the statistical probability matrix, summing the proportion of each land use type that undergoes type conversion and is converted to other land use types from each year to the next year, and calculating the proportion of each land use type in the total sum;

[0031] By multiplying the area of ​​each land use type by the proportion of the total sum, the occurrence type conversion of each land use type and the area converted to other land use types are calculated, thus realizing the historical reconstruction of the spatial distribution of land use changes including conversion information.

[0032] As one of the optional embodiments, the process of iteratively calculating the area ratio of each land use type transferred out and transferred in is as follows:

[0033] For the smallest land use type The conversion area is calculated, and the ratio of the transfer area from year y to year y+1 is as follows:

[0034] ;

[0035] in, It means that from year y to year y+1 is The proportion of the area converted to type j land use in the total area of ​​the study area, where j is the land use type number; express The proportion of type land use in the total area of ​​the study area in year y; Represents from year y to year y+1 The proportion of the area with no type conversion in land use to the total area of ​​the study area; Represents the remote sensing probability matrix The average proportion of the area converted to type j land use to the total area of ​​the study area; n is the number of land use types.

[0036] right The proportion of land use type transferred from year y to year y+1 is calculated as follows:

[0037] ;

[0038] in, Indicates that from year y to year y+1, type i land use is converted to The proportion of the area of ​​the type of land use to the total area of ​​the study area, i is the land use type number; express The proportion of type land use in the total area of ​​the study area in year y + 1; Indicates that the remote sensing probability matrix is ​​converted from year y to year y+1. The average proportion of the area of ​​the type of land use to the total area of ​​the study area; n is the number of land use types.

[0039] As one of the optional embodiments, the process of iteratively calculating the area ratio of each land use type transferred out and transferred in is as follows:

[0040] The conversion area of ​​the second smallest land use type m is calculated, and its conversion area ratio from year y to year y+1 is as follows:

[0041] ;

[0042] in, It represents the proportion of the area converted from type m land use to type j land use from year y to year y+1 to the total area of ​​the study area, where j is the land use type number; It represents the proportion of land use type m in the total area of ​​the study area in year y; It represents the proportion of the area with no type conversion of land use type m from year y to year y+1 to the total area of ​​the study area; Indicates the land use conversion of type m from year y to year y+1 The proportion of the area of ​​the type to the total area of ​​the study area; It represents the average proportion of the area of ​​land use type m converted to land use type j from year y to year y+1 in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types.

[0043] The proportion of land use area transferred from year y to year y+1 of type m is calculated as follows:

[0044] ;

[0045] in, It represents the proportion of the area converted from type i land use to type m land use from year y to year y+1 to the total area of ​​the study area, where m is the land use type number; It represents the proportion of land use type m in the total area of ​​the study area in year y+1; Represents from year y to year y+1 The proportion of area converted to type m land use in the total area of ​​the study area; It represents the average proportion of the area of ​​land use type i converted to land use type m from year y to year y+1 in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types.

[0046] The present application also provides a device for reconstructing the history of spatial distribution of land use changes including conversion information, including:

[0047] A basic data module is used to obtain a remote sensing probability matrix of land use conversion; wherein the remote sensing probability matrix includes the average proportion of the area of ​​each land use type that produces type conversion to the total area of ​​the study area;

[0048] The first conversion module is used to calculate the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, and establish a calculation rule for the area ratio without type conversion;

[0049] The second conversion module is used to integrate the long-term land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rules, determine the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculate the proportion of the area that does not undergo type conversion from each year to the next year;

[0050] The data iteration module is used to sort the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, and iteratively calculate the area ratio of each land use type transferred out and transferred in, so as to obtain the reconstructed data of the spatial distribution of land use containing conversion information.

[0051] The historical reconstruction device of the spatial distribution of land use changes containing conversion information in the embodiment of the present application calculates the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, establishes a calculation rule for the area ratio of non-type conversion, integrates the long-term series of land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rule for the area ratio of non-type conversion, determines the ratio of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculates the area ratio of non-type conversion from each year to the next year, sorts the statistical area of ​​land use in each year according to the ratio of the total area of ​​the study area, and iteratively calculates the area ratio of each land use type transferred out and transferred in, and obtains the reconstruction data of the spatial distribution of land use containing conversion information. Based on this, it not only provides a reliable land use conversion area, but also ensures the statistical reliability consistent with the survey statistical data.

[0052] At least one embodiment of the present application further provides a data control device, including:

[0053] one or more memories non-transitorily storing computer-executable instructions;

[0054] One or more processors are configured to execute computer executable instructions, wherein the computer executable instructions, when executed by the one or more processors, implement the method for historical reconstruction of the spatial distribution of land use changes containing conversion information according to any embodiment of the present application.

[0055] The above data control device calculates the proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, establishes the calculation rules for the proportion of the area without type conversion, integrates the long-term land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rules for the proportion of the area without type conversion, determines the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculates the proportion of the area without type conversion from each year to the next year, sorts the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, and iteratively calculates the proportion of the area transferred out and transferred in for each land use type, and obtains the reconstruction data of the land use spatial distribution containing conversion information. Based on this, it not only provides a reliable land use conversion area, but also ensures the statistical reliability consistent with the survey statistical data.

[0056] At least one embodiment of the present application further provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer-executable instructions, which, when executed by a processor, implement a method for historical reconstruction of the spatial distribution of land use changes including conversion information according to any embodiment of the present application.

[0057] The above non-transient computer-readable storage medium calculates the proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, establishes the calculation rules for the proportion of the area without type conversion, integrates the long-term land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rules for the proportion of the area without type conversion, determines the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculates the proportion of the area without type conversion from each year to the next year based on this, sorts the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, and iterates the proportion of the area transferred out and transferred in for each land use type to obtain the reconstruction data of the land use spatial distribution containing conversion information. Based on this, it not only provides reliable land use conversion areas, but also ensures the statistical reliability consistent with the survey statistical data. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 A flow chart of a method for historical reconstruction of spatial distribution of land use changes including conversion information provided by the present invention;

[0059] Figure 2 A flow chart of a method for historical reconstruction of spatial distribution of land use changes including conversion information provided by a preferred embodiment of the present invention;

[0060] Figure 3 A module structure diagram of a device for reconstructing the history of spatial distribution of land use changes including conversion information according to an embodiment of the present invention;

[0061] Figure 4 A schematic block diagram of a data control device provided by the present invention;

[0062] Figure 5 A schematic diagram of a non-transitory computer-readable storage medium provided by the present invention. DETAILED DESCRIPTION

[0063] In order to make the purpose, technical scheme and advantages of the embodiments of the present application clearer, the technical scheme of the embodiments of the present application will be clearly and completely described below in conjunction with the drawings of the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all of the embodiments. Based on the described embodiments of the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0064] Unless otherwise defined, the technical terms or scientific terms used in this application should be understood by people with ordinary skills in the field to which this application belongs. "First", "second" and similar words used in this application do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word cover the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0065] In order to keep the following description of the embodiments of the present application clear and concise, the present application omits detailed descriptions of some known functions and known components.

[0066] An embodiment of the present application provides a method for historical reconstruction of the spatial distribution of land use changes including conversion information.

[0067] Figure 1 FIG. 1 is a flow chart of a method for historical reconstruction of spatial distribution of land use changes including conversion information according to an embodiment of the present invention. Figure 1 As shown, a method for historical reconstruction of spatial distribution of land use changes including conversion information in one embodiment includes steps S100 to S103:

[0068] S100, obtaining a remote sensing probability matrix of land use conversion; wherein the remote sensing probability matrix includes an average proportion of the area of ​​each land use type conversion to the total area of ​​the study area;

[0069] S101, calculating the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, and establishing a calculation rule for the ratio of the area that does not undergo type conversion;

[0070] S102, integrating the long-term land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rule, determining the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and using this to calculate the proportion of the area that does not undergo type conversion from each year to the next year;

[0071] S103, sorting the land use statistics of each year by the proportion of the total area of ​​the study area, and iteratively calculating the proportion of the area transferred out and transferred in for each land use type, to obtain the reconstructed data of the land use spatial distribution containing the conversion information.

[0072] The spatial distribution of land use in the embodiment of the present application is based on the study area as a spatial reference. The land area in the study area is collected as remote sensing data, i.e., land use data. Land use data can be divided into different land use types. Any area in the study area changes from one land use type to another land use type, i.e., land use type conversion.

[0073] As one of the preferred embodiments, Figure 2 FIG. 1 is a flow chart of a method for historical reconstruction of spatial distribution of land use changes including conversion information according to a preferred embodiment of the present invention. Figure 2 As shown, the process of obtaining the remote sensing probability matrix of land use in step S100 includes steps S200 to S202:

[0074] S200, obtain the land use data covering the study area year by year and reclassify the land use types;

[0075] S201, superimposing the land use data of adjacent years, calculating the area of ​​each land use type that has not been converted and the area of ​​conversion between different land use types, and establishing a land use change conversion probability matrix for each year;

[0076] S202, calculating the average value of the land use change transfer matrix of each year to obtain the remote sensing probability matrix of land use conversion.

[0077] Preferably, the year-by-year land use data covering the study area is obtained, for example, the year-by-year land use spatial distribution data of CLCD from 1990 to 2020 is selected; the land use types are reclassified, and this embodiment divides the land use types into 6 types: cultivated land, forest land, grassland, water area, construction land and unused land.

[0078] Superimpose the land use data of adjacent years, calculate the area of ​​each land use type that has not been converted and the area that has been converted between different land use types, and based on this, establish a year-by-year land use change conversion probability matrix. This example is the year-by-year land use change conversion probability matrix from 1991 to 2020.

[0079] In this step, the average value of the land use change transfer matrix of all years is calculated, which is the remote sensing probability matrix A of land use conversion, as follows:

[0080]

[0081] Among them, A is the remote sensing probability matrix of land use conversion. Aij in the remote sensing probability matrix represents the average proportion of the area converted from type i land use to type j land use to the total area of ​​the study area. If i=j, it is the proportion of the area that has not been converted.

[0082] For example, i and j are both 6, and the remote sensing probability matrix of land use conversion is as follows:

[0083]

[0084] Among them, 1-6 represent the six types of cultivated land, forest land, grassland, water area, construction land and unused land, such as A 25 It represents the proportion of forest land converted to construction land to the total area of ​​the study area. If i=j, it is the proportion of area that has not been converted. 11 It indicates the proportion of cultivated land that has not been converted to the total area of ​​the study area.

[0085] According to the determination of the remote sensing probability matrix, the process of calculating the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix in step S101, and establishing the calculation rule of the area ratio without type conversion, is as follows:

[0086] The proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type is as follows:

[0087] ;

[0088] in, ;

[0089] Among them, k, i and j represent the numbers of land use types. It refers to the proportion of the area where the type of land use k is converted from a certain initial year to the next year to the total area of ​​the study area. represents the proportion of land use area of ​​type k in the initial year to the total area of ​​the study area; and is the average proportion of the area of ​​each land use type conversion in the remote sensing probability matrix to the total area of ​​the study area. n is the number of land use types.

[0090] In the embodiment of the present application, k, i and j can take values ​​of 1-6, representing 6 types of land use, namely, cultivated land, forest land, grassland, water area, construction land and unused land. n is the number of land use types, and n is 6 in the present embodiment.

[0091] Furthermore, a long-term series of historical land use reconstruction data covering the study area is obtained, and the statistical area of ​​each land use type in each year in the study area is counted, and the spatial resolution of the historical land use reconstruction data is used as the mapping unit size for land use conversion reconstruction in the present invention.

[0092] In step S102, the long-term land use reconstruction historical data of the study area are integrated according to the remote sensing probability matrix and the calculation rules to determine the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and the proportion of the area that does not undergo type conversion from each year to the next year is calculated as follows:

[0093] For any year, the land use type conversion calculation of the study area is performed, the statistical area of ​​each type of land use in the study area in year y and year y+1 is obtained, and its proportion of the total area of ​​the study area is calculated, which are marked as and ; Where k represents the number of the land use type;

[0094] In the embodiment of the present application, k can take values ​​of 1-6, representing six types of land, namely, cultivated land, forest land, grassland, water area, construction land and unused land.

[0095] Calculate the proportion of area where type k land use does not change from year y to year y+1 ,as follows:

[0096] ;

[0097] in, It indicates the proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type.

[0098] Preferably, if Less than , then assign it to .

[0099] As one of the preferred embodiments, Figure 2 As shown, the process of sorting the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, iteratively calculating the area ratio of each land use type transferred out and transferred in, and obtaining the reconstructed data of the land use spatial distribution containing the conversion information, includes steps S300 to S302:

[0100] S300, sorting the statistical area of ​​land use in each year by the proportion of the statistical area of ​​the land use in the total area of ​​the study area from small to large, iteratively calculating the proportion of the area transferred out and transferred in of each land use type, and obtaining a statistical probability matrix of land use conversion from each year to the next year;

[0101] S301, based on the statistical probability matrix, summing the proportion of each land use type that undergoes type conversion and is converted to other land use types from each year to the next year, and calculating the proportion of each land use type in the total sum;

[0102] S302, multiplying the area of ​​each land use type by the proportion of the total sum, calculating the occurrence type conversion of each land use type and the area converted to other land use types, and realizing the historical reconstruction of the spatial distribution of land use changes including conversion information.

[0103] Preferably, in step S301, the process of iteratively calculating the area ratio of each land use type transferred out and transferred in is as follows:

[0104] The conversion area of ​​the smallest land use type l is calculated, and the proportion of the conversion area from year y to year y+1 is as follows:

[0105] ;

[0106] in, It means that from year y to year y+1 is The proportion of the area converted to type j land use in the total area of ​​the study area, where j is the land use type number; express The proportion of type land use in the total area of ​​the study area in year y; Represents from year y to year y+1 The proportion of the area with no type conversion in land use to the total area of ​​the study area; Represents the remote sensing probability matrix The average proportion of the area converted to type j land use to the total area of ​​the study area; n is the number of land use types;

[0107] right The proportion of land use type transferred from year y to year y+1 is calculated as follows:

[0108] ;

[0109] in, Indicates that from year y to year y+1, type i land use is converted to The proportion of the area of ​​the type of land use to the total area of ​​the study area, i is the land use type number; express The proportion of type land use in the total area of ​​the study area in year y + 1; Indicates that the remote sensing probability matrix is ​​converted from year y to year y+1. The average proportion of the area of ​​the type of land use to the total area of ​​the study area; n is the number of land use types.

[0110] Preferably, j and i can take values ​​of 1-6, representing the six types of land, namely, cultivated land, forest land, grassland, water area, construction land and unused land. The number of the type land use.

[0111] Preferably, in step S301, the process of iteratively calculating the area ratio of each land use type transferred out and transferred in is as follows:

[0112] The conversion area of ​​the second smallest land use type m is calculated, and its conversion area ratio from year y to year y+1 is as follows:

[0113] ;

[0114] in, It represents the proportion of the area converted from type m land use to type j land use from year y to year y+1 to the total area of ​​the study area, where j is the land use type number; It represents the proportion of land use type m in the total area of ​​the study area in year y; It represents the proportion of area where land use type m does not undergo type conversion from year y to year y+1; Indicates the land use conversion of type m from year y to year y+1 The proportion of the area of ​​the type to the total area of ​​the study area; It represents the average proportion of the area of ​​land use type m converted to land use type j from year y to year y+1 in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types;

[0115] The proportion of land use area transferred from year y to year y+1 of type m is calculated as follows:

[0116] ;

[0117] in, It represents the proportion of the area converted from type i land use to type m land use from year y to year y+1 to the total area of ​​the study area, where m is the land use type number; It represents the proportion of land use type m in the total area of ​​the study area in year y+1; Represents from year y to year y+1 The proportion of area converted to type m land use in the total area of ​​the study area; It represents the average proportion of the area of ​​land use type i converted to land use type m from year y to year y+1 in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types.

[0118] Preferably, j and m can take values ​​of 1-6, representing the six types of land, namely, cultivated land, forest land, grassland, water area, construction land and unused land. type and m type of land use number.

[0119] The above iterative process is repeated in the order of land use area from small to large in year y, wherein the land use conversion area calculated in the previous process is no longer included in the calculation, and only the remaining land use conversion area is calculated; when the iterative process is executed to the last two types, if there is a situation where the conversion area is negative to ensure that the proportion of each converted land use area is the same as the statistical area proportion, then according to the principle of minimum conversion area, the proportion of the area that has not been converted is adjusted; finally, it is ensured that any land use type from year y to year y+1 through the transfer in and out of different land use types is the same as the land use statistical area in year y and year y+1, and the statistical probability matrix of land use conversion from year y to year y+1 is obtained.

[0120] Based on the statistical probability matrix of land use conversion, the proportion of areas converted from year y and year y+1 and converted to other land use types is summed for each land use type, and the proportion of each type in the total sum is calculated; on this basis, the area of ​​each land use type in each mapping unit is calculated; the above two are multiplied to calculate the area of ​​conversion and conversion to other land use types for each land use type in each mapping unit, and finally the historical reconstruction of the spatial distribution of land use including conversion information is achieved.

[0121] Based on this, the embodiment of the present application effectively integrates the statistical information of ground surveys and the spatial information provided by remote sensing inversion, and establishes a land use history reconstruction method containing conversion information by calculating the land use conversion matrix based on remote sensing data and designing a cyclic iterative process for year-by-year land use conversion and mutual conversion areas. Compared with traditional methods, the embodiment of the present application can not only provide reliable land use conversion areas, but also achieve the statistical reliability of the land use reconstruction area calculated based on the conversion matrix and the land use statistical area based on the survey.

[0122] The historical reconstruction method of the spatial distribution of land use changes containing conversion information in the embodiment of the present application calculates the proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, establishes a calculation rule for the proportion of the area without type conversion, and integrates the long-term series of land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rule, determines the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculates the proportion of the area that does not undergo type conversion from each year to the next year based on this, sorts the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, and iteratively calculates the proportion of the area transferred out and transferred in for each land use type, and obtains the reconstruction data of the spatial distribution of land use containing conversion information. Based on this, it not only provides a reliable land use conversion area, but also ensures the statistical reliability consistent with the survey statistical data.

[0123] The embodiment of the present application also provides a device for historically reconstructing the spatial distribution of land use changes including conversion information.

[0124] Figure 3 FIG. 1 is a block diagram of a device for reconstructing the spatial distribution of land use changes including conversion information according to an embodiment of the present invention. Figure 3 As shown, a historical reconstruction device for spatial distribution of land use changes including conversion information in one embodiment includes:

[0125] The basic data module 100 is used to obtain a remote sensing probability matrix of land use conversion; wherein the remote sensing probability matrix includes the average proportion of the area of ​​each land use type that produces type conversion to the total area of ​​the study area;

[0126] The first conversion module 101 is used to calculate the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, and establish a calculation rule for the area ratio without type conversion;

[0127] The second conversion module 102 is used to integrate the long-term land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rule, determine the proportion of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculate the proportion of the area that does not undergo type conversion from each year to the next year;

[0128] The data iteration module 103 is used to sort the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, and iteratively calculate the area ratio of each land use type transferred out and transferred in, so as to obtain the reconstructed data of the spatial distribution of land use containing conversion information.

[0129] The historical reconstruction device of the spatial distribution of land use change containing conversion information in the embodiment of the present application calculates the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix, establishes a calculation rule for the area ratio without type conversion, integrates the long-term series of land use reconstruction historical data of the study area according to the remote sensing probability matrix and the calculation rule, determines the ratio of the statistical area of ​​land use in each year to the total area of ​​the study area, and calculates the area ratio of each year to the next year without type conversion, sorts the statistical area of ​​land use in each year according to the ratio of the total area of ​​the study area, and iteratively calculates the area ratio of each land use type transferred out and transferred in, and obtains the reconstruction data of the spatial distribution of land use containing conversion information. Based on this, it not only provides a reliable land use conversion area, but also ensures the statistical reliability consistent with the survey statistical data.

[0130] At least one embodiment of the present application also provides a data control device. Figure 4 A schematic block diagram of a data control device provided in at least one embodiment of the present application. Figure 4 As shown, the data control device 20 may include one or more memories 200 and one or more processors 201. The memory 200 is used to store computer executable instructions non-transiently; the processor 201 is used to run the computer executable instructions, and when the computer executable instructions are run by the processor 201, the processor 201 may execute one or more steps in the historical reconstruction method of the spatial distribution of land use changes containing conversion information according to any embodiment of the present application.

[0131] The specific implementation and related explanation of each step of the method for historical reconstruction of spatial distribution of land use changes including conversion information can be found in the above-mentioned embodiment of the method for historical reconstruction of spatial distribution of land use changes including conversion information, and will not be repeated here. It should be noted that Figure 4 The components of the data control device 20 shown are merely exemplary and non-limiting. The data control device 20 may also have other components according to actual application requirements.

[0132] In one embodiment, the processor 201 and the memory 200 can communicate with each other directly or indirectly. For example, the processor 201 and the memory 200 can communicate via a network connection. The network may include a wireless network, a wired network, and / or any combination of a wireless network and a wired network, and the present application does not limit the type and function of the network. For another example, the processor 201 and the memory 200 can also communicate via a bus connection. The bus may be a peripheral component interconnect standard (PCI) bus or an extended industrial standard architecture (EISA) bus, etc. For example, the processor 201 and the memory 200 may be arranged at a remote data server end (cloud) or a distributed energy system end (local end), or may be arranged at a client end (for example, a mobile device such as a mobile phone). For example, the processor 201 may be a device having data processing capability and / or instruction execution capability such as a central processing unit (CPU), a tensor processor (TPU) or a graphics processor GPU, and may control other components in the data prediction device 20 to perform the desired functions. The central processing unit (CPU) may be an X86 or ARM architecture, etc.

[0133] In one embodiment, the memory 200 may include any combination of one or more computer program products, and the computer program product may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. Volatile memory may include, for example, random access memory (RAM) and / or cache memory (cache), etc. Non-volatile memory may include, for example, read-only memory (ROM), hard disk, erasable programmable read-only memory (EPROM), portable compact disk read-only memory (CD-ROM), USB memory, flash memory, etc. One or more computer executable instructions may be stored on the computer-readable storage medium, and the processor 201 may execute the computer executable instructions to implement various functions of the data prediction device 20. Various applications and various data, as well as various data used and / or generated by the application, etc. may also be stored in the memory 200.

[0134] It should be noted that the data control device 20 can achieve technical effects similar to the aforementioned method for historical reconstruction of spatial distribution of land use changes including conversion information, and the repeated parts will not be repeated.

[0135] At least one embodiment of the present application also provides a non-transitory computer-readable storage medium. Figure 5 A schematic diagram of a non-transitory computer-readable storage medium provided for at least one embodiment of the present application. Figure 5As shown, one or more computer executable instructions 301 may be non-transitory stored on a non-transitory computer readable storage medium 30. For example, when the computer executable instructions 301 are executed by a computer, the computer may execute one or more steps in the method for historical reconstruction of spatial distribution of land use changes including conversion information according to any embodiment of the present application.

[0136] In one embodiment, the non-transitory computer-readable storage medium 30 may be applied to the above-mentioned data control device 20 , for example, it may be the memory 200 in the data control device 20 .

[0137] In one embodiment, the description of the non-transitory computer-readable storage medium 30 may refer to the description of the memory 200 in the embodiment of the data control device 20, and the repeated parts will not be repeated.

[0138] It should be noted that the memory 200 stores different non-transiently stored computer executable instructions, and the data control device 20 corresponds to a firmware upgrade device. When the computer executable instructions are executed by the processor 201, the processor 201 can execute one or more steps in the historical reconstruction method of the spatial distribution of land use changes containing conversion information according to any embodiment of the present application.

[0139] There are a few points to note about this application:

[0140] (1) The drawings of the embodiments of the present application only relate to the structures related to the embodiments of the present application, and other structures may refer to the general design.

[0141] (2) For the sake of clarity, the thickness and size of layers or structures are exaggerated in the drawings used to describe the embodiments of the present invention. It is understood that when an element such as a layer, film, region, or substrate is referred to as being "on" or "under" another element, the element may be "directly" "on" or "under" the other element, or there may be intervening elements.

[0142] (3) In the absence of conflict, the embodiments of the present application and the features in the embodiments may be combined with each other to obtain new embodiments. The above are only specific implementation methods of the present application, but the protection scope of the present application is not limited thereto. The protection scope of the present application shall be based on the protection scope of the claims.

[0143] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0144] The above embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the patent application. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent application shall be subject to the attached claims.

Claims

1. A method for historical reconstruction of the spatial distribution of land use changes including conversion information, characterized in that Includes steps: Obtain remote sensing probability matrix of land use conversion; Based on the remote sensing probability matrix, the proportion of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type is calculated according to the following formula: ； ； Among them, k, i and j represent the numbers of land use types. It refers to the proportion of the area where the type of land use k is converted from a certain initial year to the next year to the total area of ​​the study area. represents the proportion of land use area of ​​type k in the initial year to the total area of ​​the study area; and is the average proportion of the area of ​​each land use type conversion in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types; For any year, the land use type conversion calculation of the study area is performed, the statistical area of ​​each type of land use in the study area in year y and year y+1 is obtained, and its proportion of the total area of ​​the study area is calculated, which are marked as and ; Where k represents the number of the land use type; The proportion of area where type k land use does not change from year y to year y+1 is calculated according to the following formula: : ； in, It represents the ratio of the theoretical non-conversion area of ​​each land use type from year y to year y+1 to the average area of ​​the corresponding land use type; The land use statistical area of ​​each year is ranked according to its proportion to the total area of ​​the study area, and the area ratio of each land use type transferred out and in is iteratively calculated to obtain the reconstructed data of the spatial distribution of land use containing conversion information.

2. The method for historical reconstruction of the spatial distribution of land use changes including conversion information according to claim 1, characterized in that: The process of obtaining the remote sensing probability matrix of land use includes the steps of: Obtain the land use data covering the study area year by year and reclassify the land use types; Superimpose the land use data of adjacent years, calculate the area of ​​each land use type that has not been converted and the area of ​​conversion between different land use types, and establish the land use change conversion probability matrix for each year; The average value of the land use change transfer matrix for each year was calculated to obtain the remote sensing probability matrix of land use conversion.

3. The method for historical reconstruction of spatial distribution of land use changes including conversion information according to claim 1, characterized in that: The process of sorting the statistical area of ​​land use in each year by the proportion of the total area of ​​the study area, iteratively calculating the proportion of the area transferred out and transferred in of each land use type, and obtaining the reconstructed data of the spatial distribution of land use containing conversion information includes the steps of: The statistical area of ​​land use in each year is sorted from small to large according to the proportion of the total area of ​​the study area, and the area ratio of each land use type transferred out and transferred in is iteratively calculated to obtain the statistical probability matrix of land use conversion from each year to the next year; Based on the statistical probability matrix, summing the proportion of each land use type that undergoes type conversion and is converted to other land use types from each year to the next year, and calculating the proportion of each land use type in the total sum; By multiplying the area of ​​each land use type by the proportion of the total sum, the occurrence type conversion of each land use type and the area converted to other land use types are calculated, thus realizing the historical reconstruction of the spatial distribution of land use changes including conversion information.

4. The method for historical reconstruction of spatial distribution of land use changes including conversion information according to claim 3, characterized in that: The process of iteratively calculating the proportion of the area transferred out and transferred in of each land use type by sorting the statistical area of ​​each year in the total area of ​​the study area from small to large is as follows: The smallest land use type is calculated according to the following formula The conversion area is calculated, and the ratio of the area transferred out from year y to year y+1 is: ； in, It means that from year y to year y+1 is The proportion of the area converted to type j land use in the total area of ​​the study area, where j is the land use type number; express The proportion of type land use in the total area of ​​the study area in year y; Represents from year y to year y+1 The proportion of the area with no type conversion in land use to the total area of ​​the study area; Represents the remote sensing probability matrix The average proportion of the area converted to type j land use to the total area of ​​the study area; n is the number of land use types; According to the following formula The conversion ratio of land use type from year y to year y+1 is calculated as: ； in, Indicates that from year y to year y+1, type i land use is converted to The proportion of the area of ​​the type of land use to the total area of ​​the study area, i is the land use type number; express The proportion of type land use in the total area of ​​the study area in year y + 1; Indicates that the remote sensing probability matrix is ​​converted from year y to year y+1. The average proportion of the area of ​​this type of land use to the total area of ​​the study area.

5. The method for historical reconstruction of spatial distribution of land use changes including conversion information according to claim 4, characterized in that: The process of iteratively calculating the proportion of the area transferred out and transferred in of each land use type by sorting the statistical area of ​​each year in the total area of ​​the study area from small to large is as follows: The conversion area of ​​the second smallest land use type m is calculated according to the following formula, and its conversion area ratio from year y to year y+1 is: ； in, It represents the proportion of the area converted from type m land use to type j land use from year y to year y+1 to the total area of ​​the study area, where j is the land use type number; It represents the proportion of land use type m in the total area of ​​the study area in year y; It represents the proportion of area where land use type m does not undergo type conversion from year y to year y+1; Indicates the land use conversion of type m from year y to year y+1 The proportion of the area of ​​the type to the total area of ​​the study area; It represents the average proportion of the area of ​​land use type m converted to land use type j from year y to year y+1 in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types; The proportion of land use area transferred from year y to year y+1 of type m is calculated according to the following formula: ； in, It represents the proportion of the area converted from type i land use to type m land use from year y to year y+1 to the total area of ​​the study area, where m is the land use type number; It represents the proportion of land use type m in the total area of ​​the study area in year y+1; Represents from year y to year y+1 The proportion of area converted to type m land use in the total area of ​​the study area; It represents the average proportion of the area of ​​land use type i converted to land use type m from year y to year y+1 in the remote sensing probability matrix to the total area of ​​the study area.

6. A device for historical reconstruction of the spatial distribution of land use changes including conversion information, characterized in that include: Basic data module, used to obtain remote sensing probability matrix of land use conversion; The first conversion module is used to calculate the ratio of the theoretical non-conversion area of ​​each land use type to the average area of ​​the corresponding land use type based on the remote sensing probability matrix according to the following formula: ； ； ； Among them, k, i and j represent the numbers of land use types. It refers to the proportion of the area where the type of land use k is converted from a certain initial year to the next year to the total area of ​​the study area. represents the proportion of land use area of ​​type k in the initial year to the total area of ​​the study area; and is the average proportion of the area of ​​each land use type conversion in the remote sensing probability matrix to the total area of ​​the study area; n is the number of land use types; The second conversion module is used to calculate the land use type conversion of the study area for any year, obtain the statistical area of ​​each type of land use in the study area in year y and year y+1, and calculate its proportion of the total area of ​​the study area, marked as and ; Where k represents the number of the land use type; Calculate the proportion of area where type k land use does not change from year y to year y+1 : ； It represents the ratio of the theoretical non-conversion area of ​​each land use type from year y to year y+1 to the average area of ​​the corresponding land use type; The data iteration module is used to sort the statistical area of ​​land use in each year according to the proportion of the total area of ​​the study area, and iteratively calculate the area ratio of each land use type transferred out and transferred in, so as to obtain the reconstructed data of the spatial distribution of land use containing conversion information.

7. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium stores computer-executable instructions, which, when executed by a processor, implement the method for historical reconstruction of spatial distribution of land use changes containing conversion information as described in any one of claims 1 to 5.

8. A data control device, characterized in that: include: one or more memories non-transitorily storing computer-executable instructions; One or more processors are configured to run computer executable instructions, wherein the computer executable instructions, when executed by the one or more processors, implement the method for historical reconstruction of the spatial distribution of land use changes containing conversion information as described in any one of claims 1 to 5.

Images