Method, system and device for calculating volume of glacial lake based on morphological and topographical constraints, and storage medium

Abstract

The application discloses a method, system and device for calculating the volume of an ice lake based on morphology and terrain constraints and a storage medium, and belongs to the technical field of ice lake volume calculation. The method comprises the following steps: obtaining the area and perimeter of an ice lake based on pre-processed remote sensing image data; obtaining the morphology coefficient of the ice lake based on the area and perimeter of the ice lake; collecting digital elevation model (DEM) data of a range above the overflow outlet of the terminal moraine of the target ice lake to obtain the fluctuation coefficient of the ice moraine around the ice lake and the morphology index of the ice lake basin; obtaining DEM data of the range of the slope of the ice moraine around the target ice lake to further obtain the terrain coefficient around the ice lake; and obtaining the volume of the target ice lake based on the area of the ice lake, the morphology coefficient of the ice lake, the fluctuation coefficient of the ice moraine around the ice lake, the terrain coefficient around the ice lake and the morphology index of the ice lake basin. The application can significantly improve the calculation accuracy and more accurately calculate the volume of the ice lake, thereby providing key basis for accurately evaluating the risk of ice lake outburst and formulating targeted prevention measures.

Description

Technical Field

[0001] This invention belongs to the field of ice lake volume calculation technology, and particularly relates to ice lake volume calculation methods, systems, devices and storage media based on morphological and topographic constraints. Background Technology

[0002] Glacial lakes are lakes formed by glacial activity, typically in high-altitude regions. Against the backdrop of global warming, glacial lakes are generally expanding, leading to an increased risk of glacial lake outburst floods. In recent years, glacial lake outburst floods have occurred frequently, seriously threatening the lives and property of people downstream. The volume of water within a glacial lake is a key element, controlling the scale and destructive power of the floodwaters following an outburst, and thus determining the degree of risk associated with such a flood.

[0003] Using instruments such as underwater sonar to conduct on-site measurements of glacial lakes can accurately obtain their volume data. However, since most glacial lakes are located in high-altitude and remote areas, reaching the site for on-site measurements is extremely difficult. Glacial lake survey and cataloging data show that there are approximately 1,680 glacial lakes in the Himalayas of my country, but currently only a very small number, accounting for less than 3% of the total, have undergone on-site measurements. Therefore, current methods mostly utilize the area and volume data of a limited number of already measured glacial lakes to fit empirical formulas relating glacial lake area to volume. Based on this, remote sensing interpretation methods are used to obtain the area of ​​unknown glacial lakes, and the aforementioned empirical formulas are used to calculate the volume of unknown glacial lakes. Examples of such empirical formulas include the Huggel formula, the Evans formula, and the O'Connor formula, which are widely used internationally.

[0004] However, the volume of a glacial lake is essentially controlled by the characteristics of its basin, and the lake's shape and surrounding topography are direct reflections of this basin. Glacial lakes of the same area can have vastly different actual volumes due to variations in shape and basin topography. Therefore, calculating a glacial lake's volume solely based on its area can lead to significant errors. Summary of the Invention

[0005] The purpose of this invention is to provide a method, system, device, and storage medium for calculating the volume of glacial lakes based on morphological and topographical constraints, in order to solve the problem that calculating the volume of glacial lakes using only the area of ​​the lake may cause large errors.

[0006] The embodiments of this application are implemented as follows: a method for calculating the volume of glacial lakes based on morphological and topographical constraints, including:

[0007] S01. Obtain preprocessed remote sensing image data, and based on the preprocessed remote sensing image data, obtain the area and perimeter of the ice lake;

[0008] S02. Obtain the morphological coefficient of the ice lake based on its area and perimeter.

[0009] S03. Collect digital elevation model (DEM) data of the area above the spillway of the final moraine dam of the target glacial lake to obtain the undulation coefficient of the moraine around the glacial lake and the morphological index of the glacial lake basin.

[0010] S04. Obtain the digital elevation model (DEM) data of the slope range of glacial moraine around the target glacial lake, and then obtain the topographic coefficients around the glacial lake.

[0011] S05. Based on the area of ​​the glacial lake, the morphological coefficient of the glacial lake, the undulation coefficient of the glacial moraine around the glacial lake, the topographic coefficient around the glacial lake, and the morphological index of the glacial lake basin, the target volume of the glacial lake is obtained.

[0012] Optionally, in some embodiments of this application, remote sensing image data of the area surrounding the target glacial lake is collected, and the remote sensing image data is preprocessed to obtain preprocessed remote sensing image data; and / or

[0013] The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was used. The area of ​​the glacial lake was delineated using polygons, and the area of ​​the polygons was calculated using ArcGIS's computational geometry functions. This area represents the glacial lake area. Calculate the perimeter of the polygon, which is the perimeter of the ice lake. P ; and / or

[0014] The formula for calculating the topographic coefficient around the glacial lake is as follows:

[0015] ω = tanα

[0016] In the formula: ω The topographic coefficient around the glacial lake, where α is the average slope of the slope range around the glacial lake; and / or

[0017] The formula for calculating the target glacial lake volume is as follows:

[0018] ;

[0019] In the formula: For the target glacial lake volume, The area of ​​the glacial lake, For the morphological coefficient of glacial lakes, The coefficient of variation of glacial moraine around the glacial lake. The topographic coefficient around the glacial lake. This is a morphological index for the glacial lake basin.

[0020] Optionally, in some embodiments of this application, the preprocessing includes geometric correction, orthorectification, and image fusion of the remote sensing image data using ENVI software.

[0021] Optionally, in some embodiments of this application, the formula for calculating the glacial lake morphology coefficient is as follows:

[0022] ;

[0023] In the formula, For the morphological coefficient of glacial lakes, For the area of ​​the glacial lake, The perimeter of the ice lake; and / or

[0024] The formula for calculating the undulation coefficient of glacial moraine around the glacial lake is as follows:

[0025] ;

[0026] In the formula, The coefficient of variation of glacial moraine around the glacial lake. The elevation difference of the glacial moraine slopes surrounding the glacial lake. This refers to the elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine.

[0027] Optionally, in some embodiments of this application, ArcGIS software is used to overlay the glacial lake extent onto the digital elevation model (DEM) data of the area above the spillway of the target glacial lake terminal moraine, and the glacial lake surface elevation is extracted using the glacial lake extent. The elevation of the glacial lake terminal moraine overflow outlet was directly read using the digital elevation model (DEM) data of the area above the target glacial lake terminal moraine overflow outlet. The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine is... for:

[0028] ;

[0029] In the formula: The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine dike; This refers to the elevation of the ice lake surface. Elevation of the spillway of the glacial lake terminal moraine;

[0030] The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was employed. Polygons were used to delineate the slopes of the lateral and terminal moraines surrounding the glacial lake, facing the lake. ArcGIS's cropping tool was used to crop the Digital Elevation Model (DEM) data of the area above the spillway of the target glacial lake's terminal moraines to the defined polygon area, thus obtaining the DEM data of the glacial moraine slopes surrounding the target glacial lake. The maximum elevation of the glacial moraine slopes around the target glacial lake was then read from this DEM data. The elevation difference of the glacial moraine slopes around the glacial lake for:

[0031] ;

[0032] In the formula: The elevation difference of the glacial moraine slopes surrounding the glacial lake; This represents the maximum elevation of the glacial moraine slope surrounding the glacial lake. This refers to the elevation of the ice lake surface.

[0033] Optionally, in some embodiments of this application, the digital elevation model (DEM) data of the area above the spillway of the target glacial lake terminal moraine is cropped using ArcGIS software's cropping tool to the watershed area above the glacial lake surface, thus obtaining the DEM data of the watershed area above the glacial lake surface. Contour lines for the watershed area above the glacial lake surface are then generated using ArcGIS software's contour line tool, and the area of ​​the watershed area above the glacial lake surface is calculated using ArcGIS software. The maximum elevation values ​​of the DEM data for the watershed above the surface of the glacial lake were extracted using ArcGIS software. and the elevation of the ice lake surface The difference between the maximum elevation of the DEM data for the watershed above the glacial lake surface and the elevation of the glacial lake surface is... for:

[0034] ;

[0035] Using DEM data of the catchment area above the glacial lake surface, the area 'a' above each contour line and the elevation difference between each contour line and the glacial lake surface elevation are extracted. h Then a i h is the area above the contour line with elevation value i. i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake.

[0036] make , ;

[0037] Multiple ( Connecting these curves yields the function curve: ;

[0038] For functions Integrating the results, we obtain the Glacier Lake Basin Morphological Index HI: ;

[0039] In the formula, HI For the morphological index of the glacial lake basin; a i Let i be the area above the contour line with elevation value i. The area of ​​the catchment above the surface of the glacial lake; h i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake. This represents the elevation difference between the maximum elevation value of the DEM data for the catchment area above the surface of the glacial lake and the elevation of the glacial lake surface.

[0040] Optionally, in some embodiments of this application, the digital elevation model (DEM) data of the glacial moraine slope range around the target glacial lake is used, and ArcGIS software is used to calculate the slope data of the glacial moraine slope range around the target glacial lake. Then, ArcGIS software is used to statistically analyze and calculate the slope data of the glacial moraine slope range around the glacial lake to obtain the average slope α of the slope range around the glacial lake.

[0041] Accordingly, embodiments of this application also provide a system for calculating the volume of glacial lakes based on morphological and topographical constraints, including:

[0042] The module for measuring the area and perimeter of a glacial lake obtains preprocessed remote sensing image data, and based on the preprocessed remote sensing image data, it obtains the area and perimeter of the glacial lake.

[0043] The ice lake morphology coefficient module obtains the ice lake morphology coefficient based on the ice lake area and ice lake perimeter;

[0044] The modules for the undulation coefficient of glacial moraine around the glacial lake and the morphological index of the glacial lake basin collect digital elevation model (DEM) data of the area above the spillway of the final moraine dike of the target glacial lake to obtain the undulation coefficient of glacial moraine around the glacial lake and the morphological index of the glacial lake basin.

[0045] The Glacier Lake Surrounding Topographic Coefficient module obtains the Digital Elevation Model (DEM) data of the glacial moraine slope range surrounding the target glacier lake, and then obtains the glacier lake surrounding topographic coefficients.

[0046] The target glacial lake volume module obtains the target glacial lake volume based on the glacial lake area, glacial lake morphology coefficient, glacial moraine undulation coefficient around the glacial lake, glacial lake topography coefficient, and glacial lake watershed morphology index.

[0047] Accordingly, embodiments of this application also provide a computer device, including a storage device and a processor, wherein the storage device stores a computer program, and when the computer program is executed by the processor, the processor performs the steps of the method described above.

[0048] Accordingly, embodiments of this application also provide a computer-readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform the steps of the method described above.

[0049] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

[0050] This invention provides a method for calculating glacial lake volume based on glacial lake morphology and surrounding topographic constraints. Compared to current methods that only estimate glacial lake volume by area, this invention introduces five key factors controlling glacial lake volume as variables: glacial lake area, glacial lake morphology coefficient, glacial moraine undulation coefficient around the glacial lake, glacial lake surrounding topographic coefficient, and glacial lake basin morphology index. This constructs a more scientific method for calculating glacial lake volume. Compared to current methods, which often result in errors of several times or even orders of magnitude, this invention significantly improves calculation accuracy, providing a more precise calculation of glacial lake volume and offering crucial evidence for accurately assessing glacial lake outburst flood risk and developing targeted prevention and control measures. Attached Figure Description

[0051] Figure 1 This is a flowchart of the method for calculating the volume of glacial lakes based on morphological and topographical constraints according to the present invention;

[0052] Figure 2 This is a graph showing the relationship between the area of ​​a glacial lake and its actual volume, collected or measured according to the present invention.

[0053] Figure 3 This is a graph showing the relationship between the measured glacial lake morphology coefficient, glacial lake topography coefficient, glacial moraine undulation coefficient around the glacial lake, glacial lake basin morphology index, and glacial lake volume in this invention.

[0054] Figure 4 This is a diagram showing the interpretation and calculation results of the area and perimeter of a certain glacial lake according to the present invention;

[0055] Figure 5 This is a diagram showing the slope range and gradient of the lateral moraines and terminal moraines surrounding a certain glacial lake, facing the glacial lake.

[0056] Figure 6 This is a statistical result of the average slope of the lateral moraines and terminal moraines surrounding a certain glacial lake on the side facing the glacial lake.

[0057] Figure 7 This is a schematic diagram showing the values ​​of the morphological index calculation factors for a certain glacial lake basin in this invention.

[0058] Figure 8 This is a graph showing the calculation results of the morphological index of a certain glacial lake basin according to the present invention;

[0059] Figure 9 This is a comparison chart of the calculation error of the present invention and the error of the commonly used calculation result using only area. Detailed Implementation

[0060] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0061] The technical solution of this application is as follows:

[0062] Firstly, please refer to Figure 1 This application provides a method for calculating the volume of glacial lakes based on morphological and topographical constraints, including:

[0063] S01. Obtain preprocessed remote sensing image data, and based on the preprocessed remote sensing image data, obtain the area and perimeter of the ice lake;

[0064] S02. Obtain the morphological coefficient of the ice lake based on its area and perimeter.

[0065] S03. Collect digital elevation model (DEM) data of the area above the spillway of the final moraine dam of the target glacial lake to obtain the undulation coefficient of the moraine around the glacial lake and the morphological index of the glacial lake basin.

[0066] S04. Obtain the digital elevation model (DEM) data of the slope range of glacial moraine around the target glacial lake, and then obtain the topographic coefficients around the glacial lake.

[0067] S05. Based on the area of ​​the glacial lake, the morphological coefficient of the glacial lake, the undulation coefficient of the glacial moraine around the glacial lake, the topographic coefficient around the glacial lake, and the morphological index of the glacial lake basin, the target volume of the glacial lake is obtained.

[0068] Glacial lake volume is a core factor in accurately assessing glacial lake outburst flood risk and its likelihood, predicting peak outburst flow, and formulating disaster prevention measures for downstream valley residents, buildings, hydropower, and transportation. Therefore, accurately calculating the volume of glacial lakes is extremely crucial. This invention provides a method for calculating glacial lake volume based on glacial lake morphology and surrounding topographic constraints. Compared to the traditional area-volume empirical formula, this invention introduces five key factors controlling glacial lake volume as variables: glacial lake area, glacial lake morphology coefficient, glacial moraine undulation coefficient around the glacial lake, glacial lake surrounding topographic coefficient, and glacial lake basin morphology index. The physical meaning of each variable corresponds to the calculation result, and the dimensions of the variables are consistent with the calculation result, thus constructing a more scientific method for calculating glacial lake volume. Compared to current methods, this invention can significantly improve calculation accuracy and more precisely calculate glacial lake volume, providing a crucial basis for accurately assessing glacial lake outburst flood risk and developing targeted prevention and control measures.

[0069] This application utilizes the results of 24 field measurements of glacial lake volume, combined with glacial lake morphology and surrounding topographic features, to determine that:

[0070] Factors influencing the volume of a glacial lake include its morphology and surrounding topographic features, specifically its area, perimeter, surface elevation, overflow elevation, maximum elevation and average slope of surrounding glacial moraine slopes, and the elevation and area of ​​its upper catchment area. Based on these factors, further calculations can yield influencing factors such as the morphological coefficient of the glacial lake, the undulation coefficient of the surrounding glacial moraine, the topographic coefficient of the glacial lake, and the morphological index of the glacial lake catchment area.

[0071] Among the factors influencing glacial lake volume, glacial lake area has a significant impact, generally showing a positive correlation. Further analysis of field-measured glacial lake data reveals a near-exponential relationship between glacial lake area and volume. However, glacial lakes with similar areas exhibit significant volume differences, indicating that glacial lake volume is influenced not only by glacial lake area but also by other factors, such as… Figure 2 As shown. However, glacial lakes of similar area exhibit significant fluctuations in volume, indicating that glacial lake volume is influenced not only by lake area but also by other factors such as lake shape and topography. Figure 3 As shown. Therefore, based on the measured volumes of 24 glacial lakes in the field, and the corresponding data on factors such as glacial lake area, glacial lake morphology coefficient, glacial moraine undulation coefficient, glacial lake topography coefficient, and glacial lake watershed morphology index, and considering that the relationship between glacial lake area and glacial lake volume is a near-exponential function, SPSS software was used to perform multivariate nonlinear fitting analysis on the glacial lake volume and factor data, resulting in the glacial lake volume calculation method based on glacial lake morphology and topography proposed in this application.

[0072] In S01:

[0073] In some embodiments, remote sensing image data of the area surrounding the target glacial lake is collected, and the remote sensing image data is preprocessed to obtain preprocessed remote sensing image data.

[0074] Furthermore, the preprocessing includes using ENVI software to perform geometric correction, orthorectification, and image fusion on the remote sensing image data.

[0075] Furthermore, geometric correction involves using control point data between remote sensing image data and standard maps in the surveying and mapping industry to obtain a mathematical model of geometric distortion, and then using this mathematical model to correct the geometric distortion of the remote sensing image data.

[0076] It is understandable that the purpose of using this geometric distortion mathematical model to correct the geometric distortion of the remote sensing image data is to reduce or eliminate the geometric distortion of the acquired remote sensing image (image) relative to the ground target caused by the influence of the attitude, altitude, speed of the aircraft acquiring the remote sensing image data and the Earth's rotation.

[0077] Furthermore, orthorectification involves selecting ground control points on the remote sensing image data, using digital elevation model (DEM) data within the range of the remote sensing image, simultaneously performing tilt correction and projection difference correction on the remote sensing image, and resampling the remote sensing image into an orthorectified image.

[0078] It is understandable that performing tilt correction and projection difference correction on the remote sensing image simultaneously aims to reduce or eliminate significant distortions caused by terrain, camera geometry, and sensor-related errors.

[0079] Furthermore, image fusion involves processing redundant or complementary remote sensing image data in space, time, and spectrum using IHS transform to obtain synthetic remote sensing image data with new spatial, temporal, and spectral characteristics.

[0080] In some embodiments, preprocessed remote sensing image data is loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation is used. The area of ​​the glacial lake is delineated using polygons, and the area of ​​the polygons is calculated using ArcGIS computational geometry functions, which represents the area of ​​the glacial lake. Calculate the perimeter of the polygon, which is the perimeter of the ice lake. P .

[0081] In S02:

[0082] In some embodiments, the formula for calculating the glacial lake morphology coefficient is as follows:

[0083] ;

[0084] In the formula, For the morphological coefficient of glacial lakes, For the area of ​​the glacial lake, This is the circumference of the ice lake.

[0085] In S03:

[0086] In some embodiments, the formula for calculating the undulation coefficient of glacial moraine around a glacial lake is as follows:

[0087] ;

[0088] In the formula, The coefficient of variation of glacial moraine around the glacial lake. The elevation difference of the glacial moraine slopes surrounding the glacial lake. This refers to the elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine.

[0089] Furthermore, using ArcGIS software, the digital elevation model (DEM) data of the glacial lake area was overlaid onto the area above the spillway of the target glacial lake's terminal moraine, and the lake surface elevation was extracted using the glacial lake area. The elevation of the glacial lake terminal moraine overflow outlet was directly read using the digital elevation model (DEM) data of the area above the target glacial lake terminal moraine overflow outlet. The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine is... for:

[0090] ;

[0091] In the formula: The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine dike; This refers to the elevation of the ice lake surface. Elevation of the spillway of the glacial lake terminal moraine;

[0092] The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was employed. Polygons were used to delineate the slopes of the lateral and terminal moraines surrounding the glacial lake, facing the lake. ArcGIS's cropping tool was used to crop the Digital Elevation Model (DEM) data of the area above the spillway of the target glacial lake's terminal moraines to the defined polygon area, thus obtaining the DEM data of the glacial moraine slopes surrounding the target glacial lake. The maximum elevation of the glacial moraine slopes around the target glacial lake was then read from this DEM data. The elevation difference of the glacial moraine slopes around the glacial lake for:

[0093] ;

[0094] In the formula: The elevation difference of the glacial moraine slopes surrounding the glacial lake; This represents the maximum elevation of the glacial moraine slope surrounding the glacial lake. This refers to the elevation of the ice lake surface.

[0095] Understandable. The elevation difference between the surface of the glacial lake and the spillway of the terminal moraine is expressed in meters (m). This refers to the elevation of the ice lake surface, in meters (m). Elevation of the spillway of the final moraine of the glacial lake, in meters.

[0096] Understandable. The elevation difference of the glacial moraine slopes around the glacial lake, in meters; The maximum elevation of the glacial moraine slope surrounding the glacial lake is shown in meters (m).

[0097] In some embodiments, the ArcGIS software clipping tool is used to clip the digital elevation model (DEM) data of the area above the spillway of the target glacial lake terminal moraine to the watershed area above the lake surface, obtaining the DEM data of the watershed area above the lake surface. The contour line tool in ArcGIS software is then used to generate contour lines for the watershed area above the lake surface. Finally, ArcGIS software is used to calculate the area of ​​the watershed area above the lake surface. The maximum elevation values ​​of the DEM data for the watershed above the surface of the glacial lake were extracted using ArcGIS software. and the elevation of the ice lake surface The difference between the maximum elevation of the DEM data for the watershed above the glacial lake surface and the elevation of the glacial lake surface is... for:

[0098] ;

[0099] Using DEM data of the catchment area above the glacial lake surface, the area 'a' above each contour line and the elevation difference between each contour line and the glacial lake surface elevation are extracted. h Then a i h is the area above the contour line with elevation value i. i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake.

[0100] make , ;

[0101] Multiple ( Connecting these curves yields the function curve: ;

[0102] For functions Integrating the results, we obtain the Glacier Lake Basin Morphological Index HI: ;

[0103] In the formula, HI For the morphological index of the glacial lake basin; a i Let i be the area above the contour line with elevation value i. The area of ​​the catchment above the surface of the glacial lake; h i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake. This represents the elevation difference between the maximum elevation value of the DEM data for the catchment area above the surface of the glacial lake and the elevation of the glacial lake surface.

[0104] Understandable. HI a is a dimensionless watershed morphological index above the surface of the glacial lake; i The area above the contour line with elevation value i, in km². 2 ; The watershed area above the surface of the glacial lake is expressed in km². 2 ; h i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake, in meters; The elevation difference, in meters, is the maximum elevation value of the DEM data for the catchment area above the surface of the glacial lake and the elevation of the glacial lake surface.

[0105] In S04:

[0106] In some embodiments, the formula for calculating the topographic coefficient around the glacial lake is as follows:

[0107] ω = tanα;

[0108] In the formula: ω α represents the topographic coefficient around the glacial lake, and α represents the average slope of the slope range around the glacial lake.

[0109] Furthermore, using the digital elevation model (DEM) data of the glacial moraine slope range around the target glacial lake, ArcGIS software was used to calculate the slope data of the glacial moraine slope range around the target glacial lake. Then, ArcGIS software was used to statistically analyze and calculate the slope data of the glacial moraine slope range around the glacial lake to obtain the average slope α of the slope range around the glacial lake.

[0110] In S05:

[0111] In some embodiments, the target glacial lake volume is calculated using the following formula:

[0112] ;

[0113] In the formula: For the target glacial lake volume, The area of ​​the glacial lake, For the morphological coefficient of glacial lakes, The coefficient of variation of glacial moraine around the glacial lake. The topographic coefficient around the glacial lake. This is a morphological index for the glacial lake basin.

[0114] Secondly, embodiments of this application provide a system for calculating the volume of glacial lakes based on morphological and topographical constraints, including:

[0115] The module for measuring the area and perimeter of a glacial lake obtains preprocessed remote sensing image data, and based on the preprocessed remote sensing image data, it obtains the area and perimeter of the glacial lake.

[0116] The ice lake morphology coefficient module obtains the ice lake morphology coefficient based on the ice lake area and ice lake perimeter;

[0117] The modules for the undulation coefficient of glacial moraine around the glacial lake and the morphological index of the glacial lake basin collect digital elevation model (DEM) data of the area above the spillway of the final moraine dike of the target glacial lake to obtain the undulation coefficient of glacial moraine around the glacial lake and the morphological index of the glacial lake basin.

[0118] The Glacier Lake Surrounding Topographic Coefficient module obtains the Digital Elevation Model (DEM) data of the glacial moraine slope range surrounding the target glacier lake, and then obtains the glacier lake surrounding topographic coefficients.

[0119] The target glacial lake volume module obtains the target glacial lake volume based on the glacial lake area, glacial lake morphology coefficient, glacial moraine undulation coefficient around the glacial lake, glacial lake topography coefficient, and glacial lake watershed morphology index.

[0120] In some embodiments of the ice lake area and ice lake perimeter module, remote sensing image data of the surrounding area of ​​the target ice lake is collected, and the remote sensing image data is preprocessed to obtain preprocessed remote sensing image data.

[0121] Furthermore, the preprocessing includes using ENVI software to perform geometric correction, orthorectification, and image fusion on the remote sensing image data.

[0122] Furthermore, geometric correction involves using control point data between remote sensing image data and standard maps in the surveying and mapping industry to obtain a mathematical model of geometric distortion, and then using this mathematical model to correct the geometric distortion of the remote sensing image data.

[0123] It is understandable that the purpose of using this geometric distortion mathematical model to correct the geometric distortion of the remote sensing image data is to reduce or eliminate the geometric distortion of the acquired remote sensing image (image) relative to the ground target caused by the influence of the attitude, altitude, speed of the aircraft acquiring the remote sensing image data and the Earth's rotation.

[0124] Furthermore, orthorectification involves selecting ground control points on the remote sensing image data, using digital elevation model (DEM) data within the range of the remote sensing image, simultaneously performing tilt correction and projection difference correction on the remote sensing image, and resampling the remote sensing image into an orthorectified image.

[0125] It is understandable that performing tilt correction and projection difference correction on the remote sensing image simultaneously aims to reduce or eliminate significant distortions caused by terrain, camera geometry, and sensor-related errors.

[0126] Furthermore, image fusion involves processing redundant or complementary remote sensing image data in space, time, and spectrum using IHS transform to obtain synthetic remote sensing image data with new spatial, temporal, and spectral characteristics.

[0127] In some embodiments, preprocessed remote sensing image data is loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation is used. The area of ​​the glacial lake is delineated using polygons, and the area of ​​the polygons is calculated using ArcGIS computational geometry functions, which represents the area of ​​the glacial lake. Calculate the perimeter of the polygon, which is the perimeter of the ice lake. P。

[0128] In the ice lake morphology coefficient module

[0129] In some embodiments, the formula for calculating the glacial lake morphology coefficient is as follows:

[0130] ;

[0131] In the formula, For the morphological coefficient of glacial lakes, For the area of ​​the glacial lake, This is the circumference of the ice lake.

[0132] In the modules of glacial moraine undulation coefficient and glacial lake watershed morphology index,

[0133] In some embodiments, the formula for calculating the undulation coefficient of glacial moraine around a glacial lake is as follows:

[0134] ;

[0135] In the formula, The coefficient of variation of glacial moraine around the glacial lake. The elevation difference of the glacial moraine slopes surrounding the glacial lake. This refers to the elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine.

[0136] Furthermore, using ArcGIS software, the digital elevation model (DEM) data of the glacial lake area was overlaid onto the area above the spillway of the target glacial lake's terminal moraine, and the lake surface elevation was extracted using the glacial lake area. The elevation of the glacial lake terminal moraine overflow outlet was directly read using the digital elevation model (DEM) data of the area above the target glacial lake terminal moraine overflow outlet. The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine is... for:

[0137] ;

[0138] In the formula: The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine dike; This refers to the elevation of the ice lake surface. Elevation of the spillway of the glacial lake terminal moraine;

[0139] The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was employed. Polygons were used to delineate the slopes of the lateral and terminal moraines surrounding the glacial lake, facing the lake. ArcGIS's cropping tool was used to crop the Digital Elevation Model (DEM) data of the area above the spillway of the target glacial lake's terminal moraines to the defined polygon area, thus obtaining the DEM data of the glacial moraine slopes surrounding the target glacial lake. The maximum elevation of the glacial moraine slopes around the target glacial lake was then read from this DEM data. The elevation difference of the glacial moraine slopes around the glacial lake for:

[0140] ;

[0141] In the formula: The elevation difference of the glacial moraine slopes surrounding the glacial lake; This represents the maximum elevation of the glacial moraine slope surrounding the glacial lake. This refers to the elevation of the ice lake surface.

[0142] Understandable. The elevation difference between the surface of the glacial lake and the spillway of the terminal moraine is expressed in meters (m). This refers to the elevation of the ice lake surface, in meters (m). Elevation of the spillway of the final moraine of the glacial lake, in meters.

[0143] Understandable. The elevation difference of the glacial moraine slopes around the glacial lake, in meters; The maximum elevation of the glacial moraine slope surrounding the glacial lake is shown in meters (m).

[0144] In some embodiments, the ArcGIS software clipping tool is used to clip the digital elevation model (DEM) data of the area above the spillway of the target glacial lake terminal moraine to the watershed area above the lake surface, obtaining the DEM data of the watershed area above the lake surface. The contour line tool in ArcGIS software is then used to generate contour lines for the watershed area above the lake surface. Finally, ArcGIS software is used to calculate the area of ​​the watershed area above the lake surface. The maximum elevation values ​​of the DEM data for the watershed above the surface of the glacial lake were extracted using ArcGIS software. and the elevation of the ice lake surface The difference between the maximum elevation of the DEM data for the watershed above the glacial lake surface and the elevation of the glacial lake surface is... for:

[0145] ;

[0146] Using DEM data of the catchment area above the glacial lake surface, the area 'a' above each contour line and the elevation difference between each contour line and the glacial lake surface elevation are extracted. h Then a i h is the area above the contour line with elevation value i. i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake.

[0147] make , ;

[0148] Multiple ( Connecting these curves yields the function curve: ;

[0149] For functions Integrating the results, we obtain the Glacier Lake Basin Morphological Index HI: ;

[0150] In the formula, HI For the morphological index of the glacial lake basin; a i Let i be the area above the contour line with elevation value i. The area of ​​the catchment above the surface of the glacial lake; h i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake. This represents the elevation difference between the maximum elevation value of the DEM data for the catchment area above the surface of the glacial lake and the elevation of the glacial lake surface.

[0151] Understandable. HI a is a dimensionless watershed morphological index above the surface of the glacial lake; i The area above the contour line with elevation value i, in km². 2 ; The watershed area above the surface of the glacial lake is expressed in km². 2 ; h i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake, in meters; The elevation difference, in meters, is the maximum elevation value of the DEM data for the catchment area above the surface of the glacial lake and the elevation of the glacial lake surface.

[0152] In the terrain coefficient module surrounding the glacial lake

[0153] In some embodiments, the formula for calculating the topographic coefficient around the glacial lake is as follows:

[0154] ω = tanα;

[0155] In the formula: ω α represents the topographic coefficient around the glacial lake, and α represents the average slope of the slope range around the glacial lake.

[0156] Furthermore, using the digital elevation model (DEM) data of the glacial moraine slope range around the target glacial lake, ArcGIS software was used to calculate the slope data of the glacial moraine slope range around the target glacial lake. Then, ArcGIS software was used to statistically analyze and calculate the slope data of the glacial moraine slope range around the glacial lake to obtain the average slope α of the slope range around the glacial lake.

[0157] In the target glacial lake volume module

[0158] In some embodiments, the target glacial lake volume is calculated using the following formula:

[0159] ;

[0160] In the formula: For the target glacial lake volume, The area of ​​the glacial lake, For the morphological coefficient of glacial lakes, The coefficient of variation of glacial moraine around the glacial lake. The topographic coefficient around the glacial lake. This is a morphological index for the glacial lake basin.

[0161] Thirdly, this application provides a computer device including a storage device and a processor, wherein the storage device stores a computer program, and when the computer program is executed by the processor, the processor performs the steps of the above-described method for calculating the volume of an ice lake based on morphological and topographical constraints.

[0162] The computer device can be a desktop computer, laptop, handheld computer, or cloud server, etc. The computer device can interact with the user via a keyboard, mouse, remote control, touchpad, or voice control.

[0163] The memory includes at least one type of readable storage medium, including flash memory, hard disk, multimedia card, card-type memory (e.g., SD or D-interface display memory), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, optical disk, etc. In some embodiments, the memory may be an internal storage unit of the computer device, such as the hard disk or memory of the computer device. In other embodiments, the memory may also be an external storage device of the computer device, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc., equipped on the computer device. Of course, the memory may include both internal storage units and external storage devices of the computer device. In this embodiment, the memory is often used to store the operating system and various application software installed on the computer device, such as the program code of the ice lake volume calculation method based on morphology and terrain constraints. In addition, the memory can also be used to temporarily store various types of data that have been output or will be output.

[0164] In some embodiments, the processor may be a central processing unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chip. The processor is typically used to control the overall operation of the computer device. In this embodiment, the processor is used to run program code stored in the memory or process data, for example, to run the program code for the morphological and topographical constraint-based ice lake volume calculation method.

[0165] Fourthly, this application provides a computer-readable storage medium storing a computer program, which, when executed by a processor, causes the processor to perform the steps of the above-described method for calculating the volume of an ice lake based on morphological and topographical constraints.

[0166] The computer-readable storage medium stores an interface display program that can be executed by at least one processor to cause the at least one processor to perform the steps of the glacial lake volume calculation method based on morphological and topographic constraints as described above.

[0167] The present invention will be further described below with reference to experimental examples.

[0168] Experimental Example 1

[0169] This experimental example provides a method for calculating the volume of glacial lakes based on morphological and topographical constraints, and selects a glacial lake in my country for calculation and comparative verification.

[0170] Publicly available remote sensing imagery and digital elevation model (DEM) data of a glacial lake were collected. ENVI software was used to perform geometric correction, orthorectification, and image fusion processing on the remote sensing imagery. Geometric correction involved obtaining a geometric distortion mathematical model using control point data between the remote sensing imagery and standard surveying and mapping maps, and then using this model to correct the geometric distortion of the remote sensing imagery. Orthorectification involved selecting ground control points on the remote sensing imagery, using the DEM data within the imagery area, and simultaneously performing tilt correction and projection difference correction on the remote sensing imagery, resampling it into an orthorectified image. Image fusion involved processing spatially, temporally, and spectrally redundant or complementary remote sensing imagery data using the commonly used IHS transform method to obtain synthetic remote sensing imagery with new spatial, temporal, and spectral characteristics.

[0171] The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of interactive and visual interpretation methods was used, and polygons were employed to delineate the extent of the glacial lake. Figure 4 As shown. Using ArcGIS's computational geometry capabilities, the area of ​​the polygon is calculated, which represents the area of ​​the glacial lake. = 1.22 km2 The perimeter of the polygon was measured using ArcGIS and found to be P = 674 m.

[0172] Glacial lake morphology coefficient ;

[0173] Using ArcGIS software, the extent of the glacial lake was overlaid onto the digital elevation model (DEM) data, such as... Figure 5 and Figure 6 As shown. The elevation of the glacial lake surface is extracted using the extent of the glacial lake. Read the elevation of the spillway of the glacial lake terminal moraine. The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine is... for:

[0174] ;

[0175] The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was employed. Polygons were used to delineate the slopes of the lateral and terminal moraines surrounding the glacial lake, facing the lake. ArcGIS's cropping tool was used to crop the Digital Elevation Model (DEM) data to the defined polygonal area, resulting in the DEM data for the glacial moraine slopes surrounding the lake. Using this DEM data, the maximum elevation of the slopes surrounding the lake was determined. The elevation difference of the glacial moraine slopes around the glacial lake for:

[0176] ;

[0177] Glacial moraine undulation coefficient around the glacial lake:

[0178] ;

[0179] Using the digital elevation model (DEM) data of the glacial moraine slopes surrounding the glacial lake, ArcGIS software was used to calculate the slope data of the glacial moraine slopes surrounding the glacial lake. Then, ArcGIS software was used to statistically analyze and calculate the slope data of the glacial moraine slopes surrounding the glacial lake, and the average slope α = 20.2° was obtained.

[0180] Calculate the topographic coefficient ω around the glacial lake:

[0181] ω = tanα=0.3679;

[0182] Using ArcGIS software's cropping tool, the Digital Elevation Model (DEM) data was cropped to cover the catchment area above the glacial lake surface. Using ArcGIS's contour line tool, contour lines were generated for the catchment area above the glacial lake surface. The catchment area above the glacial lake surface was calculated to be A = 49.15 km². 2 The maximum elevation values ​​of the DEM (Digital Elevation Model) for the watershed above the surface of the glacial lake were extracted using ArcGIS software. =7452m and the elevation of the Ice Lake surface =5483m, then the elevation difference between the maximum elevation above the surface of the glacial lake and the elevation of the glacial lake surface is... for:

[0183] ;

[0184] Using the digital elevation model (DEM) data of the watershed above the glacial lake surface, the area 'a' above each elevation line is extracted using the method described above. i (i is the elevation value of a certain contour line), and the elevation difference between the elevation value of each corresponding contour line and the elevation of the glacial lake surface. h i .like Figure 7 As shown in the example, the area of ​​a glacial lake above an elevation of 5800m is denoted as 'a'. 5800 =39.18km 2 The elevation difference between the contour line at 5800m and the glacial lake is... h 5800 =5800-5483=317m.

[0185] make , ;

[0186] but , .

[0187] Multiple ( Connecting these curves yields the function curve: ;

[0188] like Figure 8 As shown, for the function Integrating the results, we obtain the Glacier Lake Basin Morphological Index HI: ;

[0189] The volume of a glacial lake can be calculated using the formula of this invention:

[0190] ;

[0191] The calculated volume of the glacial lake is 59.14 × 10⁻⁶. 6 m 3The actual measured volume of the ice lake was 64 × 10⁻⁶. 6 m 3 The error is only 7.5%, and the calculation results are highly reliable.

[0192] Experimental Example 2

[0193] The area, morphological coefficient, moraine undulation coefficient, topographic coefficient, and watershed morphological index of 24 glacial lakes were collected. The glacial lake volume was calculated using the method of this invention, as shown in Table 1. Table 2 compares the calculation results of this invention with other formulas that only use area to estimate glacial lake volume. Statistical analysis of calculation errors is shown below. Figure 9 As shown.

[0194] Table 1. Glacial lake volume calculated by the method of the present invention

[0195]

[0196] As shown in Table 1, the calculated volume of this invention is close to the measured volume. The calculation result for 1 / 3 of the ice lake has an error of less than 5%, with a maximum error of only -25.95% and a minimum error of -0.12%.

[0197] Table 2 Comparison of calculation results of the present invention with other commonly used formulas

[0198]

[0199] As shown in Table 2, the calculation results of the present invention have an error range of -25.95% to 19.82%, with a median error of 0.11%; the calculation results of the O'Connor formula have an error range of 53.82% to 1147.37%, with a median error of 313.32%; the calculation results of the Huggel formula have an error range of -55.54% to 129.04%, with a median error of -20.53%; and the calculation results of the Evans formula have an error range of -54.58% to 115.38%, with a median error of -20.24%.

[0200] By comparison, the calculation formula of this invention considers multiple parameters that directly affect the volume of glacial lakes, rather than estimating the volume of glacial lakes solely using the area of ​​the glacial lake. The calculation results are significantly better than commonly used calculation formulas such as the O'Connor formula, Huggel formula, and Evans formula, and can provide a basis for accurately assessing the risk of glacial lake outburst.

[0201] The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for calculating the volume of glacial lakes based on morphological and topographical constraints, characterized in that, include: Obtain preprocessed remote sensing image data, and based on the preprocessed remote sensing image data, obtain the area and perimeter of the glacial lake; The morphological coefficient of the glacier lake is obtained based on the area and perimeter of the glacier lake; Collect digital elevation model (DEM) data of the area above the spillway of the terminal moraine of the target glacial lake to obtain the undulation coefficient of the moraine around the glacial lake and the morphological index of the glacial lake basin. Obtain the digital elevation model (DEM) data of the slope range of glacial moraine around the target glacial lake, and then obtain the topographic coefficients around the glacial lake; The target glacial lake volume is obtained based on the glacial lake area, glacial lake morphology coefficient, glacial moraine undulation coefficient around the glacial lake, glacial lake topography coefficient, and glacial lake watershed morphology index. Collect remote sensing image data around the target glacial lake, preprocess the remote sensing image data, and obtain preprocessed remote sensing image data; The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was used. The area of ​​the glacial lake was delineated using polygons, and the area of ​​the polygons was calculated using ArcGIS's computational geometry functions. This area represents the glacial lake area. Calculate the perimeter of the polygon, which is the perimeter of the ice lake. P ; and / or The formula for calculating the topographic coefficient around the glacial lake is as follows: ω = tanα; In the formula: ω α represents the topographic coefficient around the glacial lake, and α is the average slope of the slope range around the glacial lake. The formula for calculating the target glacial lake volume is as follows: ； In the formula: For the target glacial lake volume, The area of ​​the glacial lake, For the morphological coefficient of glacial lakes, The coefficient of variation of glacial moraine around the glacial lake. The topographic coefficient around the glacial lake. This is a morphological index for the glacial lake basin.

2. The method for calculating the volume of glacial lakes based on morphological and topographical constraints according to claim 1, characterized in that, The preprocessing includes using ENVI software to perform geometric correction, orthorectification, and image fusion on the remote sensing image data.

3. The method for calculating the volume of glacial lakes based on morphological and topographical constraints according to claim 1, characterized in that, The formula for calculating the morphological coefficient of a glacial lake is as follows: ； In the formula, For the morphological coefficient of glacial lakes, For the area of ​​the glacial lake, The perimeter of the ice lake; and / or The formula for calculating the undulation coefficient of glacial moraine around the glacial lake is as follows: ； In the formula, The coefficient of variation of glacial moraine around the glacial lake. The elevation difference of the glacial moraine slopes surrounding the glacial lake. This refers to the elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine.

4. The method for calculating the volume of glacial lakes based on morphological and topographical constraints according to claim 3, characterized in that, Using ArcGIS software, the digital elevation model (DEM) data of the glacial lake area was overlaid onto the area above the spillway of the target glacial lake's terminal moraine. The lake surface elevation was then extracted from the glacial lake area. The elevation of the glacial lake terminal moraine overflow outlet was directly read using the digital elevation model (DEM) data of the area above the target glacial lake terminal moraine overflow outlet. The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine is... for: ； In the formula: The elevation difference between the surface of the glacial lake and the overflow outlet of the terminal moraine dike; This refers to the elevation of the ice lake surface. Elevation of the spillway of the glacial lake terminal moraine; The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was employed. Polygons were used to delineate the slopes of the lateral and terminal moraines surrounding the glacial lake, facing the lake. ArcGIS's cropping tool was used to crop the Digital Elevation Model (DEM) data of the area above the spillway of the target glacial lake's terminal moraines to the defined polygon area, thus obtaining the DEM data of the glacial moraine slopes surrounding the target glacial lake. The maximum elevation of the glacial moraine slopes around the target glacial lake was then read from this DEM data. The elevation difference of the glacial moraine slopes around the glacial lake for: ； In the formula: The elevation difference of the glacial moraine slopes surrounding the glacial lake; This represents the maximum elevation of the glacial moraine slope surrounding the glacial lake. This refers to the elevation of the ice lake surface.

5. The method for calculating the volume of glacial lakes based on morphological and topographical constraints according to claim 1, characterized in that, Using ArcGIS software's cropping tool, the digital elevation model (DEM) data of the area above the spillway of the target glacial lake's terminal moraine dam was cropped to the watershed area above the lake surface, obtaining the DEM data of the watershed area above the lake surface. Then, using ArcGIS's contour line tool, contour lines were generated for the watershed area above the lake surface. Finally, ArcGIS software was used to calculate the area of ​​the watershed area above the lake surface. The maximum elevation values ​​of the DEM data for the watershed above the surface of the glacial lake were extracted using ArcGIS software. and the elevation of the ice lake surface The difference between the maximum elevation of the DEM data for the watershed above the glacial lake surface and the elevation of the glacial lake surface is... for: ； Using DEM data of the catchment area above the glacial lake surface, the area 'a' above each contour line and the elevation difference between each contour line and the glacial lake surface elevation are extracted. h Then a i h is the area above the contour line with elevation value i. i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake. make , ; Multiple ( Connecting these curves yields the function curve: ; For functions Integrating the results, we obtain the Glacier Lake Basin Morphological Index HI: ; In the formula, HI For the morphological index of the glacial lake basin; a i Let i be the area above the contour line with elevation value i. The area of ​​the catchment above the surface of the glacial lake; h i The elevation difference between the contour line with elevation value i and the surface elevation of the glacial lake. This represents the elevation difference between the maximum elevation value of the DEM data for the catchment area above the surface of the glacial lake and the elevation of the glacial lake surface.

6. The method for calculating the volume of glacial lakes based on morphological and topographical constraints according to claim 1, characterized in that, Using the digital elevation model (DEM) data of the glacial moraine slope area surrounding the target glacial lake, ArcGIS software was used to calculate the slope data of the glacial moraine slope area surrounding the target glacial lake. Then, ArcGIS software was used to statistically analyze and calculate the slope data of the glacial moraine slope area surrounding the glacial lake to obtain the average slope α of the slope area surrounding the glacial lake.

7. A system for calculating the volume of glacial lakes based on morphological and topographical constraints, characterized in that, include: The module for measuring the area and perimeter of a glacial lake obtains preprocessed remote sensing image data, and based on the preprocessed remote sensing image data, it obtains the area and perimeter of the glacial lake. The ice lake morphology coefficient module obtains the ice lake morphology coefficient based on the ice lake area and ice lake perimeter; The modules for the undulation coefficient of glacial moraine around the glacial lake and the morphological index of the glacial lake basin collect digital elevation model (DEM) data of the area above the spillway of the final moraine dike of the target glacial lake to obtain the undulation coefficient of glacial moraine around the glacial lake and the morphological index of the glacial lake basin. The Glacier Lake Surrounding Topographic Coefficient module obtains the Digital Elevation Model (DEM) data of the glacial moraine slope range surrounding the target glacier lake, and then obtains the glacier lake surrounding topographic coefficients. The target glacial lake volume module obtains the target glacial lake volume based on the glacial lake area, glacial lake morphology coefficient, glacial moraine undulation coefficient around the glacial lake, glacial lake topography coefficient, and glacial lake watershed morphology index. Collect remote sensing image data around the target glacial lake, preprocess the remote sensing image data, and obtain preprocessed remote sensing image data; The preprocessed remote sensing image data was loaded into ArcGIS software. A combination of human-computer interactive interpretation and visual interpretation was used. The area of ​​the glacial lake was delineated using polygons, and the area of ​​the polygons was calculated using ArcGIS's computational geometry functions. This area represents the glacial lake area. Calculate the perimeter of the polygon, which is the perimeter of the ice lake. P ; and / or The formula for calculating the topographic coefficient around the glacial lake is as follows: ω = tanα; In the formula: ω α represents the topographic coefficient around the glacial lake, and α is the average slope of the slope range around the glacial lake. The formula for calculating the target glacial lake volume is as follows: ； In the formula: For the target glacial lake volume, The area of ​​the glacial lake, For the morphological coefficient of glacial lakes, The coefficient of variation of glacial moraine around the glacial lake. The topographic coefficient around the glacial lake. This is a morphological index for the glacial lake basin.

8. A computer device, characterized in that, It includes a storage device and a processor, the storage device storing a computer program that, when executed by the processor, causes the processor to perform the steps of the method as described in any one of claims 1-6.

9. A computer-readable storage medium, characterized in that, The device stores a computer program that, when executed by a processor, causes the processor to perform the steps of the method as described in any one of claims 1-6.

Images