CN103364787B - A kind of multi-beam side-scan sonar image mosaic fissure removing method - Google Patents
A kind of multi-beam side-scan sonar image mosaic fissure removing method Download PDFInfo
- Publication number
- CN103364787B CN103364787B CN201310104485.4A CN201310104485A CN103364787B CN 103364787 B CN103364787 B CN 103364787B CN 201310104485 A CN201310104485 A CN 201310104485A CN 103364787 B CN103364787 B CN 103364787B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mover
- rightarrow
- mfrac
- rsqb
- 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
Abstract
The invention belongs to acoustic picture to splice field, and in particular to a kind of multi-beam side-scan sonar image mosaic fissure removing method based on gps coordinate and course angle amendment.The present invention includes:Read half-tone information;Draw original track line;Establish model;Recording angular difference calculates Δ angle [x] average per 10Ping, establishes correction model;Obtain seamless landform image.The present invention can obtain more accurately positional information, reduce multi-beam side-scan sonar image and inlay caused gap, improve syncretizing effect.
Description
Technical field
The invention belongs to acoustic picture to splice field, and in particular to a kind of more ripples based on gps coordinate and course angle amendment
Beam side-scanning sonar image mosaic fissure removing method.
Background technology
Increasingly reduced with the resource of landing field, the mankind have turned to wide ocean space for naturally exploitation nature, sea
The exploration of foreign resource, the description of marine territory, the research in exclusive economic zone are finally required for the ocean map to be in a variety of manners
According to criterion.Current marine charting technology depends on sonar detection means, but because sonar GPS location is inaccurate and abnormal
The destabilizing factors such as saltus step are present, and result in and crack is inevitably present after side-scanning sonar image is inlayed, therefore can not obtain
Accurate sea-floor relief image.
Sonar Mosaic gap is repaired, and exactly reads the gps coordinate and course angle information in original sonar data, is carried out
Corresponding processing, eliminate the change of the asynchronous and speed of a ship or plane of influence and sound sonar GPS location scans record with gathering to(for) sonar
Receive due to the flexible delay of tow strap and the influence of caused course angle mutation, finally give seamless bottom mounted sonar mosaic image.
Related product in the market is external company's exploitation, and expensive, core technology is underground, and the country there is no related production
Product, and existing sonar mosaic image gap removing method only applies the pixel in interpolation filling chink, not only without fully profit
With the data in sonar, and it is computationally intensive, display effect is bad, therefore, to the multi-beam based on gps coordinate and course angle amendment
Side-scanning sonar image mosaic fissure removing method has great importance.
The content of the invention
It is an object of the invention to provide a kind of precision based on gps coordinate and course angle amendment is higher, syncretizing effect is more
Good multi-beam side-scan sonar image mosaic fissure removing method.
The object of the present invention is achieved like this:
The present invention comprises the following steps:
(1)The half-tone information of each passage in multi-beam side-scan sonar file is read, shows sonar image;
(2)GPS latitude and longitude coordinates information, ship's speed information and course angle information are successively read, draws original track line;
(3)When track points are that Tracepoint [k] arrives Tracepoint [k+n], gps coordinate is different from collection due to positioning
Step does not change, and at Tracepoint [k+2], speed changes to V2 by V1, establishes model:
Wherein, Tracepoint is the array of storage track points information, and k is that transducer launches sound wave Ping numbers, and n is whole
Number, m=k+2, Displacement respectively before and after velocity variations, Tracepoint [m-1]~Tracepoint [m+1] are exactly weight
Track points after building;
(4)Course angle heading [x], the ship's head angle shipheading [x] of towfish are recorded, note Δ angle [x] is
Both differential seat angles, i.e.,
Δangle[x]=shipheading[x]-heading[x];
(5)Calculate Δ angle [x] average per 10Ping, be designated as averageangle, towfish between adjacent two Ping and
Ship's head angular difference value is designated as Δ 1 and Δ 2:
Δ1=heading[x]-heading[x-1]
Δ2=shipheading[x]-shipheading[x-1];
(6)Establish correction model:
Wherein, α is corrected parameter;
(7)According to revised flight path and course angle information embedding image, seamless landform image is obtained.
The beneficial effects of the present invention are:It is right for multi-beam side scan sonar system GPS and course angle data acquisition feature
The relation that velocity information changes between GPS information carries out accurate position correction;According to towfish course angle in itself and hull
Relation between course angle, course angle information is modified.So as to obtain more accurately positional information, multi-beam side is reduced
Gap caused by Sonar Mosaic is swept, improves syncretizing effect.
Brief description of the drawings
Fig. 1 is the basic flow sheet of the inventive method;
Fig. 2 is the mosaic image for not carrying out filling up in gap;
Fig. 3 is the image that picture element interpolation filling chink is utilized after inlaying;
Fig. 4 is the mosaic image after being handled using the inventive method.
Embodiment
Illustrate below in conjunction with the accompanying drawings and the present invention is described in more detail:
It is characteristic of the invention that:
1. step(3)In when being corrected to GPS position information, using velocity information as foundation, when GPS changes,
Its position distribution is not equidistant arrangement, but adds speed factor, and positional information is more defined closer to virtual condition, display
Really.
2. when multi-beam side scan sonar system carries out data acquisition, system can record two course angle informations, come respectively
From towfish and hull, the present invention utilizes this feature, carries out course angle correction.Actual acquired data comes from towfish, therefore should choose
Towfish course angle is as processing foundation, but when ship is turned, and because fluid retention acts on, its towfish angle change is jumped
Property it is larger, it is relatively stable from the change of ship's head angle, so according to step(4)、(5)、(6)Middle course angle correction model, knot
Towfish and hull course angle information is closed to be modified course angle information, so as to get course angle information it is more accurate, more meet sea
Bottom data characteristicses.
With reference to Fig. 1, the present invention is eliminated based on the multi-beam side-scan sonar image mosaic fissure of gps coordinate and course angle amendment
Method comprises the following steps that:
(1)The half-tone information of each passage in multi-beam side-scan sonar file is read, shows sonar image;
(2)GPS latitude and longitude coordinates information, ship's speed information and course angle information are successively read, draws original track line;
(3)When arriving Tracepoint [k+n] according to track points for Tracepoint [k](Remember that Tracepoint navigates for storage
The array of mark point information, k are that transducer launches sound wave Ping numbers, and n is usually 4), gps coordinate is due to GPS location and sonar data
Gather it is asynchronous do not change, but speed changes to V2 by V1 at Tracepoint [k+2], and we carry out data and shown
When, overlapping gps coordinate should carry out stretching amendment, i.e., be uniformly distributed the GPS information between this n PING therebetween, if speed around here
Information changes, then its distribution should be directly proportional to speed, therefore establishes such as drag:
Wherein, m=k+2, Displacement respectively before and after velocity variations, Tracepoint [m-1]~Tracepoint
[m+1] is exactly the track points after rebuilding;
(4)Because delay of the towfish in water acts on, when ship steady steaming, between towfish course angle and hull course angle
In the presence of certain stable difference, course angle heading [x], the ship's head angle shipheading [x] of towfish are recorded, remembers Δ
Angle [x] is both differential seat angles, i.e.,
Δangle[x]=shipheading[x]-heading[x];
(5)Δ angle [x] average per 10Ping is calculated, is designated as averageangle, is(4)In the stabilization mentioned
Course angular difference, the influence of towfish course angle mutation is excluded, has calculated towfish and ship's head angular difference value between adjacent two Ping,
Δ 1 and Δ 2 are designated as respectively:
Δ1=heading[x]-heading[x-1]
Δ2=shipheading[x]-shipheading[x-1];
(6)Establish following correction model:
(7)According to revised flight path and course angle information embedding image, seamless landform image is obtained.
(8)Figure of description 2 is the mosaic image for not carrying out filling up in gap, and Fig. 3 is to be filled out after inlaying using picture element interpolation
The image for gap of squeegeeing, it can be seen that after the existing processing using picture element interpolation method, sawtooth effect of image is obvious, and fill
Pixel has obscured display effect, and Fig. 4 is the mosaic image after being handled using the inventive method, can see by contrast,
After being handled with the inventive method, sonar data is fully used so that mosaic image is more accurate, better into scheming.
Claims (1)
1. a kind of multi-beam side-scan sonar image mosaic fissure removing method, it is characterised in that comprise the following steps:
(1) half-tone information of each passage in multi-beam side-scan sonar file is read, shows sonar image;
(2) GPS latitude and longitude coordinates information, ship's speed information and course angle information are successively read, draws original track line;
(3) when track points are that Tracepoint [k] arrives Tracepoint [k+n], gps coordinate is asynchronous simultaneously due to positioning and collection
Do not change, at Tracepoint [k+2], speed changes to V2 by V1, establishes model:
<mrow>
<mover>
<mi>P</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<mover>
<mrow>
<mi>V</mi>
<mn>1</mn>
</mrow>
<mo>&RightArrow;</mo>
</mover>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<mover>
<mrow>
<mi>V</mi>
<mn>2</mn>
</mrow>
<mo>&RightArrow;</mo>
</mover>
<mo>*</mo>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mo>(</mo>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>T</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mi>p</mi>
<mi>o</mi>
<mi>i</mi>
<mi>n</mi>
<mi>t</mi>
<mo>&lsqb;</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mfrac>
<mrow>
<mrow>
<mo>|</mo>
<mover>
<mi>&delta;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
</mrow>
<mo>*</mo>
<mover>
<mi>P</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
<mrow>
<mover>
<mi>P</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>+</mo>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
</mfrac>
<mo>+</mo>
<mi>T</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mi>p</mi>
<mi>o</mi>
<mi>int</mi>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>T</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mi>p</mi>
<mi>o</mi>
<mi>i</mi>
<mi>n</mi>
<mi>t</mi>
<mo>&lsqb;</mo>
<mi>m</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mrow>
<mrow>
<mo>|</mo>
<mover>
<mi>&delta;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
</mrow>
<mo>*</mo>
<mover>
<mi>P</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
<mrow>
<mover>
<mi>P</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>+</mo>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
</mfrac>
<mo>+</mo>
<mi>T</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mi>p</mi>
<mi>o</mi>
<mi>int</mi>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mi>T</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mi>p</mi>
<mi>o</mi>
<mi>int</mi>
<mo>&lsqb;</mo>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mfrac>
<mrow>
<mrow>
<mo>|</mo>
<mover>
<mi>&delta;</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
</mrow>
<mo>*</mo>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
<mrow>
<mover>
<mi>P</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>+</mo>
<mover>
<mi>Q</mi>
<mo>&RightArrow;</mo>
</mover>
</mrow>
</mfrac>
<mo>+</mo>
<mi>T</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mi>p</mi>
<mi>o</mi>
<mi>int</mi>
<mo>&lsqb;</mo>
<mi>m</mi>
<mo>&rsqb;</mo>
</mrow>
Wherein, Tracepoint is the array of storage track points information, and k is that transducer launches sound wave Ping numbers, and n is integer, m=
K+2,Displacement respectively before and after velocity variations, Tracepoint [m-1]~Tracepoint [m+1] are exactly after rebuilding
Track points;
(4) course angle heading [x], the ship's head angle shipheading [x] of towfish are recorded, note Δ angle [x] is both
Differential seat angle, i.e.,
Δ angle [x]=shipheading [x]-heading [x];
(5) Δ angle [x] average per 10Ping is calculated, is designated as averageangle, the course of the towfish between adjacent two Ping
Angular difference value is designated as Δ 1, and the ship's head angular difference value between adjacent two Ping is designated as Δ 2:
Δ 1=heading [x]-heading [x-1]
Δ 2=shipheading [x]-shipheading [x-1];
(6) correction model is established:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>di</mi>
<mi>n</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>s</mi>
<mi>h</mi>
<mi>i</mi>
<mi>p</mi>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>d</mi>
<mi>i</mi>
<mi>n</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>r</mi>
<mi>a</mi>
<mi>g</mi>
<mi>a</mi>
<mi>n</mi>
<mi>g</mi>
<mi>l</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mn>1</mn>
<mo><</mo>
<mi>&Delta;</mi>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>d</mi>
<mi>i</mi>
<mi>n</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>h</mi>
<mi>e</mi>
<mi>a</mi>
<mi>d</mi>
<mi>i</mi>
<mi>n</mi>
<mi>g</mi>
<mo>&lsqb;</mo>
<mi>x</mi>
<mo>-</mo>
<mn>1</mn>
<mo>&rsqb;</mo>
<mo>+</mo>
<mi>&alpha;</mi>
<mo>*</mo>
<mi>&Delta;</mi>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mn>1</mn>
<mo>></mo>
<mi>&Delta;</mi>
<mn>2</mn>
<mo>;</mo>
<mn>0</mn>
<mo><</mo>
<mi>&alpha;</mi>
<mo><</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
Wherein, α is corrected parameter;
(7) according to revised flight path and course angle information embedding image, seamless landform image is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310104485.4A CN103364787B (en) | 2013-03-28 | 2013-03-28 | A kind of multi-beam side-scan sonar image mosaic fissure removing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310104485.4A CN103364787B (en) | 2013-03-28 | 2013-03-28 | A kind of multi-beam side-scan sonar image mosaic fissure removing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103364787A CN103364787A (en) | 2013-10-23 |
CN103364787B true CN103364787B (en) | 2017-11-28 |
Family
ID=49366558
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310104485.4A Active CN103364787B (en) | 2013-03-28 | 2013-03-28 | A kind of multi-beam side-scan sonar image mosaic fissure removing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103364787B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592731A (en) * | 2009-07-09 | 2009-12-02 | 浙江大学 | A kind of side-scan sonar towfish flight path disposal route based on the track line file |
US8213740B1 (en) * | 2009-05-18 | 2012-07-03 | The United States Of America, As Represented By The Secretary Of The Navy | Coherent image correlation |
CN102622732A (en) * | 2012-03-14 | 2012-08-01 | 上海大学 | Front-scan sonar image splicing method |
CN102707289A (en) * | 2012-06-15 | 2012-10-03 | 哈尔滨工程大学 | Course angle rotation-based multi-beam side-scan sonar image real-time splicing method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008508539A (en) * | 2004-08-02 | 2008-03-21 | ジョンソン アウトドアズ インコーポレイテッド | Sonar imaging system for ship |
-
2013
- 2013-03-28 CN CN201310104485.4A patent/CN103364787B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8213740B1 (en) * | 2009-05-18 | 2012-07-03 | The United States Of America, As Represented By The Secretary Of The Navy | Coherent image correlation |
CN101592731A (en) * | 2009-07-09 | 2009-12-02 | 浙江大学 | A kind of side-scan sonar towfish flight path disposal route based on the track line file |
CN102622732A (en) * | 2012-03-14 | 2012-08-01 | 上海大学 | Front-scan sonar image splicing method |
CN102707289A (en) * | 2012-06-15 | 2012-10-03 | 哈尔滨工程大学 | Course angle rotation-based multi-beam side-scan sonar image real-time splicing method |
Non-Patent Citations (4)
Title |
---|
"A side scan sonar image denoising algorithm based on compound of fuzzy weighted average and kalman filter";X Ye等;《International conference on mechatronics and automation》;20121231;第720-724页 * |
"Sidescan sonar image processing techniques";P Cervenka等;《IEEE journal of oceanic engineering》;19931231;第18卷(第2期);第108-122页 * |
"Sidescan sonar image processing:correcting brightness variation and patching gaps";YC Chang等;《Journal of marine science and technology》;20101231;第18卷(第6期);第785-789页 * |
"一种基于航迹线文件的声纳图像裂缝处理方法";尹天鹤等;《武汉大学学报信息科学版》;20090831;第34卷(第8期);第898-901页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103364787A (en) | 2013-10-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105352476B (en) | Boat-carrying waterfront line underwater integrated measuring system integrated approach waterborne | |
CN106886024B (en) | Deep-sea multi-beam sound ray precise tracking method | |
US9651698B2 (en) | Multi-beam bathymetric chart construction method based on submarine digital depth model feature extraction | |
CN103456041B (en) | A kind of dimensional topography based on S-57 electronic chart data and radar terrain generation method | |
Bistacchi et al. | Photogrammetric digital outcrop reconstruction, visualization with textured surfaces, and three-dimensional structural analysis and modeling: Innovative methodologies applied to fault-related dolomitization (Vajont Limestone, Southern Alps, Italy) | |
CN105241457A (en) | Establishing method of three-dimensional aided navigation system for ship handling | |
US10520631B2 (en) | Magnetic field measurement via streamer cables | |
Menna et al. | A photogrammetric approach to survey floating and semi-submerged objects | |
CN105488852B (en) | A kind of 3-D image joining method based on geocoding and multidimensional calibrating | |
CN102707289B (en) | Based on the real-time joining method of multi-beam side-scan sonar image that course angle rotates | |
CN103292792A (en) | Actual measurement SVP reconstruction method suitable for submarine detection and pseudo-landform processing | |
CN109031319A (en) | A kind of side-scanning sonar image splicing system and its method | |
CN104764445A (en) | Method and device for determining coordinates of underwater object point | |
CN110095786A (en) | Three-dimensional point cloud based on a line laser radar ground drawing generating method and system | |
Chemisky et al. | Underwater survey for oil and gas industry: A review of close range optical methods | |
CN103592650A (en) | Three-dimensional sonar imaging system based on graph processor and three-dimensional image method thereof | |
CN106197398A (en) | It is applied to navigation system and the map constructing method thereof of River Channel | |
CN110706177B (en) | Method and system for equalizing gray level of side-scan sonar image | |
CN101592731B (en) | Method for processing navigation track of side scan sonar tow-fish based on navigation track line file | |
US20170219731A1 (en) | Method of visualizing and interpreting wide azimuth profile (wap) | |
KR102339500B1 (en) | Method for providing recommended sea route based on electronic nautical chart by measuring topography and depth of water | |
CN109782355A (en) | The detection method and device of OBS detection point drift | |
CN103364787B (en) | A kind of multi-beam side-scan sonar image mosaic fissure removing method | |
Bu et al. | Improved calibration method for refraction errors in multibeam bathymetries with a wider range of water depths | |
CN101290223A (en) | Sea-floor relief measurement beam angle effect correction method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |