Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength
In this study, we present a new method called <i>BathySent</i> to retrieve shallow bathymetry from space that is based on the joint measurement of ocean wave celerity (c) and wavelength (λ). We developed the method to work with Sentinel 2 data, exploiting the time lag between two Sentine...
Main Authors: | , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
MDPI AG
2021-05-01
|
Series: | Remote Sensing |
Subjects: | |
Online Access: | https://www.mdpi.com/2072-4292/13/11/2149 |
_version_ | 1827691044525834240 |
---|---|
author | Marcello de Michele Daniel Raucoules Deborah Idier Farid Smai Michael Foumelis |
author_facet | Marcello de Michele Daniel Raucoules Deborah Idier Farid Smai Michael Foumelis |
author_sort | Marcello de Michele |
collection | DOAJ |
description | In this study, we present a new method called <i>BathySent</i> to retrieve shallow bathymetry from space that is based on the joint measurement of ocean wave celerity (c) and wavelength (λ). We developed the method to work with Sentinel 2 data, exploiting the time lag between two Sentinel 2 spectral bands, acquired quasi-simultaneously, from a single satellite dataset. Our method was based on the linear dispersion law, which related water depth to wave celerity and wavelength: when the water depth was less than about half the dominant wavelength, the wave celerity and wavelength decreased due to decreasing water depth (h) as the waves propagated towards the coast. Instead of using a best weighted <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> fit with the linear dispersion relation to retrieve h, we proposed solving the linear dispersion relation for each <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> pair to find multiple h-values within the same resolution cell. Then, we calculated the weighted averaged h-value for each resolution cell. To improve the precision of the final bathymetric map, we stacked the bathymetry values from <i>N</i>-different datasets acquired from the same study area on different dates. We first tested the algorithm on a set of images representing simulated ocean waves, then we applied it to a real set of Sentinel 2 data obtained of our study area, Gâvres peninsula (France, 47°,67 lat.; −3°35 lon.). A comparison with in situ bathymetry yielded good results from the synthetic images (r<sup>2</sup> = 0.9) and promising results with the Sentinel 2 images (r<sup>2</sup> = 0.7) in the 0–16 m depth zone. |
first_indexed | 2024-03-10T10:52:31Z |
format | Article |
id | doaj.art-4561a1e092e04fa4a51aebb1bb5c440f |
institution | Directory Open Access Journal |
issn | 2072-4292 |
language | English |
last_indexed | 2024-03-10T10:52:31Z |
publishDate | 2021-05-01 |
publisher | MDPI AG |
record_format | Article |
series | Remote Sensing |
spelling | doaj.art-4561a1e092e04fa4a51aebb1bb5c440f2023-11-21T22:07:40ZengMDPI AGRemote Sensing2072-42922021-05-011311214910.3390/rs13112149Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and WavelengthMarcello de Michele0Daniel Raucoules1Deborah Idier2Farid Smai3Michael Foumelis4BRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceBRGM, French Geological Survey, Direction des Risques et de la Prevention des Risques—3 av. C. Guillemin, 45000 Orleans, FranceIn this study, we present a new method called <i>BathySent</i> to retrieve shallow bathymetry from space that is based on the joint measurement of ocean wave celerity (c) and wavelength (λ). We developed the method to work with Sentinel 2 data, exploiting the time lag between two Sentinel 2 spectral bands, acquired quasi-simultaneously, from a single satellite dataset. Our method was based on the linear dispersion law, which related water depth to wave celerity and wavelength: when the water depth was less than about half the dominant wavelength, the wave celerity and wavelength decreased due to decreasing water depth (h) as the waves propagated towards the coast. Instead of using a best weighted <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> fit with the linear dispersion relation to retrieve h, we proposed solving the linear dispersion relation for each <inline-formula><math display="inline"><semantics><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">c</mi><mo>,</mo><mi mathvariant="sans-serif">λ</mi><mo stretchy="false">)</mo></mrow></semantics></math></inline-formula> pair to find multiple h-values within the same resolution cell. Then, we calculated the weighted averaged h-value for each resolution cell. To improve the precision of the final bathymetric map, we stacked the bathymetry values from <i>N</i>-different datasets acquired from the same study area on different dates. We first tested the algorithm on a set of images representing simulated ocean waves, then we applied it to a real set of Sentinel 2 data obtained of our study area, Gâvres peninsula (France, 47°,67 lat.; −3°35 lon.). A comparison with in situ bathymetry yielded good results from the synthetic images (r<sup>2</sup> = 0.9) and promising results with the Sentinel 2 images (r<sup>2</sup> = 0.7) in the 0–16 m depth zone.https://www.mdpi.com/2072-4292/13/11/2149bathymetrySentinel 2cross correlation |
spellingShingle | Marcello de Michele Daniel Raucoules Deborah Idier Farid Smai Michael Foumelis Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength Remote Sensing bathymetry Sentinel 2 cross correlation |
title | Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength |
title_full | Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength |
title_fullStr | Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength |
title_full_unstemmed | Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength |
title_short | Shallow Bathymetry from Multiple Sentinel 2 Images via the Joint Estimation of Wave Celerity and Wavelength |
title_sort | shallow bathymetry from multiple sentinel 2 images via the joint estimation of wave celerity and wavelength |
topic | bathymetry Sentinel 2 cross correlation |
url | https://www.mdpi.com/2072-4292/13/11/2149 |
work_keys_str_mv | AT marcellodemichele shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength AT danielraucoules shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength AT deborahidier shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength AT faridsmai shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength AT michaelfoumelis shallowbathymetryfrommultiplesentinel2imagesviathejointestimationofwavecelerityandwavelength |