River inundation drawing method, device and equipment based on section water level and storage medium

Abstract

The invention discloses a section water level-based riverway submerging drawing method, which comprises the following steps of: acquiring section data of each riverway section, riverbed roughness of a riverway, upper boundary hourly flow and lower boundary water level, and determining the section water level and a riverway section water level line; calculating a first submerging distance and a plurality of second submerging distances from the first section river bank point of each river section to the central section of the target river; determining the coordinate of a first side river bank point adjacent to the first section river bank point according to the third coordinate of the river channel point closest to the river channel section, the second submerging distance and the first coordinate of the first section river bank point; coordinates of other first side river bank points and coordinates of second side river bank points; and drawing a first bank according to the coordinates of the first-side river bank point, and drawing a second bank according to the coordinates of the second-side river bank point. According to the method, the submerging range can be accurately determined, the flood submerging area can be pre-judged in advance, the flood submerging loss is evaluated, the flood control plan is made, and threats to life and property safety are reduced.

Description

Technical Field

[0001] The present invention relates to the technical field of river channel flooding range mapping, and in particular to a river channel flooding mapping method, device, equipment and storage medium based on cross-section water level. Background Art

[0002] Floods, as highly destructive natural disasters, often cause enormous loss of life and property. They can destroy homes, inundate farmland, disrupt transportation and infrastructure, severely disrupt normal social order, and even threaten human survival. Efficiently simulating flood processes and quickly and accurately determining the inundation area during a flood is crucial for impact and loss assessment, as well as for decision-making on flood prevention and disaster reduction.

[0003] Currently, one-dimensional hydrodynamic models are widely used for rapid flood simulation. Based on the principles of conservation of mass and momentum, these models simplify water motion into a one-dimensional flow. By solving the Saint-Venant equations, they can rapidly simulate the temporal evolution of water levels and flow at key sections, providing fundamental data for subsequent flood analysis.

[0004] In real-world flood scenarios, simply understanding the inundation extent of a river section is far from sufficient. The impact of a flood is not limited to the specific section, but extends to the vast areas between them. This lack of effective understanding of the inundation extent between sections makes it impossible to comprehensively and accurately assess the extent of flood impacts and losses in different areas. This hinders effective protection of life and property, and prevents the rational planning of flood control facilities. Summary of the Invention

[0005] Based on this, it is necessary to propose a river channel submergence mapping method, device, equipment and storage medium based on cross-sectional water level to address the above problems.

[0006] A method for mapping river channel inundation based on cross-sectional water levels, wherein a first bank of the river channel is composed of a plurality of first-side river bank points, and a second bank of the river channel is composed of a plurality of second-side river bank points, wherein the first-side river bank points and the second-side river bank points are located on either side of a river channel centerline, the method comprising:

[0007] Loading a Tiandi Map satellite image into an ArcGIS geographic information system, determining a target river channel in the Tiandi Map satellite image, and determining a river channel centerline of the target river channel;

[0008] Defining multiple river sections of the target river, and obtaining cross-sectional data, riverbed roughness, upper boundary hourly flow, and lower boundary water level of each river section; the cross-sectional data includes: river section lines and river section curves;

[0009] Determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level;

[0010] Determining a river section water level line according to the cross-section water level;

[0011] Obtaining a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, where the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point;

[0012] The centerline of the river between two adjacent river sections is divided into multiple river center segments at equal intervals. The connecting points of adjacent river center segments are river points, and the third coordinate of each river point is determined;

[0013] Calculate the first flooding distance from the first bank point of each river section to the target river center section in sequence, where the target river center section is the line connecting the two river points closest to the river section;

[0014] Determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are the vertical distances from each first side riverbank point to the corresponding nearest river channel center section;

[0015] determining the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and repeating this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0016] The first bank is drawn according to the coordinates of each first-side river bank point, and the second bank is drawn according to the coordinates of each second-side river bank point.

[0017] In one embodiment, determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate, and the lower boundary water level includes:

[0018] The river channel cross-sectional curve and the riverbed roughness are used as inputs of a one-dimensional hydrodynamic model, the upper boundary hourly flow rate is used as the upper boundary condition of the one-dimensional hydrodynamic model, and the lower boundary water level is used as the lower boundary condition of the one-dimensional hydrodynamic model, so that the one-dimensional hydrodynamic model outputs the cross-sectional water level.

