CN111047704B - Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm - Google Patents

Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm Download PDF

Info

Publication number
CN111047704B
CN111047704B CN201911036196.9A CN201911036196A CN111047704B CN 111047704 B CN111047704 B CN 111047704B CN 201911036196 A CN201911036196 A CN 201911036196A CN 111047704 B CN111047704 B CN 111047704B
Authority
CN
China
Prior art keywords
point
water depth
data
value
points
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
CN201911036196.9A
Other languages
Chinese (zh)
Other versions
CN111047704A (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.)
NATIONAL MARINE DATA AND INFORMATION SERVICE
Original Assignee
NATIONAL MARINE DATA AND INFORMATION SERVICE
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 NATIONAL MARINE DATA AND INFORMATION SERVICE filed Critical NATIONAL MARINE DATA AND INFORMATION SERVICE
Priority to CN201911036196.9A priority Critical patent/CN111047704B/en
Publication of CN111047704A publication Critical patent/CN111047704A/en
Application granted granted Critical
Publication of CN111047704B publication Critical patent/CN111047704B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/187Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention belongs to the field of post-processing of marine mapping data, in particular to an automatic multi-beam sounding data rough difference removing method for improving a region growing algorithm, which comprises the following steps of: analyzing data: analyzing the original multi-Beam data, outputting standardized data of Ping number, beam number, longitude, latitude and water depth, and setting a threshold value for data rejection; step 2: constructing a Delaunay triangulation network: according to longitude, latitude and water depth, adopting a Delaunay growth method to construct an irregular triangular network TIN; step 3: selecting a seed point: the first step of region growth is the selection of seed points, and the invention selects the most reliable water depth point as the seed point by taking the submarine topography as a standard when selecting the seed point; step 4: removing depth measurement data gross errors by using an improved region growing algorithm: starting to perform regional growth from the seed points, and judging whether the water depth is an abnormal value according to a judging method; step 5: outputting normal water depth points: deleting the data marked as the abnormal water depth point, and outputting the normal water depth point.

Description

Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm
Technical Field
The invention belongs to the field of marine mapping data post-processing, and particularly relates to a multi-beam sounding data rough error automatic clearing method for improving a region growing algorithm.
Background
The traditional multi-beam outlier clearing method mainly comprises two modes of manual interaction and automatic algorithm.
Mature commercial software such as Caris, fledermaus mostly adopts a manual interaction method to remove abnormal values, and the process usually needs to consume a great deal of manpower and time, so that the efficiency is low; the currently applied automatic algorithm mainly comprises median filtering, wavelet analysis, trend surface filtering and the like, a part of filtering methods are easy to delete real terrain data by mistake, most of filtering methods are based on a coordinate grid and filter along or perpendicular to the direction of a navigation path, correlation of each water depth point is not fully considered, and complexity of submarine topography cannot be well adapted.
Therefore, a method for automatically removing the data rough difference should be designed to address the above-mentioned shortcomings.
Disclosure of Invention
The invention aims at overcoming the defects of the prior art and provides a multi-beam sounding data rough difference automatic clearing method of an improved region growing algorithm applying a region growing algorithm and a sounding data point topological relation established by Delaunay.
The technical scheme adopted by the invention is as follows:
the multi-beam sounding data coarse error automatic clearing method for improving the region growing algorithm is characterized by comprising the following steps of:
step 1: analyzing data: analyzing the original multi-Beam data, outputting standardized data of Ping number, beam number, longitude, latitude and water depth, and setting a threshold value for data rejection;
step 2: constructing a Delaunay triangulation network: according to longitude, latitude and water depth, adopting a Delaunay growth method to construct an irregular triangular network TIN;
step 3: selecting a seed point: the first step of region growth is the selection of seed points, and the invention selects the most reliable water depth point as the seed point by taking the submarine topography as a standard when selecting the seed point;
step 4: removing depth measurement data gross errors by using an improved region growing algorithm: starting to perform regional growth from the seed points, and judging whether the water depth is an abnormal value according to a judging method;
step 5: outputting normal water depth points: deleting the data marked as the abnormal water depth point, and only outputting the normal water depth point.
Further, in the step 1, the threshold value for setting the data rejection is preferably 2.5% of the water depth value.
Further, in the step 2, the construction of the irregular triangular net TIN includes the following steps:
2.1 optional Point P in the discrete Water depth Point set 1
2.2 selecting distance P in a discrete Water depth Point set 1 Nearest point P 2 Forming a first datum edge;
2.3 on the right side of the first reference edge, finding a third point P according to the maximum opening angle criterion of Delaunay triangulation as the third point discriminant 3 Namely P i As shown in formula (1), the first two points are respectively connected to form a first Delaunay triangle;
2.4 searching other points capable of constructing Delaunay triangle by taking three sides of the initial triangle as reference until all discrete sampling points are traversed, namely forming the irregular triangle network TIN.
Further, in the step 3, the selecting of the seed points adopts the following method, including:
3.1, traversing all triangles of which each point participates in construction in the Delaunay triangle, and determining all vertexes of the triangles of which the point participates in construction;
3.2 calculating the standard deviation of the water depth values of the point and all vertexes of the triangle participating in construction, wherein n is the number of vertexes and x is shown in the formula (2) i Is the value of the water depth of the top point,is the average value of the water depth of the vertex. The smaller standard deviation shows that the more continuous the topography change of the point is, the point is the true water depth value and can be used as the seed point for the regional growth
Further, in the step 4, the method for judging the abnormal value of the water depth includes:
4.1 concentrating at discrete Water depth points, finding the point P closest to the seed point P 1 Forming a first reference edge;
4.2 finding the third vertex P of the triangle to be constructed on the right side of the reference edge according to the maximum opening angle criterion 2
4.3 calculating the Point P respectively 2 Water depth value of (2) and P 1 If the difference value of the water depth values exceeds the threshold value, the point is considered to be an abnormal value to be deleted when the difference value of the water depth values exceeds the threshold value, and the abnormal value is marked as a rough difference point; otherwise, the point is considered to be a normal water depth point to be reserved;
4.4 at P P 2 Or P 1 P 2 To continue the region growth for the reference, if P 2 Judging the water depth point to be a normal water depth point according to the step 4.3; if P 2 Is an abnormal water depth point, then according to P and P 1 Calculating P using inverse distance weighted interpolation 2 As shown in formula (3), wherein Z i D is the water depth value of the abnormal water depth point i Is the distance between the abnormal point and the other two vertexes of the triangle. When P P 2 Finding the third point P of the triangle for reference continued region growth 3 When according to the modified P 2 And judging whether the water depth value and the P point water depth value are abnormal water depth points or not.
4.5 repeating the steps 4.3 and 4.4 until the water depth point cannot find the point meeting the maximum opening angle criterion.
The invention has the advantages and positive effects that:
the method aims at automatically removing the abnormal value of the multi-beam, combines the topological relation constructed by the Delaunay triangulation network, fully considers the mutual relevance of all the water depth points, improves the region growing algorithm by editing the water depth value of the data gross error for the first time, and realizes the automatic removal of the multi-beam sounding data gross error. Experimental results prove that the multi-beam sounding data rough difference can be accurately identified and removed, the current situation of removing the multi-beam sounding abnormal value based on the grid and the path line is changed, and compared with the processing mode based on the grid, the correlation of each discrete point can be fully considered, and the multi-beam sounding abnormal value removing method is better suitable for the complexity of submarine topography.
Drawings
FIG. 1 is a flow chart of an automatic rough error removal method according to the present invention;
FIG. 2 is a schematic diagram of improved area growth;
FIG. 3 is a partial enlarged view of a triangular network generated by experimental data construction;
fig. 4 is a graph showing the comparison of point clouds and topography before and after the rough cleaning of experimental data in the present invention.
Detailed Description
The invention will now be further illustrated by reference to the following examples, which are intended to be illustrative, not limiting, and are not intended to limit the scope of the invention.
A multi-beam sounding data rough error automatic clearing method for improving a region growing algorithm is characterized in that: the method comprises the following steps that the flow chart is shown in fig. 1:
step 1: analyzing data: analyzing the original multi-Beam data, outputting standardized data of Ping number, beam number, longitude, latitude and water depth, and setting a threshold value for data rejection according to 2.5% of the water depth value;
step 2: constructing a Delaunay triangulation network: according to longitude, latitude and water depth, adopting a Delaunay growth method to construct an irregular triangular network TIN;
2.1 at discrete Water depth Point setsOptionally at a point P 1
2.2 selecting distance P in a discrete Water depth Point set 1 Nearest point P 2 Forming a first datum edge;
2.3 on the right side of the first reference edge, finding a third point p according to the maximum opening angle criterion of Delaunay triangulation as the third point discriminant 3 I.e. p i As shown in formula (1), the first two points are respectively connected to form a first Delaunay triangle;
2.4 searching other points capable of constructing Delaunay triangle by taking three sides of the initial triangle as reference until all discrete sampling points are traversed, namely forming the irregular triangle network TIN.
Step 3: selecting a seed point: the first step of region growth is the selection of seed points, and the invention selects the most reliable water depth point as the seed point by taking the submarine topography as a standard when selecting the seed point;
3.1, traversing all triangles of which each point participates in construction in the Delaunay triangle, and determining all vertexes of the triangles of which the point participates in construction;
3.2 calculating the standard deviation of the water depth values of the point and all vertexes of the triangle participating in construction, wherein as shown in the formula (2), the smaller the standard deviation is, the more continuous the topography change of the point is, the point is the real water depth value and can be used as the seed point for region growth
Step 4: removing depth measurement data gross errors by using an improved region growing algorithm: starting to perform regional growth from the seed points, and judging whether the water depth is an abnormal value according to a judging method;
in the step 4, the judgment of the abnormal value of the water depth adopts the following method that:
4.1 in the separationConcentrated water scattering deep points, finding the point P nearest to the seed point P 1 Forming a first reference edge;
4.2 finding the third vertex P of the triangle to be constructed on the right side of the reference edge according to the maximum opening angle criterion 2
4.3 calculating the Point P respectively 2 Water depth value of (2) and P 1 If the difference value of the water depth values exceeds the threshold value, the point is considered to be an abnormal value to be deleted when the difference value of the water depth values exceeds the threshold value, and the abnormal value is marked as a rough difference point; otherwise, the point is considered to be a normal water depth point to be reserved;
4.4 at P P 2 Or P 1 P 2 To continue the region growth for the reference, if P 2 Judging the water depth point to be a normal water depth point according to the step 4.3; if P 2 Is an abnormal water depth point, then according to P and P 1 Calculating P using inverse distance weighted interpolation 2 As shown in formula (3). When in PP 2 Finding the third point P of the triangle for reference continued region growth 3 When according to the modified P 2 And judging whether the water depth value and the P point water depth value are abnormal water depth points or not.
4.5 repeating the steps 4.3 and 4.4 until the water depth point cannot find the point meeting the maximum opening angle criterion.
Step 5: outputting normal water depth points: deleting the data marked as the abnormal water depth point, and only outputting the normal water depth point.
The working process of the invention is as follows:
the use of the present invention is described in connection with the following examples,
step 1: analyzing data: analyzing the original multi-Beam sounding data, outputting Ping number, beam number, longitude, latitude and water depth, obtaining 112912 water depth points in total, setting the water depth range as [ -1169.006, -777.965], and setting the threshold as 30m.
Step 2: constructing a Delaunay triangle network, and selecting seed points: (1) 112912 discrete water depth points construct 225799 triangles, as shown in fig. 3; (2) according to the calculation result, the standard deviation of the water depth values of all the vertices of the triangle which participate in construction of the 74996 th water depth point is 0.117, and the triangle is selected as a seed point for regional growth.
Step 3: removing rough difference of sounding data by using a region growing algorithm: according to the rough difference judging method, 138 water depth points are finally marked as rough difference points.
Step 4: outputting normal water depth points: the normal water depth point which is not marked as the rough difference is output according to Ping number, beam number, longitude, latitude and water depth.
The point clouds before and after the rough difference elimination and the topography comparison chart are shown in fig. 4.

