Urban earthquake prevention one map adaptive multi-scale visualization mapping method and system

Abstract

The application discloses a city earthquake prevention one-map adaptive multi-scale visualization mapping method and system, relates to the technical field of visualization mapping, and comprises the following steps: S1, constructing an earthquake prevention multi-source basic data set, and preprocessing the earthquake prevention multi-source basic data; S2, constructing an alternate confirmation section according to a relative post-earthquake time, calculating an object alternate confirmation constraint value, and generating object-level earthquake prevention correlation data; S3, when there is a boundary contact relationship in a boundary region of adjacent building objects or an overlapping relationship in a suspected damage patch region, calculating a locking expansion index value, and generating cross-scale version locking data; and S4, constructing different hierarchical expression base elements, and performing expression replacement, fission mapping processing and time sequence locking display to generate a city earthquake prevention one-map adaptive multi-scale visualization map. The application solves the problems that in the prior art, multi-source earthquake prevention data in a city earthquake prevention one-map system is only used for layer superposition, it is difficult to accurately correlate the same building object, and there is a lack of reliable quantitative basis for multi-scale map switching.

Description

Technical Field

[0001] This invention relates to the field of visualization mapping technology, specifically to an adaptive multi-scale visualization mapping method and system for urban earthquake prevention. Background Technology

[0002] With the continuous development of digital twin cities, urban lifeline monitoring, remote sensing interpretation, IoT sensing, and multi-source geographic information fusion technologies, thematic mapping technologies for urban governance, disaster prevention and mitigation, and risk assessment are evolving from single-layer representation to multi-source data collaboration, cross-scale linkage, image synthesis, and dynamic visualization. Especially in urban earthquake prevention scenarios, constructing a unified base map around building objects, suspected post-earthquake damage patches, on-site sensing points, and manually reported information, and achieving multi-scale map organization and temporal representation, has become an important technological development direction for urban earthquake damage monitoring, risk assessment, and decision support.

[0003] For example, application CN121685759A provides an automated administrative boundary fusion mapping method and meteorological mapping method based on the Python language. It belongs to the field of meteorological mapping technology, specifically involving the interdisciplinary technology of meteorological data processing and geographic information system visualization. Starting with the construction of multi-scale adaptive boundary fusion mapping, it addresses the problems of accuracy discontinuities across scales (e.g., from province to village) and large interpolation errors in complex terrain. It achieves seamless aggregation technology across administrative boundaries at different levels (village level → township level → county level → province level), thereby solving the problem of accuracy loss during cross-scale transformations (e.g., from province to village) in meteorological mapping.

[0004] For example, application CN109801296B discloses a nested remote sensing refined mapping method for the underlying surface of sponge city construction. It introduces advanced concepts such as "low-impact development" and "sponge city," along with classification methods like "hierarchical nesting." It proposes a novel remote sensing refined mapping technology supported by a hierarchical nested classification system that integrates multi-source urban surface data adapted to sponge city construction. This technology develops key technologies for extracting multi-attribute information from "mapped entities"—including surface permeability, relative slope, population density levels, and urban functional zoning—to achieve sponge city construction planning, design, and engineering implementation goals. It also develops multi-scale mapping and visualization methods. This system can realize intelligent calculation and statistical mapping methods to serve the identification of sponge city construction goals.

[0005] With the continuous development of digital twin cities, urban lifeline monitoring, remote sensing interpretation, IoT sensing, and multi-source geographic information fusion technologies, thematic mapping technologies for urban governance, disaster prevention and mitigation, and risk assessment are evolving from single-layer representation to multi-source data collaboration, cross-scale linkage, image synthesis, and dynamic visualization. Especially in urban earthquake prevention scenarios, constructing a unified base map around building objects, suspected post-earthquake damage patches, on-site sensing points, and manually reported information, and achieving multi-scale map organization and temporal representation, has become an important technological development direction for urban earthquake damage monitoring, risk assessment, and decision support.

[0006] Therefore, in order to address the above issues, there is an urgent need for an adaptive multi-scale visualization mapping method and system for urban earthquake prevention. Summary of the Invention

[0007] Technical problems to be solved

[0008] To address the shortcomings of existing technologies, this invention provides an adaptive multi-scale visualization mapping method and system for urban earthquake prevention maps. This solves the problems in existing urban earthquake prevention map systems where multi-source earthquake prevention data are simply overlaid on layers, making it difficult to accurately correlate the same building object, and where there is a lack of reliable quantitative basis for switching between multi-scale maps.

[0009] Technical solution

[0010] To achieve the above objectives, this invention provides the following technical solution: an adaptive multi-scale visualization mapping method for urban earthquake prevention, comprising: S1, constructing a unified spatial reference coordinate system, collecting spatial and temporal data corresponding to the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations, constructing a multi-source earthquake prevention basic dataset, and preprocessing the multi-source earthquake prevention basic data; S2, based on the preprocessed multi-source earthquake prevention basic data, constructing building object boundary regions, building core regions, building boundary buffer zones, and suspected damaged patch regions, and determining the building corresponding to the current record by combining the overlapping area value and the location falling relationship. The survey objects are numbered, and alternating confirmation segments are constructed according to the relative post-earthquake time. The alternating confirmation constraint value of the objects is calculated to generate object-level seismic defense correlation data. S3, Based on the object-level seismic defense correlation data, the locking version kernel is determined. When there is a boundary contact relationship or an overlapping relationship between suspected damaged patches in the boundary area of ​​adjacent building objects, the locking expansion index value is calculated to generate cross-scale version locking data. S4, Based on the cross-scale version locking data, different levels of expression base elements are constructed. According to the maximum value of the alternating confirmation constraint value of the objects and the relative post-earthquake time, expression replacement, fission mapping processing and temporal locking display are performed to generate an adaptive multi-scale visualization map of urban seismic defense.

[0011] Furthermore, a unified spatial reference coordinate system is constructed, and spatial and temporal data corresponding to the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations are collected to build a multi-source basic dataset for earthquake prevention. The specific steps for preprocessing the multi-source basic dataset for earthquake prevention are as follows: Select the plane coordinate reference datum and elevation reference datum corresponding to the urban earthquake prevention area. Use the plane coordinate reference datum to determine the east and north coordinate axes, and use the elevation reference datum to determine the height reference axis, thus constructing a unified spatial reference coordinate system. Then, map the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations to the unified spatial reference coordinate system, and collect the east coordinate values ​​of the building base outline boundary points, the north coordinate values ​​of the building base outline boundary points, the east coordinate values ​​of the suspected damaged patch boundary points, the north coordinate values ​​of the suspected damaged patch boundary points, and the east coordinate values ​​of the IoT sensing device installation points. The system collects the north coordinates of the IoT sensing device installation points, the east coordinates of the manually reported locations, the north coordinates of the manually reported locations, and the building survey object number. Simultaneously, it collects the earthquake event occurrence time, remote sensing acquisition time, IoT sensing sampling time, manually reported time, total building height, total building area, building tilt angle, structural crack width, building settlement, and building horizontal displacement to construct a multi-source earthquake prevention basic dataset. For the collected multi-source earthquake prevention basic data, a precise time protocol master-slave clock deviation compensation algorithm is used to perform unified processing of multi-source acquisition times; a Bursa seven-parameter coordinate transformation algorithm is used to perform unified processing of spatial location benchmarks; a median absolute deviation anomaly detection algorithm is used to perform anomaly data identification processing; a sliding median filtering algorithm is used to perform local abrupt earthquake prevention record smoothing processing; and a range normalization algorithm is used to perform numerical scale unification processing, outputting the preprocessed multi-source earthquake prevention basic data.

[0012] Furthermore, based on the preprocessed earthquake-resistant multi-source basic data, the following steps are taken to construct the building object boundary region, building core region, building boundary buffer zone, and suspected damaged patch region, and to determine the building survey object number corresponding to the current record by combining the overlapping area value and the positional relationship: Read the preprocessed earthquake-resistant multi-source basic dataset; construct the building object boundary region according to the building survey object number value; determine the inward offset distance according to the product of the total building height value and the preset offset ratio; perform inward offset on the building object boundary region to generate the building core region; and connect the building object boundary region and the building core region... The area is defined as the building boundary buffer zone; the eastward coordinates and northward coordinates of the boundary points of suspected damaged patches are collected according to the same remote sensing acquisition time value, and then connected end to end in the order of the boundary points to construct the suspected damaged patch area; the overlap area value between each suspected damaged patch area and each building core area and the overlap area value between each building boundary buffer zone are calculated respectively; and the building survey object number value that satisfies the largest patch overlap area value and the IoT sensing device installation point or the manually reported location falls within the building core area or the building boundary buffer zone is determined as the building survey object number value corresponding to the current record.

[0013] Further, the specific steps for constructing alternating confirmation segments based on relative post-earthquake times and calculating the alternating confirmation constraint values ​​for objects are as follows: Using the earthquake event occurrence time as a benchmark, convert the remote sensing acquisition time value, IoT sensing sampling time value, and manual reporting time value into relative post-earthquake times; aggregate all records corresponding to the same building object according to the building survey object number value, and arrange them in ascending order of relative post-earthquake times; when adjacent records correspond to the same building survey object number value and the adjacent records contain three types of sources—remote sensing, IoT, and manual reporting—merge the adjacent records into an alternating confirmation segment, and calculate the first time difference between the remote sensing acquisition time value and the earthquake event occurrence time value, the second time difference between the IoT sensing sampling time value and the earthquake event occurrence time value, and the third time difference between the manual reporting time value and the earthquake event occurrence time value; for For each alternating confirmation segment, the first spatial constraint term is obtained by adding one to the overlapping area value between the suspected damaged patch area and the building core area and taking the square root; the second spatial constraint term is obtained by adding the natural constant e to the overlapping area value between the suspected damaged patch area and the building boundary buffer zone and taking the natural logarithm; the second spatial constraint term is obtained by adding one to the building tilt angle value and taking the square root; the third crack width value is obtained by adding the natural constant e to the structural crack width value and taking the natural logarithm; the first time difference value, the second time difference value, and the third time difference value are added one to the sum of the building settlement value and the building horizontal displacement value, divided by the total building height value, and then added to the constant one; the spatiotemporal suppression term is obtained by multiplying the first spatial constraint term, the second spatial constraint term, the tilt angle constraint term, and the crack constraint term and dividing by the spatiotemporal suppression term.

[0014] Further, the specific steps for generating object-level seismic correlation data are as follows: Compare the alternating confirmation segments under the same building survey object number value. If the object alternating confirmation constraint value corresponding to the subsequent alternating confirmation segment is less than the object alternating confirmation constraint value corresponding to the preceding alternating confirmation segment, and at least one of the building tilt angle, structural crack width, building settlement, and building horizontal displacement values ​​corresponding to the subsequent alternating confirmation segment is greater than the corresponding values ​​of the preceding alternating confirmation segment, then the subsequent alternating confirmation segment is marked as a reverse rollback record. For alternating confirmation segments not marked as reverse rollback records, continue to be processed according to relative seismic data. The time sequence is compared. If the alternating confirmation constraint value of the object corresponding to the subsequent alternating confirmation segment is not less than the alternating confirmation constraint value of the object corresponding to the preceding alternating confirmation segment, and the time difference between the start relative post-earthquake time of the subsequent alternating confirmation segment and the end relative post-earthquake time of the preceding alternating confirmation segment is not greater than the continuous threshold, the subsequent alternating confirmation segment and the preceding alternating confirmation segment are compressed into a single object association record. When there are multiple single object association records under the same building survey object number value, the single object association record with the largest object alternating confirmation constraint value is retained to construct object-level seismic prevention association data.

[0015] Furthermore, based on object-level seismic defense correlation data, the specific steps for determining the locking version core and calculating the locking expansion index value when there is a boundary contact relationship or a suspected damage patch area overlap relationship between adjacent building objects are as follows: Read the object-level seismic defense correlation data and determine the single object association record corresponding to the number value of each building survey object as the locking version core of the current building object; determine whether there is a common boundary line segment between the boundary areas of adjacent building objects, or whether the overlapping area between the corresponding suspected damage patch areas is greater than the overlap threshold; for adjacent building objects that meet the boundary contact relationship or suspected damage patch area overlap relationship, read the overlap area value between the suspected damage patch area and the building core area, the overlap area value between the suspected damage patch area and the building boundary buffer zone, the object alternation confirmation constraint value and the relative post-earthquake time corresponding to the current locking version core, as well as the object alternation confirmation constraint value and the relative post-earthquake time corresponding to the adjacent building objects, and calculate the locking expansion index value of the current locking version core for adjacent building objects.

[0016] Further, the specific steps for calculating the locking expansion index value of the current locking version core relative to adjacent building objects are as follows: Add one to the overlap area value between the suspected damaged patch area and the building boundary buffer zone corresponding to the current locking version core, then divide by the sum of the overlap area values ​​between the suspected damaged patch area and the building boundary buffer zone, the overlap area value between the suspected damaged patch area and the building core area, and a constant one, to obtain the buffer expansion ratio term; add one to the absolute value of the difference between the alternating confirmation constraint values ​​of the objects corresponding to the current locking version core and adjacent building objects, and multiply by one the absolute value of the difference between the relative post-earthquake times corresponding to the current locking version core and adjacent building objects, to obtain the expansion suppression term; divide the buffer expansion ratio term by the expansion suppression term to obtain the locking expansion index value.

[0017] Further, the specific steps for generating cross-scale version locking data are as follows: When the locking expansion index value is greater than the locking expansion threshold, the current building object and adjacent building objects are merged into a block locking unit; the adjacent block locking units are further compared, and when there is a continuous connection chain formed by the sequential connection of the boundary contact relationship of two or more building objects between adjacent block locking units, and the relative post-earthquake time corresponding to the locking version core in the later block locking unit is later than the relative post-earthquake time corresponding to the locking version core in the previous block locking unit, the adjacent block locking units are merged into a district locking unit, and the locking version core with the largest object alternation confirmation constraint value in the district locking unit is determined as the district dominant locking core; the locking version cores, block locking units, district locking units and district dominant locking cores corresponding to the building survey object number values ​​are summarized to construct cross-scale version locking data.

[0018] Furthermore, based on cross-scale version-locked data, different levels of representational base elements are constructed. Representation replacement, fission mapping processing, and time-series locking display are performed according to the maximum value of the object's alternating confirmation constraint value and the relative post-earthquake time. The specific steps for generating an adaptive multi-scale visualization map of urban earthquake prevention are as follows: Read the cross-scale version-locked data, and use the boundary areas of building objects, the boundaries of street-level locking units, and the boundaries of district-level locking units as building-level, street-level, and district-level representational base elements, respectively; express the building-level, street-level, and district-level representational base elements respectively. Replacement occurs when the maximum value of the alternating confirmation constraint value of the object corresponding to the lower-level representation element is greater than the maximum value of the alternating confirmation constraint value of the object corresponding to the upper-level representation element, and the relative post-earthquake time corresponding to the lower-level representation element is later than the relative post-earthquake time corresponding to the upper-level representation element; when the maximum value of the alternating confirmation constraint value of the object corresponding to the upper-level representation element is greater than or equal to the maximum value of the alternating confirmation constraint value of the object corresponding to the lower-level representation element, the upper-level representation element remains unchanged; fission mapping processing is performed on each level of representation element after the representation replacement is completed. When two or more objects within the same block-level representation unit have alternating confirmation constraint values ​​greater than the fission threshold, and the corresponding relative post-earthquake time difference is greater than the temporal difference threshold, the block-level representation unit is split into multiple building-level representation units. When two or more objects within the same area-level representation unit have alternating confirmation constraint values ​​greater than the fission threshold, and the corresponding relative post-earthquake time difference is greater than the temporal difference threshold, the area-level representation unit is split into multiple block-level representation units. Temporal locking display is applied to each level of representation unit after fission mapping processing. When the requested time of the map display is not less than the relative post-earthquake time corresponding to the locked version core or the dominant locking core of the area, the suspected damaged patch area and the boundary area of ​​the building object corresponding to the current expression unit are displayed; when the requested time of the map display is less than the relative post-earthquake time corresponding to the locked version core or the dominant locking core of the area, the historical expression unit with the largest object alternation confirmation constraint value in the corresponding level is displayed; the building-level expression units, block-level expression units and area-level expression units that have completed expression replacement, fission mapping processing and time-series locking display are overlaid to generate an adaptive multi-scale visualization map of urban earthquake prevention.

[0019] The second aspect of this invention provides an adaptive multi-scale visualization mapping system for urban earthquake prevention, characterized by comprising: a data acquisition and processing module, an object verification and association module, a span version locking module, and a constraint fission mapping module, wherein: the data acquisition and processing module is used to construct a unified spatial reference coordinate system, collect spatial and temporal data corresponding to the locations of building survey objects, remote sensing patch locations, IoT sensing device installation locations, and manually reported locations, construct an earthquake prevention multi-source basic dataset, and preprocess the earthquake prevention multi-source basic data; the object verification and association module is used to construct building object boundary regions, building core regions, building boundary buffer regions, and suspected damaged patch regions based on the preprocessed earthquake prevention multi-source basic data, combining overlapping area values ​​and position... The system determines the building survey object number corresponding to the current record by setting the falling relationship, and constructs an alternating confirmation segment according to the relative post-earthquake time, calculates the alternating confirmation constraint value of the object, and generates object-level seismic defense association data; the span version locking module is used to determine the locking version kernel based on the object-level seismic defense association data, and calculates the locking expansion index value when there is a boundary contact relationship or a suspected damage patch area overlap relationship between adjacent building objects, generating cross-scale version locking data; the constraint fission mapping module is used to construct different levels of expression base elements based on the cross-scale version locking data, and performs expression replacement, fission mapping processing and temporal locking display according to the maximum value of the object alternating confirmation constraint value and the relative post-earthquake time, generating an adaptive multi-scale visualization map of urban seismic defense.

[0020] Beneficial effects

[0021] The present invention has the following beneficial effects:

[0022] (1) An adaptive multi-scale visualization mapping method and system for urban earthquake prevention maps, which unifies the mapping and association of building survey objects, suspected damaged patches, IoT sensing device installation locations and manually reported locations, so that multi-source earthquake prevention data can form a stable association relationship for the same building object, thereby improving the consistency of building object identification and attribution.

[0023] (2) An adaptive multi-scale visualization mapping method and system for urban earthquake prevention, which constructs alternating confirmation segments according to the relative post-earthquake time and calculates the alternating confirmation constraint values ​​of objects in conjunction with multiple spatial and temporal constraints, transforms the post-earthquake state confirmation from a single judgment to a continuous quantitative judgment, thereby improving the stability of post-earthquake state identification.

[0024] (3) The method and system for adaptive multi-scale visualization mapping of urban earthquake prevention maps, by determining the locking version core and combining the locking expansion index value to generate block locking units, area locking units and area dominant locking cores, enables the cross-scale locking process to adaptively expand with the evolution of object state, thereby improving the reliability of cross-scale state transmission.

[0025] (4) An adaptive multi-scale visualization mapping method and system for urban earthquake prevention maps. By constructing different levels of representation elements and performing representation replacement, fission mapping processing and time-series locking display, the multi-scale map can be dynamically adjusted according to the object status and time sequence changes, thereby improving the visualization effect of urban earthquake prevention maps. Attached Figure Description

[0026] Figure 1 Flowchart of an adaptive multi-scale visualization mapping method for urban earthquake prevention;

[0027] Figure 2 A structural diagram of an adaptive multi-scale visualization mapping system for urban earthquake prevention;

[0028] Figure 3 Alternating confirmation constraint evolution decision diagram for objects;

[0029] Figure 4 A flowchart for constraint-driven adaptive multi-scale visualization mapping. Detailed Implementation

[0030] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0031] Please see Figures 1-4This invention provides a technical solution: an adaptive multi-scale visualization mapping method for urban earthquake prevention, comprising: S1, constructing a unified spatial reference coordinate system, collecting spatial and temporal data corresponding to the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations, constructing a multi-source earthquake prevention basic dataset, and preprocessing the multi-source earthquake prevention basic data; S2, constructing building object boundary regions, building kernel regions, building boundary buffer regions, and suspected damaged patch regions based on the preprocessed multi-source earthquake prevention basic data, and determining the building survey object corresponding to the current record by combining the overlapping area value and the location falling relationship. The system generates object-level seismic defense correlation data by assigning numbers to objects and constructing alternating confirmation segments based on relative post-earthquake times. S3, based on the object-level seismic defense correlation data, it determines the locking version kernel and calculates the locking expansion index value when there is boundary contact relationship or overlapping relationship between suspected damaged patches in the boundary areas of adjacent building objects, generating cross-scale version locking data. S4, based on the cross-scale version locking data, it constructs different levels of expression base elements and performs expression replacement, fission mapping processing, and temporal locking display according to the maximum value of the object alternating confirmation constraint value and the relative post-earthquake time, generating an adaptive multi-scale visualization map of urban seismic defense.

[0032] Specifically, a unified spatial reference coordinate system is constructed, and spatial and temporal data corresponding to the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations are collected to build a multi-source basic dataset for earthquake prevention. The specific steps for preprocessing the multi-source basic data for earthquake prevention are as follows: Select the plane coordinate reference datum and elevation reference datum corresponding to the urban earthquake prevention area. Use the plane coordinate reference datum to determine the east and north coordinate axes, and use the elevation reference datum to determine the height reference axis, thus constructing a unified spatial reference coordinate system. When constructing the unified spatial reference coordinate system, first read the mapping datum parameters, projection parameters, coordinate origin parameters, scale parameters, and elevation zero point parameters corresponding to the urban earthquake prevention area. The coordinates of the reference points in the planar coordinate reference datum are written into the initial positioning tables of the east and north coordinate axes, and the zero elevation point in the elevation reference datum is written into the initial positioning table of the height reference axis. This ensures that all subsequent data entering the same technical link have a unified spatial reference starting point. The technical principle is to first fix the spatial reference datum and then perform multi-source location mapping. This avoids spatial offsets caused by inconsistent projection datums for the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations under different acquisition sources. Then, the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations are uniformly mapped to a unified spatial reference coordinate system to collect building foundation data. The coordinates of the outline boundary points (eastward), the building base outline boundary points (northward), the suspected damaged patch boundary points (eastward), the suspected damaged patch boundary points (northward), the IoT sensing device installation points (eastward), the IoT sensing device installation points (northward), the manually reported locations (eastward), the manually reported locations (northward), and the building survey object number are all mapped to a unified spatial reference coordinate system. First, the building outline nodes in the building survey object locations are converted to the eastward and northward coordinates of the building base outline boundary points. Then, the patch boundary nodes in the remote sensing patch locations are converted to the eastward and northward coordinates of the suspected damaged patch boundary points. The coordinate values ​​are then converted into the installation point positioning results of the IoT sensing device installation location and the IoT sensing device installation point east coordinate value and IoT sensing device installation point north coordinate value. The reporting point positioning results of the manually reported location are then converted into the manually reported location east coordinate value and manually reported location north coordinate value. At the same time, the unique identifier in the building census ledger is written into the building census object number value, so that each set of spatial location data can be pointed back to a unique building object. The technical principle is that the outline-type location data is kept in the form of a boundary point sequence and the point-type location data is kept in the form of a single point coordinate. This can provide direct input for the subsequent construction of building object boundary areas, construction of suspected damaged patch areas, and location fall-in determination.Simultaneously, the system collects earthquake event occurrence timestamps, remote sensing acquisition timestamps, IoT sensing sampling timestamps, manually reported timestamps, total building height, total building area, building tilt angle, structural crack width, building settlement, and building horizontal displacement values ​​to construct a multi-source basic dataset for earthquake prevention. When collecting time-based data, earthquake event occurrence timestamps are written to an event start timetable, remote sensing acquisition timestamps are written to a patch acquisition timetable, IoT sensing sampling timestamps are written to a sensing sampling timetable, and manually reported timestamps are written to a reporting timetable. When collecting state-based data, total building height is directly entered from the building height record in the building survey results into the building attribute table, and total building area is directly entered from the building survey results into the building attribute table. The area record in the data is directly written into the building attribute table. The building tilt angle value is written into the attitude monitoring table based on the tilt angle sensor sampling result. The structural crack width value is written into the crack monitoring table based on the crack measurement result. The building settlement value is written into the settlement monitoring table based on the settlement monitoring result. The building horizontal displacement value is written into the displacement monitoring table based on the displacement monitoring result. This allows spatial data, temporal data, and structural status data to form a traceable data chain under the same building survey object number value. The technical principle is to use the earthquake event occurrence time value as a unified temporal reference and the building survey object number value as a unified object reference, so that the spatial, temporal, and status data have been bound at the object level before entering the subsequent judgment. For the collected seismic multi-source foundations... The data is processed using a precise time protocol master-slave clock offset compensation algorithm to unify the time of multi-source acquisition. During execution, the original timestamps corresponding to the earthquake event occurrence time, remote sensing acquisition time, IoT sensing sampling time, and manual reporting time are first extracted. Then, the precise time protocol master-slave clock offset compensation algorithm is used to uniformly synchronize the original timestamps, outputting the earthquake event occurrence time, remote sensing acquisition time, IoT sensing sampling time, and manual reporting time under a unified time reference. The technical principle is to eliminate the time offset caused by local clock drift of different acquisition devices, ensuring a unified starting point for subsequent calculations of relative post-earthquake times. A seven-parameter coordinate transformation algorithm is used to unify the spatial location reference. During execution, the data is first read... The system takes the eastward and northward coordinates of the building base outline boundary points, the eastward and northward coordinates of the suspected damaged patch boundary points, the eastward and northward coordinates of the IoT sensing device installation points, the eastward and northward coordinates of manually reported locations, and then performs coordinate projection conversion, offset correction, and axial normalization based on the transformation parameters corresponding to the unified spatial reference coordinate system. This outputs spatial location data under a unified reference. The technical principle is to project planar locations from different sources onto the same reference frame, eliminating spatial misalignment caused by differences in coordinate references. An abnormal data identification and processing algorithm is used to detect anomalies in the median absolute deviation.During execution, the median, absolute deviation, and deviation threshold are calculated for the total building height, total building area, building tilt angle, structural crack width, building settlement, and building horizontal displacement. Data items exceeding the deviation threshold are marked as anomalies. The technical principle lies in utilizing the insensitivity of the median absolute deviation anomaly detection algorithm to extreme outliers to identify sudden anomalies in multi-source monitoring data. The deviation threshold is obtained by multiplying the median absolute deviation of the total building height, total building area, building tilt angle, structural crack width, building settlement, and building horizontal displacement within the current batch of data by a preset amplification factor, which ranges from 2 to 5. A sliding median filtering algorithm is used to smooth local abrupt earthquake records. During execution, the building survey object number is used as the grouping key, and the sampling time values ​​from IoT sensing and manual reporting are used as the grouping keys. The corresponding sequence is a sliding order, replacing the building tilt angle, structural crack width, building settlement, and building horizontal displacement values ​​with the median value within a window. The technical principle is to suppress the amplified impact of local abrupt noise on the calculation of subsequent object alternating confirmation constraint values. A range normalization algorithm is used to perform numerical scale unification processing, outputting preprocessed seismic and anti-seismic multi-source foundation data. During execution, the maximum and minimum values ​​of the total building height, total building area, building tilt angle, structural crack width, building settlement, and building horizontal displacement values ​​within the current processing batch are extracted, and then dimensionless processing is performed based on the relative position of each data point within the current processing batch. The technical principle is to eliminate the dimensional influence caused by differences in the original value ranges of different physical quantities through the range normalization algorithm, ensuring that the calculation of subsequent object alternating confirmation constraint values ​​is based on preprocessed dimensionless data with a unified scale.

[0033] In this implementation plan, a unified spatial reference coordinate system is constructed, and unified processing of multi-source acquisition time, spatial location benchmark, abnormal data identification, smoothing of local abrupt seismic records, and unified numerical scale are performed on the multi-source basic data of seismic prevention. This ensures that the east-facing coordinates of the building base outline boundary points, the north-facing coordinates of the building base outline boundary points, the east-facing coordinates of the suspected damaged patch boundary points, the north-facing coordinates of the suspected damaged patch boundary points, the east-facing coordinates of the IoT sensing device installation points, the north-facing coordinates of the IoT sensing device installation points, the east-facing coordinates of the manually reported locations, and the north-facing coordinates of the manually reported locations are all consistent. Before being incorporated into subsequent calculations, the data chain of standard values, earthquake event occurrence times, remote sensing acquisition times, IoT sensing sampling times, manual reporting times, total building height, total building area, building tilt angle, structural crack width, building settlement, and building horizontal displacement values ​​is formed, establishing a unified object benchmark, unified time benchmark, unified spatial benchmark, and unified numerical benchmark. This improves the comparability, continuity, and computational stability of multi-source earthquake prevention basic data in subsequent processes such as building object boundary area construction, suspected damaged patch area construction, alternating confirmation segment construction, and object alternating confirmation constraint value calculation.

[0034] Specifically, the steps for constructing building object boundary regions, building core regions, building boundary buffer zones, and suspected damaged patch regions based on preprocessed seismic and fire prevention multi-source basic data, and determining the building survey object number value corresponding to the current record by combining the overlapping area value and the positional relationship, are as follows: Read the preprocessed seismic and fire prevention multi-source basic dataset; construct building object boundary regions according to building survey object number values; determine the inward offset distance according to the product of the total building height value and the preset offset ratio; then perform inward offset on the building object boundary regions based on the inward offset distance to generate the building core region; when constructing building object boundary regions according to building survey object number values, first extract the same... The east and north coordinates of all building base outline boundary points corresponding to the building survey object's ID value are then used to perform a closed connection on each boundary point according to the boundary point collection order, forming a closed outline area that corresponds one-to-one with the building survey object's ID value. This ensures that subsequent spatial inclusion determination and overlapping area calculation are based on the same building object's geometric boundary. When determining the inward offset distance by multiplying the total building height value by a preset offset ratio, and then performing inward offset on the building object's boundary area based on this inward offset distance, the total building height value corresponding to the building survey object's ID value is read first, and then the total building height value is multiplied by the preset offset ratio. Multiplying this by the inward offset distance of the current building object's boundary region, synchronous contraction is performed on each side of the building object's boundary region to obtain the building's core area. This core area corresponds to a stable region closer to the main body's center within the building object's boundary region. The technical principle lies in using the total building height value and a preset offset ratio to construct an inwardly contracted boundary that matches the building's scale. This allows the overlapping area values ​​of subsequent suspected damaged patches and the building's core area to more effectively represent the degree of damage to the building's core. The preset offset ratio is pre-calibrated based on the correspondence between the typical planar dimensions of building survey objects within the urban earthquake-resistant area and the total building height value, and the preset offset ratio is set to 0. The range is from 0.05 to 0.20. The area between the boundary region of the building object and the core region of the building is defined as the building boundary buffer. When defining the building boundary buffer, the closed loop area between the outer boundary of the building object boundary region and the outer boundary of the core region of the building is extracted as the building boundary buffer. The building boundary buffer is used to characterize the extent of damage expansion near the edge of the building. The technical principle is that by dividing the boundary region of the building object into a two-level spatial structure of the core region and the building boundary buffer, the damage to the main core and the damage to the boundary expansion can be expressed separately in subsequent calculations, thereby improving the ability of the same suspected damaged patch area to distinguish different damage levels inside the building object. The eastward coordinate values ​​and northward coordinate values ​​of the suspected damaged patch boundary points are collected according to the same remote sensing acquisition time value, and then connected end to end in the order of the boundary points to construct the suspected damaged patch area.When aggregating the eastward and northward coordinates of suspected damaged patch boundary points according to the same remote sensing acquisition time value, the remote sensing acquisition time value is first used as the aggregation key. All eastward and northward coordinates of suspected damaged patch boundary points under the same remote sensing acquisition time value are written into the same patch boundary sequence. Then, the boundary points are sequentially closed to form suspected damaged patch regions, ensuring that each suspected damaged patch region corresponds to a unique remote sensing acquisition time value. The technical principle is to maintain the synchronous correspondence between the patch spatial outline and the acquisition time using the remote sensing acquisition time value, which can avoid the mixing of patch boundaries at different times, leading to subsequent errors in the attribution of building objects. The overlap area between each suspected damaged patch region and the core area of ​​each building, as well as the overlap area with the buffer zone of each building boundary, are calculated separately. When calculating the overlap area, each suspected damaged patch region is first compared with all building core areas... The core area is spatially overlaid, and then each suspected damaged patch area is spatially overlaid with all building boundary buffer zones. The area of ​​the overlapping area formed after overlay is written into the overlap area value between the suspected damaged patch area and the building core area, and the overlap area value between the suspected damaged patch area and the building boundary buffer zone. This allows the same suspected damaged patch area to form multiple sets of comparable spatial fit amounts under different building objects. The technical principle is that by calculating the overlap area of ​​the core area and the overlap area of ​​the boundary area respectively, the difference in the effect of the suspected damaged patch area on the central area of ​​the building body and the edge area of ​​the building can be simultaneously characterized, providing dual spatial basis for the subsequent determination of the building survey object number value; and the building survey object number value that meets the condition of having the largest patch overlap area value and whose IoT sensing device installation point or manual reporting location falls in the building core area or building boundary buffer zone is determined as the building survey object number value corresponding to the current record.When determining the building survey object number corresponding to the current record, the overlapping area values ​​of all suspected damaged patch areas and the building core area, as well as the overlapping area values ​​of the suspected damaged patch areas and the building boundary buffer zone, corresponding to the same suspected damaged patch area are first compared. The candidate building survey object number value with the largest patch overlapping area value is extracted. Then, the east coordinate value, north coordinate value, east coordinate value, and north coordinate value of the IoT sensing device installation point, are read. The following steps are then performed to determine whether the IoT sensing device installation point falls within the building core area, whether the IoT sensing device installation point falls within the building boundary buffer zone, and whether the manually reported location falls within the building core area. The determination of whether a manually reported location falls within the building boundary buffer zone is based on the following principle: when the candidate building survey object number value corresponds to a building core area or a building boundary buffer zone that receives at least one type of location from the IoT sensing device installation point or the manually reported location, the candidate building survey object number value is determined as the building survey object number value corresponding to the current record. The technical principle is to first complete the main spatial matching using the patch overlap area value, and then complete the location verification using the IoT sensing device installation point and the manually reported location. This avoids misattribution of building survey object numbers based solely on the coverage relationship of the outer edge of suspected damaged patch areas, ensuring that the building survey object number value corresponding to the current record simultaneously satisfies both spatial fit constraints and location fall-in constraints.

[0035] In this implementation plan, a boundary region, a core region, and a buffer zone for a building object are constructed based on the building survey object number value. The overlapping area values ​​of the suspected damaged patch area and the core region, the overlapping area values ​​of the suspected damaged patch area and the buffer zone, the east coordinate value of the IoT sensing device installation point, the north coordinate value of the IoT sensing device installation point, the east coordinate value of the manually reported location, and the north coordinate value of the manually reported location are all incorporated into the process of determining the building survey object number value corresponding to the current record. This provides a unified association basis for the object-level spatial attribution expression in subsequent image synthesis, enabling the representation results of the suspected damaged patch area on the degree of damage to the core of the building and the degree of damage to the boundary extension to form a unified attribution constraint under the same building survey object number value. This improves the accuracy and stability of the determination result of the building survey object number value corresponding to the current record and the reliability of the calculation of the subsequent object alternation confirmation constraint value.

[0036] Specifically, the steps for constructing alternating confirmation segments based on relative post-earthquake times and calculating the alternating confirmation constraint values ​​for objects are as follows: Using the earthquake event occurrence time as a benchmark, the remote sensing acquisition time, IoT sensing sampling time, and manual reporting time are converted into relative post-earthquake times. When converting relative post-earthquake times, the earthquake event occurrence time is first written into the event start time table corresponding to the current building survey object number. Then, the earthquake event occurrence time is subtracted from the remote sensing acquisition time, IoT sensing sampling time, and manual reporting time, respectively, to obtain the relative post-earthquake times corresponding to the remote sensing records, IoT records, and manual reporting records. This ensures that records from different sources all fall under the same post-earthquake time benchmark. The technical principle is... By establishing a unified zero point using the moment of earthquake events, records from different sources and with different acquisition rhythms can be converted into comparable records on the same timeline. All records corresponding to the same building object are aggregated according to the building survey object number, and then arranged in ascending order based on the relative post-earthquake time. When aggregating records by building survey object number, all remote sensing records, all IoT records, and all manually reported records under the current building survey object number are first extracted, and then sorted in ascending order using the relative post-earthquake time as the sorting key. This ensures that records under the same building survey object number form a continuous time series. The technical principle is that object-level aggregation is completed first using the building survey object number, and then time-level sorting is completed using the relative post-earthquake time, which can guarantee alternating confirmation segments. The construction process always revolves around the same building object. When adjacent records correspond to the same building survey object number and contain three types of sources: remote sensing, IoT, and manual reporting, the adjacent records are merged into an alternating confirmation segment. The first time difference between the remote sensing acquisition time value and the earthquake event occurrence time value, the second time difference between the IoT sensing sampling time value and the earthquake event occurrence time value, and the third time difference between the manual reporting time value and the earthquake event occurrence time value are calculated separately. When merging adjacent records into an alternating confirmation segment, the source composition of adjacent records under the current building survey object number value is first retrieved, and then it is confirmed that the adjacent record set simultaneously contains records corresponding to remote sensing acquisition time values, IoT sensing sampling time values, and manual reporting time values. The time value corresponds to a record, and then the set of adjacent records is written into the alternating confirmation segment sequence list. Then, the remote sensing acquisition time value minus the earthquake event occurrence time value is written into the first time difference value, the IoT sensing sampling time value minus the earthquake event occurrence time value is written into the second time difference value, and the manual reporting time value minus the earthquake event occurrence time value is written into the third time difference value. The technical principle is that by forcing the alternating confirmation segment to simultaneously include records from three sources, the subsequent object alternating confirmation constraint value can be established on the basis of multi-source synchronous confirmation, avoiding confirmation deviation caused by a single source record. For each alternating confirmation segment, based on the preprocessed dimensionless data, the overlap area value between the suspected damaged patch area and the building core area is added by one and the square root is taken to obtain the first spatial constraint term.When calculating the first spatial constraint term, the overlap area value between the suspected damaged patch area and the building core area corresponding to the current alternating confirmation segment is written into the spatial core fitting table, then incremented by one, and then a square root transformation is performed. This ensures that records with larger overlap area values ​​maintain an increasing trend, while records with smaller overlap area values ​​are not over-compressed due to their low original magnitude. The technical principle is that when the overlap area value between the suspected damaged patch area and the building core area is small, it is more sensitive to distinguishing the degree of damage to the main core. When the value is large, continued linear growth can easily lead to excessive amplification. Therefore, the square root transformation compresses the increase in high-value intervals, making the change in overlap area value maintain a smoother correspondence with the change in the perceived intensity of damage to the main core. The overlap area value between the suspected damaged patch area and the building core area is then added to the spatial core fitting table. The second spatial constraint term is obtained by adding the natural constant e to the overlapping area value of the boundary buffer and then taking the natural logarithm. When calculating the second spatial constraint term, the overlapping area value between the suspected damaged patch area corresponding to the current alternating confirmation segment and the building boundary buffer is written into the boundary expansion fitting table, added to the natural constant e, and then a natural logarithmic transformation is performed. This ensures that boundary expansion damage contributes incrementally to the calculation result without excessive amplification. The technical principle is that boundary expansion damage is more critical for identifying the outward expansion trend in the initial growth stage, while continued linear amplification after the value continues to increase weakens the effect of other constraint terms. Therefore, the natural logarithmic transformation preserves the trend information of the increasing boundary expansion area while suppressing the unilateral dominance of large-scale boundary expansion values ​​on the final result. The building tilt... The tilt constraint term is obtained by adding one to the angle value and taking the square root. When calculating the tilt constraint term, the building tilt value is written into the attitude offset table, then the addition process is performed, followed by a square root transformation. This ensures that the tilt offset continuously strengthens the alternating confirmation constraint value of the object. The technical principle is that changes in the building tilt value within the low-value range are usually sufficient to characterize the building attitude offset trend. Continued linear growth in the high-value range can easily cause individual abnormal tilt records to exert an excessive influence on the overall judgment. Therefore, the square root transformation compresses the increase in high values, suppressing the abnormal amplification of results by outlier tilt values ​​and making the building tilt value closer to the perceived damage level in manual interpretation. That is, low to medium tilt changes can clearly differentiate, while high tilt changes remain enhanced but not... Imbalance; the crack constraint term is obtained by adding the natural constant e to the crack width value and taking the natural logarithm; when calculating the crack constraint term, the crack width value is written into the crack characterization table, added to the natural constant e, and then the natural logarithmic transformation is performed so that the crack width change can participate in the confirmation calculation in an incremental manner. The technical principle is that the increase of the crack width value in a small range often corresponds to the critical stage of the damage level changing from slight to obvious. The human perception of crack width change is also closer to the nonlinear pattern of sensitivity in the early stage and slowing down in the later stage. Therefore, the natural logarithmic transformation can enhance the ability to distinguish small to medium crack width changes, while suppressing the unilateral dominance of the maximum crack width value on the joint calculation, so that the mapping relationship between crack width value and damage level is more coordinated.To obtain the spatiotemporal suppression term, first increment the first, second, and third time differences by one, then multiply them sequentially. Next, divide the sum of the building settlement and horizontal displacement values ​​by the total building height, increment by one, and add this sum to a constant. Alternatively, to calculate the spatiotemporal suppression term, first increment the first, second, and third time differences by one, then multiply them sequentially to form the time suppression base term. Then, read the building settlement, horizontal displacement, and total building height values, sum the building settlement and horizontal displacement values, divide by the total building height, and increment by one to form the deformation suppression base term. Subsequently, the deformation suppression base term is added to the constant and written into the spatiotemporal suppression term. The technical principle is to use three types of time difference values ​​to characterize the degree of deviation of the multi-source confirmation time series, and to use the ratio of building settlement value and building horizontal displacement value to the total building height value to characterize the degree of deformation disturbance. Then, the time disturbance and deformation disturbance are written into the denominator side suppression structure, so that records with a large degree of time series dispersion and strong deformation disturbance are suppressed in the final judgment. After multiplying the first spatial constraint term, the second spatial constraint term, the tilt angle constraint term and the crack constraint term, and dividing by the spatiotemporal suppression term, the corresponding object alternating confirmation constraint value is obtained. When obtaining the alternating confirmation constraint value of an object, the first spatial constraint item, the second spatial constraint item, the tilt angle constraint item, and the crack constraint item are first written into the spatial enhancement joint table. Then, a product calculation is performed to obtain the spatial enhancement base value. Finally, the spatial enhancement base value is divided by the spatiotemporal suppression item to output the alternating confirmation constraint value of the object corresponding to the current alternating confirmation segment. The technical principle is that the degree of core damage fit, the degree of boundary expansion fit, the degree of building attitude deviation, and the degree of crack development are all written into the enhancement side, and the degree of multi-source time deviation and the degree of overall building deformation disturbance are all written into the suppression side. This allows the alternating confirmation constraint value of the object to simultaneously characterize the damage confirmation intensity, temporal consistency, and deformation stability of the current alternating confirmation segment.

[0037] The specific formula for calculating the object alternating confirmation constraint value is as follows:

[0038] ;

[0039] In the formula, This indicates that the constraint values ​​are alternately confirmed by the object. This indicates the area of ​​overlap between the suspected damaged patch and the building's core area. This represents the area of ​​overlap between the suspected damaged patch and the building boundary buffer zone. Indicates the building's tilt angle value. This represents the width of the structural crack. Indicates the first time difference. Indicates the second time difference. Indicates the third time difference. Indicates the building settlement value. This represents the horizontal displacement value of the building. This indicates the total height of the building.

[0040] In this implementation plan, the timing values ​​obtained by remote sensing, IoT-based sampling, and manually reported are uniformly constrained by the timing value of the earthquake event. The overlapping area values ​​of suspected damaged patch areas and building core areas, the overlapping area values ​​of suspected damaged patch areas and building boundary buffer zones, building tilt angle, structural crack width, first time difference, second time difference, third time difference, building settlement, building horizontal displacement, and total building height are all incorporated into the calculation process of the object alternating confirmation constraint value. This allows records from different sources under the same building survey object number to form continuous confirmation results under a unified time reference and a unified damage judgment scale. This improves the stability of the alternating confirmation segment identification results, the ability of the object alternating confirmation constraint value to represent post-earthquake state changes, and the reliability of subsequent reverse backtracking record identification and single object association record retention.

[0041] Specifically, the steps for generating object-level seismic correlation data are as follows: Compare the alternating confirmation segments under the same building survey object number value. If the object alternating confirmation constraint value corresponding to the subsequent alternating confirmation segment is less than that corresponding to the preceding alternating confirmation segment, and at least one of the building tilt angle, structural crack width, building settlement, and building horizontal displacement values ​​corresponding to the subsequent alternating confirmation segment is greater than the corresponding values ​​of the preceding alternating confirmation segment, then mark the subsequent alternating confirmation segment as a reverse rollback record. When comparing the alternating confirmation segments under the same building survey object number value, first read the current building survey data in ascending order of relative post-earthquake time. Search all alternating confirmation segments under the object number value, then write the previous alternating confirmation segment into the preceding comparison bit and the next alternating confirmation segment into the following comparison bit. Next, extract the object alternating confirmation constraint value, building tilt angle value, structural crack width value, building settlement value, and building horizontal displacement value corresponding to the preceding alternating confirmation segment. Then, extract the object alternating confirmation constraint value, building tilt angle value, structural crack width value, building settlement value, and building horizontal displacement value corresponding to the following alternating confirmation segment. Perform item-by-item comparison on the two sets of data. When marking the following alternating confirmation segment as a reverse rollback record, first determine whether the object alternating confirmation constraint value corresponding to the following alternating confirmation segment is less than the preceding alternating confirmation segment. The system iterates through the alternating confirmation constraint values ​​of objects corresponding to the previous confirmation segment, then determines whether at least one of the following values—building tilt angle, structural crack width, building settlement, and building horizontal displacement—is greater than the corresponding value of the preceding alternating confirmation segment. If both determinations are met simultaneously, the subsequent alternating confirmation segment is written into the reverse rollback record table. The technical principle is that if the alternating confirmation constraint value decreases between the preceding and subsequent alternating confirmation segments, it indicates a weakening of the multi-source joint confirmation strength. If, at this time, there are still increasing values ​​among the building tilt angle, structural crack width, building settlement, and building horizontal displacement, it indicates that the building damage continues to manifest. However, the joint confirmation results weakened, and such records are not suitable to be retained as stable association results. Therefore, they need to be removed from the subsequent compressed links by marking them as reverse backtracking records. For the alternating confirmation segments that are not marked as reverse backtracking records, they are compared in order of relative post-earthquake time. When the object alternating confirmation constraint value corresponding to the subsequent alternating confirmation segment is not less than the object alternating confirmation constraint value corresponding to the preceding alternating confirmation segment, and the time difference between the starting relative post-earthquake time of the subsequent alternating confirmation segment and the ending relative post-earthquake time of the preceding alternating confirmation segment is not greater than the continuous threshold, the subsequent alternating confirmation segment and the preceding alternating confirmation segment are compressed into a single object association record.When comparing alternating confirmation segments not marked as reverse rollback records in order of relative post-earthquake time, the alternating confirmation constraint values ​​of the objects corresponding to the subsequent alternating confirmation segments are extracted first, followed by the extraction of the alternating confirmation constraint values ​​of the objects corresponding to the preceding alternating confirmation segments. The size of the alternating confirmation constraint values ​​is then judged, and the starting relative post-earthquake time of the subsequent alternating confirmation segment and the ending relative post-earthquake time of the preceding alternating confirmation segment are extracted. The time difference between the two is calculated, and this time difference is compared with a continuous threshold. The continuous threshold is pre-set based on the minimum sampling interval of remote sensing acquisition time values, IoT sensing sampling time values, and manually reported time values ​​within the current processing batch, and the continuous threshold value is 1 to 3 times the minimum sampling interval. When compressing the subsequent alternating confirmation segments and the preceding alternating confirmation segments into a single object-related record, the preceding alternating confirmation segments that satisfy the conditions of no decrease in object alternating confirmation constraint values ​​and time differences not exceeding the continuous threshold are first compared with the subsequent alternating confirmation segments. The confirmed segments are written into the same compressed merging table. The earliest starting relative post-earthquake time in the compressed merging table is used as the starting time marker for the single object association record, and the latest ending relative post-earthquake time in the compressed merging table is used as the ending time marker for the single object association record. The maximum object alternation confirmation constraint value in the compressed merging table is used as the confirmation strength marker for the single object association record. The technical principle is that if the object alternation confirmation constraint value does not decrease, it indicates that the multi-source joint confirmation strength remains stable during the time progression. If the time difference between the starting relative post-earthquake time and the ending relative post-earthquake time does not exceed the continuity threshold, it indicates that there is no obvious temporal break between adjacent alternation confirmation segments. Therefore, two alternation confirmation results can be compressed into the same continuous object association result. When there are multiple single object association records under the same building survey object number value, the single object association record with the largest object alternation confirmation constraint value is retained to construct object-level seismic prevention association data. When retaining the single object association record with the largest alternating confirmation constraint value, the alternating confirmation constraint values ​​corresponding to all single object association records under the same building survey object number are first extracted. Then, these values ​​are sorted from largest to smallest, and the single object association record with the highest sorting result is written into the object-level seismic prevention association data table. The technical principle is that, under the same building survey object number, the larger the alternating confirmation constraint value corresponding to a single object association record, the higher the multi-source joint confirmation strength, the stronger the spatiotemporal consistency, and the more stable the damage characterization matching degree. Therefore, using the single object association record with the largest alternating confirmation constraint value as the final retention result can maintain higher stability, uniqueness, and traceability of the object-level seismic prevention association data.

[0042] In this implementation plan, by performing continuous comparison, reverse rollback record marking, and single-object association record compression and retention on alternating confirmation segments under the same building survey object number value, the object alternating confirmation constraint value, building tilt angle value, structural crack width value, building settlement value, building horizontal displacement value, starting relative post-earthquake time, and ending relative post-earthquake time form a sequential screening constraint relationship in the same judgment link. This improves the accuracy of object-level seismic prevention association data in representing the evolution of continuous post-earthquake states and enhances the stability, uniqueness, and temporal continuity of the final retained records.

[0043] In this embodiment, Table 1 is a data table of alternating confirmation constraint values ​​for objects, listing exemplary values ​​used to illustrate the calculation process of alternating confirmation constraint values ​​under five alternating confirmation segments corresponding to the same building survey object number. Specifically: In alternating confirmation segment 1, the overlap area between the suspected damaged patch area and the building core area is 18.6, the overlap area between the suspected damaged patch area and the building boundary buffer zone is 7.4, the building tilt angle is 1.8, the structural crack width is 2.6, the first time difference is 2, the second time difference is 3, the third time difference is 4, the building settlement is 6.2, the building horizontal displacement is 5.4, the total building height is 28.0, and the object alternating confirmation constraint value is 0.1871; In alternating confirmation segment 2, the overlap area between the suspected damaged patch area and the building core area is 22.3, the suspected damaged... The overlap area between the damaged patch area and the building boundary buffer zone is 9.1, the building tilt angle is 2.1, the structural crack width is 3.1, the first time difference is 2, the second time difference is 2, the third time difference is 3, the building settlement is 7.0, the building horizontal displacement is 6.1, the total building height is 28.0, and the object alternating confirmation constraint is 0.3578; in alternating confirmation segment 3, the overlap area between the suspected damaged patch area and the building core area is 27.8, the overlap area between the suspected damaged patch area and the building boundary buffer zone is 11.6, the building tilt angle is 2.7, and the structural crack width is... 2.8, first time difference value is 1, second time difference value is 2, third time difference value is 2, building settlement value is 8.4, building horizontal displacement value is 7.3, total building height value is 28.0, object alternating confirmation constraint value is 0.9713; in alternating confirmation segment 4, the overlap area between the suspected damaged patch area and the building core area is 25.4, the overlap area between the suspected damaged patch area and the building boundary buffer zone is 10.2, building tilt angle is 2.4, structural crack width is 3.5, first time difference value is 3, second time difference value is 4, third time difference value is 4, building settlement value is 6.1. The building's horizontal displacement is 8.0, the total building height is 28.0, and the object alternation confirmation constraint is 0.1155. In alternation confirmation segment 5, the overlap area between the suspected damaged patch area and the building's core area is 31.2, the overlap area between the suspected damaged patch area and the building boundary buffer zone is 12.8, the building tilt angle is 3.0, the structural crack width is 4.2, the first time difference is 1, the second time difference is 1, the third time difference is 2, the building settlement is 10.3, the building's horizontal displacement is 8.7, the total building height is 28.0, and the object alternation confirmation constraint is 1.9202.As shown in Table 1, during the continuous alternating confirmation process corresponding to the same building survey object number value, there are significant differences in the overlapping area values ​​of the suspected damaged patch area and the building core area, the overlapping area values ​​of the suspected damaged patch area and the building boundary buffer zone, the building tilt angle, the structural crack width, the first time difference, the second time difference, the third time difference, the building settlement value, and the building horizontal displacement value corresponding to each alternating confirmation segment. This results in significant differences in the alternating confirmation constraint values ​​of the objects corresponding to each alternating confirmation segment.

[0044] Table 1. Object Alternate Confirmation Constraint Value Data Table

[0045] Alternating confirmation segment 1 18.6 7.4 1.8 2.6 2 3 4 6.2 5.4 28.0 0.1871 2 22.3 9.1 2.1 3.1 2 2 3 7.0 6.1 28.0 0.3578 3 27.8 11.6 2.7 2.8 1 2 2 8.4 7.3 28.0 0.9713 4 25.4 10.2 2.4 3.5 3 4 4 6.1 8.0 28.0 0.1155 5 31.2 12.8 3.0 4.2 1 1 2 10.3 8.7 28.0 1.9202

[0046] like Figure 3 As shown, this figure is a determination diagram of the evolution of alternating confirmation constraints for objects. The figure sequentially presents the evolution results of alternating confirmation constraints for alternating confirmation segments 1 to 5 according to the relative post-earthquake time sequence. Each alternating confirmation segment has a first spatial constraint term, a second spatial constraint term, a dip angle constraint term, and a crack constraint term on its upper side, corresponding to the spatial confirmation enhancement source of the numerator of the alternating confirmation constraint value. Each alternating confirmation segment also has a first time difference term, a second time difference term, a third time difference term, and a settlement displacement suppression term on its lower side, corresponding to the spatiotemporal suppression source of the denominator of the alternating confirmation constraint value. The diamond-shaped determination kernel in the middle of each alternating confirmation segment is used to characterize the alternating confirmation constraints of the corresponding segment. The value is used to represent the continuous comparison relationship between adjacent alternating confirmation segments through evolutionary bridges. Specifically, a backtracking breakpoint is set between alternating confirmation segment 3 and alternating confirmation segment 4, and a reverse backtracking is marked on alternating confirmation segment 4. This is used to indicate that the alternating confirmation constraint value of the object corresponding to alternating confirmation segment 4 is lower than that of the preceding alternating confirmation segment 3, and the structural crack width value and building horizontal displacement value corresponding to alternating confirmation segment 4 are greater than the corresponding values ​​of alternating confirmation segment 3, thus satisfying the identification conditions of the reverse backtracking record. The alternating confirmation constraint value of the object corresponding to alternating confirmation segment 5 is the highest and is marked as the final retention. This is used to indicate that after comparison and compression retention, this alternating confirmation segment is determined to be a single object associated record.

[0047] Specifically, the steps for determining the locking version core based on object-level seismic defense correlation data and calculating the locking expansion index value when there is boundary contact relationship or overlapping relationship of suspected damaged patch areas in the boundary areas of adjacent building objects are as follows: Read the object-level seismic defense correlation data and determine the single object association record corresponding to each building survey object number as the locking version core of the current building object; when reading the object-level seismic defense correlation data, first extract the building survey object number value, single object association record, building object boundary area, overlapping area value between suspected damaged patch area and building core area, overlapping area value between suspected damaged patch area and building boundary buffer area, object alternating confirmation constraint value, and relative post-earthquake time, then according to the building survey... The object number value is used to group and aggregate all single-object associated records. The single-object associated records retained under each building survey object number value are written into the corresponding object's lock candidate table. When determining the single-object associated record corresponding to each building survey object number value as the current building object's lock version core, the single-object associated records in the lock candidate table are directly written into the lock version core table, ensuring that each building survey object number value corresponds to a unique lock version core. The technical principle is that the object-level seismic fortification associated data has undergone alternating confirmation segment comparison, reverse rollback record removal, and single-object associated record compression and retention. The alternating confirmation constraint value of the object corresponding to the currently retained record can stably characterize the current confirmed intensity under that building survey object number value. Therefore, it can be directly used as the benchmark record for subsequent cross-object extension locking; determine whether there is a common boundary line segment between adjacent building object boundary areas, or whether the overlapping area between corresponding suspected damaged patch areas is greater than the overlap threshold; when determining whether there is a boundary contact relationship between adjacent building object boundary areas, first extract the outer contour of the building object boundary area corresponding to any two building survey object number values, and then perform boundary adjacency search on the two; when there is a common contact boundary between the two outer contours, write the pair of building survey object number values ​​into the boundary contact relationship table; when determining whether there is an overlapping relationship between suspected damaged patch areas, first extract the suspected damaged patch areas corresponding to two adjacent building survey object number values, and then perform boundary adjacency search on the two Spatial overlap retrieval is performed. When the spatial overlap area between two objects exceeds a preset overlap threshold, the corresponding building survey object IDs are written into the patch overlap relationship table. The technical principle is that boundary contact relationships are used to characterize the direct spatial adjacency between building objects, while suspected damaged patch area overlap relationships are used to characterize the cross-boundary continuity of post-earthquake damage between adjacent building objects. By introducing both types of relationships simultaneously, simple geometric adjacency and damage propagation continuity can be incorporated into the subsequent locking and expansion judgment link. The overlap threshold is preset based on the minimum effective identification area of ​​suspected damaged patch areas within the urban earthquake prevention area, and the overlap threshold value is 5% to 20% of the median area of ​​all suspected damaged patch areas in the current processing batch.For adjacent building objects that satisfy boundary contact relationships or suspected damaged patch area overlap relationships, the overlapping area values ​​of the suspected damaged patch area and the building core area, the overlapping area values ​​of the suspected damaged patch area and the building boundary buffer area corresponding to the current locking version kernel, the object alternation confirmation constraint value, and the relative post-earthquake time are read, as well as the object alternation confirmation constraint value and the relative post-earthquake time corresponding to the adjacent building objects. The locking expansion index value of the current locking version kernel for adjacent building objects is then calculated, thus providing a directly executable quantitative basis for the subsequent formation of block locking units.

[0048] In this implementation plan, by using the building survey object number as the main collection line, the associated records of a single object are stably mapped to the locking version core. Then, the overlapping area values ​​of the building object boundary area, the suspected damaged patch area and the building core area, the overlapping area value of the suspected damaged patch area and the building boundary buffer area, the object alternating confirmation constraint value, and the relative post-earthquake time are all included in the locking expansion index value calculation process. This makes the spatial adjacency relationship between adjacent building objects, the outward expansion trend of the damaged boundary, the degree of closeness of the confirmed intensity, and the degree of closeness of the temporal progression form a unified expansion constraint within the same judgment link. This improves the accuracy of the judgment when the locking version core expands to adjacent building objects and enhances the stability, continuity, and reliability of the subsequent block locking unit formation process.

[0049] Specifically, the steps for calculating the lock expansion index value of the current lock version core relative to adjacent building objects are as follows: Based on the preprocessed dimensionless data, add one to the overlap area value between the suspected damaged patch area and the building boundary buffer corresponding to the current lock version core, then divide by the sum of the overlap area value between the suspected damaged patch area and the building boundary buffer, the overlap area value between the suspected damaged patch area and the building core area, and a constant one to obtain the buffer expansion ratio term; when calculating the buffer expansion ratio term, first read the overlap area value between the suspected damaged patch area and the building boundary buffer corresponding to the current lock version core, then read the overlap area value between the suspected damaged patch area and the building core area corresponding to the current lock version core, and then add the overlap area value between the suspected damaged patch area and the building boundary buffer. The overlap area between the damaged patch area and the building boundary buffer is incremented by one and written into the numerator-side calculation table. The sum of the overlap area between the suspected damaged patch area and the building boundary buffer, the overlap area between the suspected damaged patch area and the building core area, and a constant is written into the denominator-side calculation table. The numerator-side result is then divided by the denominator-side result to output the buffer expansion ratio. The technical principle is that the overlap area between the suspected damaged patch area and the building boundary buffer, and the overlap area between the suspected damaged patch area and the building core area together represent the total amount of bonding within the currently locked version core. Among them, the overlap area between the suspected damaged patch area and the building boundary buffer represents the degree of boundary damage bonding, and the overlap area between the suspected damaged patch area and the building core area represents the degree of bonding within the core. The overlapping area value characterizes the degree of core damage adhesion. Adding both the overlapping area value and a constant to the denominator is equivalent to using the total adhesion within the current locked version kernel as a normalization benchmark, and then adding one to the overlapping area value between the suspected damaged patch area and the building boundary buffer zone as the numerator. This allows the buffer expansion ratio term to stably characterize the proportion of boundary adhesion in the total adhesion, thus avoiding the amplified bias caused by magnitude differences when simply using absolute area values. Adding one to the absolute value of the difference between the alternating confirmation constraint value of the object corresponding to the current locked version kernel and the alternating confirmation constraint value of the adjacent building object, and adding one to the absolute value of the difference between the relative post-earthquake time corresponding to the current locked version kernel and the relative post-earthquake time corresponding to the adjacent building object... The values ​​are multiplied by one to obtain the extended suppression term. When calculating the extended suppression term, first, the alternating confirmation constraint value of the object corresponding to the current locked version core is read, then the alternating confirmation constraint value of the adjacent building object is read, the difference between the alternating confirmation constraint values ​​is calculated, the absolute value of the difference is processed, and then one is added to form the confirmation difference suppression sub-term. Next, the relative post-earthquake time corresponding to the current locked version core is read, then the relative post-earthquake time corresponding to the adjacent building object is read, the difference between the relative post-earthquake times is calculated, the absolute value of the difference is processed, and then one is added to form the temporal difference suppression sub-term. Finally, the confirmation difference suppression sub-term and the temporal difference suppression sub-term are multiplied to calculate the extended suppression term.The technical principle is that the greater the difference between the alternating confirmation constraint value of the object corresponding to the current locked version core and the alternating confirmation constraint value of the object corresponding to the adjacent building object, the more inconsistent they are in terms of confirmation strength. The greater the difference between the relative post-earthquake time corresponding to the current locked version core and the relative post-earthquake time corresponding to the adjacent building object, the more asynchronous they are in terms of post-earthquake evolution. Therefore, by writing the confirmation difference suppression sub-item and the temporal difference suppression sub-item into the expansion suppression term, the boundary expansion capability of the current locked version core can be jointly restricted by the confirmation difference and the temporal difference. The lock expansion index value is obtained by dividing the buffer expansion ratio term by the expansion suppression term. When obtaining the lock expansion index value, the buffer expansion ratio is first written into the expansion enhancement table, and then the expansion inhibition ratio is written into the expansion constraint table. Subsequently, the buffer expansion ratio is divided by the expansion inhibition ratio to output the lock expansion index value of the current lock version core relative to adjacent building objects. The technical principle is that by writing the boundary expansion trend of the current lock version core into the enhancement side, and then writing the confirmation difference and timing difference between the current lock version core and adjacent building objects into the inhibition side, the lock expansion index value can simultaneously represent expansion intention, expansion resistance, and cross-object locking feasibility in the same calculation chain, thereby providing a directly executable quantitative judgment basis for the subsequent generation of block lock units.

[0050] The specific formula for calculating the lock-up expansion index value is as follows:

[0051] ;

[0052] In the formula, This indicates the lock extension index value corresponding to the currently locked version of the core. This represents the overlap area between the suspected damaged patch area corresponding to the currently locked version kernel and the building kernel area. This represents the overlap area between the suspected damaged patch area corresponding to the currently locked version core and the building boundary buffer zone. This indicates the alternating confirmation constraint value for the object corresponding to the currently locked version core. This indicates that the constraint values ​​for alternating confirmation of adjacent building objects are represented. This indicates the relative post-earthquake time corresponding to the currently locked version of the core. This indicates the relative post-earthquake time for adjacent building objects.

[0053] In this implementation plan, by incorporating the overlapping area values ​​of the suspected damaged patch area and the building boundary buffer zone, the overlapping area values ​​of the suspected damaged patch area and the building core area, the object alternation confirmation constraint value, and the relative post-earthquake time into the calculation process of the locking expansion index value, the boundary expansion trend of the current locking version core, the degree of confirmation strength similarity between the current locking version core and adjacent building objects, and the degree of post-earthquake evolution synchronization are formed into a unified constraint relationship in the same quantitative judgment link. This improves the ability of the locking expansion index value to characterize the feasibility of cross-object expansion and enhances the stability, accuracy, and continuity of the subsequent block locking unit generation results.

[0054] Specifically, the steps for generating cross-scale version locking data are as follows: When the locking expansion index value is greater than the locking expansion threshold, the current building object and its adjacent building objects are merged into a block locking unit; wherein, the locking expansion threshold is preset based on the quantile statistical results of all locking expansion index values ​​in the current processing batch, and the locking expansion threshold value is between the 60th and 80th quantile values ​​of all locking expansion index values; in the process of merging the current building object and its adjacent building objects into a block locking unit when the locking expansion index value is greater than the locking threshold, the building census object number value and the locking version core value corresponding to the current building object are first read. First, the building object boundary area is read. Then, the building census object number value, lock version core, and building object boundary area corresponding to adjacent building objects are read. The current building object and adjacent building objects are written into the same block merging table. Then, all building census object numbers contained in the block merging table are used as the object member set of the block locking unit, and the joint outer boundary of all building object boundary areas in the block merging table is used as the boundary of the block locking unit. When generating the block locking unit, the lock version core corresponding to each building census object number value in the block merging table is also retained in the internal index table of the block locking unit, so that subsequent block-level comparisons can still be traced back to the block. The basic principle of building-level locking is that when the locking expansion index value is greater than the locking expansion threshold, it indicates that the current building object's boundary expansion trend has been able to break through the confirmation difference constraint and the temporal difference constraint, forming a stable expansion pull on adjacent building objects. Therefore, merging the current building object and adjacent building objects into the same block locking unit can upgrade the building-level expansion relationship to a block-level continuous locking relationship. Further comparison of adjacent block locking units is then performed. When there is a continuous connection chain between adjacent block locking units formed by the sequential connection of the boundary contact relationships of two or more building objects, and the relative seismicity corresponding to the locking version core in the later block locking unit is considered... When the time of the later block is later than the relative post-earthquake time corresponding to the locking version core in the previous block locking unit, the adjacent block locking units are merged into a regional locking unit, and the locking version core with the largest alternating confirmation constraint value of objects in the regional locking unit is determined as the regional dominant locking core; when comparing adjacent block locking units, the boundary of the previous block locking unit and the boundary of the next block locking unit are extracted first, and then the boundary contact relationship between the boundary areas of all building objects in the previous block locking unit and the boundary areas of all building objects in the next block locking unit is retrieved, and the boundary contact relationship that can be connected end to end to form a chain is written into a continuous linked list;Subsequently, the relative post-earthquake times corresponding to all locking version cores within the preceding block locking unit are extracted, and then the relative post-earthquake times corresponding to all locking version cores within the following block locking unit are extracted. The relative post-earthquake time corresponding to the latest-ranked locking version core in the preceding block locking unit is used as the comparison time for the preceding block locking unit, and the relative post-earthquake time corresponding to the earliest-ranked locking version core in the following block locking unit is used as the comparison time for the following block locking unit. When the comparison time of the following block locking unit is later than that of the preceding block locking unit, the preceding and following block locking units are written into the same area merging table. When merging adjacent block locking units into an area locking unit, all block locking units included in the area merging table are used as the member set of the area locking unit, and the joint outer boundary of all block locking unit boundaries is used as the boundary of the area locking unit. Then, the object intersection corresponding to all locking version cores in the area merging table is extracted. The locking version core with the largest alternating confirmation constraint value is written into the area-dominant locking core table as the area-dominant locking core. The technical principle is that the continuous connection chain represents the continuous spatial transmission channel between multiple block locking units. The relative post-earthquake time corresponding to the locking version core in the later block locking unit is later than the relative post-earthquake time corresponding to the locking version core in the previous block locking unit, which represents the temporal progression relationship of the post-earthquake state evolution across blocks. Therefore, by constraining the area merging process through spatial continuity conditions and temporal progression conditions, the area locking unit can be established on the basis of continuous expansion links. Then, by using the locking version core with the largest alternating confirmation constraint value as the area-dominant locking core, the area-level expression can always be anchored to the building object with the highest confirmation intensity within the current area. The locking version cores, block locking units, area locking units, and area-dominant locking cores corresponding to the building survey object number values ​​are summarized to construct cross-scale version locking data. When summarizing the lock version cores, block lock units, area lock units, and area-dominant lock cores corresponding to the ID values ​​of each building survey object, the corresponding lock version cores are first written into the building-level lock table using the building survey object ID value as the index key. Then, the attribution relationship between the building survey object ID value and the block lock unit is written into the block attribution table. Next, the attribution relationship between the block lock unit and the area lock unit is written into the area attribution table. Finally, the area-dominant lock core is written into the area-dominant table. Subsequently, association processing is performed according to the building survey object ID value, block lock unit ID, and area lock unit ID, outputting cross-scale version lock data. The technical principle is that by uniformly organizing building-level, block-level, and area-level lock relationships in the same data link, a hierarchical mapping relationship can be formed between the lock version core, block lock unit attribution, area lock unit attribution, and area-dominant lock core of a single building object. This provides direct input for the subsequent construction of building-level, block-level, and area-level expression base elements.

[0055] In this implementation plan, the block locking unit is generated based on the locking expansion index value, and the area locking unit is generated based on the continuous connection chain and the relative post-earthquake time sequence. Then, the locking version core with the largest object alternation confirmation constraint value is used as the area-dominant locking core. This enables the building survey object number value, locking version core, block locking unit, area locking unit, and area-dominant locking core to form a hierarchical mapping relationship in the same data link. This improves the unified representation ability of cross-scale version locking data for building-level state evolution, block-level continuous expansion, and area-level dominant expression, and enhances the continuity, stability, and traceability of the subsequent construction process of different levels of expression base elements.

[0056] Specifically, based on cross-scale version-locked data, different levels of representational base elements are constructed. Representation replacement, fission mapping processing, and time-series locking display are performed according to the maximum value of the constraint value confirmed by the alternating objects and the relative post-earthquake time. The specific steps for generating an adaptive multi-scale visualization map of urban earthquake prevention are as follows: Figure 4As shown, cross-scale version locking data is read, and the boundary areas of building objects, street blocks, and districts are used as the base elements for building-level, street block-level, and district-level expressions, respectively. When reading the cross-scale version locking data, the boundary areas of building objects corresponding to the building survey object numbers, the street block locking unit boundaries corresponding to the street block locking units, and the district locking unit boundaries corresponding to the district locking units are extracted first. Then, the alternating confirmation constraint values ​​of objects corresponding to the locking version kernel, the alternating confirmation constraint values ​​of objects corresponding to the district-dominant locking kernel, the relative post-earthquake time corresponding to the locking version kernel, and the relative post-earthquake time corresponding to the district-dominant locking kernel are extracted. The boundary areas of building objects are then written into... The architectural-level representation base table writes the street block locking unit boundary into the street block-level representation base table, and the area locking unit boundary into the area-level representation base table. This ensures that architectural-level, street block-level, and area-level representation bases correspond to architectural-level, street block-level, and area-level boundary carriers respectively within the same hierarchical link. The technical principle lies in first splitting the boundary carriers in the cross-scale version locking data hierarchically, and then performing subsequent map organization processing. This ensures that subsequent representation replacement, fission mapping processing, and time-series locking display are all based on clearly defined boundary units. Representation replacement is then performed on architectural-level, street block-level, and area-level representation bases, with the next level representation base being... When the maximum value of the corresponding object alternation confirmation constraint is greater than the maximum value of the object alternation confirmation constraint corresponding to the previous level expression element, and the relative post-earthquake time corresponding to the next level expression element is later than the relative post-earthquake time corresponding to the previous level expression element, the next level expression element replaces the previous level expression element; when the maximum value of the object alternation confirmation constraint corresponding to the previous level expression element is greater than or equal to the maximum value of the object alternation confirmation constraint corresponding to the next level expression element, the previous level expression element remains unchanged; when replacing building-level expression elements, street-level expression elements, and district-level expression elements respectively, the district-level expression element is first used as the previous level expression element and the street-level expression element is used as the next level expression element. The first hierarchical comparison is performed, and then the second hierarchical comparison is performed with the street-level expression base element as the upper-level expression base element and the building-level expression base element as the lower-level expression base element. When performing each hierarchical comparison, the maximum value of the alternating confirmation constraint value of the object corresponding to the lower-level expression base element is extracted first, and then the maximum value of the alternating confirmation constraint value of the object corresponding to the upper-level expression base element is extracted. The size judgment is performed, and then the relative post-earthquake time corresponding to the lower-level expression base element and the relative post-earthquake time corresponding to the upper-level expression base element are extracted. The order judgment is performed. When the lower-level expression base element meets the replacement condition in both the confirmation intensity and time progression judgment dimensions, the lower-level expression base element is written into the current display table, overwriting the original upper-level expression base element.The technical principle lies in the fact that the maximum value of the alternating confirmation constraint is used to characterize the strongest confirmation intensity within the current level boundary unit, while the relative post-earthquake time is used to characterize the corresponding time progression position of the current level boundary unit. The replacement is not executed directly based on a single hierarchical relationship or a single time sequence. Instead, it is triggered by the simultaneous fulfillment of two conditions: "the maximum value of the alternating confirmation constraint is larger" and "the relative post-earthquake time is later." This ensures that the next-level representation element only enters the map when it is superior to the previous-level representation element in both confirmation intensity and time progression updates, thereby avoiding... Simply relying on hierarchical refinement can lead to low-confidence local representations prematurely replacing higher-level representation units, while also avoiding the overwriting of high-confidence higher-level representation units by representations with insufficient confirmation strength due to later timing. After representation replacement, fission mapping is performed on each level of representation unit. When two or more building-level representation units exist within the same block-level representation unit, with alternating confirmation constraint values ​​greater than the fission threshold and corresponding relative post-earthquake time differences greater than the temporal difference threshold, the block-level representation unit is split into multiple building-level representation units. The fission threshold is determined based on the alternating confirmation constraint values ​​of objects within the current level of representation unit. The distribution results are pre-set, and the fission threshold is taken as the mean of the alternating confirmation constraint values ​​of the current level objects plus one standard deviation. When there are two or more block-level expression base elements within the same area-level expression base element whose alternating confirmation constraint values ​​are greater than the fission threshold and whose corresponding relative post-earthquake time differences are greater than the temporal difference threshold, the area-level expression base element is split into multiple block-level expression units. When performing fission mapping processing, first, the alternating confirmation constraint values ​​of all building-level expression base elements within the current block-level expression base element are retrieved, then the set of building-level expression base elements whose alternating confirmation constraint values ​​are greater than the fission threshold is selected. Extract the relative post-earthquake time corresponding to each building-level expression base element within the set, and calculate the difference between the relative post-earthquake times one by one. When at least two building-level expression base elements in the set satisfy the object alternation confirmation constraint value being greater than the fission threshold and the difference between the relative post-earthquake times being greater than the temporal difference threshold, the current block-level expression base element is re-split into multiple building-level expression units according to the boundary area of ​​the building object. The temporal difference threshold is preset based on the statistical results of the difference between the relative post-earthquake times within the current level expression base element, and the temporal difference threshold is 1 to 2 times the median of the difference between the relative post-earthquake times of the current level.Next, the system retrieves the object alternation confirmation constraint values ​​corresponding to all block-level expression base elements within the current area-level expression base element. Then, it filters out the block-level expression base element sets whose object alternation confirmation constraint values ​​are greater than the fission threshold. Finally, it extracts the relative post-earthquake times corresponding to each block-level expression base element within this set and calculates the difference between the relative post-earthquake times one by one. When at least two block-level expression base elements in the block-level expression base element set satisfy the conditions that the object alternation confirmation constraint value is greater than the fission threshold and the difference between the relative post-earthquake times is greater than the temporal difference threshold, the current area-level expression base element is re-splitted into multiple block-level expression units according to the block-level locking unit boundary. The technical principle is that fission mapping processing is not based on the number of lower-level expression units. The direct triggering of the quantity is not solely based on the object alternation confirmation constraint value being greater than the fission threshold. Instead, it binds the object alternation confirmation constraint value being greater than the fission threshold and the difference between the relative post-earthquake time and the temporal difference threshold as a joint fission condition. When multiple lower-level expression units with high confirmation intensity and significantly asynchronous time progression appear simultaneously within the same block-level expression unit or the same area-level expression unit, maintaining the original aggregation boundary would compress the post-earthquake states appearing at different times into the same expression unit, causing the map to mislead the order of state evolution within the block and area. Therefore, by binding the fission triggering with the temporal difference, the aggregation boundary can be made to exist only within the actual time. The data is only split during differentiation, ensuring that the map representation can synchronously reflect the spatial and temporal differentiation structures. After the fission mapping process is completed, each level of representation unit is displayed with temporal locking. When the map display request time is not less than the relative post-earthquake time corresponding to the locking version core or the dominant locking core of the area, the suspected damaged patch area and building object boundary area corresponding to the current representation unit are displayed. When the map display request time is less than the relative post-earthquake time corresponding to the locking version core or the dominant locking core of the area, the historical representation unit with the largest alternating confirmation constraint value in the corresponding level is displayed. When performing temporal locking display, the map display request time is read first, where the map display request time can be... The client sends input when initiating a map display request; then, according to the level to which the current expression unit belongs, the corresponding locking version relative to the post-earthquake time and the area's dominant locking core relative to the post-earthquake time are extracted, and the map display request time is compared with the current level's reference time; when the map display request time is not less than the current level's reference time, the suspected damaged patch area and building object boundary area corresponding to the current expression unit are written into the current map display table; when the map display request time is less than the current level's reference time, the current level's historical expression unit set is first retrieved, then the historical expression unit with the largest object alternating confirmation constraint value is extracted, and this historical expression unit is written into the current map display table;The technical principle lies in establishing a direct comparison relationship between the map display request time and the relative post-earthquake time corresponding to the locked version core and the relative post-earthquake time corresponding to the area-dominant locking core. This ensures that the map content always remains consistent with the timeline position. When the map display request time has not yet progressed to the time position corresponding to the current locking state, historical expression units are called to replace the current expression unit, which can prevent the map from prematurely displaying state results that have not yet occurred. By overlaying the building-level expression units, street-level expression units, and area-level expression units that have completed expression replacement, fission mapping processing, and time-series locking display, an adaptive multi-scale visualization map of urban earthquake prevention is generated. When overlaying architectural, street-level, and district-level expression units that have undergone expression replacement, fission mapping, and time-series locking, the layer overlay table is first written in the order of district-level, street-level, and architectural expression units. Then, layer overlay control is executed according to layer priority, current display status, and boundary coverage relationship, placing architectural expression units on the top layer, street-level expression units on the middle layer, and district-level expression units on the bottom layer. The result of all layer overlays is then output as an adaptive multi-scale visualization map of urban earthquake prevention. The technical principle is that by first completing layer replacement, boundary fission, and time-series filtering, and then executing layer overlay in a fixed layer order, the same map can simultaneously maintain the overall situational expression capability at the district level, the continuous expansion expression capability at the street level, and the detailed object expression capability at the architectural level, thereby achieving a multi-scale visualization map that automatically adjusts with changes in object status and time progression.

[0057] In this implementation plan, by unifying the boundary areas of building objects, the boundaries of block locking units, and the boundaries of area locking units into the construction process of different levels of expression base elements, and then incorporating the maximum value of the alternating confirmation constraint value of objects, the relative post-earthquake time, the locking version core, and the area-dominant locking core into the expression replacement, fission mapping processing, temporal locking display, and image synthesis links, the building-level expression units, block-level expression units, and area-level expression units form a hierarchical expression relationship in the same map generation link that can be dynamically adjusted according to changes in confirmed intensity, time progression, and spatial aggregation status. This improves the expression accuracy, map switching continuity, and temporal display consistency of the urban earthquake prevention single map adaptive multi-scale visualization map at different post-earthquake state levels.

[0058] like Figure 2As shown, the second aspect of this invention provides an adaptive multi-scale visualization mapping system for urban earthquake prevention using a single map. The system comprises: a data acquisition and processing module, an object verification and association module, a span version locking module, and a constraint-based mapping module. The data acquisition and processing module is used to construct a unified spatial reference coordinate system, collect spatial and temporal data corresponding to the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations, construct a multi-source earthquake prevention basic dataset, and preprocess the multi-source earthquake prevention basic data. The object verification and association module is used to construct building object boundary regions, building core regions, building boundary buffer zones, and suspected damaged patch regions based on the preprocessed multi-source earthquake prevention basic data, combining overlapping area values ​​and... The location-based relationship determines the building survey object number corresponding to the current record, and constructs an alternating confirmation segment according to the relative post-earthquake time, calculates the object alternating confirmation constraint value, and generates object-level seismic defense correlation data; the span-scale version locking module is used to determine the locking version kernel based on the object-level seismic defense correlation data, and calculates the locking expansion index value when there is a boundary contact relationship or suspected damage patch area overlap relationship between adjacent building objects, generating cross-scale version locking data; the constraint fission mapping module is used to construct different levels of expression base elements based on the cross-scale version locking data, and performs expression replacement, fission mapping processing and temporal locking display according to the maximum value of the object alternating confirmation constraint value and the relative post-earthquake time, generating an adaptive multi-scale visualization map of urban seismic defense.

[0059] In this implementation plan, by organizing multi-source basic data on earthquake prevention, object-level earthquake prevention correlation data, and cross-scale version locking data sequentially in the same technical link, and performing image synthesis processing for a unified urban earthquake prevention map on this basis, the building survey object number value, building object boundary area, building kernel area, building boundary buffer area, suspected damaged patch area, object alternating confirmation constraint value, locking version kernel, locking expansion index value, and different level expression base elements can form a continuous processing relationship, thereby improving the overall carrying capacity of the adaptive multi-scale visualization map of urban earthquake prevention for unified correlation results of multi-source objects, cross-scale state locking results, and temporal evolution expression results, and enhancing the stability, continuity, and traceability of the map generation process.

[0060] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0061] The preferred embodiments of the present invention disclosed above are merely illustrative of the invention. These preferred embodiments do not exhaustively describe all details, nor do they limit the invention to the specific implementations described. Clearly, many modifications and variations can be made based on the content of this specification. This specification selects and specifically describes these embodiments to better explain the principles and practical applications of the invention, thereby enabling those skilled in the art to better understand and utilize the invention. The invention is limited only by the claims and their full scope and equivalents.

Claims

1. An adaptive multi-scale visualization mapping method for urban earthquake prevention, characterized in that: Includes the following steps: S1. Construct a unified spatial reference coordinate system, collect spatial and temporal data corresponding to the location of building survey objects, remote sensing patch locations, IoT sensing device installation locations, and manually reported locations, construct a multi-source basic dataset for earthquake prevention, and preprocess the multi-source basic data for earthquake prevention. S2, based on the preprocessed earthquake prevention multi-source basic data, construct the building object boundary area, building core area, building boundary buffer area and suspected damaged patch area, combine the overlapping area value and the position falling relationship to determine the building survey object number value corresponding to the current record, and construct the alternating confirmation segment according to the relative post-earthquake time, calculate the object alternating confirmation constraint value, and generate object-level earthquake prevention association data; S3. Based on object-level seismic and anti-seismic correlation data, the locking version kernel is determined, and when there is a boundary contact relationship or a suspected damage patch area overlap relationship in the boundary area of ​​adjacent building objects, the locking expansion index value is calculated to generate cross-scale version locking data. S4 constructs different levels of representation elements based on cross-scale version locking data, and performs representation replacement, fission mapping processing and time-series locking display based on the maximum value of the constraint value of the object alternation confirmation and the relative post-earthquake time to generate an adaptive multi-scale visualization map of urban earthquake prevention.

2. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 1, characterized in that: The specific steps for constructing a unified spatial reference coordinate system, collecting spatial and temporal data corresponding to the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations, constructing a multi-source earthquake prevention basic dataset, and preprocessing the multi-source earthquake prevention basic data are as follows: Select the plane coordinate reference datum and elevation reference datum corresponding to the urban earthquake prevention area, determine the east-west coordinate axis and north-north axis with the plane coordinate reference datum, and determine the height reference axis with the elevation reference datum to construct a unified spatial reference coordinate system; Then, the locations of building survey objects, remote sensing patches, IoT sensing device installation locations, and manually reported locations are uniformly mapped to a unified spatial reference coordinate system. The following coordinates are collected: eastward coordinates of building base outline boundary points, northward coordinates of suspected damaged patch boundary points, eastward coordinates of IoT sensing device installation points, northward coordinates of manually reported locations, and building survey object numbers. Simultaneously, the following data are collected: earthquake event occurrence time, remote sensing acquisition time, IoT sensing sampling time, manually reported time, total building height, total building area, building tilt angle, structural crack width, building settlement, and building horizontal displacement. This data is used to construct a multi-source earthquake prevention basic dataset. For the collected multi-source basic data of earthquake prevention, a master-slave clock deviation compensation algorithm based on a precise time protocol is used to perform unified processing of the multi-source acquisition time. The Bursa seven-parameter coordinate transformation algorithm is used to perform spatial location benchmark unification processing; the median absolute deviation anomaly detection algorithm is used to perform abnormal data identification processing. The sliding median filtering algorithm is used to smooth the local abrupt earthquake prevention records; the range normalization algorithm is used to unify the numerical scale, and the preprocessed earthquake prevention multi-source basic data is output.

3. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 2, characterized in that: The specific steps for constructing building object boundary regions, building core regions, building boundary buffer zones, and suspected damaged patch regions based on preprocessed earthquake-resistant multi-source basic data, and determining the building survey object number value corresponding to the current record by combining the overlapping area value and the positional relationship, are as follows: Read the preprocessed earthquake prevention multi-source basic dataset, construct the building object boundary region according to the building survey object number value; determine the inward offset distance according to the product of the total building height value and the preset offset ratio, perform inward offset on the building object boundary region, and generate the building kernel region; Define the area between the boundary region of the building object and the building kernel region as the building boundary buffer; The eastward and northward coordinates of the boundary points of suspected damaged patches are collected according to the same remote sensing acquisition time value. The boundary points are then connected in sequence to form a closed loop, thus constructing a suspected damaged patch area. The overlap area between each suspected damaged patch area and the core area of ​​each building, as well as the overlap area with the buffer zone of each building, are calculated. The building survey object number value that satisfies the condition of having the largest patch overlap area value and whose IoT sensing device installation point or manually reported location falls within the core area or buffer zone of the building is determined as the building survey object number value corresponding to the current record.

4. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 3, characterized in that: The specific steps for constructing alternating confirmation segments based on relative post-earthquake times and calculating the alternating confirmation constraint values ​​of the objects are as follows: Based on the time of the earthquake event, the remote sensing acquisition time, IoT sensing sampling time, and manual reporting time are converted into relative post-earthquake time. All records corresponding to the same building object are collected according to the building survey object number value and arranged in ascending order of relative post-earthquake time. When adjacent records correspond to the same building survey object number value and the adjacent records contain three types of sources: remote sensing, IoT, and manual reporting, the adjacent records are merged into an alternating confirmation segment, and the first time difference between the remote sensing acquisition time value and the earthquake event occurrence time value, the second time difference between the IoT sensing sampling time value and the earthquake event occurrence time value, and the third time difference between the manual reporting time value and the earthquake event occurrence time value are calculated respectively. For each alternating confirmation segment, the first spatial constraint term is obtained by adding one to the overlapping area value of the suspected damaged patch area and the building core area and taking the square root. The second spatial constraint term is obtained by adding the natural constant e to the overlapping area of ​​the suspected damaged patch area and the building boundary buffer zone, and then taking the natural logarithm. The tilt angle value is increased by one and the square root is taken to obtain the tilt angle constraint term; the structural crack width value is increased by the natural constant e and the natural logarithm is taken to obtain the crack constraint term; the first time difference value, the second time difference value, and the third time difference value are increased by one and then multiplied in sequence, and the sum of the building settlement value and the building horizontal displacement value is divided by the total building height value plus one, and then added to the constant one to obtain the spatiotemporal suppression term; the first spatial constraint term, the second spatial constraint term, the tilt angle constraint term, and the crack constraint term are multiplied and then divided by the spatiotemporal suppression term to obtain the corresponding object alternation confirmation constraint value.

5. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 4, characterized in that: The specific steps for generating object-level earthquake prevention correlation data are as follows: For alternating confirmation segments under the same building survey object number, if the object alternating confirmation constraint value of the subsequent alternating confirmation segment is less than that of the preceding alternating confirmation segment, and at least one of the building tilt angle, structural crack width, building settlement, and building horizontal displacement values ​​of the subsequent alternating confirmation segment is greater than the corresponding values ​​of the preceding alternating confirmation segment, the subsequent alternating confirmation segment is marked as a reverse regression record. For alternating confirmation segments not marked as reverse regression records, comparisons are continued in the order of relative post-earthquake time. If the object alternating confirmation constraint value of the subsequent alternating confirmation segment is not less than that of the preceding alternating confirmation segment, and the time difference between the starting relative post-earthquake time of the subsequent alternating confirmation segment and the ending relative post-earthquake time of the preceding alternating confirmation segment is not greater than a continuous threshold, the subsequent alternating confirmation segment and the preceding alternating confirmation segment are compressed into a single object association record. When there are multiple single object association records under the same building survey object number, the single object association record with the largest object alternating confirmation constraint value is retained to construct object-level seismic defense association data.

6. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 1, characterized in that: The specific steps for determining the locking version core based on object-level seismic correlation data and calculating the locking expansion index value when there is a boundary contact relationship or overlapping relationship between suspected damaged patches in the boundary area of ​​adjacent building objects are as follows: Read the object-level seismic resistance association data and determine the single object association record corresponding to the number value of each building survey object as the locked version core of the current building object; Determine whether there is a common boundary line segment between the boundary areas of adjacent building objects, or whether the overlapping area between the corresponding suspected damaged patch areas is greater than the overlap threshold; for adjacent building objects that satisfy the boundary contact relationship or the suspected damaged patch area overlap relationship, read the overlap area value between the suspected damaged patch area and the building core area, the overlap area value between the suspected damaged patch area and the building boundary buffer area, the object alternating confirmation constraint value and the relative post-earthquake time, as well as the object alternating confirmation constraint value and the relative post-earthquake time corresponding to the adjacent building objects, and calculate the locking expansion index value of the current locking version core for adjacent building objects.

7. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 6, characterized in that: The specific steps for calculating the current locking version's lock extension index value relative to adjacent building objects are as follows: Add one to the overlap area between the suspected damaged patch area and the building boundary buffer zone corresponding to the current locked version core, then divide by the sum of the overlap area between the suspected damaged patch area and the building boundary buffer zone, the overlap area between the suspected damaged patch area and the building core area, and a constant one to obtain the buffer expansion ratio term; add one to the absolute value of the difference between the alternating confirmation constraint values ​​of the current locked version core and the adjacent building objects, and multiply by one the absolute value of the difference between the relative post-earthquake times of the current locked version core and the adjacent building objects to obtain the expansion suppression term; divide the buffer expansion ratio term by the expansion suppression term to obtain the locked expansion index value.

8. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 6, characterized in that: The specific steps for generating cross-scale version locking data are as follows: When the lock expansion index value is greater than the lock expansion threshold, the current building object and the adjacent building objects are merged into a block lock unit; The adjacent block locking units are further compared. When there is a continuous connection chain between adjacent block locking units, which is formed by the sequential connection of the boundary contact relationship of two or more building objects, and the relative post-earthquake time corresponding to the locking version core in the later block locking unit is later than the relative post-earthquake time corresponding to the locking version core in the previous block locking unit, the adjacent block locking units are merged into a district locking unit, and the locking version core with the largest object alternation confirmation constraint value in the district locking unit is determined as the district dominant locking core. By summarizing the locking version core, block locking unit, area locking unit, and area dominant locking core corresponding to the number values ​​of each building survey object, cross-scale version locking data is constructed.

9. The adaptive multi-scale visualization mapping method for urban earthquake prevention as described in claim 1, characterized in that: The specific steps for constructing different levels of representational base elements based on cross-scale version locking data, and performing representation replacement, fission mapping processing, and time-series locking display based on the maximum value of the constraint value of the alternating confirmation of objects and the relative post-earthquake time to generate an adaptive multi-scale visualization map of urban earthquake prevention are as follows: Read cross-scale version locking data and use the building object boundary area, the block locking unit boundary, and the area locking unit boundary as the building-level expression base element, the block-level expression base element, and the area-level expression base element, respectively. The building-level, street-level, and district-level expression base elements are replaced respectively. When the maximum value of the alternating confirmation constraint value of the object corresponding to the lower-level expression base element is greater than the maximum value of the alternating confirmation constraint value of the object corresponding to the upper-level expression base element, and the relative post-earthquake time corresponding to the lower-level expression base element is later than the relative post-earthquake time corresponding to the upper-level expression base element, the lower-level expression base element replaces the upper-level expression base element; when the maximum value of the alternating confirmation constraint value of the object corresponding to the upper-level expression base element is greater than or equal to the maximum value of the alternating confirmation constraint value of the object corresponding to the lower-level expression base element, the upper-level expression base element remains unchanged. After the expression replacement is completed, fission mapping processing is performed on the expression base elements at all levels. When there are two or more building-level expression base elements within the same block-level expression base element whose alternating confirmation constraint value is greater than the fission threshold and whose corresponding relative post-earthquake time difference is greater than the temporal difference threshold, the block-level expression base element is split into multiple building-level expression units. When there are two or more block-level expression base elements within the same area-level expression base element whose alternating confirmation constraint value is greater than the fission threshold and whose corresponding relative post-earthquake time difference is greater than the temporal difference threshold, the area-level expression base element is split into multiple block-level expression units. After the fission mapping process is completed, the time-series locking display is performed on the expression units at all levels. When the time requested for map display is not less than the relative post-earthquake time corresponding to the locking version core or the dominant locking core of the area, the suspected damaged patch area and the boundary area of ​​the building object corresponding to the current expression unit are displayed. When the time requested for map display is less than the relative post-earthquake time corresponding to the locking version core or the dominant locking core of the area, the historical expression unit with the largest object alternation confirmation constraint value in the corresponding level is displayed. The architectural, street, and district-level expression units, which have undergone expression replacement, fission mapping, and time-series locking display, are overlaid to generate an adaptive multi-scale visualization map of urban earthquake prevention.

10. An adaptive multi-scale visualization mapping system for urban earthquake prevention, characterized in that: include: The module includes a data acquisition and processing module, an object verification and association module, a span version locking module, and a constraint fission mapping module, among which: The data acquisition and processing module is used to construct a unified spatial reference coordinate system, collect spatial and temporal data corresponding to the location of building survey objects, remote sensing patch locations, IoT sensing device installation locations, and manually reported locations, construct a multi-source basic dataset for earthquake prevention, and preprocess the multi-source basic dataset for earthquake prevention. The object verification and association module is used to construct the building object boundary area, building core area, building boundary buffer area and suspected damage patch area based on the preprocessed earthquake prevention multi-source basic data. It determines the building survey object number value corresponding to the current record by combining the overlapping area value and the position falling relationship, and constructs the alternating confirmation segment according to the relative post-earthquake time, calculates the object alternating confirmation constraint value, and generates object-level earthquake prevention association data. The cross-scale version locking module is used to determine the locking version core based on object-level seismic and mitigation correlation data, and to calculate the locking expansion index value when there is a boundary contact relationship or a suspected damage patch area overlap relationship in the boundary area of ​​adjacent building objects, thereby generating cross-scale version locking data. The constrained fission mapping module is used to construct different levels of representation base elements based on cross-scale version locked data, and to perform representation replacement, fission mapping processing and time-series locked display according to the maximum value of the constraint value confirmed by the object alternation and the relative post-earthquake time, so as to generate an adaptive multi-scale visualization map of urban earthquake prevention.

Images