[go: up one dir, main page]

CN111443344B - Automatic extraction method and device for side-scan sonar sea bottom line - Google Patents

Automatic extraction method and device for side-scan sonar sea bottom line Download PDF

Info

Publication number
CN111443344B
CN111443344B CN202010263405.XA CN202010263405A CN111443344B CN 111443344 B CN111443344 B CN 111443344B CN 202010263405 A CN202010263405 A CN 202010263405A CN 111443344 B CN111443344 B CN 111443344B
Authority
CN
China
Prior art keywords
echo data
starboard
port
scan sonar
seabed
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
CN202010263405.XA
Other languages
Chinese (zh)
Other versions
CN111443344A (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 Defense Technology Innovation Institute PLA Academy of Military Science
Original Assignee
National Defense Technology Innovation Institute PLA Academy of Military Science
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 Defense Technology Innovation Institute PLA Academy of Military Science filed Critical National Defense Technology Innovation Institute PLA Academy of Military Science
Priority to CN202010263405.XA priority Critical patent/CN111443344B/en
Publication of CN111443344A publication Critical patent/CN111443344A/en
Application granted granted Critical
Publication of CN111443344B publication Critical patent/CN111443344B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523Details of pulse systems
    • G01S7/526Receivers
    • G01S7/527Extracting wanted echo signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8902Side-looking sonar
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a method and a device for automatically extracting a side-scan sonar sea bottom line. The method comprises the following steps: acquiring a blind area boundary of echo data; extracting the port and starboard echo data of all Ping scanning lines according to the width of the blind zone and the maximum sampling number of the port and starboard echo data, and eliminating abnormal points of the port and starboard echo data; detecting the coordinates of the seabed points by using the port and starboard echo data for eliminating the abnormal points; and extracting a seabed line according to the seabed point coordinates. The method can overcome strong noise interference without manual intervention to set the initial value of the bottom line, and has high accuracy rate of the bottom line extraction and good detail continuity.

Description

Automatic extraction method and device for side-scan sonar sea bottom line
Technical Field
The invention relates to the technical field of underwater acoustic signal processing, in particular to a method and a device for automatically extracting a side-scan sonar sea bottom line.
Background
Side Scan Sonar (SSS) records the intensity of submarine backscattered echoes by using an echo depth sounding principle, and generates sonar images with different brightness according to the echo intensity to detect submarine landforms and underwater objects. In the side scan sonar image, an obvious boundary line between a water column area and an image area is called as a sea bottom line, which represents the distance from the side scan sonar to the sea bottom and is an important parameter for measuring the sea bottom landform or underwater objects, correcting the slant distance and equalizing the image gray level.
At present, the working mode of the side scan sonar is mainly embodied by adopting towed measurement information of towed fish. For example, the seabed line extraction may be performed by a maximum amplitude method and a gradient method. However, in order to improve the submarine line extraction effect under strong noise interference, the method needs manual intervention and sets the initial value of the submarine line.
In addition, aiming at the problem that the image in the side-scan sonar water column area is interfered by transmitted pulses, sea surface echoes, wake flows, large-area suspended matters and the like, a sea bottom line final peak value detection method and a sea bottom line repairing method are adopted, and the automatic tracking and extraction of the sea bottom line under the influence of complex ocean noise are realized. However, the algorithm is complex in calculation process, and the continuity of extracting the details of the sea bottom line is poor.
In addition, a comprehensive method integrating the technologies of seabed tracking, slant range correction, towed fish homing calculation, geocoding, gap filling and the like is provided. However, in the data processing of the method, the existence of a middle blind area of the sonar transducer array is ignored, and the precision is not high.
Disclosure of Invention
The invention aims to provide a method and a device for automatically extracting a submarine line by side scan sonar, so that stronger noise interference is overcome under the condition of not needing manual intervention to set an initial value of the submarine line, and the submarine line is high in extraction accuracy and good in detail continuity.
According to one aspect of the invention, the side-scan sonar sea bottom line extraction method comprises the following steps:
acquiring a blind area boundary of echo data;
extracting the port and starboard echo data of all Ping scanning lines according to the width of the blind zone and the maximum sampling number of the port and starboard echo data, and eliminating abnormal points of the port and starboard echo data;
detecting the coordinates of the seabed points by using the port and starboard echo data for eliminating the abnormal points; and
and extracting a seabed line according to the seabed point coordinates.
According to a further embodiment of this aspect, the acquiring a blind zone boundary of the echo data includes: and respectively detecting the first echo data of each Ping scanning line on the port and the starboard, which is smaller than the blind zone boundary detection threshold value, and determining the blind zone width of the echo data according to the bilateral symmetry principle.
According to a further embodiment of this aspect, median filtering is used to remove outliers of the port and starboard echo data.
Further, the eliminating the outlier of the port and starboard echo data by using median filtering includes:
calculating the signal intensity median of each echo data of each Ping scanning line and a plurality of adjacent echo data of the front and rear adjacent Ping scanning lines in the port and starboard echo data;
calculating the absolute deviation of the signal intensity of each echo data of each Ping scanning line, and determining the signal intensity value of each echo data according to the absolute deviation of the signal intensity; and
and according to the relation between the absolute deviation of the signal intensity and a set value, taking the signal intensity median value or the signal intensity sampling value as the signal intensity numerical value of the echo data of each line.
According to a further embodiment of this aspect, the detecting subsea point coordinates using port and starboard echo data with outlier points removed comprises:
matching and aligning two space-time data sequences of the port echo data and the starboard echo data of each Ping scanning line;
taking logarithm of the matched and aligned port and starboard echo data;
calculating the standard deviation of the logarithmic sequence of the port and starboard echo data of each Ping scanning line, and constructing a comprehensive sequence of the port and starboard echo data according to the standard deviation;
carrying out maximum detection on the comprehensive sequence of the port and starboard echo data to obtain an extreme point coordinate; and
and acquiring the coordinates of the port and starboard seabed points according to the extreme point coordinates.
Further, matching alignment is performed based on Euclidean distances between data points of the two spatiotemporal data sequences.
According to a further embodiment of this aspect, the extracting the seafloor line from the seafloor point coordinates comprises: and performing smooth filtering on the coordinate sequences of the port and starboard seabed points of all Ping scanning lines to extract accurate seabed lines.
According to another aspect of the present invention, there is provided a side-scan sonar sea floor line extracting apparatus, including:
a storage unit for storing computer instructions;
a processor, communicatively coupled to the memory unit, for executing the computer instructions; when a processor executes the computer instructions, the side-scan sonar undersea line extraction method according to any one of the preceding claims is executed.
According to still another aspect of the present invention, there is provided a carrier including side-scan sonar acoustic systems arranged on left and right sides, respectively, and a side-scan sonar undersea line extraction device according to the above-described aspect. Optionally, the vector is AUV.
According to the method, by considering the existence of the blind areas among the sonar transducers in practical application, firstly, the images are calculated to extract the blind area parameters and remove the high-frequency outlier noise, then, a dynamic time warping algorithm and an extreme value detection method are introduced to obtain the seabed point, and under the condition that the seabed line initial value is set without manual intervention, stronger noise interference can be overcome, the seabed line extraction accuracy is high, and the detail continuity is good.
Other features of the present invention and advantages thereof will become apparent from the following detailed description of exemplary embodiments thereof, which proceeds with reference to the accompanying drawings.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and not to limit the invention. In the drawings:
FIG. 1 is a side scan sonar waterfall plot;
FIG. 2 is a side scan sonar image taken alone for port and starboard;
FIG. 3 is a side scan sonar waterfall image main distribution parameter;
FIG. 4 is echo data for a single Ping scanline;
FIG. 5 is a flow chart of a method for automatically extracting a side-scan sonar sea bottom line according to a first embodiment of the present invention;
FIG. 6 is a side scan sonar image with no dead zones shown on port and starboard alone;
FIG. 7 is a side scan sonar image with no dead zone shown, taken separately from port and starboard after high frequency noise cancellation;
FIG. 8 is a side scan sonar image with submarine lines marked separately on the port and starboard side and without showing dead zones;
FIG. 9 is a side scan sonar waterfall plot with sea floor lines labeled and without showing dead zones;
fig. 10 is a schematic structural view of a side-scan sonar undersea line automatic extraction device according to a second embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the invention, its application, or uses. 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 invention.
The relative arrangement of the components and steps, the numerical expressions and numerical values set forth in these embodiments do not limit the scope of the present invention unless specifically stated otherwise. Techniques, methods, and apparatus known to those of ordinary skill in the relevant art may not be discussed in detail but are intended to be part of the specification where appropriate.
The side scan sonar is an active sonar system, and the working principle is that a short sound wave pulse with a very small horizontal open angle (about 1 degree) and a very large vertical open angle is emitted to one side or two sides of the vertical direction of the course of a measuring ship, and after the pulse reaches the seabed, the pulse is continuously reflected according to the distance between the seabed and a transducer, and a sonar image with uneven gray level change is generated according to the intensity degree of a reflected signal. The information such as the change of the landform of the seabed, whether the navigation obstacle, the type of the seabed substrate and the like exist can be observed from the sonar image.
When the side scan sonar emission pulse is propagated in a water body and meets a target, the target scatters sound energy to all directions, a transducer receives backscatter echoes, and the sound energy is difficult to reach behind the target, so that a blind area can be generated. The sonar array continuously advances along with the carrier, and sonar constantly launches in the advancing process, constantly receives and forms the sonar image, and the target (the strong echo signal of target) and its shade appear in the position that corresponds on the sonar image.
The side-scan sonar system mainly comprises a side-scan sonar acoustic system, a peripheral auxiliary sensor, a data real-time acquisition processing system and the like. The transducer is a core component of a side-scan sonar acoustic system, and for example, a piezoelectric ceramic structure can be adopted, when a voltage is applied to the transmitting transducer, the physical form of the transmitting transducer is changed, an oscillating electric field generated by the transmitter is converted into mechanical deformation, and the deformation is transmitted into water to generate oscillating pressure, namely acoustic pulse, in the water. Also, a receiving transducer is used to receive the echo signal and convert the pressure change into electrical energy by detecting the acoustic pressure change. Typically, side scan sonar systems use a transmit-receive unified linear array when designing the transducer.
An Autonomous Underwater Vehicle (AUV) is a widely used unmanned marine vehicle and can carry a large number of detection devices to automatically complete marine surveying and mapping tasks. Therefore, the AUV is used as a measuring ship, the side scan sonar is installed on the AUV, and data fusion can be well carried out by means of good maneuvering performance of the AUV.
For example, Sea from Marine Sonic Technology may be used
Figure BDA0002440285370000061
And carrying out offshore measurement on the ARC SCOUT Mk-II, and collecting a plurality of Ping original data. The working parameters of the equipment are set as follows: the working frequency is 900kHz, and the slant distance is 10 m.
Fig. 1 shows a waterfall diagram obtained by inverting a gray image displayed by 256 gray levels of raw data obtained by a certain measurement. Fig. 2(a), (b) are respective port and starboard side scan sonar images. Fig. 3 is a waterfall plot of several Ping echo data after the scanning in fig. 1 starts, where a1 represents a sea bottom line, a2 represents a sea surface line, A3 represents a dead zone, a4 represents a water column region, and a5 represents a scanned image region. In the figure, darker colors represent stronger echoes, whereas whiter (lighter) colors represent weaker echoes.
In the present invention, unless otherwise specified, all the figures shown are the inverse of the gray scale image in which the original data (echo data) is displayed at 256 gray scales.
Fig. 4 schematically shows a workflow of the automatic side-scan sonar undersea line extraction method according to the first embodiment of the present invention. As shown in fig. 4, the method for automatically extracting the side-scan sonar sea bottom line comprises the following steps:
step 101, acquiring a blind zone boundary of echo data;
when extracting the sea bottom line, the prior art often makes an ideal assumption about the geometric problem of a sonar image, considers that no blind area exists in the middle of a sonar transducer array, or considers that the blind area is small and ignores the blind area. However, as can be seen from fig. 3, in the actual work of the side scan sonar, the sonar transducer emits sound waves to both sides, the sonar image has a blind area, and the blind area occupies a large space and should not be ignored. Therefore, the present invention requires first calculating and extracting the blind area parameter.
Figure 5 shows echo data for a single Ping scanline. As shown in fig. 5, the blind area is located right below the sonar, and the echo data of the Ping scan line suddenly changes from the blind area to the outside, and the numerical cliff type drops, whereby the detection of the blind area boundary can be realized. In general, echoes directly below the side-scan sonar are strong positive emission waves, and therefore, the present invention optionally employs a threshold method for blind zone boundary detection.
For a selected side scan sonar, the Ping data bit number is fixed, e.g., 8, 16, or 32bit data. If the data bit number of the side scan sonar is Nbit, a blind area boundary detection threshold D can be set according to the measurement environment conditionmComprises the following steps:
Dm=2N-2
from the blind area to the outer side, the first less than the blind area boundary detection threshold D of each Ping scanning line of the port and the starboard is respectively detectedmRespectively denoted as KLeft side of、KRight side. Because the side-scan sonar port and starboard transducers are symmetrically installed relative to the central axis of the carrier, the width K of the echo data blind area of the Ping scanning line is determined according to the bilateral symmetry principlemThe determination is as follows:
Figure BDA0002440285370000071
102, extracting the port and starboard echo data of all Ping scanning lines according to the width of the blind zone and the maximum sampling number of the port and starboard echo data, and eliminating abnormal points of the port and starboard echo data;
for each Ping scanning line, the maximum sampling number of echo data of a carrier, such as AUV left and right sides, is respectively KLeft max、KRight maxThen the port echo data [ K ] of all Ping scan lines can be extractedm,KLeft max]And starboard echo data [ K ]m,KRight max]As shown in fig. 6.
As can be seen from fig. 6, in a complex environment, such as a water column region with suspended matter, etc., the echo data has many outliers, which may interfere with the submarine line extraction and may even cause false detection. The echo abnormal point is essentially high-frequency outlier noise and needs to be removed.
According to an optional embodiment of the present invention, the value filtering removes echo outliers in the present invention, and the specific process includes:
(1) calculating the signal intensity median of each echo data of each Ping scanning line and a plurality of adjacent echo data of the front and rear adjacent Ping scanning lines in the port and starboard echo data;
(2) and calculating the absolute deviation of the signal intensity of each echo data of each Ping scanning line, and determining the signal intensity value of each echo data according to the absolute deviation of the signal intensity.
And when the absolute deviation of the signal intensity is greater than a set value, taking the median of the signal intensity as the signal intensity numerical value of each echo data. Otherwise, the outlier removing operation is not performed, and the signal strength value of each echo data is the measured value (sampling value).
The set value is a calculated value based on a signal strength measurement of echo data of each Ping scanline.
For example, for port echo data [ K ]m,KLeft max]Let [ i, j)]Index identifier, I, of jth echo sample data representing ith Ping scan lineijDenotes the [ i, j ]]Signal strength value, P, of the port echo datamaxRepresenting the total number of side scan sonar Ping scans, then:
for the satisfaction of 1<i<Pmax,Km≤j<KLeft maxAll of [ i, j ] of]And combining, namely performing the following data processing:
(1) calculate [ i-1, j-1 [ ] -1]、[i-1,j]、[i-1,j+1]、[i,j-1]、[i,j]、[i,j+1]、[i+1,j-1]、[i+1,j]、[i+1,j+1]Median Z of the port echo dataijAnd calculating the absolute deviation Eij=abs(Iij-Zij);
This example selects the jth echo data of the ith Ping scan line, and its 3 adjacent echo data of the preceding, i.e., ith-1 Ping scan line, i.e., jth-1, j +1 echo data, and its 3 adjacent echo data of the following, i.e., ith +1 Ping scan line, i.e., jth-1, j +1 echo data. However, as will be understood by those skilled in the art, the number of adjacent scan lines, or the number of echo data adjacent to each adjacent scan line, can be selected according to actual needs. Similarly, the same is true for the subsequently described numerical calculation of the signal intensity of the starboard echo data.
(2) Determine absolute deviation Eij
If E isij>9.78+1.44×Iij-8.86×10-4×Iij 2+2.95×10-7×Iij 3-3.26×10-11×Iij 2Then adopt ZijAs the [ i, j ]]A new value of the signal strength of the port echo data; otherwise, no processing is done.
For port echo data [ K ]m,KRight max]Let [ i, l ] be]Denotes the ithIndex identification of the first echo sample data of the sub-Ping scanline, IilDenotes the [ i, l ]]Signal strength values for the individual starboard echo data.
For the satisfaction of 1<i<Pmax,Km≤l<KRight maxAll of [ i, l ] of]And combining, namely performing the following data processing:
(1) calculation of [ i-1, l-1 [ ]]、[i-1,l]、[i-1,l+1]、[i,l-1]、[i,l]、[i,l+1]、[i+1,l-1]、[i+1,l]、[i+1,l+1]Median Z of starboard echo datailAnd calculating the absolute deviation Eil=abs(Iil-Zil);
(2) Determine absolute deviation Eil
If E isil>9.78+1.44×Iil-8.86×10-4×Iil 2+2.95×10-7×Iil 3-3.26×10-11×Iil 2Then adopt ZilAs the (i, l)]A new value of the signal strength of the starboard echo data; otherwise, no processing is done.
Fig. 7 shows a side scan sonar image with no dead zone shown, port and starboard alone, after high frequency outlier noise cancellation. Comparing fig. 6, it can be seen that the echo abnormal point in the side-scan data water column region can be eliminated to some extent. And for the data processing, the echo abnormal point eliminating effect of the port echo data is better.
103, detecting the coordinates of the seabed points by using the port and starboard echo data for eliminating the abnormal points;
according to the sea bottom line bilateral symmetry principle, the first echoes received by the side-scan sonar port and starboard transducers are from the sea bottom right below, experience depths are the same, and the sea bottom lines extracted from the two sides are symmetrical about a navigation track line. By utilizing the principle, the seabed point detection can be better carried out.
According to an alternative embodiment of the invention, the subsea point detection comprises the steps of:
1031, matching and aligning two space-time data sequences of the port and starboard echo data of each Ping scanning line;
due to the same Ping scanningIn the line, the maximum sampling number K of the port and starboard echo data often appearsLeft max、KRight maxUnequal problem, as shown in fig. 4, results in I in the port and starboard echo data when l is jijAnd IilThe temporal and spatial information involved is different. In order to fully utilize the sea bottom line bilateral symmetry principle and realize comprehensive extraction comparison between a port side and a starboard side so as to remove the influence of suspended matters as much as possible, the problem is firstly solved.
In the invention, the seabed line extracts the data of the water column region and the data of the abrupt change part from the water column region to the image region, therefore, in order to keep the robustness of the invention, the left and right side echo data can be considered as two space-time data sequences to be matched and aligned. When matching and aligning, the consistency of the time-space information of the abrupt change from the water column area to the image area needs to be kept. According to the optional embodiment of the present invention, the matching alignment of the two sequences in space-time by using a dynamic time warping algorithm specifically includes the following steps:
(1) for each Ping scanline, two spatio-temporal data sequences, respectively, left side echo data { Iij,Km≤j≤KLeft maxAnd starboard echo data Iil,Km≤l≤KRight max}。
In each of the two sequences, a Euclidean distance D (I) between two points is calculatedij,Iil) In which K ism≤j≤KLeft max,Km≤l≤KRight max
Calculating Euclidean distances for all data points in the sequence, and constructing a Euclidean distance table, wherein the following table shows that:
D(IiKm,IiKm) D(IiKm,IiKm+1) D(IiKm,IiK Right max)
D(IiKm+1,IiKm) D(IiKm+1,IiKm+1) D(IiKm+1,IiK Right max)
D(IiK left max,IiKm) D(IiK left max,IiKm+1) D(IiK left max,IiK Right max)
(2) Searching the shortest path in the Euclidean distance table;
in the Euclidean distance table, D (I) is searched for, for example, by using a dynamic programming or greedy search algorithmiKm,IiKm) To D (I)iK left max,IiK Right max) The shortest path of (c). This shortest path needs to satisfy: if the current node is D (I)ij,Iil) Then the next node must be at D (I)ij+1,Iil)、D(Iij,Iil+1)、D(Iij+1,Iil+1) And the path must be shortest.
(3) According to the Euclidean distance table nodes passed by the shortest path, the aligned data point pairs of two discrete sequences can be obtained, and the sequencesThe lengths are the same and are marked as Kmax
Step 1032, taking logarithm of the matched and aligned port and starboard echo data;
the logarithm of the echo data after the high-frequency noise cancellation processing is performed on the port and starboard matching and aligning, for example, a common logarithm with a base of 10 is taken.
Step 1033, calculating a standard deviation of the logarithmic sequence of the port and starboard echo data of each Ping scanning line, and constructing a comprehensive sequence of the port and starboard echo data according to the standard deviation;
setting the window as W, the window should be selected as small as possible to avoid flooding the abrupt features of the data. Calculating the standard deviation of the logarithmic sequence of the port echo data and the starboard echo data of each Ping scanning line to obtain two standard deviation sequences, and respectively recording { L }iq,1≤q≤Kmax-W +1} and { L }is,1≤s≤Kmax-W +1}, and the sequence lengths are all Kmax-W + 1. Wherein the value of W can also be determined by reference to the following procedure:
Figure BDA0002440285370000111
for example, echo outliers existing in the water column region in the port echo data in fig. 5 may cause a large fluctuation in the standard deviation sequence of the water column region when calculating the standard deviation sequence, which may affect the accuracy of the subsequent extremum detection method.
Therefore, in order to fully utilize the sea bottom line bilateral symmetry principle, the invention further takes the minimum value of the starboard and starboard logarithmic sequence standard deviation of each Ping scanning line, namely: for any s is q, q is more than or equal to 1 and less than or equal to Kmax-W+1,1≤s≤Kmax-W +1, taking Lig=min(Liq,Lis) And then constructing a comprehensive sequence { L) of the port and starboard echo data by using the minimum valueig,1≤g≤Kmax-W+1}。
1034, carrying out maximum detection on the comprehensive sequence of the port and starboard echo data to obtain an extreme point coordinate;
according to the characteristics of the side-scan sonar data, sea can be knownThe bottom line is composed of the first submarine strong echo sequence from port to starboard in each Ping scanning line, so that sudden change of echo data from the water column region to the image region can result in a port and starboard comprehensive sequence { L }ig,1≤g≤Kmax-W +1 }. Thus, by synthesizing sequences { L ] for port and starboardig,1≤g≤Kmax-W +1} to obtain an extreme point coordinate gb
And 1035, obtaining the left and right ship bottom point coordinates according to the extreme point coordinates.
Using extreme point coordinates gbThe coordinate q of the standard sequence of the logarithm of the port echo data and the starboard echo data can be correspondingly obtainedbAnd sbAnd then the coordinate q of the echo data sequence after the port and starboard matching and aligning corresponding to the extreme point can be obtainedb+ W-1 and corresponding echo data Ib left side、sb+ W-1 and corresponding echo data Ib right side
For each Ping scanline, search for the port echo data sequence Iij,Km≤j≤KLeft maxFind the first one and ILeft of bEqual echo data, the echo data coordinate is the port sea bottom point coordinate jb. Likewise, starboard echo data sequence { I ] is searchedil,Km≤l≤KRight maxFind the first one and Ib right sideEqual echo data, the echo data coordinate is starboard seabed point coordinate lb
And 104, extracting a seabed line according to the seabed point coordinates.
The coordinate sequences { j of the port seabed points of all Ping scanning lines obtained by calculationbAnd starboard sea-bottom point coordinate sequence (l)bAnd performing smooth filtering, namely extracting an accurate seabed line.
The side-scan sonar image which is obtained by adopting the extraction method disclosed by the invention and is obtained by marking the sea bottom line on the left and right sides independently and does not display the blind area is shown in figure 8, and the side-scan sonar waterfall map which is marked with the sea bottom line and does not display the blind area is shown in figure 9.
According to a second embodiment of the present invention, there is also provided an automatic side-scan sonar undersea line extraction apparatus, as shown in fig. 10, including:
a storage unit 10 for storing computer instructions;
a processor 20, communicatively coupled to the memory unit 10, for executing the computer instructions. When the processor 20 executes the computer instructions, the automatic extraction method of the side-scan sonar undersea line described above is performed.
According to a third embodiment of the present invention, there is also provided a carrier for mounting a side scan sonar system.
The carrier comprises side-scan sonar acoustic systems which are respectively arranged on a left board and a right board, and a side-scan sonar sea bottom line extraction device.
The vector may be AUV.
The above examples are only intended to illustrate the technical solution of the present invention and not to limit it; although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that: modifications to the specific embodiments of the invention or equivalent substitutions for parts of the technical features may be made; without departing from the spirit of the present invention, it is intended to cover all aspects of the invention as defined by the appended claims.

Claims (8)

1. A side scan sonar sea bottom line extraction method is characterized by comprising the following steps:
acquiring a blind area boundary of echo data;
extracting the port and starboard echo data of all Ping scanning lines according to the width of the blind zone and the maximum sampling number of the port and starboard echo data, and eliminating abnormal points of the port and starboard echo data;
detecting the coordinates of the seabed points by using the port and starboard echo data for eliminating the abnormal points; and
extracting a seabed line according to the seabed point coordinates;
the method for detecting the coordinates of the seabed points by using the port and starboard echo data with the abnormal points eliminated comprises the following steps:
matching and aligning two space-time data sequences of the port and starboard echo data of each Ping scanning line; wherein the matching alignment is based on Euclidean distances between data points of the two spatiotemporal data sequences;
taking logarithm of the matched and aligned port and starboard echo data;
calculating the standard deviation of the logarithmic sequence of the port and starboard echo data of each Ping scanning line, and constructing a comprehensive sequence of the port and starboard echo data according to the standard deviation;
carrying out maximum detection on the comprehensive sequence of the port and starboard echo data to obtain extreme point coordinates; and
and acquiring the coordinates of the port and starboard seabed points according to the extreme point coordinates.
2. The side-scan sonar undersea line extraction method according to claim 1, wherein the acquiring of the dead zone boundary of the echo data includes: and respectively detecting the echo data of which the first of each Ping scanning line on the port and the starboard is smaller than a blind zone boundary detection threshold value, and determining the width of the blind zone of the echo data according to a bilateral symmetry principle.
3. The side-scan sonar sea bottom line extraction method according to claim 1, wherein outliers of the port and starboard echo data are eliminated using median filtering.
4. The side-scan sonar sea bottom line extraction method according to claim 3, wherein the removing of the outliers of the port and starboard echo data by median filtering includes:
calculating the signal intensity median of each echo data of each Ping scanning line and a plurality of adjacent echo data of the front and rear adjacent Ping scanning lines in the port and starboard echo data;
calculating the absolute deviation of the signal intensity of each echo data of each Ping scanning line, and determining the signal intensity value of each echo data according to the absolute deviation of the signal intensity; and
and according to the relation between the absolute deviation of the signal intensity and a set value, taking the signal intensity median value as a signal intensity numerical value of each echo data.
5. The side-scan sonar sea bottom line extraction method according to claim 1, wherein extracting sea bottom lines from sea bottom point coordinates includes: and performing smooth filtering on the coordinate sequences of the port and starboard seabed points of all Ping scanning lines to extract accurate seabed lines.
6. The utility model provides a side scan sonar undersea line extraction element which characterized in that, the device includes:
a storage unit for storing computer instructions;
a processor, communicatively coupled to the memory unit, for executing the computer instructions; when a processor executes the computer instructions, the side-scan sonar undersea line extraction method according to any one of claims 1-5 is performed.
7. A carrier comprising side-scan sonar acoustic systems arranged on the left and right sides, respectively, and a side-scan sonar sea-bottom line extraction device according to claim 6.
8. The vector of claim 7, wherein the vector is AUV.
CN202010263405.XA 2020-04-07 2020-04-07 Automatic extraction method and device for side-scan sonar sea bottom line Active CN111443344B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010263405.XA CN111443344B (en) 2020-04-07 2020-04-07 Automatic extraction method and device for side-scan sonar sea bottom line

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010263405.XA CN111443344B (en) 2020-04-07 2020-04-07 Automatic extraction method and device for side-scan sonar sea bottom line

Publications (2)

Publication Number Publication Date
CN111443344A CN111443344A (en) 2020-07-24
CN111443344B true CN111443344B (en) 2022-06-03

Family

ID=71651204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010263405.XA Active CN111443344B (en) 2020-04-07 2020-04-07 Automatic extraction method and device for side-scan sonar sea bottom line

Country Status (1)

Country Link
CN (1) CN111443344B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112882037B (en) * 2021-04-28 2021-07-23 北京星天科技有限公司 Side-scan sonar sea bottom line detection method and device
CN114627367B (en) * 2022-05-17 2022-07-26 国家海洋局北海海洋技术保障中心 Sea bottom line detection method for side-scan sonar image
CN116930976B (en) * 2023-06-19 2024-03-26 自然资源部第一海洋研究所 Submarine line detection method of side-scan sonar image based on wavelet mode maximum value
CN117557894B (en) * 2023-10-19 2025-02-11 中科探海(深圳)海洋科技有限责任公司 A method for detecting and tracking seabed positions based on sonar image features
CN118655554B (en) * 2024-08-21 2024-11-22 中国船舶集团有限公司第七〇七研究所 Frequency shift calculation method and system for non-beam-separated phased array Doppler speedometer

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4063212A (en) * 1976-05-19 1977-12-13 Western Marine Electronics, Inc. Side scan sonar system
CN1402018A (en) * 2001-09-13 2003-03-12 中国科学院声学研究所 High resolution submarine microgeomorphy-measuring sounding side scan sonar system and measuring method
CN104536005A (en) * 2014-11-17 2015-04-22 哈尔滨工程大学 Multi-beam side-scan sonar slant-distance distortion elimination method based on blind area correction
CN106500671A (en) * 2016-09-22 2017-03-15 天津大学 A kind of based on LM algorithms decompose laser radar waveform determine sea water advanced method
CN108320291A (en) * 2018-01-25 2018-07-24 中国科学院深海科学与工程研究所 Submarine geomorphy boundary extraction method based on Otsu algorithm
CN108872997A (en) * 2018-05-10 2018-11-23 国家海洋局第二海洋研究所 A kind of seabed line detecting method based on side scan sonar data fusion and accurate processing
CN108957462A (en) * 2018-05-22 2018-12-07 中国海洋大学 A kind of multi-beam water body data processing method based on smooth bottom
CN110675410A (en) * 2019-09-25 2020-01-10 江苏海洋大学 Side-scan sonar sunken ship target unsupervised detection method based on selective search algorithm

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4063212A (en) * 1976-05-19 1977-12-13 Western Marine Electronics, Inc. Side scan sonar system
CN1402018A (en) * 2001-09-13 2003-03-12 中国科学院声学研究所 High resolution submarine microgeomorphy-measuring sounding side scan sonar system and measuring method
CN104536005A (en) * 2014-11-17 2015-04-22 哈尔滨工程大学 Multi-beam side-scan sonar slant-distance distortion elimination method based on blind area correction
CN106500671A (en) * 2016-09-22 2017-03-15 天津大学 A kind of based on LM algorithms decompose laser radar waveform determine sea water advanced method
CN108320291A (en) * 2018-01-25 2018-07-24 中国科学院深海科学与工程研究所 Submarine geomorphy boundary extraction method based on Otsu algorithm
CN108872997A (en) * 2018-05-10 2018-11-23 国家海洋局第二海洋研究所 A kind of seabed line detecting method based on side scan sonar data fusion and accurate processing
CN108957462A (en) * 2018-05-22 2018-12-07 中国海洋大学 A kind of multi-beam water body data processing method based on smooth bottom
CN110675410A (en) * 2019-09-25 2020-01-10 江苏海洋大学 Side-scan sonar sunken ship target unsupervised detection method based on selective search algorithm

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种联合单波束测深的侧扫声呐海底线提取方法;库安邦等;《海洋通报》;20190228;第38卷(第1期);47-52 *
侧扫声呐图像海底线自动提取方法研究;赵建虎等;《武汉大学学报 信息科学版》;20171231;第42卷(第12期);1797-1803 *
基于LOG 算子的侧扫声呐海底线检测;张济博等;《海洋通报》;20100630;第29卷(第3期);324-328 *

Also Published As

Publication number Publication date
CN111443344A (en) 2020-07-24

Similar Documents

Publication Publication Date Title
CN111443344B (en) Automatic extraction method and device for side-scan sonar sea bottom line
US6438071B1 (en) Method for producing a 3D image
US20120099400A1 (en) Estimating position and orientation of an underwater vehicle relative to underwater structures
WO2003023446A1 (en) High resolution bathymetric sonar system and measuring method for measuring the physiognomy of the seabed
CN103477244A (en) Detecting structural changes to underwater structures
CN114488164A (en) Underwater vehicle synchronous positioning and mapping method and underwater vehicle
KR102263037B1 (en) A Method of Underwater Environment Mapping System using Underwater Vehicle and Underwater Acoustic Detection Equipment
JP2002168952A (en) Method of reconstituting submarine three-dimensional structure
Yu et al. Bottom Detection Method of Side‐Scan Sonar Image for AUV Missions
JP7259503B2 (en) Sonar image processing device, sonar image processing method and program
CN213813947U (en) Offshore wind power foundation submarine cable access end detection system
JP2008076294A (en) Under-bottom-of-water survey method and instrument
CN116523822A (en) Submarine cable detection and identification method based on side-scan sonar
Hellequin et al. Postprocessing and signal corrections for multibeam echosounder images
CN113777596B (en) Adaptive fish finder gain compensation method and device based on sea depth detection
US11668821B2 (en) Position correction using towed sensor
JPH0385476A (en) Sea bottom searching apparatus
JP2010145224A (en) Scanning sonar device
CN116930976B (en) Submarine line detection method of side-scan sonar image based on wavelet mode maximum value
Fietz et al. Background noise cancellation from active sonar images utilizing transmission-free recordings
US11846742B1 (en) Systems and methods for the localization of objects buried in the seabed
Chen et al. Single ping filtering of multi-beam sounding data based on alpha shapes
US11947006B1 (en) Systems and methods for the localization of objects buried in the seabed
CN119620051A (en) Detection method for side-scan sonar seabed line
Tang et al. A Multibeam Underwater Terrain Modeling Algorithm Based on Echo Energy Distribution Gradient

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