[0019] In one embodiment, determining a plurality of second submergence distances according to the two first submergence distances respectively corresponding to two adjacent river sections includes:

[0020] A plurality of second flooding distances are determined by performing difference calculation on two first flooding distances corresponding to two adjacent river sections.

[0021] In one embodiment, the first submergence distance is calculated as follows:

[0022]

[0023] Where D1 is the first flooding distance, and the coordinates of the first section bank point S1 are (x s ,y s ), the coordinates of the two river points closest to the river section are P1(x1,y1) and P2(x2,y2).

[0024] In one embodiment, the second submergence distance is calculated as follows:

[0025]

[0026] Among them, D i is the second flooding distance, D1 and D n are the two first flooding distances corresponding to two adjacent river sections respectively, i is the i-th riverbank point, and n is the n-th riverbank point.

[0027] A river channel submergence mapping device based on cross-section water level, comprising:

[0028] A loading module is used to load the Tiandi Map satellite image into the ArcGIS geographic information system, determine the target river channel in the Tiandi Map satellite image, and determine the centerline of the target river channel;

[0029] A first acquisition module is used to define multiple river sections of the target river and obtain cross-sectional data, riverbed roughness, upper boundary hourly flow rate, and lower boundary water level of each river section; the cross-sectional data includes: river section line and river section curve;

[0030] A first determination module is used to determine the cross-section water level according to the river cross-section curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level;

[0031] A second determining module is used to determine the water level line of the river section according to the water level of the section;

[0032] a second acquisition module, configured to acquire a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, wherein the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point;

[0033] A third determination module is used to split the river centerline between two adjacent river sections into multiple river center segments at equal intervals, with the connecting points of adjacent river center segments being river points, and determine the third coordinate of each river point;

[0034] a calculation module for sequentially calculating a first flooding distance from a first bank point of each river channel section to a target river channel center section, where the target river channel center section is a line connecting two river channel points closest to the river channel section;

[0035] a fourth determining module, configured to determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are vertical distances from each first side riverbank point to the corresponding nearest river channel center segment;

[0036] a fifth determining module, configured to determine the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and to repeat this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0037] The drawing module is used to draw the first bank according to the coordinates of each first-side river bank point, and to draw the second bank according to the coordinates of each second-side river bank point.

[0038] In one embodiment,

[0039] The second determination module is also used to use the river section curve and the riverbed roughness as inputs of a one-dimensional hydrodynamic model, the upper boundary hourly flow rate as the upper boundary condition of the one-dimensional hydrodynamic model, and the lower boundary water level as the lower boundary condition of the one-dimensional hydrodynamic model, so that the one-dimensional hydrodynamic model outputs the cross-sectional water level.

[0040] In one embodiment,

[0041] The fourth determination module is further configured to determine a plurality of second flooding distances by performing difference calculation on two first flooding distances corresponding to two adjacent river sections.

[0042] A computer device comprises a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor performs the following steps:

[0043] Loading a Tiandi Map satellite image into an ArcGIS geographic information system, determining a target river channel in the Tiandi Map satellite image, and determining a river channel centerline of the target river channel;

[0044] Defining multiple river sections of the target river, and obtaining cross-sectional data, riverbed roughness, upper boundary hourly flow, and lower boundary water level of each river section; the cross-sectional data includes: river section lines and river section curves;

[0045] Determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level;

[0046] Determining a river section water level line according to the cross-section water level;

[0047] Obtaining a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, where the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point;

[0048] The centerline of the river between two adjacent river sections is divided into multiple river center segments at equal intervals. The connecting points of adjacent river center segments are river points, and the third coordinate of each river point is determined;

[0049] Calculate the first flooding distance from the first bank point of each river section to the target river center section in sequence, where the target river center section is the line connecting the two river points closest to the river section;

[0050] Determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are the vertical distances from each first side riverbank point to the corresponding nearest river channel center section;

[0051] determining the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and repeating this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0052] A computer-readable storage medium stores a computer program, which, when executed by a processor, causes the processor to perform the following steps:

[0053] Loading a Tiandi Map satellite image into an ArcGIS geographic information system, determining a target river channel in the Tiandi Map satellite image, and determining a river channel centerline of the target river channel;

[0054] Defining multiple river sections of the target river, and obtaining cross-sectional data, riverbed roughness, upper boundary hourly flow, and lower boundary water level of each river section; the cross-sectional data includes: river section lines and river section curves;

