Multi-node cooperative construction land approval process whole life cycle tracing system

Abstract

This invention relates to the field of approval process traceability technology, and discloses a multi-node collaborative full lifecycle traceability system for construction land approval processes, including: constructing candidate traceability paths; calculating boundary inheritance coefficients; calculating boundary variation density; calculating the age of valid evidence; calculating boundary inheritance shear potential; and outputting full lifecycle traceability results. This invention transforms each processing node in the entire construction land approval process into event object nodes containing spatial objects, completion times, and evidence formation times. Candidate traceability paths are constructed based on spatial coverage relationships and temporal order. Combined with the recursive calculation of boundary inheritance coefficients, boundary variation density, and the age of valid evidence, path-level quantitative evaluation indicators are formed. Finally, the main traceability chain is selected, and a unified full lifecycle traceability result is output. This improves the uniformity of judgment standards during cross-departmental review processes and adapts to the practical application needs of digital supervision of natural resources.

Description

Technical Field

[0001] This invention relates to the field of approval process traceability technology, and more specifically, to a multi-node collaborative construction land approval process full lifecycle traceability system. Background Technology

[0002] With the continuous digital transformation of natural resource management and engineering construction project approval, and the deepening reforms of integrating multiple reviews and certificates, the approval of construction land is no longer a linear process within a single department. Instead, it is a collaborative process spanning multiple stages, including pre-approval and site selection, agricultural land conversion, planning permits, land supply, surveying, and registration. The same project corresponds to different business records, spatial graphics, and evidentiary materials at different stages. The approval chain is constantly lengthening, and data sources are continuously increasing. Full lifecycle process traceability has become the core foundation for natural resource supervision, compliance review, and law enforcement auditing. The industry is placing higher demands on the accuracy and automation of approval traceability.

[0003] Throughout the entire approval process, the boundaries of the same construction land object are not always constant. The site selection area, approval area, land parcels, and registration units formed in the previous stage may largely overlap, but partial adjustments may also occur. Simultaneously, the remote sensing imagery, land use control data, and survey and demarcation results cited at each node are not generated at the same time. This means that while subsequent nodes may technically inherit from previous nodes, they may not be completely synchronized in terms of spatial boundaries and the temporality of evidence, making it difficult to accurately determine the true lifecycle in retrospect. Existing traceability technologies generally use project numbers, certificate numbers, or chronological order to linearly connect approval records, focusing only on the temporal flow of the process and document associations, without incorporating changes in spatial boundaries and the temporal transmission of evidence into the core traceability system.

[0004] Existing technologies can only meet basic information query needs and struggle to simultaneously handle the combined effects of spatial coverage, boundary inheritance, morphological changes, and evidence aging. They also cannot perform unified quantitative comparisons of multiple possible traceability links. When the same project undergoes boundary refinement, partial replacement, or multi-path integration at multiple stages, existing methods often only yield a few scattered records, failing to automatically identify the main traceability chain or output a single, quantitative, full-lifecycle traceability result. This results in traceability verification still requiring significant manual intervention, leading to low efficiency and poor consistency in cross-departmental reviews, making it difficult to meet the practical application needs of intelligent supervision of the entire natural resource process. Summary of the Invention

[0005] This invention provides a multi-node collaborative construction land approval process full life cycle traceability system to solve the technical problems mentioned in the background.

[0006] This invention provides a multi-node collaborative system for tracing the entire lifecycle of construction land approval processes, including: The first module represents the approval node as an event object node containing spatial objects, completion time, and evidence formation time. For two event object nodes with sequential completion times, the coverage coefficient between the corresponding spatial objects is calculated. If the coverage coefficient is greater than zero, a directed edge is established from the event object node with the earlier completion time to the event object node with the later completion time. All directed edges are connected to form a candidate traceability path. The second module calculates the boundary inheritance coefficient between the spatial objects corresponding to the event object nodes at both ends of the directed edge. The third module calculates the boundary variation density using the corresponding spatial objects of the event object nodes at both ends of the directed edge; The fourth module extracts the completion time and evidence formation time of the event object node along the candidate tracing path, and calculates the effective evidence age of the event object node by recursively calculating the boundary inheritance coefficient according to the node. The fifth module multiplies the boundary variation density by the age of valid evidence to obtain the local shearing amount, calculates the time weight corresponding to the directed edge in combination with the completion time, weights the time weight with the local shearing amount and accumulates them along the candidate tracing path to obtain the boundary inheritance shearing potential. The sixth module maps the boundary inherited shear potential to traceability values, selects the candidate traceability path with the largest traceability value as the main traceability chain, and outputs the full lifecycle traceability result composed of the main traceability chain and the traceability value.

[0007] The beneficial effects of this invention are as follows: This invention transforms each processing node in the entire process of construction land approval into event object nodes that include spatial objects, completion times, and evidence formation times. Candidate tracing paths are constructed based on spatial coverage relationships and temporal order. Combined with recursive calculations of boundary inheritance coefficients, boundary variation density, and the age of valid evidence, path-level quantitative evaluation indicators are formed. Finally, the main tracing chain is selected, and a unified full-lifecycle tracing result is output. This invention is adaptable to construction land approval business scenarios involving multiple departments and nodes, taking into account spatial boundary changes and cross-node transmission of evidence temporality in the approval process. It provides a unified quantitative judgment basis for tracing the approval chain, reduces the workload of manual verification, improves the uniformity of judgment standards in cross-departmental review processes, and meets the practical application needs of digital supervision of natural resources. Attached Figure Description

[0008] Figure 1 This is a calculation flowchart of the multi-node collaborative construction land approval process full life cycle traceability system of the present invention. Detailed Implementation

[0009] The subject matter described herein will now be discussed with reference to exemplary embodiments. It should be understood that these embodiments are discussed only to enable those skilled in the art to better understand and implement the subject matter described herein, and changes may be made to the function and arrangement of the elements discussed without departing from the scope of this specification. Various processes or components may be omitted, substituted, or added as needed in the examples. Furthermore, features described in some examples may be combined in other examples.

