Ice lake terrain estimation method, system, device, and storage medium

Abstract

The application discloses an ice lake terrain estimation method, system, device and storage medium, and belongs to the technical field of ice lake terrain measurement and estimation. The method comprises the following steps: acquiring basic data of an ice lake; extracting morphological parameters from the basic data; obtaining a volume-area fitting formula of the ice lake; obtaining a theoretical volume of the ice lake; simultaneously solving the volume-area fitting formula of the ice lake and the theoretical volume of the ice lake to obtain a scaling factor and a relationship between the maximum depth of the ice lake and the area of the ice lake; generating ice lake contour lines; and generating an ice lake terrain. The application simplifies the complex modeling process of the ice lake terrain, improves the terrain estimation efficiency under the condition of retaining high precision, and provides technical support for regional ice lake risk management and monitoring system construction. The characteristics of the ice lake, such as the long and short axes, the central coordinates, the volume-area relationship and the like, are fully considered, which is helpful to reveal the spatial distribution law of the ice lake terrain and reflect the terrain details, thereby providing support for the regional ice lake risk management and monitoring system construction.

Description

Technical Field

[0001] This invention belongs to the field of glacial lake topography measurement and estimation technology, and particularly relates to glacial lake topography estimation methods, systems, equipment and storage media. Background Technology

[0002] Glacial lake topography is crucial foundational data in glaciology, hydrology, water resources, and disaster prevention and mitigation. Accurate estimation of glacial lake topography is essential for assessing lake stability, monitoring water volume changes, and providing early warnings of outburst risk. Traditional methods for estimating glacial lake topography primarily rely on field surveys or inversion from single remote sensing data. Field surveys are limited by terrain conditions and struggle to cover remote or hazardous glacial lake areas; methods based on single remote sensing data often ignore the nonlinear characteristics of glacial lake morphology, leading to insufficient accuracy in topography estimation, particularly in the simulation of isobath distribution and depth gradients. Summary of the Invention

[0003] The purpose of this invention is to provide a method, system, device, and storage medium for estimating the topography of glacial lakes, so as to solve the problems of insufficient accuracy and difficulty in reflecting topographic details in the prior art.

[0004] The embodiments of this application implement the glacial lake topography estimation method as follows:

[0005] Obtain basic data on the glacial lake;

[0006] Extract morphological parameters from the basic data;

[0007] Obtain the fitted formula for the volume-area of ​​the glacial lake;

[0008] To obtain the theoretical volume of the glacial lake;

[0009] By combining the volume-area fitting formula of the glacial lake with the theoretical volume of the glacial lake, we can obtain the scaling factor and the relationship between the maximum depth of the glacial lake and its area.

[0010] Based on the basic data, morphological parameters, scaling factors, and the relationship between the maximum depth and area of ​​the glacial lake, isobaths of the glacial lake are generated.

[0011] Glacial lake topography is generated based on the isobaths of the glacial lake.

[0012] Optionally, in some embodiments of this application, the basic data includes the surface area, volume, maximum depth, and outline polygon of the glacial lake; and / or

[0013] The fitted formula for the volume-area of ​​a glacial lake is:

[0014] ;

[0015] In the formula, The actual volume of the glacial lake, in units of: ;

[0016] is the area of the ice lake, unit: ;

[0017] and is the fitting coefficient, unit: dimensionless; and / or

[0018] The theoretical volume of the ice lake is obtained based on the similar scaling assumption of the contour line, and the formula of the theoretical volume of the ice lake is:

[0019] ;

[0020] In the formula, is the theoretical volume of the ice lake, unit: ;

[0021] is the area of the ice lake, unit: ;

[0022] is the depth of the ice lake, unit: ;

[0023] is the maximum depth of the ice lake, unit: ;

[0024] is the scaling factor, unit: dimensionless.

[0025] Optionally, in some embodiments of the present application, the morphological parameters are extracted from the contour polygon of the ice lake, including the length of the major axis, the length of the minor axis, the center coordinates and the direction angle of the major axis of the ice lake.

[0026] Optionally, in some embodiments of the present application, the method for extracting the morphological parameters comprises:

[0027] Obtaining the minimum rotating rectangle of the contour polygon of the ice lake;

[0028] Calculating the length of each side of the minimum rotating rectangle, defining the long side of the rectangle as the major axis and the short side of the rectangle as the minor axis, to obtain the length of the major axis and the length of the minor axis respectively;

[0029] Determining the vector direction of the major axis and calculating the direction angle of the major axis;

[0030] Taking the center of mass of the minimum rotating rectangle as the center coordinates of the ice lake.

[0031] Optionally, in some embodiments of the present application, the fitting formula of the volume-area of the ice lake is combined with the formula of the theoretical volume of the ice lake , let = The scaling factor and the relationship between the maximum depth of the ice lake and the area of the ice lake can be obtained, and the formula is as follows:

[0032] ;

[0033] In the formula, is the scaling factor, unit: dimensionless;

[0034] and are fitting coefficients, unit: dimensionless;

[0035] is the area of the ice lake, unit: ;

[0036] is the maximum depth of the ice lake, unit: .

[0037] Optionally, in some embodiments of the present application, the method for generating the ice lake contour line comprises:

[0038] obtaining the ice lake contour line interval , determining the contour line level according to the maximum depth value of the ice lake and the contour line interval , i.e. calculating the total number of contour lines , the contour polygon is the 0th layer, and the center coordinate is the n th layer;

[0039] scaling the contour polygon to the center coordinate by the scaling factor p to generate the ice lake contour line of each level.

[0040] Optionally, in some embodiments of the present application, the method for generating the ice lake terrain comprises:

[0041] determining the grid coverage range based on the ice lake contour line to generate uniformly distributed plane grid points;

[0042] determining whether each grid point is located inside the contour polygon, and the value of the grid point not inside the contour polygon is 0;

[0043] if the grid point is inside the contour polygon, according to the contour line level i where the grid point is located, the grid point is given a corresponding depth value , until all grid points are given a depth value, i.e. the initial terrain is formed;

[0044] performing Gaussian filtering on the initial terrain to generate the ice lake terrain.

[0045] Correspondingly, the embodiment of the present application also provides an ice lake terrain estimation system, comprising:

[0046] a basic data module for obtaining basic data of an ice lake;

[0047] a shape parameter module for extracting shape parameters from the basic data;

[0048] a volume-area fitting module for obtaining a volume-area fitting formula of the ice lake;

[0049] a theoretical volume module for obtaining a theoretical volume of the ice lake;

[0050] a simultaneous module for simultaneously fitting the volume-area fitting formula of the ice lake and the theoretical volume of the ice lake to obtain a scaling factor and a relationship between a maximum depth of the ice lake and an area of the ice lake;

[0051] an ice lake contour line module for generating an ice lake contour line based on the basic data, the shape parameters, and the scaling factor and the relationship between the maximum depth of the ice lake and the area of the ice lake;

[0052] an ice lake terrain module for generating an ice lake terrain based on the ice lake contour line.

[0053] Correspondingly, the embodiment of the present application also provides a computer device comprising a storage and a processor, wherein the storage stores a computer program, and the computer program is executed by the processor to make the processor execute the steps of the above method.

[0054] Correspondingly, the embodiment of the present application also provides a computer readable storage medium storing a computer program, and the computer program is executed by a processor to make the processor execute the steps of the above method.

[0055] In summary, due to the adoption of the above technical solutions, the present application has the following advantages:

[0056] Based on the lake area, volume and contour data of the ice lake, the present application adopts the similar scaling assumption of the contour line, combines the empirical formula fitting of the volume and area, and estimates the ice lake terrain through the steps of extracting the geometric parameters of the ice lake, simultaneously determining the scaling factor, generating the contour line and terrain, etc. By simplifying the complex modeling process of the ice lake terrain, the present application improves the efficiency of terrain estimation while maintaining high accuracy, and provides technical support for the construction of regional ice lake risk management and monitoring system. The present application fully considers the characteristics of the ice lake such as the major and minor axes, the center coordinates, the volume and area relationship, etc., which helps to reveal the spatial distribution law of the ice lake terrain and reflects the terrain details, thereby providing technical support for the construction of regional ice lake risk management and monitoring system.

