CN109059796B - Shallow sea water depth multispectral satellite remote sensing inversion method for water depth control point-free area - Google Patents
Shallow sea water depth multispectral satellite remote sensing inversion method for water depth control point-free area Download PDFInfo
- Publication number
- CN109059796B CN109059796B CN201810805032.7A CN201810805032A CN109059796B CN 109059796 B CN109059796 B CN 109059796B CN 201810805032 A CN201810805032 A CN 201810805032A CN 109059796 B CN109059796 B CN 109059796B
- Authority
- CN
- China
- Prior art keywords
- water depth
- remote sensing
- shallow sea
- water
- depth
- 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
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/22—Measuring arrangements characterised by the use of optical techniques for measuring depth
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
A shallow sea water depth multispectral satellite remote sensing inversion method in a no-water-depth control point area belongs to the technical field of satellite ocean remote sensing application. Based on a derived two-waveband linear shallow sea water depth inversion model, an optimal waveband rotation unit vector is solved through pixel point pairs of data sets of different seabed types at different depths, seabed parameters in the model are estimated based on typical substrate type pixel adoption of land and water boundary line positions, meanwhile, a blue-green waveband diffusion attenuation coefficient ratio is obtained through pixel data of the same substrate type at different depths, and a green waveband attenuation coefficient is calculated through a half-analysis and diffusion attenuation coefficient algorithm based on deep water area data of the nearest shallow sea area. And (4) calculating the parameters to realize shallow sea water depth remote sensing inversion of the water-depth-free control point area.
Description
Technical Field
The invention belongs to the technical field of satellite marine remote sensing application, and particularly relates to a shallow sea water depth multispectral satellite remote sensing inversion method in a water depth control point-free area.
Background
The depth of shallow sea water is used as an important parameter for safe navigation guarantee of ships and warships, offshore ecosystems and optical research, and is always an important content of ocean mapping and optical remote sensing. The conventional shallow sea water depth measurement is mainly ship-borne multi/single beam acoustic measurement, and the airborne laser depth measurement technology developed in recent years is gradually widely applied. However, for some dangerous or disputed sea areas, these techniques are time consuming, laborious and even impossible to implement. Although the detection depth and precision cannot replace the conventional ocean measurement, the satellite remote sensing technology is the only feasible method for acquiring shallow sea water depth data in the regions. Therefore, the development of the research on the shallow sea water depth satellite remote sensing inversion technology has important significance and application prospect.
The existing shallow sea water depth satellite remote sensing inversion method can be generally summarized into two categories: the high spectrum remote sensing based on the semi-analytical algorithm and the multi spectrum remote sensing based on the empirical model are basically equivalent in the current water depth inversion accuracy from the literature reports. Although the hyperspectral remote sensing has the advantages of definite physical basis, no need of actually measuring water depth points and the like, the existing hyperspectral image has the defects of low spatial resolution and less available data. In contrast, the spatial resolution of the multispectral image can reach 2m at most, and the number of available satellites is large, so that the multispectral image is more suitable for carrying out shallow sea water depth remote sensing. However, the multispectral shallow sea water depth inversion model requires a certain amount of actually measured or reliable sea chart water depth data as input to carry out model coefficient calculation. Due to the change of the water body property and the seabed type, the multispectral water depth inversion model has obvious regionality. For some remote, dangerous or controversial remote sensing inversion methods, the number of measured water depth points is small, the water depth inversion result is unreliable, and sometimes even no available water depth control point exists, so that the shallow sea water depth remote sensing inversion cannot be carried out. Therefore, it is necessary to develop a shallow sea water depth multispectral satellite remote sensing inversion method without a water depth control point area.
Disclosure of Invention
The invention aims to provide a shallow sea water depth multispectral satellite remote sensing inversion method in a region without a water depth control point by analyzing and expressing a two-waveband linear model and combining multispectral image sampling points to calculate parameters required by the shallow sea water depth model aiming at the difficult problem that the multispectral shallow sea water depth remote sensing inversion requires the water depth control point.
The invention comprises the following steps:
1) deriving a shallow sea water depth inversion formula (1) according to a vector product form of the dual-waveband linear model:
Xi＝ln[rw(λi)-rdp(λi)]i＝1,2
wherein z is the shallow sea depth to be inverted α1,α2Weighting feature vectors of blue-green wave bands; g1,g2The water body double-pass diffusion attenuation coefficient is a blue-green waveband spectrum; r isw(λi) Remote reflectance below the water surface in the ith band (cyan); r isb(λi) The seabed remote sensing reflectivity of the ith wave band (blue green); r isdp(λi) The remote sensing reflectivity under the water surface of the ith wave band (blue green) optical deepwater zone; obtaining a water depth remote sensing inversion result of the water depth point-free area;
2) selecting adjacent pixel pairs with different depths and different substrates on the image, and carrying out minimum solving on an Xi data set of the adjacent pixel pairs to obtain a group of optimal waveband rotation unit vectors [ α ]1,α2]：
Where i denotes a certain adjacent pixel pair, Δ sziA, B shows different substrate types corresponding to the pixel point pair, n is the number of the pixel point pair, and f is a minimization function;
3) selecting a plurality of typical seabed sediment image element sets at the water line of the image, and combining to obtain an optimal waveband rotation unit vector [ α ]1,α2]Obtaining sea bottom parameters by average statistics
A value;
4) using the same type of seabed substrate on the image and X on different depth positions1～X2Data set, calculating blue-green wave band double-pass diffusion attenuation coefficient ratio g1/g2；
5) On the premise of assuming uniform water body properties, calculating the green band diffusion attenuation coefficient g nearest to the optical deepwater zone of the shallow sea area by using a semi-analysis and diffusion attenuation coefficient algorithm2；
6) The coefficients obtained by the calculation of the steps comprise [ α ]1、α2]Subsea parameter
g1/g2And g2Substituting the shallow sea water depth inversion formula, and applying the shallow sea water depth inversion formula to the whole image to realize the shallow sea water depth multispectral satellite remote sensing inversion of the region without the water depth control point.
Compared with the prior art, the invention has the following advantages:
1) the method is based on a derived two-waveband linear shallow sea water depth inversion model, an optimal waveband rotation unit vector is solved through pixel point pairs of data sets of different seabed types at different depths, seabed parameters in the model are estimated based on typical substrate type pixel adoption of land and water boundary line positions, meanwhile, a blue-green waveband diffusion attenuation coefficient ratio is obtained through pixel data of the same substrate type at different depths, a deep water region data of the nearest shallow sea region is based, and a green waveband attenuation coefficient is calculated through a half-analysis and diffusion attenuation coefficient algorithm. And (4) calculating the parameters to realize shallow sea water depth remote sensing inversion of the water-depth-free control point area.
2) Compared with the existing water depth inversion method which needs the participation of water depth control points, the method can eliminate the influence of the type difference of the seabed sediment on the water depth inversion to a certain extent, and can realize the shallow sea water depth remote sensing inversion of the area without the water depth control points.
Drawings
FIG. 1 is a distribution diagram of image pixel sampling points, water depth precision verification points and the like in the embodiment of the invention;
FIG. 2 shows a blue and green band sandy seabed sampling point X in the embodiment of the invention1～X2Scatter plot and its linear fit;
FIG. 3 is a water depth remote sensing inversion result diagram in the embodiment of the invention;
FIG. 4 is a result of verifying the inversion water depth accuracy of the new method in the embodiment of the present invention.
Detailed Description
According to the method, aiming at the water depth multispectral satellite remote sensing of the water depth in the no-water-depth point area, relevant parameters required by water depth inversion are obtained through the image, and therefore shallow-sea water depth remote sensing inversion of the no-water-depth point area is achieved.
The following describes in detail the specific implementation process of the technical solution of the present invention with reference to the accompanying drawings and examples:
the method comprises the following steps: acquiring a high-resolution multispectral satellite image of a research area, carrying out atmospheric correction to obtain a reflectivity image, and converting reflectivity data into underwater remote sensing reflectivity data rw. When the satellite image positioning error is lower than 6m, geometric correction processing needs to be carried out firstly; when the image has obvious flare interference, flare correction processing is required to be carried out;
step two: through visual interpretation, selecting the pixels (figure 1) of the deep water area nearest to the shallow sea area, and counting the blue and green bands rdpValue and according to formula Xi＝ln[rw(λi)-rdp(λi)]Calculating X of blue and green wave bandiImage data;
selecting adjacent pixel point pairs of different substrate types from different distances (representing different depths) along the shore line to obtain a data set (figure 1), and calculating by using an optimal algorithm (formula 2) to obtain an optimal waveband rotation unit vector [ α ]1,α2]；
Step four: referring to a near infrared image, selecting a pixel set (figure 1) of a typical substrate type on an amphibious boundary line of the image, and calculating XiData, and the obtained optimal band rotation unit vector [ α ]1,α2]Calculating the sea floor parameters
Step five: x for selecting positions of sandy seabed and different depths on image1～X2Establishing a linear regression formula (shown as figure 2) of the data set (figure 1) and the data set by using a least square method to obtain a ratio g of the two-way diffusion attenuation coefficients of the water body in blue and green wave bands1/g2；
Step six: according to a semi-analysis and diffusion attenuation coefficient algorithm, calculating the sum g of upward and downward diffusion attenuation coefficients of green wave bands of pixels in a deep water area nearest to the shallow sea area2；
Step seven, calculating the obtained water-free deep point region [ α1、α2]、
g1/g2And g2Substituting the parameters into a formula 1, and applying the parameters to the whole image to obtain a water depth remote sensing inversion result of the water depth-free area (as shown in FIG. 3);
step eight: randomly selecting an actual measurement water depth data set (distribution is shown in figure 1) to carry out precision verification on the water depth remote sensing inversion result, drawing a scatter diagram for comparing the inversion water depth with the actual measurement water depth, and calculating a root mean square error RMSE (shown in figure 4). In this example, the RMSE error for the water depth inversion is 1.18 m.
Claims (1)
1. The shallow sea water depth multispectral satellite remote sensing inversion method of the no water depth control point area is characterized by comprising the following steps:
1) deriving a shallow sea water depth inversion formula according to a vector product form of the dual-waveband linear model:
Xi＝ln[rw(λi)-rdp(λi)]i＝1,2
wherein z is the shallow sea depth to be inverted α1,α2Weighting feature vectors of blue-green wave bands; g1,g2The water body double-pass diffusion attenuation coefficient is a blue-green waveband spectrum; r isw(λi) The remote sensing reflectivity under the water surface of the ith wave band; r isb(λi) The seabed remote sensing reflectivity is the ith wave band; r isdp(λi) The remote sensing reflectivity under the water surface of the i-band optical deep water area is obtained; obtaining a water depth remote sensing inversion result of the water depth point-free area;
2) selecting adjacent pixel pairs with different depths and different substrates on the image, and carrying out minimum solving on an Xi data set of the adjacent pixel pairs to obtain a group of optimal waveband rotation unit vectors [ α ]1,α2]：
Where i denotes a certain adjacent pixel pair, Δ sziA, B shows different substrate types corresponding to the pixel point pair, n is the number of the pixel point pair, and f is a minimization function;
3) selecting a plurality of typical seabed sediment image element sets at the water line of the image, and combining to obtain an optimal waveband rotation unit vector [ α ]1,α2]Obtaining sea bottom parameters by average statistics
A value;
4) using the same type of seabed substrate on the image and X on different depth positions1～X2Data set, calculating blue-green wave band double-pass diffusion attenuation coefficient ratio g1/g2；
5) On the premise of assuming uniform water body properties, calculating the green band diffusion attenuation coefficient g nearest to the optical deepwater zone of the shallow sea area by using a semi-analysis and diffusion attenuation coefficient algorithm2；
6) The coefficients obtained by the calculation of the steps comprise [ α ]1、α2]Subsea parameter
g1/g2And g2Substituting the shallow sea water depth inversion formula, and applying the shallow sea water depth inversion formula to the whole image to realize the shallow sea water depth multispectral satellite remote sensing inversion of the region without the water depth control point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810805032.7A CN109059796B (en) | 2018-07-20 | 2018-07-20 | Shallow sea water depth multispectral satellite remote sensing inversion method for water depth control point-free area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810805032.7A CN109059796B (en) | 2018-07-20 | 2018-07-20 | Shallow sea water depth multispectral satellite remote sensing inversion method for water depth control point-free area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109059796A CN109059796A (en) | 2018-12-21 |
| CN109059796B true CN109059796B (en) | 2020-07-31 |
Family
ID=64834858
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810805032.7A Active CN109059796B (en) | 2018-07-20 | 2018-07-20 | Shallow sea water depth multispectral satellite remote sensing inversion method for water depth control point-free area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109059796B (en) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110274858B (en) * | 2019-07-15 | 2021-08-31 | 南京吉泽信息科技有限公司 | Remote sensing method for estimating lake suspended sediment concentration by utilizing GOCI data |
| CN110823190B (en) * | 2019-09-30 | 2020-12-08 | 广州地理研究所 | Island reef shallow sea water depth prediction method based on random forest |
| CN111561916B (en) * | 2020-01-19 | 2021-09-28 | 自然资源部第二海洋研究所 | Shallow sea water depth uncontrolled extraction method based on four-waveband multispectral remote sensing image |
| CN111474122B (en) * | 2020-04-21 | 2022-12-20 | 自然资源部第二海洋研究所 | Remote sensing extraction method for shallow seabed material reflectivity |
| CN111651707B (en) * | 2020-05-28 | 2023-04-25 | 广西大学 | Tidal level inversion method based on optical shallow water region satellite remote sensing image |
| CN113140000A (en) * | 2021-03-26 | 2021-07-20 | 中国科学院东北地理与农业生态研究所 | Water body information estimation method based on satellite spectrum |
| CN113793374B (en) * | 2021-09-01 | 2023-12-22 | 自然资源部第二海洋研究所 | Method for inverting water depth based on water quality inversion result by improved four-band remote sensing image QAA algorithm |
| CN113960625B (en) * | 2021-10-22 | 2024-06-11 | 自然资源部第二海洋研究所 | Water depth inversion method based on satellite-borne single-photon laser active and passive remote sensing fusion |
| CN114459438B (en) * | 2022-01-10 | 2024-02-02 | 山东科技大学 | Method for judging validity of high-resolution multispectral water depth inversion data |
| CN114199827B (en) * | 2022-02-21 | 2022-05-10 | 中国石油大学(华东) | Remote sensing data-based method for inverting vertical change of PAR diffuse attenuation coefficient |
| CN114594503B (en) * | 2022-03-02 | 2024-06-25 | 中南大学 | Shallow sea topography inversion method, computer equipment and storage medium |
| CN115235431B (en) * | 2022-05-19 | 2024-05-14 | 南京大学 | Shallow sea water depth inversion method and system based on spectrum layering |
| CN114758254B (en) * | 2022-06-15 | 2022-11-18 | 中国地质大学（武汉） | Dual-band unsupervised water depth inversion method and system |
| CN115640670A (en) * | 2022-09-05 | 2023-01-24 | 浙江水利水电学院 | Terrain self-adaptive water depth model partition weighting fusion method |
| CN115422981B (en) * | 2022-11-04 | 2023-04-07 | 自然资源部第一海洋研究所 | Land and water classification method and system for single-frequency airborne laser sounding data and application |
| CN115797760B (en) * | 2023-01-29 | 2023-05-12 | 水利部交通运输部国家能源局南京水利科学研究院 | Active and passive fusion type water quality three-dimensional remote sensing inversion method, system and storage medium |
| CN117152636B (en) * | 2023-10-29 | 2024-03-15 | 自然资源部第二海洋研究所 | Shallow sea substrate reflectivity remote sensing monitoring method based on dual-band relation |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7042802B2 (en) * | 2003-09-18 | 2006-05-09 | Schlumberger Technology Corporation | Determination of stress characteristics of earth formations |
| CN104181515B (en) * | 2013-05-21 | 2016-03-09 | 时春雨 | A kind of shallow water depth inversion method based on blue-yellow wave band high-spectral data |
| CN105651263B (en) * | 2015-12-23 | 2018-02-23 | 国家海洋局第一海洋研究所 | Shallow water depth multi-source remote sensing merges inversion method |
| CN105627997A (en) * | 2015-12-23 | 2016-06-01 | 国家海洋局第一海洋研究所 | Multi-angle remote sensing water depth decision fusion inversion method |
| CN105865424B (en) * | 2016-04-13 | 2019-09-17 | 中测新图(北京)遥感技术有限责任公司 | A kind of multispectral remote sensing inversion method and device based on nonlinear model |
-
2018
- 2018-07-20 CN CN201810805032.7A patent/CN109059796B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109059796A (en) | 2018-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109059796B (en) | Shallow sea water depth multispectral satellite remote sensing inversion method for water depth control point-free area | |
| CN105445751B (en) | A kind of shallow water area depth of water ratio remote sensing inversion method | |
| CN110749568B (en) | MODIS remote sensing inversion method for phytoplankton absorption characteristics of high-turbidity eutrophic lakes | |
| CN105865424B (en) | A kind of multispectral remote sensing inversion method and device based on nonlinear model | |
| CN102176001B (en) | Permeable band ratio factor-based water depth inversion method | |
| CN112013822A (en) | Multispectral remote sensing water depth inversion method based on improved GWR model | |
| CN105651263B (en) | Shallow water depth multi-source remote sensing merges inversion method | |
| CN111781146B (en) | Wave parameter inversion method using high-resolution satellite optical image | |
| Kanno et al. | Shallow water bathymetry from multispectral satellite images: Extensions of Lyzenga's method for improving accuracy | |
| CN109781073B (en) | Shallow sea water depth remote sensing extraction method integrating sea wave characteristics and spectral characteristics | |
| CN104156629B (en) | A kind of pathfinder image inverting wind direction of ocean surface method based on relative detector calibration | |
| CN111651707A (en) | Tidal level inversion method based on optical shallow water satellite remote sensing image | |
| CN114201732A (en) | Sentinel-2A image-based shallow sea water depth inversion method | |
| CN104197902A (en) | Method for extracting shallow sea terrain by single-shot high-resolution optical remote sensing image | |
| CN115979972B (en) | Real-time monitoring method and system for hyperspectral of crude oil film on sea surface | |
| CN114297938A (en) | Optical shallow water bottom depth inversion method based on neural network | |
| CN111947628B (en) | Linear water depth inversion method based on inherent optical parameters | |
| CN113255144A (en) | Shallow sea remote sensing water depth inversion method based on FUI partition and Randac | |
| Zhu et al. | An APMLP deep learning model for bathymetry retrieval using adjacent pixels | |
| CN117274831B (en) | Offshore turbid water body depth inversion method based on machine learning and hyperspectral satellite remote sensing image | |
| CN114973010A (en) | Remote sensing water depth inversion model migration method based on no prior data | |
| Ashphaq et al. | Evaluation and performance of satellite-derived bathymetry algorithms in turbid coastal water: a case study of Vengurla rocks | |
| CN110097562B (en) | Sea surface oil spill area image detection method | |
| CN115235431B (en) | Shallow sea water depth inversion method and system based on spectrum layering | |
| Liu et al. | Bathymetric ability of SPOT-5 multi-spectral image in shallow coastal water |
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 | ||
| CB02 | Change of applicant information |
Address after: 361005 Fujian Province University of Siming District of Xiamen City Road No. 178 Applicant after: THIRD INSTITUTE OF OCEANOGRAPHY, MINISTRY OF NATURAL RESOURCES Address before: 361005 Fujian Province University of Siming District of Xiamen City Road No. 178 Applicant before: Third Institute of Oceanography, State Oceanic Administration |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |