Spatial visualization display method and system combined with earthquake disaster database

Abstract

The present invention provides a spatial visualization display method and system combined with an earthquake disaster database. First, an earthquake disaster database containing multiple types of disaster record units with geographic coordinate marks is obtained, and then disaster feature extraction processing is performed on the earthquake disaster database to obtain disaster distribution association features and attribute association features of the disaster record units. A pre-built visualization mapping model is called to perform visualization rule mapping processing on the disaster distribution association features and attribute association features to generate a visualization rule set containing spatial position mapping rules and attribute value mapping rules. Based on this, a disaster spatial visualization display model is constructed for converting the disaster record units into visualization elements with spatial position identifiers and attribute value identifiers. Finally, by outputting earthquake disaster spatial visualization display results containing disaster distribution graphics of different geographic coordinate positions and association annotation information of corresponding attribute values, the intuitiveness and scientificity of earthquake disaster information display are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of earthquake disaster visualization, and in particular to a spatial visualization display method and system combined with an earthquake disaster database. Background Art

[0002] In the field of earthquake disaster response and research, accurately and intuitively presenting earthquake disaster information is crucial for rescue decision-making, post-disaster assessment, and scientific research. Traditional methods of presenting earthquake disaster information mainly rely on text reports, simple statistical charts, and two-dimensional map annotations.

[0003] Although written reports can record detailed information about the earthquake, such as source parameters, damage scope, and emergency response, the information is relatively scattered and lacks intuitiveness. Decision makers and researchers need to spend a lot of time and energy to sort out and analyze key information, making it difficult to quickly grasp the overall situation and spatial distribution characteristics of the earthquake disaster.

[0004] Simple statistical charts, such as bar charts and line charts, are mainly used to display certain statistical data on earthquake disasters, such as the number of people affected and economic losses in different regions. However, these charts cannot intuitively reflect the geographical distribution of earthquake disasters, nor can they clearly show the relationship between the disaster situation and the geographical location. They have limited guiding role in the rational deployment of rescue forces and post-disaster reconstruction planning.

[0005] Although two-dimensional map annotation can show the geographical location of earthquake disasters to a certain extent, it usually only simply marks the epicenter and the approximate scope of the disaster area. It lacks effective means to display detailed information about the disaster, such as the degree of damage at different locations and emergency response measures, making it difficult to meet the complex needs of earthquake disaster analysis and decision-making.

[0006] Furthermore, existing methods for displaying earthquake disaster data often overlook the correlations between different attributes within earthquake disaster data. There is a close, intrinsic connection between earthquake source information, damage severity, and emergency response. For example, the depth and intensity of the earthquake source directly influence the extent of damage in different regions, while the extent of damage, in turn, determines the level and measures of emergency response. Traditional display methods fail to intuitively present this interconnected information, hindering a comprehensive understanding and in-depth analysis of earthquake disasters. Summary of the Invention

[0007] In view of the above-mentioned problems, in combination with the first aspect of the present invention, an embodiment of the present invention provides a spatial visualization display method combined with an earthquake disaster database, the method comprising:

[0008] Acquire an earthquake disaster database, wherein the earthquake disaster database includes multiple types of disaster record units with geographic coordinate marks, and the disaster record units include earthquake source information records, damage extent records, and emergency response records;

[0009] Performing disaster feature extraction processing on the earthquake disaster database to obtain disaster distribution association features and attribute association features of the disaster record units, wherein the disaster distribution association features represent the correlation between disaster distributions at different geographic coordinate locations, and the attribute association features represent the attribute correlation between earthquake source information, damage degree, and emergency response;

[0010] Calling a pre-built visualization mapping model to perform visualization rule mapping processing on the disaster distribution association features and the attribute association features, and generating a visualization rule set including spatial position mapping rules and attribute value mapping rules;

[0011] Constructing a disaster situation spatial visualization display model based on the visualization rule set, wherein the disaster situation spatial visualization display model is used to convert disaster situation record units in the earthquake disaster situation database into visualization elements with spatial location identifiers and attribute value identifiers;

[0012] The spatial visualization display model for disaster conditions is used to output spatial visualization display results of earthquake disaster conditions, wherein the spatial visualization display results include disaster distribution graphics at different geographic coordinate positions and associated annotation information of corresponding attribute values.

[0013] On the other hand, an embodiment of the present invention also provides a spatial visualization display system combined with an earthquake disaster database, including a processor and a machine-readable storage medium, the machine-readable storage medium is connected to the processor, the machine-readable storage medium is used to store programs, instructions or codes, and the processor is used to execute the programs, instructions or codes in the machine-readable storage medium to implement the above method.

[0014] Based on the above aspects, by obtaining an earthquake disaster database containing multiple types of disaster record units with geographic coordinate marks, and comprehensively integrating key information such as earthquake source information records, damage degree records and emergency response records, the earthquake disaster database is subjected to disaster feature extraction and processing, and the disaster distribution correlation characteristics and attribute correlation characteristics of the disaster record units are obtained. The distribution pattern of earthquake disasters in geographic space and the intrinsic connection between earthquake source information, damage degree and emergency response are deeply explored.

[0015] The pre-built visualization mapping model is called to perform visualization rule mapping processing on the disaster distribution association characteristics and attribute association characteristics, generating a visualization rule set including spatial position mapping rules and attribute value mapping rules, realizing the formulation of conversion rules from abstract disaster characteristics to intuitive visualization elements, ensuring that the visualization display can accurately and effectively present the spatial distribution and attribute information of the earthquake disaster.

[0016] Based on a set of visualization rules, a disaster spatial visualization display model is constructed, which converts the disaster record units in the earthquake disaster database into visualization elements with spatial location identifiers and attribute value identifiers. It effectively integrates spatial information and attribute information, allowing users to fully understand all aspects of the earthquake disaster on a unified platform.

[0017] The spatial visualization model outputs a spatial visualization of earthquake disasters, including disaster distribution graphics at different geographic coordinate locations and associated annotation information for corresponding attribute values. This intuitive and clear presentation presents both the overall situation and detailed characteristics of the earthquake disaster. Decision-makers and researchers can quickly and accurately grasp the spatial distribution of earthquake disasters, understand the extent of damage in different regions, and assess emergency response measures. This enables efficient utilization and intuitive presentation of earthquake disaster information, significantly improving the scientific nature and efficiency of earthquake disaster analysis and decision-making. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] Figure 1 It is a schematic diagram of the execution flow of the spatial visualization display method combined with the earthquake disaster database provided by an embodiment of the present invention.

[0019] Figure 2 Schematic diagram of exemplary hardware and software components of a spatial visualization display system combined with an earthquake disaster database provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0020] The present invention will be described in detail below with reference to the accompanying drawings. Figure 1 This is a flow chart of a spatial visualization display method combined with an earthquake disaster database provided by an embodiment of the present invention. The spatial visualization display method combined with an earthquake disaster database is introduced in detail below.

[0021] Step S110: Acquire an earthquake disaster database, wherein the earthquake disaster database includes multiple types of disaster record units with geographic coordinate marks, and the disaster record units include earthquake source information records, damage extent records, and emergency response records.

[0022] In scenarios related to earthquake disaster response, a predefined architecture is used to collect and integrate various earthquake-related information. Data is continuously output from multiple geographically dispersed earthquake monitoring stations, equipped with appropriate seismic monitoring instruments that can capture various parameters generated by earthquakes in real time. Simultaneously, information is fed back from the affected areas through various channels, with on-site personnel using tools such as handheld smart devices and vehicle-mounted survey equipment to survey the affected areas.

[0023] Data collected by earthquake monitoring stations constitute the bulk of earthquake source information records. These records precisely identify the geographic coordinates of the earthquake source, acquired through a combination of calibrated satellite positioning and ground-based measurement techniques. They also include key data such as focal depth and magnitude. All data is collected in strict accordance with international and industry standards and specifications.

[0024] In the information provided by on-site disaster-stricken areas, the damage extent records detail the damage to buildings in different locations. This includes structural types (frame, brick-concrete, steel), as well as the specific locations of damage, such as the location and extent of wall cracks and the extent of roof collapse. In addition, casualties are broken down by age group and gender. Emergency response records detail the precise time the emergency response was initiated, accurate to the second, and the deployment of rescue teams during the response, including the number of specialized rescue teams, their arrival times, and the types, quantities, and locations of supplies distributed.

[0025] The above information from different channels is integrated into the earthquake disaster database. Each disaster record unit in the earthquake disaster database is clearly marked with geographical coordinates.

[0026] Step S120: Perform disaster feature extraction processing on the earthquake disaster database to obtain the disaster distribution association features and attribute association features of the disaster record unit. The disaster distribution association features represent the disaster distribution correlation at different geographic coordinate positions, and the attribute association features represent the attribute correlation between earthquake source information, degree of damage and emergency response.

[0027] Step S121: performing geographic coordinate standardization processing on the disaster record units in the earthquake disaster database, unifying the geographic coordinate representation formats of different disaster record units, and obtaining standardized disaster record units with a unified coordinate reference, wherein the geographic coordinate representation formats include latitude and longitude formats and plane rectangular coordinate formats.

[0028] Step S1211: Identify the geographic coordinate representation format identifier of each disaster record unit in the earthquake disaster database, where the representation format identifier includes a latitude and longitude identifier and a plane rectangular coordinate identifier.

[0029] In this embodiment, for each disaster record unit in the earthquake disaster database, specific data fields or data structure features within the record unit are analyzed based on pre-set format recognition rules. If a specific field in a disaster record unit presents a data pattern that conforms to longitude and latitude identification, such as containing clear latitude and longitude symbols, and the value range is consistent with the conventional longitude and latitude value range, then the geographic coordinate representation format of the unit is determined to be longitude and latitude format; if the data structure presents the form of coordinate values ​​in a rectangular coordinate system, has clear X and Y coordinate identifications, and the coordinate units conform to the common settings of a plane rectangular coordinate system, then it is determined to be in plane rectangular coordinate format.

[0030] Step S1212: For the disaster record unit identified as being in the latitude and longitude format, extract its latitude and longitude values, and call the geographic coordinate conversion function to convert the latitude and longitude values ​​into coordinate values ​​of a plane rectangular coordinate system.

[0031] Once a disaster record unit is determined to be in latitude and longitude format, the latitude and longitude values ​​of the unit are accurately extracted according to the storage structure and data format of the latitude and longitude in the record unit. Assuming that the latitude and longitude values ​​are stored in a specific field order, they are retrieved in this order.

[0032] After obtaining the latitude and longitude values, a geographic coordinate conversion function based on map projection principles is called. This geographic coordinate conversion function takes into account the shape of the Earth's ellipsoid and the parameter settings of different map projection methods. Taking the Gauss-Krüger projection as an example, the relevant parameters on the ellipsoid are first calculated based on the Earth's semi-major and semi-minor axes, combined with the obtained longitude and latitude values. Then, through a series of mathematical transformations such as coordinate translation, rotation, and scaling, the longitude and latitude values ​​are gradually converted to coordinate values ​​in a plane rectangular coordinate system. After multiple calculation steps and parameter adjustments, the accurate coordinate values ​​in the plane rectangular coordinate system are guaranteed.

[0033] Step S1213: extracting the original coordinate values ​​of the disaster record units marked as being in plane rectangular coordinate format, checking whether the coordinate units are in a unified unit system, and performing unit conversion processing on the coordinate values ​​in a non-unified unit system.

[0034] For disaster record units identified as plane rectangular coordinates, the original coordinate values ​​are extracted based on their storage structure. After extraction, the unit identification corresponding to each coordinate value is compared with the unified unit system standard according to the unit inspection rules. If the units of some coordinate values ​​are found to be inconsistent with the standard unit system, such as some in meters and some in kilometers, the non-unified unit system coordinate values ​​are converted based on the meter-kilometer conversion relationship (1 kilometer = 1000 meters). The coordinate values ​​in kilometers are multiplied by 1000 to convert them to meters, achieving unified units for all plane rectangular coordinates.

[0035] Step S1214: performing coordinate system one processing on the converted plane rectangular coordinates and the original plane rectangular coordinates, selecting the center position of the area covered by the earthquake disaster database as the reference origin, and constructing a local plane rectangular coordinate system.

[0036] For example, step S1214-1: calculate the center position of all disaster record units in the latitude and longitude format in the earthquake disaster database, and the center position is determined by the statistical characteristics of the latitude and longitude values.

[0037] The latitude values ​​of all disaster record units in the longitude and latitude format in the earthquake disaster data database are summarized, and the central trend in the latitude direction is preliminarily determined by calculating the average of these latitude values. Similarly, all longitude values ​​are summarized and the average is calculated to determine the central trend in the longitude direction. The longitude and latitude corresponding to these two average values ​​are the preliminarily determined center position. However, considering the possible uneven distribution of the data, it is necessary to further analyze the degree of dispersion of the longitude and latitude values, such as calculating the variance of the latitude and longitude values. If the variance is large, it means that the data distribution is relatively dispersed. At this time, it may be necessary to recalculate the center position using a weighted average method. The weight can be determined based on the distance between the data point and the preliminarily calculated center position. The closer the distance, the higher the weight. Finally, the longitude and latitude values ​​that can more accurately represent the central trend of the database coverage area are obtained, and the center position is determined by this.

[0038] Step S1214 - 2 : Calling a geographic coordinate conversion function to convert the center position into a reference coordinate value of a global plane rectangular coordinate system.

[0039] The previously used geographic coordinate conversion function for converting longitude and latitude to rectangular coordinates is used as input, taking the longitude and latitude values ​​of the determined center location. The geographic coordinate conversion function performs complex calculations based on knowledge of Earth shape models and map projection principles. For example, if the Gauss-Krüger projection is used, the function will perform a series of mathematical transformations on the input longitude and latitude values, including translation, rotation, and scaling, taking into account parameters such as the Earth's semi-major and semi-minor axes. The final output is the corresponding base coordinate value in a global rectangular coordinate system.

[0040] Step S1214 - 3 : performing a translation transformation on the converted plane rectangular coordinates and the original plane rectangular coordinates, subtracting the reference coordinate value from the coordinate value of each coordinate point to obtain the local plane rectangular coordinate relative to the center position.

[0041] For all plane rectangular coordinates, whether converted from longitude and latitude or originally in plane rectangular coordinate format, subtract the X component of the reference coordinate value from its X coordinate value, and subtract the Y component of the reference coordinate value from its Y coordinate value. For example, if the plane rectangular coordinates of a coordinate point are (X1, Y1) and the reference coordinate value is (X0, Y0), then the local plane rectangular coordinates obtained after the translation transformation are (X1-X0, Y1-Y0). In this way, all coordinate points are converted to a relative coordinate system with the center position as the reference and the center position as the origin.

[0042] Step S1214 - 4 : Determine the coordinate axis direction of the local plane rectangular coordinate system, where the coordinate axis direction is determined by the directional characteristics of the geographic space.

[0043] Based on the directional characteristics of geographic space, the Y-axis direction is determined by the north-south direction of geographic space, with north as the positive direction of the Y-axis and south as the negative direction of the Y-axis. The X-axis direction is determined by the east-west direction of geographic space, with east as the positive direction of the X-axis and west as the negative direction of the X-axis. In this way, based on the natural directional characteristics of geographic space, the coordinate axis directions of the local plane rectangular coordinate system are clarified.

[0044] Step S1214-5: Perform sign correction on the local plane rectangular coordinates so that the coordinate value in the target direction is positive and the opposite direction is negative.

[0045] Check the local plane rectangular coordinate values ​​after the translation transformation. For the X coordinate, if the actual direction corresponding to the coordinate value is west, and the specified positive direction is east, then the X coordinate value is reversed to make its sign consistent with the specified positive direction identification. For the Y coordinate, if the actual direction corresponding to the coordinate value is south, and the specified positive direction is north, then the Y coordinate value is also reversed. Through this sign correction, ensure that the signs of all coordinate values ​​are consistent with the determined positive direction of the coordinate axis.

[0046] Step S1214-6: Perform unit unification processing on the corrected local plane rectangular coordinates, convert the units of all coordinate values ​​into a unified unit system, and generate a local plane rectangular coordinate system with a unified reference origin, unified coordinate axis direction and unified unit.

[0047] Check the units of the corrected local plane rectangular coordinate values. If there are multiple units, such as some coordinate values ​​in meters and some in kilometers, convert all coordinate values ​​to a unified unit based on the unit conversion relationship, for example, convert them all to meters. If other units are involved, such as centimeters, convert them according to the corresponding conversion ratio. By unifying the units, the local plane rectangular coordinate system not only has a unified reference origin and coordinate axis directions, but also has unified coordinate units, forming a complete and standardized local plane rectangular coordinate system.

[0048] Step S1214-7: converting the geographic coordinates of all disaster record units in the earthquake disaster database into coordinate values ​​of the local plane rectangular coordinate system to generate standardized disaster record units with a unified coordinate reference.

[0049] Traverse all disaster record units in the earthquake disaster database. Regardless of whether their original geographic coordinates are in latitude and longitude format or plane rectangular coordinate format, they are converted according to the conversion rules of the local plane rectangular coordinate system constructed previously. For those in latitude and longitude format, first convert them into plane rectangular coordinates according to the previous conversion process, and then perform operations such as translation, symbol correction, and unit unification; for those that are originally in plane rectangular coordinate format, directly perform translation, symbol correction, and unit unification. After these steps, the geographic coordinates of all disaster record units are converted into coordinate values ​​under the local plane rectangular coordinate system, thereby generating standardized disaster record units with a unified coordinate base.

[0050] Step S1215: performing a translation transformation on the plane rectangular coordinates of all disaster record units relative to the reference origin, adjusting the offset of the coordinate value, and making the coordinate value of the reference origin the reference point.

[0051] After constructing the local rectangular coordinate system, the rectangular coordinates of all disaster record units are further processed. For each disaster record unit's rectangular coordinate, its X coordinate value is subtracted from the X coordinate value of the reference origin in the global rectangular coordinate system, and its Y coordinate value is subtracted from the Y coordinate value of the reference origin in the global rectangular coordinate system. This operation adjusts the coordinate value offset so that the coordinate value of the reference origin in the newly constructed local rectangular coordinate system becomes the set reference point, such as setting the reference point coordinates to (0, 0).

[0052] Step S1216: performing precision uniformity processing on the coordinate values ​​after the translation transformation, retaining the coordinate values ​​to a uniform precision digit, and generating a standardized disaster record unit with a uniform coordinate reference, uniform unit, and uniform precision.

[0053] The coordinate values ​​after translation transformation are processed for uniform accuracy according to pre-set accuracy rules. A certain number of decimal places is set, and all coordinate values ​​are rounded to a uniform precision. This generates standardized disaster record units with a unified coordinate base, unified units, and unified precision. This ensures that the geographic coordinates of all disaster record units are unified in format, unit, and precision, facilitating subsequent analysis and processing.

[0054] Step S122: performing spatial neighborhood analysis on the standardized disaster record unit, calculating the overlapping relationship between the disaster record coverage area of ​​each geographic coordinate position and the adjacent geographic coordinate position, and generating a neighborhood overlap feature representing the spatial adjacency strength.

[0055] A spatial neighborhood analysis is performed on the disaster record units after coordinate standardization. Taking the geographic coordinate location of each standardized disaster record unit as the core, the adjacent geographic coordinate locations are determined based on pre-defined neighborhood range rules. For example, a certain coordinate location is set as the center, and other coordinate locations within a set distance range are considered adjacent locations.

[0056] After determining adjacent locations, conventional regional analysis methods from related technologies are used to determine the boundaries of the disaster record coverage area corresponding to each geographic coordinate location. The boundaries of the disaster record coverage areas at adjacent coordinate locations are then compared to determine overlap, such as whether any areas overlap. If so, the extent of the overlap is determined.

[0057] Based on the overlapping relationship, conventional computational logic in related technologies is used to generate a neighborhood overlap feature representing the strength of spatial adjacency. This computational logic may involve calculating the ratio of the area of ​​the overlapping region to the area of ​​the individual regions to quantify the strength of spatial adjacency. For example, if the overlapping area of ​​a region with an adjacent region accounts for a large proportion of its own area, it indicates that the spatial adjacency strength of the region with the adjacent region is high, and the generated neighborhood overlap feature value is correspondingly large.

[0058] Step S123: performing time series alignment processing on the standardized disaster record unit, mapping the occurrence time of the earthquake source information record, the assessment time of the damage extent record, and the start time of the emergency response record to the same time axis, and generating a time series aligned record unit.

[0059] Time series alignment is performed on the time information contained in the earthquake source information records, damage extent records, and emergency response records in the standardized disaster record unit. The time information in each record is first clarified, and can be represented in different time dimensions such as year, month, day, hour, minute, and second.

[0060] Then, based on a unified time standard, the occurrence time of earthquake source information records, the assessment time of damage extent records, and the initiation time of emergency response records are aligned and mapped. For example, all time information is converted into a timestamp calculated from a specific starting point in time. This method places the time information of different types of records on the same timeline, generating time-aligned record units, which facilitates subsequent comparison and analysis of the relationships between different types of records in the temporal dimension.

[0061] Step S124: performing attribute co-occurrence analysis on the time-series aligned record unit, counting the co-occurrence frequency of the characteristic parameters of the earthquake source information record and the characteristic parameters of the damage degree record, and generating a co-occurrence frequency feature representing the attribute association strength.

[0062] After generating time-aligned record units, attribute co-occurrence analysis was performed to determine the characteristic parameters contained in the earthquake source information records, such as magnitude and focal depth, and the characteristic parameters contained in the damage extent records, such as the type of damaged buildings and the number of casualties.

[0063] Traverse all time-aligned recording units and count the number of co-occurrences of characteristic parameters recorded for earthquake source information and characteristic parameters recorded for damage severity. For example, if a certain magnitude earthquake and a certain type of building damage occur simultaneously in a recording unit, a co-occurrence is recorded.

[0064] The co-occurrence counts, combined with the total number of recorded units, generate a co-occurrence frequency feature representing the strength of attribute association. For example, the co-occurrence frequency is calculated by dividing the co-occurrence count by the total number of recorded units. A higher co-occurrence frequency indicates a stronger association between the earthquake source information record and the damage severity record.

[0065] Step S125: Perform feature weighted fusion processing on the neighborhood overlapping features and the co-occurrence frequency features to obtain the disaster distribution association features and attribute association features of the disaster record unit, wherein the disaster distribution association features include the weighted results of the neighborhood overlapping features, and the attribute association features include the weighted results of the co-occurrence frequency features, wherein the weight of the neighborhood overlapping features is determined by the distribution density of the disaster records at the geographic coordinate location, and the weight of the co-occurrence frequency features is determined by the statistical significance of the attribute co-occurrence frequency.

[0066] After obtaining the neighborhood overlap and co-occurrence frequency features, a weighted fusion process is performed on these features. For example, the weight of the neighborhood overlap feature is first determined by analyzing the distribution density of disaster records at the geographic coordinate location. For example, the number of disaster records within a set range of a certain geographic coordinate location is counted. If the number is large, i.e., the density of disaster records is high, then the neighborhood overlap feature is assigned a higher weight; otherwise, it is assigned a lower weight.

[0067] The weight of the co-occurrence frequency feature is determined based on the statistical significance of the attribute co-occurrence frequency. For example, through statistical analysis, determine whether the co-occurrence frequency is significant in the overall data. If the co-occurrence frequency varies stably and with little variation across different record units, indicating high statistical significance, the co-occurrence frequency feature is given a higher weight; otherwise, it is given a lower weight.

[0068] Based on the determined weights, the neighborhood overlap features and co-occurrence frequency features are weighted. For example, the neighborhood overlap feature value is multiplied by its weight, and the co-occurrence frequency feature value is multiplied by its weight to obtain a weighted result. The weighted neighborhood overlap feature result is used as part of the disaster distribution association feature, and the weighted co-occurrence frequency feature result is used as part of the attribute association feature, thereby obtaining the disaster distribution association feature and attribute association feature of the disaster record unit.

[0069] Step S130: calling a pre-built visualization mapping model to perform visualization rule mapping processing on the disaster distribution association features and the attribute association features, and generating a visualization rule set including spatial position mapping rules and attribute value mapping rules.

[0070] Step S131: Input the disaster distribution association characteristics into the spatial mapping submodule of the visualization mapping model. The spatial mapping submodule generates a spatial position mapping rule representing the correspondence between the geographic coordinate position and the center coordinate of the visualization graphic by learning the historical mapping relationship between the geographic coordinate position and the visualization interface layout.

[0071] The pre-built visualization mapping model includes a spatial mapping submodule. The previously derived disaster distribution correlation features are input into this submodule. This submodule stores historical mapping data between geographic coordinate locations and visualization interface layouts. This historical mapping data records how different geographic coordinate locations were mapped to visualization interface layouts in similar scenarios in the past.

[0072] The spatial mapping submodule analyzes and learns from this historical mapping data, for example, using data analysis algorithms to uncover potential relationships between geographic coordinates and the center coordinates of visual graphics. After this analysis and learning, it generates spatial mapping rules that represent the correspondence between geographic coordinates and the center coordinates of visual graphics. For example, a predefined functional relationship or correspondence table is determined so that, given a geographic coordinate position, the corresponding center coordinate position in the visual graphics can be determined according to this spatial mapping rule.

[0073] Step S132: Input the attribute association features into the attribute mapping submodule of the visualization mapping model. The attribute mapping submodule generates attribute value mapping rules representing the correspondence between the degree of damage and the graphic color features, and the correspondence between the emergency response start time and the graphic border features by learning the historical association between the attribute values ​​and the visualization graphic styles.

[0074] The attribute association features are input into the attribute mapping submodule of the visualization mapping model. The attribute mapping submodule stores historical association data between attribute values ​​and visualization graphic styles. These historical association data record how different attribute values ​​were associated with visualization graphic styles in the past.

[0075] The attribute mapping submodule analyzes this historical correlation data, applying data mining techniques, for example, to explore the inherent connections between attribute values ​​and visual graphic styles. Through this analysis and learning, it generates attribute value mapping rules that represent the correspondence between damage severity levels and graphic color features, and between emergency response initiation times and graphic border features. For example, it determines which colors correspond to different damage severity levels, and which graphic border features correspond to different emergency response initiation times, forming clear mapping rules.

[0076] Step S133: performing density balancing adjustment processing on the spatial position mapping rule, calculating the graphic distribution density of each area in the visualization interface, and adjusting the direction of the graphic center coordinates of the area whose density exceeds the preset balancing range. The adjustment direction is determined by the regional density distribution trend.

[0077] Step S133 - 1 : Divide the visual interface into grid-like density calculation areas, where the size of each density calculation area is determined by the overall size of the interface and the density calculation accuracy.

[0078] Consider the overall width and height of the visualization interface and divide the grid according to the pre-set density calculation accuracy requirements. For example, if you want to calculate the distribution density of the graphic more accurately, you can divide the visualization interface into smaller grid regions; if the accuracy requirements are relatively low, the grid region size can be appropriately increased. Assuming the width of the visualization interface is W and the height is H, based on the accuracy settings, divide the width into N equal parts and the height into M equal parts. This divides the visualization interface into N×M grid-like density calculation regions, each with a width of W / N and a height of H / M, forming a regular grid structure for subsequent density calculations.

[0079] Step S133 - 2 : Count the number of visual graphics in each density calculation area, and calculate the ratio of the number of graphics to the area of ​​the area as the area density parameter.

[0080] For each grid-like density calculation region, we traverse all visual shapes and determine whether the center coordinates of each shape fall within the region. If the center coordinates of a shape fall within the region, the number of shapes in that region is incremented by 1. After counting the number of shapes in each region, we calculate the region's area, which is the product of the region's width and height (for the previously divided regions, the area is (W / N) × (H / M)). The ratio of the number of shapes in each region divided by the region's area is the density parameter for that region, which quantifies the density of the visual shapes within each region.

[0081] Step S133 - 3 : Identify a first density region whose regional density parameter exceeds a density balancing range, and a second density region outside the first density region.

[0082] A density balance range is pre-set, for example, with a lower limit of Dmin and an upper limit of Dmax. The density parameters of each density calculation area are compared with this balance range. If the density parameter of an area is greater than Dmax or less than Dmin, the area is identified as a first density area, i.e., an abnormal density area. Areas with a density parameter between Dmin and Dmax are identified as second density areas, i.e., relatively balanced density areas. This method allows the various areas of the visualization interface to be categorized by density, allowing for subsequent implementation of different treatment measures for different density areas.

[0083] Step S133 - 4 : sorting the visualized graphics within the first density area to obtain sorted graphics, wherein the sorting rule is determined by the spatial distribution characteristics of the center coordinates of the graphics.

[0084] For the visualized graphics within the first density region, analyze the spatial distribution of the center coordinates of the graphics. If the center coordinates of the graphics exhibit a linear distribution trend, such as gradually becoming denser from left to right or from top to bottom, sort them according to that trend. For example, if the graphics are vertically concentrated and gradually denser from top to bottom, sort them by the Y value of the center coordinates from smallest to largest. If the graphics are horizontally concentrated and gradually denser from left to right, sort them by the X value of the center coordinates from smallest to largest.

[0085] Step S133-5: Directionally adjust the sorted graphics in sequence. The adjustment direction is determined by the spatial distribution trend of the first density region, and the adjustment amount is determined by the difference between the regional density parameter and the balanced range.

[0086] Observe the spatial distribution trend of the first density area. If the graphics in this area gradually become denser from left to right and the density exceeds the upper limit, adjust the center coordinates of the graphics to the right to balance the density. For each sorted graphic, determine the adjustment amount based on the difference between the density parameter of the area and the density balance upper limit Dmax. The larger the difference, the larger the adjustment amount. For example, if the density parameter of a first density area is D1, and D1-Dmax=ΔD, an adjustment coefficient k can be set, and the adjustment amount is k×ΔD. According to this rule, the direction of the center coordinates of the sorted graphics is adjusted in sequence to make the distribution of graphics in the area more uniform.