[0010] It should be noted that, unless otherwise defined, the technical or scientific terms used in one or more embodiments of the present invention should have the ordinary meaning understood by one of ordinary skill in the art to which this invention pertains. The terms "first," "second," and similar terms used in one or more embodiments of the present invention do not indicate any order, quantity, or importance, but are merely used to distinguish different components. Terms such as "comprising" or "including" indicate that the element or object preceding the term encompasses the elements or objects listed following the term and their equivalents, without excluding other elements or objects. Terms such as "connected" or "linked" are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. Terms such as "upper," "lower," "left," and "right" are used only to indicate relative positional relationships; when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0011] like Figure 1 As shown, the multi-node collaborative construction land approval process full lifecycle traceability system includes: The first module represents the approval node as an event object node containing spatial objects, completion time, and evidence formation time. For two event object nodes with sequential completion times, the coverage coefficient between the corresponding spatial objects is calculated. If the coverage coefficient is greater than zero, a directed edge is established from the event object node with the earlier completion time to the event object node with the later completion time. All directed edges are connected to form a candidate traceability path. The second module calculates the boundary inheritance coefficient between the spatial objects corresponding to the event object nodes at both ends of the directed edge. The third module calculates the boundary variation density using the corresponding spatial objects of the event object nodes at both ends of the directed edge; The fourth module extracts the completion time and evidence formation time of the event object node along the candidate tracing path, and calculates the effective evidence age of the event object node by recursively calculating the boundary inheritance coefficient according to the node. The fifth module multiplies the boundary variation density by the age of valid evidence to obtain the local shearing amount, calculates the time weight corresponding to the directed edge in combination with the completion time, weights the time weight with the local shearing amount and accumulates them along the candidate tracing path to obtain the boundary inheritance shearing potential. The sixth module maps the boundary inherited shear potential to traceability values, selects the candidate traceability path with the largest traceability value as the main traceability chain, and outputs the full lifecycle traceability result composed of the main traceability chain and the traceability value.

[0012] In one embodiment of the present invention, the approval node is represented as an event object node containing a spatial object, a completion time, and an evidence formation time. For two event object nodes with sequential completion times, a coverage coefficient between their corresponding spatial objects is calculated. If the coverage coefficient is greater than zero, a directed edge is established from the event object node with the earlier completion time to the event object node with the later completion time. All directed edges are connected to form a candidate tracing path, including: The definition of an event object node is as follows: ; In the formula, For event object nodes, The space object owned by the event object node. The completion time of the event object node. The time when the evidence was formed for the event object node; The formula for calculating the coverage factor is as follows: ; In the formula, For coverage factor, For area operation operators, To complete the space object owned by the earliest event object node, To complete the space object owned by the event object node with the later time, The intersecting region is generated when the spatial objects of the event object node that completed earlier overlap with the spatial objects of the event object node that completed later. exist Under the condition of establishing a directed edge from the event object node with the earlier completion time to the event object node with the later completion time. ; by all directed edges Connections constitute candidate tracing paths .

[0013] It should be noted that an approval node is a business processing unit in the construction land approval process that can be independently identified, processed, and whose status changes can be independently recorded. It can be collected through the engineering construction project approval management system, the natural resources approval business system, the electronic certificate system, and the item flow log. An event object node is a structured computation node that is a unified transformation of the approval node. It is used to bind spatial objects, completion times, and evidence formation times to the same record, serving as the basic carrier for subsequent edge construction, recursion, and path filtering. Spatial objects are spatial geometric entities such as the land parcel scope, site selection scope, land supply scope, surveying scope, or registration scope corresponding to the approval node. They can be collected through vector graphics, surveying and demarcation results, joint surveying results, and real estate registration graphic results in the national land spatial information platform. The completion time is the time marker for the approval node to complete acceptance, review, approval, permitting, surveying confirmation, or registration. It describes the node's temporal position throughout its entire lifecycle and can be collected through the approval system's completion time, approval issuance time, electronic certificate generation time, and registration time. Evidence formation time is a time marker for the formation, collection, updating, or version solidification of spatial evidence on which the approval node relies. It is used to measure the age of the evidence cited at that node and can be collected through remote sensing image metadata, land use control database version records, surveying and mapping results entry records, and registration results generation records.

[0014] It should be noted that the intersecting region is the area where spatial objects with earlier completion times and those with later completion times overlap on the same coordinate system. It is used to determine whether the spatial extent of the preceding node covers the following node. The area of ​​the intersecting region is the numerical value of the area occupied by the intersecting region, representing the degree of actual spatial overlap between the two spatial objects. The complete coverage area is the numerical value of the entire coverage area of ​​the spatial object with the later completion time, used as the denominator of the coverage coefficient and reflecting the overall scale of the following node object. The coverage coefficient is the proportional value obtained by dividing the area of ​​the intersecting region by the complete coverage area, used to determine the degree to which the preceding node spatial object covers the following node spatial object. A directed edge is the directional connection relationship between event object nodes with earlier completion times and event object nodes with later completion times, used to represent node relationships with temporal sequence and spatial continuity in the lifecycle. A candidate traceability path is the set of paths formed by continuously connecting all directed edges in their directions, used to carry all possible traceability main lines from the front-end approval node to the back-end registration node.

