Freeze-thaw geological disaster meteorological risk early warning method based on temperature and ground temperature changes

Abstract

The invention discloses a freeze-thaw geological disaster meteorological risk early warning method based on temperature and ground temperature changes, and relates to the technical field of computer processing, and the method comprises the following steps: S01, obtaining vector data of a target area; s02, converting the collected vector data into original raster data based on ArcGIS software, and endowing the converted original raster data with corresponding attribute values; s03, endowing corresponding attribute values to the original raster data through a Kriging interpolation method for processing so as to generate soil ground temperature contour surface raster data of different depths; and S04, superposing the soil ground temperature contour surface grid data and the original grid data according to a preset matrix discrimination relation to generate a geological disaster freeze-thaw induction grade grid map and a freeze-thaw geological disaster susceptible grade grid map. Through real-time monitoring and early warning, geological disaster risk information can be accurately provided, and casualties and property loss can be reduced.

Description

Technical Field

[0001] The present invention relates to the field of computer processing technologies, and particularly relates to a meteorological risk warning method for freeze-thaw geological disasters based on air temperature and ground temperature changes. Background Art

[0002] Since China spans a relatively large latitude, many provinces and cities span three climate zones, so the climate changes greatly, the freeze-thaw areas are widely distributed, and the area of high-risk areas for geological disasters caused by freeze-thaw is 16,000 km2, the medium-risk area is 40,000 km2, and the low-risk area is 76,000 km2. During the freeze-thaw cycle, the changes in the physical and mechanical properties of loess are affected by many factors such as soil moisture content, number of freeze-thaw cycles, freezing temperature, and freezing speed, which will cause drastic changes in the physical and mechanical properties of soil samples, thus triggering catastrophic damages caused by freeze-thaw such as frost heaving, thaw settlement, and landslides. According to statistics, from 2012 to 2023, 88 freeze-thaw geological disasters occurred in Shaanxi Province, causing 28 deaths, 9 injuries, and direct economic losses of 24.2339 million yuan.

[0003] Currently, the meteorological risk warning for geological disasters mainly targets geological disasters induced by rainfall. For example, in Chinese Patent, publication number CN107679167A, a meteorological risk assessment method and analysis and warning platform based on gridded meteorological data is disclosed, including: a data acquisition module, a data cleaning module, a data fusion module, and a weather risk assessment module. This meteorological risk assessment method and analysis and warning platform fully cleans the historical meteorological data obtained through multiple channels, corrects the errors in the historical meteorological data to ensure the rigor of the data source of the basic database. After interpolation of the cleaned historical meteorological data, combined with refined topographic data and satellite data inversion verification, gridded meteorological data is formed, and all relevant data is unified into the same physical space, and fused through a spatial field algorithm to obtain the value closest to the actual situation. This meteorological risk assessment method and analysis and warning platform incorporates the existing meteorological forecast data into the meteorological risk assessment model, realizes the real-time broadcast of meteorological disasters, and can accurately assess the risk level of the target meteorological disaster.

[0004] In the prior art including the above patent, it can be seen that the existing warning systems all detect rain and snow weather, and judge geological disasters by comparing the rainfall with the historical same-position or similar rainfall, but lack a warning mechanism for meteorological warnings of geological disasters caused by freeze-thaw cycles. Summary of the Invention

[0005] The purpose of the present invention is to provide a meteorological risk warning method for freeze-thaw geological disasters based on air temperature and ground temperature changes to solve the above problems.

[0006] To achieve the above purpose, the present invention provides the following technical solutions:

[0007] A meteorological risk warning method for freeze-thaw geological disasters based on air temperature and ground temperature changes, comprising the following steps:

[0008] S01. Obtain the vector data of the target area;

[0009] S02. Based on the ArcGIS software, convert the collected vector data into original raster data, and assign corresponding attribute values to the converted original raster data;

[0010] S03. Process the original raster data with corresponding attribute values assigned by the Kriging interpolation method to generate soil ground temperature isosurface raster data at different depths;

[0011] S04. Superimpose the soil ground temperature isosurface raster data and the original raster data according to a preset matrix discrimination relationship to generate a geological disaster freeze-thaw induced grade raster map and a freeze-thaw geological disaster prone grade raster map;

[0012] S05. Superimpose the geological disaster freeze-thaw induced grade raster map and the freeze-thaw geological disaster prone grade raster map according to a preset risk discrimination matrix, and render them into a five-level warning map of red, orange, yellow, blue, and white.

[0013] Preferably, the vector data includes a freeze-thaw zoning map, a geological disaster prone degree zoning map, and a forecast air temperature distribution map for the next 24 hours.

[0014] Preferably, the freeze-thaw zoning map is obtained based on the data collected by a satellite equipped with a thermal infrared sensor;

[0015] The acquisition of the geological disaster prone degree zoning map includes:

[0016] S11. Establish a database according to the disaster types that occurred in different terrains under different freeze-thaw parameters collected historically;

[0017] S12. Obtain the current freeze-thaw zoning map, extract the freeze-thaw parameters corresponding to the terrain, and match them with the established database;

[0018] Preferably, the judgment results include high prone, medium prone, low prone, and non-prone.

[0019] Preferably, in step S01, the grid cell size in the original raster data is 3*3 km, and the coordinate system adopts the GCS_WGS_1984 geographic coordinate system.

[0020] Preferably, in step S02, the processing of the vector data by the ArcGIS software includes:

[0021] S21. Use ArcGIS software to import the collected vector data, which contains the geospatial features to be converted.

[0022] S22. In the ArcToolbox of ArcGIS, select the Raster Conversion sub-menu under the Conversion Tools, and perform polygon to raster according to the predetermined conversion features of the vector data to obtain raster data.

[0023] S23. Use the Raster Calculator tool to assign new attribute values to each pixel of the raster data to obtain the original raster data.

[0024] The original raster data is the topographic feature of the target area.

[0025] Preferably, in step S03, obtaining the raster data of the soil geothermal isotherm surfaces at different depths includes:

[0026] S31. Based on the data of 981 monitoring stations arranged in the target area according to a predetermined specification;

[0027] S32. Clean the obtained data to remove missing values or outliers;

[0028] S33. Construct a Kriging interpolation model based on the Kriging interpolation method to predict the soil temperature values at the monitoring stations;

[0029] S34. According to the temperature value intervals of the 981 monitoring stations, derive the temperature differences between adjacent two monitoring stations, and the temperature decreasing data of the linear distance between adjacent two monitoring stations under the temperature differences.

[0030] Preferably, in step S04, superimposing and calculating the soil geothermal isotherm surface raster data and the original raster data according to a preset matrix discrimination relationship includes:

[0031] S41. Superimpose and calculate the geothermal raster data and the forecast air temperature raster data according to a preset matrix discrimination relationship to generate a raster map of the freezing and thawing induced grade of geological disasters;

[0032] S42. Superimpose and calculate the raster map of the prone degree zoning of geological disasters and the raster map of the freezing and thawing zoning according to a preset risk discrimination matrix relationship to generate a raster map of the prone grade of freezing and thawing geological disasters.

[0033] Preferably, the risk levels of red, orange, yellow, blue, and white in step S05 are sorted from high to low.

[0034] In the above technical solution, a method for meteorological risk warning of freezing and thawing geological disasters based on the changes in air temperature and ground temperature provided by the present invention has the following beneficial effects:

[0035] 1. By integrating multiple data sources, including freeze-thaw zoning maps, geological hazard susceptibility zoning maps, predicted air temperature distribution maps for the next 24 hours, and advanced geographic information technology, a comprehensive and integrated assessment of the meteorological risks of freeze-thaw geological disasters is achieved, enabling more accurate capture of potential geological disaster risks and improving the accuracy and timeliness of early warnings.

[0036] 2. By converting vector data into raster data using ArcGIS software and assigning corresponding attribute values, a refined description of the spatial characteristics of the target area is achieved. The raster cell size is 3*3 km, ensuring high spatial resolution of the early warning results and facilitating the identification of more detailed geological disaster risk areas.

[0037] 3. By processing the original raster data using Kriging interpolation method, raster data of soil geothermal isotherm surfaces at different depths is generated, combining the key factor of soil temperature with depth information, providing more comprehensive data support for assessing geological disaster risks.

[0038] 4. By overlaying and calculating multiple factors through preset matrix discrimination relationships and risk discrimination matrices, a raster map of the freeze-thaw induced grade of geological disasters and a raster map of the susceptibility grade of freeze-thaw geological disasters are finally generated, and further rendered into a five-level early warning map of red, orange, yellow, blue, and white, which helps decision-makers take corresponding disaster prevention and mitigation measures according to the early warning level. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] To more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in the present invention, and those of ordinary skill in the art can also obtain other drawings based on these drawings.

[0040] Figure 1 Flow chart provided for Embodiment 1 of the present invention;

[0041] Figure 2 Flow chart of S01 provided for Embodiment 1 of the present invention;

[0042] Figure 3 Flow chart of S02 provided for Embodiment 1 of the present invention;

[0043] Figure 4 Flow chart of S03 provided for Embodiment 1 of the present invention;

[0044] Figure 5 Schematic diagram of the surface raster data provided for Embodiment 1 of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0045] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0046] Embodiment 1

[0047] As Figures 1-5 shown, a meteorological risk early warning method for freeze-thaw geological disasters based on air temperature and ground temperature changes includes the following steps:

[0048] S01. Obtain the vector data of the target area, where the vector data includes a freeze-thaw zoning map, a geological disaster susceptibility zoning map, and a forecast air temperature distribution map for the next 24 hours;

[0049] S02. Based on the ArcGIS software, convert the collected vector data into original raster data, and assign corresponding attribute values to the converted original raster data;

[0050] S03. Process the original raster data with corresponding attribute values through Kriging interpolation to generate soil ground temperature isosurface raster data at different depths;

[0051] S04. Superimpose the soil ground temperature isosurface raster data and the original raster data according to a preset matrix discrimination relationship to generate a geological disaster freeze-thaw induced grade raster map and a freeze-thaw geological disaster susceptibility grade raster map;

[0052] S05. Superimpose the geological disaster freeze-thaw induced grade raster map and the freeze-thaw geological disaster susceptibility grade raster map according to a preset risk discrimination matrix, and render them into a five-level early warning map of red, orange, yellow, blue, and white.

[0053] Specifically, the grid cell size in the above-mentioned original raster data is 3*3 km, and the coordinate system adopts the GCS_WGS_1984 geographic coordinate system ( Figure 5 ).

[0054] Furthermore, it is displayed by grouping the attribute values of grid cells into different categories, and corresponding colors are assigned to these categories) to generate a five-level early warning map of red, orange, yellow, blue, and white with very high risk (the grid color is rendered red), high risk (the grid color is rendered orange), relatively high risk (the grid color is rendered yellow), certain risk (the grid color is rendered blue), and low risk (the grid color is rendered white).

[0055] By integrating multiple data sources, including freeze-thaw zoning maps, geological hazard susceptibility zoning maps, and predicted temperature distribution maps for the next 24 hours, and advanced geographic information technology, a comprehensive and integrated assessment of the meteorological risks of freeze-thaw geological hazards is achieved, enabling more accurate capture of potential geological hazard risks and improving the accuracy and timeliness of early warnings.

[0056] As a further embodiment provided by the present invention, as Figure 2 shown, the acquisition of the freeze-thaw zoning map is based on the data collected by a satellite equipped with a thermal infrared sensor;

[0057] The acquisition of the geological hazard susceptibility zoning map includes:

[0058] S11. Establish a database according to the historical data of the types of disasters that occurred in different terrains under different freeze-thaw parameters;

[0059] S12. Obtain the current freeze-thaw zoning map, extract the freeze-thaw parameters corresponding to the terrain, and match them with the established database;

[0060] S13. Generate a judgment result for the corresponding terrain in the new mask of the freeze-thaw zoning map according to the matching result. The judgment result includes high susceptibility, medium susceptibility, low susceptibility, and non-susceptibility.

[0061] Specifically, by collecting data through a satellite equipped with a thermal infrared sensor, rapid and accurate monitoring of a large area can be achieved, thereby generating a high-precision freeze-thaw zoning map. This method has higher spatio-temporal resolution and wider coverage compared to traditional ground observations. And based on the statistical analysis of historical data and disaster types, a database is established, and the geological hazard susceptibility zoning map is generated according to the current freeze-thaw parameters.

[0062] As yet another embodiment provided by the present invention, as Figure 3 shown, the processing of vector data by ArcGIS software in step S02 includes:

[0063] S21. Use ArcGIS software to import the collected vector data, which contains geospatial elements to be converted;

[0064] S22. In the ArcToolbox of ArcGIS, select the Rasterize sub-menu under the Conversion Tools, and perform polygon to raster according to the predetermined conversion elements of the vector data to obtain raster data;

[0065] S23. Use the Raster Calculator tool to assign new attribute values to each pixel of the raster data to obtain the original raster data;

[0066] The original raster data is the terrain feature of the target area.

[0067] Specifically, the use of ArcGIS software for data import, conversion, and attribute value assignment has automated the data processing and improved work efficiency. Further processing through professional geographic information system software ensures the accuracy and reliability of the data.

[0068] As yet another embodiment further provided by the present invention, as Figure 4 shown, in step S03, obtaining soil geothermal isothermal surface raster data at different depths includes;

[0069] S31. Based on the data of 981 monitoring stations arranged in the target area according to a predetermined specification;

[0070] S32. Cleaning the obtained data to remove missing values or outliers;

[0071] S33. Constructing a Kriging interpolation model based on the Kriging interpolation method to predict the soil temperature values at the monitoring stations;

[0072] S34. According to the temperature value spacing between 981 monitoring stations, deriving the temperature difference between two adjacent monitoring stations and the temperature decreasing data of the linear spacing between two adjacent monitoring stations under the temperature difference.

[0073] Specifically, arranging 981 monitoring stations in the target area ensures the wide coverage and representativeness of the data. Cleaning and interpolating the obtained data removes missing values and outliers, improving the accuracy and integrity of the data. At the same time, an interpolation model is constructed through the Kriging interpolation method to achieve the prediction of soil temperature and the generation of isothermal surface raster data.

[0074] As yet another embodiment further provided by the present invention, in step S04, superimposing and calculating the soil geothermal isothermal surface raster data and the original raster data according to a preset matrix discrimination relationship includes:

[0075] S41. Superimposing and calculating the geothermal raster data and the predicted air temperature raster data according to a preset matrix discrimination relationship to generate a raster map of the freeze-thaw induced grade of geological disasters;

[0076] S42. Superimposing and calculating the raster map of the geological disaster prone degree zoning and the freeze-thaw zoning raster map according to a preset risk discrimination matrix relationship to generate a raster map of the freeze-thaw geological disaster prone grade.

[0077] Specifically, superimposing and calculating the geothermal raster data and the predicted air temperature raster data, and the raster map of the geological disaster prone degree zoning and the freeze-thaw zoning raster map realizes a comprehensive and integrated assessment of the geological disaster risk. And through the preset matrix discrimination relationship and risk discrimination matrix, a raster map of the freeze-thaw induced grade of geological disasters and a raster map of the freeze-thaw geological disaster prone grade are generated, providing strong support for scientific early warning.

[0078] The above-mentioned matrix discrimination relationship and risk discrimination matrix are respectively the daily maximum air temperature and the daily maximum ground temperature. By setting the range of classification thresholds, for example: the daily maximum air temperature of 0-5 degrees, 5-10 degrees, 10-15 degrees, three thresholds are set respectively; the daily maximum ground temperature of 15-25 degrees, 25-35 degrees, above 35°, three thresholds are set respectively. Or more threshold types can be refined and more detailed classifications can be set, so as to generate a grid map of the freeze-thaw induced grade of geological disasters and a grid map of the prone grade of freeze-thaw geological disasters, and more can be set.

[0079] In the first embodiment, the vector data is converted into raster data through ArcGIS software, and corresponding attribute values are assigned, realizing the refined description of the spatial characteristics of the target area. The raster cell size is 3*3km, ensuring the high resolution of the warning result in space and helping to identify more detailed geological disaster risk areas. Then, the original raster data is processed by the Kriging interpolation method to generate raster data of soil ground temperature isosurfaces at different depths, combining the key factor of soil temperature with depth information, providing more comprehensive data support for evaluating geological disaster risks. Finally, through the preset matrix discrimination relationship and risk discrimination matrix, multiple factors are superimposed and calculated, and finally a grid map of the freeze-thaw induced grade of geological disasters and a grid map of the prone grade of freeze-thaw geological disasters are generated, and further rendered into a five-level warning map of red, orange, yellow, blue, and white, which helps decision-makers take corresponding disaster prevention and mitigation measures according to the warning level.

[0080] Embodiment 2

[0081] The embodiment of the present invention provides a non-transitory computer-readable storage medium, in which at least one instruction or at least one segment of program is stored, and at least one instruction or at least one segment of program is loaded and executed by a processor to implement the steps:

[0082] Obtain the vector data of the target area;

[0083] Based on ArcGIS software, convert the collected vector data into original raster data, and assign corresponding attribute values to the converted original raster data;

[0084] Process the original raster data with corresponding attribute values by the Kriging interpolation method to generate raster data of soil ground temperature isosurfaces at different depths;

[0085] Superimpose the raster data of soil ground temperature isosurfaces and the original raster data according to the preset matrix discrimination relationship to generate a grid map of the freeze-thaw induced grade of geological disasters and a grid map of the prone grade of freeze-thaw geological disasters;

[0086] Overlay the grid map of the freeze-thaw induced grade of geological disasters and the grid map of the susceptibility grade of freeze-thaw geological disasters according to a preset risk discrimination matrix, and render it into a five-level warning map of red, orange, yellow, blue, and white.

[0087] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, storage, database, or other medium used in the various embodiments provided in the present application can include non-volatile and / or volatile memories. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and Rambus dynamic RAM (RDRAM), etc.

[0088] Those skilled in the art can clearly understand that for the convenience and simplicity of description, only the above-mentioned division of each functional unit and module is used as an example. In practical applications, the above functions can be allocated to different functional units and modules according to needs, that is, the internal structure of the device is divided into different functional units or modules to complete all or part of the functions described above.

[0089] Embodiment 3

[0090] The embodiment of the present invention provides an electronic device, including a processor and a memory. At least one instruction or at least one program segment is stored in the memory, and the at least one instruction or the at least one program segment is loaded and executed by the processor to implement the steps:

[0091] Obtain the vector data of the target area;

[0092] Based on the ArcGIS software, convert the collected vector data into original raster data, and assign corresponding attribute values to the converted original raster data;

[0093] The original raster data is processed by Kriging interpolation method to assign corresponding attribute values to generate raster data of soil ground temperature isosurface at different depths;

[0094] The raster data of soil ground temperature isosurface is superimposed on the original raster data according to the preset matrix discrimination relationship to generate a raster map of the freeze-thaw induced grade of geological disasters and a raster map of the susceptibility grade of freeze-thaw geological disasters;

[0095] The raster map of the freeze-thaw induced grade of geological disasters and the raster map of the susceptibility grade of freeze-thaw geological disasters are superimposed according to the preset risk discrimination matrix and rendered into a five-level early warning map of red, orange, yellow, blue and white.

[0096] The above are only the preferred embodiments of the present invention and do not impose any formal limitations on the present invention. Although the present invention has been disclosed above with the preferred embodiments, it is not intended to limit the present invention. Any person skilled in the art can make some changes or modifications to the equivalent embodiments with equivalent changes within the scope of the technical solution of the present invention. However, as long as it does not depart from the content of the technical solution of the present invention, any simple modification, equivalent change and modification made to the above embodiments according to the technical essence of the present invention still fall within the scope of the technical solution of the present invention.

Claims

1. A meteorological risk early warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature, characterized in that, It includes the following steps: S01. Obtain the vector data of the target area; S02. Based on the ArcGIS software, convert the collected vector data into original raster data, and assign corresponding attribute values to the converted original raster data; S03. Process the original raster data with corresponding attribute values assigned through the Kriging interpolation method to generate soil geothermal isothermal surface raster data at different depths; S04. Superimpose the soil geothermal isothermal surface raster data and the original raster data according to a preset matrix discrimination relationship to generate a raster map of the freeze-thaw induced grade of geological disasters and a raster map of the prone grade of freeze-thaw geological disasters; S05. Superimpose the raster map of the freeze-thaw induced grade of geological disasters and the raster map of the prone grade of freeze-thaw geological disasters according to a preset risk discrimination matrix, and render it into a five-level warning map of red, orange, yellow, blue, and white.

2. The meteorological risk warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature according to claim 1, wherein The vector data includes a freeze-thaw zoning map, a zoning map of the prone degree of geological disasters, and a forecast air temperature distribution map for the next 24 hours.

3. The meteorological risk early warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature according to claim 2, characterized in that, The acquisition of the freeze-thaw zoning map is based on the data collected by a satellite equipped with a thermal infrared sensor; The acquisition of the zoning map of the prone degree of geological disasters includes: S11. Establish a database according to the disaster types that occurred in different terrains under different freeze-thaw parameters collected historically; S12. Obtain the current freeze-thaw zoning map, extract the freeze-thaw parameters corresponding to the terrain, and match them with the established database; S13. Generate a new mask in the freeze-thaw zoning map and generate a judgment result for the corresponding terrain according to the matching result. The judgment result includes high prone, medium prone, low prone, and non-prone.

4. The meteorological risk warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature according to claim 1, wherein In step S01, the grid cell size in the original raster data is 3*3 km, and the coordinate system uses the GCS_WGS_1984 geographic coordinate system.

5. The meteorological risk early warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature according to claim 1, wherein In step S02, the processing of the vector data by the ArcGIS software includes: S21. Use the ArcGIS software to import the collected vector data, which contains the geographical spatial elements to be converted; S22. In the ArcToolbox of ArcGIS, select the convert to raster sub-menu under the conversion tool, and execute the polygon to raster according to the predetermined conversion elements of the vector data to obtain raster data; S23. Use the raster calculator tool to assign new attribute values to each pixel of the raster data to obtain the original raster data; The original raster data is the terrain feature of the target area.

6. The meteorological risk warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature according to claim 1, characterized in that In step S03, obtaining the soil geothermal isothermal surface raster data at different depths includes: S31. Arrange the data of 981 monitoring stations in the target area according to a predetermined specification; S32. Clean the obtained data to remove missing values or outliers; S33. Build a Kriging interpolation model based on the Kriging interpolation method to predict the soil temperature values of the monitoring stations; S34. According to the temperature value spacing of the 981 monitoring stations, deduce the temperature difference between two adjacent monitoring stations and the temperature decreasing data of the linear distance between two adjacent monitoring stations under the temperature difference.

7. The meteorological risk warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature according to claim 1, wherein In step S04, superimposing and calculating the soil geothermal isothermal surface raster data and the original raster data according to a preset matrix discrimination relationship includes: S41. Superimposing and calculating the geothermal raster data and the predicted air temperature raster data according to a preset matrix discrimination relationship to generate a raster map of the induced level of geological disaster freeze-thaw; S42. Superimposing and calculating the raster map of the zoning of the susceptibility degree of geological disasters and the raster map of the freeze-thaw zoning according to a preset risk discrimination matrix relationship to generate a raster map of the susceptibility level of freeze-thaw geological disasters.

8. The meteorological risk warning method for freeze-thaw geological disasters based on the changes in air temperature and ground temperature according to claim 1, wherein, The risk levels of red, orange, yellow, blue, and white in step S05 are sorted from high to low.

9. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the steps of the freeze-thaw geological disaster meteorological risk warning method based on the changes in air temperature and geothermal temperature described in any one of claims 1 to 8.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the freeze-thaw geological disaster meteorological risk warning method based on the changes in air temperature and geothermal temperature described in any one of claims 1 to 8.

Images