[0057] In the real world, the isobath from the shore to the lake center of the ice lake often has certain geometric similarity. The ice lake terrain is estimated by the assumption of isobath similarity scaling in the application, which can make the result closer to the true situation, avoid unreasonable depth distribution and terrain characteristics, and improve the stability and reliability of the estimation.

[0058] The application combines experience fitting and geometric similarity principle, and can effectively estimate the terrain characteristics of the ice lake under different conditions through systematic steps calculation, thereby providing scientific basis for ice lake change monitoring, water resource assessment and disaster prevention and mitigation.

[0059] According to the ice lake terrain estimation method proposed in the application, the scaling factor, depth-area relationship, isobath and terrain data and other key parameters of the ice lake can be obtained, and the corresponding department can make good ice lake related resource management and disaster response measures in advance according to the estimation result. BRIEF DESCRIPTION OF DRAWINGS

[0060] Figure 1 The flow chart of the ice lake terrain estimation method of the application is shown in the figure;

[0061] Figure 2 The ice lake volume-area empirical formula fitting of the application example is shown in the figure;

[0062] Figure 3 The ice lake filled isobath chart of the application example is shown in the figure;

[0063] Figure 4 The 3D lake bottom terrain chart of the application example is shown in the figure;

[0064] Figure 5 The long axis profile chart of the application example is shown in the figure;

[0065] Figure 6 The short axis profile chart of the application example is shown in the figure. DETAILED DESCRIPTION

[0066] In order to make the purpose, technical scheme and advantages of the application more clear and understandable, the application will be further described in detail below in combination with the drawings and examples. It should be understood that the specific examples described herein are only used to explain the application, and are not used to limit the application.

[0067] The technical scheme of the application is as follows:

[0068] Please refer to Figure 1 The ice lake terrain estimation method provided by the application embodiment comprises the following steps:

[0069] S01, obtaining the basic data of the ice lake;

[0070] S02, extracting the morphological parameters from the basic data;

[0071] S03, obtaining the fitting formula of volume-area of ice lake;

[0072] S04, obtaining the theoretical volume of ice lake;

[0073] S05, obtaining the scaling factor and the relationship between the maximum depth of ice lake and the area of ice lake by combining the fitting formula of volume-area of ice lake and the theoretical volume of ice lake;

[0074] S06, generating the ice lake contour based on the basic data, the morphological parameters, and the scaling factor and the relationship between the maximum depth of ice lake and the area of ice lake;

[0075] S07, generating the ice lake terrain based on the ice lake contour.

[0076] In the S01,

[0077] In some embodiments, the basic data includes the surface area, volume, maximum depth and contour polygon of the ice lake.

[0078] In the S02,

[0079] Further, the morphological parameters are extracted from the contour polygon of the ice lake, including the length of the major axis, the length of the minor axis, the center coordinates and the direction angle of the major axis of the ice lake.

[0080] Further, the method for extracting the morphological parameters includes:

[0081] S021, obtaining the minimum rotating rectangle of the contour polygon of the ice lake;

[0082] S022, calculating the length of each side of the minimum rotating rectangle, defining the long side of the rectangle as the major axis and the short side of the rectangle as the minor axis to obtain the length of the major axis and the length of the minor axis respectively;

[0083] S023, determining the vector direction of the major axis and calculating the direction angle of the major axis;

[0084] S024, taking the center of mass of the minimum rotating rectangle as the center coordinates of the ice lake.

[0085] In the S03,

[0086] It can be understood that by plotting the scatter plot of the volume and area of the ice lake, it is found that the volume and area of the ice lake conform to the power function relationship, so the fitting formula of the volume-area of the ice lake adopts the power function.

[0087] Further, the fitting formula of the volume-area of the ice lake is:

[0088] ;

[0089] In the formula, is the actual volume of the ice lake, unit: ;

[0090] The area of ​​the glacial lake, in units of: ;

[0091] and The fitting coefficients are dimensionless.

[0092] Understandable. and It was obtained by least squares fitting.

[0093] In S04:

[0094] Furthermore, based on the similarity scaling assumption of isobaths, the theoretical volume of the glacial lake is obtained. The formula for the theoretical volume of the glacial lake is:

[0095] ;

[0096] In the formula, The theoretical volume of the glacial lake, unit: ;

[0097] The area of ​​the glacial lake, in units of: ;

[0098] The depth of the glacial lake, in units of: ;

[0099] The maximum depth of the glacial lake, in units of: ;

[0100] Scaling factor, unit: dimensionless.

[0101] In S05:

[0102] Furthermore, the fitting formula for the volume-area of ​​the glacial lake is established. Formula for the theoretical volume of a glacial lake ,make = This allows us to obtain the scaling factor and the relationship between the maximum depth of the glacial lake and its area, as shown in the following formula:

[0103] ;

[0104] In the formula, Scaling factor, unit: dimensionless;

[0105] and These are the fitting coefficients, in dimensionless form.

[0106] is the area of the ice lake, unit: ;

[0107] is the maximum depth of the ice lake, unit: .

[0108] In the S06:

[0109] Further, the method for generating the ice lake contour line comprises:

[0110] S061, obtaining the contour line interval of the ice lake , according to the maximum depth value of the ice lake and the contour line interval , determine the contour line level, that is, calculate the total number of contour lines , the contour polygon is the 0th layer, and the center coordinate is the n layer;

[0111] S062, the contour polygon is proportionally scaled to the center coordinate according to the scaling factor p , to generate the ice lake contour line of each level.

[0112] Exemplarily, the contour line interval is 0.1m.

[0113] It can be understood that the contour polygon is the 0th layer of the contour line, and the depth is 0. For the i layer, the depth is , i 0, 1, 2,..., n; to the nth layer, the depth is .

[0114] It can be understood that for the i layer contour line, the depth is ( ).

[0115] It can be understood that and are known data, and the scaling factor can be directly calculated. .

[0116] In the S07:

[0117] Further, the method for generating the ice lake terrain comprises:

[0118] S071, determining the grid coverage range based on the ice lake contour line, and generating uniformly distributed plane grid points;

[0119] S072, determining whether each grid point is located inside the contour polygon, the value of the grid point not inside the contour polygon is 0;

[0120] S073, if the grid point is inside the contour polygon, according to the depth contour level where the grid point is located i , the corresponding depth value of the grid point is given , until all grid points are given depth values, that is, the initial terrain is formed;

[0121] S074, Gaussian filtering is performed on the initial terrain to generate the ice lake terrain.

[0122] It can be understood that Gaussian filtering is performed on the initial terrain to make the generated terrain smoother.

[0123] It can be understood that the grid coverage range is determined based on the polygon boundary of the ice lake depth contour.

[0124] In a second aspect, embodiments of the present application provide an ice lake terrain estimation system, comprising:

[0125] a basic data module for obtaining basic data of an ice lake;

[0126] a shape parameter module for extracting a shape parameter from the basic data;

[0127] a volume-area fitting module for obtaining a volume-area fitting formula of the ice lake;

[0128] a theoretical volume module for obtaining a theoretical volume of the ice lake;

[0129] a simultaneous module for simultaneously solving the volume-area fitting formula of the ice lake and the theoretical volume of the ice lake to obtain a scaling factor and a relationship between the maximum depth of the ice lake and the area of the ice lake;

[0130] an ice lake depth contour module for generating an ice lake depth contour based on the basic data, the shape parameter, and the relationship between the scaling factor and the maximum depth of the ice lake and the area of the ice lake;

[0131] an ice lake terrain module for generating an ice lake terrain based on the ice lake depth contour.

[0132] In the basic data module:

[0133] In some embodiments, the basic data includes the lake surface area, volume, maximum depth and contour polygon of the ice lake.

[0134] In the shape parameter module;

[0135] Further, the shape parameter is extracted from the contour polygon of the ice lake, and the shape parameter includes the length of the major axis, the length of the minor axis, the center coordinates and the direction angle of the major axis of the ice lake.

[0136] Further, the method for extracting the shape parameters comprises:

[0137] S021, obtaining a minimum rotating rectangle of the ice lake contour polygon;

[0138] S022, calculating the length of each side of the minimum rotating rectangle, defining the long side of the rectangle as the long axis and the short side of the rectangle as the short axis, and obtaining the length of the long axis and the length of the short axis respectively;

[0139] S023, determining the vector direction of the long axis and calculating the direction angle of the long axis;

[0140] S024, taking the center of mass of the minimum rotating rectangle as the center coordinates of the ice lake.

[0141] In the volume-area fitting module:

[0142] It can be understood that by plotting the scatter plot of the volume and area of the ice lake, it is found that the volume and area of the ice lake conform to the power function relationship, so the volume-area fitting formula of the ice lake adopts the power function.

[0143] Further, the volume-area fitting formula of the ice lake is:

[0144] ;

[0145] In the formula, is the actual volume of the ice lake, unit: ;

[0146] is the area of the ice lake, unit: ;

[0147] and is the fitting coefficient, unit: dimensionless.

[0148] It can be understood that and are obtained by least square fitting.

[0149] In the theoretical volume module:

[0150] Further, the theoretical volume of the ice lake is obtained based on the similar scaling assumption of the contour line, and the theoretical volume formula of the ice lake is:

[0151] ;

[0152] In the formula, is the theoretical volume of the ice lake, unit: ;

[0153] is the area of the ice lake, unit: ;

[0154] is the depth of the ice lake, unit: ;

[0155] is the maximum depth of the ice lake, unit: ;

[0156] is the scaling factor, unit: dimensionless.

[0157] In the simultaneous module:

[0158] Further, the fitting formula of the volume-area of the ice lake and the theoretical volume formula of the ice lake , let = The scaling factor and the relationship between the maximum depth of the ice lake and the area of the ice lake can be obtained, and the formula is as follows:

[0159] ;

[0160] In the formula, is the scaling factor, unit: dimensionless;

[0161] and are fitting coefficients, unit: dimensionless;

[0162] is the area of the ice lake, unit: ;

[0163] is the maximum depth of the ice lake, unit: .

[0164] In the ice lake contour module:

[0165] Further, the method for generating ice lake contours comprises:

[0166] S061, obtain the ice lake contour interval , according to the maximum depth value of the ice lake and the contour interval , determine the contour level, that is, calculate the total number of contour layers , the contour polygon is the 0th layer, and the center coordinates are the n th layer;

[0167] S062, scale the contour polygon to the center coordinates according to the scaling factor p , to generate the ice lake contours of each level.

[0168] Exemplarily, the contour polygon is the 0th layer of the contour line, the depth is 0, and for the 1st layer, the depth is .

[0169] It can be understood that the contour polygon is the 0th layer of the contour line, the depth is 0, and for the 1st layer, the depth is i . , i is 0, 1, 2,..., n; and for the nth layer, the depth is .

[0170] It can be understood that for the 1st layer of the contour line, the depth is i . .

[0171] It can be understood that and are known data, and the scaling factor can be directly calculated. .

[0172] In the ice lake terrain module, the following steps are performed:

[0173] Further, the method for generating the ice lake terrain comprises the following steps:

[0174] S071, determining the grid coverage range based on the contour line of the ice lake, and generating the uniformly distributed plane grid points;

[0175] S072, determining whether each grid point is located inside the contour polygon, and the value of the grid point not inside the contour polygon is 0;

[0176] S073, if the grid point is inside the contour polygon, according to the layer level of the contour line where the grid point is located i , the corresponding depth value is given to the grid point, until all the grid points are given the depth value, that is, the initial terrain is formed;

[0177] S074, performing Gaussian filtering on the initial terrain to generate the ice lake terrain.

[0178] It can be understood that the Gaussian filtering on the initial terrain makes the generated terrain smoother.

[0179] It can be understood that the grid coverage range is determined based on the polygon boundary of the contour line of the ice lake.

[0180] In a third aspect, the present application provides a computer device, comprising a storage and a processor, wherein the storage stores a computer program, and the computer program is executed by the processor to make the processor perform the steps of the ice lake terrain estimation method as described above.

[0181] ​The computer device can be a desktop computer, a notebook computer, a palm computer, a cloud server, or the like. The computer device can interact with a user through a keyboard, a mouse, a remote controller, a touchpad, a voice control device, or the like.

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