[0015] It should be noted that since the approval of construction land use is not a single document flow, but involves changes in spatial scope, processing time, and evidence version simultaneously, representing approval nodes as event object nodes that include spatial objects, completion time, and evidence formation time allows business actions, spatial geometry, and evidence temporality to be uniformly incorporated into the same calculation unit. This enables subsequent continuous execution of edge construction, recursion, and path selection within the same node system. Since a traceable succession relationship can only be formed between two approval nodes if they are sequential in time and have real spatial overlap, establishing directed edges from the event object node with earlier completion time to the event object node with later completion time, under the condition that the coverage coefficient is greater than zero, excludes nodes that are purely temporally adjacent but spatially unrelated, thus limiting the traceability relationship to node pairs with a spatial inheritance basis. Furthermore, since the entire process of construction land use approval often involves multiple intermediate nodes and multiple possible succession routes, connecting all successfully established directed edges according to direction to form candidate traceability paths allows all possible links from front-end approval to back-end registration to be retained at once, providing a complete set of paths for subsequent comparison of the cumulative degree of evidence aging and boundary variation of different paths.

[0016] It should be noted that the extraction method for approval nodes is as follows: First, the item number, item name, processing status, completion time, associated project number, associated land parcel number, and evidence document number are extracted from the construction land approval business system, engineering construction project approval management system, electronic certificate system, and registration system. Then, these are merged according to the project number and spatial object number, removing withdrawn, invalid, and test records, and retaining valid processing records that can enter the lifecycle chain, ultimately generating an approval node list. The construction method for event object nodes is as follows: Based on a single approval node, the spatial object graphic, completion time, and evidence formation time bound to that node are called, and the three types of fields are written into the same structured record; when an approval node corresponds to multiple spatial objects, an event object node is formed for each spatial object; when multiple approval nodes correspond to the same spatial object but have different completion times or different evidence formation times, they are retained as different event object nodes. The method for determining the order of completion times is as follows: First, unify the time zone, date format, and second-level precision of all time fields, and then compare the completion times of two event object nodes; if the earlier time is earlier than the later time, it is determined that there is an order; if the two times are the same, no order relationship is established; if a node is missing a completion time, it is first supplemented based on the approval and issuance time, certificate generation time, or system completion time, and if it still cannot be supplemented, the node is removed.

[0017] It should be noted that the method for extracting intersecting regions is as follows: First, transform the two spatial objects to the same projected coordinate system, perform geometric validity repair and gap elimination, and then superimpose the two face objects to find their intersection. If the intersection result is multiple discrete faces, they are merged into the same set of intersecting regions; if the intersection result is empty, it is recorded as a non-intersecting region. The method for extracting the complete coverage area is as follows: Use the spatial objects owned by the event object node with the later completion time as the area calculation object, and calculate the area covered by the entire outer contour of the object under the unified projected coordinate system; when the object contains multiple discrete faces, the areas of each discrete face are accumulated; when the object has holes, only the effective coverage part is calculated, and the area of ​​the holes is not included in the complete coverage area. The method for establishing directed edges is as follows: For any two event object nodes whose time order has been determined, first calculate the coverage coefficient; when the coverage coefficient is greater than 0, generate a directed edge in the graph structure with the starting point of the event object node with the earlier completion time and the ending point of the event object node with the later completion time, and write the coverage coefficient into the edge attribute; when the coverage coefficient is equal to 0, no directed edge is generated. The candidate tracing path is constructed as follows: the event object node with an in-degree of 0 is used as the path source node, and the event object node with no successor directed edge is used as the path end node. A depth-first traversal or a breadth-first traversal is performed according to the direction of the directed edge to enumerate all directed paths between the source node and the end node; each non-repeating sequence of directed nodes is recorded as a candidate tracing path.

[0018] In one embodiment of the present invention, calculating the boundary inheritance coefficient between the spatial objects corresponding to the event object nodes at both ends of a directed edge includes: ; In the formula, For boundary inheritance coefficients, For length operation operators, To define the spatial object boundaries of the earliest event object node, To define the spatial object boundaries for event object nodes that occur later in time, The boundary sharing segment is formed by the overlap between the spatial object boundaries of the event object node that was completed earlier and the spatial object boundaries of the event object node that was completed later.

[0019] It should be noted that the spatial object boundary is the set of outer contours of the spatial object and inner contour lines that are legally required to be included in the boundary calculation, used to describe the linear range of the object's geometry at the edge position. The boundary shared segment is the set of line segments formed by the overlapping geometric contours of the spatial object boundaries completed earlier and later, used to describe the portion of the later node boundary inherited from the earlier node boundary. The boundary shared segment length is the cumulative value of all line segment lengths in the boundary shared segment, used to measure the actual scale of the inherited boundary portion. The total boundary length is the cumulative value of all line segment lengths of the spatial object boundary completed later, used as the denominator of the boundary inheritance coefficient. The boundary inheritance coefficient is a proportional value obtained by dividing the boundary shared segment length by the total boundary length, used to measure the degree of inheritance of the earlier spatial object boundary by the later-completed spatial object boundary.

[0020] It should be noted that, since the inheritance relationship in the approval of construction land mainly reflects whether the boundary is continued by subsequent nodes, rather than simply whether the area still overlaps, extracting the boundary-shared segment formed by the geometrical overlap of spatial object boundaries can identify the truly inherited boundary segments from the overall boundary, transforming boundary inheritance from an area concept to a linear geometric concept. Because the degree of inheritance of subsequent node boundaries from preceding node boundaries needs to be expressed by a comparable proportion, dividing the length of the boundary-shared segment by the total boundary length to calculate the boundary inheritance coefficient allows for comparison of boundary objects of different scales and lengths under the same proportional benchmark, thus reflecting how much of the subsequent node boundary directly originates from the preceding node boundary.

[0021] It should be noted that the method for extracting the boundaries of spatial objects is as follows: The outer contour lines are extracted from the planar geometry of the spatial object, and whether to include hole boundaries is determined according to implementation needs; for spatial objects composed of multiple discrete surfaces, the outer contour lines of each discrete surface are extracted separately and then merged into a boundary line set; after extraction, line topology cleanup is performed to eliminate duplicate line segments and hanging short lines. The method for determining shared boundary segments is as follows: the boundaries of spatial objects with earlier and later completion times are converted to the same coordinate reference, a unified line coincidence tolerance is set, and then the intersection of the two sets of boundaries is performed; the set of coincident line segments obtained from the intersection is the shared boundary segment; when two boundaries are only locally close and do not meet the tolerance requirements, they are not considered a shared boundary segment. The method for calculating the total boundary length is as follows: the length of all valid line segments in the boundary line set of the spatial object with the later completion time is calculated segment by segment and accumulated; when the spatial object is composed of multiple discrete surfaces, the boundary lengths of each discrete surface are summed separately; when there are hole boundaries on the boundary and the implementation plan requires them to be included in the calculation, the length of the hole boundary is also included in the total boundary length.

