Methods, equipment, and storage media for vertical mesh generation in complex terrain

Abstract

This application discloses a method, device, and storage medium for generating vertical grids suitable for complex terrain, relating to the field of fluid numerical computation technology. The key technical points are: the method includes: acquiring topographic elevation data of the target terrain region; generating planar triangular or quadrilateral grids based on the topographic elevation data; defining the control water depth and grid layer number for several terrain feature control points; determining the critical water depth between deep and shallow water areas; calculating the characteristic values ​​and elevation values ​​of each grid point; performing local vertical subtraction to generate a series of characteristic values ​​and elevation values; and generating the vertical grid of the target terrain region based on the series of characteristic values ​​and elevation values ​​of each grid point. This generation method can adapt to complex underwater terrain and allocate appropriate vertical grid layers in deep and shallow water areas respectively, satisfying the requirements of vertical resolution while improving model computational efficiency and accuracy.

Description

Technical Field

[0001] This invention relates to the field of fluid numerical computation technology, specifically to a method, device, and storage medium for generating vertical grids suitable for complex terrain. Background Technology

[0002] This section is intended to provide background or context for the embodiments set forth in the claims. The description herein is not an admission that it is prior art simply because it is included in this section.

[0003] The construction of the model mesh includes the construction of planar and vertical mesh scales. The determination of the planar mesh mainly refers to the terrain plane, while the division of the vertical mesh is the basis for characterizing the longitudinal heterogeneity of terrain features and needs to fully consider the variation law of water depth.

[0004] In existing technologies, three common vertical coordinate systems used in 3D fluid numerical computation are z-coordinate, S-coordinate, and hybrid SZ-coordinate. The z-coordinate is simple, with each grid layer parallel to the others in the vertical direction, eliminating the need for secondary processing of the vertical coordinates. However, it can produce stepped terrain in complex terrain conditions. The S-coordinate, on the other hand, has grid layers that are not parallel to each other in the vertical direction, and the surfaces containing each grid layer are often curved. The S-coordinate fits terrain and free water surfaces well. However, when there are significant variations in water depth, it can lead to overly dense vertical layering at shallower depths and insufficient vertical coordinate resolution at deeper depths. Furthermore, for more complex underwater terrain, such as steep slopes, deep pits, and submerged breakwaters, all three vertical coordinate systems have limitations. Summary of the Invention

[0005] To address the aforementioned technical problems, this invention proposes a method, device, and storage medium for generating vertical meshes suitable for complex terrains. This method can adapt to complex underwater terrains and allocate appropriate vertical mesh layers in deep and shallow water areas, thus meeting the requirements for vertical resolution while improving model computation efficiency and accuracy.

[0006] To solve the above-mentioned technical problems, the technical solution adopted by the present invention includes three aspects.

[0007] Firstly, a method for generating vertical meshes suitable for complex terrain is provided, comprising the following steps:

[0008] Obtain the topographic elevation data of the target terrain area;

[0009] A planar triangular or quadrilateral grid is generated based on the terrain elevation data;

[0010] Based on the terrain elevation data, define the control water depth and the number of grid layers for several terrain feature control points;

[0011] Determine the critical water depth at the transition point between deep and shallow water areas in the target topographic region;

[0012] Calculate the characteristic value and elevation value of each grid point based on the actual water depth and critical water depth of each grid point within the target terrain area;

[0013] Based on the actual water depth of each grid point and the two adjacent control water depths close to that actual water depth, a series of characteristic values ​​and elevation values ​​are generated for each grid point by performing local vertical difference calculations.

[0014] The vertical grid of the target terrain region is generated based on a series of feature values ​​and elevation values ​​of each grid point in the target terrain region.

[0015] In some embodiments, calculating the characteristic value and elevation value of each grid point based on the actual water depth and critical water depth of each grid point within the target terrain area includes:

[0016] Compare the actual water depth at each grid point with the critical water depth.

[0017] When the actual water depth at a grid point is greater than or equal to the critical water depth, the characteristic value and elevation value at that grid point are calculated using the first set of preset formulas.

[0018] When the actual water depth at a grid point is less than the critical water depth, the characteristic value and elevation value at that grid point are calculated using the preset second set of formulas.

[0019] In some embodiments, the first set of formulas includes a first formula for calculating elevation values ​​and a second formula for calculating characteristic values; the first formula is:

[0020]

[0021] In the formula, σ is the eigenvalue. h c Where η is the critical water depth, h is the actual water depth, and η is the reference water level.