[0055] Determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level;

[0056] Determining a river section water level line according to the cross-section water level;

[0057] Obtaining a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, where the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point;

[0058] The centerline of the river between two adjacent river sections is divided into multiple river center segments at equal intervals. The connecting points of adjacent river center segments are river points, and the third coordinate of each river point is determined;

[0059] Calculate the first flooding distance from the first bank point of each river section to the target river center section in sequence, where the target river center section is the line connecting the two river points closest to the river section;

[0060] Determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are the vertical distances from each first side riverbank point to the corresponding nearest river channel center section;

[0061] determining the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and repeating this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0062] The first bank is drawn according to the coordinates of each first-side river bank point, and the second bank is drawn according to the coordinates of each second-side river bank point.

[0063] The above-mentioned river flooding range drawing method provided in this application can generate a two-dimensional flooding range consistent with the river channel trend based on the cross-sectional water level calculated by the one-dimensional hydrodynamic model, which helps relevant departments to predict in advance the areas that may be flooded by floods, assess the losses caused by floods, and formulate more scientific and reasonable flood control plans, thereby reducing the threat of floods to people's lives and property safety. BRIEF DESCRIPTION OF THE DRAWINGS

[0064] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0065] in:

[0066] Figure 1 A diagram illustrating an application environment of a river channel inundation mapping method based on cross-section water levels in one embodiment;

[0067] Figure 2 A flow chart of a method for drawing river channel inundation based on cross-section water levels in one embodiment;

[0068] Figure 3 A schematic diagram of a cross-sectional water level in one embodiment;

[0069] Figure 4 A schematic diagram of river channel points and river bank points on a river channel in one embodiment;

[0070] Figure 5 This is a rendering of the river flooding range in one embodiment;

[0071] Figure 6 A structural block diagram of a river channel submergence drawing device based on cross-section water level in one embodiment;

[0072] Figure 7 FIG. 1 is a structural block diagram of a computer device in one embodiment. DETAILED DESCRIPTION

[0073] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

[0074] Figure 1 This is an application environment diagram of a river channel flooding drawing method based on cross-section water level in one embodiment. Figure 1This cross-sectional water level-based river channel inundation mapping method is applied to a cross-sectional water level-based river channel inundation mapping system. This cross-sectional water level-based river channel inundation mapping system includes a terminal 110 and a server 120. Terminal 110 and server 120 are connected via a network. Terminal 110 can be a desktop terminal or a mobile terminal. The mobile terminal can be at least one of a mobile phone, tablet computer, and laptop computer. Server 120 can be implemented as a standalone server or a server cluster consisting of multiple servers. The terminal 110 is used to: load the Tiandi Map satellite image in the ArcGIS geographic information system, determine the target river in the Tiandi Map satellite image, and determine the centerline of the target river; define multiple river sections of the target river, and obtain the section data, riverbed roughness, upper boundary hourly flow and lower boundary water level of each river section; the section data includes: a river section line and a river section curve; determine the section water level according to the river section curve, the riverbed roughness, upper boundary hourly flow and lower boundary water level; determine the river section water level line according to the section water level; the server 120 is used to obtain the first coordinate of the first section riverbank point and the second coordinate of the second section riverbank point, the two intersection points of the river section water level line and the river section curve are the first section riverbank point and the second section riverbank point of the river section, the first section riverbank point belongs to the first side riverbank point, and the second section riverbank point belongs to the second side riverbank point; two adjacent rivers are separated The center line of the river channel between the channel sections is divided into multiple river channel center sections at equal intervals, the connection points of adjacent river channel center sections are river channel points, and the third coordinate of each river channel point is determined; the first inundation distance from the first section river bank point of each river channel section to the target river channel center section is calculated in sequence, and the target river channel center section is the line connecting the two river channel points closest to the river channel section; multiple second inundation distances are determined according to the two first inundation distances corresponding to two adjacent river channel sections, and the second inundation distance is the vertical distance from each first side river bank point to the corresponding nearest river channel center section; the coordinates of the first side river bank point adjacent to the first section river bank point are determined according to the third coordinate of the river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section river bank point; and this step is repeated to calculate the coordinates of other first side river bank points and the coordinates of the second side river bank points; the first bank is drawn according to the coordinates of each first side river bank point, and the second bank is drawn according to the coordinates of each second side river bank point.