Claims (3)

1. A multi-beam sounding data rough error automatic clearing method for improving a region growing algorithm is characterized in that: the method comprises the following steps:
step 1: analyzing data: analyzing the original multi-Beam data, outputting standardized data of Ping number, beam number, longitude, latitude and water depth, and setting a threshold value for data rejection;
step 2: constructing a Delaunay triangulation network: according to longitude, latitude and water depth, adopting a Delaunay growth method to construct an irregular triangular network TIN;
step 3: selecting a seed point: the first step of the region growth is the selection of seed points, the most reliable water depth point is selected as the seed point by taking the seabed topography as a standard when the seed point is selected, and the seed point is selected by adopting the following method which comprises the following steps:
3.1, traversing all triangles of which each point participates in construction in the Delaunay triangle, and determining all vertexes of the triangles of which the point participates in construction;
3.2, calculating the standard deviation of the water depth values of the point and all vertexes of the triangle participating in construction, wherein as shown in the following formula, the smaller the standard deviation is, the more continuous the topography change of the point is, and the point is a real water depth value and can be used as a seed point for regional growth;
wherein n is the number of vertexes, x i Is the value of the water depth of the top point,is the average value of the water depth of the vertex;
step 4: removing depth measurement data gross errors by using an improved region growing algorithm: starting to perform regional growth from the seed points, and judging whether the water depth is an abnormal value according to a judging method;
the judgment of the abnormal value of the water depth adopts the following method that:
4.1 concentrating at discrete Water depth points, finding the point P closest to the seed point P 1 Forming a first reference edge;
4.2 finding the third vertex P of the triangle to be constructed on the right side of the reference edge according to the maximum opening angle criterion 2
4.3 calculating the Point P respectively 2 Water depth value of (2) and P 1 If the difference value of the water depth values exceeds the threshold value, the point is considered to be an abnormal value to be deleted when the difference value of the water depth values exceeds the threshold value, and the abnormal water depth point is marked; otherwise, the point is considered to be a normal water depth point to be reserved;
4.4 use of PP 2 Or P 1 P 2 To continue the region growth for the reference, if P 2 Judging the water depth point to be a normal water depth point according to the step 4.3; if P 2 Is an abnormal water depth point, then according to P and P 1 Calculating P using inverse distance weighted interpolation 2 As shown in formula (3), when the water depth is represented by PP 2 Finding the third point P3 of the triangle for the reference continued region growth according to the modified P 2 Judging whether the water depth value and the P point water depth value are abnormal water depth points or not;
4.5 repeating the steps 4.3 and 4.4 until the water depth point cannot find the point meeting the maximum opening angle criterion;
step 5: outputting normal water depth points: deleting the data marked as the abnormal water depth point, and only outputting the normal water depth point.
2. The method for automatically removing coarse and fine data of multi-beam sounding with improved area growing algorithm according to claim 1, wherein the method comprises the steps of: in the step 1, the threshold value for data rejection is set to be 2.5% of the water depth value.
3. The method for automatically removing coarse and fine data of multi-beam sounding with improved area growing algorithm according to claim 1, wherein the method comprises the steps of: in the step 2, the construction of the irregular triangular net TIN includes the following steps:
2.1 optional Point P in the discrete Water depth Point set 1
2.2 selecting distance P in a discrete Water depth Point set 1 Nearest point P 2 Forming a first datum edge;
2.3 on the right side of the first reference edge, finding a third point P according to the maximum opening angle criterion of Delaunay triangulation as the third point discriminant 3 Namely P i As shown in formula (1), the two points are respectively connected with each other to form an initial triangle;
2.4 searching other points capable of constructing Delaunay triangle by taking three sides of the initial triangle as reference until all discrete sampling points are traversed, namely forming the irregular triangle network TIN.
CN201911036196.9A 2019-10-29 2019-10-29 Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm Active CN111047704B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911036196.9A CN111047704B (en) 2019-10-29 2019-10-29 Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911036196.9A CN111047704B (en) 2019-10-29 2019-10-29 Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm

Publications (2)

Publication Number Publication Date
CN111047704A CN111047704A (en) 2020-04-21
CN111047704B true CN111047704B (en) 2023-12-01

Family

ID=70232022

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911036196.9A Active CN111047704B (en) 2019-10-29 2019-10-29 Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm

Country Status (1)

Country Link
CN (1) CN111047704B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112035430B (en) * 2020-07-13 2023-12-19 国家海洋信息中心 Boundary constraint self-adaptive submarine topography water depth point detection and screening method
CN112907615B (en) * 2021-01-08 2022-07-26 中国石油大学(华东) Submarine landform unit contour and detail identification method based on region growing
CN113129328B (en) * 2021-04-22 2022-05-17 中国电子科技集团公司第二十九研究所 Target hotspot area fine analysis method
CN114509050B (en) * 2021-12-31 2023-07-21 广州南方卫星导航仪器有限公司 Water depth data processing method and system
CN115291182B (en) * 2022-07-29 2023-06-06 连云港港口工程设计研究院有限公司 Filtering method for channel single-beam dense mining water depth

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102446367A (en) * 2011-09-19 2012-05-09 哈尔滨工程大学 Method for constructing three-dimensional terrain vector model based on multi-beam sonar submarine measurement data
CN104180873A (en) * 2014-09-02 2014-12-03 长江航道测量中心 Single-wave-beam depth finder water depth gross error detection and correction method and system
CN107229594A (en) * 2017-05-05 2017-10-03 深圳市建设综合勘察设计院有限公司 A kind of multibeam bathymetric data trend surface filtering method and system
CN108447125A (en) * 2018-02-28 2018-08-24 哈尔滨工程大学 A kind of Delaunay triangulation network rapid generation based on multibeam bathymetric data
CN110111375A (en) * 2019-03-29 2019-08-09 中国地质大学(武汉) A kind of Image Matching elimination of rough difference method and device under Delaunay triangulation network constraint

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102446367A (en) * 2011-09-19 2012-05-09 哈尔滨工程大学 Method for constructing three-dimensional terrain vector model based on multi-beam sonar submarine measurement data
CN104180873A (en) * 2014-09-02 2014-12-03 长江航道测量中心 Single-wave-beam depth finder water depth gross error detection and correction method and system
CN107229594A (en) * 2017-05-05 2017-10-03 深圳市建设综合勘察设计院有限公司 A kind of multibeam bathymetric data trend surface filtering method and system
CN108447125A (en) * 2018-02-28 2018-08-24 哈尔滨工程大学 A kind of Delaunay triangulation network rapid generation based on multibeam bathymetric data
CN110111375A (en) * 2019-03-29 2019-08-09 中国地质大学(武汉) A kind of Image Matching elimination of rough difference method and device under Delaunay triangulation network constraint

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Xiaodong Cui,etc.A method for multibeam seafloor terrain classification based on self-adaptive geographic classification unit.《Applied Acoustics》.2019,全文. *
张雪洋等.基于三角网滤波检验的KD 树粗差剔除法.《测绘与空间地理信息》.2017,全文. *
彭彤 ; 邹亚未 ; 周平 ; 卢凯乐 ; .基于趋势面的多波束测深粗差剔除研究.江西测绘.2016,(第02期),全文. *