[0022] The second formula is:

[0023]

[0024] In the formula, σ is the eigenvalue, and θ b and θ f Define and set coefficients to control the density of the target terrain area at the bottom and water surface in the grid;

[0025] When θ f When the coordinates approach zero, the vertical coordinate becomes the characteristic value, and the elevation value is calculated using the following formula:

[0026]

[0027] In the formula, h c η is the critical water depth, h is the actual water depth, η is the reference water level, and σ is the characteristic value.

[0028] In some embodiments, the second set of formulas includes a third formula for calculating elevation values ​​and a fourth formula for calculating characteristic values; the third formula is:

[0029] z k =h i σk=1,2,3......,N i ;

[0030] In the formula, h i The control water depth for the topographic feature control point of this grid point is σ, where σ is the characteristic value.

[0031] The fourth formula is: σ=ασ 2 +(1+α)σ;

[0032] In the formula, the magnitude of σ is evenly distributed between -1 and 0, where -1 represents the bottom and 0 represents the water surface; α is a defined set coefficient.

[0033] In some embodiments, generating a vertical grid for the target terrain region based on a series of feature values ​​and elevation values ​​of each grid point in the target terrain region includes:

[0034] For grid points in the vertical grid where the actual water depth is greater than or equal to the critical water depth, when its θ f When the value is greater than or equal to 1, the vertical grid at the boundary of the grid point is refined;

[0035] For grid points in the vertical grid where the actual water depth is greater than or equal to the critical water depth, when its θ b When θ approaches 0, the vertical grid at the point near the water surface is refined; when its θ... b When the value approaches 1, the vertical grid on both the surface and bottom of the grid point is refined.

[0036] In some embodiments, generating a vertical grid for the target terrain region based on a series of feature values ​​and elevation values ​​of each grid point in the target terrain region includes:

[0037] For grid points in the vertical grid where the actual water depth is less than the critical water depth, when α = 0, the traditional eigenvalue is obtained through the second set of formulas; when α < 0, the vertical grid near the water surface of the grid point is densified in layers; when α > 0, the vertical grid near the bottom of the grid point is densified in layers.

[0038] In some embodiments, the step of generating a series of feature values ​​and elevation values ​​for each grid point by performing local vertical difference calculations based on the actual water depth of each grid point and two adjacent control water depths close to the actual water depth includes:

[0039] When the thickness of the bottommost grid layer of the interpolated grid point is less than the preset minimum grid layer thickness, the grid layer is merged with the grid layer above it to preserve the actual water depth of the grid point.

[0040] In some embodiments, generating a vertical grid for the target terrain region based on a series of feature values ​​and elevation values ​​of each grid point in the target terrain region includes:

[0041] When adjacent grid points have different vertical grid layers, the vertical grid at the bottom junction is represented by trapezoidal elements.

[0042] In a second aspect, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the aforementioned vertical mesh generation method.

[0043] Thirdly, a computer-readable storage medium is provided, wherein when the computer program is executed by a processor, it implements the steps of the aforementioned vertical mesh generation method.

[0044] Compared with the prior art, one or more embodiments of the above solutions may have the following advantages or beneficial effects:

[0045] This application provides a method, device, and storage medium for generating vertical meshes suitable for complex terrain. It can define the control water depth and the number of mesh layers according to terrain changes to represent terrain features and ensure that the transition of the vertical mesh is smooth enough when the water depth difference between adjacent nodes is large. Each mesh point on the plane can have different numbers of layers in the vertical direction according to different water depths, which can adapt to complex underwater terrain. At the same time, appropriate vertical mesh layers are allocated in deep water areas and shallow water areas respectively, which not only meets the requirements of vertical resolution, but also improves the model calculation efficiency and accuracy. Attached Figure Description

[0046] The present application will be described in more detail below based on embodiments and with reference to the accompanying drawings;

[0047] Figure 1 This is a flowchart illustrating a vertical mesh generation method suitable for complex terrain provided in Embodiment 1 of the present invention.

[0048] Figure 2 In the embodiments of the present invention, corresponding to Figure 1 An exemplary flowchart of step S5 shown;

[0049] Figure 3 In the embodiments of the present invention, corresponding to Figure 1 An exemplary flowchart of step S7 shown;

[0050] Figure 4 This is an example diagram of the vertical mesh generated in Embodiment 1 of the present invention;

[0051] Figure 5 This is an example diagram of the vertical mesh generated in Embodiment 2 of the present invention;

