CN113535672A - Turbidity data processing method based on autonomous underwater robot platform - Google Patents
Turbidity data processing method based on autonomous underwater robot platform Download PDFInfo
- Publication number
- CN113535672A CN113535672A CN202010296975.9A CN202010296975A CN113535672A CN 113535672 A CN113535672 A CN 113535672A CN 202010296975 A CN202010296975 A CN 202010296975A CN 113535672 A CN113535672 A CN 113535672A
- Authority
- CN
- China
- Prior art keywords
- turbidity
- underwater robot
- autonomous underwater
- turbidity data
- coordinate
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 10
- 239000013535 sea water Substances 0.000 claims abstract description 6
- 238000013499 data model Methods 0.000 claims description 22
- 238000000034 method Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 230000008447 perception Effects 0.000 claims description 2
- 230000001568 sexual effect Effects 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000011835 investigation Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000033116 oxidation-reduction process Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/211—Schema design and management
- G06F16/212—Schema design and management with details for data modelling support
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
- B63G—OFFENSIVE OR DEFENSIVE ARRANGEMENTS ON VESSELS; MINE-LAYING; MINE-SWEEPING; SUBMARINES; AIRCRAFT CARRIERS
- B63G8/00—Underwater vessels, e.g. submarines; Equipment specially adapted therefor
- B63G8/001—Underwater vessels adapted for special purposes, e.g. unmanned underwater vessels; Equipment specially adapted therefor, e.g. docking stations
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/04—Control of altitude or depth
- G05D1/06—Rate of change of altitude or depth
- G05D1/0692—Rate of change of altitude or depth specially adapted for under-water vehicles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mechanical Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Automation & Control Theory (AREA)
- Operations Research (AREA)
- Remote Sensing (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Radar, Positioning & Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The invention belongs to the field of turbidity data processing of an ocean motion platform, and particularly relates to a turbidity data processing method based on an autonomous underwater robot platform. The invention comprises the following steps: carrying out detection movement on a planned line by an autonomous underwater robot carrying water quality instrument equipment, and measuring and recording turbidity data of the water quality instrument on the line; and carrying out digital modeling according to the turbidity data of the water quality instrument, predicting the turbidity of an undetected area according to the linear relation of the turbidity data of the detected water quality instrument, and constructing a ground turbidity data map according to the detection motion of the autonomous underwater robot. The invention fuses a large amount of discrete turbidity data into a digital model, so that the turbidity data meets the linear relation, and more obvious evidence is provided for the research of seawater turbidity.
Description
Technical Field
The invention belongs to the field of turbidity data processing of an ocean motion platform, and particularly relates to a turbidity data processing method based on an autonomous underwater robot platform.
Background
With the rapid development of autonomous underwater robots, the scientific value of ocean exploration is increasingly obvious. Parameters such as sea water temperature, turbidity, oxidation-reduction potential and the like of the ocean are important indexes for ocean science investigation, and the scientific and efficient processing of the data is particularly important. The autonomous underwater robot has the characteristics of large detection range, long operation time and flexibility due to the characteristics of unmanned and cable-free operation. The autonomous underwater robot can obtain water body data of a vertical plane and can also obtain large-area plane water body data. The autonomous underwater robot for each diving operation can obtain rich water body detection data, particularly turbidity information can provide direct evidence for detecting a submarine plume hydrothermal solution area to a great extent, a traditional spatial interpolation mode aims at site discrete sampling observation, is not suitable for continuous time observation of an autonomous underwater robot platform, occupies valuable surgical investigation time, and is difficult to exert scientific investigation value maximization. The method is particularly important for fusing and predicting the turbidity data of the seabed area by establishing a digital model by utilizing the time and space dependence relationship of the motion direction of the autonomous underwater robot in the turbidity data field.
Disclosure of Invention
The invention relates to the field of processing turbidity data of hydrothermal fluid searching in the ocean floor of an autonomous underwater robot, and aims at the defects of a turbidity data processing method of a traditional autonomous underwater robot platform.
The technical scheme adopted by the invention for realizing the purpose is as follows:
a turbidity data processing method based on an autonomous underwater robot platform comprises the following steps:
carrying out detection movement on a planned line by an autonomous underwater robot carrying water quality instrument equipment, and measuring and recording turbidity data of the water quality instrument on the line;
and carrying out digital modeling according to the turbidity data of the water quality instrument, predicting the turbidity of an undetected area according to the linear relation of the turbidity data of the detected water quality instrument, and constructing a ground turbidity data map according to the detection motion of the autonomous underwater robot.
The method comprises the following steps of carrying out digital modeling according to turbidity data of a water quality instrument, predicting the turbidity of an undetected area according to the linear relation of the turbidity data of the detected water quality instrument, and constructing a ground turbidity data map according to the detection motion of the autonomous underwater robot:
1) constructing a linear relation between adjacent turbidity values by using turbidity data of the water quality instrument and the time dependence relation of forward motion of the autonomous underwater robot;
2) constructing a turbidity data model of the autonomous underwater robot in the forward motion direction and a turbidity data model of the autonomous underwater robot in the auxiliary motion direction, namely the right direction, by utilizing the linear relation between the adjacent turbidity values;
3) constructing a turbidity data prediction model according to a turbidity data model of the autonomous underwater robot in the forward motion direction and a turbidity data model of the autonomous underwater robot in the right motion direction;
4) repeating the steps 1) to 3) according to a plurality of detection motion lines of the autonomous underwater robot to obtain a plurality of turbidity data prediction models;
5) and fusing a plurality of turbidity data prediction models of the autonomous underwater robot to construct a geodetic turbidity data map based on autonomous underwater machine data perception.
The step 1) comprises the following steps:
establishing a motion coordinate system by taking the autonomous underwater robot as an origin, the forward direction of the autonomous underwater robot as the positive direction of an X axis, and the right direction of the autonomous underwater robot as the positive direction of a Y axis, wherein a turbidity z is a function of a two-dimensional space position (X, Y), namely z (X, Y) ═ f (X, Y), f represents a mapping from the two-dimensional space position (X, Y) to the turbidity z, and on the axis of the X direction, Y equals 0, z (X,0) is abbreviated as z (X), namely f (X,0) is abbreviated as f (X), and A is defined to represent a coordinate set of a measuring line of the autonomous underwater robot; current intervalOn the line-of-sight trajectory A of the autonomous underwater robot, i.e.When, definePoints representing the interval AAndthe function values are all known measurement values; defining pointsHas a gradient value ofThey are the coefficients of undetermined, m1,m2Are all constants.
1.1) constructing gradient model of key discrete points
In thatRespectively weighing delta n coordinate points (delta n usually takes a value between 20 and 50 in engineering, a user can set the value by himself, the suggested value is 30) on the left and right sides of the vicinity, and calculating the pointsGradient value ofThe method comprises the following steps:
where T is the matrix transpose operator,representing coordinate pointsThe increase in turbidity between its adjacent coordinate points,is a coordinate pointThe position increment between the adjacent coordinate points is an intermediate variable, and the calculation method is as follows:
whereinIs a coordinate pointThe position increment between its neighboring coordinate points,representing coordinate pointsTurbidity between adjacent coordinate pointsDegree increments, which are intermediate variables, are calculated as follows:
1.2) utilizing the Gaussian relation of the spatial data to serialize the discrete data and construct the linear relation between adjacent turbidity values:
the turbidity data model of the autonomous underwater robot in the forward motion direction and the turbidity data model of the autonomous underwater robot in the auxiliary motion direction, namely the right direction, are as follows:
where Δ Y is a position increment in the Y direction of the two-dimensional position coordinate (x,0), and is a set value.
The turbidity data prediction model is as follows:
where Δ Y is a position increment in the Y direction of the two-dimensional position coordinate (x,0), and is a set value.
Taking the linear interval in the X direction satisfying the z (X, y) mathematical modelWherein Defining a coordinate set A representing a measuring line of the autonomous underwater robot; definition B denotes autonomous underwater vehicleThe coordinate of a measuring line of the robot and the coordinate set of the unknown area in the prediction range are integrated, delta L represents the prediction range and is a set value, and L represents the maximum value of the distance of the unknown area. Calculating any point in the sexual interval according to the step 1) to obtain the turbidity of the measuring pointAnda value of, andandvalue of (A) and haze valueThereby predicting the seawater turbidity z of the adjacent unknown area outside the autonomous underwater robot survey line.
Step 5) comprises the following steps:
xN=xcosθ-ysinθ
yE=xsinθ+ycosθ
Wherein xNIs true north coordinate of the earth, yEIs the geodetic east coordinate. Theta is an autonomous underwater robot track angle;
establishing a lattice (x, y, z) of the motion coordinate system of the autonomous underwater robot by a polynomial z (x, y), and then converting the (x, y, z) into a geodetic coordinate system (x) by a solved projection formulaN,yEZ) to obtain a turbidity value of the actual seafloor coordinate.
The invention has the following beneficial effects and advantages:
1. the invention fuses a large amount of discrete turbidity data into a digital model, so that the turbidity data meets the linear relation, and more obvious evidence is provided for the research of seawater turbidity.
2. The turbidity data of the adjacent undetected area is predicted according to the linear relation of the known measurement data, so that submarine hydrothermal solution searching channels are widened, and effective data support is provided for the submarine turbidity detection in the later period.
3. The application range is wide. The turbidity processing method can be applied to the sailing type ocean detection equipment which can carry a water quality meter or a turbidity meter, such as an underwater glider, a cabled underwater robot and the like.
Drawings
FIG. 1 is a schematic view of the turbidity data processing flow of the water quality meter of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and examples.
The method comprises the steps of constructing an ocean turbidity data model of the autonomous underwater robot in the forward motion direction (the forward motion is the main motion direction of the autonomous underwater robot) by utilizing the time dependence relationship between the ocean turbidity data and the forward motion of the autonomous underwater robot; constructing a turbidity data model of the autonomous underwater robot in the right motion direction (the right motion is the auxiliary motion direction of the autonomous underwater robot) by utilizing the spatial dependence of the ocean turbidity data and the gradient information of the forward motion; model prediction is carried out on the turbidity of an area, which is close to an undetected area, of the autonomous underwater robot track on the basis of a turbidity data model of the autonomous underwater robot in the motion direction; and finally, fusing the data of a plurality of sampling strips of the autonomous underwater robot, and constructing a geodetic turbidity data map based on autonomous underwater machine data sensing.
According to the invention, an ocean turbidity data model which is established by taking the motion direction of the autonomous underwater robot as the leading factor by utilizing the time and space dependence relationship of turbidity data is used, so that the continuous time observation of the autonomous underwater robot platform is realized, and the method is more efficient than the traditional site discrete sampling observation in a space interpolation mode. The invention can digitally analyze the linear relation of the known measured turbidity data, predict the turbidity data of the adjacent undetected area according to the linear relation and solve the problem that the turbidity characteristics of the local area cannot be effectively analyzed by analyzing the discrete turbidity data in a single way. The digitalized turbidity data model greatly improves the efficiency of processing a large amount of discrete turbidity data, and the spatial predictability of the model can bring further scientific research value. The method has stronger practicability and can be expanded to the fields of underwater gliders and cabled underwater robots.
The invention aims to efficiently process turbidity data in a water quality instrument by using a water quality instrument data processing method carried on the basis of an autonomous underwater robot.
The technical scheme adopted by the invention for realizing the purpose is as follows: a turbidity data processing method of a water quality instrument based on an autonomous underwater robot comprises the following operation steps:
(1) the autonomous underwater robot carries effective water quality instrument equipment to enter a comb-shaped detection planning route;
(2) the autonomous underwater robot navigates a long enough distance on the comb-shaped detection measuring line and measures the turbidity data of the water quality instrument at the distance;
(3) uploading turbidity data of the water quality instrument or acquiring the turbidity data on line after the latent time is finished;
(4) and selecting turbidity data in the turbidity data of the water quality instrument to carry out digital modeling treatment, predicting the turbidity of a nearby undetected area according to a linear relation obtained by the detected turbidity data, and constructing a turbidity map under the geodetic system according to the detection motion of the autonomous underwater robot.
As shown in fig. 1, the operation steps of the autonomous underwater robot (4) are as follows:
actual turbidity data is a large amount of discrete data, and now a spatial temporal relationship is constructed using a continuous mathematical model description. Based on the space motion relation of the autonomous underwater robot, the motion of the autonomous underwater robot is decoupled into two motion directions which do not interfere with each other.
And assuming that x is the forward motion coordinate of the autonomous underwater robot and y is the right motion coordinate of the autonomous underwater robot. And z (x, y) is the turbidity value of the (x, y) point in the coordinate system of the autonomous underwater robot, wherein x is a forward motion coordinate, and y is a right motion coordinate.
Step 1, constructing an ocean turbidity data model of the autonomous underwater robot in the forward motion direction (the forward motion is the main motion direction of the autonomous underwater robot) by utilizing the time dependence relationship between the ocean turbidity data and the forward motion of the autonomous underwater robot.
Turbidity z is a function of two-dimensional spatial position (x, y), i.e. z (x, y) ═ f (x, y), f denotes a mapping from two-dimensional spatial position (x, y) to turbidity z. In the x-direction axis, y is 0, and z (x,0) is abbreviated as z (x), i.e., f (x,0) is abbreviated as f (x). Defining a coordinate set A representing a measuring line of the autonomous underwater robot; current intervalOn the line-of-sight trajectory A of the autonomous underwater robot, i.e.When, define Points representing the interval AAndthe function values are all known measurement values; defining points Has a gradient value ofThey are the coefficients of undetermined, m1,m2Are all constants.
1) And constructing a gradient model of the key discrete points.
In thatRespectively weighing delta n coordinate points (delta n usually takes a value between 20 and 50 in engineering, a user can set the value by himself, the suggested value is 30) on the left and right sides of the vicinity, and calculating the pointsGradient value ofThe method comprises the following steps:
where T is the matrix transpose operator,representing coordinate pointsThe increase in turbidity between its adjacent coordinate points,is a coordinate pointThe position increment between the adjacent coordinate points is an intermediate variable, and the calculation method is as follows:
whereinIs a coordinate pointThe position increment between its neighboring coordinate points,representing coordinate pointsThe turbidity increase between adjacent coordinate points, which are intermediate variables, is calculated as follows:
2) discrete data are serialized by using the Gaussian relation of spatial data.
Adjacent turbidity values in water quality meter equipment generally satisfy a linear relationship, satisfying a polynomial as follows:
and 2, constructing a turbidity data model of the autonomous underwater robot in the right motion direction (the right motion is the auxiliary motion direction of the autonomous underwater robot) by using the spatial dependence of the ocean turbidity data and the gradient information of the forward motion.
Because the right motion of the autonomous underwater robot is only the auxiliary motion direction, the data volume in the Y direction is insufficient, and the gradient in the plane is approximated based on the gradient in the X direction. Where Δ y is any distance that satisfies the y direction.
z(x,y)=z(x,0)+G(x,y)Δy········(2)
G (x, y) is a polynomial for gradient of z (x) in the x direction.
G (x, y) ═ G (x,0) ═ G (x) · z' (x) · z (3)
Step 1 and step 2, decoupling in X and Y directions respectively, and constructing a planar continuous model as follows:
where Δ Y is a position increment in the Y direction of the two-dimensional position coordinate (x,0), and is a set value.
And 3, performing model prediction on the turbidity of the area, adjacent to the undetected area, of the autonomous underwater robot track on the basis of a turbidity data model of the autonomous underwater robot in the motion direction.
The planar linear model of the motion measuring line is known as equation (5), the measured discrete turbidity data has a spatial linear relation, a smooth data model is constructed according to the internal data dependency relation, and the turbidity data of the adjacent position area after the end of each measuring line is predicted by the smooth data model. The prediction model is as follows:
where Δ Y is a position increment in the Y direction of the two-dimensional position coordinate (x,0), and is a set value.
Taking the linear interval in the X direction satisfying the z (X, y) mathematical modelWherein (A represents a coordinate set of a measuring line of the autonomous underwater robot; B represents a coordinate set of a measuring line of the autonomous underwater robot and a coordinate set of an adjacent unknown region; Delta L represents a prediction range which is a set value; L represents a maximum value of a distance from the unknown region), and turbidity of a measuring point can be obtained by calculating the formula (1)Andthe value of (c). Obtaining from the formula (6)Andvalue of (A) and haze valueTherefore, the prediction of the seawater turbidity z of an unknown area outside the measuring line area of the autonomous underwater robot is realized.
And 4, fusing the data of a plurality of sampling strips of the autonomous underwater robot, and constructing a geodetic turbidity data map based on autonomous underwater machine data sensing.
Since the autonomous underwater robot generally adopts a geodetic coordinate system in practical scientific research application, a motion coordinate system (x is a forward motion coordinate and y is a right motion coordinate) of the autonomous underwater robot needs to be projected into the geodetic coordinate system.
Wherein xNIs true north coordinate of the earth, yEIs the geodetic east coordinate. ThetaIs the autonomous underwater robot track angle. Any survey line direction of the autonomous underwater robot has a determined track angle relative to the geodetic coordinates. Establishing a lattice (x, y, z) of an autonomous underwater robot motion coordinate system (x is a forward motion coordinate, and y is a right motion coordinate) by a polynomial z (x, y), and then converting the lattice (x, y, z) into a geodetic coordinate system (x, y, z) by a formula (7)N,yEZ) the turbidity value of the actual seafloor coordinate is determined according to the method.
The algorithm for projecting the motion coordinate system of the autonomous underwater robot to the geodetic coordinate system is as follows:
TABLE 1 Algorithm for projecting motion coordinate system of autonomous underwater robot to geodetic coordinate system
Wherein i represents the ith detection measuring line executed by the autonomous underwater robot, namely the ith turbidity strip under the motion coordinate system of the autonomous underwater robot; and N represents the total detection line number of the autonomous underwater robot, namely the total turbidity strip number under the motion coordinate system of the autonomous underwater robot.
Claims (6)
1. A turbidity data processing method based on an autonomous underwater robot platform is characterized by comprising the following steps:
carrying out detection movement on a planned line by an autonomous underwater robot carrying water quality instrument equipment, and measuring and recording turbidity data of the water quality instrument on the line;
and carrying out digital modeling according to the turbidity data of the water quality instrument, predicting the turbidity of an undetected area according to the linear relation of the turbidity data of the detected water quality instrument, and constructing a ground turbidity data map according to the detection motion of the autonomous underwater robot.
2. The turbidity data processing method based on the autonomous underwater robot platform as claimed in claim 1, wherein the steps of performing digital modeling according to the turbidity data of the water quality instrument, predicting the turbidity of an undetected area according to the linear relation of the turbidity data of the detected water quality instrument, and constructing a map of the turbidity data of the ground system according to the detected motion of the autonomous underwater robot comprise:
1) constructing a linear relation between adjacent turbidity values by using turbidity data of the water quality instrument and the time dependence relation of forward motion of the autonomous underwater robot;
2) constructing a turbidity data model of the autonomous underwater robot in the forward motion direction and a turbidity data model of the autonomous underwater robot in the auxiliary motion direction, namely the right direction, by utilizing the linear relation between the adjacent turbidity values;
3) constructing a turbidity data prediction model according to a turbidity data model of the autonomous underwater robot in the forward motion direction and a turbidity data model of the autonomous underwater robot in the right motion direction;
4) repeating the steps 1) to 3) according to a plurality of detection motion lines of the autonomous underwater robot to obtain a plurality of turbidity data prediction models;
5) and fusing a plurality of turbidity data prediction models of the autonomous underwater robot to construct a geodetic turbidity data map based on autonomous underwater machine data perception.
3. The method for processing the turbidity data based on the autonomous underwater robot platform as claimed in claim 2, wherein the step 1) comprises:
establishing a motion coordinate system by taking the autonomous underwater robot as an origin, the forward direction of the autonomous underwater robot as the positive direction of an X axis, and the right direction of the autonomous underwater robot as the positive direction of a Y axis, wherein a turbidity z is a function of a two-dimensional space position (X, Y), namely z (X, Y) ═ f (X, Y), f represents a mapping from the two-dimensional space position (X, Y) to the turbidity z, and on the axis of the X direction, Y equals 0, z (X,0) is abbreviated as z (X), namely f (X,0) is abbreviated as f (X), and A is defined to represent a coordinate set of a measuring line of the autonomous underwater robot; current intervalIn the autonomyOn the line-measuring trajectory A of the underwater robot, i.e.When, definePoints representing the interval AAndthe function values are all known measurement values; defining pointsHas a gradient value ofThey are the coefficients of undetermined, m1,m2Are all constants:
1.1) constructing gradient model of key discrete points
In thatRespectively weighing delta n coordinate points around, and calculating pointsGradient value ofThe method comprises the following steps:
where T is the matrix transpose operator,representing coordinate pointsThe increase in turbidity between its adjacent coordinate points,is a coordinate pointThe position increment between the adjacent coordinate points is an intermediate variable, and the calculation method is as follows:
whereinIs a coordinate pointThe position increment between its neighboring coordinate points,representing coordinate pointsThe turbidity increase between adjacent coordinate points, which are intermediate variables, is calculated as follows:
1.2) utilizing the Gaussian relation of the spatial data to serialize the discrete data and construct the linear relation between adjacent turbidity values:
4. the method for processing the turbidity data based on the autonomous underwater robot platform as claimed in claim 2, wherein the turbidity data model of the autonomous underwater robot in the forward motion direction and the turbidity data model of the autonomous underwater robot in the auxiliary motion direction, i.e. the right direction, are as follows:
where Δ Y is a position increment in the Y direction of the two-dimensional position coordinate (x,0), and is a set value.
5. The method for processing the turbidity data based on the autonomous underwater robot platform as claimed in claim 2, wherein the turbidity data prediction model is:
wherein, Δ Y is a position increment of the two-dimensional position coordinate (x,0) in the Y direction and is a set value;
taking the number satisfying z (x, y)Linear interval of learning model in X directionWherein Defining a coordinate set A representing a measuring line of the autonomous underwater robot; defining B to represent the coordinate of a measuring line of the autonomous underwater robot and an unknown area coordinate set in a prediction range, wherein delta L represents the prediction range and is a set value, and L represents the maximum value of the distance of the unknown area. Calculating any point in the sexual interval according to the step 1) to obtain the turbidity of the measuring pointAnda value of, andandvalue of (A) and haze valueThereby predicting the seawater turbidity z of the adjacent unknown area outside the autonomous underwater robot survey line.
6. The method for processing the turbidity data based on the autonomous underwater robot platform as claimed in claim 2, wherein the step 5) comprises:
xN=xcosθ-ysinθ
yE=xsinθ+ycosθ
Wherein xNIs true north coordinate of the earth, yEIs the geodetic east coordinate. Theta is an autonomous underwater robot track angle;
establishing a lattice (x, y, z) of the motion coordinate system of the autonomous underwater robot by a polynomial z (x, y), and then converting the (x, y, z) into a geodetic coordinate system (x) by a solved projection formulaN,yEZ) to obtain a turbidity value of the actual seafloor coordinate.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010296975.9A CN113535672B (en) | 2020-04-16 | 2020-04-16 | Turbidity data processing method based on autonomous underwater robot platform |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010296975.9A CN113535672B (en) | 2020-04-16 | 2020-04-16 | Turbidity data processing method based on autonomous underwater robot platform |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113535672A true CN113535672A (en) | 2021-10-22 |
CN113535672B CN113535672B (en) | 2024-03-12 |
Family
ID=78088299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010296975.9A Active CN113535672B (en) | 2020-04-16 | 2020-04-16 | Turbidity data processing method based on autonomous underwater robot platform |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113535672B (en) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN203053954U (en) * | 2012-08-13 | 2013-07-10 | 河海大学 | Water-quality multi-parameter underwater three-dimensional distribution detection, marking and data analysis system |
KR101339401B1 (en) * | 2012-12-14 | 2013-12-09 | 경운대학교 산학협력단 | Underwater probing robot using usn |
CN103796912A (en) * | 2011-05-17 | 2014-05-14 | 艾尼股份公司 | Autonomous underwater system for 4D environmental monitoring |
US20160335291A1 (en) * | 2015-05-13 | 2016-11-17 | U.S.A. Represented By The Administrator Of The National Aeronautics And Space Administration | System and method for providing a modern-era retrospective analysis for research and applications (merra) data analytic service |
RU2609618C1 (en) * | 2015-11-05 | 2017-02-02 | Владимир Васильевич Чернявец | Underwater robot system |
KR20170075928A (en) * | 2015-12-23 | 2017-07-04 | 포항공과대학교 산학협력단 | Method and apparatus for controlling position of underwater robot and underwater robot having the same |
US20180004217A1 (en) * | 2014-12-16 | 2018-01-04 | Robert Bosch Gmbh | Method for Mapping a Processing Area for Autonomous Robot Vehicles |
WO2018071416A1 (en) * | 2016-10-11 | 2018-04-19 | Kaarta, Inc. | Laser scanner with real-time, online ego-motion estimation |
CN108008099A (en) * | 2017-11-01 | 2018-05-08 | 佛山科学技术学院 | A kind of pollution sources localization method |
BR102017012861A2 (en) * | 2017-06-14 | 2019-01-15 | Universidade De São Paulo - Usp | autonomous oceanographic variables acquisition system |
CN109343350A (en) * | 2018-11-20 | 2019-02-15 | 清华大学 | A kind of underwater robot path tracking control method based on Model Predictive Control |
CN109901205A (en) * | 2019-03-14 | 2019-06-18 | 南京信息工程大学 | A kind of underwater robot Multi-sensor Fusion and motion profile prediction technique |
CN110045079A (en) * | 2019-05-13 | 2019-07-23 | 自然资源部第二海洋研究所 | A kind of Hydrothermal plumes autonomous classification sampling apparatus and method for autonomous underwater robot |
CN110794855A (en) * | 2019-11-26 | 2020-02-14 | 南方电网调峰调频发电有限公司 | Comprehensive control system and method for underwater robot |
-
2020
- 2020-04-16 CN CN202010296975.9A patent/CN113535672B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103796912A (en) * | 2011-05-17 | 2014-05-14 | 艾尼股份公司 | Autonomous underwater system for 4D environmental monitoring |
CN203053954U (en) * | 2012-08-13 | 2013-07-10 | 河海大学 | Water-quality multi-parameter underwater three-dimensional distribution detection, marking and data analysis system |
KR101339401B1 (en) * | 2012-12-14 | 2013-12-09 | 경운대학교 산학협력단 | Underwater probing robot using usn |
US20180004217A1 (en) * | 2014-12-16 | 2018-01-04 | Robert Bosch Gmbh | Method for Mapping a Processing Area for Autonomous Robot Vehicles |
US20160335291A1 (en) * | 2015-05-13 | 2016-11-17 | U.S.A. Represented By The Administrator Of The National Aeronautics And Space Administration | System and method for providing a modern-era retrospective analysis for research and applications (merra) data analytic service |
RU2609618C1 (en) * | 2015-11-05 | 2017-02-02 | Владимир Васильевич Чернявец | Underwater robot system |
KR20170075928A (en) * | 2015-12-23 | 2017-07-04 | 포항공과대학교 산학협력단 | Method and apparatus for controlling position of underwater robot and underwater robot having the same |
WO2018071416A1 (en) * | 2016-10-11 | 2018-04-19 | Kaarta, Inc. | Laser scanner with real-time, online ego-motion estimation |
BR102017012861A2 (en) * | 2017-06-14 | 2019-01-15 | Universidade De São Paulo - Usp | autonomous oceanographic variables acquisition system |
CN108008099A (en) * | 2017-11-01 | 2018-05-08 | 佛山科学技术学院 | A kind of pollution sources localization method |
CN109343350A (en) * | 2018-11-20 | 2019-02-15 | 清华大学 | A kind of underwater robot path tracking control method based on Model Predictive Control |
CN109901205A (en) * | 2019-03-14 | 2019-06-18 | 南京信息工程大学 | A kind of underwater robot Multi-sensor Fusion and motion profile prediction technique |
CN110045079A (en) * | 2019-05-13 | 2019-07-23 | 自然资源部第二海洋研究所 | A kind of Hydrothermal plumes autonomous classification sampling apparatus and method for autonomous underwater robot |
CN110794855A (en) * | 2019-11-26 | 2020-02-14 | 南方电网调峰调频发电有限公司 | Comprehensive control system and method for underwater robot |
Non-Patent Citations (2)
Title |
---|
张少伟;俞建成;张艾群;赵文涛;: "多水下机器人自主海洋特征场跟踪研究", 科学通报, no. 2 * |
潘爽;施建礼;聂永芳;王希彬;: "自主水下航行器同时定位与制图技术研究", 舰船科学技术, no. 11 * |
Also Published As
Publication number | Publication date |
---|---|
CN113535672B (en) | 2024-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108362281B (en) | Long-baseline underwater submarine matching navigation method and system | |
CN105823480B (en) | Underwater moving target location algorithm based on single beacon | |
Kularatne et al. | Time and Energy Optimal Path Planning in General Flows. | |
Lermusiaux et al. | Optimal planning and sampling predictions for autonomous and Lagrangian platforms and sensors in the northern Arabian Sea | |
Michini et al. | Robotic tracking of coherent structures in flows | |
CN111695299B (en) | Mesoscale vortex trajectory prediction method | |
CN105352496A (en) | AUV navigation method based on sonar-assisted autonomous navigation | |
de Fommervault et al. | Seaexplorer underwater glider: A new tool to measure depth-resolved water currents profiles | |
Zhang et al. | Ocean current-aided localization and navigation for underwater gliders with information matching algorithm | |
Wu et al. | Glider CT: Reconstructing flow fields from predicted motion of underwater gliders | |
CN113535672B (en) | Turbidity data processing method based on autonomous underwater robot platform | |
CN112215395A (en) | Underwater equipment adaptability information guarantee system based on ocean big data | |
Suzuki et al. | SLAM using ICP and graph optimization considering physical properties of environment | |
Demin et al. | A numerical method of calculation of currents and sea surface topography in multiply connected domains of the ocean | |
Yang et al. | Energy-optimal path planning with active flow perception for autonomous underwater vehicles | |
Tuohy et al. | Map based navigation for autonomous underwater vehicles | |
Padial et al. | Correlation of imaging sonar acoustic shadows and bathymetry for ROV terrain-relative localization | |
CN114417601B (en) | Method for quickly estimating mesoscale eddy underwater sound velocity field based on satellite altimeter data | |
CN113190025B (en) | Motion control method suitable for variable-structure underwater robot | |
Yiğit | Path planning methods for autonomous underwater vehicles | |
CN114114438A (en) | Quasi-three-dimensional inversion method for loop source ground-air transient electromagnetic data | |
Lin et al. | Underwater geomagnetic navigation based on ICP algorithm | |
Shi et al. | Cooperative estimation of a three-dimensional flow field using multiple AUVs and local measurements | |
Mugnai et al. | Developing affordable bathymetric analysis techniques using non-conventional payload for cultural heritage inspections | |
CN115047406B (en) | Reconstruction method of ground-air link propagation attenuation region |
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 |