A land survey change analysis method and system that eliminates false changes

Abstract

The application discloses a land survey change analysis method and system for eliminating false changes and relates to the technical field of land space data processing. The method comprises the following steps: obtaining a superimposed patch set by spatially superimposing and analyzing two-adjust patch data and three-adjust patch data; screening, from the superimposed patch set, superimposed patches corresponding to information in which source layer information is linear features and land class change information is that land class has not changed, as a reference anchor patch; constructing a local search window based on the reference anchor patch and the width of the linear features; determining, according to a preset search condition and the local search window, forward change patches and reverse change patches corresponding to each land class change information from the superimposed patch set; and combining a preset offset rule to offset land class change flow, and determining a target land class change flow statistical result. The application eliminates false changes caused by differences in linear and planar expression through anchor point screening and area offset mechanism, and avoids the misdeletion of real changes caused by excessive filtering.

Description

Technical Field

[0001] This invention belongs to the field of land spatial data processing technology, and in particular relates to a method and system for analyzing land survey changes to eliminate spurious changes. Background Technology

[0002] The Second National Land Survey (hereinafter referred to as the Second Survey) and the Third National Land Survey (hereinafter referred to as the Third Survey) are important land surveys in my country. The two surveys differ significantly in their technical standards. Specifically, the Second Survey used linear features to represent linear land features with widths less than a certain threshold (such as rural roads, ditches, and forest belts); while the Third Survey uniformly used isometric land feature patches to represent all land types.

[0003] This difference in standards leads to numerous spurious changes when directly overlaying second and third national land survey data. For example, a rural road might be captured as a line (crossing farmland) in the second land survey, based on its centerline and average width from higher-resolution remote sensing imagery. However, in the third land survey, the boundary of the land parcel is captured based on the actual land area covered by higher-resolution remote sensing imagery, resulting in a planar land parcel. Because of these inconsistencies in data collection methods and precision, overlaying the data might statistically show a "linear road" from the second survey becoming "farmland" in the third, and vice versa, creating a bidirectional change flow. In reality, however, the road hasn't undergone any actual land use change; it's merely a difference in representation and collection precision. These land use change flows caused by differences in representation are spurious changes.

[0004] Currently, conventional change analysis methods (i.e., direct superposition statistical methods that do not consider pseudo-change issues) include all the aforementioned pseudo-changes in change transfer flows, resulting in the inclusion of pseudo-change components in the statistical results. This leads to significant errors in the statistical results and makes it impossible to truly reflect the actual land use changes. Summary of the Invention

[0005] To address the aforementioned problems in the existing technology, this invention provides a method and system for analyzing land survey changes by eliminating spurious changes.

[0006] The technical problem to be solved by this invention is achieved through the following technical solution: In a first aspect, the present invention provides a method for analyzing land survey changes to eliminate spurious changes, comprising: Obtain the patch data from the Second National Land Survey and the Third National Land Survey; the patch data from the Second National Land Survey includes source layer information and the width of the source being linear features; Spatial overlay analysis was performed on the second and third land survey parcel data to obtain an overlay parcel set. The overlay parcel set includes multiple overlay parcels and corresponding land use change information, which includes information on land use changes and information on land use unchanged land use. From the set of overlay patches, select the overlay patches whose source layer information is linear features and whose land use change information is that the land use has not changed, and use them as the benchmark anchor patches. Based on the width of the benchmark anchored patch and the corresponding linear feature, a local search window is constructed based on the benchmark anchored patch; According to the preset search conditions and the local search window based on the benchmark anchored patches, the positive change patches and negative change patches corresponding to the land use change information of each land use change are determined from the set of overlaid patches; the land use change information of the positive change patches indicates that the land use type changes from a certain linear feature to a certain non-linear feature; the land use change information of the negative change patches indicates that the land use type changes from a certain non-linear feature to a certain linear feature. Based on preset offsetting rules, positive and negative change patches corresponding to land type change information are used to offset land type change flows and determine the statistical results of target land type change flows.