Also Published As

Publication number Publication date
CN111047704A (en) 2020-04-21

Similar Documents

Publication Publication Date Title
CN111047704B (en) Multi-beam sounding data rough error automatic clearing method for improving region growing algorithm
CN106529469B (en) Unmanned aerial vehicle-mounted LiDAR point cloud filtering method based on self-adaptive gradient
CN110136072B (en) Point cloud noise removing method, denoising system, computer device and storage medium
CN111986115A (en) Accurate elimination method for laser point cloud noise and redundant data
CN109900338A (en) A kind of road surface pit slot volume measuring method and device
CN111598780B (en) Terrain adaptive interpolation filtering method suitable for airborne LiDAR point cloud
CN109816664A (en) A kind of three-dimensional point cloud dividing method and device
CN111667574B (en) Method for automatically reconstructing regular facade three-dimensional model of building from oblique photography model
CN109410265B (en) A kind of TIN filtering innovatory algorithm based on previous DEM auxiliary
CN102419436A (en) Multibeam data processing method based on total propagation error filter
CN109919955A (en) The tunnel axis of ground formula laser radar point cloud extracts and dividing method
CN108519587A (en) A kind of real-time aerial target recognizing model of movement and method for parameter estimation
CN112669333A (en) Single tree information extraction method
CN114417951A (en) Unsupervised machine learning-based automatic subdivision optimization method for ocean unstructured grid
CN115131231A (en) Complex terrain area point cloud hierarchical filtering method assisted by multi-feature clustering
CN110490800B (en) Shallow sea terrain reconstruction method
CN113703438A (en) AUV autonomous navigation path planning method for water delivery tunnel inspection
Wang et al. Bottom Tracking Method Based on LOG/Canny and the Threshold Method for Side-scan Sonar.
CN110675501A (en) Efficient conversion method from TIN to Grid model under DEM framework
CN113567968B (en) Underwater target real-time segmentation method based on shallow water multi-beam water depth data and application thereof
CN116524127A (en) Vegetation removing method for high-gradient area
Li et al. Generating 3D model of slope eroded gully based on photo reconstruction technique
CN113763280A (en) Region growing algorithm based on spatial hierarchical topological relation for point cloud denoising
CN107798339B (en) Multi-beam abnormal value data processing algorithm based on Thiessen polygon map
Yan et al. Multi-beam Data Automatic Filtering Technology

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