[0075] Floods, as extremely destructive natural disasters, often cause enormous loss of life and property. They can destroy houses, inundate farmland, disrupt transportation and infrastructure, severely disrupt normal social order, and even threaten human survival. Efficiently simulating flood processes and quickly and accurately determining the inundation extent during floods is crucial for flood impact and loss assessment, as well as for flood prevention and mitigation decision-making. Currently, one-dimensional hydrodynamic models are widely used for rapid flood simulation. Based on the principles of conservation of mass and momentum, one-dimensional hydrodynamic models simplify water flow into a one-dimensional flow. By solving the Saint-Venant equations, they rapidly simulate the temporal evolution of water levels and flow at key sections, providing fundamental data for subsequent flood analysis. In real-world flood scenarios, simply understanding the inundation extent of a river section is far from sufficient. The impact of a flood is not limited to the section itself, but extends to the vast areas between sections. The lack of effective understanding of the inundation extent between sections hinders comprehensive and accurate assessment of the impact of floods on different regions during flood impact and loss assessments. It is impossible to effectively protect the safety of people's lives and property and it is impossible to reasonably plan flood control facilities. In order to solve the above technical problems, this application provides a river channel flooding mapping method based on cross-sectional water level. The first bank of the river channel is composed of a number of first side river bank points, and the second bank of the river channel is composed of a number of second side river bank points. The first side river bank points and the second side river bank points are located on both sides of the center line of the river channel, such as Figure 2 As shown, the method includes:

[0076] S10: loading a Tiandi Map satellite image into an ArcGIS geographic information system, determining a target river channel in the Tiandi Map satellite image, and determining a river channel centerline of the target river channel;

[0077] S20: defining multiple river sections of the target river, and obtaining cross-sectional data of each river section, riverbed roughness, upper boundary hourly flow rate, and lower boundary water level; the cross-sectional data includes: a river section line and a river section curve;

[0078] S30: determining a cross-section water level according to the river channel cross-section curve, the riverbed roughness, the upper boundary hourly flow rate, and the lower boundary water level;

[0079] S40: determining a river section water level line according to the cross-section water level;

[0080] S50: Obtain the first coordinates of the first-section riverbank point S1 (i.e., the left-side riverbank point S1) and the second coordinates of the second-section riverbank point S2. The two intersection points of the river section water level line and the river section curve are the first-section riverbank point S1 and the second-section riverbank point S2 of the river section, respectively. The first-section riverbank point S1 belongs to the first-side riverbank point, and the second-section riverbank point S2 belongs to the second-side riverbank point.

[0081] S60: Split the river centerline between two adjacent river sections into multiple river center segments at equal intervals, with the connecting points of adjacent river center segments as river points, and determine the third coordinate of each river point P;

[0082] S70: sequentially calculating a first flooding distance D1 from a first bank point S1 of each river channel section to a target river channel center section, where the target river channel center section is a line connecting two river channel points closest to the river channel section;

[0083] S80: Determine multiple second flooding distances D according to the two first flooding distances D1 corresponding to two adjacent river sections. i , the second flooding distance D i The vertical distance from each first-side riverbank point to the corresponding nearest river channel center segment;

[0084] S90: According to the third coordinate of the river point closest to the river section, the second flooding distance D i , determining the coordinates of the first side riverbank point adjacent to the first section riverbank point S1 using the first coordinate of the first section riverbank point S1; and repeating this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0085] S100: Draw a first bank according to the coordinates of each first-side riverbank point, and draw a second bank according to the coordinates of each second-side riverbank point.

[0086] The above-mentioned river flooding range drawing method provided in this application can generate a two-dimensional flooding range consistent with the river channel trend based on the cross-sectional water level calculated by the one-dimensional hydrodynamic model, which helps relevant departments to predict in advance the areas that may be flooded by floods, assess the losses caused by floods, and formulate more scientific and reasonable flood control plans, thereby reducing the threat of floods to people's lives and property safety.

[0087] In one embodiment, determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate, and the lower boundary water level in step S30 includes:

[0088] S301: The river channel cross-sectional curve and the riverbed roughness are used as inputs of a one-dimensional hydrodynamic model, the upper boundary hourly flow rate is used as the upper boundary condition of the one-dimensional hydrodynamic model (HEC-RAS model), and the lower boundary water level is used as the lower boundary condition of the one-dimensional hydrodynamic model (HEC-RAS model), so that the one-dimensional hydrodynamic model (HEC-RAS model) outputs the cross-sectional water level.

