CN110849311A - Estimation method for sea ice output area flux of polar region key channel - Google Patents
Estimation method for sea ice output area flux of polar region key channel Download PDFInfo
- Publication number
- CN110849311A CN110849311A CN201911132406.4A CN201911132406A CN110849311A CN 110849311 A CN110849311 A CN 110849311A CN 201911132406 A CN201911132406 A CN 201911132406A CN 110849311 A CN110849311 A CN 110849311A
- Authority
- CN
- China
- Prior art keywords
- sea ice
- channel
- grid
- vector
- drift velocity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000004907 flux Effects 0.000 title claims abstract description 34
- 101000703681 Homo sapiens Single-minded homolog 1 Proteins 0.000 claims description 6
- 102100031980 Single-minded homolog 1 Human genes 0.000 claims description 6
- 241000764238 Isis Species 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 238000009827 uniform distribution Methods 0.000 claims 1
- 230000005012 migration Effects 0.000 abstract description 4
- 238000013508 migration Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 101000616761 Homo sapiens Single-minded homolog 2 Proteins 0.000 description 1
- 102100021825 Single-minded homolog 2 Human genes 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008278 dynamic mechanism Effects 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/28—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring areas
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G06T3/06—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20068—Projection on vertical or horizontal image axis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30192—Weather; Meteorology
-
- 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/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to a method for estimating sea ice output area flux of a polar region key channel. Obtaining Sea Ice migration Speed (SIM) data and Sea Ice Concentration (SIC) data by using satellite-based remote sensing data, determining grids covered by channels according to channel positions and data resolution, calculating SIM speed components perpendicular to the channels at each grid, and accumulating products of SIM speeds, time, grid side lengths and Sea Ice Concentration perpendicular to the channels of all the grids to obtain Sea Ice area output flux of a specific channel. The method solves the problems of low accuracy of basic data, lack of long time sequence results and the like of estimation of sea ice area output flux of the key channel in the polar region, and is beneficial to understanding of the material balance of the sea ice in the two polar regions in a new climate environment. The method is high in precision, and the sea ice output area flux of any channel in the polar region can be calculated according to requirements.
Description
Technical Field
The invention relates to a method for estimating sea ice output area flux of a polar region key channel, in particular to a flow algorithm for estimating sea ice area output of the key channel.
Background
Under new climatic conditions, the bipolar sea ice plays an important role in indicating global climate change and is widely concerned in various aspects. The sea ice output through each channel is also an important dynamic mechanism influencing the change of the quality of the arctic sea ice, plays an important role in maintaining the mass balance of the north ocean, and in addition, the fresh water generated by the sea ice output can change the strength of the radial overturning flow of the north ocean sea area, thereby influencing the global ocean temperature and salt circulation. The early sea ice migration output observation mode is mainly based on limited field buoy and submerged buoy measurement data, and a convincing quantitative result is difficult to obtain due to the fact that field observation data samples are rare and time sequences are not coordinated. The data obtained by the satellite sensor has the advantages of high time resolution and historical data accumulation, and the characteristics lay a foundation for obtaining the output of the long-time sequence arctic sea ice area. Besides making up for the data deficiency of the existing estimation method in the aspect of basic data, the estimation method has clear logic, higher scientificity and accuracy and easy realization, and can acquire sea ice output area flux of any channel by modifying the channel position.
Disclosure of Invention
In view of the technical deficiencies, the invention aims to provide a method for estimating sea ice output area flux of a polar region key channel. According to the method, sea ice migration speed data and sea ice density data extracted based on satellite remote sensing data are used, a series of data grids covered by a channel are determined according to the position of the channel and the resolution of the data, sea ice area output flux at each grid is calculated, and the sea ice area output flux of a specific channel is obtained by accumulating all the grids. The method solves the problems of low precision of basic data, lack of long time sequence results and the like of sea ice area output flux estimation of critical channels (straits, gulfs and the like) in polar regions, is simple and convenient to operate and is easy to realize.
The technical scheme adopted by the invention for solving the technical problems is as follows: a method for estimating sea ice output area flux of a polar region key channel comprises the following steps:
1) acquiring position coordinates of two ends of the channel under a longitude and latitude coordinate system;
2) respectively establishing projection coordinate systems corresponding to the sea ice drift speed and the sea ice density data;
3) respectively converting the position coordinates of the two ends of the channel under the longitude and latitude coordinate systems into end point coordinates under the two projection coordinate systems; respectively dispersing the channel into a plurality of uniformly distributed discrete points under two projection coordinate systems, and calculating the projection coordinate of each discrete point under the projection coordinate system;
4) for the projection coordinates of each discrete point under the two projection coordinate systems, searching a grid which is closest to the surrounding grid, and taking the attribute data of the grid as the attribute value of the discrete point under the current projection coordinate system; the attribute data of the sea ice drift velocity projection coordinate system is a sea ice drift velocity vector, and the attribute data of the sea ice density data projection coordinate system is a sea ice density numerical value;
5) under the sea ice drift velocity projection coordinate system, calculating the position relation between the sea ice drift velocity vector of the current grid and the channel vector to judge whether the sea ice flows into or flows out of the channel, and calculating the sea ice drift velocity component of the sea ice at the grid, which is vertical to the channel;
6) under the sea ice density data projection coordinate system, reading a sea ice density numerical value of a current grid, and taking the product of the size of a sea ice drift velocity component perpendicular to a channel, the limiting time, the grid side length and the sea ice density numerical value as the sea ice area output flux in the limited time interval of the grid;
7) and repeating the steps 5) -6), carrying out iterative processing on all grids, and accumulating the sea ice area output flux in the limited time interval of each grid to obtain the sea ice area output flux in the limited time interval of the whole channel.
The projection coordinate systems respectively corresponding to the sea ice drift velocity and the sea ice density data are specifically as follows: and setting a projection type, an ellipsoid, a central meridian and a unit, and establishing a required coordinate system.
The two ends of the channel are in longitude and latitude coordinatesThe position coordinates under the system are respectively converted into the endpoint coordinates under the two projection coordinate systems, which are specifically as follows: according to the two established projection coordinate systems, the longitude and latitude coordinates (lambda,) Converting into projection coordinate (x) under corresponding projection coordinate system0,y0) And (x)1,y1) It is a projection conversion process.
The method is characterized in that the channel is dispersed into a plurality of uniformly distributed discrete points by adopting the interval resolution r, and the method specifically comprises the following steps:
suppose that the coordinates of the end points of the two ends of the channel under the current projection coordinate system are A (x)0,y0)、B(x1,y1);
Projection coordinate (x) of any discrete point j on the channelj,yj) Comprises the following steps:
xj=x0+r×cosα×j;
yj=y0+r×sinα×j。
the searching for the grid closest to the surrounding grid specifically includes: the distance between the discrete point and the center of the surrounding grid is calculated.
The calculation of the position relation between the sea ice drift velocity vector and the channel vector of the current grid is realized by calculating the cross product of the sea ice drift velocity vector and the channel vector:
after the corresponding data grid position is determined, the sea ice drift velocity SIM (u) at the current grid position is read1,v1) And calculating:wherein the channelThe method comprises the following steps of (1) knowing;
then: if it isIs negative to explainIn the clockwise direction of the channel, the sea ice at the current grid flows out of the channel, and vice versa represents the sea ice flowing into the channel.
The calculation of the sea ice drift velocity component of the sea ice at the grid, which is perpendicular to the channel, is specifically as follows:
and calculating the projection vector of the sea ice drift velocity on the channel vector according to a vector algorithm, and further calculating the size of the sea ice drift velocity component perpendicular to the channel at the grid.
The projection vector of the sea ice drift velocity on the channel vector calculated according to the vector algorithm is as follows:
calculate sea ice drift velocity vector SIM ((u) at the grid1,v1) Projection vector SIM1 (u) on channelm,vm):
Among them, SIM1 (u)m,vm) For the sea ice drift velocity vector SIM (u) at the grid1,v1) In the channelIs a projection vector of, theta is a vectorAndthe included angle of (A);
the calculation of the size of the sea ice drift velocity component perpendicular to the channel at the grid is as follows:
wherein u ispSIM2 (u) is the sea ice drift velocity component perpendicular to the channel at the current gridn,vn) The size of (2).
The method is implemented in IDL programming language.
The invention has the following beneficial effects and advantages:
1. the invention relates to a method for estimating sea ice output area flux of a polar region key channel, which can provide a whole set of flow algorithm from data reading, data processing and iterative operation.
2. The basic data adopted by the method is verified in a large quantity, so that the precision and the accuracy are higher.
3. The method can be completed by adopting various programming languages, and is simple, convenient and easy to operate.
4. The method has high flexibility, and users can obtain the estimation results of the sea ice area output flux at different channels only by changing the positions of the channels.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is an example channel position and pattern;
FIG. 3 is a geometric schematic of calculating the projected coordinates of discrete points covered by a channel;
FIG. 4 is a diagram of data grid locations corresponding to discrete points covered by the Davis channel;
FIG. 5 is a geometric schematic of the sea ice drift velocity component perpendicular to the channel at the computational grid.
Detailed Description
The present invention will be described in further detail with reference to examples. The method steps are explained with reference to the attached drawings. The channel in the following embodiments is exemplified by the davis channel in the arctic ocean, and is only for the purpose of illustrating the present invention, and is not intended to limit the scope of the present invention.
As shown in fig. 1, firstly, the channel position is determined, the projection coordinate systems corresponding to the sea ice drift velocity and the sea ice density data are respectively established, and the positions of the points covered by the channel in the SIC and SIM data coordinate systems are searched. Reading SIC and SIM data of corresponding grids, calculating a sea ice drift velocity component perpendicular to a channel at the grid, taking the product of the sea ice drift velocity component perpendicular to the channel, time, grid side length and sea ice concentration as a sea ice area output flux within a specific time interval of the grid, and accumulating the value of each grid to obtain the sea ice area output flux of the whole channel. The channel in the present invention refers to straits, bays and the like. The data used in the method is Sea Ice migration velocity (SIM) data and Sea Ice Concentration (SIC) data extracted based on satellite remote sensing data. The method mainly comprises the following steps:
as shown in fig. 2, the channel position is determined according to the longitude and latitude coordinates of the two ends of the channel, in the figure, the solid line represents the channel position, and the arrow represents the direction of the channel, which is used for the vector calculation in the subsequent step. And respectively establishing projection coordinate systems corresponding to the sea ice drift speed and the sea ice density data, wherein the coordinate system corresponding to the SIM data is an equal-area extensible grid, the coordinate system corresponding to the SIC data is polar orthographic projection (a coordinate system used for drawing in figure 1), and parameters such as a projection type, an ellipsoid, a central meridian, a unit and the like are set according to grid configuration to establish a required coordinate system.
As shown in fig. 3, the projected coordinates of the two end points of the channel in the sea ice drift velocity data coordinate system are calculated, which essentially are the geographic coordinates (λ,) And converting into projection coordinates (x, y) in a corresponding coordinate system. Let the coordinates of the starting end point of the channel be A (x)0,y0) The coordinate of the end point of the channel is B (x)1,y1) The angle α of the channel can then be calculated by the following equation:
the number n of occupied grids is calculated from the length of the channel and the resolution r (25km) of the data:
the formula II is as follows:
according to the coordinates (x) of the starting end point of the channel0,y0) Angle α of channel and number n of grids, discretizing a continuous straight line into points spaced at data resolution, and obtaining projection coordinates (x) of all discrete points covered by the channeli,yi) Where i is 1,2, … … n, for example, in fig. 3, point F is the jth point in the channel, and the coordinate of point F (x)j,yj) The calculation method comprises the following steps:
the formula III is as follows: x is the number ofj=x0+r×cosα×j
The formula four is as follows: y isj=y0+r×sinα×j
For each discrete point (x) covered by a channeli,yi) And determining the nearest grid in the sea ice drift data grid as a corresponding data grid, recording row and column values, and acquiring the position of the sea ice drift speed data of all grids covered by the channel. FIG. 4 represents data grid locations corresponding to discrete points covered by Davis strait channels, the channels corresponding to sea ice density data grids being searched in a manner similar to that described above.
After determining the corresponding data grid position, reading the sea ice drift velocity SIM (u) at one grid position1,v1) And judging the position relation between the sea ice motion velocity vector and the channel according to a formula five:
if it isIs negative to explainIn the clockwise direction of the channel, the ice flows out of the davis channel, whereas the ice flows into the channel as indicated by the light green arrows in fig. 5.
According to the relation between the channel and the sea ice drift velocity vector, firstly calculating the sea ice drift velocity vector SIM (u)1,v1) Projection vector SIM1 (u) on channelm,vm):
calculating to obtain sea ice drift velocity vector SIM (u)1,v1) Projection vector SIM1 (u) on channelm,vm) Thereafter, the sea ice drift velocity component SIM2 (u) perpendicular to the channel at this grid is calculated by the following formulan,vn) And the size u thereofp。
the formula is nine: u. ofn=u1-um,vn=v1-vm
Formula ten:
Then, reading the sea ice density data c at the grid position, and measuring the size u of the sea ice drift velocity component perpendicular to the channelpAnd the product of the grid side length r and the sea ice density c under the time t and SIC coordinate system is taken as the sea ice area output flux in the specific time interval of the grid, and the value of each grid is accumulated to obtain the sea ice area output flux of the whole channel (formula eleven).
wherein i is a grid serial number.
The method for estimating sea ice output area flux of the polar region key channel is realized by IDL programming language.
Claims (10)
1. A method for estimating sea ice output area flux of a polar region key channel is characterized by comprising the following steps:
1) acquiring position coordinates of two ends of the channel under a longitude and latitude coordinate system;
2) respectively establishing projection coordinate systems corresponding to the sea ice drift speed and the sea ice density data;
3) respectively converting the position coordinates of the two ends of the channel under the longitude and latitude coordinate systems into end point coordinates under the two projection coordinate systems; respectively dispersing the channel into a plurality of uniformly distributed discrete points under two projection coordinate systems, and calculating the projection coordinate of each discrete point under the projection coordinate system;
4) for the projection coordinates of each discrete point under the two projection coordinate systems, searching a grid which is closest to the surrounding grid, and taking the attribute data of the grid as the attribute value of the discrete point under the current projection coordinate system; the attribute data of the sea ice drift velocity projection coordinate system is a sea ice drift velocity vector, and the attribute data of the sea ice density data projection coordinate system is a sea ice density numerical value;
5) under the sea ice drift velocity projection coordinate system, calculating the position relation between the sea ice drift velocity vector of the current grid and the channel vector to judge whether the sea ice flows into or flows out of the channel, and calculating the sea ice drift velocity component of the sea ice at the grid, which is vertical to the channel;
6) under the sea ice density data projection coordinate system, reading a sea ice density numerical value of a current grid, and taking the product of the size of a sea ice drift velocity component perpendicular to a channel, the limiting time, the grid side length and the sea ice density numerical value as the sea ice area output flux in the limited time interval of the grid;
7) and repeating the steps 5) -6), carrying out iterative processing on all grids, and accumulating the sea ice area output flux in the limited time interval of each grid to obtain the sea ice area output flux in the limited time interval of the whole channel.
2. The method for estimating sea ice output area flux of a polar region key channel according to claim 1, wherein the projection coordinate systems respectively corresponding to the sea ice drift velocity and the sea ice density data are specifically: and setting a projection type, an ellipsoid, a central meridian and a unit, and establishing a required coordinate system.
3. The method for estimating sea ice output area flux of the polar region key channel according to claim 1, wherein the step of respectively converting the position coordinates of the two ends of the channel under the longitude and latitude coordinate systems into the end point coordinates under the two projection coordinate systems comprises the following steps: according to the two established projection coordinate systems, the longitude and latitude coordinates are converted into the coordinate systemConverting into projection coordinate (x) under corresponding projection coordinate system0,y0) And (x)1,y1) It is a projection conversion process.
4. The method for estimating sea ice output area flux of polar region key channel according to claim 3, wherein the discrete points for dispersing the channel into uniform distribution by using the interval resolution r are specifically: suppose that the coordinates of the end points of the two ends of the channel under the current projection coordinate system are A (x)0,y0)、B(x1,y1);
Projection coordinate (x) of any discrete point j on the channelj,yj) Comprises the following steps:
xj=x0+r×cosα×j;
yj=y0+r×sinα×j。
5. the method for estimating sea ice output area flux of polar region key channel according to claim 1, wherein the step of searching the grid closest to the surrounding grid is specifically as follows: the distance between the discrete point and the center of the surrounding grid is calculated.
6. The method for estimating sea ice output area flux of polar region key channel according to claim 1, wherein the calculating the position relation of sea ice drift velocity vector and channel vector of current grid is realized by calculating cross product of two:
after the corresponding data grid position is determined, the sea ice drift velocity SIM (u) at the current grid position is read1,v1) And calculating:wherein the channelThe method comprises the following steps of (1) knowing;
7. The method for estimating sea ice output area flux of polar region key channel according to claim 1, wherein said calculating the sea ice drift velocity component of the sea ice at the grid perpendicular to the channel is specifically:
and calculating the projection vector of the sea ice drift velocity on the channel vector according to a vector algorithm, and further calculating the size of the sea ice drift velocity component perpendicular to the channel at the grid.
8. The method for estimating sea ice output area flux of polar critical path according to claim 7, wherein the vector algorithm for calculating the projection vector of sea ice drift velocity on the path vector is:
calculate sea ice drift velocity vector SIM ((u) at the grid1,v1) Projection vector SIM1 (u) on channelm,vm):
Among them, SIM1 (u)m,vm) For the sea ice drift velocity vector SIM (u) at the grid1,v1) In the channelIs a projection vector of, theta is a vectorAndthe included angle of (A);
9. the method of claim 8, wherein the step of calculating the drift velocity component of the ice perpendicular to the channel at the grid comprises:
un=u1-um,vn=v1-vm
wherein u ispSIM2 (u) is the sea ice drift velocity component perpendicular to the channel at the current gridn,vn) The size of (2).
10. The method of claim 1, wherein the method is implemented in IDL programming language.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911132406.4A CN110849311B (en) | 2019-11-19 | 2019-11-19 | Estimation method for sea ice output area flux of polar region key channel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911132406.4A CN110849311B (en) | 2019-11-19 | 2019-11-19 | Estimation method for sea ice output area flux of polar region key channel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110849311A true CN110849311A (en) | 2020-02-28 |
CN110849311B CN110849311B (en) | 2021-03-26 |
Family
ID=69602584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911132406.4A Active CN110849311B (en) | 2019-11-19 | 2019-11-19 | Estimation method for sea ice output area flux of polar region key channel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110849311B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113689041A (en) * | 2021-08-25 | 2021-11-23 | 国家海洋环境预报中心 | Ice fixing condition prediction method and device, electronic equipment and storage medium |
CN116597319A (en) * | 2023-04-19 | 2023-08-15 | 北京师范大学 | Sea ice daily data comprehensive processing method and device and computing equipment |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107871323A (en) * | 2016-09-23 | 2018-04-03 | 国家海洋环境监测中心 | A kind of sea ice motion information acquisition method based on all solid state radar |
CN109254767A (en) * | 2018-08-28 | 2019-01-22 | 中国科学院海洋研究所 | A kind of polar region Sea Ice Drift vector visualized algorithm |
CN109447993A (en) * | 2018-10-25 | 2019-03-08 | 哈尔滨工程大学 | A kind of sea ice image partition method based on mixing true and false sample strategy |
CN110442665B (en) * | 2019-08-02 | 2022-09-30 | 上海海事大学 | Polar region sea area ice condition fusion display method |
-
2019
- 2019-11-19 CN CN201911132406.4A patent/CN110849311B/en active Active
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113689041A (en) * | 2021-08-25 | 2021-11-23 | 国家海洋环境预报中心 | Ice fixing condition prediction method and device, electronic equipment and storage medium |
CN113689041B (en) * | 2021-08-25 | 2022-04-29 | 国家海洋环境预报中心 | Ice fixing condition prediction method and device, electronic equipment and storage medium |
CN116597319A (en) * | 2023-04-19 | 2023-08-15 | 北京师范大学 | Sea ice daily data comprehensive processing method and device and computing equipment |
CN116597319B (en) * | 2023-04-19 | 2024-02-02 | 北京师范大学 | Sea ice daily data comprehensive processing method and device and computing equipment |
Also Published As
Publication number | Publication date |
---|---|
CN110849311B (en) | 2021-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108534779B (en) | Indoor positioning map construction method based on track correction and fingerprint improvement | |
CN110146895B (en) | Acoustic velocity profile inversion method based on inverted multi-beam echometer | |
CN106871880B (en) | The method that rib position, half-breadth and height carry out ship closure is calculated using three-dimensional coordinate | |
CN110849311B (en) | Estimation method for sea ice output area flux of polar region key channel | |
CN102004856B (en) | Rapid collective Kalman filtering assimilating method for real-time data of high-frequency observation data | |
CN109978275B (en) | Extreme strong wind speed prediction method and system based on mixed CFD and deep learning | |
JPWO2006104087A1 (en) | Depth distribution prediction method for predetermined water temperature range, migratory fish fishing ground prediction method, and migratory fish fishing ground prediction information distribution system | |
CN109543356A (en) | Consider the ocean interior temperature-salinity structure remote sensing inversion method of Space atmosphere | |
CN108896040B (en) | Inertia/gravity combined navigation method and system for sky-sea integrated underwater vehicle | |
Yang et al. | Marine surveying and mapping system based on Cloud Computing and Internet of Things | |
CN115307714B (en) | Rapid sound velocity distribution estimation method based on cross-space-time sound velocity profile clustering | |
CN106501774A (en) | A kind of underwater acoustic sensor network node positioning method | |
CN106570936B (en) | A kind of equidistant weight interpolation encryption method based on gridded DEM data | |
CN110441760B (en) | Wide-range seabed topographic map expansion composition method based on prior topographic map | |
CN109738902B (en) | High-precision autonomous acoustic navigation method for underwater high-speed target based on synchronous beacon mode | |
CN104507097A (en) | Semi-supervised training method based on WiFi (wireless fidelity) position fingerprints | |
US10122796B2 (en) | System for monitoring temperature and other conditions in water bodies | |
CN106488526B (en) | Mobile multi-hop underwater acoustic network dynamic method for self-locating based on layering | |
CN111123345A (en) | GNSS measurement-based empirical ionosphere model data driving method | |
CN109085656B (en) | Feature-oriented high-precision gravity graph construction and interpolation method | |
CN115423955B (en) | Multi-source data-based optimal depth reference surface geodetic height model construction method | |
CN113009531A (en) | Small-scale high-precision low-altitude troposphere atmosphere refractive index model | |
CN111579011A (en) | Walking type ADCP river flow testing method suitable for ultralow flow speed condition | |
CN113254850B (en) | River channel sludge amount monitoring and calculating method | |
CN108924734A (en) | A kind of three-dimension sensor node positioning method and system |
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 |