[0007] Secondly, the present invention provides a land survey change analysis system for eliminating spurious changes, comprising: The data acquisition module is used to acquire the patch data from the Second National Land Survey and the Third National Land Survey; the patch data from the Second National Land Survey includes source layer information and the width of the source being linear features; The overlay analysis module is used to perform spatial overlay analysis on the second and third land survey parcel data to obtain an overlay parcel set. The overlay parcel set includes multiple overlay parcels and corresponding land use change information, which includes information on land use changes and information on land use unchanged land use. The anchoring patch determination module is used to select overlay patches from the set of overlay patches whose source layer information is linear features and whose land use change information is that the land use has not changed, as the benchmark anchoring patches; The search window determination module is used to construct a local search window based on the width of the benchmark anchored patch and the corresponding linear feature. The forward and reverse change flow determination module is used to determine the forward and reverse change patches corresponding to the land use change information from the overlay patch set according to preset search conditions and a local search window based on the benchmark anchored patch. The land use change information of the forward change patch indicates that the land use type changes from a linear feature to a non-linear feature; the land use change information of the reverse change patch indicates that the land use type changes from a non-linear feature to a linear feature. The change flow offset module is used to offset land use change flow based on preset offset rules, positive change patches and negative change patches corresponding to land use change information, and to determine the statistical results of target land use change flow.

[0008] Thirdly, the present invention provides an electronic device, including a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, used to store computer programs; When a processor executes a computer program stored in memory, it implements the steps of a land survey change analysis method for eliminating spurious changes, as described in the first aspect.

[0009] Fourthly, the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of a land survey change analysis method for eliminating spurious changes, as described in the first aspect.

[0010] This invention provides a method and system for analyzing land survey changes to eliminate spurious changes. Through anchor point screening and area offsetting mechanisms, it identifies and deducts spurious changes caused by differences in line and surface representation. Furthermore, since the identified spurious changes are essentially "suspected spurious changes," this invention uses offsetting deduction rather than direct zeroing to avoid the potential for over-filtering and the misdeletion of genuine changes. It employs a benchmark anchored patch as a spatial reference to construct a local search window, making spurious change identification more targeted. Using the width of the corresponding linear feature as a buffer radius reflects the application of geological principles, namely that the spatial positional deviation of the same physical feature in two surveys usually does not exceed its own width. This gives the invention adaptive capabilities, allowing for different search ranges for linear features of different widths. Because the search only needs to be conducted locally near the benchmark anchored patch, and the width of the linear feature is typically several meters, the search range is small, the computational load is controllable, and the computational efficiency is high.

[0011] The present invention will now be described in further detail with reference to the accompanying drawings. Attached Figure Description

[0012] Figure 1 This is a flowchart illustrating a land survey change analysis method for eliminating spurious changes provided in an embodiment of the present invention; Figure 2 This is a schematic diagram of the spatial overlay analysis process of a land survey change analysis method for eliminating spurious changes provided in an embodiment of the present invention; Figure 3 This is a structural block diagram of a land survey change analysis system for eliminating spurious changes provided in an embodiment of the present invention; Figure 4This is a schematic diagram of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0013] The present invention will be further described in detail below with reference to specific embodiments, but the implementation of the present invention is not limited thereto.

[0014] This invention provides a method for analyzing land survey changes to eliminate spurious changes. See also... Figure 1 The method includes the following steps: S10. Obtain the patch data from the Second National Land Survey and the Third National Land Survey.

[0015] The second-level map patch data includes source layer information and the width of the source being linear features.

[0016] For example, in the original land use data (including linear features and land use patches) obtained from the Second National Land Survey, linear features are those with a width less than a certain threshold (such as rural roads, ditches, and forest belts), which are represented by linear elements in the original land use data of the Second Survey. The land use patch data of the Second Survey in this embodiment refers to the data obtained after converting the linear features in the original land use data of the Second Survey into area patches through buffer analysis, and then trimming and merging them with the surrounding land use patches to form a complete full-pattern result. This conversion can be implemented using the technical solution in Chinese Patent Application No. 202010770183.0, entitled "A Method for Land Use Patch Transformation of Linear Features Based on Total Quantity Control in Land Surveys," which will not be elaborated upon in this embodiment.