[0089] In one embodiment, in step S80, the second flooding distances D are determined based on the two first flooding distances D1 corresponding to the two adjacent river sections. i include:

[0090] S801: Determine multiple second flooding distances D by performing difference calculation on the two first flooding distances D1 corresponding to two adjacent river sections. i .

[0091] In one embodiment, the calculation method of the first flooding distance D1 in step S70 is as follows:

[0092]

[0093] Where D1 is the first flooding distance, and the coordinates of the first section bank point S1 are (x s ,y s ), the coordinates of the two river points closest to the river section are P1(x1,y1) and P2(x2,y2);

[0094] In one embodiment, for the second flooding distance D in step S801 i The calculation method is as follows:

[0095]

[0096] Among them, D i is the second flooding distance, D1 and D n are the two first flooding distances corresponding to two adjacent river sections respectively, i is the i-th riverbank point, and n is the n-th riverbank point.

[0097] Specifically, we use actual cases to illustrate the solution of this application:

[0098] At a specific moment, in the geographic information system, the river section water level line is determined based on the cross-section water level. The two intersection points of the river section water level line and the river section curve are the first section river bank point S1 and the second section river bank point S2 of the river section. The first section river bank point S1 belongs to the first side river bank point, and the second section river bank point S2 belongs to the second side river bank point. The first coordinate of the first section river bank point S1 and the second coordinate of the second section river bank point S2 are determined. The first section river bank point S1 and the second section river bank point S2 are on both sides of the river centerline. Figure 3 shown.

[0099] The center line of the river between two adjacent river sections is divided into multiple river center segments at equal intervals. The connection points of adjacent river center segments are river points, and the coordinates of each river point are determined, such as Figure 4 As shown in the figure, the river points are represented as P1 to P7. It is ensured that the winding shape of the river can be reflected after the river points are connected in sequence.

[0100] Calculate the first flooding distance D1 from the first bank point S1 of each river channel section to the river channel center section in sequence, where the target river channel center section is the line connecting the two river channel points closest to the river channel section;

[0101] from Figure 4 It can be seen that there is a first section riverbank point S1 on one of the river sections S1S2, which is located on the left side of the river centerline. ， The distance between every two channel points on the channel centerline is equal. The vertical distance between the first section bank point S1 and the line connecting the first two channel points P1 and P2 is the first inundation distance D1. Calculate the vertical distance D7 between the line connecting the left bank point E1 and the last two channel points P6 and P7 of the other channel section E1E2. Based on the first inundation distance D1 and the vertical distance D7, interpolate the inundation distances of each channel point between the two sections to calculate the left inundation distance of each channel point.

[0102] Riverbank point S1(x s ,y s ) to the line connecting the river point P1 (x1, y1) and the river point P2 (x2, y2) as follows:

[0103]

[0104] Riverbank point E1(x e ,y e ) to the line connecting river point P6 (x3, y3) and river point P7 (x4, y4) as follows:

[0105]

[0106] Interpolate each river point between the two sections using the linear interpolation method. The left flooding distance of the i-th river point is calculated as follows:

[0107]

[0108] For each river point P between the two sections i , according to its relationship with the previous river point P i-1 Vertical direction of the connecting line, flooding distance D on the left i , and the previous riverbank point b i-1 Coordinates, determine the left bank point b corresponding to each river point i The coordinates of the river bank point b1 are the river bank point S1, and the river bank point b7 are the river bank point E1.

[0109] According to the river point P i-1 (x1,y1), river point P i (x2,y2), river bank point b i-1 (x m ,y m ), left flooding distance D i To calculate b i (x n ,y n ) is as follows:

[0110]

[0111] dot=(x1-x m )(x2-x n )+(y1-y m )(y2-y n ) (8)

[0112]

[0113] Where x' n ,y' n for b i The coordinates of the trial point, dot represents and When dot>0, it means and The angle is less than 90°, b i Trial point and b i-1 On the straight line P i-1 ipsilateral to P, b i Coordinates are taken from the trial point coordinates (x' n ,y' n ); when dot<0, it means and The angle is greater than 90°, b iwith b i-1 On the straight line P i-1 On the opposite side of P, b needs to be adjusted i Coordinate calculation formula, get b i Coordinate (x n ,y n ).

