Incrementally Detecting Change Types of Spatial Area Object: A Hierarchical Matching Method Considering Change Process

: Detecting and extracting the change types of spatial area objects can track area objects’ spatiotemporal change pattern and provide the change backtracking mechanism for incrementally updating spatial datasets. To respond to the problems of high complexity of detection methods, high redundancy rate of detection factors, and the low automation degree during incrementally update process, we take into account the change process of area objects in an integrated way and propose a hierarchical matching method to detect the nine types of changes of area objects, while minimizing the complexity of the algorithm and the redundancy rate of detection factors. We illustrate in details the identiﬁcation, extraction, and database entry of change types, and how we achieve a close connection and organic coupling of incremental information extraction and object type-of-change detection so as to characterize the whole change process. The experimental results show that this method can successfully detect incremental information about area objects in practical applications, with the overall accuracy reaching above 90%, which is much higher than the existing weighted matching method, making it quite feasible and applicable. It helps establish the corresponding relation between new-version and old-version objects, and facilitate the linked update processing and quality control of spatial data.


Introduction
Spatial databases are mainly updated incrementally, and one difficulty with the incremental updating of spatial databases is the extraction of the said incremental (i.e., changed) information [1][2][3].
The detection and extraction of information about object changes is the basis for the incremental updating of spatial databases [4,5].Recording complete process about object changes should not only include the changed object itself, but should also contain information about the types of changes that trigger the physical change of the object [6][7][8].Recording the changed object is the key to collecting and storing incremental update information, while further detection of the types of changes influences the linked update processing, quality control, and release of the incremental change information [9].Information about object change can be applied in fields such as spatial object update, lifecycle tracking, historical data backtracking, statistics, and analysis of change information, as well as prediction of the spatiotemporal transmission patterns and change trend of spatial objects [10][11][12].In updating practices, if incremental updating only record those old object in history database, with paying little attention to change process of geographical objects, then multiple updates will hide inevitably the corresponding connection between the old and new data.
On the one hand, the nature of incremental extraction is the matching of spatial objects, and feature matching is often used to detect the changed objects.This is because since single-geometric-feature (e.g., centroids) matching usually cannot accurately measure the similarity among different objects or determining spatial relationships of polygon objects.Although there indeed exist many other alternative approaches (e.g., MBR-Minimum Bounding Rectangle; Hausdorff, Voronoi, etc.) that solve these problems, most researchers prefer integrating a number of operators reflecting the different geometric characteristics of the objects from the perspective of reducing the complexity of the algorithm [10,[13][14][15][16][17].One of the typical methods conducts a weighted processing of several operators reflecting different geometric features of the objects, obtaining a comprehensive similarity index to judge whether there is a geometric matching relationship between those objects [11,[18][19][20].This method avoids the selection of single operator thresholds by increasing the complexity and cost of the algorithm, and can better describe the overall geometric features of the objects, greatly improving the matching integrity and accuracy.However, both the weighted method and the single operator matching method ignore an important fact: geometric matching requires that similarity in position, shape, and size should be satisfied simultaneously, and similarity in only one or two geometric features is not enough [12].They are liable to produce matching redundancy, and to affect the matching accuracy due to different weight assignment methods.
On the other hand, as change types are concerned, changes are usually divided into nine types in existing studies: vanish; appearance; reappearance; attribute change; expansion; shrink; translation; rotation, and; deformation [21][22][23][24].For the detection and extraction of incremental information, scholars mainly conduct studies from two perspectives: geographic entity change (change event-driven) and feature change (change information-driven).Of these, most researchers update a spatial database with information about the types of changes in spatial objects based on event-driven incremental updates [25][26][27], which alleviates the difficulty of detecting the types of changes during the incremental information extraction process.However, this method requires information about the types of changes ahead of time and conduct incremental updates with that information (geographic change event), which means that the information needs to be stored in the database before the update.In essence, this method requires manually adding the type-of-change information to the database; however, in most cases, when updating a large amount of spatial data, it is impossible to have information ahead of time about the types of changes experienced by these objects.In addition, the topological relations among objects, or non-topology attributes such as size, shape, and position, are also used to determine the types of changes [19,[27][28][29].Some scholars apply the snapshot difference triple-descriptive model [30] or the quartet model [31] to detect the types of changes.Compared with methods based on the change event, these methods (which determine the types of changes using object information) have made remarkable progress towards automation, but they face problems like high complexity in their detection method, high redundancy rates of detection factors, and rough, inaccurate inference of type-of-change information.Moreover, most of the existing literature either studies the extraction of changed objects or the detection of the types of changes, and there is little research that combines both the extraction of incremental information about area objects and the determination of the types of changes in order to completely detect the change information and to record the change process.The link between updating data and history data remains a difficult issue to date [8,11].
To solve the above problems, this study considers the organic unity of the extraction of change information about the objects and the detection of the types of changes, and proposes a method to completely detect the incremental information about area objects.Taking into account the process of object change, we propose an algorithm which implements the identification, extraction, and database entry of the types of changes of area objects during the incremental information extraction process.
While ensuring the accuracy of the algorithm, both its complexity and the redundancy of the detection factors involved are reduced to the greatest possible extent.

Incremental Detection of Change Information of Area Objects
We propose a new, universal method that is based on object matching features and detects both the changing objects and their types of changes.This method first establishes a hierarchical matching model to extract incremental information datasets of area objects in both the old and new versions.Subsequently, based on the attributes, geometric position, size, shape and directional information of the area objects, this method defines appropriate matching operators to define the rules to automatically detect the types of changes.Then, the method detect and extract the nine types of changes of individual area objects.It uses concise detection conditions and operational operators to ensure the accuracy of detection, reduce algorithmic complexity, and achieve effective tracking and management of various changing tracks within the lifecycle of area objects.The incremental information extraction process and its connection with the extraction of incremental objects are shown in Figure 1.

Incremental Detection of Change Information of Area Objects
We propose a new, universal method that is based on object matching features and detects both the changing objects and their types of changes.This method first establishes a hierarchical matching model to extract incremental information datasets of area objects in both the old and new versions.Subsequently, based on the attributes, geometric position, size, shape and directional information of the area objects, this method defines appropriate matching operators to define the rules to automatically detect the types of changes.Then, the method detect and extract the nine types of changes of individual area objects.It uses concise detection conditions and operational operators to ensure the accuracy of detection, reduce algorithmic complexity, and achieve effective tracking and management of various changing tracks within the lifecycle of area objects.The incremental information extraction process and its connection with the extraction of incremental objects are shown in Figure 1.

Basic Idea
The identification and extraction of change information depends largely on matching objects in the new-version database to the corresponding objects in the old-version database, and the matching accuracy will directly affect the accuracy of the data update.The matching of spatial objects includes both geometric matching and attribute matching.The geometrical characteristics of area objects include their position, size, shape, and direction [21], whose influence on feature matching should be fully considered.We use M to represent the matching model, with M = (Attribute, Position, Size, Shape, Direction).
To detect the change information, we propose to use matching operators reflecting information such as attribute, position, size, shape, and direction in the matching model M to identify the changes of the area object and discriminate the types of changes.When the attribute data are not included in the dataset, other geometric information apart from attribute can also be used to conduct geometric matching and discriminate the types of changes [32].Therefore, this study only focuses on geometric matching.The matching operators are screened and combined according to their complexity from the lowest to the highest, reducing the calculation workload while ensuring accuracy.For any two matching objects in the old and new databases under the same measuring scale, if any of the operators matches successfully, the condition judgment value will be denoted as T; otherwise, it will be denoted as F. Therefore, when M = (T, T, T, T, T), the match succeeds.After traversing all the objects, those that cannot be successfully matched are the changing ones.
In addition, the types of changes in the objects are detected.Vanish, appearance, and reappearance are a combination of delete and add, while the other six types of changes (attribute change, expansion, shrink, translation, rotation, and deformation) may be inferred using the combination of operators in Model M shown in Figure 2.For vanish and appearance, we only need to find out whether any one area object in a database intersects with another area object in the other database within the same geographical range.If there is object A in the new database while there is

Basic Idea
The identification and extraction of change information depends largely on matching objects in the new-version database to the corresponding objects in the old-version database, and the matching accuracy will directly affect the accuracy of the data update.The matching of spatial objects includes both geometric matching and attribute matching.The geometrical characteristics of area objects include their position, size, shape, and direction [21], whose influence on feature matching should be fully considered.We use M to represent the matching model, with M = (Attribute, Position, Size, Shape, Direction).
To detect the change information, we propose to use matching operators reflecting information such as attribute, position, size, shape, and direction in the matching model M to identify the changes of the area object and discriminate the types of changes.When the attribute data are not included in the dataset, other geometric information apart from attribute can also be used to conduct geometric matching and discriminate the types of changes [32].Therefore, this study only focuses on geometric matching.The matching operators are screened and combined according to their complexity from the lowest to the highest, reducing the calculation workload while ensuring accuracy.For any two matching objects in the old and new databases under the same measuring scale, if any of the operators matches successfully, the condition judgment value will be denoted as T; otherwise, it will be denoted as F. Therefore, when M = (T, T, T, T, T), the match succeeds.After traversing all the objects, those that cannot be successfully matched are the changing ones.
In addition, the types of changes in the objects are detected.Vanish, appearance, and reappearance are a combination of delete and add, while the other six types of changes (attribute change, expansion, shrink, translation, rotation, and deformation) may be inferred using the combination of operators in Model M shown in Figure 2.For vanish and appearance, we only need to find out whether any one area object in a database intersects with another area object in the other database within the same geographical range.If there is object A in the new database while there is no corresponding object of A in the old database, then such change can be classified as appearance, otherwise it belongs to vanish.For reappearance, it is necessary to identify objects in the databases over three periods, in which it is the succession of vanish and appearance.The specific solution and reasoning are shown in the following sections.no corresponding object of A in the old database, then such change can be classified as appearance, otherwise it belongs to vanish.For reappearance, it is necessary to identify objects in the databases over three periods, in which it is the succession of vanish and appearance.The specific solution and reasoning are shown in the following sections.

Selection of Hierarchical Matching Operators
Once the changing objects are discriminated, providing attribute data for each, only the attribute fields of two area objects need to be directly matched.However, when selecting geometric discriminant operators, their integrity must be guaranteed.Moreover, the objects' shape, size, and position should be considered as well to ensure the accuracy of the matching.In addition, the complexity of the algorithm should be taken into account.In the end, the detection model formed by all the operators should not be affected by the data types and the measurement scale, and it should be able to accurately describe both the overall and individual geometric features of the objects.
Among the discriminant operators, the simple and accurate centroid distance operator measures the position similarity of area objects and is used by many researchers [8,33].As an alternative, some

Selection of Hierarchical Matching Operators
Once the changing objects are discriminated, providing attribute data for each, only the attribute fields of two area objects need to be directly matched.However, when selecting geometric discriminant operators, their integrity must be guaranteed.Moreover, the objects' shape, size, and position should be considered as well to ensure the accuracy of the matching.In addition, the complexity of the algorithm should be taken into account.In the end, the detection model formed by all the operators should not be affected by the data types and the measurement scale, and it should be able to accurately describe both the overall and individual geometric features of the objects.
Among the discriminant operators, the simple and accurate centroid distance operator measures the position similarity of area objects and is used by many researchers [8,33].As an alternative, some scholars use the Hausdorff distance model, the generalized Hausdorff distance model, or the median Hausdorff distance model to measure the position similarity [25,34].This method not only measures the position difference of spatial objects, but also reflects differences between objects in terms of shape and overall distribution.However, only as far as the measurement of position similarity is concerned, the Hausdorff methods are not as accurate as the centroid matching method, lying in that they are highly complicated with heavy computational burden to a great extent.
As far as shape similarity is concerned, there are also various operators to measure shape similarity [20,28,31].Among them, the steering angle function matching algorithm can accurately describe the local, detailed features of complex graphics on their shape and direction, and it is easier to operate and more accurate than other matching algorithms.In addition, those other operators, for example, the fractal dimension and compactness index, the function describing the distance from edge feature point to the centroid of the area object, the Fourier transform method, the symmetric difference index, and the operator to measure the shape of tangent space, etc. can also measure the shape similarity of area objects [21,25,27].Although these operators can reflect the shape similarity between two area objects to a certain extent, they either fail to describe the details of the shape or involve complicated calculation.
Additionally, most researchers use the degree of area overlapping as a criterion to measure the size similarity between area objects [20].In addition, Hausdorff distance and symmetric difference could also reflect the size similarity of the objects to a certain extent [11].However, measuring the size similarity using the degree of area overlapping is not only simple to operate, but is also relatively high accuracy in nature.
In summary, after a systematic analysis and modification of the operators that come from the existing studies, this study uses the centroid discriminant operator to measure the position features of area objects, the area overlapping operator to measure their size features, and the steering angle cumulative function matching operator to measure their shape and direction characteristics.The specific calculation method is as follows: (1) Centroid discriminant operator.This is an index that reflects the distance (position) between objects.The smaller the centroid distance, the higher the position similarity between two area objects.Here, F 1 and F 2 denote the two area objects to be matched, D(F 1 , F 2 ) denotes the centroid distance between the two area objects, and r 1 and r 2 represent the half of the diagonal of the minimum enclosing rectangle of F 1 and F 2 , respectively, used so as not to be affected by the actual size of the objects.The formula is as follows: (2) Area discriminant operator.Assume that S(F 1 ) and S(F 2 ) represent the area of F 1 and that of F 2 , respectively.Formula (2), which discriminates the size similarity between area objects with area difference, is shown as follows: (3) Steering angle cumulative function discriminant operator.This operator measures the shape similarity between area objects and can measure the direction similarity between them to a certain extent.The principle is as follows: First, a shape feature point in a specific direction in the vertices of the two objects is chosen as a reference point.Giving consideration to information about the direction of the area objects, the reference point is selected per the following standard: the sine value of the angle between the line formed by connecting the reference point to the object centroid and the abscissa axis of the coordinate system should be the maximum, and the X axis coordinate of the reference point should be larger than that of the object centroid.The angle between the counterclockwise segmental arc of the reference point and the X axis is then recorded, along with the normalized length (the ratio of arc length to perimeter) of each counterclockwise arc.Then, the normalized length is taken as the X axis, and the accumulated value of each point along the peripheral steering angle is taken as the Y axis.Assuming that F 1 (x) and F 2 (x) are the steering angle cumulative functions of F 1 and F 2 , then the steering angle cumulative function discriminant operator formula can be expressed as follows: Although the operator can simultaneously measure the shape and direction information of area objects, certain shape or direction differences may cause the discriminant operator to be inadequate, in which case the combination of the above three operators cannot discriminate the type of rotation change.Therefore, in the case where the centroid discriminant operator and the area difference discriminant operator are adequate, while the steering angle function discriminant operator is inadequate, the auxiliary discriminant operator is introduced in this paper to discriminate the type of rotation change.If the computed result of the auxiliary operator is larger than the set threshold value, it means that the two objects have the same shape but different directions, and the type of changes can be identified as rotation.The discriminant operator measures the shape difference with the similarity of form factor [8]. Assume that S(F 1 ), S(F 2 ), L(F 1 ) and L(F 2 ) represent the areas and perimeters of the two area objects, F 1 and F 2 , respectively, then the formula is as follows:

Extraction of Incremental Information Based on Hierarchical Matching
This study proposes a geometric matching method for area object based on hierarchical matching, with the hierarchy determined by the accuracy of object matching, the increase of operator complexity, and the decrease of the number of filtered objects.Using the matching operators described above, the algorithm filters the three operators with increasing computational complexity layer by layer to reduce the calculational burden and ensure the accuracy of the matching.
In this study, the similarity threshold of the four geometric matching operators (namely O centroid , O area O angle and O Form ) are set as ϕ 1 , ϕ 2 , ϕ 3 and ϕ 4 respectively.Let IsAttributeMatch denote whether the attributes of the two area objects match or not.If they match, it will be denoted by 1 (True), otherwise by 0 (False).For any A and B, which are area objects from the new-version and old-version databases, respectively, the detailed matching method to extract the change information based on hierarchical feature matching is as follows, also as shown in Figure 2.
The first layer: Matching of centroid distance.This layer determines which pairs of candidate objects meet the 1:1 criterion of centroid distance matching, with the remaining mismatched ones regarded as changed objects.Centroid distance matching is conducted among all objects of the same semantic theme in the new and old databases in pairs.If the centroid distance between two objects is smaller than the given matching threshold ϕ 1 , the ID number of the pair of objects will be stored in the ID set IDCollection 1 for the matching calculation at the next layer.
The second layer: Matching of area overlapping degree.The algorithm's first layer returns a one-to-one relationship of the objects, indicating identical positions.Next, a one-to-one relationship at the second layer indicates that the two objects are not only matched in terms of centroid distance but also area similarity, with mismatched objects regarded as the changed ones.In IDCollection 1 , each pair of objects satisfying centroid distance matching conditions are traversed, calculating the area overlapping degree of the pair.If the overlapping degree is greater than the given area matching threshold ϕ 2 , then the two objects meet the criterion of the overlapping degree and will be re-saved to the new ID set (IDCollection 2 ).Pairs of objects ultimately retained are those matched to each other in both centroid distance and area overlapping degree.
The third layer: Matching of steering angle function.The algorithm proceeds only on those pairs matching in the first two layers.If the two objects are matched at the third layer as well, they are successfully matched geometrically; otherwise, they are changed objects.The matching process is similar to that of the second layer.In IDCollection 2 , each pair of objects is traversed, calculating their steering function matching.If it is greater than the matching threshold (ϕ 3 ), then the two objects are successfully matched in terms of shape.
After the above hierarchical matching, the finally retained pairs of objects are considered successfully matched, which can be formalized as follows.
MatchingRule: if ((IsAttributeMatch = True) and (PositionResult ≤ ϕ 1 ) and (AreaResult ≥ ϕ 2 ) and (ShapeDirectionResult ≥ ϕ 3 )) then Matching (A, B).From the above, the matching result set IDCollection (M) is obtained.Unsuccessfully matched objects in the old and new databases are considered to be the changed objects and are recorded as ChangeCollection (N) and ChangeCollection (O), respectively, to be used in the subsequent detection of the types of changes.
Let the numbers of objects in the new and old databases be m and n, respectively.We compare the hierarchical matching algorithm proposed in this study and the existing weighted matching algorithm from the perspective of algorithm complexity.Between the two methods, the algorithm complexity of the three operators, namely distance, area, and shape, is the same.If the number of pairs of objects meeting the distance matching is k 1 , and that of pairs of objects meeting the distance and area matching at the same time is k 2 , then Min(m, n) > k 1 > k 2 .For each operator, O centroid , O area , and O angle , the times of calculations conducted by the hierarchical matching algorithm is m * n, k 1 , k 2 , respectively, much less than the number of calculations conducted by the weighted algorithm, which is m*n, as shown in Table 1.Therefore, when the time complexity and spatial complexity of each operator are the same, the hierarchical matching algorithm is much less complicated than the weighted matching algorithm.

Detection of Change Types Based on Matching Operators
Based on the above extracted change datasets, that is, ChangeCollection (N) and ChangeCollection (O), assume that ChangeType denotes the detected results about the types of changes.Meanwhile, the corresponding dataset (i.e., GeoDatabase (P)), earlier than GeoDatabase (O), is introduced to help detect the change type of Reappearance.For F n , F o , and F p , which are any objects derived from ChangeCollection (N), ChangeCollection (O), and GeoDatabase (P), respectively, the rules to infer the nine basic types of changes of single area objects are defined as follows, also shown in Figure 3.
ChangeCollection (O), assume that ChangeType denotes the detected results about the types of changes.Meanwhile, the corresponding dataset (i.e., GeoDatabase (P)), earlier than GeoDatabase (O), is introduced to help detect the change type of Reappearance.For Fn, Fo, and Fp, which are any objects derived from ChangeCollection (N), ChangeCollection (O), and GeoDatabase (P), respectively, the rules to infer the nine basic types of changes of single area objects are defined as follows, also shown in Figure 3.

Experimental Test
To verify the feasibility of the method proposed in this study, C# programming was combined with ArcGIS Engine to conduct an experiment.Large-scale thematic maps of buildings in the same area of Fangshan District, Beijing City during different periods of time were selected to carry out the discriminating experiment.In consideration that few cases in real world contain the two change types of rotation and translation, while the type of reappearance is based on the change data in the three periods.Therefore, to fully prove the feasibility of this method, simulated buildings with the translation and rotation types of changes were added to the extracted set of changed objects in the experiment.
In the algorithm implementation, the threshold settings of the discriminant operators will directly affect the discriminating effect of the method.Since the OTSU (maximization of interclass variance) algorithm exists as one of popular methods on automatic threshold selection, it has a very wide range of applications in the image segmentation for its simplicity and intuitive.In case that the object and the background in image obey normal distributions, the selection of the global threshold could be realized by statistics of the histogram of the whole image [3,35].In this study, we found sampling values of each operators fit in with the principle of the OTSU method.Therefore, to get the optimal threshold, the adaptive threshold calculation method of geometric matching based on OSTU, was used in this study.Using this method, the thresholds of the four operators (O centroid , O area , O angle , and O Form ) in the algorithm to determine the types of changes were 0.85, 0.85, 0.71, and 0.95, respectively.A hierarchical matching experiment was carried out using the above thresholds, and a weighted matching was performed simultaneously among objects meeting the centroid matching conditions on the basis of centroid matching.
Additionally, in order to better measure and compare the experimental results, integrity (R i ) and accuracy (R a ) were chosen as the indices to evaluate the performance of the algorithm.The number of objects with correct results detected by this algorithm and the number of objects in the old and new databases with actual corresponding relationships with each other (i.e., matching relationships or the actual types of changes) are represented by t and T respectively, with R denoting the number of objects with corresponding relationships determined by the algorithm.The calculation method is as follows:

Incremental Information Extraction
From comparative analysis of the experimental results that are shown in Figure 4 and Table 2, the area object matching method based on hierarchical matching can not only achieve high accuracy and stability, but is also superior to the weighted algorithm when it comes to integrity and accuracy of matching.From the theoretical point of view, there are cases where some indices of the weighted matching method are too large or too small, but the final weighted results satisfy the set threshold, which was also verified by the experiment.
accuracy (Ra) were chosen as the indices to evaluate the performance of the algorithm.The number of objects with correct results detected by this algorithm and the number of objects in the old and new databases with actual corresponding relationships with each other (i.e., matching relationships or the actual types of changes) are represented by t and T respectively, with R denoting the number of objects with corresponding relationships determined by the algorithm.The calculation method is as follows:

Incremental Information Extraction
From comparative analysis of the experimental results that are shown in Figure 4 and Table 2, the area object matching method based on hierarchical matching can not only achieve high accuracy and stability, but is also superior to the weighted algorithm when it comes to integrity and accuracy of matching.From the theoretical point of view, there are cases where some indices of the weighted matching method are too large or too small, but the final weighted results satisfy the set threshold, which was also verified by the experiment.An further analysis was carried out with MD1_32 (meaning that the object was in the old-version candidate dataset MD1 and the ID was 32, and so on) and MD2_30 (meaning that the object was in the new-version candidate dataset MD2 and the ID was 30, and so on), MD1_33 and MD2_32, and MD1_65 and MD2_47 as the cases, as shown in Table 3, the analysis indicates that when two indices among the three indices are large while the other is smaller or two results are small and the other is larger, the weighted matching method can match two objects with each other, while the hierarchical matching method does not.In fact, these objects are mismatched.The first two pairs of matching objects contain an object with shape mutation, and there is a significant difference in positional distance when it comes to the latter as shown in Figure 5.This reason is that, the results obtained by weighted matching cannot reliably satisfy the overall matching probability of the geometric features of the area objects, when there are objects that match only some geometric features.The above analysis not only explains the reason why the hierarchical matching method is superior to the weighted matching method, but also points out the shortcomings of the weighted method.
To verify the reliability of the method, other four datasets of different study areas were selected to carry out the matching experiment.As shown in Table 4, as the number of matching objects increases in the study area, the integrity and accuracy of the hierarchical matching method and the weighted matching method remain within a certain range, and the result of the hierarchical matching method is better than that of the weighted algorithm in both integrity and accuracy.This indicates that the two methods both have good matching ability; however, the hierarchical matching method is better than the weighted matching method.An further analysis was carried out with MD1_32 (meaning that the object was in the old-version candidate dataset MD1 and the ID was 32, and so on) and MD2_30 (meaning that the object was in the new-version candidate dataset MD2 and the ID was 30, and so on), MD1_33 and MD2_32, and MD1_65 and MD2_47 as the cases, as shown in Table 3, the analysis indicates that when two indices among the three indices are large while the other is smaller or two results are small and the other is larger, the weighted matching method can match two objects with each other, while the hierarchical matching method does not.In fact, these objects are mismatched.The first two pairs of matching objects contain an object with shape mutation, and there is a significant difference in positional distance when it comes to the latter as shown in Figure 5.This reason is that, the results obtained by weighted matching cannot reliably satisfy the overall matching probability of the geometric features of the area objects, when there are objects that match only some geometric features.The above analysis not only explains the reason why the hierarchical matching method is superior to the weighted matching method, but also points out the shortcomings of the weighted method.To verify the reliability of the method, other four datasets of different study areas were selected to carry out the matching experiment.As shown in Table 4, as the number of matching objects increases in the study area, the integrity and accuracy of the hierarchical matching method and the weighted matching method remain within a certain range, and the result of the hierarchical matching method is better than that of the weighted algorithm in both integrity and accuracy.This indicates that the two methods both have good matching ability; however, the hierarchical matching method is better than the weighted matching method.

Detection of Change Types
The detection results obtained by the algorithm were compared with the actual types of changes, accompanied by a precision analysis.As implied by the data samples in Figure 6 and the results in Table 5, except that the accuracy and integrity to discriminate rotation and deformation are maintained between 89% and 93%, the accuracies of detecting the other changes mostly reach 100%.

Detection of Change Types
The detection results obtained by the algorithm were compared with the actual types of changes, accompanied by a precision analysis.As implied by the data samples in Figure 6 and the results in Table 5, except that the accuracy and integrity to discriminate rotation and deformation are maintained between 89% and 93%, the accuracies of detecting the other changes mostly reach 100%.Similarly, in order to verify the reliability and robustness of the method, other four datasets from different scales were also selected to carry out the type-of-change detection experiment.As shown in Table 6, the results of which show that both the integrity (R i ) and accuracy (R a ) of this method are relatively high, reaching above 90%, indicating that the detection method proposed in this study can accurately discriminate the types of changes.

Discussion and Conclusions
The detection and extraction of information about object change provides the basis for incrementally updating spatial databases.Considering spatial objects' change process, complete information about object change includes extraction of the changed objects and detection of their types of changes.To synthetically detect the change types of spatial area objects while extracting incremental information, so as to track object' spatiotemporal change pattern and provide the change backtracking mechanism for incrementally updating spatial datasets, as well as to solve problems like difficulty in detecting the types of changes, high complexity of detection methods, high redundancy rate of detection factors, and low degree of automation during the identification process during the incremental update of the spatial database, we proposed a method based on hierarchical matching model to detect change types of area objects in an integrated way.Using hierarchical matching operators reflecting different geometric features of area objects, this method filtered all the operators according to their complexity from low to high, creating hierarchical matching rules to extract incremental information.On this basis, we established the rules to judge the types of changes of objects based on the relevant matching operators, and designed the process to detect the nine types of changes of area objects, while minimizing the complexity of the algorithm and the redundancy rate of detection factors.Considering to track spatial object's change process, we discussed the identification, extraction, and database entry of the types of changes of area objects during the incremental information extraction process, achieved a close connection and organic coupling of incremental information extraction and object type-of-change detection, and ensured the accuracy and integrity of incremental information detection.The experimental results show that this method can successfully detect incremental information about area objects in practical applications, with the overall accuracy reaching above 90%, which is much higher than the existing matching method, making it quite feasible and applicable.It helps establish the corresponding relation between new and old objects, and facilitates the management of backtracking historical data and maintenance of spatial data.
However, in this study, only the detection on the types of changes of simple area objects at the same scale was explored.It is necessary to further study adaptive thresholds in the algorithm implementation, finding a method to examine the types of changes of complex objects under different scales, and more and varied test datasets are needed to demonstrate the robustness of the method in this study.Meanwhile, due to the limits of objective conditions and academic ability, we have not compared it with other methods except for the existing weighted matching method, nor performed with non-convex polygons or multi-part ones, it is inevitable for this paper that sensitivity to these issues needs to be addressed as some omissions or shortcoming during the study.Moreover, due to uncertainty and maybe fuzzy definitions, there may be an intersection with another object of the same kind even if an object has vanished.Therefore, the existing definitions of rules on extracting incremental information are not robust enough for real-world data, which should be addressed to be further studied in the future.

Figure 1 .
Figure 1.The incremental information extraction process.

Figure 1 .
Figure 1.The incremental information extraction process.

Figure 2 .
Figure 2. Illustration of extracting the change information based on hierarchical feature matching.

Figure 2 .
Figure 2. Illustration of extracting the change information based on hierarchical feature matching.

Figure 3 .
Figure 3. Illustration of type-of-change detection.Figure 3. Illustration of type-of-change detection.

Figure 3 .
Figure 3. Illustration of type-of-change detection.Figure 3. Illustration of type-of-change detection.
∀F n ∈ ChangeCollection (N), ∃F o ∈ ChangeCollection (O), if (F o ∩F n = Ø) then ChangeType (F o ) ← Vanish.Rule 2: If there is no intersection between the area object F n and any object F o within the same geographical theme, the change of F o is classified as Vanish.∀F o ∈ ChangeCollection (O), ∃F n ∈ ChangeCollection (N), if (F n ∩F o = Ø) then ChangeType (F 0 ) ← Vanish.Rule 3: If F p and F n are matched, and the types of changes of F p and F n are Vanish and Appearance, then the change of F n is classified as Reappearance.if ((Matching (F o , F n ) = Ture)and (ChangeType (F o ) = Vanish) and (ChangeType (F n ) = Appearance)) then ChangeType (F n ) ← Reappearance.Rule 4: If geometrical features of F o and F n are determined to be unchanged and the attribute information has changed, then the change of the object is classified as AttributeChange.if ((IsAttributeMatch = False) and (PositionResult ≥ ϕ 1 ) and (AreaResult ≥ ϕ 2 ) and (ShapeDirectionResult ≥ ϕ 3 )) then ChangeType (F o → F n ) ← AttributeChange.Rule 5: If area geometrical features of F n and F o are determined to be unchanged and the area of F n decreases from that of F o , while the other discriminant indexes are unchanged, then the change of the object is classified as Shrink.if ((IsAttributeMatch = True) and (PositionResult ≥ ϕ 1 ) and (AreaResult < ϕ 2 ) and (ShapeDirectionResult ≥ ϕ 3 ) and (S(F n ) < S(F o ))) then ChangeType (F o → F n ) ← Shrink.Rule 6: If geometrical features of the areas of F n and F o are determined to be changed and the area of F n increases from that of F o , while the other discriminant indexes are unchanged, then the change of the object is classified as Expansion.if ((IsAttributeMatch = True) and (PositionResult ≥ ϕ 1 ) and (AreaResult < ϕ 2 ) and (ShapeDirectionResult ≥ ϕ 3 ) and (S(F n ) > S(F o ))) then ChangeType (F o → F n ) ← Expansion.Rule 7: If geometrical features of the positions of F n and F o are determined to be changed while the other discriminant indexes are unchanged, then the change of the object is classified as Translation.if ((IsAttributeMatch = True) and (PositionResult < ϕ 1 ) and (AreaResult ≥ ϕ 2 ) and (ShapeDirectionResult ≥ ϕ 3 )) then ChangeType (F o → F n ) ← Translation.Rule 8: If geometrical features of the directions of F n and F o are determined to be changed while the other discriminant indexes are deemed as unchanged, then the types of changes of the object is classified as Rotation.if ((IsAttributeMatch = True) and (PositionResult ≥ ϕ 1 ) and (AreaResult ≥ ϕ 2 ) and (ShapeDirectionResult < ϕ 3 ) and (AssistResult ≥ ϕ 4 )) then ChangeType (F o → F n ) ← Rotation.Rule 9: The types of changes other than the types of changes stated in Rule 1-Rule 8 are defined as Deformation.∀E 1 ∈ ChangeCollection (O), F n ∈ ChangeCollection (N), if ((ChangeType / ∈ (Rule 1∪Rule 2∪Rule 3∪Rule 4∪Rule 5∪Rule 6∪Rule 7∪Rule 8)) then ChangeType (F o → F n ) ← Deformation.

Figure 4 .
Figure 4. Illustration of hierarchical matching result (a) and weighted matching result (b), respectively (the highlighted represents matched objects).

Figure 4 .
Figure 4. Illustration of hierarchical matching result (a) and weighted matching result (b), respectively (the highlighted represents matched objects).

Figure 5 .
Figure 5. Case comparison of matching objects from hierarchical matching and weighted algorithm.

Figure 5 .
Figure 5. Case comparison of matching objects from hierarchical matching and weighted algorithm.

Figure 6 .
Figure 6.Illustration of change samples from both old and new dataset.

Figure 6 .
Figure 6.Illustration of change samples from both old and new dataset.

Table 1 .
Comparison of calculation times of sub-operators from two algorithms.

Table 2 .
Comparison of experimental results from study area.

Table 3 .
Matching analysis of object instances.

Table 3 .
Matching analysis of object instances.
area O angle O

Table 4 .
Matching results in different study areas. No.

of Study Area Number of Samples Hierarchical Matching Weighted Matching R
i (%) R a (%) R i (%) R a (%)

Table 4 .
Matching results in different study areas.

Table 5 .
Detection result of change type in study area.

Table 5 .
Detection result of change type in study area.

Table 6 .
Detection results of change type in different study areas.