[0022] In one embodiment of the present invention, the boundary variation density is calculated using the corresponding spatial objects of the event object nodes at both ends of a directed edge, including: ; In the formula, For boundary variation density, For area operation operators, The symmetrical difference region is formed by the combination of spatial objects owned by event object nodes with earlier completion times and spatial objects owned by event object nodes with later completion times. The union region is formed by combining the spatial objects owned by the event object node with the event object node with the later completion time. For length operation operators, To define the spatial object boundaries for event object nodes that occur later in time, Pi is a constant. The space object owned by the event object node that completes later.

[0023] It should be noted that the symmetry difference region is the set of non-overlapping areas formed by merging spatial objects completed earlier and later, used to characterize the areas where the actual spatial forms changed between the two. The area of ​​the symmetry difference region is the cumulative value of the total area of ​​the symmetry difference region, used to quantify the scale of area change between the two spatial objects. The union region is the overall coverage area formed by combining the spatial objects completed earlier and later, used to represent the total spatial range after the merger. The area of ​​the union region is the area value of the union region, used as the denominator of the area change ratio. The area change ratio is the proportional value obtained by dividing the area of ​​the symmetry difference region by the area of ​​the union region, used to reflect the overall degree of change between the two spatial objects. The total coverage area is the area value of the spatial object completed later, used for subsequent conversion to the circumference of an ideal circle of equal area. The pi constant is a fixed mathematical constant used to calculate the circumference of an ideal circle of equal area from the area, with a preferred value of 3.1415926. This value already meets the accuracy requirements for length and area conversion in the spatial calculation of construction land approval and is consistent with the default calculation accuracy of commonly used spatial databases and surveying software. The square root value is an intermediate calculated value obtained by multiplying the total area covered by the constant pi and then taking the square root. It is used to convert the half-circumference of an ideal circle with the same area as the total area covered. The constant 2 is a fixed multiplier used when converting the square root value to the circumference of an ideal circle with the same area; it is 2. The conversion relationship between the circumference and area dictates that this multiplier must be fixed at 2. Using other values ​​will lead to distortion of the circumference of the ideal circle with the same area. The circumference of the ideal circle with the same area is the circumference value corresponding to the ideal circle with the same total area covered by the spatial object completed later. It is used as a standard reference for constructing the complexity of the boundary. The actual boundary length is the cumulative value of the true length of the boundary of the spatial object completed later in a unified coordinate system. It is used to reflect the actual curvature of the object's boundary. The shape irregularity is a proportional value obtained by dividing the actual boundary length by the circumference of the ideal circle with the same area. It is used to measure the complexity of the spatial object's boundary relative to a regular circle. The boundary variation density is a composite value obtained by multiplying the area change ratio by the shape irregularity. It is used to simultaneously express the scale of boundary change and the complexity of the boundary.

[0024] It should be noted that since the actual changes between spatial objects are not in the intersecting areas but in the non-overlapping parts, calculating the symmetry difference region and using the area of ​​the symmetry difference region and the area of ​​the union region to form the area change ratio can limit the area change to the actual spatial range of change, and normalize different object scales through the union region. Since the actual boundary length alone is insufficient to determine the complexity of a boundary, a rule-based reference corresponding to the area scale is needed. Therefore, converting the total coverage area into the circumference of an ideal circle of equal area provides a unified benchmark for subsequent comparisons of actual boundaries and ideal regular boundaries, allowing the boundary complexity of objects at different area scales to be measured using the same standard. Because the same area change has different impacts on the difficulty of tracing on smooth and zigzag boundaries, dividing the actual boundary length by the circumference of an ideal circle of equal area calculates the shape irregularity, which can transform factors such as boundary zigzags, fragmentation, and increased angles into a quantifiable indicator of complexity. Furthermore, since the scale of boundary changes and the complexity of the boundary jointly determine the strength of geometric disturbances in the tracing chain, multiplying the area change ratio by the shape irregularity to calculate the boundary variation density can simultaneously incorporate the magnitude and complexity of the change into a composite quantity, so that subsequent calculations no longer treat regular changes and complex changes as the same.

[0025] It should be noted that the calculation method for the symmetric difference region is as follows: First, perform a union operation on the spatial objects with earlier and later completion times, then remove the intersecting parts from the union result; the remaining region is the symmetric difference region. When the result consists of multiple discrete surfaces, all discrete surfaces are treated as the same set of symmetric difference regions. The calculation method for the union region is as follows: Perform a geometric merge on the spatial objects with earlier and later completion times in the same coordinate system to obtain a planar result that covers the entire range of both; if multiple discrete surfaces appear after the merge, all discrete surfaces together form the union region. The calculation method for the area change ratio is as follows: First, calculate the area of ​​the symmetric difference region, then calculate the area of ​​the union region; when the area of ​​the union region is greater than 0, divide the area of ​​the symmetric difference region by the area of ​​the union region to obtain the area change ratio; when the area of ​​the union region is equal to 0, record the area change ratio of the node pair as 0 and mark it as an abnormal geometric pair, which will not be used as a valid edge in the main chain selection. The calculation method for the circumference of an ideal circle with equal area is as follows: First, extract the total area covered by the spatial object completed later. Multiply the total area by the constant pi and take the square root. Then, multiply the result by the constant two to obtain the circumference of the ideal circle with the same total area. The calculation method for shape irregularity is as follows: First, calculate the actual boundary length of the spatial object completed later. Then, calculate the circumference of the ideal circle with the same area. Divide the actual boundary length by the circumference of the ideal circle with equal area. When the result equals 1, it indicates that the boundary is close to a regular circle. When the result is greater than 1, it indicates that the boundary has varying degrees of irregularity.

[0026] In one embodiment of the present invention, the completion time and evidence formation time of the event object node are extracted along the candidate tracing path, and the effective evidence age of the event object node is calculated recursively by node based on the boundary inheritance coefficient, including: The formula for calculating the valid evidence age of the starting event object node is as follows: ; In the formula, The age of valid evidence for the starting event object node. The completion time of the starting event object node. The time when evidence was formed for the starting event object node; The formula for calculating the effective evidence age of an event object node by recursively calculating the boundary inheritance coefficient is as follows: ; In the formula, To determine the valid evidence age for later event object nodes, To determine the completion time of event object nodes with later completion times, To complete the evidence formation time for later event object nodes, For boundary inheritance coefficients, The age of valid evidence for the earliest event object node.

[0027] It should be noted that the starting event object node is the first event object node at the source of the candidate tracing path network, used as the starting node for recursively calculating the age of valid evidence. The age of valid evidence is a comprehensive time difference that combines the time difference between the node's own completion time and evidence formation time with the aging amount transmitted from upstream nodes via the boundary inheritance coefficient; it describes the actual degree of evidence aging carried by the current node. The time difference of the current node is the time difference obtained by subtracting the evidence formation time of the event object node with a later completion time; it reflects the newness of the evidence at the current node. The inherited aging time difference is the time difference obtained by multiplying the boundary inheritance coefficient by the age of valid evidence of the event object node with an earlier completion time; it represents the amount of influence of evidence aging from upstream nodes transmitted to downstream nodes via boundary inheritance.

[0028] It should be noted that since subsequent nodes do not always regenerate all boundaries, but rather inherit a portion of the boundaries from preceding nodes, the effective evidence age of event object nodes can be calculated recursively by combining the boundary inheritance coefficient. This allows the aging of evidence from preceding nodes to be passed on to subsequent nodes through the inherited boundaries, thus reflecting the actual time lag of comprehensive evidence carried by the current node. Furthermore, since the impact of upstream node evidence aging on downstream nodes is not entirely transmitted, but rather proportionally according to the boundary inheritance ratio, multiplying the boundary inheritance coefficient by the effective evidence age of the event object node with the earliest completion time to calculate the inheritance aging time difference can limit the cross-node impact of evidence aging to the range of the inherited boundaries.

[0029] It should be noted that the starting event object node is determined as follows: within a single candidate tracing path, find the event object node with an in-degree of 0 located at the beginning of the path, and determine it as the starting event object node; if multiple source points enter the same successor node in parallel on a candidate tracing path, they are retained as different candidate tracing paths, and the starting event object node is determined separately in each path. The evidence formation time is selected according to the principle that whoever actually supports the spatial judgment of the node is the evidence time, selecting the time most directly corresponding to the node from the remote sensing image date, land use control data version time, surveying and mapping result generation time, surveying and demarcation result completion time, or registration result formation time; if multiple pieces of evidence are used for the same node, the evidence time with the latest formation time and actually written into the approval basis is taken. The recursive method for the effective evidence age is as follows: first, determine the starting event object node on each candidate tracing path, and subtract its evidence formation time from the completion time of the starting event object node to obtain the initial effective evidence age; then, recursively calculate along the path direction, using the current node's time difference plus the inheritance aging time difference to obtain the effective evidence age of subsequent nodes, until all nodes on the path have been calculated.

[0030] In one embodiment of the present invention, the local shearing amount is obtained by multiplying the boundary variation density by the age of valid evidence, and the time weight corresponding to the directed edge is calculated by combining the completion time. The time weight is weighted with the local shearing amount and accumulated along the candidate tracing path to obtain the boundary inheritance shearing potential, including: The formula for calculating local shear volume is as follows: ; In the formula, This is the local shear amount. For boundary variation density, To determine the age of valid evidence for event object nodes that are completed later; The formula for calculating time weight is as follows: ; In the formula, As time weight, To determine the completion time of event object nodes with later completion times, To determine the completion time of the event object node that was completed earlier, As candidate tracing paths, These are independent directed edges within the candidate tracing path. The completion time of the event object node for an independent directed edge terminal. The completion time of the event object node at the beginning of an independent directed edge; The formula for calculating the boundary inherited shear potential is as follows: ; In the formula, Inherited shear potential at the boundary, It is a directed edge.

[0031] It should be noted that the local shearing amount is a local intensity value formed by multiplying the boundary variation density by the age of valid evidence, used to describe the degree of combined effect of geometric changes and evidence aging on a single directed edge. The local interval time is the time difference obtained by subtracting the completion time of the earlier-completion event object node from the completion time of the later-completion event object node, used to represent the lifecycle duration corresponding to a single directed edge. The global total duration is the total duration obtained by summing all local interval times within the same candidate tracing path, used as a normalization benchmark for time weights. The time weight is a proportional value obtained by dividing the local interval time by the global total duration, used to characterize the time proportion of a certain directed edge in the entire candidate tracing path. The time-weighted local shearing amount is the result of multiplying the time weight by the local shearing amount, used to weight the local shearing amount according to the lifecycle duration. The boundary inheritance shear potential is a path-level value formed by summing all time-weighted local shearing amounts within the candidate tracing path, used to comprehensively represent the cumulative intensity of boundary inheritance, geometric variation, and evidence aging propagation on that path.

