CN109141376B - Monocular vision-based wave direction detection method - Google Patents
Monocular vision-based wave direction detection method Download PDFInfo
- Publication number
- CN109141376B CN109141376B CN201810883025.9A CN201810883025A CN109141376B CN 109141376 B CN109141376 B CN 109141376B CN 201810883025 A CN201810883025 A CN 201810883025A CN 109141376 B CN109141376 B CN 109141376B
- Authority
- CN
- China
- Prior art keywords
- image
- straight line
- pixel
- point
- position point
- 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
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention relates to a method for estimating wave direction by utilizing monocular images of waves. Firstly, setting a position point path traversed by an image, selecting straight lines with different slopes at each position point, fitting a straight line according to pixel value data of pixel points superposed in each straight line and the image, and then calculating the size of a pixel value change period by using the pixel value data and the fitted straight lines. And then, calculating the slope of the straight line corresponding to the minimum pixel period obtained by calculation in all the position points to obtain an average value, wherein the average value is the wave slope of the waves in the image. The method provided by the invention can obtain better accuracy in wave direction estimation and has the advantages of small influence of an image shooting angle, and can improve the water surface environment perception capability of a ship when the ship sails in a complex water area.
Description
Technical Field
The invention belongs to the technical field of water surface environment detection, and particularly relates to a monocular vision wave direction detection method.
Background
When a ship sails in the wind and waves, the sailing resistance encountered is increased. When the ship is sailed by top waves, because the collision period of the surge and the ship is shortened, the collision frequency of the surge is increased, the collision degree is intensified, and the damage to the ship body is increased; when the ship sails perpendicular to the direction of the surge propagation, the generated rolling can cause the life raft, the life boat, the anchor and other decks to be easily bound, loosened and separated by animals, and can fall into the sea when serious; surge resistance can also overload and slow down the marine main engine, can load up electromechanical equipment, and can cause abnormal conditions to occur in such equipment. The above situations all bring difficulties to the ship control and affect the navigation safety.
In order to solve the above problems, the ship needs to know some hydrological information of its own water area in real time to deal with the interference caused by the wind waves, that is, the ship needs to have a certain sensing capability of the water surface environment, and the wave direction detection is an indispensable part of the sensing capability of the water surface environment of the ship. Common wave direction detection methods comprise a threshold segmentation method, a Fourier transform method and a Radon transform method, and the methods have several defects, namely, due to the complex algorithm and large calculation amount, the instantaneity is insufficient; secondly, the method is easily influenced by illumination or shooting angles, and the identification effect is not accurate enough; thirdly, the image to be detected needs high resolution and is high in cost. Therefore, in view of the above disadvantages, there is a need for a wave direction detection method with strong real-time performance and anti-interference of illumination and shooting angle.
Disclosure of Invention
The technical problem to be solved by the present invention is to provide a wave direction detection method based on monocular vision by using a computer image processing technology, according to the phenomenon that in a monocular image, a wave has a strong texture characteristic, the pixel gray value of a wave peak relative to a wave trough is high, any one straight line direction in the image must pass through a plurality of wave peaks and wave troughs, the pixel gray value from the wave peak to the wave trough is periodically changed from high to low in each straight line direction, and the period of the change is shortest in the wave direction.
In order to solve the above technical problems, the present invention provides a wave direction detection method based on monocular vision, which is characterized by comprising the following steps:
step 1: set point of position
Let the size of the image F shot by the monocular camera be M lines, N columns, the ith position point biPosition in the image (x)i,yi) Comprises the following steps:
wherein i is 1,2, … n; n is N/5, and the total number of the position points is obtained.
Step 2: recording pixel values
Let the current position point be biTaking the position point as the center, 36 straight lines with an angle interval of 5 degrees are selected, and the jth line passes through the position point bi(xi,yi) Straight line Y ofjThe equation of (a) is:
Yj=tan(5j)(Xj-xi)+yi (2)
wherein j is 1,2, … 36; xj=1,2,…N;
Straight line YjPixel value H of pixel point coinciding with image FwComprises the following steps:
Hw=F(Xj,Yj′) (3)
wherein w is 1,2, … wj,wjIs a straight line YjThe number of pixels coincident with the image F; y isj' is YjAn integer rounded up; f (x, y) is the pixel value of image F at the (x, y) location.
And step 3: fitting straight line and calculating pixel period
Using least squares method, according to HwFitting a straight line by the pixel value data in the step (b), finding all the zero-crossing points by an up-crossing zero point method (namely, the point which is intersected with the fitted straight line when the pixel value rises is taken as the up-crossing zero point), and recording the difference value of corresponding abscissas of every two adjacent zero-crossing points as tpWherein p is 1,2, …, pj;pjThe number of the upper zero crossing points is reduced by 1.
All t are comparedpSorting the values in descending order, discarding the value a percent after the sorting value, and leaving tpIs denoted by tq(ii) a Wherein q is 1,2, …, qj;qj=pj(100-a); and a is a fault tolerance ratio coefficient, and a is 1-50, and an empirical value is selected. Then at biThe pixel period in the jth straight line direction at the position point is as follows:
wherein i is the ith position point, i is 1,2, …, n; j is the jth line, j being 1,2, …, 36. Step 4: drawing wave-direction lines
Recording the slope of the line with the minimum pixel period in all the lines at the ith position point as KiThe slope of the wave lineComprises the following steps:
in the image F taken by the monocular camera, the wave direction line equation representing the wave direction is:
compared with the prior art, the method has higher real-time performance, can quickly acquire the wave direction of the target water area, and is slightly influenced by the shooting angle of the camera; wave direction interference information on an appointed navigation path can be timely and effectively provided for the ship in navigation; the method provided by the invention can obtain better accuracy in wave direction estimation and can improve the water surface environment perception capability of the ship.
The invention also has the advantages of low detection cost and low dependence of the detection device.
Description of the drawings:
fig. 1 is a wave direction detection flow chart based on monocular vision according to the present invention.
The specific implementation mode is as follows:
the following is a preferred embodiment of the present invention and is further described with reference to the accompanying drawings, but the present invention is not limited to the embodiment.
As shown in fig. 1, the monocular vision-based wave direction detection method of the present invention includes the following steps:
step 1: set point of position
The size of the image F shot by monocular vision is M rows, N columns, the ith position point biPosition in the image (x)i,yi) Comprises the following steps:
wherein i is 1,2, … n; n is N/5, and the total number of the position points is obtained.
Step 2: recording pixel values
Let the current position point be biTaking the position point as the center, 36 straight lines with an angle interval of 5 degrees are selected, and the jth line passes through the position point bi(xi,yi) Straight line Y ofjThe equation of (a) is:
Yj=tan(5j)(Xj-xi)+yi (2)
wherein j is 1,2, … 36; xj=1,2,…N;
Straight line YjPixel value H of pixel point coinciding with image FwComprises the following steps:
Hw=F(Xj,Yj′) (3)
wherein w is 1,2, … wj,wjIs a straight line YjThe number of pixels coincident with the image F; y isj' is YjAn integer rounded up; f (x, y) is the pixel value of image F at the (x, y) location.
And step 3: fitting straight line and calculating pixel period
Using least squares method, according to HwFitting a straight line by the pixel value data in the step (b), finding all the zero-crossing points by an up-crossing zero point method (namely, the point which is intersected with the fitted straight line when the pixel value rises is taken as the up-crossing zero point), and recording the difference value of corresponding abscissas of every two adjacent zero-crossing points as tpWherein p is 1,2, …, pj;pjThe number of the upper zero crossing points is reduced by 1.
All t are comparedpSorting the values in descending order, discarding the value a percent after the sorting value, and leaving tpIs denoted by tq(ii) a Wherein q is 1,2, …, qj;qj=pj(100-a); and a is a fault tolerance ratio coefficient, and a is 1-50, and an empirical value is selected. Then at biThe pixel period in the jth straight line direction at the position point is as follows:
wherein i is the ith position point, i is 1,2, …, n; j is the jth line, j is 1,2, …, 36;
and 4, step 4: drawing wave-direction lines
Recording the slope of the line with the minimum pixel period in all the lines at the ith position point as KiThe slope of the wave lineComprises the following steps:
in the image F taken by the monocular camera, the wave direction line equation representing the wave direction is:
Claims (1)
1. a wave direction detection method based on monocular vision is characterized by comprising the following steps:
step 1: set point of position
Let the size of the image F shot by the monocular camera be M lines, N columns, the ith position point biPosition in the image (x)i,yi) Comprises the following steps:
wherein i is 1,2, … n; n is N/5 and is the total number of the position points;
step 2: recording pixel values
Let the current position point be biTaking the position point as the center, 36 straight lines with an angle interval of 5 degrees are selected, and the jth line passes through the position point bi(xi,yi) Straight line Y ofjThe equation of (a) is:
Yj=tan(5j)(Xj-xi)+yi (2)
wherein j is 1,2, … 36; xj1,2, … N; straight line YjPixel value H of pixel point coinciding with image FwComprises the following steps:
Hw=F(Xj,Yj′) (3)
wherein w is 1,2, … wj,wjIs a straight line YjThe number of pixels coincident with the image F; y isj' is YjAn integer rounded up; f (x, y) is the pixel value of image F at the (x, y) location;
and step 3: fitting straight line and calculating pixel period
Using least squares method, according to HwFitting a straight line by the pixel value data in the step (2), setting the point crossed with the fitted straight line when the pixel value rises as an ascending zero point, finding all the ascending zero points, and recording the difference value of the corresponding abscissa of every two adjacent ascending zero points as tpWherein p is 1,2, …, pj;pjThe number of the upper span zero points is reduced by 1;
all t are comparedpSorting the values in descending order, discarding the value a percent after the sorting value, and leaving tpIs denoted by tq(ii) a Wherein q is 1,2, …, qj;qj=pj(100-a); a is a fault-tolerant ratio coefficient, a is 1-50, and an empirical value is selected; then at biThe pixel period in the jth straight line direction at the position point is as follows:
wherein i is the ith position point, i is 1,2, …, n; j is the jth line, j is 1,2, …, 36;
and 4, step 4: drawing wave-direction lines
Recording the slope of the line with the minimum pixel period in all the lines at the ith position point as KiThe slope of the wave lineComprises the following steps:
in the image F taken by the monocular camera, the wave direction line equation representing the wave direction is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810883025.9A CN109141376B (en) | 2018-08-06 | 2018-08-06 | Monocular vision-based wave direction detection method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810883025.9A CN109141376B (en) | 2018-08-06 | 2018-08-06 | Monocular vision-based wave direction detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109141376A CN109141376A (en) | 2019-01-04 |
CN109141376B true CN109141376B (en) | 2021-02-26 |
Family
ID=64791571
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810883025.9A Active CN109141376B (en) | 2018-08-06 | 2018-08-06 | Monocular vision-based wave direction detection method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109141376B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101586956A (en) * | 2009-06-18 | 2009-11-25 | 上海交通大学 | River water level monitoring method based on monocular camera |
CN201748922U (en) * | 2010-07-07 | 2011-02-16 | 南京信息工程大学 | Wind wave element value live-action monitoring system |
CN102721411A (en) * | 2012-06-28 | 2012-10-10 | 上海海事大学 | Wave scale monitoring method based on water wave image |
CN102878985A (en) * | 2012-06-27 | 2013-01-16 | 上海海事大学 | Water surface wave scale monitoring method based on image texture features |
WO2017179344A1 (en) * | 2016-04-11 | 2017-10-19 | 古野電気株式会社 | Wave height calculating device, radar device, and wave height calculating method |
-
2018
- 2018-08-06 CN CN201810883025.9A patent/CN109141376B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101586956A (en) * | 2009-06-18 | 2009-11-25 | 上海交通大学 | River water level monitoring method based on monocular camera |
CN201748922U (en) * | 2010-07-07 | 2011-02-16 | 南京信息工程大学 | Wind wave element value live-action monitoring system |
CN102878985A (en) * | 2012-06-27 | 2013-01-16 | 上海海事大学 | Water surface wave scale monitoring method based on image texture features |
CN102721411A (en) * | 2012-06-28 | 2012-10-10 | 上海海事大学 | Wave scale monitoring method based on water wave image |
WO2017179344A1 (en) * | 2016-04-11 | 2017-10-19 | 古野電気株式会社 | Wave height calculating device, radar device, and wave height calculating method |
Non-Patent Citations (2)
Title |
---|
"关于依上跨零点波高和依下跨零点波高的统计分析";范顺庭等;《海洋科学》;19841231(第4期);1-7 * |
"基于单目视觉的波高检测方法研究";倪文军,王晨;《水道港口》;20151231;第36卷(第6期);591-595 * |
Also Published As
Publication number | Publication date |
---|---|
CN109141376A (en) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108765458B (en) | Sea surface target scale self-adaptive tracking method of high-sea-condition unmanned ship based on correlation filtering | |
CN106384344B (en) | A kind of remote sensing image surface vessel target detection and extracting method | |
Chen et al. | Ship detection from coastal surveillance videos via an ensemble Canny-Gaussian-morphology framework | |
US9031285B2 (en) | Detection of floating objects in maritime video using a mobile camera | |
CN114782905B (en) | Ship draft detection method based on video monitoring | |
CN108509919B (en) | Method for detecting and identifying waterline in video or picture based on deep learning | |
CN104217215A (en) | Classification and identification method for foggy water surface image and clear water surface image | |
CN113112540B (en) | Method for positioning ship image target by using AIS (automatic identification system) Calibration CCTV (CCTV) camera in VTS (video tape server) system | |
CN113177593B (en) | Fusion method of radar point cloud and image data in water traffic environment | |
CN110619328A (en) | Intelligent ship water gauge reading identification method based on image processing and deep learning | |
CN111292371A (en) | Intelligent water level detection equipment, intelligent water level monitoring management device and detection method | |
CN116168028B (en) | High-speed rail original image processing method and system based on edge filtering under low visibility | |
CN116824570B (en) | Draught detection method based on deep learning | |
CN111738071B (en) | Inverse perspective transformation method based on motion change of monocular camera | |
CN109141376B (en) | Monocular vision-based wave direction detection method | |
CN111105419B (en) | Vehicle and ship detection method and device based on polarized SAR image | |
Shi et al. | Object detection based on saliency and sea-sky line for USV vision | |
CN117630911A (en) | Unmanned ship radar photoelectric fusion method based on Singer model | |
Kaur et al. | Sea situational awareness (seasaw) dataset | |
Wei et al. | Automatic water line detection for an USV system | |
CN112115891B (en) | SAR image port detection method based on one-dimensional multidirectional extraction | |
CN107941220B (en) | Unmanned ship sea antenna detection and navigation method and system based on vision | |
CN114994705A (en) | Unmanned ship obstacle avoidance method based on laser radar angle correction | |
CN108198214A (en) | The extraction of sea drift rubbish and area computation method based on video monitoring | |
CN105785990B (en) | Ship mooring system and obstacle recognition method based on panoramic looking-around |
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 |