[0183] The processor can be a central processing unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chips in some embodiments. The processor is usually used to control the overall operation of the computer device. In this embodiment, the processor is used to run program codes or process data stored in the memory, such as program codes of the ice lake terrain estimation method.

[0184] In a fourth aspect, the present application provides a computer readable storage medium storing a computer program. When the computer program is executed by a processor, the processor performs the steps of the ice lake terrain estimation method as described above.

[0185] The computer readable storage medium stores an interface display program. The interface display program can be executed by at least one processor to make the at least one processor perform the steps of the ice lake terrain estimation method as described above.

[0186] Through the description of the above embodiments, those skilled in the art can clearly understand that the above-mentioned embodiment method can be realized by means of software and the necessary general hardware platform, of course, it can also be realized by hardware, but in many cases, the former is a better embodiment. Based on such understanding, the technical solutions of the present application can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, magnetic disk, optical disk), and includes a plurality of instructions for making a terminal device (which can be a mobile phone, computer, server or network device, etc.) execute the ice lake terrain estimation method described in the embodiments of the present application.

[0187] The application will be further described below in combination with application examples.

[0188] Application Example

[0189] In a certain domestic ice lake, in this application example, the ice lake data is obtained through remote sensing and actual measurement, including lake contour, maximum depth (8.07 m), area (7195.27 m²), and actual volume (27100 m³); when extracting parameters, the grid resolution is set to 150x150, and the contour interval is 0.1 m; when generating the contour, the lake geometric center is taken as the scaling origin, and the volume-area empirical formula (a=0.0445, b=1.5) can be obtained by fitting the existing ice lake area and ice lake volume data in this area, and the fitting image is as shown in Figure 2 The scaling factor p is calculated to be 0.8781 based on the fitted volume-area empirical formula and the contour similarity scaling assumption; in the terrain generation, the Gaussian filter with sigma=1 is used to eliminate discrete errors. Under the above parameter conditions, the estimation of the ice lake terrain is successfully realized, and the terrain estimation result is as shown in Figures 3 to 6 Figure 3 The filled contour map is shown in blue gradient to intuitively present the depth distribution (light blue represents shallow water area, and dark blue represents deep water area), and the red solid line marks the long axis, and the orange dashed line marks the short axis; Figure 4 The 3D lake bottom terrain map clearly shows the undulating shape of the lake bottom; Figure 5 The long axis profile and Figure 6 The short axis profile respectively reflect the depth variation law along the axis direction, and the lake surface height is set to 0 m in the figure, and the depth is a negative value indicating the corresponding vertical downward distance from the lake surface.

[0190] Through calculation, the contour volume estimation value is 26798.34 m³, and the relative error with the actual volume is 1.11 %, which verifies the accuracy of the ice lake terrain estimation method of the present application.

[0191] ​The above merely describes preferred embodiments of the present application, and is not used to limit the present application, any modification, equivalent replacement and improvement within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims

1. A method of estimating ice lake terrain, characterized by, The method comprises: obtaining basic data of the ice lake; extracting morphological parameters from the basic data; obtaining a volume-area fitting formula of the ice lake; obtaining a theoretical volume of the ice lake; obtaining a scaling factor and a relationship between the maximum depth of the ice lake and the area of the ice lake by simultaneously solving the volume-area fitting formula of the ice lake and the theoretical volume of the ice lake; generating ice lake contour lines based on the basic data, the morphological parameters, and the scaling factor and the relationship between the maximum depth of the ice lake and the area of the ice lake; generating an ice lake terrain based on the ice lake contour lines; The simultaneous fitting of the volume-area formula of the ice lake The theoretical volume formula of the ice lake , let = The scaling factor and the relationship between the maximum depth of the ice lake and the area of the ice lake can be obtained, and the formula is as follows: ； wherein is a scaling factor, dimensionless; and = fitting coefficient, dimensionless; Vice, in units of: ; Theoretical volume of the ice lake, units: ; Area of ice lake, unit: ; Depth of the ice lake, units: ; Maximum depth of the ice lake, unit: .

2. The ice-lake terrain estimation method according to claim 1, characterized by, the basic data comprises the surface area, volume, maximum depth, and contour polygon of the ice lake; and / or the volume-area fitting formula of the ice lake is: ； wherein Vice, in m3, of the ice lake ; Area of ice lake, unit: ; and are fitting coefficients, dimensionless; and / or the theoretical volume of the ice lake is obtained based on the similar scaling assumption of the contour lines, and the formula of the theoretical volume of the ice lake is: ； wherein Vice, in m3, of the ice lake ; Area of ice lake, unit: ; Depth of the ice lake, units: ; is the maximum depth of the ice lake, in units of ; Scaling factor, dimensionless.

3. The ice-lake terrain estimation method according to claim 2, characterized by, the morphological parameters are extracted from the contour polygon of the ice lake, and the morphological parameters comprise the length of the major axis, the length of the minor axis, the center coordinates, and the direction angle of the major axis of the ice lake.

4. The ice-lake terrain estimation method according to claim 3, characterized by, The method for extracting the morphological parameters comprises: obtaining a minimum rotating rectangle of the contour polygon of the ice lake; calculating the length of each side of the minimum rotating rectangle, defining the long side of the rectangle as the major axis and the short side of the rectangle as the minor axis to obtain the length of the major axis and the length of the minor axis respectively; determining the vector direction of the major axis and calculating the direction angle of the major axis; taking the center of mass of the minimum rotating rectangle as the center coordinates of the ice lake.

5. The ice-lake terrain estimation method according to claim 1, characterized by, The method for generating the ice lake contour lines comprises: obtaining the isobath interval of the ice lake , determining the isobath level according to the maximum depth value of the ice lake and the isobath interval , that is, calculating the total number of isobaths , the contour polygon is the 0th layer, and the center coordinate is the 1st layer n ; According to the scaling factor p The contour polygons are scaled to the center coordinates proportionally to generate the ice lake isobaths of each level.

6. The ice lake terrain estimation method according to claim 5, characterized by, The method for generating the ice lake terrain comprises: determining the grid coverage range based on the ice lake contour lines to generate uniformly distributed plane grid points; determining whether each grid point is located inside the contour polygon, and the value of the grid point not inside the contour polygon is 0; If the grid point is inside the outline polygon, it depends on the contour level where the grid point is located. i Assign the depth value corresponding to the grid point. This continues until all grid points have been assigned depth values, meaning the initial terrain has been formed. performing Gaussian filtering on the initial terrain to generate the ice lake terrain.

7. A system for estimating ice lake terrain, characterized by The method comprises: a basic data module for obtaining basic data of the ice lake; a morphological parameter module for extracting morphological parameters from the basic data; a volume-area fitting module for obtaining a volume-area fitting formula of the ice lake; a theoretical volume module for obtaining a theoretical volume of the ice lake; a simultaneous solving module for obtaining a scaling factor and a relationship between the maximum depth of the ice lake and the area of the ice lake by simultaneously solving the volume-area fitting formula of the ice lake and the theoretical volume of the ice lake; an ice lake contour line module for generating ice lake contour lines based on the basic data, the morphological parameters, and the scaling factor and the relationship between the maximum depth of the ice lake and the area of the ice lake; an ice lake terrain module for generating an ice lake terrain based on the ice lake contour lines; The simultaneous fitting of the volume-area formula of the ice lake The theoretical volume formula of the ice lake Let = The scaling factor and the relationship between the maximum depth of the ice lake and the area of the ice lake can be obtained, and the formula is as follows: ； wherein is a scaling factor, dimensionless; and = fitting coefficient, dimensionless; Vice, in units of: ; Theoretical volume of the ice lake, units: ; Area of ice lake, unit: ; Depth of ice lake, units: ; Depth of ice lake, unit: .

8. Computer device, characterized in that a computer program is stored in a storage and executed by a processor, and the computer program causes the processor to execute the steps of the method according to any one of claims 1-6 when executed by the processor.

9. A computer readable storage medium, characterized in that, A computer program is stored in a storage and executed by a processor, and the computer program causes the processor to execute the steps of the method according to any one of claims 1-6 when executed by the processor.

Images