[0032] It should be noted that since the retrospective distortion on a single directed edge depends on both the strength of the boundary change and the aging of the evidence carried by that edge, multiplying the boundary variation density by the age of the valid evidence to obtain the local shearing amount can couple geometric change and temporal aging into the same local intensity quantity. Furthermore, since different directed edges have varying durations throughout the entire candidate retrospective path, the influence of long-term edges should not be treated equally with that of short-term edges. Therefore, weighting the time weight with the local shearing amount and accumulating it along the candidate retrospective path to obtain the boundary inherited shear potential can form a path-level index reflecting the cumulative distortion of the entire path.

[0033] It should be noted that the total global duration is calculated as follows: the local interval times of all directed edges within the same candidate tracing path are summed one by one to obtain the total duration of the path from the start to the end. When a path contains only one node and no directed edges, the total global duration is not calculated, and the path is excluded from the boundary inheritance shear potential calculation. The time weight is calculated as follows: first, the local interval time of a single directed edge is calculated, and then this local interval time is divided by the total global duration of the candidate tracing path; the result is recorded as the time weight of the directed edge; the sum of the time weights of all directed edges within the same path should equal 1. The accumulation method of the boundary inheritance shear potential is as follows: first, the local shear amount is calculated for each directed edge, then the local shear amount is multiplied by the time weight corresponding to the edge to obtain the time-weighted local shear amount, and finally, all time-weighted local shear amounts are accumulated according to the edge order of the candidate tracing path to obtain the boundary inheritance shear potential of a single candidate tracing path.

[0034] In one embodiment of the present invention, the boundary inherited shear potential is mapped to a traceability value, the candidate traceability path with the largest traceability value is selected as the main traceability chain, and the full lifecycle traceability result composed of the main traceability chain and the traceability value is output, including: The formula for calculating the traceability value is as follows: ; In the formula, For traceability value, For the natural base exponentiation operator, Inherit the shear potential at the boundary; It should be noted that, in order to eliminate the influence of the time dimension on the stability of the exponential operation, a characteristic time constant can be introduced for standardization. The resulting negative exponential parameter is a dimensionless input quantity obtained by dividing the boundary inherited shear potential by the characteristic time constant and assigning it a negative sign. This input quantity is used to map the cumulative intensity to the decaying traceability value. The traceability value is the result of performing a natural base exponential operation on this negative exponential parameter, and is used to quantify the traceability of the candidate traceability path. Specifically, the formula for calculating the traceability value is as follows: ,in This represents a preset characteristic time constant used for unifying the time dimension. It can be set according to the statutory maximum processing time limit for approval business (e.g., set to the same time value corresponding to 1 standard calendar month or 30 working days). The extraction formula for the main traceability chain includes: ; ; ; In the formula, The maximum traceability value indicates the path number. To find operators for maximal numerical numbering, Main traceability chain, The maximum traceability value indicates the candidate traceability path bound to the path number. The maximum traceability value, The maximum traceability value indicates the traceability value bound to the path number; The formula for synthesizing the full lifecycle traceability results is as follows: ; In the formula, This provides results for full lifecycle traceability.

[0035] It should be noted that the negative boundary inherited shear potential is the input quantity obtained by assigning a negative sign to the boundary inherited shear potential, used to map the cumulative intensity to the decaying traceability value. The traceability value is the result of performing a natural base exponential operation on the negative boundary inherited shear potential as the exponential parameter, used to quantify the traceability of candidate traceability paths. The maximum traceability value is the highest value retrieved from all candidate traceability path traceability values, used to determine the optimal traceability result. The maximum traceability value indicator path number is the path identifier corresponding to the maximum traceability value among all candidate traceability paths, used to locate the main traceability chain. The main traceability chain is the candidate traceability path bound to the maximum traceability value, used to represent the main lifecycle undertaking path finally confirmed by the system. The full lifecycle traceability result is the final output result formed by the combination of the main traceability chain and the maximum traceability value, used to provide a single and interpretable traceability conclusion to the review, supervision, and verification stages.

[0036] It should be noted that since a larger boundary inherited shear potential indicates a stronger cumulative path distortion, and the system ultimately needs to output a single-value result that is easy to compare, mapping the boundary inherited shear potential to a traceability value can convert the path-level cumulative intensity into a traceability evaluation quantity that is easier to understand with a larger numerical value. Since there is usually more than one candidate traceability path, the system must select the one that best represents the true lifecycle continuity relationship from all possible links. Therefore, selecting the candidate traceability path with the largest traceability value as the main traceability chain can transform path selection from empirical judgment to automatic selection under a unified numerical criterion. Furthermore, since the review, supervision, and verification stages need to know both which lifecycle link is ultimately identified and the traceability degree of that link, outputting a full lifecycle traceability result composed of the main traceability chain and the traceability value can simultaneously provide path conclusions and quantitative conclusions, facilitating subsequent review.

[0037] It should be noted that the mapping method for the traceability value is as follows: first, a negative sign is assigned to the boundary inheritance shear potential, and then the negative value is used as the exponent parameter to perform natural base exponentiation; the calculation result is the traceability value; the larger the boundary inheritance shear potential, the closer the traceability value is to 0; the smaller the boundary inheritance shear potential, the closer the traceability value is to 1. When the maximum traceability values ​​are tied, the method for determining the main traceability chain is as follows: first, compare the boundary inheritance shear potentials among the tied paths, and select the path with the smaller boundary inheritance shear potential; if the boundary inheritance shear potentials are still the same, then compare the cumulative coverage coefficient values, and select the path with the larger cumulative coverage coefficient value; if they are still the same, then select the path with the later completion time as the main traceability chain. The output method for the full lifecycle traceability results is as follows: the node sequence, directed edge sequence, maximum traceability value, boundary inheritance shear potential, and the effective evidence age of each node corresponding to the main traceability chain are combined into a unified result object, and simultaneously output as the system interface display result and the structured data result; the system interface display result is for manual viewing, and the structured data result is for auditing and subsequent review.

