CN110390687B - A method for measuring erosion and sedimentation of river channels based on 3D laser scanning - Google Patents

A method for measuring erosion and sedimentation of river channels based on 3D laser scanning Download PDF

Info

Publication number
CN110390687B
CN110390687B CN201910686756.9A CN201910686756A CN110390687B CN 110390687 B CN110390687 B CN 110390687B CN 201910686756 A CN201910686756 A CN 201910686756A CN 110390687 B CN110390687 B CN 110390687B
Authority
CN
China
Prior art keywords
point cloud
registration
point
river
phase
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.)
Active
Application number
CN201910686756.9A
Other languages
Chinese (zh)
Other versions
CN110390687A (en
Inventor
陈骎
周家文
李海波
杨兴国
戚顺超
范刚
鲁功达
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan University
Original Assignee
Sichuan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan University filed Critical Sichuan University
Priority to CN201910686756.9A priority Critical patent/CN110390687B/en
Publication of CN110390687A publication Critical patent/CN110390687A/en
Application granted granted Critical
Publication of CN110390687B publication Critical patent/CN110390687B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B5/00Measuring arrangements characterised by the use of mechanical techniques
    • G01B5/18Measuring arrangements characterised by the use of mechanical techniques for measuring depth
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M10/00Hydrodynamic testing; Arrangements in or on ship-testing tanks or water tunnels
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4023Scaling of whole images or parts thereof, e.g. expanding or contracting based on decimating pixels or lines of pixels; based on inserting pixels or lines of pixels
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/64Analysis of geometric attributes of convexity or concavity

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Fluid Mechanics (AREA)
  • Geometry (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于三维激光扫描的河道冲淤测量方法,包括以下步骤:S1.河道地形测量:根据测量区域的几何形状,确定扫描仪架设站点和靶标的平面布置,获取能够满足冲淤变形测量要求和完整覆盖扫描区域的点云数据;S2.点云数据配准:进行第一期的多站点点云配准,通过靶标进行初步配准,再通过最近点迭代算法减小配准误差;接着,进行多时期的点云配准,以第一期点云坐标系为工程坐标系,通过最近点迭代算法对第二期点云进行配准,导出第二期点云;S3.河道冲淤变形计算及可视化:通过对比不同时期平面坐标相同点的高度变化值,得到河道的冲淤变化量。本发明结合三维激光扫描技术实现河道冲淤测量,具有效率高、工作量小、精度高、低风险等优势。

Figure 201910686756

The invention discloses a method for measuring scouring and silting of a river channel based on three-dimensional laser scanning. Deformation measurement requirements and point cloud data that completely covers the scanning area; S2. Point cloud data registration: carry out the first phase of multi-site point cloud registration, conduct preliminary registration through the target, and then reduce the registration through the closest point iteration algorithm Then, perform multi-period point cloud registration, take the first phase point cloud coordinate system as the engineering coordinate system, register the second phase point cloud through the closest point iteration algorithm, and derive the second phase point cloud; S3. Calculation and visualization of river erosion and deposition deformation: By comparing the height change values of the same point on the plane coordinates in different periods, the erosion and deposition change of the river channel can be obtained. The invention combines the three-dimensional laser scanning technology to realize the measurement of river erosion and sedimentation, and has the advantages of high efficiency, small workload, high precision, low risk and the like.

Figure 201910686756

Description

Three-dimensional laser scanning-based river channel erosion and deposition measuring method
Technical Field
The invention relates to river channel erosion and deposition measurement, in particular to a river channel erosion and deposition measurement method based on three-dimensional laser scanning.
Background
In the process of river development and utilization by human beings, a large amount of hydraulic engineering is built in a river channel, so that the water and sand conditions of the river are changed, and a series of river bed evolution at the downstream is triggered. The erosion and deposition evolution of the river channel can have great influence on the safety of hydraulic buildings, the operation mode of a reservoir, shipping engineering and the like, so that the research on the erosion and deposition rule of the river channel is of great significance, wherein the direct measurement of the deformation of the river bed is an effective means for revealing the evolution rule of the river channel. The traditional river channel measurement is carried out single-point measurement through a total station instrument, a sounding rod, a sounding hammer and the like, then the integral terrain is obtained through interpolation, and the deformation of a riverbed is calculated by comparing data in different periods. The patent with the application number of CN201610354918.5 discloses a river scouring depth detection device and a detection method, and has the advantages of low cost, simple method operation, capability of accurately measuring and calculating scouring deformation before and after flood and the like.
However, the device has limited application range, the outer sleeve needs to be driven into the bedrock, the device is not suitable for the situation of deeper covering layer, only can measure the erosion and deposition deformation of the sandy riverbed, and is not suitable for the pebble riverbed; the stability of the device is not high, only the influence of the dynamic water pressure is considered, and the device is possibly damaged by floaters or pushed stones in the water; the observation range is limited, and only the erosion deformation at the bed sheet point of the river can be obtained; the measuring accuracy is not high, and the measuring rod is contacted with fine sand in the pipe during measurement, so that the volume of the fine sand is compressed, and further errors are generated.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, provides a river channel erosion and deposition measuring method based on three-dimensional laser scanning, realizes river channel erosion and deposition measurement by combining a three-dimensional laser scanning technology, and has the advantages of high efficiency, small workload, high precision, low risk and the like.
The purpose of the invention is realized by the following technical scheme: a river channel erosion and deposition measuring method based on three-dimensional laser scanning comprises the following steps:
s1, river terrain measurement: determining the plane arrangement of a scanner erection station and a target according to the geometric shape of a measurement area, and acquiring point cloud data which can meet the erosion-deposition deformation measurement requirement and completely cover the scanning area;
s2, point cloud data registration: firstly, multi-site point cloud registration of a first period is carried out, preliminary registration is carried out through a target, registration errors are further reduced through a closest point iterative algorithm, and after completion, point cloud is exported and named by date; then, performing multi-period point cloud registration, using a first-period point cloud coordinate system as an engineering coordinate system, registering a second-period point cloud through a closest point iterative algorithm, and finishing deriving the second-period point cloud, wherein the second-period point cloud is named by date;
s3, calculating and visualizing erosion-deposition deformation of the river channel: through the altitude variation value of the same point of the plane coordinate in different periods of contrast, the erosion and deposition variable quantity of the river channel is obtained: because plane coordinates of point clouds in different periods are not coincident, difference calculation cannot be directly carried out, firstly, deleting the point clouds in regions outside an original point cloud river channel, and extracting a river channel point cloud boundary; then, generating a plane coordinate matrix with equal spacing, and deleting points outside the boundary to obtain a new matrix; then, obtaining a height coordinate corresponding to each point through interpolation, and storing the height coordinates in a matrix form; and finally, calculating the height difference of the corresponding points, and displaying the junction in the form of a deformed cloud picture or a contour map.
Wherein the step S1 includes the following substeps:
s101, determining an effective measurement distance rho of a scanner according to the following formula:
N=1/(ρ×h×((tan(θ+dθ)-tanθ)×dα);
wherein θ is arccos (h/ρ), N is the minimum point cloud density, h is the height from the center of the instrument to the riverbed, ρ is the effective measurement distance of the instrument, and da and do are the horizontal resolution and the vertical resolution of the scanner;
s102, acquiring the geometric shape of the measured river reach by using the satellite image, measuring the data of the width and the length of the river reach, and calculating the maximum distance L between scanner standsmax
Figure BDA0002146551660000021
Wherein d is 1/2 of the river reach width;
arranging the erection stations on the image according to the maximum distance between the erection stations, and alternately arranging the stations on two banks to avoid instrument scanning blind areas;
s103, erecting an instrument at each preset station position, and arranging at least 3 targets between two different testing stations; and then scanning at different stations until the measurement of the whole area is completed, and acquiring point cloud data which meets the erosion and deposition deformation measurement requirement and completely covers the scanning area.
Wherein the step S2 includes the following substeps:
s201, importing first-stage point cloud data, and carrying out point cloud denoising and vegetation removal; taking the scanner coordinate system of the first station as an engineering coordinate system, taking the target coordinate in the first station as an engineering coordinate, and firstly, calculating a coordinate conversion parameter by solving the following formula:
A=RA'+ΔT
in the formula:
Figure BDA0002146551660000022
r comprises three independent parameters alpha, beta,
Figure BDA0002146551660000023
Rotation angles around x, y, z axes, respectively;
Figure BDA0002146551660000031
Δ X, Δ Y, Δ Z are three translation parameters;
Figure BDA0002146551660000032
engineering coordinates for the target;
Figure BDA0002146551660000033
scanner coordinates for the target;
the simultaneous solution is carried out to obtain six parameters of alpha, beta,
Figure BDA0002146551660000034
Delta X, delta Y and delta Z, completing the initial registration of the second station to the engineering coordinate system, and then reducing the registration error by adopting a closest point iterative algorithm to realize the final registration; and repeating the registration of the 3 rd, 4 th,.., n stations;
s202, importing second-stage point clouds, denoising and removing vegetation, importing the first-stage registered point clouds in the same project, deleting the point clouds in the river channel, and reserving points of which the external topography of the river channel is not changed; secondly, registering the clouds of the first station in the second period to the clouds of the first station by using the external area of the riverway overlapped by the point clouds in the two periods by adopting a closest point iterative algorithm; finishing the registration of the remaining point clouds in the second period according to the same method;
s203, deleting the point clouds in the region outside the river channel from the two-stage point clouds after registration, reserving the rest point clouds according to the convex polygon region, respectively exporting point cloud coordinate information in a txt format, and naming by using dates.
The invention has the beneficial effects that: the river channel erosion and deposition measurement method is combined with a three-dimensional laser scanning technology to achieve river channel erosion and deposition measurement, and has the advantages of being fast in data acquisition, high in measurement accuracy, simple and convenient to operate and small in workload.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a deformed cloud diagram of erosion and deposition of a river channel downstream of a water gate in the embodiment;
FIG. 3 is a contour map of erosion deformation of a river channel downstream of a water gate in an embodiment;
Detailed Description
The technical solutions of the present invention are further described in detail below with reference to the accompanying drawings, but the scope of the present invention is not limited to the following.
The invention adopts a three-dimensional laser scanner to obtain the river channel point cloud information. The measurement objects are exposed riverbed and bank slope. Because the emitted laser is refracted on the water surface, the intensity of the return signal received by the instrument is very weak, and the terrain below the water surface is difficult to measure, so the measurement is recommended to be carried out in a dry season or when the exposed area of a riverbed is large in the specific implementation process.
As shown in fig. 1, a river channel erosion and deposition measuring method based on three-dimensional laser scanning includes the following steps:
s1, river terrain measurement: according to the geometric shape of the measurement area, determining the plane arrangement of the scanner erection station and the target, and acquiring point cloud data which can meet the erosion-deposition deformation measurement requirement and completely cover the scanning area: in the embodiment of the application, the measurement area is selected as a downstream river channel of a certain sluice engineering, and the downstream river channel is straight. The measuring instrument adopts a Ruigger VZ-2000i laser scanner, the distance measurement reaches 2000m furthest, the speed measurement reaches 3000 ten thousand points/minute, the maximum horizontal angle resolution is 0.0007 degrees, and the maximum vertical angle resolution is 0.0015 degrees;
the step S1 includes the following sub-steps:
s101, considering that the farther the relative distance between the scanner and the scanner is, the smaller the point cloud density is, the effective measurement distance of the scanner is determined according to the minimum point cloud density required by deformation observation before measurement. For a river channel with large surface relief, a larger point cloud density can be selected, for a river channel with a flat surface, a smaller point cloud density can be selected, and the effective measurement distance rho of the scanner is determined according to the following formula:
N=1/(ρ×h×((tan(θ+dθ)-tanθ)×dα);
wherein θ is arccos (h/ρ), N is the minimum point cloud density, h is the height from the center of the instrument to the riverbed, ρ is the effective measurement distance of the instrument, and da and do are the horizontal resolution and the vertical resolution of the scanner;
in the above embodiment, the minimum point cloud density required for observation is 1000 points/m2. The height h from the center of the instrument to the river bed surface is the height of the embankment plus the height of a tripod, and is about 20 m. The horizontal resolution of the scanner was set to 0.004 °, and the vertical resolution was set to 0.002 °. Converting the angle into radian when calculating the effective measuring distance; for the convenience of solution, firstly, a rho value is assumed to obtain a point cloud density N, and if the point cloud density is close to the set minimum point cloud density, the assumed rho value is used as an approximate value of the effective observation distance. Assuming that ρ is 200m, N is calculated to be about 1100 points/m2. Therefore, the effective measuring distance is 200 m;
s102, acquiring the measured river by using satellite imagesMeasuring the geometric shape of the section, measuring the width and length data of the river section, and calculating the maximum spacing L of the scanner erecting stationsmax
Figure BDA0002146551660000041
Wherein d is 1/2 of the river reach width;
arranging the erection stations on the image according to the maximum distance between the erection stations, and alternately arranging the stations on two banks to avoid instrument scanning blind areas;
in the above embodiment, the river reach width is about 300m, the length is about 650m, etc. measured from the satellite image, Lmax is calculated to be 132m, 12 stations are arranged at equal intervals of 130m on both banks, and 6 stations are respectively arranged on the right bank and the left bank;
s103, erecting an instrument at each preset station position, and arranging at least 3 targets between two different testing stations; and then scanning at different stations until the measurement of the whole area is completed, and acquiring point cloud data which meets the erosion and deposition deformation measurement requirement and completely covers the scanning area.
In the above embodiment, 3 targets were placed between two different stations. For example, when scanning is performed at the second station, 3 targets are arranged between the first and second stations and between the second and third stations. And after the scanning is finished, if the missing area is found, immediately performing the additional measurement. Likewise, the second stage scanning is done at the same site location.
S2, point cloud data registration: firstly, multi-site point cloud registration of a first period is carried out, preliminary registration is carried out through a target, registration errors are further reduced through a closest point Iteration (ICP) algorithm, and point cloud is exported and named by date after completion; then, multi-period point cloud registration is carried out, a first-period point cloud coordinate system is used as an engineering coordinate system, a second-period point cloud is registered through a closest point iterative algorithm, the second-period point cloud is derived, and the method is named by date:
specifically, in the embodiment of the present application, the step S2 includes the following sub-steps:
s201, a project is newly built in point cloud processing software, point cloud data in the same period are imported, point cloud denoising is carried out, and vegetation is removed. Taking the scanner coordinate system of the first station as an engineering coordinate system, taking the target coordinate in the first station as an engineering coordinate, and firstly, calculating a coordinate conversion parameter by solving the following formula:
A=RA'+ΔT
in the formula:
Figure BDA0002146551660000051
r comprises three independent parameters alpha, beta,
Figure BDA0002146551660000052
Rotation angles around x, y, z axes, respectively;
Figure BDA0002146551660000053
Δ X, Δ Y, Δ Z are three translation parameters;
Figure BDA0002146551660000054
engineering coordinates for the target;
Figure BDA0002146551660000055
scanner coordinates for the target;
the simultaneous solution is carried out to obtain six parameters of alpha, beta,
Figure BDA0002146551660000056
Delta X, delta Y and delta Z, completing the initial registration of the second station to the engineering coordinate system, and then reducing the registration error by adopting a closest point iterative algorithm to realize the final registration; and repeating the registration of the 3 rd, 4 th,.., n stations; in the above embodiment, n is 12.
S202, building a new project, importing second-stage point clouds, denoising and removing vegetation, importing the first-stage registered point clouds in the same project, deleting the point clouds in the river channel, and reserving points of which the external topography of the river channel is not changed; then, registering the clouds of the first station in the second stage to the clouds of the first stage by using an external region of the riverway overlapped by the clouds of the two stages by adopting a closest point Iteration (ICP) algorithm; finishing the registration of the remaining point clouds in the second period according to the same method;
s203, deleting the point clouds in the region outside the river channel from the two-stage point clouds after registration, reserving the rest point clouds according to the convex polygon region, respectively exporting point cloud coordinate information in a txt format, and naming by using dates.
S3, calculating and visualizing erosion-deposition deformation of the river channel: through the altitude variation value of the same point of the plane coordinate in different periods of contrast, the erosion and deposition variable quantity of the river channel is obtained: because plane coordinates of point clouds in different periods are not coincident, difference calculation cannot be directly carried out, firstly, deleting the point clouds in regions outside an original point cloud river channel, and extracting a river channel point cloud boundary; then, generating a plane coordinate matrix with equal spacing, and deleting points outside the boundary to obtain a new matrix; then, obtaining a height coordinate corresponding to each point through interpolation, and storing the height coordinates in a matrix form; and finally, calculating the height difference of the corresponding points, and displaying the junction in the form of a deformed cloud picture or a contour map.
Specifically, in the embodiment of the present application, the step S3 includes:
s301: and reading txt format point cloud in matlab software, storing the point cloud in a 4-row matrix, and naming the point cloud as pointclosed _1 and pointclosed _2, wherein the first row is the serial number of the point, and the second, third and fourth rows are the x, y and z coordinates of the point. And then, a delaunay algorithm is used for generating a triangular net on the basis of the first-stage point cloud, and the side and the triangular surface are numbered respectively. All the edges are stored in a 3-column matrix line _1, the first column is the serial number of the edge, and the second and third columns are the serial numbers of two points included in the edge. All triangle faces are stored in a 4-column matrix tri _1, the first column is the number of the triangle face, and the second, third and fourth are the serial numbers of three sides of each triangle face. And counting the serial numbers of the edges in the matrix for storing the triangular surface information, wherein the serial number is the inner edge of the triangular network when the serial number appears for 2 times, and the serial number is the edge on the boundary of the triangular network when the serial number appears for 1 time. And extracting point cloud coordinates of points forming the outer boundary edge, storing the point cloud coordinates as a matrix bound _1, and reordering the point sequence in the matrix counterclockwise.
S302: firstly, a meshgrid function is utilized to generate grid coordinate matrixes x _ coord and y _ coord, x-direction coordinates and y-direction coordinates of grid points are stored respectively, and coordinate increment is 0.03. And then, interpolating the vertical coordinates of the grid points on the basis of the point cloud pointclose _1 and pointclose _2 by using a gritdata function, wherein the interpolation method can be a bilinear interpolation method, a cubic spline interpolation method, a bicubic interpolation method and the like. And obtaining interpolated z-coordinate matrixes z1_ coord and z2_ coord, and subtracting the former from the latter to obtain a z-coordinate change value matrix delz _ coord.
S303: the reshape function is then applied to straighten the x _ coord, y _ coord, delz _ coord into column vectors, which are combined and stored as a matrix gridcoord. And then, judging the positions of the points in the gridcoord and the boundary bound _1 by using an inpolygon function, and storing the points inside the boundary in a matrix in _ gridcoord.
S304: and exporting the coordinates in the in _ gridcoord as a txt file, and drawing a erosion deformation cloud picture or a contour map through importing drawing software, as shown in fig. 2 and 3.
The sorting method for S032 comprises the following steps: firstly, calculating barycentric coordinates according to plane coordinates of boundary points, then generating vectors from the barycenter to all the boundary points, calculating included angles between all the vectors and the x direction, and finally sequencing the boundary points according to the sequence of the angles from small to large.
The invention is based on three-dimensional laser scanning, thus improving the data acquisition speed; the measurement precision is high, the error generated by contact measurement is reduced, and the precision of laser scanning can reach mm level; compared with the traditional single-point measurement, the deformation value of any point of the river channel plane can be obtained; the method is simple and convenient to operate, has small workload, and can be realized by using a point cloud registration function of software and a simple matlab statement; the visualization of results can be realized through a software drawing function, and the distribution condition and the deformation amount of the erosion and deposition deformation area of the river channel are vividly displayed.
The foregoing is a preferred embodiment of the present invention, it is to be understood that the invention is not limited to the form disclosed herein, but is not to be construed as excluding other embodiments, and is capable of other combinations, modifications, and environments and is capable of changes within the scope of the inventive concept as expressed herein, commensurate with the above teachings, or the skill or knowledge of the relevant art. And that modifications and variations may be effected by those skilled in the art without departing from the spirit and scope of the invention as defined by the appended claims.

Claims (2)

1.一种基于三维激光扫描的河道冲淤测量方法,其特征在于:包括以下步骤:1. a river scouring and silting measurement method based on three-dimensional laser scanning, is characterized in that: comprise the following steps: S1.河道地形测量:根据测量区域的几何形状,确定扫描仪架设站点和靶标的平面布置,获取能够满足冲淤变形测量要求和完整覆盖扫描区域的点云数据;S1. River topography measurement: According to the geometric shape of the measurement area, determine the layout of the scanner erection site and the target, and obtain point cloud data that can meet the measurement requirements of erosion and deposition deformation and completely cover the scanning area; S2.点云数据配准:首先进行第一期的多站点点云配准,通过靶标进行初步配准,再通过最近点迭代算法进一步减小配准误差,完成后将点云导出,以日期命名;接着,进行多时期的点云配准,以第一期点云坐标系为工程坐标系,通过最近点迭代算法对第二期点云进行配准,完成导出第二期点云,以日期命名;S2. Point cloud data registration: First, perform the first phase of multi-site point cloud registration, perform preliminary registration through the target, and then further reduce the registration error through the closest point iteration algorithm. After completion, export the point cloud and use the date Name; then, perform multi-period point cloud registration, take the first phase point cloud coordinate system as the engineering coordinate system, register the second phase point cloud through the closest point iterative algorithm, and complete the export of the second phase point cloud. date naming; S3.河道冲淤变形计算及可视化:通过对比不同时期平面坐标相同点的高度变化值,得到河道的冲淤变化量:由于不同期点云平面坐标不重合,无法直接进行差值计算,首先,对原始点云河道以外的区域点云进行删除,对河道点云边界进行提取;接着,生成等间距的平面坐标矩阵,对边界外的点进行删除,得到一个新的矩阵;接着,通过插值得到每一点对应的高度坐标,以矩阵的形式存放;最后,计算对应点的高度差,以变形云图或等高线图的形式展现结;S3. Calculation and visualization of the erosion and deposition deformation of the river channel: By comparing the height changes of the same point on the plane coordinates in different periods, the amount of erosion and deposition changes in the river channel can be obtained: Since the plane coordinates of the point clouds in different periods do not overlap, the difference calculation cannot be performed directly. First, Delete the point cloud of the original point cloud outside the river channel, and extract the boundary of the river point cloud; then, generate an equal-spaced plane coordinate matrix, delete the points outside the boundary, and obtain a new matrix; then, obtain a new matrix by interpolation The height coordinates corresponding to each point are stored in the form of a matrix; finally, the height difference of the corresponding points is calculated, and the knot is displayed in the form of a deformed cloud map or contour map; 所述步骤S1包括以下子步骤:The step S1 includes the following sub-steps: S101.按照下式确定扫描仪有效测量距离ρ:S101. Determine the effective measurement distance ρ of the scanner according to the following formula: N=1/(ρ×h×((tan(θ+dθ)-tanθ)×dα);N=1/(ρ×h×((tan(θ+dθ)-tanθ)×dα); 其中,θ=arccos(h/ρ),N为最小点云密度,h为仪器中心距河床的高度,ρ为仪器的有效测量距离,dα,dθ分别为设置的扫描仪水平分辨率和竖直分辨率;Among them, θ=arccos(h/ρ), N is the minimum point cloud density, h is the height from the center of the instrument to the river bed, ρ is the effective measurement distance of the instrument, dα, dθ are the horizontal resolution and vertical resolution of the scanner, respectively. resolution; S102.利用卫星图像获取所测河段的几何形状,量取河段宽度、长度数据,计算扫描仪架站的最大间距LmaxS102. Use satellite images to obtain the geometric shape of the measured river reach, measure the width and length data of the river reach, and calculate the maximum distance L max of the scanner stand;
Figure FDA0003023067520000011
Figure FDA0003023067520000011
其中,d为河段宽度的1/2;Among them, d is 1/2 of the width of the reach; 依据架站最大间距在图像上进行架站布置,站点在两岸交替布设,以避免仪器扫描盲区;According to the maximum distance between the racks and the stations, the racks are arranged on the image, and the stations are alternately arranged on both sides to avoid the blind spot of instrument scanning; S103.在各个预设的站点位置上架设仪器,不同两测站间布置最少3个标靶;然后在不同站点进行扫描,直到完成整个区域的测量,获取满足冲淤变形测量要求和完整覆盖扫描区域的点云数据。S103. Set up instruments at each preset site position, and arrange at least 3 targets between different two stations; then scan at different stations until the measurement of the entire area is completed, and the requirements for scouring and deposition deformation measurement and complete coverage scanning are obtained. The point cloud data of the area.
2.根据权利要求1所述的一种基于三维激光扫描的河道冲淤测量方法,其特征在于:所述步骤S2包括以下子步骤:2. The method for measuring scour and sedimentation of a river channel based on three-dimensional laser scanning according to claim 1, wherein the step S2 comprises the following sub-steps: S201.将第一期点云数据导入,进行点云去噪,植被的去除;以第一站的扫描仪坐标系为工程坐标系,第一站中的靶标坐标为工程坐标,首先,通过求解下式计算坐标转换参数:S201. Import the point cloud data of the first phase, and perform point cloud denoising and vegetation removal; take the scanner coordinate system of the first station as the engineering coordinate system, and the target coordinates in the first station as the engineering coordinates. First, by solving The following formula calculates the coordinate transformation parameters: A=RA'+ΔTA=RA'+ΔT 式中:where:
Figure FDA0003023067520000021
Figure FDA0003023067520000021
R包含三个独立参数α、β、
Figure FDA0003023067520000022
分别为绕x,y,z轴的旋转角;
R contains three independent parameters α, β,
Figure FDA0003023067520000022
are the rotation angles around the x, y, and z axes, respectively;
Figure FDA0003023067520000023
ΔX、ΔY、ΔZ为三个平移参数;
Figure FDA0003023067520000023
ΔX, ΔY, ΔZ are three translation parameters;
Figure FDA0003023067520000024
为靶标的工程坐标;
Figure FDA0003023067520000024
is the engineering coordinates of the target;
Figure FDA0003023067520000025
为靶标的扫描仪坐标;
Figure FDA0003023067520000025
is the scanner coordinates of the target;
联立求解得到六参数α、β、
Figure FDA0003023067520000026
ΔX、ΔY、ΔZ,完成第二站到工程坐标系的初配准,接着采用最近点迭代算法减小配准误差,实现最终配准;并重复进行第3,4,...,n站的配准;
Simultaneously solve the six parameters α, β,
Figure FDA0003023067520000026
ΔX, ΔY, ΔZ, complete the initial registration of the second station to the engineering coordinate system, and then use the closest point iterative algorithm to reduce the registration error to achieve the final registration; and repeat the 3rd, 4th, ..., n stations registration;
S202.导入第二期点云,进行去噪与去植被,在同一工程中导入第一期配准好的点云,删去河道内的点云,保留河道外部地形未发生变化的点;接着,采用最近点迭代算法,利用两期点云重合的河道外部区域,将第二期第一站点云向第一期点云进行配准;按照同样的方法完成第二期剩下点云的配准;S202. Import the point cloud of the second phase, perform denoising and devegetation, import the point cloud registered in the first phase in the same project, delete the point cloud in the river channel, and keep the points where the terrain outside the river channel has not changed; then , using the closest point iteration algorithm to register the first site cloud of the second phase to the point cloud of the first phase by using the outer area of the river where the point clouds of the two phases overlap; complete the registration of the remaining point clouds of the second phase in the same way allow; S203.将配准后的两期点云,删去河道外区域点云,并将剩余点云按凸多边形区域保留,以txt格式分别导出点云坐标信息,用日期命名。S203. Delete the point cloud in the area outside the river channel from the registered two-phase point cloud, and keep the remaining point cloud as a convex polygon area, export the point cloud coordinate information in txt format, and name them with a date.
CN201910686756.9A 2019-07-29 2019-07-29 A method for measuring erosion and sedimentation of river channels based on 3D laser scanning Active CN110390687B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910686756.9A CN110390687B (en) 2019-07-29 2019-07-29 A method for measuring erosion and sedimentation of river channels based on 3D laser scanning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910686756.9A CN110390687B (en) 2019-07-29 2019-07-29 A method for measuring erosion and sedimentation of river channels based on 3D laser scanning

Publications (2)

Publication Number Publication Date
CN110390687A CN110390687A (en) 2019-10-29
CN110390687B true CN110390687B (en) 2021-06-01

Family

ID=68287747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910686756.9A Active CN110390687B (en) 2019-07-29 2019-07-29 A method for measuring erosion and sedimentation of river channels based on 3D laser scanning

Country Status (1)

Country Link
CN (1) CN110390687B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7385445B2 (en) * 2019-11-27 2023-11-22 東急建設株式会社 Concrete surface evaluation system
CN111833401B (en) * 2020-06-12 2022-05-27 华中科技大学 Rapid ranging method and system based on double-sided point cloud information
CN112461496B (en) * 2020-09-09 2023-08-18 福建省水利水电勘测设计研究院有限公司 Visual processing method for hydraulic and river model dredging test results
CN112762904B (en) * 2020-12-25 2023-09-19 长沙理工大学 Rapid measurement equipment and method for shoreline dredging evolution in beach repair test
CN112880582B (en) * 2021-01-19 2023-01-17 河海大学 A monitoring device for local scour-silting deformation of tidal flat bed in intertidal zone
CN112685918B (en) * 2021-02-04 2024-09-27 中建一局华江建设有限公司 River channel dredging method based on unmanned aerial vehicle inclination measurement technology
CN113160374B (en) * 2021-05-25 2022-07-15 西安科技大学 Three-dimensional calculation method for volume change of gully based on terrain point cloud
CN113465524B (en) * 2021-05-26 2023-07-28 中国水利水电第七工程局有限公司 A real-time dam surface deformation monitoring method for earth-rockfill dam filling based on point cloud data
CN115014701B (en) * 2022-05-10 2023-05-05 浙江省水利河口研究院(浙江省海洋规划设计研究院) Integrated measuring device and method for underwater topography in indoor flushing test
CN118031904B (en) * 2024-04-12 2024-06-25 四川省公路规划勘察设计研究院有限公司 Expressway tunnel clearance measurement method and device based on vehicle-mounted laser point cloud

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607447A (en) * 2012-03-08 2012-07-25 北京北科安地科技发展有限公司 Method for rapidly monitoring deformation by aid of ground-based three-dimensional laser scanner
CN103852025A (en) * 2014-03-19 2014-06-11 北京工业大学 Method for monitoring vertical deformation in rail way underlying substratum by applying 3D laser scanning technology
CN105205855A (en) * 2015-09-14 2015-12-30 西北农林科技大学 Method for measuring water and soil loss conditions
CN106225707A (en) * 2016-08-01 2016-12-14 三峡大学 A kind of method for the deformation of fast monitored high CFRD extrusion side wall
CN106679579A (en) * 2016-12-02 2017-05-17 中国电建集团昆明勘测设计研究院有限公司 Mobile landslide mass deformation monitoring device and method
CN107316299A (en) * 2017-07-13 2017-11-03 云南数云信息科技有限公司 Ancient architecture deformed monitoring method and system based on three-dimensional point cloud technology

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9619691B2 (en) * 2014-03-07 2017-04-11 University Of Southern California Multi-view 3D object recognition from a point cloud and change detection
US10319101B2 (en) * 2016-02-24 2019-06-11 Quantum Spatial, Inc. Systems and methods for deriving spatial attributes for imaged objects utilizing three-dimensional information
CN109375281A (en) * 2018-09-05 2019-02-22 上海华测导航技术股份有限公司 Goaf surveying method based on 3 D laser scanning

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102607447A (en) * 2012-03-08 2012-07-25 北京北科安地科技发展有限公司 Method for rapidly monitoring deformation by aid of ground-based three-dimensional laser scanner
CN103852025A (en) * 2014-03-19 2014-06-11 北京工业大学 Method for monitoring vertical deformation in rail way underlying substratum by applying 3D laser scanning technology
CN105205855A (en) * 2015-09-14 2015-12-30 西北农林科技大学 Method for measuring water and soil loss conditions
CN106225707A (en) * 2016-08-01 2016-12-14 三峡大学 A kind of method for the deformation of fast monitored high CFRD extrusion side wall
CN106679579A (en) * 2016-12-02 2017-05-17 中国电建集团昆明勘测设计研究院有限公司 Mobile landslide mass deformation monitoring device and method
CN107316299A (en) * 2017-07-13 2017-11-03 云南数云信息科技有限公司 Ancient architecture deformed monitoring method and system based on three-dimensional point cloud technology

Also Published As

Publication number Publication date
CN110390687A (en) 2019-10-29

Similar Documents

Publication Publication Date Title
CN110390687B (en) A method for measuring erosion and sedimentation of river channels based on 3D laser scanning
Mukupa et al. A review of the use of terrestrial laser scanning application for change detection and deformation monitoring of structures
CN101788666B (en) Underwater three dimensional terrain reconstruction method based on multi-beam sonar data
CN109960838B (en) An automatic generation method of river terrain reflecting the basic characteristics of rivers
CN113465524B (en) A real-time dam surface deformation monitoring method for earth-rockfill dam filling based on point cloud data
CN102494611B (en) Method for rapidly measuring volume of object
CN113160374B (en) Three-dimensional calculation method for volume change of gully based on terrain point cloud
CN102609940A (en) Method for processing errors generated by point cloud registration in process of surface reconstruction of measuring object by using ground laser scanning technique
CN103278115A (en) Method and system for calculating deposition volume of check dam based on DEM (digital elevation model)
CN111879300A (en) Method for monitoring collapse erosion development based on three-dimensional laser scanning technology
CN109631786B (en) Three-dimensional laser scanning method for surface deformation of similar materials in underground engineering simulation test
CN115861527A (en) Method and device for constructing live-action three-dimensional model, electronic equipment and storage medium
CN104331927A (en) Method and system for simulating digital outcrop surface of reservoir
CN111257870A (en) An underwater terrain inversion method for coal mining subsidence ponded area using InSAR monitoring data
CN107122522B (en) Digital elevation model precision evaluation method based on fill-cut analysis
CN115909091A (en) Earth volume calculation method based on unmanned aerial vehicle three-dimensional scanning live-action modeling
Cheng et al. Study on reverse engineering of historical architecture based on 3D laser scanner
CN104613945A (en) Reconstruction method for terrain of shallow-sea large-sized complicated sand wave area
CN112734929B (en) Method for calculating excavation volume of complex earth and rockfill dam earth stock ground based on grid subdivision algorithm
Lo Brutto et al. Image-based and range-based 3D modelling of archaeological cultural heritage: the Telamon of the Temple of Olympian Zeus in Agrigento (Italy)
CN117576334B (en) Underwater live-action three-dimensional data base plate acquisition method based on digital twin technology
Alba et al. Geometric modelling of a large dam by terrestrial laser scanning
CN110287560A (en) A Calculation Method of Excavation Quantity for Complex Shape Soil Field
CN102426010B (en) Method for generating variable resolution digital elevation model
US12038277B2 (en) High-precision waterway reconstruction method based on multi-satellite source information coupling

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
GR01 Patent grant
GR01 Patent grant