[0017] Furthermore, this embodiment adds source layer information and the width of linear features as the source to the second land survey patch data. The source layer information includes linear features and land cover patches. Linear features represent patches obtained by patching linear features from the original land cover data of the second land survey, while land cover patches represent the corresponding land cover patches from the original land cover data of the second land survey. The assignment rule is: source_flag = 1 indicates that the source layer information is linear features; source_flag = 0 indicates that the source layer information is land cover patches. Additionally, for patches whose source layer information is linear features, the width of the corresponding linear feature is recorded in the second land survey patch data.

[0018] S20. Perform spatial overlay analysis on the second and third land survey patch data to obtain an overlay patch set.

[0019] The overlay set of map features includes multiple overlay map features and corresponding land use change information. The land use change information includes information on land use changes and information on land use unchanged land use.

[0020] For example, refer to Figure 2 Each overlay patch in the overlay patch set can include a set of corresponding second-level and third-level land survey patch data. The corresponding land use change information is used to characterize the land use change between the second-level and third-level land survey patch data in the same set. For example, for a set of second and third national land survey patch data where the land type is paddy field, the land type has not changed, and the corresponding land type change information can be "paddy field → paddy field"; for a set of second and third national land survey patch data where the land type is rural road, the land type has not changed, and the corresponding land type change information can be "rural road → rural road"; for a set of second national land survey patch data where the land type is rural road and third national land survey patch data where the land type is paddy field, the land type changes from rural road to paddy field, and the corresponding land type change information can be "rural road → paddy field"; for a set of second national land survey patch data where the land type is paddy field and third national land survey patch data where the land type is rural road, the land type changes from paddy field to rural road, and the corresponding land type change information can be "paddy field → rural road".

[0021] S30. From the set of overlay patches, select the overlay patches whose source layer information is linear features and whose land use change information is that the land use has not changed, and use them as the reference anchor patches.

[0022] For example, refer to Figure 2 If, in a certain overlay patch, the source layer information of the second land survey patch data is a linear feature (source_flag = 1), such as a rural road, and its corresponding land type change information is information where the land type has not changed, such as "rural road → rural road", then this overlay patch is used as a reference anchor patch. Such patches represent the same physical feature identified as the same land type in both surveys and can be used as spatial reference points. One or more reference anchor patches can be obtained from the set of overlay patches.

[0023] S40. Based on the width of the benchmark anchored patch and the corresponding linear feature, construct a local search window based on the benchmark anchored patch.

[0024] Optionally, step S40 may specifically include: The boundary of the benchmark anchored patch is extended outward along the width direction to the width of the corresponding linear feature, and the extended area is used as a local search window based on the benchmark anchored patch.

[0025] For example, refer to Figure 2 Using the geometric boundary of the baseline anchored patch as the source element, the width of the corresponding linear feature is extended outward on both sides of the baseline anchored patch along its width direction, so that an extension area can be obtained on each side of the baseline anchored patch, and the extension area is used as a local search window.

[0026] Understandably, due to geometric transformation errors during the second land survey's linear feature patch conversion process, and the accuracy deviations in the boundary interpretation of the third land survey's isometric land cover patches, different representations of the same physical feature may exhibit slight spatial shifts. Therefore, using the width of the linear feature corresponding to the benchmark anchor patch as the extension distance of the local search window can effectively cover this shift range and improve the accuracy of subsequent searches for positively and negatively changing patches.

[0027] S50. Based on the preset search conditions and the local search window of the benchmark anchored map pieces, determine the positive change map pieces and negative change map pieces corresponding to the land cover change information from the set of overlaid map pieces.

[0028] Among them, the land use change information of positive change patches indicates that the land use type changes from a linear feature to a non-linear feature; the land use change information of negative change patches indicates that the land use type changes from a non-linear feature to a linear feature.