[0038] It should be noted that this invention is suitable for deployment within the existing engineering construction project approval management platforms, territorial spatial information platforms, and real estate registration business platforms of provincial, municipal, and county natural resources authorities. During deployment, a unified project identifier, a unified land parcel identifier, and a unified coordinate benchmark are first established at the data layer. Then, acceptance records, pre-approval and site selection records, agricultural land conversion approval records, construction land planning permit records, land supply parcel records, joint surveying records, and real estate registration records are connected to the same traceability database. After each business record enters the traceability database, the spatial object, completion time, and evidence formation time must be synchronously written. The spatial object can come from survey and demarcation results, parcel maps, joint surveying results, and real estate unit graphics. The completion time can come from the system completion time, approval issuance time, and registration entry time. The evidence formation time can come from the remote sensing image date, land use control data version time, surveying result generation time, and registration result generation time. After the system completes data access, it first automatically generates event object nodes, then automatically calculates the coverage coefficient and establishes directed edges, followed by calculating the boundary inheritance coefficient, boundary variation density, and effective evidence age, and finally forms the boundary inheritance shear potential and tracing value for all candidate tracing paths. Specifically, all data involved in this invention has been explicitly authorized and consented to by the user before collection, and strictly adheres to relevant regulations on data security and privacy protection.

[0039] It should be noted that in actual operation, staff do not need to manually compare the graphics node by node, nor do they need to rely on experience to judge which link is closer to the actual connection relationship. After receiving a query request for a construction land project, the system will first list all possible links from the candidate traceability paths, and then automatically filter out the main traceability chain. For example, an industrial project has 8 event object nodes, corresponding to acceptance, pre-approval and site selection, agricultural land conversion approval, provincial review, construction land planning permit, land supply, joint surveying and mapping, and registration. The system forms 3 candidate traceability paths. The boundary inheritance shear potential of the first path is 0.1824, and the traceability value is 0.8333. The boundary inheritance shear potential of the second path is 0.3271, and the traceability value is 0.7210. The boundary inheritance shear potential of the third path is 0.4516, and the traceability value is 0.6365. Since the first path has the highest traceability value, the system identifies the first path as the main traceability chain and outputs the node order, the age of valid evidence for each node, the inheritance coefficient of each edge, and the final maximum traceability value along with the path.

[0040] It should be noted that the final output results can be divided into two categories. The first category is a results page for human readers, which displays the node order of the main traceability chain, the spatial object corresponding to each node, the completion time, the evidence formation time, and the age of valid evidence, as well as the maximum traceability value. The second category is a structured result for system calls, which includes at least the main traceability chain identifier, node set, edge set, maximum traceability value, boundary inheritance shear potential, and timestamp. This result can be directly used for review and verification, internal audit, law enforcement evidence collection, project review, and historical responsibility tracing. Its practical application value lies in its ability to reorganize the previously scattered relationships across multiple systems, stages, and sets of graphs into a readable, calculable, and verifiable main chain, reducing the workload of manual verification, improving the consistency of cross-departmental verification, and providing a unified traceability basis for subsequent supervision, spot checks, and dispute resolution. After deployment, it can also be run in batches by project number, land parcel number, registration unit number, or administrative division to sort all projects within the same region by traceability value, prioritizing items with lower traceability values ​​to assist management departments in conducting spot checks and supplementary certification work, which will not be elaborated here.

[0041] It should be noted that the interval and threshold sizes are set for ease of comparison. The size of the threshold depends on the amount of sample data and the base number set by those skilled in the art for each set of sample data, as long as it does not affect the proportional relationship between the parameter and the quantized value. Furthermore, the above formulas are all dimensionless calculations, and the formulas are derived from software simulations using a large amount of collected data to obtain the most recent real-world results. The preset parameters in the formulas are set by those skilled in the art according to the actual situation.

[0042] The embodiments of this example have been described above. However, this example is not limited to the specific implementation methods described above. The specific implementation methods described above are merely illustrative and not restrictive. Those skilled in the art can make many other forms based on the guidance of this example, and all of them are within the protection scope of this example.

Claims

1. A multi-node cooperative construction land approval process whole life cycle tracing system, characterized in that, include: The first module represents the approval node as an event object node containing spatial objects, completion time, and evidence formation time. For two event object nodes with sequential completion times, the coverage coefficient between the corresponding spatial objects is calculated. If the coverage coefficient is greater than zero, a directed edge is established from the event object node with the earlier completion time to the event object node with the later completion time. All directed edges are connected to form a candidate traceability path. The second module calculates the boundary inheritance coefficient between the spatial objects corresponding to the event object nodes at both ends of the directed edge. The calculation formula of the boundary inheritance coefficient is as follows: ; In the formula, For boundary inheritance coefficients, For length operation operators, To define the spatial object boundaries of the earliest event object node, To define the spatial object boundaries for event object nodes that occur later in time, The boundary sharing segment is formed by the overlap between the spatial object boundaries of the event object node that was completed earlier and the spatial object boundaries of the event object node that was completed later; The third module calculates the boundary variation density using the corresponding spatial objects of the event object nodes at both ends of the directed edge; The calculation formula of the boundary variation density is as follows: ; In the formula, For boundary variation density, For area operation operators, The symmetrical difference region is formed by the combination of spatial objects owned by event object nodes with earlier completion times and spatial objects owned by event object nodes with later completion times. The union region is formed by combining the spatial objects owned by the event object node with the event object node with the later completion time; The fourth module extracts the completion time and evidence formation time of the event object node along the candidate tracing path, and calculates the effective evidence age of the event object node by recursively calculating the boundary inheritance coefficient according to the node. The formula for calculating the valid evidence age of the start event object node is as follows: ; In the formula, is the valid evidence age of the start event object node, is the completion time of the start event object node, is the evidence formation time of the start event object node; The calculation formula of the effective evidence age of the event object node according to the node recursive calculation of the combination boundary inheritance coefficient is as follows: ; In the formula, To determine the valid evidence age for later event object nodes, To determine the completion time of event object nodes with later completion times, To complete the evidence formation time for later event object nodes, To determine the age of valid evidence for event object nodes with earlier completion times; The fifth module multiplies the boundary variation density by the age of valid evidence to obtain the local shearing amount, calculates the time weight corresponding to the directed edge in combination with the completion time, weights the time weight with the local shearing amount and accumulates them along the candidate tracing path to obtain the boundary inheritance shearing potential. The sixth module maps the boundary inherited shear potential to a traceability value, selects the candidate traceability path with the largest traceability value as the main traceability chain, and outputs the full lifecycle traceability result composed of the main traceability chain and the traceability value.