[0114] Connect the left bank point S1 of the first river section, each of the left bank points b2 to b6, and the left bank point E1 of the second river section to form the left bank line. Calculate the vertical distance between the right bank point S2 of the first river section and the line connecting the first two river points P1 and P2. Calculate the vertical distance between the right bank point E2 of the second river section and the line connecting the last two river points P6 and P7. Based on these two vertical distances, interpolate the inundation distances of each river point between the two sections to calculate the right inundation distance of each river point.

[0115] For each river point P between the two sections i , according to its relationship with the previous river point P i-1 The vertical direction of the connecting line, the right-side inundation distance, and the coordinates of the previous riverbank point are used to determine the coordinates of the right-side riverbank point corresponding to each river channel point.

[0116] The right bank point S2 of the first river section, each right bank point, and the right bank point E2 of the second river section are connected in sequence to form a right bank line.

[0117] Connect the left riverbank and the right riverbank at the beginning and end to form a closed curve, which serves as the flooding range between the two sections.

[0118] Repeat the above operation for every two adjacent sections to obtain the flooding range of the entire river channel. Figure 5 shown.

[0119] The above operation is performed at each calculation moment to obtain the change process of the flooding range of the entire river channel.

[0120] This application also provides a river channel flooding mapping device based on cross-sectional water level, such as Figure 6 Shown, including:

[0121] The loading module 10 is used to load the Tiandi Map satellite image into the ArcGIS geographic information system, determine the target river channel in the Tiandi Map satellite image, and determine the centerline of the target river channel;

[0122] The first acquisition module 20 is used to define multiple river sections of the target river and obtain cross-sectional data, riverbed roughness, upper boundary hourly flow rate and lower boundary water level of each river section; the cross-sectional data includes: river section line and river section curve;

[0123] A first determining module 30 is configured to determine a cross-sectional water level based on the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate, and the lower boundary water level;

[0124] A second determining module 40 is used to determine the water level line of the river section according to the water level of the section;

[0125] a second acquisition module 50, configured to acquire a first coordinate of the first section riverbank point and a second coordinate of the second section riverbank point, wherein the two intersection points of the river section water level line and the river section curve are the first section riverbank point and the second section riverbank point of the river section, respectively; the first section riverbank point belongs to a first side riverbank point, and the second section riverbank point belongs to a second side riverbank point;

[0126] A third determination module 60 is configured to split the channel centerline between two adjacent channel sections into a plurality of channel center segments at equal intervals, wherein the connecting points of adjacent channel center segments are channel points, and determine the third coordinate of each channel point;

[0127] A calculation module 70 is configured to sequentially calculate a first flooding distance from a first bank point of each river channel section to a target river channel center segment, where the target river channel center segment is a line connecting two river channel points closest to the river channel section;

[0128] A fourth determining module 80 is configured to determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are vertical distances from each first side riverbank point to the corresponding nearest river channel center segment;

[0129] a fifth determining module 90 configured to determine the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of the river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and to repeat this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0130] The drawing module 100 is used to draw the first bank according to the coordinates of each first-side riverbank point, and to draw the second bank according to the coordinates of each second-side riverbank point.

[0131] In one embodiment, the second determination module is further used to use the river section curve and the riverbed roughness as inputs of a one-dimensional hydrodynamic model, use the upper boundary hourly flow rate as the upper boundary condition of the one-dimensional hydrodynamic model (HEC-RAS model), and use the lower boundary water level as the lower boundary condition of the one-dimensional hydrodynamic model (HEC-RAS model), so that the one-dimensional hydrodynamic model (HEC-RAS model) outputs the cross-sectional water level.

[0132] In one embodiment, the fourth determining module is further configured to determine a plurality of second flooding distances by performing difference calculation on two first flooding distances corresponding to two adjacent river sections.

[0133] In one embodiment, a computer device is provided, comprising a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor performs the following steps:

[0134] Loading a Tiandi Map satellite image into an ArcGIS geographic information system, determining a target river channel in the Tiandi Map satellite image, and determining a river channel centerline of the target river channel;

[0135] Defining multiple river sections of the target river, and obtaining cross-sectional data, riverbed roughness, upper boundary hourly flow, and lower boundary water level of each river section; the cross-sectional data includes: river section lines and river section curves;