[0029] Optionally, step S50 may specifically include: Based on the land cover change information, from the set of overlay patches, we select overlay patches whose land cover changes from the same linear feature to the same non-linear feature and whose overlap area with the local search window exceeds a preset threshold, and obtain the positive change patches of the corresponding land cover change information. From the set of overlay patches, filter the overlay patches whose land use type changes from the same non-linear feature to the same linear feature, and whose overlap area with the local search window exceeds a preset threshold, to obtain the reverse change patches of the corresponding land use type change information.

[0030] For example, the preset threshold can be flexibly set by those skilled in the art according to the actual situation; for instance, the preset threshold can be 80%. For example, if a local search window contains an overlay patch (land use change information: rural road → dry land) with a total area of ​​50 square meters, and the overlap area between this patch and the local search window is 45 square meters, the overlap ratio is 90% (>80%), then it meets the selection criteria for a positive change patch of land use change information. Another overlay patch (land use change information: rural road → dry land) with a total area of ​​50 square meters overlaps with the local search window with an area of ​​30 square meters, the overlap ratio is 60% (<80%), which does not meet the selection criteria for a positive change patch and therefore will not participate in the subsequent offsetting statistics. By setting an 80% overlap threshold, it ensures that only overlay patches that are spatially highly correlated participate in the offsetting, avoiding the erroneous inclusion of changes with large positional deviations in pseudo-change processing.

[0031] S60. Based on the preset offsetting rules, the positive and negative change patches corresponding to the land type change information, offset the land type change flow, and determine the statistical results of the target land type change flow.

[0032] Optionally, step S60 may specifically include: Based on positive and negative change patches, the total positive and negative change areas corresponding to the land cover change information are statistically obtained. The minimum of the total positive change area and the total negative change area corresponding to each land type change information is determined as the offset amount corresponding to each land type change information; The corresponding offset amount is deducted from the total positive and negative change areas corresponding to the land type change information to obtain the statistical results of the target land type change flow.

[0033] For example, in step S50, multiple positive change patches and multiple negative change patches corresponding to each land type change information may be selected, or only positive change patches or negative change patches may be selected. The areas of all positive change patches and all negative change patches corresponding to each land type change information are counted respectively to obtain the total positive change area A_forward and the total negative change area A_reverse corresponding to each land type change information.

[0034] If A_forward > 0 and A_reverse > 0, then area offsetting is performed. The offset amount C = min(A_forward, A_reverse) corresponding to each land cover change information is calculated. Then, in the initial land cover change flow statistics table obtained from the spatial overlay analysis of the second and third land cover changes, the corresponding offset amount C is subtracted from the forward change flow (i.e., the total forward change area A_forward) and reverse change flow (i.e., the total reverse change area A_reverse) corresponding to each land cover change information, respectively. This yields the forward change retention amount R_forward = A_forward – C, and the reverse change retention amount R_reverse = A_reverse – C. If A_forward > 0 and A_reverse = 0, i.e., offset amount C = 0, then all forward change flows are retained; if A_forward = 0 and A_reverse > 0, i.e., offset amount C = 0, then all reverse change flows are retained. Specific details are shown in Table 1 below. Table 1. Area cancellation results for A_forward and A_reverse under different conditions.

[0035] For example, if the total positive change area (rural road → paddy field) for a certain land use change is A_forward = 100 square meters, and the total negative change area (paddy field → rural road) is A_reverse = 60 square meters, then the corresponding offset amount C = min(100, 60) = 60 square meters; the positive change retention amount R_forward = 100 - 60 = 40 square meters, and the negative change retention amount R_reverse = 60 - 60 = 0 square meters. By statistically analyzing the positive and negative change retention amounts corresponding to all land use change information, the statistical results of the target land use change flow are obtained. A comparison with the ordinary change analysis method is shown in Table 2 below: Table 2 Comparison of results between the present invention and conventional change analysis methods

