CN111444561A - River channel terrain model establishing method and device - Google Patents
River channel terrain model establishing method and device Download PDFInfo
- Publication number
- CN111444561A CN111444561A CN202010141812.3A CN202010141812A CN111444561A CN 111444561 A CN111444561 A CN 111444561A CN 202010141812 A CN202010141812 A CN 202010141812A CN 111444561 A CN111444561 A CN 111444561A
- Authority
- CN
- China
- Prior art keywords
- target
- grid
- data
- boundary
- river channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000004590 computer program Methods 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 1
- 238000012876 topography Methods 0.000 description 12
- 238000004364 calculation method Methods 0.000 description 7
- 230000001154 acute effect Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
The application discloses a river terrain model establishing method and device, and relates to the technical field of river terrain. The method for establishing the river terrain model comprises the following steps: acquiring boundary data and key section data of a target river channel, and obtaining a target data group through interpolation based on the boundary data and the key section data; generating a two-dimensional orthogonal curve grid based on the target data set, and interpolating to obtain a Z coordinate of each target grid node based on the X coordinate and the Y coordinate of each target node in the two-dimensional orthogonal curve grid; interpolating to obtain Z coordinates of all other grid nodes in the two-dimensional orthogonal curve grid based on the X coordinate, the Y coordinate and the Z coordinate of each target grid node and the boundary data of the target river channel; therefore, a terrain three-dimensional model of the target river channel is established based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid. The method and the device for establishing the river channel terrain model can reduce the economic cost required by establishing the river channel terrain model and improve the efficiency of establishing the river channel terrain model.
Description
Technical Field
The application relates to the technical field of river terrain, in particular to a river terrain model building method and device.
Background
The river terrain is an important characteristic of a water system watershed and is also a basis for carrying out water flow data simulation calculation and risk prevention and control effect simulation on rivers.
In the prior art, the river terrain is obtained through manual field surveying, and a large amount of manpower, material resources, financial resources and time are consumed in the actual measurement process. The prior art has the defects of high economic cost and low acquisition efficiency in acquiring the river terrain through manual field surveying.
Disclosure of Invention
The application provides a river channel terrain model building method and device, which can reduce the economic cost for building a river channel terrain model and improve the river channel terrain model building efficiency.
In order to achieve the above technical effect, a first aspect of the present application provides a method for establishing a river terrain model, including:
acquiring boundary data and key section data of a target river channel, wherein the boundary data is a three-dimensional coordinate of a boundary discrete point of the target river channel, and the key section data is a three-dimensional coordinate of a key section discrete point of the target river channel;
obtaining a target data group through interpolation based on the boundary data and the key section data of the target river channel, wherein the target data group comprises a target boundary data group and a target section data group;
generating a two-dimensional orthogonal curve grid based on the target data set;
obtaining a Z coordinate of each target grid node through interpolation based on the X coordinate and the Y coordinate of each target grid node in the two-dimensional orthogonal curve grid, wherein the target grid node is a grid node corresponding to the target section data set;
obtaining Z coordinates of all grid nodes except the target grid node in the two-dimensional orthogonal curve grid through interpolation based on the X coordinate, the Y coordinate and the Z coordinate of each target grid node and the boundary data of the target river channel;
and establishing a terrain three-dimensional model of the target river channel based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid.
Optionally, the acquiring boundary data and key section data of the target river specifically includes: and acquiring boundary data and key section data of the target river channel from the existing standard terrain file.
Optionally, the acquiring boundary data and key section data of the target river specifically includes:
calling the existing electronic map software to respectively acquire longitude and latitude coordinates and elevation information of the boundary discrete point and the key section discrete point of the target river channel;
and acquiring boundary data and key section data of the target river channel based on longitude and latitude coordinates and elevation information of the boundary discrete points and the key section discrete points.
Optionally, after obtaining the target data group through interpolation, the method further includes: drawing a boundary line of the target river channel based on the target data group;
the generating of the two-dimensional orthogonal curve grid based on the target data set specifically includes: and generating a two-dimensional orthogonal curve grid based on the target data group and the boundary line.
The second aspect of the present application provides a river terrain model building apparatus, including:
the data acquisition module is used for acquiring boundary data and key section data of a target river channel, wherein the boundary data is a three-dimensional coordinate of a boundary discrete point of the target river channel, and the key section data is a three-dimensional coordinate of a key section discrete point of the target river channel;
the first interpolation module is used for obtaining a target data set through interpolation based on the boundary data and the key section data of the target river channel, wherein the target data set comprises a target boundary data set and a target section data set;
a grid generating module for generating a two-dimensional orthogonal curve grid based on the target data set;
a second interpolation module, configured to obtain a Z coordinate of each target grid node through interpolation based on an X coordinate and a Y coordinate of each target grid node in the two-dimensional orthogonal curve grid, where the target grid node is a grid node corresponding to the target section data set;
a third interpolation module, configured to obtain, based on the X coordinate, the Y coordinate, and the Z coordinate of each target grid node and the boundary data of the target river, Z coordinates of all grid nodes except the target grid node in the two-dimensional orthogonal curve grid through interpolation;
and the model establishing module is used for establishing a terrain three-dimensional model of the target river channel based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid.
Optionally, the data obtaining module is specifically configured to:
and acquiring boundary data and key section data of the target river channel from the existing standard terrain file.
Optionally, the data obtaining module includes a map calling sub-module and a coordinate conversion sub-module:
the map calling sub-module is used for calling the existing electronic map software to respectively acquire longitude and latitude coordinates and elevation information of the boundary discrete point and the key section discrete point of the target river channel;
and the coordinate conversion sub-module is used for acquiring boundary data and key section data of the target river channel based on longitude and latitude coordinates and elevation information of the boundary discrete points and key section discrete points.
Optionally, the river terrain model building device further includes a boundary drawing module;
the boundary drawing module is used for drawing a boundary line of the target river channel based on the target data group;
the mesh generation module is configured to generate a two-dimensional orthogonal curve mesh based on the target data set and the boundary line.
A third aspect of the present application provides a computer device, which includes a memory and a processor, where the memory stores a computer program, and the processor implements the steps of the method for establishing a river terrain model when executing the computer program.
A fourth aspect of the present application provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the steps of the method for creating a river terrain model.
According to the technical scheme, the target data group is obtained through interpolation by acquiring boundary data and key section data of the target river channel and based on the boundary data and the key section data; generating a two-dimensional orthogonal curve grid based on the target data set, and obtaining a Z coordinate of each target grid node through interpolation based on the X coordinate and the Y coordinate of each target node in the two-dimensional orthogonal curve grid; interpolating to obtain Z coordinates of all other grid nodes in the two-dimensional orthogonal curve grid based on the X coordinate, the Y coordinate and the Z coordinate of each target grid node and the boundary data of the target river channel; therefore, a terrain three-dimensional model of the target river channel is established based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid. Because this scheme obtains the required data point of establishing river course topography model through the interpolation after obtaining the boundary data and the key section data of target river course, need not carry out a large amount of manual work and surveys on the spot, consequently, for the scheme that obtains the river course topography through artificial wetland survey among the prior art, the required economic cost of establishing river course topography model has greatly been reduced and the efficiency of establishing river course topography model has been improved to this application.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the embodiments or the prior art descriptions will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive exercise.
Fig. 1 is a flowchart of a method for establishing a river terrain model according to an embodiment of the present disclosure;
fig. 2 is a structural block diagram of a river terrain model building apparatus according to an embodiment of the present application.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system structures, techniques, etc. in order to provide a thorough understanding of the embodiments of the invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.
It will be understood that the terms "comprises" and/or "comprising," when used in this specification and the appended claims, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
It is also to be understood that the terminology used in the description of the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used in the specification of the present application and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
It should be further understood that the term "and/or" as used in this specification and the appended claims refers to and includes any and all possible combinations of one or more of the associated listed items.
As used in this specification and the appended claims, the term "if" may be interpreted contextually as "when …" or "upon" or "in response to a determination" or "in response to a detection". Similarly, the phrase "if it is determined" or "if a [ described condition or event ] is detected" may be interpreted depending on the context to mean "upon determining" or "in response to determining" or "upon detecting [ described condition or event ]" or "in response to detecting [ described condition or event ]".
The technical solutions in the embodiments of the present application are clearly and completely described below with reference to the drawings of the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present application, but the present application may be practiced in other ways than those described herein, and it will be apparent to those of ordinary skill in the art that the present application is not limited by the specific embodiments disclosed below.
Fig. 1 shows a river terrain model building method provided in an embodiment of the present application, which is detailed as follows:
The boundary data is three-dimensional coordinates of boundary discrete points of the target river channel, and the key section data is three-dimensional coordinates of key section discrete points of the target river channel.
In an application scenario, the acquiring boundary data and key section data of the target river specifically includes: and acquiring boundary data and key section data of the target river channel from the existing standard terrain file. Specifically, step 101 includes: and acquiring a standard terrain file containing data of a target river channel, and acquiring longitude and latitude coordinates and elevation information of boundary discrete points of the target river channel and longitude and latitude coordinates and elevation information of key section discrete points of the target river channel from the standard terrain file. Respectively converting longitude and latitude coordinates and elevation information of the boundary discrete point of the target river channel and the key section discrete point of the target river channel into three-dimensional coordinates required by establishment of a river channel terrain model through a formula (1):
the method comprises the following steps of establishing a river channel terrain model, establishing a discrete point on the basis of a discrete point coordinate system, establishing a river channel terrain model, establishing a discrete point coordinate system on the basis of the discrete point coordinate system, and solving the river channel terrain model, wherein X is an abscissa of the discrete point when the river channel terrain model is established, Y is a ordinate of the discrete point when the river channel terrain model is established, Z is a vertical coordinate of the discrete point when the river channel terrain model is established, B is a latitude of the discrete point, L is a longitude of the discrete point, H is a geodetic height of the discrete point:
wherein, a is the length of the earth's major semi-axis, and b is the length of the earth's minor semi-axis.
Optionally, the key section is a national assessment section in the target river channel region.
In another application scenario, the acquiring boundary data and key section data of the target river includes: calling the existing electronic map software to respectively acquire longitude and latitude coordinates and elevation information of the boundary discrete point and the key section discrete point of the target river channel; and acquiring boundary data and key section data of the target river channel based on longitude and latitude coordinates and elevation information of the boundary discrete points and the key section discrete points. Optionally, after the existing electronic map software is called, a target river channel is found and a river channel boundary and a key section of the target river channel are defined, so that boundary discrete points and key section discrete points of the target river channel, and longitude and latitude coordinates and elevation information corresponding to each discrete point are obtained. And respectively converting longitude and latitude coordinates and elevation information of the boundary discrete point of the target river channel and the key section discrete point of the target river channel into three-dimensional coordinates required by establishment of a river channel terrain model through the formula (1), the formula (2) and the formula (3).
And 102, obtaining a target data group through interpolation based on the boundary data and the key section data of the target river channel.
Wherein the target data set comprises a target boundary data set and a target section data set.
Optionally, based on the boundary data and the key section data of the target river channel, a target data group required by the establishment of the river channel terrain model is obtained through interpolation by an inverse distance weighted method. In the interpolation process, the topographic fluctuation change among the discrete points is small, so that the reverse distance weight interpolation method which is relatively simple and convenient to calculate is adopted, and the calculation efficiency is improved.
And 103, generating a two-dimensional orthogonal curve grid based on the target data set.
Optionally, after the target data group is obtained through interpolation, the method further includes: drawing a boundary line of the target river channel based on the target data group; the generating of the two-dimensional orthogonal curve grid based on the target data set specifically includes: and generating a two-dimensional orthogonal curve grid based on the target data group and the boundary line. Specifically, the generated two-dimensional orthogonal curve grid is close to the boundary line of the target river channel as much as possible.
In an application scenario, after the generating the two-dimensional orthogonal curve grid, the method further includes detecting whether a first target sub-grid exists in the two-dimensional orthogonal curve grid, if so, adjusting the two-dimensional orthogonal curve grid to improve orthogonality of the first target sub-grid, and returning to re-detect whether the first target sub-grid still exists in the orthogonal curve grid until the first target sub-grid does not exist in the two-dimensional orthogonal curve grid. The first target sub-grid is a grid in which the cosine value of the minimum acute angle with the grid node as the vertex in the grid node of the two-dimensional orthogonal curve grid is greater than a preset first cosine threshold; the improving the orthogonality of the first target sub-mesh is specifically achieved by topology optimization, deleting a long and narrow mesh, or adjusting the positions of mesh vertices.
In another application scenario, after the generating the two-dimensional orthogonal curve grid, the method further includes detecting whether a second target sub-grid or a third target sub-grid exists in the two-dimensional orthogonal curve grid, if so, adjusting the two-dimensional orthogonal curve grid to improve orthogonality of the second target sub-grid and the third target sub-grid, and returning to re-detect whether the second target sub-grid or the third target sub-grid still exists in the orthogonal curve grid until the second target sub-grid and the third target sub-grid do not exist in the two-dimensional orthogonal curve grid. The second target sub-grid is a grid in which the cosine value of the minimum acute angle with the grid node as the vertex at the grid node in the two-dimensional orthogonal curve grid is greater than a preset first cosine threshold value and the terrain does not belong to a preset first terrain range; the second target sub-grid is a grid in which the cosine value of the minimum acute angle with the grid node as the vertex in the grid node of the two-dimensional orthogonal curve grid is greater than a preset second cosine threshold value, and the terrain in which the grid node is located belongs to a preset first terrain range. In this embodiment, the first cosine threshold is 0.02, and the second cosine threshold is 0.05. Optionally, the preset first terrain range includes one or more of a shoal terrain, a shoal terrain and a river boundary terrain in water, and is not specifically limited herein.
Optionally, after the two-dimensional orthogonal grid is generated, some river channel branches and the river center continent terrain in the target river channel are processed in a manner of adding or deleting part of the grid, so that the calculation amount for establishing the river channel terrain model is reduced, and the establishment efficiency of the river channel terrain model is improved.
And 104, obtaining the Z coordinate of each target grid node through interpolation based on the X coordinate and the Y coordinate of each target grid node in the two-dimensional orthogonal curve grid.
And the target grid node is a grid node corresponding to the target section data set.
And 105, obtaining the Z coordinates of all grid nodes except the target grid node in the two-dimensional orthogonal curve grid through interpolation based on the X coordinate, the Y coordinate and the Z coordinate of each target grid node and the boundary data of the target river channel.
Optionally, in the step 104 and the step 105, a common kriging interpolation method is used for interpolation, and Z coordinates of each target node and all grid nodes except the target grid node in the two-dimensional orthogonal curve grid are respectively obtained. Because the topographic relief change between each scattered point of the key section is large, the common kriging interpolation method is adopted, the calculation precision is improved, and the accuracy of the river channel topographic model is improved.
And 106, establishing a terrain three-dimensional model of the target river channel based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid.
Optionally, after the establishing of the three-dimensional terrain model of the target river channel, the method further includes outputting and displaying the three-dimensional terrain model of the target river channel.
Optionally, after the establishing of the three-dimensional terrain model of the target river channel, obtaining terrain parameters of the target river channel, and marking and displaying the terrain parameters on the three-dimensional terrain model of the target river channel. The topographic parameters include one or more of river length, river width, roughness, river course shoreline coordinates and emphasis section elevation distribution, and may further include other topographic parameters, which are not specifically limited herein.
As can be seen from the above, the method for establishing a river terrain model provided in the embodiment of the present application obtains a target data group by obtaining boundary data and key section data of a target river and interpolating based on the boundary data and the key section data; generating a two-dimensional orthogonal curve grid based on the target data set, and obtaining a Z coordinate of each target grid node through interpolation based on the X coordinate and the Y coordinate of each target node in the two-dimensional orthogonal curve grid; interpolating to obtain Z coordinates of all other grid nodes in the two-dimensional orthogonal curve grid based on the X coordinate, the Y coordinate and the Z coordinate of each target grid node and the boundary data of the target river channel; therefore, a terrain three-dimensional model of the target river channel is established based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid. Because this scheme obtains the required data point of establishing river course topography model through the interpolation after obtaining the boundary data and the key section data of target river course, need not carry out a large amount of manual work and surveys on the spot, consequently, for the scheme that obtains the river course topography through artificial wetland survey among the prior art, the required economic cost of establishing river course topography model has greatly been reduced and the efficiency of establishing river course topography model has been improved to this application.
Corresponding to the method for establishing a river terrain model provided in the foregoing embodiment, fig. 2 illustrates a river terrain model establishing apparatus provided in the second aspect of the embodiment of the present application. For convenience of explanation, only the portions related to the present embodiment are shown. Unless the present embodiment clearly indicates otherwise, what is not explicitly stated in the present embodiment corresponds to the above-mentioned method for establishing a river terrain model.
Referring to fig. 2, the river terrain model building apparatus includes:
a data obtaining module 201, configured to obtain boundary data and key section data of a target river, where the boundary data is a three-dimensional coordinate of a boundary discrete point of the target river, and the key section data is a three-dimensional coordinate of a key section discrete point of the target river;
a first interpolation module 202, configured to obtain a target data set through interpolation based on the boundary data and the key section data of the target river, where the target data set includes a target boundary data set and a target section data set;
a grid generating module 203, configured to generate a two-dimensional orthogonal curve grid based on the target data set;
a second interpolation module 204, configured to obtain a Z coordinate of each target grid node through interpolation based on an X coordinate and a Y coordinate of each target grid node in the two-dimensional orthogonal curve grid, where the target grid node is a grid node corresponding to the target section data set;
a third interpolation module 205, configured to obtain, through interpolation, Z coordinates of all grid nodes except the target grid node in the two-dimensional orthogonal curve grid based on the X coordinate, the Y coordinate, and the Z coordinate of each target grid node and the boundary data of the target river;
and a model establishing module 206, configured to establish a three-dimensional model of the terrain of the target river based on the two-dimensional orthogonal curve grid and the Z coordinate of each grid node in the two-dimensional orthogonal curve grid.
In an application scenario, the data obtaining module 201 is specifically configured to: and acquiring boundary data and key section data of the target river channel from the existing standard terrain file. Specifically, a standard terrain file containing data of a target river channel is obtained, and longitude and latitude coordinates and elevation information of boundary discrete points of the target river channel and longitude and latitude coordinates and elevation information of key section discrete points of the target river channel are obtained from the standard terrain file. And (3) respectively converting longitude and latitude coordinates and elevation information of the boundary discrete points of the target river channel and the key section discrete points of the target river channel into three-dimensional coordinates required by establishment of a river channel terrain model through a formula (4):
the method comprises the following steps of establishing a river channel terrain model, establishing a discrete point on the basis of a discrete point coordinate system, establishing a river channel terrain model, establishing a discrete point coordinate system on the basis of the discrete point coordinate system, and solving the river channel terrain model, wherein X is an abscissa of the discrete point when the river channel terrain model is established, Y is a ordinate of the discrete point when the river channel terrain model is established, Z is a vertical coordinate of the discrete point when the river channel terrain model is established, B is a latitude of the discrete point, L is a longitude of the discrete point, H is a geodetic height of the discrete point:
wherein, a is the length of the earth's major semi-axis, and b is the length of the earth's minor semi-axis.
Optionally, the key section is a national assessment section in the target river channel region.
In another application scenario, the data obtaining module 201 is configured to: calling the existing electronic map software to respectively acquire longitude and latitude coordinates and elevation information of the boundary discrete point and the key section discrete point of the target river channel; and acquiring boundary data and key section data of the target river channel based on longitude and latitude coordinates and elevation information of the boundary discrete points and the key section discrete points. Optionally, after the existing electronic map software is called, a target river channel is found and a river channel boundary and a key section of the target river channel are defined, so that boundary discrete points and key section discrete points of the target river channel, and longitude and latitude coordinates and elevation information corresponding to each discrete point are obtained. And respectively converting longitude and latitude coordinates and elevation information of the boundary discrete point of the target river channel and the key section discrete point of the target river channel into three-dimensional coordinates required by establishment of a river channel terrain model through the formula (4), the formula (5) and the formula (6).
Optionally, the first interpolation module is specifically configured to: and based on the boundary data and the key section data of the target river channel, obtaining a target data group required by establishing a river channel terrain model through interpolation by an inverse distance weight method. In the interpolation process, the topographic fluctuation change among the discrete points is small, so that the reverse distance weight interpolation method which is relatively simple and convenient to calculate is adopted, and the calculation efficiency is improved.
Optionally, the river terrain model building apparatus further includes a boundary drawing module (not shown in the figure), and the boundary drawing module is configured to draw a boundary line of the target river based on the target data set. The grid generating module is specifically configured to generate a two-dimensional orthogonal curve grid based on the target data group and the boundary line, so that the generated two-dimensional orthogonal curve grid is as close as possible to the boundary line of the target river.
In an application scenario, the above river terrain model building apparatus further includes a grid verification module (not shown in the figure) configured to: and detecting whether a first target sub-grid exists in the two-dimensional orthogonal curve grid, if so, adjusting the two-dimensional orthogonal curve grid to improve the orthogonality of the first target sub-grid, and returning to re-detect whether the first target sub-grid still exists in the orthogonal curve grid until the first target sub-grid does not exist in the two-dimensional orthogonal curve grid. The first target sub-grid is a grid in which the cosine value of the minimum acute angle with the grid node as the vertex in the grid node of the two-dimensional orthogonal curve grid is greater than a preset first cosine threshold; the improving the orthogonality of the first target sub-mesh is specifically achieved by topology optimization, deleting a long and narrow mesh, or adjusting the positions of mesh vertices.
In another application scenario, the river terrain model building apparatus further includes a grid verification module (not shown in the figure) configured to: and detecting whether a second target sub-grid or a third target sub-grid exists in the two-dimensional orthogonal curve grid, if so, adjusting the two-dimensional orthogonal curve grid to improve the orthogonality of the second target sub-grid and the third target sub-grid, and returning to detect whether the second target sub-grid or the third target sub-grid still exists in the orthogonal curve grid again until the second target sub-grid and the third target sub-grid do not exist in the two-dimensional orthogonal curve grid. The second target sub-grid is a grid in which the cosine value of the minimum acute angle with the grid node as the vertex at the grid node in the two-dimensional orthogonal curve grid is greater than a preset first cosine threshold value and the terrain does not belong to a preset first terrain range; the second target sub-grid is a grid in which the cosine value of the minimum acute angle with the grid node as the vertex in the grid node of the two-dimensional orthogonal curve grid is greater than a preset second cosine threshold value, and the terrain in which the grid node is located belongs to a preset first terrain range. In this embodiment, the first cosine threshold is 0.02, and the second cosine threshold is 0.05. Optionally, the preset first terrain range includes one or more of a shoal terrain, a shoal terrain and a river boundary terrain in water, and is not specifically limited herein.
Optionally, the river terrain model building apparatus further includes a grid correction module (not shown in the figure), configured to process some river branches and river center continent terrains in the target river in a manner of adding or deleting part of grids after the grid generation module 203 generates the two-dimensional orthogonal grid, so as to reduce the amount of calculation for building the river terrain model and improve the river terrain model building efficiency.
Optionally, the second interpolation module 204 and the third interpolation module 205 perform interpolation specifically by using a common kriging interpolation method, so as to improve the calculation precision, thereby improving the accuracy of the river terrain model.
Optionally, the river terrain model building apparatus further includes a display module (not shown in the figure), configured to obtain the terrain three-dimensional model of the target river, which is built by the model generating module 206, and output and display the terrain three-dimensional model of the target river.
Optionally, the display module is further configured to obtain a terrain parameter of the target river channel, and mark and display the terrain parameter on the terrain three-dimensional model of the target river channel. The topographic parameters include one or more of river length, river width, roughness, river course shoreline coordinates and emphasis section elevation distribution, and may further include other topographic parameters, which are not specifically limited herein.
As can be seen from the above, the river terrain model establishing apparatus provided in the embodiment of the present application obtains boundary data and key section data of a target river through the data obtaining module 201, and obtains a target data group through interpolation based on the boundary data and the key section data through the first interpolation module 202; generating a two-dimensional orthogonal curve grid based on the target data set through a grid generating module 203, and interpolating through a second interpolation module 204 to obtain a Z coordinate of each target grid node based on an X coordinate and a Y coordinate of each target node in the two-dimensional orthogonal curve grid; then, the third interpolation module 205 interpolates to obtain Z coordinates of all other grid nodes in the two-dimensional orthogonal curve grid based on the X coordinate, the Y coordinate, and the Z coordinate of each target grid node and the boundary data of the target river; therefore, a terrain three-dimensional model of the target river channel is established through the model generation module 206 based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid. Because this scheme obtains the required data point of establishing river course topography model through the interpolation after obtaining the boundary data and the key section data of target river course, need not carry out a large amount of manual work and surveys on the spot, consequently, for the scheme that obtains the river course topography through artificial wetland survey among the prior art, the required economic cost of establishing river course topography model has greatly been reduced and the efficiency of establishing river course topography model has been improved to this application.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-mentioned division of the functional units and modules is illustrated, and in practical applications, the above-mentioned functions may be distributed as different functional units and modules according to needs, that is, the internal structure of the apparatus may be divided into different functional units or modules to implement all or part of the above-mentioned functions. Each functional unit and module in the embodiments may be integrated in one processing unit, or each unit may exist alone physically, or two or more units are integrated in one unit, and the integrated unit may be implemented in a form of hardware, or in a form of software functional unit. In addition, specific names of the functional units and modules are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working processes of the units and modules in the system may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the above embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to the related descriptions of other embodiments for parts that are not described or illustrated in a certain embodiment.
Corresponding to the above embodiments, the present application also provides a computer device and a computer-readable storage medium. The computer device includes a memory and a processor, where the memory stores a computer program, and the processor implements the steps of the method for establishing a river terrain model according to the embodiment when executing the computer program. The computer readable storage medium stores a computer program, and the computer program, when executed by a processor, implements the steps of the method for establishing a river terrain model provided by the above embodiments.
Those of ordinary skill in the art would appreciate that the elements and algorithm steps of the examples described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus/terminal device and method may be implemented in other ways. For example, the above-described embodiments of the apparatus/terminal device are merely illustrative, and for example, the division of the above modules or units is only one logical division, and the actual implementation may be implemented by another division, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed.
The integrated modules/units described above, if implemented in the form of software functional units and sold or used as separate products, may be stored in a computer readable storage medium. Based on such understanding, all or part of the flow in the method of the embodiments described above may be implemented by a computer program, which may be stored in a computer-readable storage medium and can implement the steps of the embodiments of the methods described above when the computer program is executed by a processor. The computer program includes computer program code, and the computer program code may be in a source code form, an object code form, an executable file or some intermediate form. The computer readable medium may include: any entity or device capable of carrying the above-mentioned computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier signal, telecommunication signal, software distribution medium, etc. It should be noted that the contents contained in the computer-readable storage medium can be increased or decreased as required by legislation and patent practice in the jurisdiction.
The above-mentioned embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those skilled in the art; the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the embodiments of the present application, and they should be construed as being included therein.
Claims (10)
1. A river terrain model building method is characterized by comprising the following steps:
acquiring boundary data and key section data of a target river channel, wherein the boundary data are three-dimensional coordinates of boundary discrete points of the target river channel, and the key section data are three-dimensional coordinates of key section discrete points of the target river channel;
obtaining a target data group through interpolation based on the boundary data and the key section data of the target river channel, wherein the target data group comprises a target boundary data group and a target section data group;
generating a two-dimensional orthogonal curve grid based on the target data set;
obtaining a Z coordinate of each target grid node through interpolation based on the X coordinate and the Y coordinate of each target grid node in the two-dimensional orthogonal curve grid, wherein the target grid node is a grid node corresponding to the target section data set;
obtaining Z coordinates of all grid nodes except the target grid node in the two-dimensional orthogonal curve grid through interpolation based on the X coordinate, the Y coordinate and the Z coordinate of each target grid node and the boundary data of the target river channel;
and establishing a terrain three-dimensional model of the target river channel based on the two-dimensional orthogonal curve grid and the Z coordinates of grid nodes in the two-dimensional orthogonal curve grid.
2. The method for establishing a river terrain model according to claim 1, wherein the acquiring of the boundary data and the key section data of the target river specifically comprises: and acquiring boundary data and key section data of the target river channel from the existing standard terrain file.
3. The method for establishing a river terrain model according to claim 1, wherein the acquiring of the boundary data and the key section data of the target river specifically comprises:
calling the existing electronic map software to respectively acquire longitude and latitude coordinates and elevation information of boundary discrete points and key section discrete points of the target river channel;
and acquiring boundary data and key section data of the target river channel based on longitude and latitude coordinates and elevation information of the boundary discrete points and the key section discrete points.
4. A method as claimed in any one of claims 1 to 3, wherein the obtaining of the target data set by interpolation further comprises: drawing a boundary line of the target river channel based on the target data group;
the generating of the two-dimensional orthogonal curve grid based on the target data group is specifically: a two-dimensional orthogonal curve grid is generated based on the target data set and the boundary line.
5. A riverway terrain model building device is characterized by comprising:
the data acquisition module is used for acquiring boundary data and key section data of a target river channel, wherein the boundary data is a three-dimensional coordinate of a boundary discrete point of the target river channel, and the key section data is a three-dimensional coordinate of a key section discrete point of the target river channel;
the first interpolation module is used for obtaining a target data group through interpolation based on the boundary data and the key section data of the target river channel, wherein the target data group comprises a target boundary data group and a target section data group;
a grid generation module for generating a two-dimensional orthogonal curve grid based on the target data set;
a second interpolation module, configured to obtain a Z coordinate of each target grid node through interpolation based on an X coordinate and a Y coordinate of each target grid node in the two-dimensional orthogonal curve grid, where the target grid node is a grid node corresponding to the target section data set;
a third interpolation module, configured to obtain, based on the X coordinate, the Y coordinate, and the Z coordinate of each target grid node and the boundary data of the target river, Z coordinates of all grid nodes except the target grid node in the two-dimensional orthogonal curve grid through interpolation;
and the model establishing module is used for establishing a terrain three-dimensional model of the target river channel based on the two-dimensional orthogonal curve grid and the Z coordinates of each grid node in the two-dimensional orthogonal curve grid.
6. The riverway terrain model building device according to claim 5, wherein the data acquisition module is specifically configured to:
and acquiring boundary data and key section data of the target river channel from the existing standard terrain file.
7. The riverway terrain model building apparatus of claim 5, wherein the data acquisition module comprises a map calling sub-module and a coordinate transformation sub-module:
the map calling sub-module is used for calling the existing electronic map software to respectively acquire longitude and latitude coordinates and elevation information of boundary discrete points and key section discrete points of the target river channel;
and the coordinate conversion sub-module is used for acquiring boundary data and key section data of the target river channel based on longitude and latitude coordinates and elevation information of the boundary discrete points and key section discrete points.
8. The river terrain model building apparatus of any one of claims 5 to 7, further comprising a boundary drawing module;
the boundary drawing module is used for drawing a boundary line of the target river channel based on the target data group;
the grid generating module is used for generating a two-dimensional orthogonal curve grid based on the target data group and the boundary line.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 6 when executing the computer program.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010141812.3A CN111444561A (en) | 2020-03-04 | 2020-03-04 | River channel terrain model establishing method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010141812.3A CN111444561A (en) | 2020-03-04 | 2020-03-04 | River channel terrain model establishing method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111444561A true CN111444561A (en) | 2020-07-24 |
Family
ID=71652851
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010141812.3A Pending CN111444561A (en) | 2020-03-04 | 2020-03-04 | River channel terrain model establishing method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111444561A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112182814A (en) * | 2020-09-11 | 2021-01-05 | 河海大学 | River course underwater terrain modeling method based on sparse cross section point data |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030023412A1 (en) * | 2001-02-14 | 2003-01-30 | Rappaport Theodore S. | Method and system for modeling and managing terrain, buildings, and infrastructure |
| CN102496168A (en) * | 2011-11-22 | 2012-06-13 | 南京大学 | Complicated river channel gridding method used for simulation of river channel hydrological numerical value |
| CN103425856A (en) * | 2012-05-24 | 2013-12-04 | 泰瑞数创科技(北京)有限公司 | Method of automatic generation of river terrain according to river cross section data |
| CN104949628A (en) * | 2015-04-14 | 2015-09-30 | 上海大学 | Method for reconstructing complex morphology of flexible platy structure based on two-dimensional orthogonal curvature |
| CN105844709A (en) * | 2016-03-25 | 2016-08-10 | 中国水利水电科学研究院 | Submerged line tracking method for flood evolution simulation in complicated river channel landform area |
| US20180010936A1 (en) * | 2016-07-07 | 2018-01-11 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | River discharge and depth estimation |
| CN108010103A (en) * | 2017-11-24 | 2018-05-08 | 武汉大学 | The quick fine generation method of river with complicated landform |
| CN108563837A (en) * | 2018-03-21 | 2018-09-21 | 中国水利水电科学研究院 | A kind of the model parameter real-time correction method and system of alluvial stream Water-sand model |
| CN108986222A (en) * | 2018-08-15 | 2018-12-11 | 四川大学 | Without branch of a river river digital terrain generation method |
-
2020
- 2020-03-04 CN CN202010141812.3A patent/CN111444561A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030023412A1 (en) * | 2001-02-14 | 2003-01-30 | Rappaport Theodore S. | Method and system for modeling and managing terrain, buildings, and infrastructure |
| CN102496168A (en) * | 2011-11-22 | 2012-06-13 | 南京大学 | Complicated river channel gridding method used for simulation of river channel hydrological numerical value |
| CN103425856A (en) * | 2012-05-24 | 2013-12-04 | 泰瑞数创科技(北京)有限公司 | Method of automatic generation of river terrain according to river cross section data |
| CN104949628A (en) * | 2015-04-14 | 2015-09-30 | 上海大学 | Method for reconstructing complex morphology of flexible platy structure based on two-dimensional orthogonal curvature |
| CN105844709A (en) * | 2016-03-25 | 2016-08-10 | 中国水利水电科学研究院 | Submerged line tracking method for flood evolution simulation in complicated river channel landform area |
| US20180010936A1 (en) * | 2016-07-07 | 2018-01-11 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | River discharge and depth estimation |
| CN108010103A (en) * | 2017-11-24 | 2018-05-08 | 武汉大学 | The quick fine generation method of river with complicated landform |
| CN108563837A (en) * | 2018-03-21 | 2018-09-21 | 中国水利水电科学研究院 | A kind of the model parameter real-time correction method and system of alluvial stream Water-sand model |
| CN108986222A (en) * | 2018-08-15 | 2018-12-11 | 四川大学 | Without branch of a river river digital terrain generation method |
Non-Patent Citations (2)
| Title |
|---|
| 耿庆斋等: "基于MO的二维正交曲线网格生成技术及其应用", 《水电能源科学》 * |
| 耿庆斋等: "基于MO的二维正交曲线网格生成技术及其应用", 《水电能源科学》, vol. 24, no. 03, 30 June 2006 (2006-06-30), pages 18 - 21 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112182814A (en) * | 2020-09-11 | 2021-01-05 | 河海大学 | River course underwater terrain modeling method based on sparse cross section point data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Klinkenberg et al. | The fractal properties of topography: a comparison of methods | |
| Erdogan | A comparision of interpolation methods for producing digital elevation models at the field scale | |
| Hutchinson et al. | ANUSPLIN version 4.4 user guide | |
| CN108763825B (en) | Numerical simulation method for simulating wind field of complex terrain | |
| CN108648271A (en) | A kind of interpolation method generating complicated landform grid model based on GIS data | |
| CN115861527A (en) | Method and device for constructing live-action three-dimensional model, electronic equipment and storage medium | |
| Ibrahim et al. | Comparison between inverse distance weighted (IDW) and Kriging | |
| CN114219917B (en) | Vulnerability restoration method for laser radar point cloud data | |
| CN111444561A (en) | River channel terrain model establishing method and device | |
| Xie et al. | Using contour lines to generate digital elevation models for steep slope areas: a case study of the Loess Plateau in North China | |
| Hou et al. | Automatic Extraction of Drainage Networks from DEMs Base on Heuristic Search. | |
| Hapep et al. | Comparison of Different DEM Generation Methods based on Open Source Datasets. | |
| CN109752723B (en) | Underwater topographic map drawing method and device | |
| CN111583406A (en) | Pole tower foot base point coordinate calculation method and device and terminal equipment | |
| Liu et al. | Evaluation of the performance of DEM interpolation algorithms for LiDAR data | |
| Xiao-Ping et al. | An algorithm for generation of DEMs from contour lines considering geomorphic features | |
| Badora et al. | Effect of DTM resolution on the determination of slope values in an upland catchment using different computational algorithms | |
| CN110853489B (en) | Water system basin basic data acquisition method and device | |
| Saksena | Investigating the role of DEM resolution and accuracy on flood inundation mapping | |
| EP1953730B1 (en) | Determining elevation values in a geocoding system | |
| Bakuła | Reduction of DTM obtained from LiDAR data for flood modeling | |
| Wang et al. | A new method of generating grid DEM from contour lines | |
| Al-Taee et al. | Comparison of two deterministic interpolation methods for predicting ground water level in Baghdad | |
| Shang et al. | Study on the improved IDW in the maanshan reach of the yangtze river modeling | |
| CN118736152A (en) | A digital riverbed reconstruction method, device and medium based on DEM |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200724 |
|
| RJ01 | Rejection of invention patent application after publication |