[0087] Step S133-6: Supplement the position of the visualized graphics in the second density area, where the supplement position is determined by the spatial distribution characteristics of the second density area.

[0088] Analyze the spatial distribution of the second-density area. If you find that graphics are relatively sparse in a certain area, such as a corner or near a border, determine the location of additional graphics based on the shape of the area and the distribution of existing graphics. You can choose to add additional graphics in the center of the sparse area or near the border without affecting the overall layout. For example, if the second-density area is rectangular and the upper right corner is sparsely populated, you can determine a supplementary location near the center of the upper right corner to increase the graphics density in this area and make the graphics distribution more balanced across the entire visualization interface.

[0089] Step S133-7: recalculate the regional density parameters after direction adjustment and position supplementation, and repeat the above adjustment process for the regions still outside the density balance range until the density parameters of all regions are within the density balance range.

[0090] After adjusting the orientation of the first density region graphic and repositioning the second density region graphic, the method of step S133-2 is repeated to count the number of visual graphics within each density calculation region and calculate the ratio of the number of graphics to the region's area to obtain new regional density parameters. These new density parameters are compared with the density balance range (Dmin, Dmax). If there are still areas where the density parameters exceed this range, the adjustment process of steps S133-4 to S133-6 is repeated for these areas. This process continues, continuously optimizing the distribution of graphics, until the density parameters of all density calculation regions are within the pre-set density balance range, ensuring that the distribution of graphics in each area on the visualization interface is relatively balanced.

[0091] Step S133 - 8 : Record the final adjusted positions of all graphics, and generate adjusted spatial position mapping rules, wherein the spatial position mapping rules include a correspondence between the original center coordinates and the final adjusted coordinates.

[0092] Traverse all visual graphics and record the final center coordinate position of each graphic after the above adjustment. Associate the original center coordinates of each graphic with the final adjusted center coordinates to form a one-to-one correspondence. For example, the original center coordinates of graphic A are (X1, Y1), and the final center coordinates after adjustment are (X2, Y2). In the spatial position mapping rule, record this correspondence (X1, Y1) → (X2, Y2). By recording the above correspondences for all graphics, the adjusted spatial position mapping rule is generated.

[0093] Step S134: performing conflict resolution processing on the attribute value mapping rules, identifying attribute value combinations with overlapping color features or confused border features, and reallocating styles by adjusting color hue parameters or border line type parameters to make the visualization styles of different attribute values ​​distinguishable.

[0094] Perform conflict resolution on the generated attribute value mapping rules. Check the visualization styles corresponding to different attribute value combinations in the attribute value mapping rules, particularly color and border features. Identify any overlap in color features (i.e., different attribute values ​​corresponding to the same or similar colors) or border feature confusion (i.e., different attribute values ​​corresponding to indistinguishable border line types).

[0095] For attribute value combinations with overlapping color features, styles are reallocated by adjusting the color hue parameters. For example, similar colors can be visually differentiated by changing their hue. For attribute value combinations with confusing border features, styles are reallocated by adjusting border line parameters, such as changing a solid border to a dashed one or changing the spacing between dashed lines. This ensures that the visualization styles of different attribute values ​​are distinguishable and avoids confusion during visualization.

[0096] Step S135: Verify the consistency of the adjusted spatial position mapping rules and attribute value mapping rules, check whether the correspondence between the geographic coordinate position and the graphic center coordinate covers all disaster record units, verify whether the correspondence between the attribute value and the graphic style completely maps all attribute categories, and generate the final visualization rule set.

[0097] Verify the consistency of the spatial location mapping rules that have undergone density equalization adjustment and the attribute value mapping rules that have undergone conflict resolution. Check the correspondence between the geographic coordinates and the center coordinates of the graphics in the spatial location mapping rules to confirm whether all disaster record units in the earthquake disaster database are covered. That is, for each disaster record unit's geographic coordinate position, the corresponding center coordinates of the graphics can be found in this correspondence.

[0098] At the same time, the mapping rules for attribute values ​​and graphic styles were verified to ensure that they fully mapped all attribute categories. This ensured that every possible attribute value, whether in earthquake source information records, damage extent records, or emergency response records, had a clear and unique corresponding visualization graphic style. After verification, the adjusted spatial location mapping rules and attribute value mapping rules were integrated to generate the final set of visualization rules.

[0099] Step S140: constructing a disaster spatial visualization display model based on the visualization rule set, wherein the disaster spatial visualization display model is used to convert the disaster record units in the earthquake disaster database into visualization elements with spatial location identifiers and attribute value identifiers.

[0100] Step S141: Construct a visualization element generation module, which is used to receive disaster record units in the earthquake disaster database, extract the graphic center coordinates corresponding to the geographic coordinate position according to the spatial position mapping rules in the visualization rule set, extract the color features corresponding to the degree of damage and the border features corresponding to the emergency response initiation time according to the attribute value mapping rules, and generate basic visualization element information including center coordinates, color features and border features.

[0101] Disaster record units are read from the earthquake disaster database. These disaster record units carry a wealth of earthquake-related information, among which geographic coordinates are one of the key pieces of information. According to the spatial position mapping rules generated previously, for the geographic coordinates of each disaster record unit, the corresponding relationship defined in the rules is searched to determine the center coordinates that should be presented in the visualization. For example, if the spatial position mapping rule is established based on a certain functional relationship, the geographic coordinate value is substituted into the function, and after a series of operations (similar to first performing a certain transformation on the longitude, and then performing another operation in combination with the latitude, the specific operation depends on the rules), the corresponding center coordinate value of the graphic is obtained. This process requires dimensional consistency. If the geographic coordinates themselves have been dimensionally unified, there is no need to consider the dimension again when calculating the center coordinate of the graphic.

[0102] At the same time, according to the attribute value mapping rules, the corresponding color features are found for the damage level information in the disaster record unit. The attribute value mapping rules specify the correspondence between different damage levels and colors. For example, a specific damage level A corresponds to a specific color in the red series. The color features are obtained by searching this correspondence. Similarly, for the emergency response initiation time, the corresponding border features are determined according to the attribute value mapping rules. For example, different emergency response initiation time ranges correspond to different border line types, such as a solid line border in one interval and a dashed line border in another interval.

[0103] The obtained graphic center coordinates, color features, and border features are integrated together to form the basic information of the visualization element. Each basic information set of visualization elements corresponds to a disaster record unit and contains the key visual feature information required for the visualization of the unit.

[0104] Step S142: Construct a visualization element layout module, which is used to receive basic visualization element information, and perform secondary fine-tuning on the center coordinates of adjacent visualization elements based on the density balance adjustment result in the spatial position mapping rule to maintain a recognizable minimum spacing between elements, thereby generating visualization element position information after layout optimization.

[0105] The visualization element layout module obtains basic visualization element information provided by the visualization element generation module. This basic information includes the center coordinates of each visualization element. However, in the initial state, due to the distribution characteristics of different disaster record units, the visualization elements may be unevenly distributed in space. Some areas may be too densely packed with elements, which may affect the visualization effect and information recognition.

[0106] At this point, the result after density balancing adjustment in the reference spatial position mapping rule is used. This result records the relevant information of the adjustment of the graphic distribution density in different areas, including which areas' graphic center coordinates have been adjusted, as well as the direction and magnitude of the adjustment. For adjacent visual elements, secondary fine-tuning is performed based on this adjustment information. For example, if a certain area is identified as having too high density during density balancing adjustment, and the adjustment direction is to move the graphic center coordinates a certain distance to the right, then for the adjacent visual elements in the area, their center coordinates are fine-tuned accordingly based on their relative positional relationship.

[0107] During the fine-tuning process, ensure that the minimum recognizable spacing is maintained between elements. The determination of this minimum spacing is based on a comprehensive consideration of multiple factors, including the resolution of the visualization interface, the characteristics of the display device, and the visual recognition ability of the human eye. For example, through testing of visual effects at different resolutions and analysis of a large amount of user visual recognition feedback data, it is determined that in the current visualization scenario, elements must maintain at least a certain pixel distance in order to be clearly identified. By fine-tuning the center coordinates of adjacent visualization elements, each visualization element can be reasonably distributed in space to avoid mutual occlusion or being too close to be clearly distinguishable. Finally, the visualization element position information after layout optimization is generated. This visualization element position information accurately records the exact position of each visualization element after layout optimization.

[0108] Step S143: Construct a visualization element annotation module, which is used to receive the basic information of the visualization element and the location information after layout optimization, extract the location characteristics of the earthquake source information record and the descriptive characteristics of the damage degree record according to the attribute value mapping rules, and generate associated annotation information including the location annotation and the descriptive annotation.

[0109] Step S143-1: Determine the correspondence between the attribute value and the annotation content in the attribute value mapping rule, the correspondence includes the correspondence between the characteristic parameters of the earthquake source information record and the annotation title content, the correspondence between the characteristic parameters of the damage degree record and the annotation text content, and the correspondence between the characteristic parameters of the emergency response record and the annotation time content.

[0110] By parsing the attribute value mapping rule document, the corresponding regulations for various characteristic parameters and annotation content for earthquake source information records, damage extent records, and emergency response records are determined. For earthquake source information records, the corresponding annotation title content for characteristic parameters such as earthquake focal depth and magnitude is clearly defined. For example, the rule may stipulate that earthquake magnitude within a certain range corresponds to "Strong Earthquake" as the annotation title, and different earthquake focal depth ranges correspond to different depth descriptions as supplementary title information. For damage extent records, the corresponding relationship between characteristic parameters such as building damage type and damage percentage and the annotation content is determined. For example, the damage type of a certain building type corresponds to a specific damage description as the main content, and the damage percentage corresponds to a corresponding quantitative description added to the main text. For emergency response records, the corresponding relationship between characteristic parameters such as the emergency response initiation time and the annotation time content is determined. For example, the initiation time is accurate to the minute and is used as the annotation time content according to a set format (such as "Emergency Response Initiated at XX hours and XX minutes").

[0111] Step S143 - 2 : extracting earthquake source characteristic parameters, damage degree characteristic parameters, and emergency response characteristic parameters from the disaster record unit corresponding to the visualization element as attribute values ​​to be labeled.

[0112] Based on the association between visualization elements and disaster record units, the specific disaster record unit corresponding to each visualization element is found. From this disaster record unit, characteristic parameters of the earthquake source information record, such as magnitude and focal depth, are extracted; characteristic parameters of the damage extent record, such as the type of damaged buildings and the damage ratio, are extracted; and characteristic parameters of the emergency response record, such as the emergency response initiation time, are extracted. These extracted characteristic parameters are the attribute values ​​to be annotated.

[0113] Step S143-3: Generate title annotation content representing basic information of the earthquake source based on the correspondence between the earthquake source characteristic parameters and the annotation title content.

[0114] Based on the correspondence determined in step S143-1, the extracted earthquake source characteristic parameters are substituted into the title. If the extracted earthquake magnitude is a specific value, the corresponding relationship is searched to determine the corresponding annotation title body content, such as "Moderate Earthquake." Combined with other characteristic parameters such as the earthquake source depth, the title content is further supplemented, such as "Earthquake Source Depth XX kilometers." These contents are combined to form a title annotation content that represents the basic information of the earthquake source, clearly conveying the key information of the earthquake source to the user.

[0115] Step S143 - 4 : Based on the correspondence between the damage degree characteristic parameters and the annotation text content, generate text annotation content indicating the specific situation of the damage degree.

[0116] Similarly, damage severity characteristic parameters are processed based on this correspondence. For building damage types, such as "frame structure damage," combined with damage percentages, such as "30% damage," detailed annotations are generated using the descriptions specified in the correspondence. For example, "30% of frame structures in this area have been damaged, primarily manifesting as wall cracks," this provides a comprehensive and detailed description of the actual extent of damage, allowing users to intuitively understand the disaster situation.

[0117] Step S143 - 5 : Generate time annotation content representing emergency response time information according to the correspondence between the emergency response characteristic parameters and the annotation time content.

[0118] Extracted characteristic parameters, such as the emergency response initiation time, are converted according to corresponding relationships. If the emergency response initiation time is a specific moment, a time annotation is generated according to the specified format, such as "Emergency response initiated at XX month XX day XX hour XX minute." This accurately annotates the emergency response time information and helps users understand the timeliness of the emergency response.

[0119] Step S143 - 6 : Determine the size of the graphic body of the visualization element, where the size of the graphic body is determined by the display accuracy and information carrying requirements of the visualization interface.

[0120] Consider the resolution and display requirements of the visualization interface, as well as the richness of information required within the visualization element. If the visualization interface has a high resolution and you wish to display more detailed information within the main body of the graphic, the main body size can be appropriately increased. If the resolution is limited or the information is concise, the main body size can be reduced accordingly. For example, through analysis of the visualization interface and user needs research, it was determined that in the current scenario, to clearly display information and ensure a coherent overall layout, the main body size of the visualization element should be set to X pixels wide and Y pixels high.

[0121] Step S143-7: Calculate the size parameters of the annotation content, where the size parameters of the annotation content are determined by the font characteristics and content length of the annotation text.

[0122] Analyze the length of the title, body, and time annotations, taking into account the selected font characteristics, such as font size and font style. Different font sizes correspond to different character widths and heights. For example, at a certain font size, each character is a pixels wide and b pixels high. Based on the number of characters in the annotation, calculate the width of the title annotation as the number of characters in the title × a, and the height as b. Similarly, calculate the size parameters for the body and time annotations. Combine these size parameters to determine the overall size range of the annotation content.

[0123] Step S143 - 8 : A marking area is reserved in a preset direction of the graphic body, and the size of the marking area is determined by the size parameters of the marking content and the size of the graphic body.

[0124] Based on visualization design specifications and user reading habits, reserve a callout area in a preset direction within the main graphic, such as below or to the right of the main graphic. Determine the specific size of the callout area based on the dimensions of the callout content and the main graphic size. If the overall callout width is W1 and height is H1, and the main graphic width is W and height is H, to ensure the callout information is clear and does not affect the main graphic display, reserve a callout area below the main graphic with a width of W and a height of H1+ΔH (ΔH is the appropriate spacing) to ensure that the callout content can be completely and reasonably placed within this area.

[0125] Step S143-9: Arrange the title annotation content, text annotation content and time annotation content in the annotation area in a preset order, generate associated annotation information containing title, text and time information, and calibrate the position of the associated annotation information to generate final associated annotation information.

[0126] Arrange the content in the annotation area in a predetermined order, such as first the title, then the body, and finally the time. When arranging content, ensure that appropriate spacing is maintained between each section to enhance readability. For example, a certain distance between the title and the body, and a corresponding distance between the body and the time.

[0127] After the arrangement is completed, the position of the entire associated annotation information is calibrated. Check whether the position of the annotation information in the annotation area is centered or meets the set alignment requirements. If the annotation content is not centered in the horizontal direction, calculate the offset that needs to be adjusted based on the width of the annotation area and the total width of the annotation content, and move the entire annotation information to one side by the corresponding pixels to center it in the horizontal direction. Similarly, perform a similar calibration operation in the vertical direction to ensure that the position of the annotation information in the annotation area is accurate and the layout is beautiful. After position calibration, the final associated annotation information is generated. These associated annotation information can clearly and accurately provide key information annotations such as relevant earthquake sources, degree of damage and emergency response time for visualization elements.

[0128] Step S144: Construct a visualization rendering module, which is used to receive the position information and associated annotation information after layout optimization, convert the color features into a color representation recognizable by the visualization interface, convert the border features into a border representation recognizable by the visualization interface, convert the center coordinates into interface pixel coordinates, and generate a visualization intermediate result including the graphic body, border and associated annotations.

[0129] The visualization rendering module obtains the location information of the visualization elements and the associated annotation information after layout optimization. For the color features in the basic information of the visualization elements, since different visualization interfaces may use different color representations, such as RGB color mode, CMYK color mode, etc., the color features need to be converted into the color representation used by the current visualization interface. For example, if the color feature was originally represented by an abstract name or code (such as the specific damage level defined in the attribute value mapping rule corresponds to the color code A), it is necessary to convert the code into a specific RGB value (such as R=255, G=0, B=0 represents red) according to the requirements of the visualization interface. This conversion process is performed based on the color mapping table or conversion algorithm of the visualization interface to ensure that colors are accurately converted between different representations.

[0130] Border features also need to be converted into a recognizable border representation within the visualization interface. If a border feature is initially described as a textual border line type (for example, defining a dashed border corresponding to a certain interval for the emergency response activation time defined in the attribute value mapping rule), it is necessary to convert it into a drawable border style based on the visualization interface's drawing rules. This may involve setting parameters such as the border line width and dashed line interval to accurately render the desired border effect within the visualization interface.

[0131] The center coordinates need to be converted to interface pixel coordinates. The visual interface displays in pixels, and the center coordinates obtained previously may be based on a geographic coordinate system or other relative coordinate system. A coordinate conversion algorithm is used to convert the center coordinates from the original coordinate system to the pixel coordinate system of the visual interface. For example, based on the size and scale of the visual interface, as well as the mapping between the original coordinate system and the pixel coordinate system, the specific position of the center coordinate in the pixel coordinate system is calculated (e.g., the original center coordinate X is multiplied by the scale and added with the offset to obtain the pixel coordinate X_pixel).

[0132] The converted color representation, border representation, and pixel coordinates are integrated together, along with associated annotation information, to generate an intermediate visualization result that includes the graphic body (filled with the color representation), the border (drawn with a specific border representation), and the associated annotations. This intermediate visualization result already has the basic elements for preliminary display on a visualization interface.

[0133] Step S145: Construct a visualization correction module, which is used to receive the intermediate visualization results, check the position overlap between the graphic body and the associated annotations, adjust the position of the overlapping annotation information, and adjust the direction according to the spatial distribution of the overlapping area to generate the final disaster spatial visualization display model.

[0134] The visualization correction module obtains the intermediate visualization results generated by the visualization rendering module. In this intermediate result, due to some slight deviations between the main graphic and the associated annotations during the layout and conversion process, the main graphic and the associated annotations may overlap, affecting the clear display of information.

[0135] Therefore, we examine the intermediate visualization results and determine the bounding ranges of the main graphic and associated annotations using conventional spatial position analysis methods. For example, we calculate the bounding range of the main graphic based on its center coordinates and dimensions (such as radius and side length). For associated annotations, we determine their bounding ranges based on their starting coordinates and content size. We then compare the bounding ranges of the main graphic and associated annotations to determine if there is any overlap.

[0136] If overlap is detected, the direction of adjustment is determined based on the spatial distribution of the overlapping area. For example, if the overlapping area is primarily to the right of the main graphic, the associated annotation information will be repositioned to the left. The extent of the adjustment will be determined by the size of the overlapping area and the layout requirements of the visualization interface. During the adjustment process, ensure that the associated annotation information remains logically associated with the corresponding main graphic and does not overlap with other visual elements.

[0137] For adjacent regions of multiple associated annotations, the spacing between them is also checked. If the spacing is found to be less than the minimum recognizable spacing, the horizontal or vertical position of the adjacent annotations is adjusted based on their relative positions. For example, if the vertical spacing between two adjacent annotations is too small, one of the annotations is moved vertically up or down a certain distance to ensure sufficient recognizable spacing between them.

[0138] After processing the overlap between the main graphic and the associated annotations and adjusting the spacing between adjacent annotations, the final disaster spatial visualization model was generated. This model ensures the accuracy and clarity of the visualization and effectively presents earthquake disaster information to users in an intuitive manner.

[0139] Step S150: Outputting a spatial visualization display result of the earthquake disaster through the disaster spatial visualization display model, wherein the spatial visualization display result includes disaster distribution graphics at different geographic coordinate positions and associated annotation information of corresponding attribute values.

[0140] Step S151: selecting disaster record units within a target display range from the earthquake disaster database, where the target display range is determined by geographic area boundary information input by a user.

[0141] The user inputs geographic area boundary information through the interactive interface, such as defining a rectangular area in the form of longitude and latitude ranges, or defining an irregular area in the form of polygon vertex coordinates. The system filters the earthquake disaster database based on the boundary information input by the user. For each disaster record unit in the database, its geographic coordinate information is extracted. The geographic coordinates of the disaster record unit are compared with the geographic area boundary input by the user through the spatial position judgment algorithm. For example, if the user inputs the longitude and latitude boundaries of a rectangular area, it is determined whether the longitude and latitude of the disaster record unit are within the longitude and latitude range of the rectangle; if the input is the vertex coordinates of a polygon, the point-to-polygon relationship judgment algorithm (such as the ray method, etc.) is used to determine whether the coordinate point of the disaster record unit is within the polygon area. Only disaster record units that meet the user-defined geographic area boundary conditions are selected as disaster record units in the target display range.

[0142] Step S152: inputting the disaster record units in the target display range into the visualization element generation module of the disaster space visualization display model to generate basic visualization element information corresponding to the geographic coordinate position.

[0143] The selected disaster record units within the target display range are sequentially input into the visualization element generation module of the disaster spatial visualization model. This module processes each disaster record unit according to the spatial location mapping rules and attribute value mapping rules in the visualization rule set. For each disaster record unit's geographic coordinate location, its center coordinates within the visualization graph are determined according to the spatial location mapping rules. For example, this can be achieved by searching a pre-established mapping table or executing conventional coordinate conversion algorithms (e.g., calculating the proportional relationship between geographic coordinates and visualization interface coordinates, offsets, etc.) to obtain the corresponding center coordinates.

[0144] At the same time, according to the attribute value mapping rules, information such as the degree of damage and the emergency response initiation time is extracted from the disaster record unit, and the corresponding color features and border features are determined. For example, if the degree of damage in the disaster record unit is a specific level, the color corresponding to that level is searched in the attribute value mapping rules to obtain the color features. For the emergency response initiation time, the corresponding border features are determined according to the same rules. The obtained graphic center coordinates, color features, and border features are combined to generate the basic information of the visualization element corresponding to the geographic coordinate location. Each disaster record unit generates a corresponding set of basic visualization element information.

[0145] Step S153: Optimize the layout of the basic information of the visualization elements through the visualization element layout module of the disaster space visualization display model, calculate the center coordinate spacing of adjacent visualization elements, adjust the direction of the element pairs whose center coordinate spacing is less than the minimum recognizable spacing, and generate layout-optimized position information.

[0146] The visual element layout module receives the basic visual element information output by the visual element generation module. For each visual element in this basic information, it calculates the center coordinate spacing between adjacent visual elements. It determines the distance between adjacent elements using a coordinate distance calculation formula (such as the Euclidean distance formula).

[0147] The calculated center coordinate spacing is compared with the minimum recognizable spacing. The minimum recognizable spacing is a fixed value pre-set based on the display characteristics of the visualization interface and user visual recognition requirements. For element pairs whose center coordinate spacing is less than the minimum recognizable spacing, a directional adjustment strategy is determined based on the spatial distribution of the element area. For example, if multiple elements in a certain area are closely spaced and exhibit a clustering trend, the direction of this clustering trend is analyzed (such as horizontal or vertical density), and then the pairs of elements with closely spaced spacing are adjusted in a direction perpendicular to the clustering trend. The magnitude of the adjustment is determined based on the difference between the spacing and the minimum spacing and the distribution of surrounding elements to ensure that the spacing between elements after adjustment exceeds the minimum recognizable spacing while maintaining the overall layout structure. After checking and adjusting the spacing of all adjacent elements, optimized layout position information is generated, which accurately records the exact position of each visual element after the optimized layout.

[0148] Step S154: adding associated annotation information to the visualization elements after the layout optimization through the visualization element annotation module of the disaster space visualization display model, wherein the associated annotation information includes an earthquake source location annotation and a damage degree description annotation.

[0149] The visualization element annotation module obtains the optimized visualization element information and the corresponding disaster record unit information from the earthquake disaster database. Based on attribute value mapping rules, it extracts location features from the earthquake source information records of the disaster record units to generate the earthquake source location annotation. For example, the longitude and latitude of the earthquake source may be obtained from the earthquake source information records and converted (e.g., to a coordinate representation that is easily understood on the visualization interface or scaled according to the visualization interface) to form the earthquake source location annotation content.

[0150] Similarly, descriptive features are extracted from the damage extent records within the disaster record unit to generate damage extent annotations. For example, damage extent records contain information such as the type of building damage and the percentage of damage. Based on attribute value mapping rules, this information is converted into concise and clear text descriptions, such as "X% of the building's walls are cracked, primarily in brick-concrete structures," which serve as the damage extent annotation content.

[0151] The generated earthquake source location and damage description labels are added to the corresponding optimized visualization elements. The label placement is determined based on the location of the visualization element and the layout specifications of the visualization interface. For example, labels are typically placed near the visualization element and do not overlap with other elements or labels to ensure clear display of information. By adding this associated label information to the visualization elements, the visualization can more comprehensively and accurately convey relevant information about the earthquake disaster.

[0152] Step S155: The visualization rendering module of the disaster spatial visualization display model converts the layout optimized location information and associated annotation information into interface pixel coordinates, draws the fill color corresponding to the color feature and the border line corresponding to the border feature, superimposes the associated annotation information, and generates an earthquake disaster spatial visualization display result including spatial distribution graphics and attribute annotation information.

[0153] The visualization rendering module obtains the optimized layout location information and associated annotation information. First, the center coordinates in the optimized layout visualization element location information are converted into interface pixel coordinates. This conversion process is performed based on the mapping relationship between the visualization interface coordinate system and the previously used geographic coordinate system or other relative coordinate system. For example, if the size of the visualization interface and its proportional relationship with the geographic coordinate system are known, the center coordinates can be converted to interface pixel coordinates by multiplying the center coordinates by the scale and adding an offset (the specific operation is determined by the actual mapping relationship).

[0154] For the color features of visual elements, the color features are converted to corresponding fill colors based on the color representation supported by the visualization interface. For example, if the visualization interface uses the RGB color model and the color feature is originally represented by a code or name, it is converted to a specific RGB value by looking up the color map, such as (255, 165, 0) representing orange. This RGB value is then used to draw the fill color on the visualization interface.

[0155] Similarly, for border features, they are converted into border lines according to the visualization interface's drawing rules. For example, if the border feature is a dashed border, the visualization interface's drawing function and parameter settings are used to determine parameters such as the dashed line spacing and line width, and then the corresponding border line is drawn.

[0156] Finally, the associated annotation information is overlaid onto the drawn graphic body and border. The display format and position of the annotation information are adjusted according to the design specifications of the visualization interface to ensure that the annotation information is clear and readable and does not obscure the graphic body and border. Through these operations, a spatial visualization of the earthquake disaster situation is generated, including a spatial distribution graphic (composed of fill color and border line) and attribute annotation information. This intuitively presents the distribution of earthquake disasters at different geographic coordinate locations and related attribute information.

[0157] Step S156: The visualization correction module of the disaster spatial visualization display model performs a position overlap check on the display result, and adjusts the annotation position of the area where the annotation information overlaps with the graphic body. The adjustment direction is determined by the spatial distribution of the overlapping area.

[0158] Step S156 - 1 : Obtaining the graphic body position information and the annotation information position information in the visual display result, wherein the graphic body position information includes center coordinates and size parameters, and the annotation information position information includes start coordinates and content size parameters.

[0159] Extract the position-related information of the graphic body and annotation information from the data structure of the visualization result. For the graphic body, extract its center coordinate value, which determines the core position of the graphic in the visualization interface. At the same time, obtain the size parameters of the graphic body, such as the radius of circular graphics and the length and width of rectangular graphics. These size parameters define the boundary range of the graphic body. For the annotation information, extract its starting coordinates, that is, the coordinates of the position where the annotation content begins to be displayed on the visualization interface, and the size parameters of the annotation content, such as the length and height of the annotation text, to determine the spatial range occupied by the annotation information on the visualization interface.

[0160] Step S156 - 2 : Calculate the boundary area of ​​the graphic body according to the graphic body position information, and calculate the boundary area of ​​the annotation information according to the starting coordinates and content size parameters.

[0161] For the main body of the graphic, if it is a circle, take the center coordinate as the center of the circle, and calculate its boundary area according to the radius, that is, determine a circular area with the center coordinate as the center and a radius of R as the boundary of the main body of the graphic; if it is a rectangle, calculate the coordinates of the four vertices of the rectangle according to the center coordinates, length and width, so as to determine the boundary area of ​​the rectangle. For the annotation information, its boundary is determined according to the starting coordinates and content size parameters. For example, if the starting coordinates are (X0, Y0), the width of the annotation content is W, and the height is H, then the coordinates of the four boundary points of the annotation information can be calculated to determine the boundary area of ​​the annotation information. By calculating these boundary areas, the spatial range of the main body of the graphic and the annotation information on the visual interface can be clearly defined, which facilitates the subsequent comparison of whether the two overlap.

[0162] Step S156-3: Determine whether the projections of the boundary area of ​​the graphic body and the boundary area of ​​the annotation information in the spatial dimension intersect to determine whether there is overlap. If there is overlap, calculate the intersection range of the projections of the boundary area of ​​the graphic body and the boundary area of ​​the annotation information to obtain the spatial overlap amount of the overlapping area.

[0163] Project the boundary area of ​​the graphic body and the boundary area of ​​the annotation information onto a two-dimensional plane (i.e., the plane where the visualization interface is located). By comparing the positional relationship between the two projected areas, determine whether there is an intersection. For example, for two rectangular areas, determine whether there is an overlapping part by comparing their coordinate ranges in the X-axis and Y-axis directions. If there is an overlap, further calculate the range of the overlapping area. For rectangular areas, by determining the starting and ending coordinates of the overlapping part in the X-axis and Y-axis directions, calculate the area of ​​the overlapping area and use this as the spatial overlap amount; for irregular combinations such as circles and rectangles, determine the range of the overlapping area and quantify it as the spatial overlap amount through conventional geometric calculation methods in related technologies.

[0164] Step S156 - 4 : adjusting the starting coordinates of the annotation information in the opposite direction of the overlapping direction according to the spatial overlap amount, so that the boundary area of ​​the annotation information has no projection intersection with the boundary area of ​​the graphic body.

[0165] Based on the calculated spatial overlap and overlap direction, the starting coordinates of the annotation information are adjusted. If the overlapping area is mainly above the main body of the graphic, that is, the overlapping direction is downward, then the starting coordinates of the annotation information are moved upward by a certain distance in the Y-axis direction. The moving distance is determined according to the spatial overlap and the margin of the visual interface layout to ensure that the boundary area of ​​the annotation information and the boundary area of ​​the main body of the graphic no longer intersect in the projection after the movement. Similarly, if the overlapping direction is other directions, such as left, right or oblique, adjust the X-axis or Y-axis component of the starting coordinate in the corresponding opposite direction to avoid overlap between the annotation information and the main body of the graphic.

[0166] Step S156 - 5 : recalculating the boundary area of ​​the adjusted annotation information, and checking the overlap with the boundary area of ​​the graphic body again until no overlap exists.

[0167] After adjusting the starting coordinates of the annotation information, the boundary area of ​​the annotation information is recalculated based on the adjusted starting coordinates and the annotation content size parameters. This boundary area is then compared again with the boundary area of ​​the main body of the graphic to check whether there is still overlap. If there is still overlap, steps S156-3 and S156-4 are repeated, i.e., the spatial overlap amount is recalculated and the starting coordinates of the annotation information are adjusted again based on the new overlap direction until the projections of the boundary area of ​​the annotation information and the boundary area of ​​the main body of the graphic in the spatial dimension no longer intersect, ensuring a clear spatial distinction between the annotation information and the main body of the graphic.

[0168] Step S156-6: For areas where multiple annotation information are adjacent, check the spatial spacing between the annotation information, and adjust the position of the annotation information whose spatial spacing is less than the minimum recognizable spacing in the horizontal or vertical direction, thereby generating a corrected earthquake disaster spatial visualization display result, which includes a non-overlapping graphic body and annotation information layout.

[0169] After completing the overlap check and adjustment of the annotation information and the main body of the graphic, the adjacent areas of multiple annotation information are checked. The spatial spacing between adjacent annotation information can be determined by calculating the shortest distance between the boundary points of adjacent annotation information. This spacing is compared with the pre-set minimum recognizable spacing. If it is found that the spacing between two adjacent annotation information is less than the minimum recognizable spacing, the position is adjusted according to the arrangement direction of the annotation information. If the annotation information is arranged horizontally, one of the annotation information is moved a certain distance to the left or right in the horizontal direction so that the spacing between them reaches or exceeds the minimum recognizable spacing; if the annotation information is arranged vertically, it is moved up or down in the vertical direction. By checking and adjusting the spacing of all adjacent annotation information, a corrected earthquake disaster spatial visualization display result is generated. This result ensures that there is a clear spatial layout between the main body of the graphic and the annotation information, as well as between the annotation information, without overlap or spacing that is too small to be distinguished, providing users with a clear and accurate earthquake disaster visualization information display.

[0170] Figure 2 A schematic diagram illustrates exemplary hardware and software components of a spatial visualization display system 100 incorporating an earthquake disaster database, provided in some embodiments of the present application, that can implement the concepts of the present application. For example, a processor 120 can be used in the spatial visualization display system 100 incorporating an earthquake disaster database and to perform the functions described in the present application.

[0171] The spatial visualization display system 100 for an earthquake disaster database can be a general-purpose server or a special-purpose server, both of which can be used to implement the spatial visualization display method for an earthquake disaster database in this application. Although only one server is shown in this application, for convenience, the functions described in this application can be implemented in a distributed manner on multiple similar platforms to balance the processing load.

[0172] For example, the spatial visualization display system 100 in combination with the earthquake disaster database may include a network port 110 connected to the network, one or more processors 120 for executing program instructions, a communication bus 130, and storage media 140 in different forms, such as a disk, ROM, or RAM, or any combination thereof. Exemplarily, the spatial visualization display system 100 in combination with the earthquake disaster database may also include program instructions stored in ROM, RAM, or other types of non-transitory storage media, or any combination thereof. The method of the present application can be implemented according to these program instructions. The spatial visualization display system 100 in combination with the earthquake disaster database also includes an I / O interface 150 between the computer and other input and output devices.

[0173] For ease of explanation, only one processor is described in the spatial visualization display system 100 combined with the earthquake disaster database. However, it should be noted that the spatial visualization display system 100 combined with the earthquake disaster database in the present application may also include multiple processors, so the steps performed by one processor described in the present application may also be performed jointly or individually by multiple processors. For example, if the processor of the spatial visualization display system 100 combined with the earthquake disaster database executes step A and step B, it should be understood that step A and step B may also be performed jointly by two different processors or individually in one processor. For example, the first processor executes step A, the second processor executes step B, or the first processor and the second processor execute steps A and B together.

[0174] In addition, an embodiment of the present invention further provides a readable storage medium, in which computer-executable instructions are preset. When a processor executes the computer-executable instructions, the spatial visualization display method combined with the earthquake disaster database as described above is implemented.

[0175] It should be noted that in order to simplify the description of the present invention and thus help understand one or more embodiments of the invention, in the foregoing description of the embodiments of the present invention, multiple features are sometimes combined into one embodiment, figure or description thereof.

Claims

1. A spatial visualization display method combined with an earthquake disaster database, characterized in that: The method comprises: Acquire an earthquake disaster database, wherein the earthquake disaster database includes multiple types of disaster record units with geographic coordinate marks, and the disaster record units include earthquake source information records, damage extent records, and emergency response records; Performing disaster feature extraction processing on the earthquake disaster database to obtain disaster distribution association features and attribute association features of the disaster record units, wherein the disaster distribution association features represent the correlation between disaster distributions at different geographic coordinate locations, and the attribute association features represent the attribute correlation between earthquake source information, damage degree, and emergency response; Calling a pre-built visualization mapping model to perform visualization rule mapping processing on the disaster distribution association features and the attribute association features, and generating a visualization rule set including spatial position mapping rules and attribute value mapping rules; Constructing a disaster situation spatial visualization display model based on the visualization rule set, wherein the disaster situation spatial visualization display model is used to convert disaster situation record units in the earthquake disaster situation database into visualization elements with spatial location identifiers and attribute value identifiers; The spatial visualization display model for disaster conditions is used to output spatial visualization display results of earthquake disaster conditions, wherein the spatial visualization display results include disaster distribution graphics at different geographic coordinate positions and associated annotation information of corresponding attribute values.

2. The spatial visualization display method combined with the earthquake disaster database according to claim 1 is characterized in that: The extracting of disaster characteristics from the earthquake disaster database to obtain the disaster distribution correlation characteristics and attribute correlation characteristics of the disaster record unit includes: Performing geographic coordinate standardization processing on the disaster record units in the earthquake disaster database to unify the geographic coordinate representation formats of different disaster record units to obtain standardized disaster record units with a unified coordinate reference, wherein the geographic coordinate representation formats include latitude and longitude formats and plane rectangular coordinate formats; Performing spatial neighborhood analysis on the standardized disaster record units, calculating the overlapping relationship between the disaster record coverage area of ​​each geographic coordinate position and the adjacent geographic coordinate position, and generating a neighborhood overlap feature representing the spatial adjacency strength; Performing time series alignment processing on the standardized disaster record unit, mapping the occurrence time of the earthquake source information record, the assessment time of the damage extent record, and the start time of the emergency response record to the same time axis, and generating a time series aligned record unit; Performing attribute co-occurrence analysis on the time-series aligned record units, counting the co-occurrence frequencies of characteristic parameters of earthquake source information records and characteristic parameters of damage degree records, and generating co-occurrence frequency features representing attribute association strength; The neighborhood overlapping features and the co-occurrence frequency features are subjected to feature weighted fusion processing to obtain the disaster distribution association features and attribute association features of the disaster record unit, wherein the disaster distribution association features include the weighted results of the neighborhood overlapping features, and the attribute association features include the weighted results of the co-occurrence frequency features, wherein the weight of the neighborhood overlapping features is determined by the distribution density of the disaster records at the geographic coordinate location, and the weight of the co-occurrence frequency features is determined by the statistical significance of the attribute co-occurrence frequency.

3. The spatial visualization display method combined with the earthquake disaster database according to claim 2 is characterized in that: The process of performing geographic coordinate standardization on the disaster record units in the earthquake disaster database to unify the geographic coordinate representation formats of different disaster record units to obtain standardized disaster record units with a unified coordinate reference includes: Identify the geographic coordinate representation format identifier of each disaster record unit in the earthquake disaster database, wherein the representation format identifier includes a latitude and longitude identifier and a plane rectangular coordinate identifier; For the disaster record unit marked as longitude and latitude format, extract its latitude and longitude values, and call the geographic coordinate conversion function to convert the longitude and latitude values ​​into coordinate values ​​of the plane rectangular coordinate system; For disaster record units marked as plane rectangular coordinate format, extract their original coordinate values, check whether the coordinate units are unified units, and perform unit conversion on coordinate values ​​that are not unified units; The converted plane rectangular coordinates and the original plane rectangular coordinates are processed into a coordinate system, the center position of the earthquake disaster database coverage area is selected as the reference origin, and a local plane rectangular coordinate system is constructed; Performing a translation transformation on the plane rectangular coordinates of all disaster record units relative to the reference origin, adjusting the offset of the coordinate values ​​so that the coordinate value of the reference origin is the reference point; The coordinate values ​​after translation transformation are processed with unified precision, and the coordinate values ​​are retained to a unified number of digits of precision to generate a standardized disaster record unit with a unified coordinate reference, unified unit and unified precision.

4. The spatial visualization display method combined with the earthquake disaster database according to claim 1 is characterized in that: The calling of the pre-built visualization mapping model performs visualization rule mapping processing on the disaster distribution association features and the attribute association features to generate a visualization rule set including spatial location mapping rules and attribute value mapping rules, including: Inputting the disaster distribution correlation features into the spatial mapping submodule of the visualization mapping model, the spatial mapping submodule generates a spatial position mapping rule representing the correspondence between the geographic coordinate position and the center coordinate of the visualization graphic by learning the historical mapping relationship between the geographic coordinate position and the visualization interface layout; Inputting the attribute association features into the attribute mapping submodule of the visualization mapping model, the attribute mapping submodule generates attribute value mapping rules representing the correspondence between the damage degree level and the graphic color feature, and the correspondence between the emergency response initiation time and the graphic border feature by learning the historical association relationship between the attribute value and the visualization graphic style; Performing density balancing adjustment processing on the spatial position mapping rule, calculating the graphic distribution density of each area in the visualization interface, and adjusting the direction of the graphic center coordinates of the area where the density exceeds the preset balancing range, where the adjustment direction is determined by the regional density distribution trend; Performing conflict resolution on the attribute value mapping rules to identify attribute value combinations with overlapping color features or confusing border features, and reallocating styles by adjusting color hue parameters or border line type parameters to make visualization styles of different attribute values ​​distinguishable; The adjusted spatial position mapping rules and attribute value mapping rules are verified for consistency, and the correspondence between the geographic coordinate position and the graphic center coordinate is checked to see whether it covers all disaster record units. The correspondence between the attribute value and the graphic style is verified to see whether it fully maps all attribute categories, and the final visualization rule set is generated.

5. The spatial visualization display method combined with the earthquake disaster database according to claim 4 is characterized in that: The density balancing adjustment process is performed on the spatial position mapping rule, the graphic distribution density of each area in the visualization interface is calculated, and the direction of the graphic center coordinates of the area where the density exceeds the preset balancing range is adjusted, including: The visualization interface is divided into grid-like density calculation areas. The size of each density calculation area is determined by the overall size of the interface and the density calculation accuracy. Count the number of visual graphics in each density calculation area, and calculate the ratio of the number of graphics to the area of ​​the area as the regional density parameter; identifying a first density region having a regional density parameter outside a density equilibrium range, and a second density region outside the first density region; Sorting the visualized graphics within the first density region to obtain sorted graphics, wherein the sorting rule is determined by the spatial distribution characteristics of the center coordinates of the graphics; The sorted graphics are adjusted in direction in sequence. The adjustment direction is determined by the spatial distribution trend of the first density area, and the adjustment amount is determined by the difference between the regional density parameter and the equilibrium range. Supplementing the position of the visual graphics within the second density area, where the supplemented position is determined by the spatial distribution characteristics of the second density area; Recalculating the density parameters of the regions after direction adjustment and position supplementation, and repeating the above adjustment process for regions still outside the density balance range until the density parameters of all regions are within the density balance range; The final adjusted positions of all graphics are recorded, and an adjusted spatial position mapping rule is generated, wherein the spatial position mapping rule includes a correspondence between the original center coordinates and the final adjusted coordinates.

6. The spatial visualization display method combined with the earthquake disaster database according to claim 1 is characterized in that: The step of constructing a disaster situation spatial visualization display model based on the visualization rule set includes: Constructing a visualization element generation module, the visualization element generation module is used to receive the disaster record unit in the earthquake disaster database, extract the center coordinates of the graphic corresponding to the geographic coordinate position according to the spatial position mapping rule, extract the color features corresponding to the damage degree level and the border features corresponding to the emergency response initiation time according to the attribute value mapping rule, and generate basic visualization element information including the center coordinates, color features and border features; Constructing a visualization element layout module, the visualization element layout module is used to receive the visualization element basic information, and according to the density balance adjustment result in the spatial position mapping rule, perform a secondary fine-tuning on the center coordinates of adjacent visualization elements so as to maintain a recognizable minimum spacing between elements, thereby generating visualization element position information after layout optimization; Constructing a visualization element annotation module, the visualization element annotation module is used to receive the visualization element basic information and the location information after layout optimization, extract the location characteristics of the earthquake source information record and the descriptive characteristics of the damage extent record according to the attribute value mapping rule, and generate associated annotation information including the location annotation and the descriptive annotation; Constructing a visualization rendering module, which is used to receive the layout-optimized position information and associated annotation information, convert the color features into a color representation recognizable by the visualization interface, convert the border features into a border representation recognizable by the visualization interface, convert the center coordinates into interface pixel coordinates, and generate a visualization intermediate result including the graphic body, border, and associated annotations; A visualization correction module is constructed. The visualization correction module is used to receive the intermediate visualization results, check the position overlap between the graphic body and the associated annotations, adjust the position of the overlapping annotation information, and adjust the direction according to the spatial distribution of the overlapping area to generate the final disaster spatial visualization display model.

7. The spatial visualization display method combined with the earthquake disaster database according to claim 6 is characterized in that: The constructing of a visualization element annotation module, wherein the visualization element annotation module is configured to add associated annotation information to the visualization element according to the attribute value mapping rule, includes: Determining the correspondence between the attribute value and the annotation content in the attribute value mapping rule, the correspondence including the correspondence between the characteristic parameters of the earthquake source information record and the annotation title content, the correspondence between the characteristic parameters of the damage degree record and the annotation text content, and the correspondence between the characteristic parameters of the emergency response record and the annotation time content; Extracting earthquake source characteristic parameters, damage degree characteristic parameters, and emergency response characteristic parameters from the disaster record unit corresponding to the visualization element as attribute values ​​to be labeled; Generate title annotation content representing basic information of the earthquake source according to the corresponding relationship between the earthquake source characteristic parameters and the annotation title content; Generate text annotation content indicating the specific situation of the damage degree according to the corresponding relationship between the damage degree characteristic parameters and the annotation text content; Generate time annotation content representing emergency response time information based on the correspondence between emergency response characteristic parameters and annotation time content; Determine the size of the main graphic body of the visualization element, wherein the main graphic body size is determined by the display accuracy and information carrying requirements of the visualization interface; Calculating size parameters of the annotation content, where the size parameters of the annotation content are determined by font characteristics and content length of the annotation text; A marking area is reserved in the preset direction of the graphic body. The size of the marking area is determined by the size parameters of the marking content and the size of the graphic body. The title annotation content, text annotation content and time annotation content are arranged in a preset order in the annotation area to generate associated annotation information including the title, text and time information, and the position of the associated annotation information is calibrated to generate the final associated annotation information.

8. The spatial visualization display method combined with the earthquake disaster database according to claim 1 is characterized in that: Outputting the spatial visualization display result of the earthquake disaster using the disaster spatial visualization display model includes: Selecting disaster record units within a target display range from the earthquake disaster database, wherein the target display range is determined by geographic area boundary information input by a user; Inputting the disaster record units in the target display range into the visualization element generation module of the disaster space visualization display model to generate basic visualization element information corresponding to the geographic coordinate position; The visualization element layout module of the disaster spatial visualization display model performs layout optimization processing on the basic information of the visualization elements, calculates the center coordinate spacing of adjacent visualization elements, adjusts the orientation of element pairs whose center coordinate spacing is less than the minimum recognizable spacing, and generates layout-optimized position information; Adding associated annotation information to the layout-optimized visualization elements through the visualization element annotation module of the disaster space visualization display model, wherein the associated annotation information includes an earthquake source location annotation and a damage degree description annotation; The visualization rendering module of the disaster spatial visualization display model converts the layout optimized location information and associated annotation information into interface pixel coordinates, draws the fill color corresponding to the color feature and the border line corresponding to the border feature, and superimposes the associated annotation information to generate an earthquake disaster spatial visualization display result including spatial distribution graphics and attribute annotation information; The visualization correction module of the disaster spatial visualization display model is used to check the position overlap of the display results, and the annotation position of the area where the annotation information overlaps with the graphic body is adjusted. The adjustment direction is determined by the spatial distribution of the overlapping area.

9. The spatial visualization display method combined with the earthquake disaster database according to claim 8, characterized in that: The visual correction module of the disaster space visual display model performs a position overlap check on the display result and adjusts the position of the annotation information in the area where the annotation information overlaps with the main body of the graphic, including: Obtaining the graphic body position information and the annotation information position information in the visualization display result, wherein the graphic body position information includes the center coordinates and size parameters, and the annotation information position information includes the starting coordinates and content size parameters; Calculating a boundary area of ​​the graphic body according to the graphic body position information, and calculating a boundary area of ​​the annotation information according to the starting coordinates and content size parameters; Determining whether projections of the boundary area of ​​the graphic body and the boundary area of ​​the annotation information in the spatial dimension intersect to determine whether there is overlap; if there is overlap, calculating the intersection range of the projections of the boundary area of ​​the graphic body and the boundary area of ​​the annotation information to obtain a spatial overlap amount of the overlapping area; Adjusting the starting coordinates of the annotation information in a direction opposite to the overlapping direction according to the spatial overlap amount, so that a boundary area of ​​the annotation information has no projection intersection with a boundary area of ​​the graphic body; Recalculating the boundary area of ​​the adjusted annotation information and checking the overlap with the boundary area of ​​the graphic body again until no overlap exists; For areas where multiple annotation information are adjacent, the spatial spacing between the annotation information is checked. For annotation information whose spatial spacing is less than the minimum recognizable spacing, the position is adjusted horizontally or vertically, thereby generating a corrected earthquake disaster spatial visualization display result, which includes a non-overlapping graphic body and annotation information layout.

10. A spatial visualization display system combined with an earthquake disaster database, characterized in that: It includes a processor and a memory, the memory is connected to the processor, the memory is used to store programs, instructions or codes, and the processor is used to execute the programs, instructions or codes in the memory to implement the spatial visualization display method combined with the earthquake disaster database as described in any one of claims 1 to 9 above.

Images