[0036] In this embodiment, the identified spurious changes are "suspected spurious changes," which cannot be 100% confirmed as spurious changes. Therefore, this embodiment uses a deduction method instead of direct zeroing for processing. When the forward and reverse areas are unequal, the difference is retained as a suspected real change, thereby avoiding the accidental deletion of possible real changes due to over-filtering and improving the accuracy of the target land category change flow statistics.

[0037] The positive and negative change retention amounts after offsetting are summarized, and the statistical results of target land category change flow are output.

[0038] Optionally, the statistical results of the target land use change flow are presented in tabular form. Rows in the table correspond to land use types in the second national land use survey patch data, columns correspond to land use types in the third national land use survey patch data, and values ​​correspond to the area values ​​of changes in land use types between the second and third national land use survey patch data. Compared to existing technologies, the changed area in this result largely eliminates the influence of spurious changes and more closely approximates the actual land use change situation.

[0039] The present invention provides a land survey change analysis method that eliminates spurious changes. Compared with ordinary change analysis methods (direct superposition statistical methods that do not consider spurious changes), the present invention has the following advantages: 1. This invention can eliminate the impact of spurious changes at the statistical level while avoiding over-filtering: Unlike ordinary change analysis methods that directly include all superimposed results in the statistics, this invention identifies and deducts spurious changes caused by differences in line and surface representation through anchor point screening and area cancellation mechanisms. Furthermore, since the identified spurious changes are essentially "suspected spurious changes," this invention uses cancellation deduction rather than direct zeroing, retaining the difference when the positive and negative areas are unequal, thus avoiding the erroneous deletion of genuine changes that might result from over-filtering.

[0040] 2. Using benchmark anchored patches as the baseline provides high spatial reliability: Unlike ordinary change analysis methods that do not distinguish the source of patches, this invention selects benchmark anchored patches with "linear source + consistent land type" as the spatial baseline. These patches represent the same physical feature that was correctly identified in both surveys. Based on this, a local search window is constructed, making the identification of false changes more targeted.

[0041] 3. Using the original line feature width as the buffer radius has clear physical significance: Unlike ordinary change analysis methods that do not involve local searches, this invention uses the width of the corresponding linear feature as the buffer radius, reflecting the application of geological laws—the spatial position deviation of the same physical feature in two surveys usually does not exceed its own width. This scheme has adaptive capabilities, employing different search ranges for linear features of different widths.

[0042] 4. Improved matching accuracy through overlap condition: When searching for positive and negative change patches, this invention requires that the overlap area between the patch and the local search window is greater than 80%, ensuring that only change patches that are highly correlated in spatial location participate in cancellation, further improving the accuracy of false change identification.

[0043] 5. High computational efficiency: Since only a local search is needed near the benchmark anchored patch, and the width of linear features is usually a few meters, the search range is small and the computational load is controllable.

[0044] Secondly, referring to Figure 3 This invention also provides a land survey change analysis system for eliminating spurious changes, comprising: The data acquisition module 301 is used to acquire the second national land survey patch data and the third national land survey patch data; the second national land survey patch data includes source layer information and the width of the source being linear features; The overlay analysis module 302 is used to perform spatial overlay analysis on the second and third land survey parcel data to obtain an overlay parcel set. The overlay parcel set includes multiple overlay parcels and corresponding land use change information, which includes information on land use changes and information on land use unchanged. The anchoring patch determination module 303 is used to select from the set of overlay patches information that the source layer information is linear features and the land type change information is that the land type has not changed, as the reference anchoring patches; The search window determination module 304 is used to construct a local search window based on the reference anchor patch and the width of the corresponding linear feature. The forward and reverse change flow determination module 305 is used to determine the forward change patches and reverse change patches corresponding to the land use change information from the set of overlaid patches according to preset search conditions and a local search window based on the benchmark anchored patches; the land use change information of the forward change patches indicates that the land use type changes from a linear feature to a non-linear feature; the land use change information of the reverse change patches indicates that the land use type changes from a non-linear feature to a linear feature. The change flow offset module 306 is used to offset land use change flow based on preset offset rules, positive change patches and negative change patches corresponding to land use change information, and to determine the statistical results of target land use change flow.