[0136] Determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level;

[0137] Determining a river section water level line according to the cross-section water level;

[0138] Obtaining a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, where the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point;

[0139] The centerline of the river between two adjacent river sections is divided into multiple river center segments at equal intervals. The connecting points of adjacent river center segments are river points, and the third coordinate of each river point is determined;

[0140] Calculate the first flooding distance from the first bank point of each river section to the target river center section in sequence, where the target river center section is the line connecting the two river points closest to the river section;

[0141] Determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are the vertical distances from each first side riverbank point to the corresponding nearest river channel center section;

[0142] determining the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and repeating this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0143] The first bank is drawn according to the coordinates of each first-side river bank point, and the second bank is drawn according to the coordinates of each second-side river bank point.

[0144] In one embodiment, a computer-readable storage medium is provided, storing a computer program. When the computer program is executed by a processor, the processor performs the following steps:

[0145] Loading a Tiandi Map satellite image into an ArcGIS geographic information system, determining a target river channel in the Tiandi Map satellite image, and determining a river channel centerline of the target river channel;

[0146] Defining multiple river sections of the target river, and obtaining cross-sectional data, riverbed roughness, upper boundary hourly flow, and lower boundary water level of each river section; the cross-sectional data includes: river section lines and river section curves;

[0147] Determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level;

[0148] Determining a river section water level line according to the cross-section water level;

[0149] Obtaining a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, where the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point;

[0150] The centerline of the river between two adjacent river sections is divided into multiple river center segments at equal intervals. The connecting points of adjacent river center segments are river points, and the third coordinate of each river point is determined;

[0151] Calculate the first flooding distance from the first bank point of each river section to the target river center section in sequence, where the target river center section is the line connecting the two river points closest to the river section;

[0152] Determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are the vertical distances from each first side riverbank point to the corresponding nearest river channel center section;

[0153] determining the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and repeating this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points;

[0154] The first bank is drawn according to the coordinates of each first-side river bank point, and the second bank is drawn according to the coordinates of each second-side river bank point.

[0155] Figure 7 FIG1 shows an internal structure diagram of a computer device in an embodiment. The computer device can be a terminal or a server. Figure 7 As shown, the computer device includes a processor, a memory, and a network interface connected via a system bus. The memory includes a non-volatile storage medium and an internal memory. The non-volatile storage medium of the computer device stores an operating system and may also store a computer program. When the computer program is executed by the processor, the processor can implement a river channel submergence mapping method based on a cross-sectional water level. The internal memory may also store a computer program. When the computer program is executed by the processor, the processor can implement a river channel submergence mapping method based on a cross-sectional water level. Those skilled in the art will understand that Figure 7 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.

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

[0157] The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0158] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the present invention shall be determined by the appended claims.

Claims

1. A method for mapping river channel inundation based on cross-sectional water levels, wherein a first bank of a river channel is composed of a plurality of first bank points, and a second bank of a river channel is composed of a plurality of second bank points, wherein the first bank points and the second bank points are located on either side of a river channel centerline, characterized in that: The method comprises: Loading a Tiandi Map satellite image into an ArcGIS geographic information system, determining a target river channel in the Tiandi Map satellite image, and determining a river channel centerline of the target river channel; Defining multiple river sections of the target river, and obtaining cross-sectional data, riverbed roughness, upper boundary hourly flow, and lower boundary water level of each river section; the cross-sectional data includes: river section lines and river section curves; Determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level; Determining a river section water level line according to the cross-section water level; Obtaining a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, where the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point; The centerline of the river between two adjacent river sections is divided into multiple river center segments at equal intervals. The connecting points of adjacent river center segments are river points, and the third coordinate of each river point is determined; Calculate the first flooding distance from the first bank point of each river section to the target river center section in sequence, where the target river center section is the line connecting the two river points closest to the river section; Determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are the vertical distances from each first side riverbank point to the corresponding nearest river channel center section; determining the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and repeating this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points; The first bank is drawn according to the coordinates of each first-side river bank point, and the second bank is drawn according to the coordinates of each second-side river bank point.