[0052] Figure 6 This is an example diagram of the vertical mesh generated in Embodiment 3 of the present invention;

[0053] Figure 7 This is an example diagram of the vertical mesh generated in Embodiment 4 of the present invention;

[0054] Figure 8 This is a schematic block diagram of an electronic device provided in Embodiment 5 of the present invention;

[0055] Figure 9 This is a schematic diagram of a computer-readable storage medium provided in Embodiment Six of the present invention.

[0056] In the accompanying drawings, the same parts are referred to by the same reference numerals, and the drawings are not drawn to scale. Detailed Implementation

[0057] The present disclosure will be further described below with reference to the embodiments shown in the accompanying drawings.

[0058] Example 1:

[0059] This application discloses a method for generating vertical meshes suitable for complex terrain, such as... Figure 1As shown, the method includes the following steps: acquiring topographic elevation data of the target terrain region; generating a planar triangular or quadrilateral grid based on the topographic elevation data; defining the control water depth and grid layer number of several terrain feature control points based on the topographic elevation data; determining the critical water depth at the transition point between deep and shallow water areas in the target terrain region; calculating the characteristic value and elevation value of each grid point based on the actual water depth and critical water depth of each grid point in the target terrain region; generating a series of characteristic values ​​and elevation values ​​for each grid point by performing local vertical difference calculation based on the actual water depth of each grid point and two adjacent control water depths close to the actual water depth; and generating a vertical grid for the target terrain region based on the series of characteristic values ​​and elevation values ​​of each grid point in the target terrain region. This method defines the control water depth based on topographic features according to topographic changes, which can represent topographic change features and ensure that the transition of the vertical grid is smooth enough when the water depth difference between adjacent nodes is large. Each grid point on the plane can have different numbers of layers in the vertical direction according to different water depths, which can adapt to complex underwater topography. At the same time, appropriate vertical grid layers are allocated in deep water and shallow water areas to meet the vertical resolution requirements and improve the model's computational efficiency and accuracy.

[0060] This disclosure provides at least one embodiment of a vertical mesh generation method suitable for complex terrain. The mesh generation method can be implemented in software, hardware, firmware, or any combination thereof, and can be loaded and executed by a processor in a device such as a mobile phone, tablet computer, laptop computer, desktop computer, or network server. This allows it to adapt to complex underwater terrain, while allocating an appropriate number of vertical mesh layers in deep and shallow water areas to meet the requirements of vertical resolution and improve the efficiency and accuracy of model calculation.

[0061] The following is for reference. Figure 1 As shown, a vertical grid generation method suitable for complex terrain provided by at least one embodiment of the present disclosure is described, the vertical grid generation method including steps S1 to S7.

[0062] S1. Obtain the terrain elevation data of the target terrain area;

[0063] S2. Generate a planar triangular or quadrilateral grid based on the terrain elevation data;

[0064] S3. Define the control water depth and grid layer number of several terrain feature control points based on the terrain elevation data.

[0065] In some embodiments, the terrain elevation data can reflect the terrain changes of the target terrain area. Specifically, the control water depth of each terrain feature control point and the corresponding number of grid layers are defined according to the terrain changes. The control water depth of the terrain feature control point can represent the terrain change characteristics, such as the location of abrupt changes such as steep slopes, deep pits, and submerged dikes. On the other hand, it can ensure that the transition of the vertical grid is smooth enough when the water depth difference between adjacent nodes is large.

[0066] In some embodiments, when generating a grid for a target terrain area, the number of terrain feature control points and the number of terrain feature control point grid layers can be set according to the calculation target of the vertical grid to meet the calculation needs under different circumstances.

[0067] S4. Determine the critical water depth at the transition point between deep water and shallow water areas in the target terrain region;

[0068] S5. Calculate the characteristic value and elevation value of each grid point based on the actual water depth and critical water depth of each grid point within the target terrain area.

[0069] In some embodiments, step S5, such as Figure 2 As shown, steps S51 to S53 are included:

[0070] S51. Compare the actual water depth at each grid point with the critical water depth.

[0071] In some embodiments, the actual water depth of each grid point within the target terrain area is determined from the elevation data of the target terrain area. After comparison with the critical water depth, it is possible to distinguish between deep water areas and shallow water areas, thereby allocating appropriate vertical grid layers to the deep water areas and shallow water areas respectively, meeting the vertical grid resolution requirements, and improving the model's computational efficiency and accuracy.