2. The multi-node collaborative construction land approval process full lifecycle traceability system according to claim 1, characterized in that, Extract all approval nodes included in the construction land approval process and represent all approval nodes as event object nodes that include spatial objects, completion time, and evidence formation time; Find two event object nodes whose completion times are sequential; Get the space objects owned by the event object node with the earlier completion time and the space objects owned by the event object node with the later completion time; Extract the intersecting region formed by the spatial objects owned by the event object node with the earlier completion time and the spatial objects owned by the event object node with the later completion time, and calculate the area of ​​the intersecting region. 3.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 2, characterized in that, Extract the complete coverage area occupied by the space object of the event object node with the latest completion time; The coverage factor is calculated by dividing the area of ​​the intersecting region by the area of ​​the complete coverage. When the coverage coefficient has a numerical attribute greater than zero, establish a directed edge whose starting point belongs to the event object node with an earlier completion time and whose ending point belongs to the event object node with a later completion time. All successfully established directed edges are connected and combined according to their directions to form candidate tracing paths. 4.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 1, characterized in that, Extract the spatial object boundary of the event object node with the earlier completion time contained within the event object nodes at both ends of the directed edge; Extract the spatial object boundary of the event object node with the later completion time contained within the event object nodes at both ends of the directed edge; Extract the boundary shared segment formed by the overlapping geometric contours of the spatial object boundaries of the event object node with the earlier completion time and the event object node with the later completion time; Calculate the length of the boundary sharing segment; Calculate the total boundary length of the spatial object boundary owned by the event object node that completes later; The boundary inheritance coefficient is calculated by dividing the length of the shared boundary segment by the total boundary length. 5.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 1, characterized in that, Extract the space object owned by the event object node with the earlier completion time contained within the event object nodes at both ends of the directed edge; Extract the space objects owned by the event object node with the later completion time contained within the event object nodes at both ends of the directed edge; Calculate the symmetrical difference region formed by merging the spatial objects owned by the event object node with the earlier completion time and the spatial objects owned by the event object node with the later completion time, and calculate the area of ​​the symmetrical difference region. The calculation is performed by combining the spatial objects owned by the event object node with the spatial objects owned by the event object node with the later completion time. Calculate the area of ​​the union region, and divide the area of ​​the symmetric difference region by the area of ​​the union region to obtain the area change ratio. 6.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 5, characterized in that, Extract the total area covered by the space objects owned by the event object node with the latest completion time; Multiply the total area covered by pi and take the square root to generate the square root value. Multiply the square root value by a constant two to calculate the circumference of the ideal circle with equal area. Extract the actual boundary length of the spatial object boundary of the event object node with the later completion time; The shape irregularity is calculated by dividing the actual boundary length by the circumference of an ideal circle of equal area. The boundary variation density is calculated by multiplying the area change ratio by the shape irregularity. 7.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 1, characterized in that, Extract the starting event object node from the source location of the candidate tracing path network, and extract the completion time of the starting event object node and the evidence formation time of the starting event object node; The effective evidence age of the starting event object node is calculated by subtracting the evidence formation time of the starting event object node from the completion time of the starting event object node. Extract event object nodes with later completion times connected by directed edges along the candidate tracing path; Extract the completion time of event object nodes with later completion times and the evidence formation time of event object nodes with later completion times; The time difference of the current node is calculated by subtracting the evidence formation time of the event object node with the later completion time from the completion time of the event object node with the later completion time. 8.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 7, characterized in that, Extract the boundary inheritance coefficients of directed edge associations and the valid evidence age of event object nodes with earlier completion times; The inheritance aging time difference is calculated by multiplying the boundary inheritance coefficient by the effective evidence age of the event object node with the earlier completion time. Add the inheritance aging time difference to the current node's time difference to obtain the valid evidence age of the event object node with the later completion time; The valid evidence age of all event object nodes is calculated recursively for all nodes along the candidate tracing path. 9.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 1, characterized in that, Extract the boundary variation density attached to the directed edge, and extract the effective evidence age of the event object node with the later completion time pointed to by the end of the directed edge; multiply the boundary variation density by the effective evidence age to calculate the local shearing amount. Extract the completion time of event object nodes with later completion times and the completion time of event object nodes with earlier completion times; The local interval time is calculated by subtracting the completion time of the event object node with the earlier completion time from the completion time of the event object node with the later completion time. The total global duration is obtained by summing all the local interval times contained within the candidate tracing path. The time weight is calculated by dividing the local interval time by the total global duration. Multiply the time weight by the local shearing amount to obtain the time-weighted local shearing amount; The boundary inherited shear potential is calculated by summing up all the time-weighted local shear values ​​of the candidate tracing paths. 10.The multi-node coordinated construction land approval process whole life cycle tracing system according to claim 1, characterized in that, Extract the boundary inherited shear potential, assign a negative sign to generate a negative boundary inherited shear potential; Perform an exponential operation based on the natural base and with the negative boundary inherited shear potential as the exponential parameter to calculate the retrospective value; Search among all trace values ​​generated by the operation to find the largest trace value that occupies the state with the largest numerical value; The candidate traceability path bound to the maximum traceability value is retrieved to determine the main traceability chain; The combined assembly of the main traceability chain and the output of the maximum traceability value constitute a successful full lifecycle traceability result.

Images