[0045] For details on the system, please refer to the steps of the land survey change analysis method for eliminating spurious changes provided in the first aspect; they will not be repeated here.

[0046] Thirdly, embodiments of the present invention also provide an electronic device, such as... Figure 4 As shown, it includes a processor 401, a communication interface 402, a memory 403, and a communication bus 404, wherein the processor 401, the communication interface 402, and the memory 403 communicate with each other through the communication bus 404. Memory 403 is used to store computer programs; The processor 401, when executing the program stored in the memory 403, implements a land survey change analysis method for eliminating spurious changes, as described in the first aspect.

[0047] The communication bus mentioned in the above electronic devices can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. This communication bus can be divided into address bus, data bus, control bus, etc. For ease of illustration, only one thick line is used to represent it in the diagram, but this does not mean that there is only one bus or one type of bus.

[0048] The communication interface is used for communication between the aforementioned electronic devices and other devices.

[0049] The memory may include random access memory (RAM) or non-volatile memory (NVM), such as at least one disk storage device. Optionally, the memory may also be at least one storage device located remotely from the aforementioned processor.

[0050] The processors mentioned above can be general-purpose processors, including central processing units (CPUs), network processors (NPs), etc.; they can also be digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0051] The method provided in this invention can be applied to electronic devices. Specifically, the electronic device can be a desktop computer, a portable computer, a smart mobile terminal, a server, etc. No limitation is made herein; any electronic device that can implement this invention falls within the protection scope of this invention.

[0052] Fourthly, the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of a land survey change analysis method for eliminating spurious changes, as described in the first aspect.

[0053] For system / electronic device / storage medium embodiments, since they are basically similar to method embodiments, the description is relatively simple, and relevant parts can be referred to in the description of the method embodiments.

[0054] It should be noted that the system, electronic device, and storage medium in the embodiments of the present invention are respectively the apparatus, electronic device, and storage medium for applying the above-mentioned method for analyzing land survey changes by eliminating false changes. Therefore, all embodiments of the above-mentioned method for analyzing land survey changes by eliminating false changes are applicable to the system, electronic device, and storage medium, and can achieve the same or similar beneficial effects.

[0055] It should be noted that the terms "first," "second," etc., are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present invention. Rather, they are merely examples of apparatuses and methods consistent with some aspects of the invention.

[0056] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Furthermore, those skilled in the art can combine and integrate the different embodiments or examples described in this specification.

[0057] Although the invention has been described herein in conjunction with various embodiments, those skilled in the art will understand and implement other variations of the disclosed embodiments by reviewing the accompanying drawings and the disclosure in carrying out the claimed invention. In the description of the invention, the word "comprising" does not exclude other components or steps, "a" or "an" does not exclude a plurality, and "a plurality" means two or more, unless otherwise explicitly specified. Furthermore, while different embodiments may describe certain measures, this does not mean that these measures cannot be combined to produce good results.

[0058] The above description, in conjunction with specific preferred embodiments, provides a further detailed explanation of the present invention. It should not be construed that the specific implementation of the present invention is limited to these descriptions. For those skilled in the art, various simple deductions or substitutions can be made without departing from the concept of the present invention, and all such modifications and substitutions should be considered within the scope of protection of the present invention.

Claims

