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 PDFInfo
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/523—Details of pulse systems
- G01S7/526—Receivers
- G01S7/527—Extracting wanted echo signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8902—Side-looking sonar
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment 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
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 usedAnd 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:
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:
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.
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)
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)
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 |
-
2020
- 2020-04-07 CN CN202010263405.XA patent/CN111443344B/en active Active
Patent Citations (8)
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)
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 |