[0072] S52. When the actual water depth of a grid point is greater than or equal to the critical water depth, the characteristic value and elevation value of that grid point are calculated respectively using the preset first set of formulas.

[0073] In some embodiments, the first set of formulas includes a first formula for calculating elevation values ​​and a second formula for calculating characteristic values; the first formula is:

[0074]

[0075] In the formula, σ is the eigenvalue. h c Where η is the critical water depth, h is the actual water depth, and η is the reference water level.

[0076] The second formula is:

[0077]

[0078] In the formula, σ is the eigenvalue, and θ b and θ f Define and set coefficients to control the density of the target terrain area at the bottom and water surface in the grid;

[0079] The θ b and θ f Used to control the density of the bottom and water surface of the target terrain area in the grid;

[0080] When θ f When the coordinates approach zero, the vertical coordinate of the grid point becomes the characteristic value, and the elevation value is calculated using the following formula:

[0081]

[0082] In the formula, h c η is the critical water depth, h is the actual water depth, η is the reference water level, and σ is the characteristic value.

[0083] S53. When the actual water depth at a grid point is less than the critical water depth, the characteristic value and elevation value at that grid point are calculated using the preset second set of formulas.

[0084] In some embodiments, the second set of formulas includes a third formula for calculating elevation values ​​and a fourth formula for calculating characteristic values; the third formula is:

[0085] z k =h i σk=1,2,3......,N i ;

[0086] In the formula, h i The control water depth for the topographic features of this grid point is σ, where σ is the characteristic value.

[0087] The fourth formula is: σ=ασ 2 +(1+α)σ;

[0088] In the formula, the magnitude of σ is evenly distributed between -1 and 0, where -1 represents the bottom and 0 represents the water surface; α is a defined set coefficient.

[0089] S6. Based on the actual water depth of each grid point and the two adjacent control water depths close to the actual water depth, perform local vertical difference calculations to generate a series of characteristic values ​​and elevation values ​​for each grid point.

[0090] In some embodiments, when the thickness of the bottommost grid layer of the interpolated grid point is less than the preset minimum grid layer thickness, the grid layer is merged with the grid layer above it to preserve the actual water depth of the grid point.

[0091] S7. Generate the vertical grid of the target terrain region based on a series of feature values ​​and elevation values ​​of each grid point in the target terrain region.

[0092] In some embodiments, step S7, such as Figure 3 As shown, it includes:

[0093] S71. For grid points in the vertical grid where the actual water depth is greater than or equal to the critical water depth, when its θ f When the value is greater than or equal to 1, the vertical grid at the boundary of the grid point is refined;

[0094] S72. For grid points in the vertical grid where the actual water depth is greater than or equal to the critical water depth, when its θ b When θ approaches 0, the vertical grid at the point near the water surface is refined; when its θ... b When the value approaches 1, the vertical grid on both the surface and bottom of the grid point is refined.

[0095] In some embodiments, step S7, such as Figure 3 As shown, it also includes:

[0096] S73. For grid points in the vertical grid where the actual water depth is less than the critical water depth, when α = 0, the traditional characteristic value is obtained through the second formula group; when α < 0, the vertical grid near the water surface of the grid point is densified in layers; when α > 0, the vertical grid near the bottom of the grid point is densified in layers.

[0097] In some embodiments, step S7, such as Figure 3 As shown, it also includes:

[0098] S74. When adjacent grid points have different vertical grid layers, the vertical grid at the bottom connection is represented by trapezoidal elements.

[0099] Using trapezoidal elements for connection makes the transition between adjacent grid points smoother and more seamless, improving the accuracy of the vertical grid.

[0100] In some specific embodiments, eight terrain feature control points are set within the target terrain area, θ f =1,θ b =1, the number of terrain feature control points is 10, thus generating a vertical grid as follows: Figure 4 As shown in the figure, vertically uniform grids are generated in both shallow and deep water areas, and the grids can smoothly transition between the deep and shallow water areas.

[0101] The vertical mesh generation method for complex terrain provided by the embodiments of this disclosure can define the control water depth of terrain features according to terrain changes, represent terrain change features, and ensure that the transition of the vertical mesh is smooth enough when the water depth difference between adjacent nodes is large. Each grid point on the plane can have different numbers of layers in the vertical direction according to different water depths, which can adapt to complex underwater terrain. At the same time, appropriate vertical mesh layers are allocated at different water depths to meet the vertical resolution requirements and improve the model calculation efficiency and accuracy.

[0102] Example 2:

[0103] The difference between this embodiment and Embodiment 1 is that the number of terrain feature control points set in the target terrain area is 8, θ f =1,θ b =1, the number of terrain feature control points is 5, thus generating a vertical grid as follows: Figure 5 As shown in the figure, a vertically uniform grid is generated in both shallow and deep water areas, and the grid transitions smoothly between the deep and shallow water areas. Compared with the grid generated by using 10 layers of terrain feature control points in Example 1, the number of vertical layers in the calculation area is reduced.

[0104] Example 3:

[0105] The difference between this embodiment and Embodiment 1 is that the number of terrain feature control points set in the target terrain area is 8, θ f =1,θ b =5, the number of terrain feature control points is 10, thus generating a vertical grid as follows: Figure 6 As shown in the figure, a denser vertical mesh is generated at the surface and bottom of the water body, which is suitable for fluid calculations that require high accuracy at the surface and boundary.

[0106] Example 4:

[0107] The difference between this embodiment and Embodiment 1 is that the number of terrain feature control points set in the target terrain area is 8, θ f =0, θ b =1, the number of terrain feature control points is 10, thus generating a vertical grid as follows: Figure 7 As shown in the figure, a denser vertical grid is generated on the surface of the water body, which is suitable for fluid calculations requiring high accuracy at the surface.

[0108] Example 5:

[0109] At least some embodiments of this disclosure also provide an electronic device, such as Figure 8As shown, the electronic device includes a memory 21 and a processor 22. The memory 21 stores a computer program that, when executed by the processor, performs the steps of the vertical mesh generation method provided in any embodiment of this disclosure.

[0110] In some embodiments, processor 22 is configured to perform all or part of the steps in the vertical grid generation method based on a knowledge graph attention network, as described in any embodiment of this disclosure. Memory 21 is configured to store various types of data, which may include, for example, instructions for any application or method in an electronic device, as well as application-related data.

[0111] The processor 22 may be implemented as an application-specific integrated circuit (ASIC), a digital signal processor (DSP), a digital signal processing device (DSPD), a programmable logic device (PLD), a field-programmable gate array (FPGA), a controller, a microcontroller, a microprocessor, or other electronic components, and is used to execute the vertical mesh generation method in Embodiment 1 above.

[0112] The memory 21 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk.

[0113] Example 6:

[0114] At least some embodiments of this disclosure also provide a computer-readable storage medium, such as Figure 9As shown, the readable storage medium stores a computer program 31, which, when executed by a processor, implements the steps of the vertical mesh generation method provided in any embodiment of this disclosure.

[0115] In some embodiments, the storage medium may be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, apparatus, or device. A machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. Machine-readable media may include, but are not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0116] To provide interaction with a user, the systems and techniques described herein can be implemented on a computer having: a display device for displaying information to the user (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor); and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the computer. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0117] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as a data server), or computing systems that include middleware components (e.g., an application server), or computing systems that include frontend components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with embodiments of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., a communication network). Examples of communication networks include local area networks (LANs), wide area networks (WANs), and the Internet.

[0118] Computer systems can include clients and servers. Clients and servers are generally located far apart and typically interact via communication networks. Client-server relationships are created by computer programs running on the respective computers and having a client-server relationship with each other. Servers can be cloud servers, servers in distributed systems, or servers incorporating blockchain technology.

[0119] This invention also provides a computer program product, including a computer program / instructions that, when executed by a processor, implement the steps of the vertical mesh generation method provided in any embodiment of this disclosure.

[0120] In summary, the vertical mesh generation method, device, and storage medium provided in this application, applicable to complex terrain, can define the control water depth of terrain features according to terrain changes to represent terrain change features, and ensure that the transition of the vertical mesh is smooth enough when the water depth difference between adjacent nodes is large. Each grid point on the plane can have different numbers of layers in the vertical direction according to different water depths, which can adapt to complex underwater terrain. At the same time, it can allocate appropriate vertical mesh layers at different water depths to meet the vertical resolution requirements and improve the model's computational efficiency and accuracy.

[0121] The various embodiments in this disclosure are described in a progressive manner. The same or similar parts between the various embodiments can be referred to each other. Each embodiment focuses on describing the differences from other embodiments.

[0122] The scope of protection of this disclosure is not limited to the embodiments described above. Obviously, those skilled in the art can make various modifications and variations to this disclosure without departing from its scope and spirit. If such modifications and variations fall within the scope of the claims of this disclosure and their equivalents, then the intent of this disclosure also includes such modifications and variations.