1. A method for analyzing land survey changes to eliminate spurious changes, characterized in that, include: Acquire the patch data from the Second National Land Survey and the Third National Land Survey; the patch data from the Second National Land Survey includes source layer information and the width of the source being linear features; Spatial overlay analysis is performed on the second-level and third-level land parcel data to obtain an overlay parcel set; the overlay parcel set includes multiple overlay parcels and corresponding land use change information, the land use change information includes information on land use changes and information on land use unchanged; From the set of overlaid patches, overlaid patches corresponding to information where the source layer information is linear features and the land use change information is that the land use has not changed are selected as reference anchor patches; Based on the width of the reference anchored patch and the corresponding linear feature, a local search window is constructed based on the reference anchored patch; Based on the land cover change information, from the set of overlay patches, overlay patches whose land cover changes from the same linear feature to the same non-linear feature and whose overlap area with the local search window of the reference anchor patch exceeds a preset threshold are selected to obtain the positive change patches of the corresponding land cover change information. From the set of overlaid patches, patches whose land use type changes from the same non-linear feature to the same linear feature, and whose overlap area with the local search window exceeds the preset threshold, are selected to obtain reverse change patches with corresponding land use change information; the land use change information of the positive change patches indicates that the land use type changes from a certain linear feature to a certain non-linear feature; the land use change information of the reverse change patches indicates that the land use type changes from a certain non-linear feature to a certain linear feature. Based on preset offsetting rules and the positive and negative change patches corresponding to the land type change information, land type change flow offsetting is performed to determine the target land type change flow statistical results.

2. The method for analyzing land survey changes to eliminate spurious changes according to claim 1, characterized in that, The step of constructing a local search window based on the reference anchored patch and the width of the corresponding linear feature includes: Using the boundary of the reference anchor patch as a reference, the width of the corresponding linear features is extended outward from both sides of the reference anchor patch along the width direction, and the extended area is used as a local search window based on the reference anchor patch.

3. The method for analyzing land survey changes to eliminate spurious changes according to claim 1, characterized in that, The process of offsetting land use change flows based on preset offsetting rules, the positive change patches, and the negative change patches to determine the statistical results of target land use change flows includes: Based on the positive change patches and the negative change patches, the total positive change area and the total negative change area corresponding to the land cover change information are statistically obtained; The minimum value between the total positive change area and the total negative change area corresponding to each land type change information is determined as the offset amount corresponding to each land type change information; The corresponding offset amount is deducted from the total positive change area and the total negative change area corresponding to the land type change information of each region to obtain the statistical results of the target land type change flow.

4. The method for analyzing land survey changes to eliminate spurious changes according to claim 1, characterized in that, The statistical results of the target land category change flow are presented in tabular form. The rows in the table correspond to the land categories in the second land survey patch data, the columns correspond to the land categories in the third land survey patch data, and the values ​​correspond to the area values ​​of the changes in land categories between the second and third land survey patch data.

5. A land survey change analysis system for eliminating spurious changes, characterized in that, include: The data acquisition module is used to acquire the second national land survey patch data and the third national land survey patch data; the second national land survey patch data includes source layer information and the width of the source being linear features; The overlay analysis module is used to perform spatial overlay analysis on the second-level land parcel data and the third-level land parcel data to obtain an overlay parcel set; the overlay parcel set includes multiple overlay parcels and corresponding land use change information, the land use change information includes information on land use changes and information on land use unchanged; The anchoring patch determination module is used to select from the set of overlay patches information that the source layer information is a linear feature and the land use change information is that the land use has not changed, as the reference anchoring patches; The search window determination module is used to construct a local search window based on the reference anchor patch and the width of the corresponding linear feature. The forward and reverse change flow determination module is used to filter overlay patches from the overlay patch set based on the land type change information, where the land type changes from the same linear feature to the same non-linear feature, and the overlap area ratio with the local search window based on the reference anchor patch exceeds a preset threshold, to obtain the forward change patch corresponding to the land type change information. From the set of overlaid patches, patches whose land use type changes from the same non-linear feature to the same linear feature, and whose overlap area with the local search window exceeds the preset threshold, are selected to obtain reverse change patches with corresponding land use change information; the land use change information of the positive change patches indicates that the land use type changes from a certain linear feature to a certain non-linear feature; the land use change information of the reverse change patches indicates that the land use type changes from a certain non-linear feature to a certain linear feature. The change flow offset module is used to offset land use change flow based on preset offset rules and the positive and negative change patches corresponding to land use change information, and to determine the target land use change flow statistical results.

Images