2. The method for mapping river channel submergence based on cross-sectional water level according to claim 1, characterized in that: Determining the cross-sectional water level according to the river channel cross-sectional curve, the riverbed roughness, the upper boundary hourly flow rate, and the lower boundary water level includes: The river channel cross-sectional curve and the riverbed roughness are used as inputs of a one-dimensional hydrodynamic model, the upper boundary hourly flow rate is used as the upper boundary condition of the one-dimensional hydrodynamic model, and the lower boundary water level is used as the lower boundary condition of the one-dimensional hydrodynamic model, so that the one-dimensional hydrodynamic model outputs the cross-sectional water level.

3. The method for mapping river channel submergence based on cross-sectional water level according to claim 1, characterized in that: The determining of a plurality of second submergence distances according to the two first submergence distances respectively corresponding to two adjacent river sections comprises: A plurality of second flooding distances are determined by performing difference calculation on two first flooding distances respectively corresponding to two adjacent river sections.

4. The method for mapping river channel submergence based on cross-sectional water level according to claim 1, characterized in that: The calculation method of the first flooding distance is as follows: Where D1 is the first flooding distance, and the coordinates of the first section bank point S1 are (x s ,y s ), the coordinates of the two river points closest to the river section are P1(x1,y1) and P2(x2,y2).

5. The method for mapping river channel submergence based on cross-sectional water level according to claim 4, characterized in that: The calculation method of the second flooding distance is as follows: Among them, D i is the second flooding distance, D1 and D n are the two first flooding distances corresponding to two adjacent river sections respectively, i is the i-th riverbank point, and n is the n-th riverbank point.

6. A river channel flooding mapping device based on cross-sectional water level, characterized in that: include: A loading module is used to load the Tiandi Map satellite image into the ArcGIS geographic information system, determine the target river channel in the Tiandi Map satellite image, and determine the river channel centerline of the target river channel; A first acquisition module is used to define multiple river sections of the target river and obtain cross-sectional data, riverbed roughness, upper boundary hourly flow rate, and lower boundary water level of each river section; the cross-sectional data includes: river section line and river section curve; A first determination module is used to determine the cross-section water level according to the river cross-section curve, the riverbed roughness, the upper boundary hourly flow rate and the lower boundary water level; A second determining module is used to determine the water level line of the river section according to the water level of the section; a second acquisition module, configured to acquire a first coordinate of the first-section riverbank point and a second coordinate of the second-section riverbank point, wherein the two intersection points of the river section water level line and the river section curve are the first-section riverbank point and the second-section riverbank point of the river section, respectively; the first-section riverbank point belongs to a first-side riverbank point, and the second-section riverbank point belongs to a second-side riverbank point; A third determination module is used to split the river centerline between two adjacent river sections into multiple river center segments at equal intervals, with the connecting points of adjacent river center segments being river points, and determine the third coordinate of each river point; a calculation module for sequentially calculating a first flooding distance from a first bank point of each river channel section to a target river channel center section, where the target river channel center section is a line connecting two river channel points closest to the river channel section; a fourth determining module, configured to determine a plurality of second inundation distances based on the two first inundation distances corresponding to two adjacent river sections, wherein the second inundation distances are vertical distances from each first side riverbank point to the corresponding nearest river channel center segment; a fifth determining module, configured to determine the coordinates of a first side riverbank point adjacent to the first section riverbank point based on the third coordinate of a river channel point closest to the river channel section, the second inundation distance, and the first coordinate of the first section riverbank point; and to repeat this step to calculate the coordinates of other first side riverbank points and the coordinates of the second side riverbank points; The drawing module is used to draw the first bank according to the coordinates of each first-side river bank point, and to draw the second bank according to the coordinates of each second-side river bank point.

7. The river channel submergence mapping device based on cross-sectional water level according to claim 1, characterized in that: The second determination module is also used to use the river section curve and the riverbed roughness as inputs of a one-dimensional hydrodynamic model, the upper boundary hourly flow rate as the upper boundary condition of the one-dimensional hydrodynamic model, and the lower boundary water level as the lower boundary condition of the one-dimensional hydrodynamic model, so that the one-dimensional hydrodynamic model outputs the cross-sectional water level.

8. The river channel submergence mapping device based on cross-sectional water level according to claim 1, characterized in that: The fourth determination module is further configured to determine a plurality of second flooding distances by performing difference calculation on two first flooding distances corresponding to two adjacent river sections.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor is caused to perform the steps of the method according to any one of claims 1 to 5.

10. A computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor is caused to perform the steps of the method according to any one of claims 1 to 5.

Images