Claims

1. A method for generating vertical meshes suitable for complex terrain, characterized in that, Includes the following steps: Obtain the topographic elevation data of the target terrain area; A planar triangular or quadrilateral grid is generated based on the terrain elevation data; Based on the terrain elevation data, define the control water depth and the number of grid layers for several terrain feature control points; Determine the critical water depth at the transition point between deep and shallow water areas in the target topographic region; Calculate the characteristic value and elevation value of each grid point based on the actual water depth and critical water depth of each grid point within the target terrain area; Local vertical interpolation is performed based on the actual water depth of each grid point and two adjacent control water depths close to the actual water depth to generate a series of characteristic values ​​and elevation values ​​for each grid point; The vertical grid of the target terrain region is generated based on a series of feature values ​​and elevation values ​​of each grid point in the target terrain region. Specifically, when adjacent grid points have different numbers of vertical grid layers, the vertical grid at the bottom connection is represented by trapezoidal units to make the transition between adjacent grid points smoother and more seamless. The step of calculating the characteristic value and elevation value of each grid point based on the actual water depth and critical water depth within the target terrain area includes: Compare the actual water depth at each grid point with the critical water depth. When the actual water depth at a grid point is greater than or equal to the critical water depth, the characteristic value and elevation value at that grid point are calculated using the first set of preset formulas. When the actual water depth at a grid point is less than the critical water depth, the characteristic value and elevation value at that grid point are calculated using the preset second set of formulas. The first set of formulas includes a first formula for calculating elevation values ​​and a second formula for calculating characteristic values; the first formula is: In the formula, σ is the eigenvalue. h c Where η is the critical water depth, h is the actual water depth, and η is the reference water level. The second formula is: In the formula, θ b and θ f Define and set coefficients to control the density of the target terrain area at the bottom and water surface in the grid; The second set of formulas includes a third formula for calculating elevation values ​​and a fourth formula for calculating characteristic values; the third formula is: ； In the formula, hi is the control water depth of the topographic feature control point of the grid point, and σ is the feature value; The fourth formula is: ; In the formula, the magnitude of σ is evenly distributed between -1 and 0, where -1 represents the bottom and 0 represents the water surface; α is a defined set coefficient.

2. The method for generating vertical meshes suitable for complex terrain according to claim 1, characterized in that, When θf approaches zero, the vertical coordinate becomes the characteristic value, and the elevation value is calculated using the following formula: ； In the formula, , h c η is the critical water depth, h is the actual water depth, η is the reference water level, and σ is the characteristic value.

3. The method for generating vertical grids suitable for complex terrain according to claim 1, characterized in that, The process of generating a vertical grid for the target terrain region based on a series of feature values ​​and elevation values ​​of each grid point includes: For grid points in the vertical grid where the actual water depth is greater than or equal to the critical water depth, when its θ f When the value is greater than or equal to 1, the vertical grid at the boundary of the grid point is refined; For grid points in the vertical grid where the actual water depth is greater than or equal to the critical water depth, when its θ b When θ approaches 0, the vertical grid at the point near the water surface is refined; when its θ... b When the value approaches 1, the vertical grid on both the surface and bottom of the grid point is refined.

4. The method for generating vertical meshes suitable for complex terrain according to claim 1, characterized in that, The process of generating a vertical grid for the target terrain region based on a series of feature values ​​and elevation values ​​of each grid point includes: For grid points in the vertical grid where the actual water depth is less than the critical water depth, when α=0, the traditional eigenvalue is obtained through the second set of formulas; when α<0, the vertical grid near the water surface of the grid point is densified in layers; when α>0, the vertical grid near the bottom of the grid point is densified in layers.

5. The method for generating vertical grids suitable for complex terrain according to claim 1, characterized in that, The process involves local vertical interpolation based on the actual water depth at each grid point and two adjacent control water depths close to that actual water depth, generating a series of characteristic values ​​and elevation values ​​for each grid point, including: When the thickness of the bottommost grid layer of the interpolated grid point is less than the preset minimum grid layer thickness, the grid layer is merged with the grid layer above it to preserve the actual water depth of the grid point.

6. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the vertical mesh generation method according to any one of claims 1 to 5.

7. A computer-readable storage medium, characterized in that, When the computer program is executed by a processor, it implements the steps of the vertical mesh generation method as described in any one of claims 1 to 5.

Images